[论文解读] Revisit of Semi-Implicit Schemes for Phase-Field Equations

Tao Tang|arXiv (Cornell University)|Jun 12, 2020

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

一句话总结

本文重新探討了具有通用勢函數的 Allen-Cahn 方程的半隱式格式，證明在時間步長限制下，離散最大值原理與能量穩定性均成立。該格式對拉普拉斯項採用隱式處理，對非線性項採用顯式處理，確保解有界且能量單調衰減，無需凸分裂或穩定化項。

ABSTRACT

It is a very common practice to use semi-implicit schemes in various computations, which treat selected linear terms implicitly and the nonlinear terms explicitly. For phase-field equations, the principal elliptic operator is treated implicitly to reduce the associated stability constraints while the nonlinear terms are still treated explicitly to avoid the expensive process of solving nonlinear equations at each time step. However, very few recent numerical analysis is relevant to semi-implicit schemes, while "stabilized" schemes have become very popular. In this work, we will consider semi-implicit schemes for the Allen-Cahn equation with {\em general potential} function. It will be demonstrated that the maximum principle is valid and the energy stability also holds for the numerical solutions. This paper extends the result of Tang \& Yang (J. Comput. Math., 34(5):471--481, 2016) which studies the semi-implicit scheme for the Allen-Cahn equation with {\em polynomial potentials}.

研究动机与目标

• 建立相場模型中半隱式格式的嚴謹穩定性保障，特別是針對具有一般非線性勢函數的 Allen-Cahn 方程。
• 證明在不使用凸分裂或穩定化項的情況下，最大值原理與能量穩定性在半隱式格式中成立。
• 將先前針對多項式勢函數的結果推廣至一般光滑勢函數，提升現有穩定性框架的通用性。
• 為半隱式格式在長時間模擬中的有效性提供理論基礎，其中穩定性與有界性至關重要。

提出的方法

• 採用半隱式時間離散化，其中拉普拉斯項以隱式方式處理，非線性項以顯式方式處理，並在每個時間步長求解相應的線性系統。
• 應用單調格式論證與數學歸納法，在涉及勢函數導數的時間步長條件下證明離散最大值原理。
• 運用離散分部積分與泰勒展開分析能量演化，並建立能量穩定性。
• 根據勢函數導數的最大值推導時間步長限制，以確保能量衰減。
• 考慮週期性邊界條件，並使用中心差分法對拉普拉斯項進行有限差分空間離散化。
• 依賴由最大值原理保證的解有界性，以控制能量演化。

实验结果

研究问题

• RQ1 對於具有通用勢函數的 Allen-Cahn 方程，其半隱式格式是否滿足離散最大值原理？
• RQ2 是否可以在不使用凸分裂或穩定化項的情況下，證明半隱式格式的能量穩定性？
• RQ3 基於勢函數導數的時間步長限制是否足夠確保最大值原理與能量穩定性同時成立？
• RQ4 針對多項式勢函數的穩定性結果能否推廣至一般光滑勢函數？
• RQ5 在半隱式格式中，能量衰減與解有界性的必要與充分條件是什麼？

主要发现

• 若時間步長滿足 Δt ≤ 1 / max|f′(u)|（在解的範圍內），則半隱式格式的離散最大值原理成立，確保 γ₋ ≤ φⱼⁿ ≤ γ₊。
• 能量穩定性已得證：在相同時間步長條件下，對所有 n 有 Eₕ(φⁿ⁺¹) ≤ Eₕ(φⁿ)。
• 證明依賴於由最大值原理保證的解有界性，以及離散分部積分與泰勒展開的運用。
• 能量衰減通過涉及勢函數二階導數的界來建立，時間步長限制確保能量非增。
• 該結果將先前針對多項式勢函數的工作推廣至一般光滑勢函數，擴展了半隱式格式的適用範圍。
• 該框架提示，類似能量穩定性可能也適用於具有對數自由能的 Cahn-Hilliard 方程，但需進一步深入分析。