[论文解读] Finite difference schemes for multi-term time-fractional mixed diffusion-wave equations

Zhaopeng Hao, Guang Lin|arXiv (Cornell University)|Jul 24, 2016

Fractional Differential Equations Solutions参考文献 12被引用 7

一句话总结

本文提出了一种针对具有阶数在 (1,2) 范围内的 Caputo 导数的多术语时间分数阶混合扩散-波动方程的一阶精度 L2 型有限差分格式。通过结合改进的 L2 近似与紧凑差分格式，并采用能量法，作者证明了该方法在时间上具有第一阶精度、在空间上具有第四阶精度的无条件收敛性，数值例子验证了该方法的有效性，并设计了一种时间复杂度为线对数级别的快速求解器。

ABSTRACT

The multi-term time-fractional mixed diffusion-wave equations (TFMDWEs) are considered and the numerical method with its error analysis is presented in this paper. First, a $L2$ approximation is proved with first order accuracy to the Caputo fractional derivative of order $β\in (1,2).$ Then the approximation is applied to solve a one-dimensional TFMDWE and an implicit, compact difference scheme is constructed. Next, a rigorous error analysis of the proposed scheme is carried out by employing the energy method, and it is proved to be convergent with first order accuracy in time and fourth order in space, respectively. In addition, some results for the distributed order and two-dimensional extensions are also reported in this work. Subsequently, a practical fast solver with linearithmic complexity is provided with partial diagonalization technique. Finally, several numerical examples are given to demonstrate the accuracy and efficiency of proposed schemes.

研究动机与目标

• 为阶数在 (1,2) 范围内的多术语时间分数阶混合扩散-波动方程（TFMDWEs）开发一种稳健的数值方法，尽管这些方程在模拟异常扩散方面具有重要意义，但目前研究仍显不足。
• 通过构建稳定且精确的有限差分格式并辅以严格的误差分析，解决分数阶导数带来的非局部性与内存需求挑战。
• 将该方法扩展至分布阶方程和二维问题，确保其广泛适用性。
• 设计一种基于部分对角化的快速求解器，以线对数时间复杂度降低计算成本。
• 利用能量法提供完整的误差分析，建立无条件收敛性与稳定性。

提出的方法

• 为阶数 β ∈ (1,2) 的 Caputo 分数阶导数推导出一种改进的 L2 近似，实现时间方向的一阶精度。
• 将改进的 L2 近似与低阶导数的 L1 近似相结合，形成完全隐式的紧凑有限差分格式。
• 将能量法应用于离散系统，利用一种新颖的离散分数阶 Poincaré-Friedrichs-Sobolev 不等式，建立稳定性与收敛性。
• 基于部分对角化方法开发了一种快速求解器，以处理线性系统的 Toeplitz 结构，实现 O(N log N) 的时间复杂度。
• 将该方法扩展至阶数在 (0,2) 范围内的分布阶方程以及二维问题，同时保持高阶精度。
• 通过数值实验验证了理论收敛速率，并展示了在 1D 和 2D 问题上的高效性。

实验结果

研究问题

• RQ1 能否为阶数在 (1,2) 范围内的 Caputo 分数阶导数严格推导出一种一阶精度的 L2 型近似？
• RQ2 能否构建一种紧凑有限差分格式，并证明其在多术语时间分数阶混合扩散-波动方程上的无条件收敛性？
• RQ3 该格式在时间和空间上的收敛速率是多少？能否利用能量法予以证明？
• RQ4 能否为紧凑格式所产生的线性系统设计一种时间复杂度为线对数级别的快速求解器？
• RQ5 该格式在分布阶方程和二维问题上的表现如何？

主要发现

• 证明了改进的 L2 近似对阶数 β ∈ (1,2) 的 Caputo 导数具有第一阶精度，并给出了严格的误差界。
• 利用能量法证明了紧凑有限差分格式具有无条件收敛性，时间方向为第一阶精度，空间方向为第四阶精度。
• 建立了离散分数阶 Poincaré-Friedrichs-Sobolev 不等式，该不等式在收敛性证明中起到了关键作用。
• 数值结果验证了 1D 和 2D 问题在时间和空间方向上的理论收敛速率。
• 成功实现了时间复杂度为 O(N log N) 的快速求解器，与标准求解器相比显著降低了计算成本。
• 该格式在阶数位于 (0,2) 范围内的分布阶方程中也表现出稳健性能，证实了其广泛适用性。