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.

Analysis of a temporal discretization for a semilinear fractional diffusion equation

Error analysis of an L2-type method on graded meshes for semilinear subdiffusion equations

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

High-accuracy time discretization of stochastic fractional diffusion equation

Error analysis of the L1 method on graded and uniform meshes for a fractional-derivative problem in two and three dimensions

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

An -robust analysis of finite element method for space-time fractional diffusion equation

Analysis of the L1 scheme for fractional wave equations with nonsmooth data

Second‐order difference scheme for the time fractional mixed diffusion‐wave equation with initial weak regularity

α -robust analysis of fast and novel two-grid FEM with nonuniform L1 scheme for semilinear time-fractional variable coefficient diffusion equations

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

Numerical analysis of linear and nonlinear time-fractional subdiffusion equations

Some Grönwall inequalities for a class of discretizations of time fractional equations on nonuniform meshes

Two Mixed Finite Element Methods for Time-Fractional Diffusion Equations

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

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

Regularity Analysis and High-Order Time Stepping Scheme for Quasilinear Subdiffusion

Local convergence analysis of L1/finite element scheme for a constant delay reaction-subdiffusion equation with uniform time mesh

A Corrected L1 Method for a Time-Fractional Subdiffusion Equation

Local analysis of L1-finite difference method on graded meshes for multi-term two-dimensional time-fractional initial-boundary value problem with Neumann boundary conditions