Unconditional convergence of a fast two-level linearized algorithm for semilinear subdiffusion equations

Abstract

A fast two-level linearized scheme with unequal time-steps is constructed and analyzed for an initial-boundary-value problem of semilinear subdiffusion equations. The two-level fast L1 formula of the Caputo derivative is derived based on the sum-of-exponentials technique. The resulting fast algorithm is computationally efficient in long-time simulations because it significantly reduces the computational cost $O(M{N}^2)$ and storage $O(MN)$ for the standard L1 formula to $O(MN\log N)$ and $O(M\log N)$, respectively, for $M$ grid points in space and $N$ levels in time. The nonuniform time mesh would be graded to handle the typical singularity of the solution near the time $t=0$, and Newton linearization is used to approximate the nonlinearity term. Our analysis relies on three tools: a new discrete fractional Grönwall inequality, a global consistency analysis and a discrete $H^2$ energy method. A sharp error estimate reflecting the regularity of solution is established without any restriction on the relative diameters of the temporal and spatial mesh sizes. Numerical examples are provided to demonstrate the effectiveness of our approach and the sharpness of error analysis.

arXiv comment23 pages, 5 figures

References in corpus (3)

• A discrete Grönwall inequality with application to numerical schemes for subdiffusion problems
• Unconditional convergence and optimal error estimates of a Galerkin-mixed FEM for incompressible miscible flow in porous media
• Fast summation by interval clustering for an evolution equation with memory

Cited by in corpus (10)

• A second-order and nonuniform time-stepping maximum-principle preserving scheme for time-fractional Allen-Cahn equations
• An energy stable and maximum bound preserving scheme with variable time steps for time fractional Allen-Cahn equation
• Adaptive second-order Crank-Nicolson time-stepping schemes for time fractional molecular beam epitaxial growth models
• A preconditioning technique for an all-at-once system from Volterra subdiffusion equations with graded time steps
• Simple maximum-principle preserving time-stepping methods for time-fractional Allen-Cahn equation
• Adaptive linear second-order energy stable schemes for time-fractional Allen-Cahn equation with volume constraint
• The variable-step L1 scheme preserving a compatible energy law for time-fractional Allen-Cahn equation
• A Non-uniform Time-stepping Convex Splitting Scheme for the Time-fractional Cahn-Hilliard Equation
• Compatible $L^2$ norm convergence of variable-step L1 scheme for the time-fractional MBE mobel with slope selection
• An implicit integration factor method for a kind of spatial fractional diffusion equations