The transport dynamics in complex systems governing by anomalous diffusion modelled with Riesz fractional partial differential equations

In this paper, numerical solutions of fractional Fokker–Planck equations with Riesz space fractional derivatives have been developed. Here, the fractional Fokker–Planck equations have been considered in a finite domain. In order to deal with the Riesz fractional derivative operator, shifted Grünwald approximation and fractional centred difference approaches have been used. The explicit finite difference method and Crank–Nicolson implicit method have been applied to obtain the numerical solutions of fractional diffusion equation and fractional Fokker–Planck equations, respectively. Numerical results are presented to demonstrate the accuracy and effectiveness of the proposed numerical solution techniques. Copyright © 2016 John Wiley & Sons, Ltd.     

Gauss‐Lobatto‐Legendre‐Birkhoff pseudospectral approximations for the multi‐term time fractional diffusion‐wave equation with Neumann boundaryconditions

A second‐order finite difference/pseudospectral scheme is proposed for numerical approximation of multi‐term time fractional diffusion‐wave equation with Neumann boundary conditions. The scheme is based upon the weighted and shifted Grünwald difference operators approximation of the time fractional calculus and Gauss‐Lobatto‐Legendre‐Birkhoff (GLLB) pseudospectral method for spatial discretization. The unconditionally stability and convergence of the scheme are rigorously proved. Numerical examples are carried out to verify theoretical results.     

An efficient numerical treatment of fourth‐order fractional diffusion‐wave problems

In this paper, we consider the numerical treatment of a fourth‐order fractional diffusion‐wave problem. Our proposed method includes the use of parametric quintic spline in the spatial dimension and the weighted shifted Grünwald‐Letnikov approximation of fractional integral. The solvability, stability, and convergence of the numerical scheme are rigorously proved. It is shown that the theoretical convergence order improves those of earlier work. Simulation is further carried out to demonstrate the numerical efficiency of the proposed scheme and to compare with other methods.     

An efficient difference scheme for the coupled nonlinear fractional Ginzburg–Landau equations with the fractional Laplacian

In this article, an efficient difference scheme for the coupled fractional Ginzburg–Landau equations with the fractional Laplacian is studied. We construct the discrete scheme based on the implicit midpoint method in time and a weighted and shifted Grünwald difference method in space. Then, we prove that the scheme is uniquely solvable, and the numerical solutions are bounded and unconditionally convergent in the norm. Finally, numerical tests are given to confirm the theoretical results and show the effectiveness of the scheme.     

Numerical methods for fractional partial differential equations with Riesz space fractional derivatives

In this paper, we consider the numerical solution of a fractional partial differential equation with Riesz space fractional derivatives (FPDE-RSFD) on a finite domain. Two types of FPDE-RSFD are considered: the Riesz fractional diffusion equation (RFDE) and the Riesz fractional advection–dispersion equation (RFADE). The RFDE is obtained from the standard diffusion equation by replacing the second-order space derivative with the Riesz fractional derivative of order α(1,2]. The RFADE is obtained from the standard advection–dispersion equation by replacing the first-order and second-order space derivatives with the Riesz fractional derivatives of order β(0,1) and of order α(1,2], respectively. Firstly, analytic solutions of both the RFDE and RFADE are derived. Secondly, three numerical methods are provided to deal with the Riesz space fractional derivatives, namely, the L1/L2-approximation method, the standard/shifted Grünwald method, and the matrix transform method (MTM). Thirdly, the RFDE and RFADE are transformed into a system of ordinary differential equations, which is then solved by the method of lines. Finally, numerical results are given, which demonstrate the effectiveness and convergence of the three numerical methods.     

Finite difference methods for fractional dispersion equations

In this paper, the fractional weighted average finite difference method for space-fractional advection–dispersion equation is proposed, which is based on shifted Grünwald formula. This method is unconditionally stable, consistent and convergent. A numerical example is given, and the numerical results verify the theoretical conclusions.     

Finite element method for Grwünwald–Letnikov time-fractional partial differential equation

In this article, we consider the finite element methods (FEM) for Grwünwald–Letnikov time-fractional diffusion equation, which is obtained from the standard two-dimensional diffusion equation by replacing the first-order time derivative with a fractional derivative (of order α, with 0?h ^r+1?+?τ^2-α), where h, τ and r are the space step size, time step size and polynomial degree, respectively. A numerical example is presented to verify the order of convergence.     

A fractional order epidemic model for the simulation of outbreaks of influenza A(H1N1)

