[论文解读] Arbitrarily high-order exponential cut-off methods for preserving maximum principle of parabolic equations

Buyang Li, Jiang Yang+1 more|arXiv (Cornell University)|Oct 17, 2020

Differential Equations and Numerical Methods参考文献 33被引用 5

一句话总结

该论文提出了一种针对半线性抛物方程的高阶最大值原理保持数值方法，结合了时间方向的 $k$-阶多步指数积分格式与使用 $r$-阶多项式及高斯-洛巴托积分的集中质量有限元法。通过节点截断操作强制实现离散最大值原理，实现了 $O(\tau^k + h^r)$ 的收敛速率，且无需时间步长限制，从而在相场问题（如 Allen–Cahn 方程）中实现任意高阶精度。

ABSTRACT

A new class of high-order maximum principle preserving numerical methods is proposed for solving parabolic equations, with application to the semilinear Allen--Cahn equation. The proposed method consists of a $k$th-order multistep exponential integrator in time, and a lumped mass finite element method in space with piecewise $r$th-order polynomials and Gauss--Lobatto quadrature. At every time level, the extra values violating the maximum principle are eliminated at the finite element nodal points by a cut-off operation. The remaining values at the nodal points satisfy the maximum principle and are proved to be convergent with an error bound of $O(τ^k+h^r)$. The accuracy can be made arbitrarily high-order by choosing large $k$ and $r$. Extensive numerical results are provided to illustrate the accuracy of the proposed method and the effectiveness in capturing the pattern of phase-field problems.

研究动机与目标

• 开发适用于半线性抛物方程（特别是 Allen–Cahn 方程）的高阶数值方法，以保持最大值原理。
• 克服现有方法的局限性：要么限制时间步长，要么在保持最大值原理的同时无法实现高阶精度。
• 设计一种鲁棒的高阶格式，通过节点截断操作消除违反最大值原理的虚假值，从而保持物理一致性。
• 建立与时间步长无关的收敛速率，通过选择较大的 $k$ 和 $r$ 实现任意高阶精度。
• 展示该方法在使用大时间步长时捕捉相场模型中尖锐界面的能力。

提出的方法

• 采用 $k$-阶多步指数积分格式进行时间方向的半离散化，实现时间方向的高阶精度。
• 采用基于分段 $r$-阶多项式及高斯-洛巴托积分的集中质量有限元法进行空间半离散化。
• 在每个时间步施加节点截断操作，以消除违反最大值原理的数值解，确保离散解的边界约束。
• 截断操作仅作用于有限元节点点，保留指数积分格式的结构，同时强制实现物理边界约束。
• 避免使用稳定化项，减少人工数值耗散，提升大时间步长下的数值鲁棒性。
• 利用快速傅里叶变换（FFT）高效计算指数积分格式，提升谱效率。

实验结果

研究问题

• RQ1 能否构造一种高阶数值方法，在无需时间步长限制的条件下，保持半线性抛物方程的最大值原理？
• RQ2 在保持最大值原理的背景下，截断操作如何影响高阶指数积分格式的收敛速率？
• RQ3 所提方法能否在时间与空间方向均实现任意高阶精度，同时保持离散最大值原理？
• RQ4 与稳定化或稳定化-指数格式相比，该方法在捕捉相场问题中尖锐界面方面的表现如何？
• RQ5 截断操作对理论收敛阶的影响如何，特别是与数值实验中观察到的 $r+1$-阶空间收敛相比？

主要发现

• 所提方法在时间和空间方向均实现 $O(\tau^k + h^r)$ 的收敛速率，且对时间步长 $\tau$ 或网格尺寸 $h$ 无限制。
• 数值结果表明，该方法通过消除超出物理范围的节点值，成功维持了最大值原理，确保了解的物理一致性。
• 对于 Allen–Cahn 方程，即使在大时间步长下，该方法仍能成功捕捉圆形界面形态，而稳定化 ETD-RK2 方法在相同条件下失效。
• 在大时间步长下，该方法优于稳定化 ETD-RK2 格式，因为后者在 $\tau$ 不小时引入了显著的人工误差。
• 尽管理论上空间方向的收敛阶为 $r$，但数值实验中观察到 $r+1$-阶收敛，表明截断操作在实际中并未降低收敛速率。
• 该方法在大时间步长下表现出鲁棒性与高精度，表明去除稳定化项可增强数值稳定性和保真度。