[论文解读] A second-order accurate structure-preserving scheme for the Cahn-Hilliard equation with a dynamic boundary condition

Makoto Okumura, T. Fukao+2 more|arXiv (Cornell University)|Jul 16, 2020

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

一句话总结

本文提出了一种针对带有动态边界条件的Cahn–Hilliard方程的二阶精度、结构保持的有限差分格式，采用离散变分导数法（DVDM）。通过重构能量离散化并应用分部求和公式，该格式实现了空间二阶精度，保持了质量守恒与能量耗散特性，并确保了解的稳定性、存在性与唯一性，数值例子验证了其有效性。

ABSTRACT

We propose a structure-preserving finite difference scheme for the Cahn-Hilliard equation with a dynamic boundary condition using the discrete variational derivative method (DVDM). In this approach, it is important and essential how to discretize the energy which characterizes the equation. By modifying the conventional manner and using an appropriate summation-by-parts formula, we can use a standard central difference operator as an approximation of an outward normal derivative on the discrete boundary condition of the scheme. We show that our proposed scheme is second-order accurate in space, although the previous structure-preserving scheme by Fukao-Yoshikawa-Wada (Commun. Pure Appl. Anal. 16 (2017), 1915-1938) is first-order accurate in space. Also, we show the stability, the existence, and the uniqueness of the solution for the proposed scheme. Computation examples demonstrate the effectiveness of the proposed scheme. Especially through computation examples, we confirm that numerical solutions can be stably obtained by our proposed scheme.

研究动机与目标

• 开发一种针对带有动态边界条件的Cahn–Hilliard方程的高阶精度数值格式，以保持关键的物理结构。
• 克服先前结构保持格式在空间精度上仅达一阶的局限性。
• 通过结构保持的离散化方法，确保数值格式的稳定性、解的存在性与唯一性。
• 通过长时间模拟的数值例子验证该格式的性能，展示其稳定且符合物理规律的行为。

提出的方法

• 该格式采用离散变分导数法（DVDM）以保持连续问题的能量耗散与质量守恒特性。
• 引入一种改进的能量离散化方式，结合特定的分部求和公式，实现空间上的二阶精度。
• 通过仔细处理离散边界条件，利用标准中心差分算子近似边界上的外法向导数。
• 该格式以半隐式形式构建，其中化学势 $ P $ 通过自由能的凸 splitting 方法进行更新。
• 离散边界条件对 $ u $ 和 $ p $ 均施加齐次Neumann条件，确保与连续问题的一致性。
• 通过能量估计，证明了该格式能够保持质量守恒并实现能量耗散，且解的稳定性、存在性与唯一性得以确立。

实验结果

研究问题

• RQ1 能否针对带有动态边界条件的Cahn–Hilliard方程，设计出一种结构保持的有限差分格式，实现空间二阶精度？
• RQ2 如何设计离散边界处理方式，以在保持物理结构的同时实现更高阶精度？
• RQ3 所提出的格式在长时间模拟中是否仍能保持质量守恒与能量耗散？
• RQ4 与标准Neumann边界条件相比，动态边界条件对解的长时间行为有何影响？

主要发现

• 所提出的格式实现了空间二阶精度，相较于Fukao–Yoshikawa–Wada提出的先前结构保持格式的一阶精度有显著提升。
• 在数值实验中，质量守恒被保持在14个数量级以内，充分证实了其结构保持特性。
• 能量耗散在数值上被有效保持，离散能量泛函随时间单调递减。
• 在长时间模拟中，解保持稳定，未观察到数值发散或振荡现象。
• 数值结果表明，动态边界条件导致与标准Neumann边界条件不同的长时间行为。
• 离散能量泛函 $ A_{ m d}^{(n)} $ 被保持在9个数量级以内，进一步验证了格式的结构保持性。