[论文解读] Stability and Error estimates of the SAV Fourier-spectral method for the Phase Field Crystal Equation

Xiaoli Li, Jie Shen|arXiv (Cornell University)|Jul 17, 2019

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

一句话总结

本文提出了一种用于相场晶体（PFC）方程的稳定化标量辅助变量（SAV）傅里叶谱方法，结合时间方向的SAV方法与谱空间离散化，实现了无条件能量稳定且线性的格式。该工作首次为全离散线性格式提供了严格的最优误差估计，表明时间方向二阶收敛、空间方向谱精度，其中稳定化对于在大时间步长下保持精度至关重要，尽管其并非稳定性的必要条件。

ABSTRACT

We consider fully discrete schemes based on the scalar auxiliary variable (SAV) approach and stabilized SAV approach in time and the Fourier-spectral method in space for the phase field crystal (PFC) equation. Unconditionally energy stability is established for both first- and second-order fully discrete schemes. In addition to the stability, we also provide a rigorous error estimate which shows that our second-order in time with Fourier-spectral method in space converges with order $O(Δt^2+N^{-m})$, where $Δt$, $N$ and $m$ are time step size, number of Fourier modes in each direction, and regularity index in space, respectively. We also present numerical experiments to verify our theoretical results and demonstrate the robustness and accuracy of the schemes.

研究动机与目标

• 为相场晶体（PFC）方程开发无条件能量稳定且线性的数值格式，以克服非线性格式的复杂性。
• 解决现有PFC数值方法中缺乏收敛性分析与误差估计的问题，特别是针对线性格式。
• 基于SAV框架与谱空间离散化，建立全离散格式的严格误差估计。
• 研究稳定化在保持大时间步长下精度的同时不损害稳定性的作用。
• 通过数值实验验证所提格式的鲁棒性与精度。

提出的方法

• 应用标量辅助变量（SAV）方法将PFC方程重构为线性系统，实现能量稳定性并简化时间离散化。
• 采用傅里叶谱方法进行空间离散化，以实现空间方向的谱精度，利用周期性边界条件。
• 通过引入稳定化项，提出稳定化SAV（S-SAV）格式，以提升大时间步长下的精度。
• 采用Crank-Nicolson或后向差分公式进行时间离散，构造全离散格式，确保无条件能量稳定性。
• 基于能量稳定性推导出离散解的统一有界性，从而支持严格的误差分析。
• 基于离散能量稳定性和正则性假设进行误差分析，得出最优收敛速率。

实验结果

研究问题

• RQ1 SAV方法能否扩展至PFC方程，以获得具有最优收敛速率的线性、无条件能量稳定格式？
• RQ2 S-SAV格式中的稳定化项在实现大时间步长下高精度结果中起到何种作用？
• RQ3 鉴于该问题具有高阶与非线性特性，能否为应用于PFC方程的全离散线性格式推导出严格的误差估计？
• RQ4 时间步长的选择如何影响SAV-傅里叶谱格式的精度，特别是在是否应用稳定化的情况下？
• RQ5 谱空间离散化在整体收敛速率与数值效率方面起到了多大程度的贡献？

主要发现

• 所提出的二阶SAV-傅里叶谱格式在时间与空间方向分别实现最优收敛率 $ O( au^2 + N^{-m}) $，其中 $ au $ 为时间步长，$ N $ 为傅里叶模态数。
• 已证明一阶与二阶全离散格式均具有无条件能量稳定性，确保对任意时间步长均具有鲁棒性。
• S-SAV格式中的稳定化项虽非稳定或收敛所必需，但对在大时间步长下实现合理精度至关重要。
• 数值实验验证了时间方向二阶收敛与空间方向谱收敛，证实了理论误差估计的正确性。
• 当应用稳定化时，该方法即使在大时间步长下仍能准确捕捉复杂的晶体生长动力学，包括树枝状结构、位错与界面运动。
• 该方法在模拟过冷液体中的相变与外延生长过程中表现出良好的鲁棒性与效率，且数值耗散极低。