[论文解读] On Linear and unconditionally energy stable Algorithms for Variable Mobility Cahn-Hilliard Type Equation with Logarithmic Flory-Huggins Potential

Xiaofeng Yang, Jia Zhao|arXiv (Cornell University)|Jan 25, 2017

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

一句话总结

本论文提出了一种针对具有对数 Flory-Huggins 势的可变迁移率 Cahn-Hilliard 方程的线性、无条件能量稳定时间积分格式，采用不变能量二次化（IEQ）方法。通过变量变换将非线性势能转化为二次形式，该方法实现了对非线性项的半显式处理，从而得到线性、对称正定系统，具有无条件能量稳定性，且可利用迭代求解器高效求解。这是首次实现该类复杂相场模型的线性、高阶、且可证明无条件能量稳定的格式。

ABSTRACT

In this paper, we consider the numerical approximations for the fourth order Cahn-Hilliard equation with concentration dependent mobility, and the logarithmic Flory-Huggins potential. One challenge in solving such a diffusive system numerically is how to develop proper temporal discretization for nonlinear terms in order to preserve the energy stability at the time-discrete level. We resolve this issue by developing a set of the first and second order time marching schemes based on a novel, called "Invariant Energy Quadratization" approach. Its novelty is that the proposed scheme is linear and symmetric positive definite because all nonlinear terms are treated semi-explicitly. We further prove all proposed schemes are unconditionally energy stable rigorously. Various 2D and 3D numerical simulations are presented to demonstrate the stability, accuracy and efficiency of the proposed schemes thereafter.

研究动机与目标

• 为广泛应用于聚合物与相分离系统的可变迁移率 Cahn-Hilliard 方程（具有对数 Flory-Huggins 势）开发高效且精确的时间积分格式。
• 克服现有格式的局限性，这些格式要么是非线性的（需要迭代求解），要么是条件稳定的（需要小时间步长），尤其在具有薄界面的刚性系统中表现不佳。
• 在每个时间步保持线性系统的同时实现无条件能量稳定性，从而实现高效且鲁棒的数值模拟。
• 将 IEQ 方法的应用扩展到复杂、有界的非线性势能（如对数 Flory-Huggins 势），这类势能在热力学一致模型中极为常见。
• 为具有复杂非线性自由能项的梯度流问题，提供一种构建高阶、线性、无条件稳定格式的通用框架。

提出的方法

• 采用不变能量二次化（IEQ）方法，通过引入新的辅助变量，将对数 Flory-Huggins 势转化为二次形式，从而实现对非线性项的半显式处理。
• 重构后的系统保持原始的能量耗散律，确保能量稳定性在离散层面得以保持。
• 利用 Crank-Nicolson 时间格式，构建一阶和二阶时间离散格式（LS1-CN 和 LS2-CN），其中非线性项通过 IEQ 重构实现半隐式处理。
• 所得线性系统为对称正定矩阵，可借助带质量凝聚预条件器的 Krylov 子空间方法高效求解。
• 在二维和三维空间离散中采用标准 Galerkin 有限元法进行求解，变分弱形式确保了格式的一致性。
• 理论分析证明，所有所提格式均具有无条件能量稳定性，且与时间步长无关。

实验结果

研究问题

• RQ1 能否为具有对数 Flory-Huggins 势的可变迁移率 Cahn-Hilliard 方程构造一种线性、高阶时间格式，实现无条件能量稳定性？
• RQ2 IEQ 方法是否能在传统凸分裂法或稳定显式方法失效或计算成本过高的情况下，提供稳定、高效且精确的数值格式？
• RQ3 基于 IEQ 的格式能否在避免非线性求解的同时，保持无条件能量稳定性，并确保对称正定线性系统？
• RQ4 所提格式在二维和三维相分离动力学的长时间模拟中表现如何，特别是针对不同初始浓度的情况？
• RQ5 IEQ 框架是否可推广至其他具有复杂、有界非线性势能的梯度流模型？

主要发现

• 所提出的 LS1-CN 和 LS2-CN 格式通过基于变分形式的能量稳定性分析，被严格证明为无条件能量稳定。
• 这些格式为线性且对称正定，可借助带质量凝聚预条件器的 Krylov 子空间求解器实现高效求解。
• 二维和三维数值模拟表明系统具有长时间稳定性，自由能随时间单调递减，证实了能量耗散特性。
• 三维相分离动力学显示：当 φ₀ = 0.3 和 φ₀ = 0.7 时形成球形液滴，当 φ₀ = 0.5 时形成双层结构，与热力学预期一致。
• 二阶格式 LS2-CN 在 δt = 0.001 时达到高精度，长时间积分中未出现不稳定或发散现象。
• 本工作首次为具有对数 Flory-Huggins 势的可变迁移率 Cahn-Hilliard 系统提出了线性、高阶且无条件能量稳定的格式，填补了该领域的重要空白。