Abstract
AbstractWe investigate a second-order accurate time-stepping scheme for solving a time-fractional diffusion equation with a Caputo derivative of order $$\alpha \in (0,1)$$
α
∈
(
0
,
1
)
. 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 $$\alpha \rightarrow 1$$
α
→
1
. Using a novel approach, we show that the proposed scheme is $$\alpha $$
α
-robust and second-order accurate in the $$L^2(L^2)$$
L
2
(
L
2
)
-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.
Funder
Australian Research Council
Publisher
Springer Science and Business Media LLC
Reference36 articles.
1. Alikhanov, A.A.: A new difference scheme for the time fractional diffusion equation. J. Comput. Phys. 280, 424–438 (2015)
2. Ashyralyev, A.: A note on fractional derivatives and fractional powers of operators. J. Math. Anal. Appl. 357, 232–236 (2009)
3. Banjai, L., Makridakis, C.G.: A posteriori error analysis for approximations of time-fractional subdiffusion problems. Math. Comput. 91, 1711–1737 (2022)
4. Bergh, J., Löfström, J.: Interpolation Spaces: An Introduction. Springer, Berlin (1976)
5. Cen, Z., Huang, J., Le, A., Xu, A.: An efficient numerical method for a time-fractional diffusion equation, preprint. arxiv:1810.07935 (2018)