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.

A Comparative Study of Efficient Numerical Schemes for Time-Fractional Subdiffusion Equation Involving Singularity

Novel numerical methods based on graded, adaptive and uniform meshes for a time-fractional advection-diffusion equation subjected to weakly singular solution

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

High order approximation to Caputo derivative on graded mesh and time-fractional diffusion equation for non-smooth solutions

Efficient numerical schemes based on the cubic B-spline collocation method for time-fractional partial integro-differential equations of Volterra type

Fast compact finite difference schemes on graded meshes for fourth-order multi-term fractional sub-diffusion equations with the first Dirichlet boundary conditions

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

An Efficient Numerical Scheme for a Time-Fractional Burgers Equation with Caputo-Prabhakar Derivative

A high order numerical scheme for solving a class of non‐homogeneous time‐fractional reaction diffusion equation

A fast modified finite difference method for time fractional diffusion equations with weakly singular solution

Numerical investigation of two fractional operators for time fractional delay differential equation

An efficient computational technique for solving a time-fractional reaction-subdiffusion model in 2D space

Two Finite Difference Schemes for Multi-Dimensional Fractional Wave Equations with Weakly Singular Solutions

A Highly Efficient Numerical Method for the Time-Fractional Diffusion Equation on Unbounded Domains

A high-order compact finite difference method on nonuniform time meshes for variable coefficient reaction–subdiffusion problems with a weak initial singularity

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

L1 scheme on graded mesh for subdiffusion equation with nonlocal diffusion term

A numerical study on solving a fractional time-space diffusion equation via the finite difference method

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

An efficient numerical method for a time-fractional diffusion equation

Error analysis of a high-order fully discrete method for two-dimensional time-fractional convection-diffusion equations exhibiting weak initial singularity