[论文解读] An Energy Stable Finite-Difference Scheme for Functionalized Cahn-Hilliard Equation and its Convergence Analysis

Wenqiang Feng, Zhen Guan+3 more|arXiv (Cornell University)|Oct 8, 2016

Solidification and crystal growth phenomena参考文献 6被引用 3

一句话总结

该论文提出了一种无条件能量稳定且收敛的有限差分格式，用于建模两亲系统相分离行为的六阶非线性抛物型PDE——功能化Cahn-Hilliard（FCH）方程。通过引入辅助项以凸化非凸能量泛函，作者采用凸分裂方法，确保了唯一可解性与无条件能量稳定性，并通过全局 $H_{\rm per}^{2}$ 稳定性及高效的预处理最速下降求解器，证明了对 $H^{-1}$-梯度流系统的收敛性。

ABSTRACT

We present and analyze an unconditionally energy stable and convergent finite difference scheme for the Functionalized Cahn-Hilliard equation. One key difficulty associated with the energy stability is based on the fact that one nonlinear energy functional term in the expansion appears as non-convex, non-concave. To overcome this subtle difficulty, we add two auxiliary terms to make the combined term convex, which in turns yields a convex-concave decomposition of the physical energy. As a result, an application of the convex splitting methodology assures both the unique solvability and the unconditional energy stability of the proposed numerical scheme. To deal with a 4-Laplacian solver in an $H^{-1}$ gradient flow at each time step, we apply an efficient preconditioned steepest descent algorithm to solve the corresponding nonlinear systems. In addition, a global in time $H_{ m per}^2$ stability of the numerical scheme is established at a theoretical level, which in turn ensures the full order convergence analysis of the scheme. A few numerical results are presented, which confirm the stability and accuracy of the proposed numerical scheme.

研究动机与目标

• 为两亲相分离中出现的高非线性、六阶功能化Cahn-Hilliard（FCH）方程开发一种数值稳定且收敛的格式。
• 克服能量泛函非凸性带来的挑战，该挑战使能量稳定性分析与数值实现变得复杂。
• 确保数值格式在无时间步长限制下实现无条件能量稳定与唯一可解性。
• 设计一种高效的非线性求解器，用于求解每个时间步长中出现的 $H^{-1}$-梯度流系统，以支持大规模模拟。
• 基于数值格式的全局 $H_{\rm per}^{2}$ 稳定性，建立严格的收敛性分析。

提出的方法

• 在能量泛函中引入两个辅助项，以实现凸-凹分解，从而支持凸分裂方法的应用。
• 在时间方向采用全隐式有限差分格式，在空间方向采用中心差分格式，确保无条件能量稳定性。
• 在每个时间步长中，应用预处理最速下降（PSD）算法求解由 $H^{-1}$-梯度流产生的非线性系统。
• 通过变分形式从能量泛函推导化学势 $\mu$，该泛函包含四阶拉普拉斯项与高阶非线性项。
• 建立数值解在时间上的全局 $H_{\rm per}^{2}$ 稳定性，为全阶收敛性分析提供理论基础。
• 实施周期性边界条件，并采用半隐式时间离散化，以保持质量守恒与能量耗散。

实验结果

研究问题

• RQ1 尽管能量泛函具有非凸性，是否能够构建针对FCH方程的无条件能量稳定有限差分格式？
• RQ2 如何重构非凸能量项，以支持凸分裂方法并确保唯一可解性？
• RQ3 所提格式的收敛阶是多少？能否通过稳定性估计严格证明？
• RQ4 在实际应用中，如何高效求解由 $H^{-1}$-梯度流产生的高阶非线性系统？
• RQ5 该数值格式在长时间模拟中是否能保持能量耗散与质量守恒等关键物理特性？

主要发现

• 通过辅助项实现的能量泛函凸-凹分解，使所提格式实现了无条件能量稳定与唯一可解性。
• 格式在时间方向呈现一阶收敛，空间方向观测到约2.0–2.5的收敛率，与理论预测一致。
• 预处理最速下降求解器表现出稳健的收敛性，每步迭代中残差均以近似恒定因子减少。
• 数值模拟结果证实了长时间演化过程中的能量耗散与质量守恒，符合物理预期。
• 该格式成功捕捉了复杂的形态结构，如胶束网络、双层-圆柱体连接处及Y型连接结构，与已知实验与模拟结果一致。
• 在理论层面建立了格式的全局 $H_{\rm per}^{2}$ 稳定性，从而支持全阶收敛性分析。