[论文解读] Characterizing the stabilization size for semi-implicit Fourier-spectral method to phase field equations

Dong Li, Zhonghua Qiao+1 more|arXiv (Cornell University)|Oct 28, 2014

Solidification and crystal growth phenomena参考文献 20被引用 7

一句话总结

本文为应用于相场方程（如Cahn-Hilliard方程和分子束外延模型）的半隐式傅里叶谱格式建立了无条件能量稳定性，消除了对非线性项的Lipschitz条件限制。证明了稳定化项的大小仅依赖于初始能量和扰动参数，而不依赖于时间步长，从而可在保持能量稳定性的同时实现大时间步长计算，并在最小假设下实现误差分析。

ABSTRACT

Recent results in the literature provide computational evidence that stabilized semi-implicit time-stepping method can efficiently simulate phase field problems involving fourth-order nonlinear dif- fusion, with typical examples like the Cahn-Hilliard equation and the thin film type equation. The up-to-date theoretical explanation of the numerical stability relies on the assumption that the deriva- tive of the nonlinear potential function satisfies a Lipschitz type condition, which in a rigorous sense, implies the boundedness of the numerical solution. In this work we remove the Lipschitz assumption on the nonlinearity and prove unconditional energy stability for the stabilized semi-implicit time-stepping methods. It is shown that the size of stabilization term depends on the initial energy and the perturba- tion parameter but is independent of the time step. The corresponding error analysis is also established under minimal nonlinearity and regularity assumptions.

研究动机与目标

• 消除现有半隐式傅里叶谱方法能量稳定性分析中对非线性势函数的限制性Lipschitz条件。
• 在最小非线性与正则性假设下，为Cahn-Hilliard方程和分子束外延模型等相场方程建立无条件能量稳定性。
• 将稳定化项的大小表征为与时间步长无关，从而在保持稳定性的同时支持大时间步长计算。
• 在最小假设下提供严格的误差分析，扩展第四阶相场问题高效数值模拟的理论基础。

提出的方法

• 为四阶相场方程提出一种带稳定化项的稳定半隐式傅里叶谱时间推进格式，其中稳定化项为O(Δt)阶。
• 采用变分形式与能量方法，证明无条件能量稳定性，且无需要求非线性势函数满足Lipschitz条件。
• 推导稳定化项大小作为初始能量与扰动参数的函数，证明其与时间步长Δt无关。
• 应用能量方法推导先验估计并建立在最小正则性与非线性假设下的误差界。
• 考虑周期域Ω = T²，并使用傅里叶谱方法实现空间谱精度离散化。
• 引入投影算子Π_N以确保高阶格式的一致性，并将分析扩展至更高阶时间推进格式。

实验结果

研究问题

• RQ1 能否在不假设非线性势函数满足Lipschitz连续性的前提下，为半隐式傅里叶谱格式证明无条件能量稳定性？
• RQ2 在该类格式中，稳定化项的大小对初始能量、扰动参数和时间步长的依赖关系如何？
• RQ3 在最小正则性与非线性假设下，如何为这些稳定化格式建立误差估计？
• RQ4 能否使稳定化项与时间步长无关，同时保持能量稳定性与精度？
• RQ5 该理论框架能否扩展至更高阶时间推进格式及更一般的相场模型？

主要发现

• 证明了在不需非线性势函数满足Lipschitz条件的前提下，稳定化半隐式傅里叶谱方法具有无条件能量稳定性。
• 稳定化项的大小仅依赖于初始能量与扰动参数，与时间步长Δt无关。
• 能量衰减单调且无条件，意味着在所导出的稳定化条件下，任意时间步长均能保持稳定性。
• 在最小非线性与正则性假设下建立了误差分析，为该格式提供了收敛性保证。
• 理论结果支持在模拟中使用大时间步长，与先前文献的计算结果一致。
• 该框架可扩展至更高阶时间推进格式及更一般的相场模型，包括六阶与向量型模型。