In this paper, we propose a nonlinear fractional order model in order to explain and understand the outbreaks of influenza A(H1N1). In the fractional model, the next state depends not only upon its current state but also upon all of its historical states. Thus, the fractional model is more general than the classical epidemic models. In order to deal with the fractional derivatives of the model, we rely on the Caputo operator and on the Grünwald–Letnikov method to numerically approximate the fractional derivatives. We conclude that the nonlinear fractional order epidemic model is well suited to provide numerical results that agree very well with real data of influenza A(H1N1) at the level population. In addition, the proposed model can provide useful information for the understanding, prediction, and control of the transmission of different epidemics worldwide. Copyright © 2013 John Wiley & Sons, Ltd.     

High order schemes for the tempered fractional diffusion equations

Lévy flight models whose jumps have infinite moments are mathematically used to describe the superdiffusion in complex systems. Exponentially tempering Lévy measure of Lévy flights leads to the tempered stable Lévy processes which combine both the α-stable and Gaussian trends; and the very large jumps are unlikely and all their moments exist. The probability density functions of the tempered stable Lévy processes solve the tempered fractional diffusion equation. This paper focuses on designing the high order difference schemes for the tempered fractional diffusion equation on bounded domain. The high order difference approximations, called the tempered and weighted and shifted Grünwald difference (tempered-WSGD) operators, in space are obtained by using the properties of the tempered fractional calculus and weighting and shifting their first order Grünwald type difference approximations. And the Crank-Nicolson discretization is used in the time direction. The stability and convergence of the presented numerical schemes are established; and the numerical experiments are performed to confirm the theoretical results and testify the effectiveness of the schemes.     

A series of high‐order quasi‐compact schemes for space fractional diffusion equations based on the superconvergent approximations for fractional derivatives

Based on the superconvergent approximation at some point (depending on the fractional order α, but not belonging to the mesh points) for Grünwald discretization to fractional derivative, we develop a series of high‐order quasi‐compact schemes for space fractional diffusion equations. Because of the quasi‐compactness of the derived schemes, no points beyond the domain are used for all the high‐order schemes including second‐order, third‐order, fourth‐order, and even higher‐order schemes; moreover, the algebraic equations for all the high‐order schemes have the completely same matrix structure. The stability and convergence analysis for some typical schemes are made; the techniques of treating the fractional derivatives with nonhomogeneous boundaries are introduced; and extensive numerical experiments are performed to confirm the theoretical analysis or verify the convergence orders. © 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 1345–1381, 2015     

Numerical solutions of nonlinear fractional Schrödinger equations using nonstandard discretizations

In this article, numerical study for both nonlinear space‐fractional Schrödinger equation and the coupled nonlinear space‐fractional Schrödinger system is presented. We offer here the weighted average nonstandard finite difference method (WANSFDM) as a novel numerical technique to study such kinds of partial differential equations. The space fractional derivative is described in the sense of the quantum Riesz‐Feller definition. Stability analysis of the proposed method is studied. To show that this method is reliable and computationally efficient different numerical examples are provided. We expect that the proposed schemes can be applicable to different systems of fractional partial differential equations. © 2016 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 1399–1419, 2017     

Maximal ℓp‐regularity of multiterm fractional equations with delay

We provide a characterization for the existence and uniqueness of solutions in the space of vector‐valued sequences ${?}_p(?,X)$ for the multiterm fractional delayed model in the form $${\mathrm{\Delta }}^{\alpha }u(n)+\lambda {\mathrm{\Delta }}^{\beta }u(n)=Au(n)+u(n?\tau )+f(n),n\in ?,\alpha ,\beta \in {?}_{+},\tau \in ?,\lambda \in ?,$$ where X is a Banach space, A is a closed linear operator with domain D(A) defined on X, $f\in {?}_p(?,X)$ and Δ^Γ denotes the Grünwald–Letkinov fractional derivative of order Γ > 0. We also give some conditions to ensure the existence of solutions when adding nonlinearities. Finally, we illustrate our results with an example given by a general abstract nonlinear model that includes the fractional Fisher equation with delay.     

Finite difference spectral approximation for the time‐space fractional telegraph equation and its parameter estimation

The object of this paper is to present the numerical solution of the time‐space fractional telegraph equation. The proposed method is based on the finite difference scheme in temporal direction and Fourier spectral method in spatial direction. The fast Fourier transform (FFT) technique is applied to practical computation. The stability and convergence analysis are strictly proven, which shows that this method is stable and convergent with (2?α) order accuracy in time and spectral accuracy in space. Moreover, the Levenberg‐Marquardt (L‐M) iterative method is employed for the parameter estimation. Finally, some numerical examples are given to confirm the theoretical analysis.     

Respectively scaled HSS iteration methods for solving discretized spatial fractional diffusion equations

