[论文解读] The operator-splitting method for Cahn-Hilliard is stable

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

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

一句话总结

本文首次为在二维周期环面上应用的Cahn-Hilliard方程的算子分裂方法建立了无条件能量稳定性。通过结合半群估计、Sobolev嵌入及归纳法，作者证明了在时间步长 $ \tau < \tau_* $ 的小时间步长条件下，$ H^{k_0} $ 范数的统一有界性与能量衰减，确保了在初始数据假设最小条件下的长期稳定性和数值格式的收敛性。

ABSTRACT

We prove energy stability of a standard operator-splitting method for the Cahn-Hilliard equation. We establish uniform bound of Sobolev norms of the numerical solution and convergence of the splitting approximation. This is the first unconditional energy stability result for the operator-splitting method for the Cahn-Hilliard equation. Our analysis can be extended to many other models.

研究动机与目标

• 为应用于Cahn-Hilliard方程的算子分裂方法建立首个无条件能量稳定性结果。
• 在初始数据正则性假设最小的条件下，证明数值解在所有时间步长下于 $ H^{k_0} $ 范数下的统一有界性。
• 通过建立稳定性与收敛性的通用框架，将分析推广至其他梯度流模型。
• 解决算子分裂格式在Cahn-Hilliard方程背景下能量稳定性的开放问题。

提出的方法

• 该方法采用标准的Strang分裂：首先通过半群 $ S_L(\tau) = e^{-\tau \nu \Delta^2} $ 求解线性扩散部分 $ \partial_t u = -\nu \Delta^2 u $，然后通过 $ S_N(\tau)w = w + \tau \Delta(w^3 - w) $ 求解非线性部分 $ \partial_t u = \Delta(u^3 - u) $。
• 数值格式定义为 $ u^{n+1} = S_L(\tau) S_N(\tau) u^n $，时间步长 $ \tau > 0 $，并在二维周期环面 $ \mathbb{T}^2 $ 上进行分析。
• 关键估计依赖于半群 $ e^{-\tau \nu \Delta^2} $ 的光滑性质，包括 $ L^p $ 有界性以及通过谱乘子实现的高频模态衰减。
• 证明通过归纳法与Bootstrap方法控制解的 $ L^\infty $、$ H^1 $ 及更高阶Sobolev范数，利用嵌入定理与能量型不等式。
• 一个关键步骤涉及对离散能量 $ E_1(u^{n+1}) $ 的控制，其定义为 $ \||\nabla|^{-1}(e^{\tau \nu \Delta^2} - 1)^{1/2} u^{n+1}\|_{L^2}^2 $ 与势能 $ \int ((u^{n+1})^2 - 1)^2 dx $ 的组合。
• 通过递归展开格式并应用热半群的光滑估计，将分析扩展至 $ H^{k_0} $ 正则性。

实验结果

研究问题

• RQ1 能否在不施加对数值解的限制性假设下，证明Cahn-Hilliard方程算子分裂方法的无条件能量稳定性？
• RQ2 在初始数据属于 $ H^{k_0} $ 的条件下，数值解是否对所有时间步长在 $ H^{k_0} $ 范数下保持统一有界？
• RQ3 保证分裂格式稳定性和收敛性的最小时间步长 $ \tau_* $ 是多少？
• RQ4 该框架能否推广至具有类似结构的其他梯度流模型？
• RQ5 在分裂格式下，数值解的 $ L^\infty $ 与 $ H^1 $ 范数行为如何？

主要发现

• 本文证明了在时间步长限制下，Cahn-Hilliard方程算子分裂方法的无条件能量稳定性，即对所有 $ n \geq 1 $，有 $ E_1(u^{n+1}) \leq E_1(u^n) $。
• 对于初始数据属于 $ H^{k_0}(\mathbb{T}^2) $（$ k_0 \geq 2 $）的情形，建立了统一有界性：$ \sup_{n \geq 0} \|u^n\|_{H^{k_0}} \leq B_1 < \infty $，其中 $ B_1 $ 依赖于 $ \|u^0\|_{H^{k_0}} $、$ \nu $ 与 $ k_0 $。
• 时间步长限制为 $ \tau < \tau_* = c \cdot \min\{ \alpha^{-8}, \alpha^{-8/3} \} \nu^3 $，其中 $ \alpha $ 依赖于初始数据与 $ \nu $，$ c > 0 $ 为一很小的绝对常数。
• 导出了一个 $ L^\infty $ 有界性：$ \sup_n \|u^n\|_{\infty} \leq \alpha (\nu \tau)^{-1/8} + \alpha \tau (\nu \tau)^{-7/8} $，表明在小 $ \tau $ 条件下解的增长受到控制。
• 证明了格式在 $ H^1 $ 与 $ L^\infty $ 范数下的收敛性，收敛速率由统一有界性与证明中的误差估计所隐含。
• 该分析具有鲁棒性，并可推广至其他具有类似梯度流结构的模型，特别是涉及 $ \Delta^2 $ 与双阱势项的模型。