[论文解读] An analysis of the L1 Scheme for the subdiffusion equation with nonsmooth data

Bangti Jin, Raytcho Lazarov+1 more|arXiv (Cornell University)|Jan 1, 2015

Fractional Differential Equations Solutions参考文献 24被引用 57

一句话总结

本文重新研究了带有阶数 $\alpha \in (0,1)$ 的Caputo分数阶导数的亚扩散方程的L1格式，证明了在光滑与非光滑初值下均具有 $O(\tau)$ 的收敛速率——即使解在时间方向不具有 $C^2$ 正则性也成立。该分析进一步推广至一般扇形算子及时空分数阶扩散方程，数值实验验证了误差估计的精确性与鲁棒性。

ABSTRACT

The subdiffusion equation with a Caputo fractional derivative of order $α\in(0,1)$ in time arises in a wide variety of practical applications, and it is often adopted to model anomalous subdiffusion processes in heterogeneous media. The L1 scheme is one of the most popular and successful numerical methods for discretizing the Caputo fractional derivative in time. The scheme was analyzed earlier independently by Lin and Xu (2007) and Sun and Wu (2006), and an $O(τ^{2-α})$ convergence rate was established, under the assumption that the solution is twice continuously differentiable in time. However, in view of the smoothing property of the subdiffusion equation, this regularity condition is restrictive, since it does not hold even for the homogeneous problem with a smooth initial data. In this work, we revisit the error analysis of the scheme, and establish an $O(τ)$ convergence rate for both smooth and nonsmooth initial data. The analysis is valid for more general sectorial operators. In particular, the L1 scheme is applied to one-dimensional space-time fractional diffusion equations, which involves also a Riemann-Liouville derivative of order $β\in(3/2,2)$ in space, and error estimates are provided for the fully discrete scheme. Numerical experiments are provided to verify the sharpness of the error estimates, and robustness of the scheme with respect to data regularity.

研究动机与目标

• 为解决现有L1格式收敛性分析中假设时间方向具有 $C^2$ 正则性的局限性，该假设在非光滑初值下不成立。
• 在初值数据正则性假设最小化的前提下，建立L1格式的严格 $O(\tau)$ 收敛速率。
• 将误差分析推广至一般扇形算子，包括时空分数阶扩散方程中出现的算子。
• 通过数值实验验证理论结果，展示误差估计的精确性与鲁棒性。

提出的方法

• 利用时间导数的分段线性逼近，重新表述Caputo分数阶导数的L1格式。
• 基于拉普拉斯变换技术，采用生成函数方法分析局部截断误差。
• 利用扇形算子的性质及分数阶扩散半群的光滑效应，推导误差界。
• 将分析方法应用于具有Riemann-Liouville空间导数的时空分数阶扩散方程的全离散有限元或有限差分格式。
• 采用能量方法及 $L^2$-范数下的稳定性估计，控制时间与空间方向的误差。
• 通过在均匀与非均匀时间网格上的数值实验，验证理论收敛速率。

实验结果

研究问题

• RQ1 当经典 $C^2$ 正则性假设在非光滑初值下失效时，L1格式能否在亚扩散问题中实现一阶收敛？
• RQ2 当解在 $t=0$ 处出现奇点时，L1格式的最优收敛速率是什么？
• RQ3 误差行为如何依赖于初值的正则性，特别是在 $t \to 0$ 的极限下？
• RQ4 该误差分析能否推广至一般扇形算子，包括时空分数阶扩散方程中的算子？
• RQ5 在时空分数阶模型中，$O(\tau)$ 收敛速率对 $\alpha$ 与 $\beta$ 的不同取值是否均具有鲁棒性？

主要发现

• L1格式在 $L^2$-范数下对光滑与非光滑初值均实现 $O(\tau)$ 的收敛速率，即使解在时间方向不具有 $C^2$ 正则性也成立。
• 该收敛速率对 $\alpha \in (0,1)$ 具有鲁棒性，数值实验在 $t=0.1$、$t=0.01$ 与 $t=0.001$ 处均确认了一阶收敛。
• 对于 $\beta \in (3/2,2)$ 的时空分数阶扩散方程，全离散格式在适当正则性条件下达到 $O(\tau + h^2)$ 的误差估计。
• 非光滑初值下的误差在 $t \to 0$ 时表现为 $t^{\alpha}$，证实了定理4.1中理论估计的精确性。
• 数值结果表明，无论 $\alpha$、$\beta$ 或初值类型如何，收敛速率始终近似为1.0，显示出良好的鲁棒性。
• 理论分析适用于一般扇形算子，将L1格式的应用范围从标准椭圆算子拓展至更广泛的情形。