[论文解读] Stability analysis for the Implicit-Explicit discretization of the Cahn-Hilliard equation

Dong Li, Tao Tang|arXiv (Cornell University)|Aug 9, 2020

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

一句话总结

本文针对Cahn-Hilliard方程的隐式-显式（IMEX）时间半离散化方法，提出了一项严格的稳定性分析，未添加稳定化项或修改问题结构。在时间步长约束最优的条件下，建立了能量稳定性，该约束自然地与能量和粘性系数相关联，通过引入一种新颖的‘以能量换L∞’（TEFL）框架，解决了相场模拟中长期存在的开放性难题，从而处理了最大值原理缺失及小迁移率系数带来的刚性问题。

ABSTRACT

Implicit-Explicit methods have been widely used for the efficient numerical simulation of phase field problems such as the Cahn-Hilliard equation or thin film type equations. Due to the lack of maximum principle and stiffness caused by the effect of small dissipation coefficient, most existing theoretical analysis relies on adding additional stabilization terms, mollifying the nonlinearity or introducing auxiliary variables which implicitly either changes the structure of the problem or trades accuracy for stability in a subtle way. In this work we introduce a robust theoretical framework to analyze directly the stability of the standard implicit-explicit approach without stabilization or any other modification. We take the Cahn-Hilliard equation as a model case and prove energy stability under natural time step constraints which are optimal with respect to energy scaling. These settle several questions which have been open since the work of Chen and Shen \cite{CS98}.

研究动机与目标

• 解决长期存在的挑战：为应用于Cahn-Hilliard方程的标准IMEX格式严格证明稳定性。
• 在理论分析中消除对人工稳定化、平滑化或辅助变量的依赖。
• 推导出在小迁移率系数背景下与能量标度最优且具有物理意义的时间步长约束。
• 建立一个稳健的理论框架——TEFL，实现在自然条件下的直接稳定性分析。
• 确认IMEX方法在长时间、大规模相场模拟中的鲁棒性。

提出的方法

• 作者提出一种新颖的‘以能量换L∞’（TEFL）框架，以在不修改方程或添加稳定化项的前提下控制解的L∞范数。
• 他们分析了Chen与Shen（2002）提出的标准IMEX格式，将非线性项显式处理，线性扩散项隐式处理。
• 通过离散能量耗散论证证明能量稳定性，利用Cahn-Hilliard方程的梯度流结构。
• 关键估计涉及使用L2和H1范数来有界控制解的L∞范数，通过标度分析仔细控制时间步长。
• 时间步长约束通过平衡能量衰减与L∞增长推导得出，所得约束在能量标度上为最优。
• 分析在二维和三维周期性区域中进行，采用基于傅里叶的技术，并对时间步长进行归纳法。

实验结果

研究问题

• RQ1 能否在不添加额外稳定化项的前提下，严格证明标准IMEX格式在Cahn-Hilliard方程上的能量稳定性？
• RQ2 在小迁移率系数存在的情况下，能量稳定性的最优时间步长约束是什么？
• RQ3 在理论分析中，如何处理最大值原理缺失及小粘性系数区域的刚性问题？
• RQ4 能否在不引入辅助变量或修改非线性项的前提下进行直接分析？
• RQ5 TEFL框架是否可推广至其他具有类似刚性和非线性的相场模型？

主要发现

• 本文在与能量标度最优的时间步长约束下，建立了IMEX格式在Cahn-Hilliard方程上的条件能量稳定性。
• 时间步长约束被推导为 τ ≤ 0.0007 min{β₃⁻⁸, β₄⁻⁸ᐟ³, β₅⁻⁸}，其中 β₃, β₄, β₅ 依赖于初始能量和L²范数。
• 该分析避免了稳定化、平滑化或辅助变量的使用，提供了直接且保持结构的证明。
• TEFL框架成功地利用能量和H¹范数有界控制了解的L∞范数，使得在较宽松时间步长限制下实现稳定性成为可能。
• 所推导的时间步长约束在能量标度上几乎最优，与相分离动力学的物理预期一致。
• 该方法具有可推广性，为高阶IMEX格式及其他相场模型的精细化稳定性分析铺平了道路。