For the discrete linear systems resulted from the discretization of the one‐dimensional anisotropic spatial fractional diffusion equations of variable coefficients with the shifted finite‐difference formulas of the Grünwald–Letnikov type, we propose a class of respectively scaled Hermitian and skew‐Hermitian splitting iteration method and establish its asymptotic convergence theory. The corresponding induced matrix splitting preconditioner, through further replacements of the involved Toeplitz matrices with certain circulant matrices, leads to an economic variant that can be executed by fast Fourier transforms. Both theoretical analysis and numerical implementations show that this fast respectively scaled Hermitian and skew‐Hermitian splitting preconditioner can significantly improve the computational efficiency of the Krylov subspace iteration methods employed as effective linear solvers for the target discrete linear systems.     

ON THE RATE OF WEIGHTED Lp CONVERGENCE BY THE GRUENWALD INTERPOLATORY OPERATORS

This paper gives the weighted Lp convergence rate estimations of the Gruenwald interpolatory polynomials based on the zeros of Chebyshev polymomials of the first kind,and proves that the order of the estimations is optimal for p≥1.     

An inverse problem for space‐fractional backward diffusion problem

In this paper, an inverse problem for space‐fractional backward diffusion equation, which is highly ill‐posed, is considered. This problem is obtained from the classical diffusion equation by replacing the second‐order space derivative with a Riesz–Feller derivative of order α ∈ (0,2]. We show that such a problem is severely ill‐posed, and further present a simplified Tikhonov regularization method to deal with this problem. Convergence estimate is presented under a priori choice of regularization parameter. Numerical experiments are given to illustrate the accuracy and efficiency of the proposed method. Copyright © 2013 John Wiley & Sons, Ltd.     

An explicit/implicit Galerkin domain decomposition procedure for parabolic integro-differential equations

In this paper, a two-dimensional fractional advection-dispersion equation (2D-FADE) with variable coefficients on a finite domain is considered. We use a new technique of combination of the Alternating Directions Implicit-Euler method (ADI-Euler), the unshifted Grünwald formula for the advection term, the right-shifted Grünwald formula for the diffusion term, and a Richardson extrapolation to establish an unconditionally stable second order accurate difference method. Stability, consistency and convergence of the ADI-Euler method for 2D-FADE are examined. A numerical example with known exact solution is also presented, and the behavior of the error is analyzed to verify the order of convergence of the ADI-Euler method and the extrapolated ADI-Euler method.     

High‐order compact difference schemes for the modified anomalous subdiffusion equation

In this article, two kinds of high‐order compact finite difference schemes for second‐order derivative are developed. Then a second‐order numerical scheme for a Riemann–Liouvile derivative is established based on a fractional centered difference operator. We apply these methods to a fractional anomalous subdiffusion equation to construct two kinds of novel numerical schemes. The solvability, stability, and convergence analysis of these difference schemes are studied by using Fourier method. The convergence orders of these numerical schemes are and , respectively. Finally, numerical experiments are displayed which are in line with the theoretical analysis. © 2015 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 32: 213–242, 2016     

Efficient and accurate spectral method for the time-fractional dual-phase-lag heat transfer model and its parameter estimation

In the current paper, a heat transfer model is suggested based on a time-fractional dual-phase-lag (DPL) model. We discuss the model in two parts, the direct problem and the inverse problem. Firstly, for solving it, the finite difference/Legendre spectral method is constructed. In the temporal direction, we employ the weighted and shifted Grünwald approximation, which can achieve second order convergence. In the spatial direction, the Legendre spectral method is used, it can obtain spectral accuracy. The stability and convergence are theoretically analyzed. For the inverse problem, the Bayesian method is used to construct an algorithm to estimate the four parameters for the model, namely, the time-fractional order α, the time-fractional order β, the delay time τ_T, and the relaxation time τ_q. Next, numerical experiments are provided to demonstrate the effectiveness of our scheme, with the values of τ_q and τ_T for processed meat employed. We also make a comparison with another method. The data obtained for the direct problem are used in the parameter estimation. The paper provides an accurate and efficient numerical method for the time-fractional DPL model.     

Gegenbauer spectral method for time‐fractional convection–diffusion equations with variable coefficients

In this paper, we study the numerical solution to time‐fractional partial differential equations with variable coefficients that involve temporal Caputo derivative. A spectral method based on Gegenbauer polynomials is taken for approximating the solution of the given time‐fractional partial differential equation in time and a collocation method in space. The suggested method reduces this type of equation to the solution of a linear algebraic system. Finally, some numerical examples are presented to illustrate the efficiency and accuracy of the proposed method. Copyright © 2014 John Wiley & Sons, Ltd.