Quenching of combustion explosion model with balanced space-fractional derivative

Balanced space-fractional derivative is usually applied in modelling the state-dependence, isotropy, and anisotropy in diffusion phenomena. In this paper, we introduce a class of space-fractional reaction-diffusion model with singular source term arising in combustion process. The fractional derivative employed in this model is defined in the sum of left-sided and right-sided Riemann-Liouville fractional derivatives. With assistance of Kaplan's first eigenvalue method, we prove that the classic solution of this model may not be globally well-defined, and the heat conduction governed by this model depends on the order of fractional derivative, the parameters in the equation, and the length of spatial interval. Finite difference method is implemented to solve this model, and an adaptive strategy is applied to improve the computational efficiency. The positivity, monotonicity, and stability of the numerical scheme are discussed. Numerical simulation and observation of the quenching and stationary solutions coincide the theoretical studies.     

Recovering the initial distribution for space-fractional diffusion equation by a logarithmic regularization method

In this paper, the backward problem for space-fractional diffusion equation is investigated. We proposed a so-called logarithmic regularization method to solve it. Based on the conditional stability and an a posteriori regularization parameter choice rule, the convergence rate estimates are given under a-priori bound assumption for the exact solution.     

The locally extrapolated exponential splitting scheme for multi-dimensional nonlinear space-fractional Schrödinger equations

An efficient local extrapolation of the exponential operator splitting scheme is introduced to solve the multi-dimensional space-fractional nonlinear Schrödinger equations. Stability of the scheme is examined by investigating its amplification factor and by plotting the boundaries of the stability regions. Empirical convergence analysis and calculation of the local truncation error exhibit the second-order accuracy of the proposed scheme. The performance and reliability of the proposed scheme are tested by implementing it on two- and three-dimensional space-fractional nonlinear Schrödinger equations including the space-fractional Gross-Pitaevskii equation, which is used to model optical solitons in graded-index fibers.     

Fractional diffusion-type equations with exponential and logarithmic differential operators

We deal with some extensions of the space-fractional diffusion equation, which is satisfied by the density of a stable process (see Mainardi et al. (2001)): the first equation considered here is obtained by adding an exponential differential (or shift) operator expressed in terms of the Riesz–Feller derivative. We prove that this produces a random component in the time-argument of the corresponding stable process, which is represented by the so-called Poisson process with drift. Analogously, if we add, to the space-fractional diffusion equation, a logarithmic differential operator involving the Riesz-derivative, we obtain, as a solution, the transition semigroup of a stable process subordinated by an independent gamma subordinator with drift. Finally, we show that an extension of the space-fractional diffusion equation, containing both the fractional shift operator and the Feller integral, is satisfied by the transition density of the process obtained by time-changing the stable process with an independent linear birth process with drift.     

A linear,symmetric and energy-conservative scheme for the space-fractional Klein–Gordon–Schrödinger equations

In this paper, we propose an efficient numerical scheme for the space-fractional Klein–Gordon–Schrödinger (SFKGS) equations. Motivated by the “Invariant Energy Quadratization” (IEQ) approach, we introduce two auxiliary variables to transform the SFKGS system into a new equivalent system in which the time derivative is discretized by the Crank–Nicolson method, and the space discretization is based on the Fourier spectral method. Consequently, the numerical scheme shares two good features. The first feature is that the nonlinear terms are treated semi-explicitly and a linear symmetric system is solved at each time step. The second feature is the energy conservation at the discrete level. These two advantages are proved by the theoretical analysis and illustrated by a given numerical example.     

Large-Time Behavior in Population Dynamics with Offspring Production at Fixed Ages,Child Care,and Spatial Dispersal

We present a density-dependent population dynamics model with age-dependence, child care, and spatial dispersal. The population consists of the young (under maternal care), juvenile, and adult (producing offsprings at fixed ages or of post-reproductive age) classes. Death moduli of the juvenile and adult individuals are decomposed into two-term sums. The first sum represents the death rate by natural causes and by those that do not depend on the population spatial density, while the other one describes the environmental influence depending on the spatial density of the juvenile and adult individuals. The steady-state and a class of separable solutions are considered, and the large-time behavior of separable solutions is analyzed for the stationary vital rates. The asymptotic behavior of nondispersing semelparous species is also examined.     

Differential quadrature method for space-fractional diffusion equations on 2D irregular domains

