[論文レビュー] On the energy stability of Strang-splitting for Cahn-Hilliard

Dong Li, Chaoyu Quan|arXiv (Cornell University)|Jul 12, 2021

Solidification and crystal growth phenomena参考文献 27被引用数 5

ひとこと要約

本稿は、Cahn-Hilliard方程式に適用された2次精度Strang分割スキームについて、最初の均一エネルギー安定性および高次のSobolevノルムの保存性を確立する。著者らは、新しい理論的枠組みを導入し、小さな時間刻み条件下で、すべての時間ステップにおいて数値解が $ H^k $ ノルムで有界であることを証明した。これにより、長年の未解決問題であった、この方程式に対する作用素分割法におけるエネルギー安定性に関する問題が解決された。

ABSTRACT

We consider a Strang-type second order operator-splitting discretization for the Cahn-Hilliard equation. We introduce a new theoretical framework and prove uniform energy stability of the numerical solution and persistence of all higher Sobolev norms. This is the first strong stability result for second order operator-splitting methods for the Cahn-Hilliard equation. In particular we settle several long-standing open issues in the work of Cheng, Kurganov, Qu and Tang \cite{Tang15}.

研究の動機と目的

• Cahn-Hilliard方程式に対する2次精度作用素分割法のエネルギー安定性に関する長年の未解決問題を解消すること。
• 小さな時間刻み条件下で、すべての時間ステップにおける数値解の $ H^k $ ノルムにおける一様有界性を確立すること。
• 非線形および線形部分問題を含むStrang分割スキームの安定性を解析するための厳密な理論的枠組みを提供すること。
• 非線形時間発展ステップが $ L^p $ 收縮性を欠くにもかかわらず、数値スキームにおいて高次のSobolevノルムの保存およびエネルギー減衰を証明すること。

提案手法

• Cahn-Hilliard方程式のStrang分割のエネルギー安定性を解析するための新しい理論的枠組みを導入する。
• Strang分割スキームを $ S_L( au/2)S_N( au)S_N(\tau/2) $ として定義する。ここで $ S_L $ は線形部分問題を、$ S_N $ は非線形部分問題を解く。
• 滑らかさ推定とエネルギーに基づく議論を用いて、長時間にわたり $ H^{40} $ ノルムの数値解を制御する。
• エネルギーがしきい値を超える最初の時間ステップを仮定して、矛盾による一様エネルギー境界を証明する。
• 生成子のテイラー展開による一貫性解析を用い、数値スキームと正確なPDE解との差を $ O(\tau^3) $ まで比較する。
• 初期時間ステップにおける一様 $ H^1 $ 制御を確立し、滑らかさおよび補間を用いて高階ノルムへ拡張する。

実験結果

リサーチクエスチョン

• RQ1 Cahn-Hilliard方程式に適用された2次精度Strang分割スキームについて、均一エネルギー安定性を証明できるか？
• RQ2 小さな時間刻み条件下で、長時間にわたり数値解が高次のSobolevノルムを保存するか？
• RQ3 非線形時間発展において $ L^p $ 收縮性が欠如する場合でも、安定性解析でこれを克服できるか？
• RQ4 時間刻み $ \tau $ に依存しない、数値解に対する事前境界を確立することは可能か？
• RQ5 数値解のエネルギーはどのように変化するか？また、このスキームにおいてエネルギーが非増加であることを示せるか？

主な発見

• すべての $ k \geq 1 $ に対して、数値解は $ H^k $ ノルムで一様有界であり、その境界は初期データおよびパrameter $ \nu $, $ \gamma_0 $ のみに依存する。
• スキームは均一エネルギー安定性を示す：$ \sup_{n \geq 1} E(u^n) \leq G_{\nu,\gamma_0} $ であり、$ G $ は $ \nu $ および初期 $ H^1 $ ノルムにのみ依存する。
• 著者らは $ \|u^n\|_{H^{40}} \leq C_{\nu,G} $ を $ n $ に関して一様に成立させることを証明し、高正則性の保存を確立した。
• エネルギー減衰が確認された：エネルギーがしきい値を超えると仮定すれば、$ E(u^{n_0+1}) < E(u^{n_0}) $ が成り立ち、その逆が成り立つと矛盾が生じる。
• スキームはPDEと $ O(\tau^3) $ まで一貫しており、一貫性誤差はテイラー展開および滑らかさ推定を用いて制御された。
• この結果により、Chengら[25]が提起した、Strang分割によるCahn-Hilliard方程式の安定性および収束に関する未解決問題が解決された。