On an explicit finite difference method for fractional diffusion equations

Abstract

A numerical method to solve the fractional diffusion equation, which could also be easily extended to many other fractional dynamics equations, is considered. These fractional equations have been proposed in order to describe anomalous transport characterized by non-Markovian kinetics and the breakdown of Fick's law. In this paper we combine the forward time centered space (FTCS) method, well known for the numerical integration of ordinary diffusion equations, with the Grunwald-Letnikov definition of the fractional derivative operator to obtain an explicit fractional FTCS scheme for solving the fractional diffusion equation. The resulting method is amenable to a stability analysis a la von Neumann. We show that the analytical stability bounds are in excellent agreement with numerical tests. Comparison between exact analytical solutions and numerical predictions are made.

22 pages, 6 figures

Cited by in corpus (31)

• Compact difference schemes for the modified anomalous fractional sub-diffusion equation and the fractional diffusion-wave equation
• Numerical methods for time-fractional evolution equations with nonsmooth data: a concise overview
• Trapezoidal methods for fractional differential equations: theoretical and computational aspects
• Superconvergence of a discontinuous Galerkin method for fractional diffusion and wave equations
• On reaction-subdiffusion equations
• A Finite Difference Method with Non-uniform Timesteps for Fractional Diffusion Equations
• High-order numerical algorithms for Riesz derivatives via constructing new generating functions
• High-order Numerical Methods for Riesz Space Fractional Turbulent Diffusion Equation
• Efficient computation of the Grunwald-Letnikov fractional diffusion derivative using adaptive time step memory
• Well-posedness of time-fractional, advection-diffusion-reaction equations
• Mesoscopic description of reactions under anomalous diffusion: A case study
• High-order Compact Difference Schemes for the Modified Anomalous Subdiffusion Equation
• Positivity and Boundedness Preserving Schemes for Space-Time Fractional Predator-Prey Reaction-Diffusion Model
• A Gauss-Jacobi Kernel Compression Scheme for Fractional Differential Equations
• Numerical method for the time-fractional porous medium equation
• Invariant subspace method to the initial and boundary value problem of the higher dimensional nonlinear time-fractional PDEs
• A contour method for time-fractional PDEs and an application to fractional viscoelastic beam equations
• A new second-order midpoint approximation formula for Riemann-Liouville derivative: algorithm and its application
• The BEM and DRBEM schemes for the numerical solution of the two-dimensional time-fractional diffusion-wave equations
• Subdiffusive discrete time random walks via Monte Carlo and subordination
• Positivity and boundedness preserving schemes for the fractional reaction-diffusion equation
• A linear Galerkin numerical method for a quasilinear subdiffusion equation
• Analysis of the L1 scheme for fractional wave equations with nonsmooth data
• Second order scheme for self-similar solutions of a time-fractional porous medium equation on the half-line
• Analysis and implementation of collocation methods for fractional differential equations
• Effective numerical treatment of sub-diffusion equation with non-smooth solution
• Numerical methods for the two-dimensional Fokker-Planck equation governing the probability density function of the tempered fractional Brownian motion
• High-order Time Stepping Schemes for Semilinear Subdiffusion Equations
• Fast Stencil Computations using Fast Fourier Transforms
• Fractional Crank-Nicolson-Galerkin finite element methods for nonlinear time fractional parabolic problems with time delay
• A linearized second-order scheme for nonlinear time fractional Klein-Gordon type equations