[论文解读] Compact ADI method for solving two-dimensional Riesz space fractional diffusion equation

Sohrab Valizadeh, Alaeddin Malek+1 more|arXiv (Cornell University)|Feb 6, 2018

Fractional Differential Equations Solutions参考文献 33被引用 3

一句话总结

本文提出了一种用于求解二维瑞利空间分数阶扩散方程的紧凑交替方向隐式（ADI）方法，具有二阶空间精度和无条件稳定性。通过将紧凑有限差分格式与克拉克-尼科尔森方法及算子分裂技术相结合，该方法实现了高精度与稳定性，其有效性通过理论分析与数值实验得到验证。

ABSTRACT

In this paper, a compact alternating direction implicit (ADI) method has been developed for solving two-dimensional Riesz space fractional diffusion equation. The precision of the discretization method used in spatial directions is twice the order of the corresponding fractional derivatives. It is proved that the proposed method is unconditionally stable via the matrix analysis method and the maximum error in achieving convergence is discussed. Several numerical examples are considered aiming to demonstrate the validity and applicability of the proposed technique.

研究动机与目标

• 开发一种用于求解二维瑞利空间分数阶扩散方程的高阶、无条件稳定数值方法。
• 解决在具有瑞利导数的分数阶扩散问题中实现二阶空间精度的挑战。
• 将紧凑差分格式与交替方向隐式（ADI）方法相结合，以提高计算效率与稳定性。
• 利用矩阵分析对所提出格式的无条件稳定性与收敛阶进行理论分析。
• 通过一个包含定量误差分析的数值实例，展示该方法的有效性与适用性。

提出的方法

• 采用紧凑有限差分格式对瑞利空间分数阶导数进行离散化，实现空间上的二阶精度。
• 在时间方向应用克拉克-尼科尔森格式，以确保二阶时间精度与无条件稳定性。
• 通过AD1方法采用算子分裂技术，将二维问题分解为一系列一维问题，降低计算复杂度。
• 使用三对角矩阵 $ A_x, B_x, A_y, B_y $ 将空间离散化算子表示为矩阵形式，从而通过对角化实现高效求解。
• 利用空间算子的特征值分解，将分数阶导数矩阵表示为 $ S_x = \frac{k_t}{2}\mathcal{C}_x P \mathbf{D}_x^{\alpha/2} P^{-1} $ 和 $ T_y = \frac{k_t}{2}\mathcal{C}_y Q \mathbf{D}_y^{\beta/2} Q^{-1} $。
• 确保矩阵 $ S_x $ 与 $ T_y $ 对称且可交换，从而在每时间步高效稳定地求解线性系统。

实验结果

研究问题

• RQ1 紧凑AD1方法能否实现对二维瑞利空间分数阶扩散方程的二阶空间精度？
• RQ2 在给定的公式设定下，所提出的紧凑AD1格式是否具有无条件稳定性？
• RQ3 该数值格式在空间与时间方向上的收敛阶是多少？
• RQ4 与现有方法相比，该方法在最大误差与计算效率方面表现如何？
• RQ5 该方法能否有效处理瑞利分数阶导数的非局部性，同时保持高精度与稳定性？

主要发现

• 所提出的紧凑AD1方法实现了二阶空间精度，因为紧凑差分格式将截断误差降低至 $ \mathcal{O}(h^4) $。
• 通过矩阵分析证明了方法的无条件稳定性，其依据为空间算子的正定性与可交换性。
• 理论分析确认了时间与空间方向均为二阶收敛，最大误差被限制在 $ \mathcal{O}(k_t^2 + h_x^4 + h_y^4) $ 以内。
• 数值实验验证了方法的有效性与适用性，显示出一致的收敛速率与较低的误差水平。
• 空间算子 $ -A_x^{-1}B_x $ 与 $ -A_y^{-1}B_y $ 的特征值为正且互异，确保了稳定性与可对角化性。
• 矩阵 $ S_x $ 与 $ T_y $ 对称且可交换，从而可通过对角化在每时间步高效求解线性系统。