- Journal Home
- Volume 36 - 2024
- Volume 35 - 2024
- Volume 34 - 2023
- Volume 33 - 2023
- Volume 32 - 2022
- Volume 31 - 2022
- Volume 30 - 2021
- Volume 29 - 2021
- Volume 28 - 2020
- Volume 27 - 2020
- Volume 26 - 2019
- Volume 25 - 2019
- Volume 24 - 2018
- Volume 23 - 2018
- Volume 22 - 2017
- Volume 21 - 2017
- Volume 20 - 2016
- Volume 19 - 2016
- Volume 18 - 2015
- Volume 17 - 2015
- Volume 16 - 2014
- Volume 15 - 2014
- Volume 14 - 2013
- Volume 13 - 2013
- Volume 12 - 2012
- Volume 11 - 2012
- Volume 10 - 2011
- Volume 9 - 2011
- Volume 8 - 2010
- Volume 7 - 2010
- Volume 6 - 2009
- Volume 5 - 2009
- Volume 4 - 2008
- Volume 3 - 2008
- Volume 2 - 2007
- Volume 1 - 2006
Commun. Comput. Phys., 30 (2021), pp. 567-601.
Published online: 2021-05
Cited by
- BibTex
- RIS
- TXT
It is reasonable to assume that a discrete convolution structure dominates the local truncation error of any numerical Caputo formula because the fractional time derivative and its discrete approximation have the same convolutional form. We suggest an error convolution structure (ECS) analysis for a class of interpolation-type approximations to the Caputo fractional derivative. Our assumptions permit the use of adaptive time steps, such as is appropriate for accurately resolving the initial singularity of the solution and also certain complex behavior away from the initial time. The ECS analysis of numerical approximations has two advantages: (i) to localize (and simplify) the analysis of the approximation error of a discrete convolution formula on general nonuniform time grids; and (ii) to reveal the error distribution information in the long-time integration via the global consistency error. The core result in this paper is an ECS bound and a global consistency analysis of the nonuniform Alikhanov approximation, which is constructed at an offset point by using linear and quadratic polynomial interpolation. Using this result, we derive a sharp $L^2$-norm error estimate of a second-order Crank-Nicolson-like scheme for linear reaction-subdiffusion problems. An example is presented to show the sharpness of our analysis.
}, issn = {1991-7120}, doi = {https://doi.org/10.4208/cicp.OA-2020-0124}, url = {http://global-sci.org/intro/article_detail/cicp/19126.html} }It is reasonable to assume that a discrete convolution structure dominates the local truncation error of any numerical Caputo formula because the fractional time derivative and its discrete approximation have the same convolutional form. We suggest an error convolution structure (ECS) analysis for a class of interpolation-type approximations to the Caputo fractional derivative. Our assumptions permit the use of adaptive time steps, such as is appropriate for accurately resolving the initial singularity of the solution and also certain complex behavior away from the initial time. The ECS analysis of numerical approximations has two advantages: (i) to localize (and simplify) the analysis of the approximation error of a discrete convolution formula on general nonuniform time grids; and (ii) to reveal the error distribution information in the long-time integration via the global consistency error. The core result in this paper is an ECS bound and a global consistency analysis of the nonuniform Alikhanov approximation, which is constructed at an offset point by using linear and quadratic polynomial interpolation. Using this result, we derive a sharp $L^2$-norm error estimate of a second-order Crank-Nicolson-like scheme for linear reaction-subdiffusion problems. An example is presented to show the sharpness of our analysis.