In mathematical physics, the space-fractional diffusion equations are of particular interest in the studies of physical phenomena modelled by Lévy processes, which are sometimes called super-diffusion equations. In this article, we develop the differential quadrature (DQ) methods for solving the 2D space-fractional diffusion equations on irregular domains. The methods in presence reduce the original equation into a set of ordinary differential equations (ODEs) by introducing valid DQ formulations to fractional directional derivatives based on the functional values at scattered nodal points on problem domain. The required weighted coefficients are calculated by using radial basis functions (RBFs) as trial functions, and the resultant ODEs are discretized by the Crank-Nicolson scheme. The main advantages of our methods lie in their flexibility and applicability to arbitrary domains. A series of illustrated examples are finally provided to support these points.     

Synchronism in population networks with non linear coupling

A quite general spatially explicit metapopulation model featuring density-dependent dispersal is proposed. A study of the stability of synchronized attractors is done based on an explicit calculation of the transverse Lyapunov number of these attractors. An analytic expression for the transverse Lyapunov number depending on the Lyapunov number of the synchronous trajectory, the eigenvalues of the network configuration matrix and the function modelling the density-dependent dispersal (the number of migrants as a function of the local density) is presented. Numerical simulations are also performed upon selecting biologically relevant features, such as local dynamics, network topology and density-dependent dispersal mechanism.     

A new reliable algorithm based on the sinc function for the time fractional diffusion equation

This paper deals with the numerical solution of time fractional diffusion equation. In this work, we consider the fractional derivative in the sense of Riemann-Liouville. At first, the time fractional derivative is discretized by integrating both sides of the equation with respect to the time variable and we arrive at a semi–discrete scheme. The stability and convergence of time discretized scheme are proven by using the energy method. Also we show that the convergence order of this scheme is O(τ2?α). Then we use the sinc collocation method to approximate the solution of semi–discrete scheme and show that the problem is reduced to a Sylvester matrix equation. Besides by performing some theorems, the exponential convergence rate of sinc method is illustrated. The numerical experiments are presented to show the excellent behavior and high accuracy of the proposed hybrid method in comparison with some other well known methods.     

The fractional Tikhonov regularization method for simultaneous inversion of the source term and initial value in a space-fractional Allen-Cahn equation

In this paper, we consider the inverse problem for identifying the source term and initial value simultaneously in a space-fractional Allen-Cahn equation. This problem is ill-posed, i.e., the solution of this problem does not depend continuously on the data. The fractional Tikhonov method is used to solve this problem. Under the a priori and the a posteriori regularization parameter choice rules, the error estimates between the regularization solutions and the exact solutions are obtained, respectively. Different numerical examples are presented to illustrate the validity and effectiveness of our method.     

Stability and Convergence of an Implicit Difference Approximation for the Space Riesz Fractional Reaction-Dispersion Equation

In this paper,we consider a Riesz space-fractional reaction-dispersion equation (RSFRDE).The RSFRDE is obtained from the classical reaction-dispersion equation by replacing the second-order space derivative with a Riesz derivative of orderβ∈(1,2]. We propose an implicit finite difference approximation for RSFRDE.The stability and convergence of the finite difference approximations are analyzed.Numerical results are found in good agreement with the theoretical analysis.     

Time-space fractional stochastic Ginzburg-Landau equation driven by Gaussian white noise

The current article is devoted to the time-spatial regularity of the nonlocal stochastic convolution for Caputo-type time fractional nonlocal Ornstein-Ulenbeck equations. The dependence of the order of time-fractional derivative, the order of the space-fractional derivative, and the regularity of the initial data are revealed. The global existence and uniqueness of the mild solutions for time-space fractional complex Ginzburg-Landau equation driven by Gaussian white noise are established.     

M-ary signal detection via a bistable system in the presence of Lévy noise

We study the problem of the M-ary signal detection via a bistable detector in the presence of Lévy noise. Based on the numerical solution of the space-fractional Fokker–Planck equation, the theoretical bit error rate is defined and used in the optimal detector design. The accuracy of the theoretical results are verified by the Monte Carlo simulations. It is shown that, with the same noise intensity, the optimal bistable detector performs better with the decreasing Lévy index α. Therefore, Lévy noise plays a more positive role in the nonlinear M-ary signal detection problem, compared to Gaussian noise.     

A semidefinite programming approach for solving fractional optimal control problems

This paper presents a numerical scheme for solving fractional optimal control. The fractional derivative in this problem is in the Riemann–Liouville sense. The proposed method, based upon the method of moments, converts the fractional optimal control problem to a semidefinite optimization problem; namely, the nonlinear optimal control problem is converted to a convex optimization problem. The Grunwald–Letnikov formula is also used as an approximation for fractional derivative. The solution of fractional optimal control problem is found by solving the semidefinite optimization problem. Finally, numerical examples are presented to show the performance of the method.     

A New Efficient Transform Mechanism with Convergence Analysis of the Space-Fractional Telegraph Equations

