[论文解读] On the energy stability of Strang-splitting for Cahn-Hilliard

Dong Li, Chaoyu Quan|arXiv (Cornell University)|Jul 12, 2021

Solidification and crystal growth phenomena参考文献 27被引用 5

一句话总结

本文首次为应用于Cahn-Hilliard方程的二阶Strang分裂格式建立了统一的能量稳定性以及高阶Sobolev范数的持久性。通过引入一种新颖的理论框架，作者证明了在小时间步长条件下，数值解在所有时间步上均保持在$ H^k $范数下的有界性，从而解决了该方程在算子分裂方法中长期存在的开放性问题。

ABSTRACT

We consider a Strang-type second order operator-splitting discretization for the Cahn-Hilliard equation. We introduce a new theoretical framework and prove uniform energy stability of the numerical solution and persistence of all higher Sobolev norms. This is the first strong stability result for second order operator-splitting methods for the Cahn-Hilliard equation. In particular we settle several long-standing open issues in the work of Cheng, Kurganov, Qu and Tang \cite{Tang15}.

研究动机与目标

• 解决关于Cahn-Hilliard方程二阶算子分裂方法能量稳定性的长期开放问题。
• 在小时间步长条件下，建立数值解在所有时间步上关于$ H^k $范数的统一有界性。
• 为分析具有非线性和线性子步骤的Strang分裂格式中的稳定性，提供严格的理论框架。
• 证明尽管非线性演化步骤缺乏$ L^p $收缩性，数值格式中仍能保持高阶Sobolev范数的持久性与能量衰减。

提出的方法

• 引入一种新理论框架，以分析Cahn-Hilliard方程中Strang分裂的能量稳定性。
• 将Strang分裂格式定义为$ S_L( au/2)S_N( au)S_N( au/2) $，其中$ S_L $求解线性子问题，$ S_N $求解非线性子问题。
• 利用光滑性估计与基于能量的论证方法，控制长时间区间内数值解的$ H^{40} $范数。
• 通过反证法证明统一能量界，假设存在首个能量超过阈值的时间步。
• 通过传播子的泰勒展开进行一致性分析，将数值格式与精确PDE解在$ O( au^3) $范围内进行比较。
• 建立早期时间步的统一$ H^1 $控制，并利用光滑性与插值技术将其推广至高阶范数。

实验结果

研究问题

• RQ1 能否为应用于Cahn-Hilliard方程的二阶Strang分裂格式证明统一能量稳定性？
• RQ2 在小时间步长条件下，数值解是否能在长时间区间内保持高阶Sobolev范数？
• RQ3 在稳定性分析中，能否克服非线性演化中缺乏$ L^p $收缩性的问题？
• RQ4 能否建立与时间步长$ \tau $无关的数值解先验界？
• RQ5 数值解的能量如何演化？能否证明其在该格式下非增？

主要发现

• 数值解在所有$ k \geq 1 $下均保持在$ H^k $范数下的统一有界性，界仅依赖于初始数据及参数$ \nu $、$ \gamma_0 $。
• 该格式表现出统一能量稳定性：$ \sup_{n \geq 1} E(u^n) \leq G_{\nu,\gamma_0} $，其中$ G $仅依赖于$ \nu $与初始$ H^1 $-范数。
• 作者证明了$ \|u^n\|_{H^{40}} \leq C_{\nu,G} $对所有$ n $一致成立，从而确立了高正则性的持久性。
• 能量衰减得到证实：若能量超过阈值，则$ E(u^{n_0+1}) < E(u^{n_0}) $，若违反则导致矛盾。
• 该格式与PDE在$ O(\tau^3) $范围内一致，且通过泰勒展开与光滑性估计控制了一致性误差。
• 该结果解决了Cheng等人[25]关于Cahn-Hilliard方程Strang分裂稳定性与收敛性的开放问题。