[论文解读] Efficient schemes with unconditionally energy stability for the anisotropic Cahn-Hilliard Equation using the stabilized-Scalar Augmented Variable (S-SAV) approach

Xiaofeng Yang|arXiv (Cornell University)|Apr 8, 2018

Solidification and crystal growth phenomena参考文献 12被引用 10

一句话总结

本论文提出两种针对各向异性Cahn-Hilliard方程的二阶、无条件能量稳定的全隐时间格式，采用稳定化标量增广变量（S-SAV）方法。通过将SAV方法与三个线性稳定项结合，该方法以半显式方式处理所有非线性项，并在每时间步仅求解三个解耦的常系数线性方程，从而在保持大时间步长的同时消除强各向异性引起的振荡。该方法是首个针对此模型实现二阶精度与高效率的、可严格证明无条件能量稳定的方案。

ABSTRACT

In this paper, we consider numerical approximations for the anisotropic Cahn-Hilliard equation. The main challenge of constructing numerical schemes with unconditional energy stabilities for this model is how to design proper temporal discretizations for the nonlinear terms with the strong anisotropy. We propose two, second order time marching schemes by combining the recently developed SAV approach with the linear stabilization approach, where three linear stabilization terms are added. These terms are shown to be crucial to remove the oscillations caused by the anisotropic coefficients, numerically. The novelty of the proposed schemes is that all nonlinear terms can be treated semi-explicitly, and one only needs to solve three decoupled linear equations with constant coefficients at each time step. We further prove the unconditional energy stabilities rigorously, and present various 2D and 3D numerical simulations to demonstrate the stability and accuracy.

研究动机与目标

• 为各向异性Cahn-Hilliard方程开发高效且无条件能量稳定的数值格式，该问题因强各向异性引起的数值振荡而具有挑战性。
• 通过结合SAV方法与线性稳定项，克服现有格式的局限性——如缺乏可证明的能量稳定性或计算成本过高。
• 通过在每个时间步仅求解线性常系数方程，确保时间方向二阶精度的同时保持计算效率。
• 证明所提格式即使在强各向异性条件下，也能允许大时间步长而不产生不稳定现象。
• 为所提格式提供无条件能量稳定性的严格数学证明，这是对先前方法的关键性改进。

提出的方法

• 采用SAV（标量增广变量）方法，通过引入新的标量辅助变量，将能量泛函重写为二次型形式，从而实现对非线性项的半显式处理。
• 引入三个线性稳定项（S1、S2、S3），以抑制由各向异性系数γ(∇φ/|∇φ|)引起的高频振荡，尤其在各向异性较强时效果显著。
• 构建两种二阶BDF2型时间离散格式（Schemes BDF2-S 和 BDF2-W），无条件保持能量耗散律。
• 在每个时间步将系统解耦为三个独立的常系数线性椭圆方程，可通过标准线性求解器求解。
• 采用与试函数空间相同的变分形式，推导弱形式，并通过能量估计确保稳定性。
• 引入Willmore正则化项（βG(φ)）以改善数值条件并减少振荡，尤其在三维模拟中效果显著。

实验结果

研究问题

• RQ1 能否为各向异性Cahn-Hilliard方程构造一种二阶时间格式，即使在强各向异性条件下也能保证无条件能量稳定性？
• RQ2 如何高效且稳定地以半显式方式处理涉及各向异性系数γ(∇φ/|∇φ|)的非线性项？
• RQ3 额外的线性稳定项在抑制由各向异性系数引起的数值振荡中起到何种作用？
• RQ4 所提格式是否能在强各向异性条件下保持高精度与稳定性，同时支持大时间步长？
• RQ5 是否可能在每个时间步仅求解线性常系数方程的前提下，实现可证明的无条件能量稳定性？

主要发现

• 所提S-SAV格式实现了无条件能量稳定性，其严格性通过基于变分形式的能量估计得到证明。
• 数值实验表明，S1 = S2 = 2 且 S3 = 0.001 的组合能有效抑制高频振荡；而 S1 = S2 = 0 且 S3 = 0.001 则导致非物理振荡。
• 该格式允许使用大时间步长（δt = 1e−4）以保持精度，δt = 1e−2、1e−3 和 1e−4 的结果在相场分布上几乎完全一致。
• 增大各向异性参数α（从0.1增至0.3）会使平衡态形貌的金字塔角更加尖锐，与物理预期一致。
• 减小Willmore正则化参数β可使平衡形状的角更加尖锐，且能量随β减小而更快下降，时间演化曲线已证实此现象。
• 该方法具有通用性，可推广至其他具有复杂非线性的梯度流问题，因其证明框架依赖于具有相容Galerkin近似的变分形式。