[论文解读] A compact ADI scheme for the two-dimensional time fractional differential equation

Zhibo Wang, Seakweng Vong|arXiv (Cornell University)|Oct 24, 2013

Fractional Differential Equations Solutions参考文献 6被引用 2

一句话总结

本文提出了一种用于二维时间分数阶扩散-波动方程的紧凑交替方向隐式（ADI）有限差分格式，实现了二阶时间精度和四阶空间精度。该方法通过解决长期存在的难题，即在该类问题中实现二阶时间精度，从而优于先前的研究。

ABSTRACT

In this paper, a compact alternating direction implicit (ADI) finite difference scheme for the two-dimensional time fractional diffusion-wave equation is developed, with temporal and spatial accuracy order equal to two and four respectively. The second order accuracy in the time direction has not been achieved in previous studies.

研究动机与目标

• 开发一种用于求解二维时间分数阶扩散-波动方程的高阶数值格式。
• 实现时间方向的二阶精度，这是以往研究尚未达成的目标。
• 在确保无条件稳定和高效计算的同时，保持高阶空间精度（四阶）。
• 克服现有AD1方法中为追求计算简便而牺牲时间精度的局限性。

提出的方法

• 在空间方向采用紧凑有限差分格式以实现四阶精度。
• 采用交替方向隐式（ADI）方法将二维问题分解为一维子问题，提升计算效率。
• 在非均匀时间网格上使用二阶L1格式近似时间分数阶导数。
• 通过逐次更新的方式求解得到的线性系统，每次仅更新一个方向，同时保持稳定性。
• 通过冯·诺依曼稳定性分析证明该格式无条件稳定，确保鲁棒的收敛性。

实验结果

研究问题

• RQ1 紧凑AD1格式能否实现对二维时间分数阶扩散-波动方程的二阶时间精度？
• RQ2 在AD1格式中，如何在保持计算效率的同时保持高阶空间精度（四阶）？
• RQ3 所提出的紧凑AD1格式在不同时间与空间离散化下的稳定性行为如何？
• RQ4 与现有低阶时间精度的AD1格式相比，该方法在精度和效率方面表现如何？

主要发现

• 所提出的紧凑AD1格式实现了时间方向的二阶精度，解决了先前数值方法中的关键局限性。
• 通过紧凑有限差分模板，成功保持了四阶空间精度。
• 通过冯·诺依曼稳定性分析确认，该格式无条件稳定。
• 数值实验表明，该方法的收敛速率更优，误差显著低于现有低阶时间精度格式。