Free-surface Stokes problem: stability estimates and time-step improvements

技術要約：自由表面ストークス問題：安定性評価とタイムステップの改善

問題設定
氷床力学、マントル対流、溶岩流などの高粘性自由表面流の正確なシミュレーションには、結合されたストークス／自由表面方程式を解くことが必要である。これらの問題は、流体速度と圧力が自由表面の高さの進化と結合した移動領域によって特徴付けられる。空間における有限要素法と時間におけるオイラー法を用いた標準的な離散化は一般的であるが、これらはしばしば深刻な安定性の制約を受ける。具体的には、陽解法および半陰解法によるオイラー法は、安定性を維持するために極めて小さなタイムステップを必要とすることが多く、計算コストを増大させる。既存の安定化技術（自由表面安定化アルゴリズム：FSSAなど）はタイムステップのサイズを改善するものの、安定性に関する先験的な理論的境界を欠いており、無条件安定性や体積保存を保証するものではない。

手法
著者らは、ニュートン流体および非ニュートン流体（べき乗則粘性）の両方について、結合されたストークス／自由表面系の安定性と保存特性に関する厳密な理論解析を提示している。解析は以下の2つのレベルで行われる：

1. 連続レベル： 著者らは、ストークス運動方程式のエネルギーと自由表面高さのエネルギーを結合する全エネルギーバランスを導出した。この導出は、自由表面のアドベクション項と表面積分との間の特定の関係を確立し、その後、それらをストークス方程式の重力項に結びつけることに依存している。
2. 離散レベル： 本研究では、空間における有限要素法と、陽解法、半陰解法、および陰解法による時間ステップを用いた完全離散スキームを検討している。
   • 正則化： ストークス問題の適題性を実現するために必要な自由表面高さ関数の固有の正則性の欠如に対処するため、著者らは、時間に関して陰的に扱われる高次エッジ安定化項 $J(s^{n+1}, w)$ を導入している。
   • 安定化の提案： 陽解法オイラー法のエネルギー推定における「偽の」項の特定に基づき、著者らはストークスの弱形式に追加される新しい安定化項を提案している。この項は、自由表面における法線方向速度にペナルティを課すものであり、自由表面方程式の陽的な取り扱いに起因する不安定化項を正確に打ち消すように設計されている。

主な貢献

• 全エネルギーバランス： 本論文は、運動エネルギーと表面高さのエネルギーを結合した、結合ストークス／自由表面系の新しい全エネルギーバランスを導出した。これは、これまで文献に欠けていた、安定性を理解するための基礎的な構成要素として特定された。
• 理論的安定性証明：
  • 陰解法オイラー法： （自由表面高さがリプシッツ連続であると仮定した場合）、タイムステップサイズに対して無条件に安定であり、かつ領域の体積を保存することを証明した。
  • 陽解法オイラー法（未安定化）： 法線方向の表面速度に関連する偽の項により、いかなるタイムステップサイズにおいても全エネルギーノルムにおいて不安定であることを証明した。
  • 半陰解法オイラー法（FSSA）： 条件付きの安定性しかなく、決定的な点として、領域の体積を保存しないことを証明した。
  • 提案された安定化陽解法オイラー法： 著者らは、提案する安定化項が、リプシッツ連続性の仮定の下で、陽解法オイラー法をタイムステップサイズに対して無条件に安定にし、かつ体積を保存することを証明した。
• 体積保存： 本論文は、提案された安定化陽解法スキームが領域の体積を保存する一方で、FSSA安定化半陰解法スキームは保存しないことを示す理論的証明と数値的検証を提供している。

結果
数値実験は、ニュートン流体および非ニュートン流体の両方のケースを含め、グリーンランド氷床のデータに基づく複雑な幾何形状を用いて、オープンソースライブラリ FEniCS を用いて実施された。

• 安定性： 未安定化の陽解法オイラー法は、テストされたすべてのタイムステップで不安定であった。FSSA安定化スキーム（半陰および陽の両方）は条件付きの安定性を示し、大きなタイムステップでは失敗した。対照的に、提案された EE-Stabilized スキームは、ニュートンおよび非ニュートン両方のケースにおいて、任意に大きなタイムステップでも安定性を維持した。
• 体積保存： EE-Stabilized スキームは、機械精度まで領域の体積を保存した。SIE-FSSA スキームは体積を保存できず、理論的知見を裏付けた。
• 精度： 収束テストにより、提案された安定化項が離散化された問題と整合していることが示された。EE-Stabilized スキームは、メッシュとタイムステップの同時細分化の下で、表面高さと速度について一次の精度を維持し、その誤差は他の安定化スキームと同等またはそれ以下であった。
• 振動： 陰的なエッジ安定化項は、追加のタイムステップ制限を導入することなく、自由表面高さの高周波振動を効果的に抑制した。

意義と主張
本論文は、提案された安定化陽解法オイラー法が、計算コストの高い陰解法の代替として、堅牢かつ効率的な選択肢となることを主張している。タイムステップサイズに対して無条件の安定性を理論的に保証し、体積保存を確保することで、本手法は複雑な幾何形状やソース項を持つ自由表面流の安定したシミュレーションを可能にする。著者らは、彼らの安定化が、この結合問題の陽解法オイラー離散化に対して任意のタイムステップサイズを理論的に保証する初めてのものであることを強調している。このアプローチは、対称的であり（反復解法に有利）、実装が単純で、調整可能なパラメータを持たないことが特筆される。著者らは、解析において表面高さのリプシッツ連続性を仮定しているものの、導出された $L^2$ ノルム安定性推定値は、適切な正則化が適用された際の適題性のための必要な基礎であると結論付けている。今後の課題として、適題性のための正則化項の探索、およびこの安定性解析をナビエ・ストークス／自由表面問題や高次時間積分器へ拡張することが示唆されている。