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.

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

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

Error analysis for a fractional-derivative parabolic problem on quasi-graded meshes using barrier functions

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

Error analysis of a collocation method on graded meshes for a fractional Laplacian problem

α-Robust Error Analysis of L2-1σ Scheme on Graded Mesh for Time-Fractional Nonlocal Diffusion Equation

Analysis of a temporal discretization for a semilinear fractional diffusion equation

$α$-Robust Error Analysis of $L2$-$1_σ$ Scheme on Graded Mesh for Time-fractional Nonlocal Diffusion Equation

Error analysis of a collocation method on graded meshes for nonlocal diffusion problems with weakly singular kernels

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

Local error analysis of L1 scheme for time-fractional diffusion equation on a star-shaped pipe network

A Compact Difference Scheme on Graded Meshes for the Nonlinear Fractional Integro-differential Equation with Non-smooth Solutions

Optimal error analysis of a non-uniform IMEX-L1 finite element method for time fractional PDEs and PIDEs

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

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 Methods for Caputo–Hadamard Fractional Differential Equations with Graded and Non-Uniform Meshes

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

Error estimate of GL‐ADI scheme for 2D multiterm nonlinear time‐fractional subdiffusion equation

Analysis of stability and convergence for L-type formulas combined with a spatial finite element method for solving subdiffusion problems

Analysis and numerical solution of a Riemann-Liouville fractional derivative two-point boundary value problem

$L2$-$1_σ$ Scheme on a Graded Mesh for a Multi-term Time-fractional Nonlocal Parabolic Problem