[论文解读] A fourth-order maximum principle preserving operator splitting scheme for three-dimensional fractional Allen-Cahn equations

Dongdong He, Kejia Pan+1 more|arXiv (Cornell University)|Apr 19, 2018

Solidification and crystal growth phenomena参考文献 1被引用 4

一句话总结

该论文提出了一种用于三维空间分数阶Allen-Cahn方程的四阶最大值原则保持的算子分裂格式。通过结合Strang分裂、Crank-Nicolson交替方向隐式法（ADI）处理分数阶扩散步骤、四阶有限差分格式以及Richardson外推法，该方法在时间步长约束下实现了无条件稳定性，并保持了离散最大值原则，数值实验验证了其在二维和三维测试中均具有四阶精度和最大值原则的满足性。

ABSTRACT

In this paper, by using Strang's second-order splitting method, the numerical procedure for the three-dimensional (3D) space fractional Allen-Cahn equation can be divided into three steps. The first and third steps involve an ordinary differential equation, which can be solved analytically. The intermediate step involves a 3D linear fractional diffusion equation, which is solved by the Crank-Nicolson alternating directional implicit (ADI) method. The ADI technique can convert the multidimensional problem into a series of one-dimensional problems, which greatly reduces the computational cost. A fourth-order difference scheme is adopted for discretization of the space fractional derivatives. Finally, Richardson extrapolation is exploited to increase the temporal accuracy. The proposed method is shown to be unconditionally stable by Fourier analysis. Another contribution of this paper is to show that the numerical solutions satisfy the discrete maximum principle under reasonable time step constraint. For fabricated smooth solutions, numerical results show that the proposed method is unconditionally stable and fourth-order accurate in both time and space variables. In addition, the discrete maximum principle is also numerically verified.

研究动机与目标

• 为三维空间分数阶Allen-Cahn方程开发一种高阶数值格式，以保持离散最大值原则。
• 在保持无条件稳定性的前提下，实现时间与空间方向的四阶精度。
• 确保在合理的时间步长约束下，数值解始终被限制在[−1, 1]范围内。
• 提供一种线性、高效且无条件稳定的格式，避免使用非线性求解器或在能量稳定性上进行权衡。
• 通过数值方法验证离散最大值原则的保持性以及格式的高阶收敛性。

提出的方法

• 应用Strang的二阶算子分裂方法，将分数阶Allen-Cahn方程分解为三个子步骤：其中两个子步骤涉及解析的常微分方程，另一个子步骤涉及三维线性分数阶扩散方程。
• 采用Crank-Nicolson交替方向隐式（ADI）方法求解三维分数阶扩散子步骤，将多维问题转化为一系列一维求解问题。
• 采用四阶有限差分格式对空间分数阶导数进行离散化，以保证高阶空间精度。
• 对时间积分应用Richardson外推法，将时间精度提升至四阶。
• 利用傅里叶分析证明了该整体格式对光滑解具有无条件稳定性。
• 推导出时间步长约束，以确保在数值上保持离散最大值原则。

实验结果

研究问题

• RQ1 能否为三维空间分数阶Allen-Cahn方程构造一种保持离散最大值原则的高阶（四阶）数值格式？
• RQ2 在存在分数阶导数的情况下，是否可能同时实现高阶时间与空间精度并保持无条件稳定性？
• RQ3 所提出的线性算子分裂格式是否能在合理的时间步长约束下保持解的边界限制[−1, 1]？
• RQ4 分数阶指数α如何影响相分离模拟中的界面厚度与粗化动力学？
• RQ5 Richardson外推法是否能有效提升分裂格式的时间精度，而不会损害稳定性或最大值原则的保持性？

主要发现

• 所提出的格式经傅里叶分析证明，对三维分数阶Allen-Cahn方程的光滑解具有无条件稳定性。
• 通过使用构造的光滑解进行收敛性测试，验证了该方法在时间和空间变量上均达到四阶精度。
• 离散最大值原则通过数值方法得到验证：当时间步长满足推导出的约束条件时（例如，当α=1.7时，Δt ≤ 0.4945），解被保持在[−1, 1]范围内。
• 当时间步长远超稳定性约束（例如，Δt = 2）时，解的最大值超过1，违反了离散最大值原则。
• 降低分数阶指数α会导致界面变薄且相分离粗化过程变慢，其中α = 1.2产生最不均匀且最精细的相结构。
• 该格式成功捕捉了复杂的三维相分离动力学，包括界面演化与粗化过程，具有高保真度且解保持有界。