This article"s goal is to investigate the space-fractional telegraph equation using an effective method called the Adomian natural decomposition method (ANDM), which is a combination of the Adomian decomposition method (ADM) and the natural transform method (NTM). Using the Banach fixed point theorem, we explore proofs for the existence and uniqueness theorems applying it to a nonlinear differential equation. Using our method, exact solutions of the space-fractional telegraph equation and time-fractional diffusion problems have been obtained. To demonstrate the effectiveness of the suggested scheme, four examples are provided.     

Bifurcation analysis in a diffusive Logistic population model with two delayed density-dependent feedback terms

The present paper is concerned with a diffusive population model of Logistic type with an instantaneous density-dependent term and two delayed density-dependent terms and subject to the zero-Dirichlet boundary condition. By regarding the delay as the bifurcation parameter and analyzing in detail the associated eigenvalue problem, the local asymptotic stability and the existence of Hopf bifurcation for the sufficiently small positive steady state solution are shown. It is found that under the suitable condition, the positive steady state solution of the model will become ultimately unstable after a single stability switch (or change) at a certain critical value of delay through a Hopf bifurcation. However, under the other condition, the positive steady state solution of the model will become ultimately unstable after multiple stability switches at some certain critical values of delay through Hopf bifurcations. In addition, the direction of the above Hopf bifurcations and the stability of the bifurcating periodic solutions are analyzed by means of the center manifold theory and normal form method for partial functional differential equations. Finally, in order to illustrate the correction of the obtained theoretical results, some numerical simulations are also carried out.     

A mathematical analysis of an age-space-structured autosomal diploid population dynamics model with random mating and females' pregnancy

We discuss a deterministic model of the age-structured autosomal polylocal multiallelic diploid population dynamics that takes into account random mating of sexes, females' pregnancy, and its dispersal in the whole space. This model generalizes the previous one by taking into account the spatial dispersal whose mechanism is described by the general linear elliptic differential operator of the second order. The population consists of male, single (nonfertilized) female, and fertilized female subclasses. Using the fundamental solution method for the uniformly parabolic second-order differential operator with bounded Hölder continuous coefficients, we prove the existence and uniqueness theorem for the classic solution of the Cauchy problem for this model. In the case where dispersal moduli of fertilized females do not depend on the age of mated males, we analyze the population growth and decay. Mutation is not consisdered in this paper.     

Stable numerical differentiation for the second order derivatives

Consider a numerical differential problem, which aims to compute the second order derivative of a function stably from its given noisy data. For this ill-posed problem, we introduce the Lavrent′ev regularization scheme by reformulating this differentiation problem as an integral equation of the first kind. The advantage of this proposed scheme is that we can give the regularizing solution by an explicit integral expression, therefore it is easy to be implemented. The a-priori and a-posterior choice strategies for the regularization parameter are considered, with convergence analysis and error estimate of the regularizing solution for noisy data based on the integral operator decomposition. The validity of the proposed scheme is shown by several numerical examples.     

On reflecting boundary conditions for space-fractional equations on a finite interval: Proof of the matrix transfer technique

Even in the one-dimensional case, dealing with the analysis of space-fractional differential equations on finite domains is a difficult issue. On a finite interval coupled with zero flux boundary conditions, different approaches have been proposed to define a space-fractional differential operator and to compute the solution to the corresponding fractional problem, but to the best of our knowledge, a clear relationship between these strategies is yet to be established. Here, by using the theory of α-stable symmetric Lévy flights and the master equation, we derive a discrete representation of the non-local operator embedding in its definition the concept of reflecting boundary conditions. We refer to this discrete operator as the reflection matrix and provide (and prove) a theorem on the analytic expression of its eigenvalues and eigenvectors. We then use this result to compare the reflection matrix to the discrete operator defined via the matrix transfer technique, and establish the validity of the latter technique in producing the correct solution to a space-fractional differential equation on a finite interval with reflecting boundary conditions. We finally discuss and emphasize the challenges in the generalisation of the proposed result to more than one spatial dimension.     

Solving a final value fractional diffusion problem by boundary condition regularization

The evolution process of fractional order describes some phenomenon of anomalous diffusion and transport dynamics in complex system. The equation containing fractional derivatives provides a suitable mathematical model for describing such a process. The initial boundary value problem is hard to solve due to the nonlocal property of the fractional order derivative. We consider a final value problem in a bounded domain for fractional evolution process with respect to time, which means to recover the initial state for some slow diffusion process from its present status. For this ill-posed problem, we construct a regularizing solution using quasi-reversible method. The well-posedness of the regularizing solution as well as the convergence property is rigorously analyzed. The advantage of the proposed scheme is that the regularizing solution is of the explicit analytic solution and therefore is easy to be implemented. Numerical examples are presented to show the validity of the proposed scheme.