[论文解读] The generalized scalar auxiliary variable approach (G-SAV) for gradient flows

Qing Cheng|arXiv (Cornell University)|Feb 1, 2020

Solidification and crystal growth phenomena参考文献 35被引用 7

一句话总结

本文提出了广义标量辅助变量（G-SAV）方法，这是一种用于构建梯度流无条件能量稳定数值格式的新型框架。通过用能量泛函的任意可逆函数替代传统SAV方法中受限的平方根形式辅助变量，G-SAV在每时间步实现仅含常系数的线性求解，同时保持能量稳定性。该方法克服了能量有界性限制，显著提升了刚性系统（如耦合Cahn-Hilliard模型）的鲁棒性。

ABSTRACT

We establish a general framework for developing, efficient energy stable numerical schemes for gradient flows and develop three classes of generalized scalar auxiliary variable approaches (G-SAV). Numerical schemes based on the G-SAV approaches are as efficient as the original SAV schemes \cite{SXY19,cheng2018multiple} for gradient flows, i.e., only require solving linear equations with constant coefficients at each time step, can be unconditionally energy stable. But G-SAV approaches remove the definition restriction that auxiliary variables can only be square root function. The definition form of auxiliary variable is applicable to any reversible function for G-SAV approaches . Ample numerical results for phase field models are presented to validate the effectiveness and accuracy of the proposed G-SAV numerical schemes.

研究动机与目标

• 开发一种梯度流无条件能量稳定时间离散化的通用框架，以克服现有SAV和IEQ方法中辅助变量受限于平方根形式的局限。
• 实现自由能不从下方有界的系统中能量稳定格式，这是传统SAV方法难以克服的限制。
• 通过在每个时间步仅需求解常系数线性方程，保持计算效率，与原始SAV格式相当。
• 提供多种G-SAV变体，以在鲁棒性、实现简便性与保持原始能量耗散律之间取得平衡。
• 通过Allen-Cahn和Cahn-Hilliard相场模型（包括耦合BCP模型）的大量数值模拟验证该方法。

提出的方法

• G-SAV框架引入一个辅助变量，定义为能量泛函的任意可逆函数，替代经典SAV方法中使用的平方根形式。
• 通过引入广义辅助变量，将原始梯度流转化为等价系统，基于变分形式实现能量稳定性。
• 提出三种不同的G-SAV方法：(i) 对辅助变量采用隐式处理，每步仅需代价可忽略的非线性代数求解；(ii) 对辅助变量采用显式处理，实现完全显式、非迭代格式；(iii) 一种保持原始能量耗散律而非修改后能量的变体。
• 所提格式具有无条件能量稳定性，其证明基于与试函数属于同一函数空间的变分形式的离散能量耗散分析。
• 时间离散采用二阶BDF格式（BDF2），空间离散与任意Galerkin型方法兼容，如有限元法或谱方法。
• 采用稳定化格式与自适应时间步长策略，以提升长时间模拟中的精度与效率。

实验结果

研究问题

• RQ1 能否在不将辅助变量限制为平方根函数的前提下，构造梯度流的能量稳定数值格式？
• RQ2 用于辅助变量的可逆函数选择如何影响数值格式的稳定性、精度与鲁棒性？
• RQ3 广义辅助变量框架能否保持原始能量耗散律，而非如IEQ与SAV方法中使用的修改后能量？
• RQ4 与经典SAV和ETDRK2格式相比，G-SAV方法在刚性系统（如耦合Cahn-Hilliard（BCP）模型）上的表现如何？
• RQ5 G-SAV框架能否扩展至多组分系统，同时保持能量稳定性和计算效率？

主要发现

• G-SAV方法通过使用任意可逆函数定义辅助变量，实现了梯度流的无条件能量稳定性，消除了对能量从下方有界的依赖。
• 基于G-SAV的数值格式在每个时间步仅需求解常系数线性方程，其计算效率与原始SAV格式相当。
• 第一种G-SAV方法（辅助变量隐式处理）在如耦合Cahn-Hilliard模型等刚性系统中，比第二种方法（显式辅助变量）更具鲁棒性。
• 第三种G-SAV方法保持了原始能量耗散律，而不同于IEQ与SAV方法中使用的修改后能量，且与文献[12]中的拉格朗日乘子法等价。
• BCP模型的模拟表明，采用$ G = e^{x/c} $和$ G = \tilde{\sqrt{x+c}} $的G-SAV格式产生几乎无法区分的形态演化模式，证实了该方法对辅助函数选择的鲁棒性。
• G-SAV格式在精度上可与ETDRK2格式及原始SAV方法相媲美，且在长时间模拟中表现出更优的稳定性与效率。