Abstract:The subdiffusion equation with a Caputo fractional derivative of order $\alpha\in(0,1)$ in time arises in a wide variety of practical applications, and it is often adopted to model anomalous subdiffusion processes in heterogeneous media. The L1 scheme is one of the most popular and successful numerical methods for discretizing the Caputo fractional derivative in time. The scheme was analyzed earlier independently by Lin and Xu (2007) and Sun and Wu (2006), and an $O(\tau^{2-\alpha})$ convergence rate was established, under the assumption that the solution is twice continuously differentiable in time. However, in view of the smoothing property of the subdiffusion equation, this regularity condition is restrictive, since it does not hold even for the homogeneous problem with a smooth initial data. In this work, we revisit the error analysis of the scheme, and establish an $O(\tau)$ convergence rate for both smooth and nonsmooth initial data. The analysis is valid for more general sectorial operators. In particular, the L1 scheme is applied to one-dimensional space-time fractional diffusion equations, which involves also a Riemann-Liouville derivative of order $\beta\in(3/2,2)$ in space, and error estimates are provided for the fully discrete scheme. Numerical experiments are provided to verify the sharpness of the error estimates, and robustness of the scheme with respect to data regularity.

The Stability and Convergence Analysis of Finite Difference Methods for the Fractional Neutron Diffusion Equation

The H1-error Analysis of the Finite Element Method for Solving the Fractional Diffusion Equation

Stability and convergence of finite difference schemes for a class of time-fractional sub-diffusion equations based on certain superconvergence

Finite Difference Scheme on Graded Meshes to the Time-Fractional Neutron Diffusion Equation with Non-Smooth Solutions

A New Analytical Technique of the L-type Difference Schemes for Time Fractional Mixed Sub-Diffusion and Diffusion-Wave Equations.

Uniformly Stable Explicitly Solvable Finite Difference Method for Fractional Diffusion Equations

Efficient Numerical Schemes for Fractional Sub-Diffusion Equation with the Spatially Variable Coefficient

Finite Difference Methods for the Time Fractional Diffusion Equation on Non-Uniform Meshes.

Two Mixed Finite Element Methods for Time-Fractional Diffusion Equations

Compact Difference Scheme for the Fractional Sub-Diffusion Equation with Neumann Boundary Conditions

The Finite Difference/finite Volume Method for Solving the Fractional Diffusion Equation

Numerical Algorithm With High Spatial Accuracy for the Fractional Diffusion-Wave Equation With Neumann Boundary Conditions

Numerical Methods for Time-Fractional Convection-Diffusion Problems with High-Order Accuracy

A Compact Difference Scheme for Time Fractional Diffusion Equation with Neumann Boundary Conditions.

An analysis of the L1 Scheme for the subdiffusion equation with nonsmooth data

Investigations on Several High-Order ADI Methods for Time-Space Fractional Diffusion Equation

An efficient numerical method for a time-fractional diffusion equation

Finite Difference Scheme for Multi-Term Variable-Order Fractional Diffusion Equation

Finite Difference Schemes for the Variable Coefficients Single and Multi-Term Time-Fractional Diffusion Equations with Non-Smooth Solutions on Graded and Uniform Meshes

Stability and Convergence of an Effective Numerical Method for the Time-Space Fractional Fokker-Planck Equation with A Nonlinear Source Term

Efficient L1-ADI finite difference method for the two-dimensional nonlinear time-fractional diffusion equation