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.

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

Local convergence analysis of L1-ADI scheme for two-dimensional reaction-subdiffusion equation

Efficient second-order ADI difference schemes for three-dimensional Riesz space-fractional diffusion equations

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

Convergence analysis of high-order compact alternating direction implicit schemes for the two-dimensional time fractional diffusion equation

Compact ADI Difference Scheme for the 2D Time Fractional Nonlinear Schrödinger Equation

An ADI Difference Scheme Based on Fractional Trapezoidal Rule for Fractional Integro-Differential Equation with a Weakly Singular Kernel.

Convergence analysis of a L1â€�ADI scheme for twoâ€�dimensional multiterm reactionâ€�subdiffusion equation

Finite Difference/Collocation Method for Two-Dimensional Sub-Diffusion Equation with Generalized Time Fractional Derivative

Efficient Numerical Algorithms for Three-Dimensional Fractional Partial Differential Equations

Three-Point Combined Compact Difference Schemes for Time-Fractional Advection-Diffusion Equations with Smooth Solutions

A fast finite difference method for 2D time variable fractional mobile/immobile equation

Fitted L1-ADI Scheme for Improving Convergence of Two-Dimensional Delay Fractional Equations

Finite Difference and Chebyshev Collocation for Time-Fractional and Riesz Space Distributed-Order Advection–Diffusion Equation with Time-Delay

An efficient numerical method for a time-fractional diffusion equation

Two Mixed Finite Element Methods for Time-Fractional Diffusion Equations

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

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

A high-order compact ADI finite difference scheme on uniform meshes for a weakly singular integro-differential equation in three space dimensions

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

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