An $\alpha$-robust and second-order accurate scheme for a subdiffusion equation
Numerical Analysis
2024-07-10 v2 Numerical Analysis
Abstract
We investigate a second-order accurate time-stepping scheme for solving a time-fractional diffusion equation with a Caputo derivative of order~. The basic idea of our scheme is based on local integration followed by linear interpolation. It reduces to the standard Crank--Nicolson scheme in the classical diffusion case, that is, as . Using a novel approach, we show that the proposed scheme is -robust and second-order accurate in the -norm, assuming a suitable time-graded mesh. For completeness, we use the Galerkin finite element method for the spatial discretization and discuss the error analysis under reasonable regularity assumptions on the given data. Some numerical results are presented at the end.
Cite
@article{arxiv.2401.04946,
title = {An $\alpha$-robust and second-order accurate scheme for a subdiffusion equation},
author = {Kassem Mustapha and William McLean and Josef Dick},
journal= {arXiv preprint arXiv:2401.04946},
year = {2024}
}