We 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)$. 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\to 1$. Using a novel approach, we show that the proposed scheme is $\alpha$-robust and second-order accurate in the $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.
翻译:本文研究了一种用于求解含阶数$\alpha \in (0,1)$的Caputo导数的分数阶扩散方程的二阶精度时间步进格式。该格式的基本思想基于局部积分后线性插值。在经典扩散情形下(即当$\alpha\to 1$时),该格式退化为标准Crank-Nicolson格式。通过一种新颖的方法,我们证明了在适当的时间分级网格下,该格式在$L^2(L^2)$范数下具有α-鲁棒性和二阶精度。为完整性考虑,我们采用Galerkin有限元方法进行空间离散,并在给定数据的合理正则性假设下讨论了误差分析。最后给出了若干数值结果。