A posteriori error analysis for approximations of time-fractional subdiffusion problems

Abstract

In this paper we consider a sub-diffusion problem where the fractional time derivative is approximated either by the L1 scheme or by Convolution Quadrature. We propose new interpretations of the numerical schemes which lead to a posteriori error estimates. Our approach is based on appropriate pointwise representations of the numerical schemes as perturbed evolution equations and on stability estimates for the evolution equation. A posteriori error estimates in $L^2(H)$ and $L^{\infty }(H)$ norms of optimal order are derived. Extensive numerical experiments indicate the reliability and the optimality of the estimators for the schemes considered, as well as their efficiency as error indicators driving adaptive mesh selection locating singularities of the problem.

arXiv commentaccepted for publication in Math. Comp. (Dec 2021)

References in corpus (5)

• Elliptic reconstruction and a posteriori error estimates for fully discrete linear parabolic problems
• An analysis of the L1 Scheme for the subdiffusion equation with nonsmooth data
• Fast summation by interval clustering for an evolution equation with memory
• A posteriori error estimates for leap-frog and cosine methods for second order evolution problems
• Well-posedness of the fractional Zener wave equation for heterogenous viscoelastic materials

Cited by in corpus (2)

• Fractional Dissipative PDEs
• An $\alpha$-robust and second-order accurate scheme for a subdiffusion equation