[論文レビュー] 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法で制限的とされる補助変数の平方根形を、エネルギー関数 functional に対して任意の可逆関数に置き換えることで、各時間ステップで定数係数の線形解法が可能になりながらもエネルギー安定性を維持できる。この手法はエネルギーの下界性の制限を克服し、カップルド・キャーン・ヒルヤール・モデルのような剛性の高い系に対しても安定性を向上させる。

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バリエーションを提供することを目的とする。
• アレン＝キャーンおよびキャーン＝ヒルヤール系の相場モデル、特にカップルドBCPモデルを用いた広範な数値シミュレーションにより、手法の有効性を検証することを目的とする。

提案手法

• G-SAVフレームワークは、従来のSAV法で用いられる平方根形ではなく、エネルギー関数 functional に対して任意の可逆関数を用いた補助変数を導入する。
• 補助変数を一般化することで、変分的定式化に基づいて元の勾配流れを等価な系に変換し、エネルギー安定性を実現する。
• 3つの異なるG-SAVアプローチを提案する：(i) 補助変数を陰的に取り扱い、計算コストが最小限に抑えられる非線形代数的解法を実施；(ii) 補助変数を陽的に取り扱い、完全に陽的かつ反復の不要なスキームを実現；(iii) 修正されたエネルギーではなく、元のエネルギー散逸法則を保つアプローチ。
• これらのスキームは、試行関数と同じ関数空間に属する試験関数を用いた変分的定式化に基づく離散的エネルギー散逸解析により、無条件エネルギー安定性が証明される。
• 時間積分には2次精度のBDF法（BDF2）を用い、空間離散化は有限要素法やスペクトル法を含む任意のガラーキン型手法と互換性を持つ。
• 安定化定式化と適応的時間刻み幅制御を用いて実装し、長時間シミュレーションにおける精度と効率を向上させる。

実験結果

リサーチクエスチョン

• RQ1 補助変数を平方根関数に制限せずに、勾配流れのエネルギー安定な数値スキームを構築できるか？
• RQ2 補助変数に用いる可逆関数の選択が、数値スキームの安定性、精度、耐障害性に与える影響は何か？
• RQ3 一般化された補助変数フレームワークは、IEQやSAV法が使用する修正エネルギーではなく、元のエネルギー散逸法則を保つことができるか？
• RQ4 カップルド・キャーン＝ヒルヤール（BCP）モデルのような剛性の高い系において、G-SAVアプローチは従来のSAVおよびETDRK2スキームと比較してどのように性能を発揮するか？
• RQ5 G-SAVフレームワークは、エネルギー安定性と計算効率を維持したまま、多成分系へと拡張可能か？

主な発見

• G-SAVアプローチは、補助変数の定義に任意の可逆関数を用いることで、エネルギーが下界にない場合でも無条件エネルギー安定性を達成する。
• G-SAVに基づく数値スキームは、各時間ステップで定数係数の線形方程式のみを解けばよく、元のSAVスキームと同等の計算効率を達成する。
• 補助変数を陰的に取り扱う最初のG-SAVアプローチは、カップルド・キャーン＝ヒルヤール・モデルのような剛性の高い系において、補助変数を陽的に扱う2番目のアプローチよりもより安定である。
• 3番目のG-SAVアプローチは、IEQやSAV法が使用する修正エネルギーではなく、元のエネルギー散逸法則を保つ。これは[12]のラグランジュ乗数法と同等である。
• BCPモデルのシミュレーションでは、$ G = e^{x/c} $ および $ G = \tilde{\sqrt{x+c}} $ を用いたG-SAVスキームが、見分けがつかないような形の相分離パターンを生成し、補助関数の選択に対して本手法の高い耐障害性を確認した。
• G-SAVスキームはETDRK2スキームおよび元のSAV法と同等の精度を達成し、長時間シミュレーションにおいてより優れた安定性と効率性を示した。