[論文レビュー] Selective decay for the rotating shallow-water equations with a structure-preserving discretization

Rüdiger Brecht, Werner Bauer+3 more|arXiv (Cornell University)|Jul 5, 2021

Numerical methods for differential equations参考文献 41被引用数 5

ひとこと要約

本稿では、エネルギーを保存しつつ、離散的カシミール散逸枠組みを用いてポテンシャル・エンストロフィーを選択的に散逸させる、回転流し浅い水方程式の構造保存型変分積分法を提案する。変分離散化と連続的選択的減衰理論を組み合わせることで、エネルギー損失なしに長期的安定性を実現し、物理的忠実性を保ちながら小スケール構造を解像する点で、標準的なバイハーモニック散乱と比較して優れている。

ABSTRACT

International audience

研究の動機と目的

• 長期的安定性を確保しつつ物理的忠実性を維持する数値的ウェザーモデルおよび気候モデルの課題に対処する。
• 双曲型保存則ソルバーにおけるエネルギー保存と不自然な小スケールノイズの抑制の間にある矛盾を克服する。
• 連続的選択的減衰を模倣する離散的枠組みを構築し、エネルギーを保存しながらポテンシャル・エンストロフィーを散逸させる。
• 平面上および球面上の浅い水方程式の長期的統合に適した変分積分法を安定化する。
• エネルギー保存性が劣化する標準的散乱法（例：バイハーモニック粘性）の物理的に一貫した代替案を提供する。

提案手法

• 離散的変分原理に基づく離散的オイラー＝ラグランジュ方程式を出発とする変分積分法フレームワークを採用する。
• 連続的リー＝ポワソン枠組みを模倣した離散的選択的減衰メカニズムを実装し、特にポテンシャル・エンストロフィーの散逸を目的とする。
• パヴロフらの構造保存型離散化法を用いて、離散レベルでのミメティック性および保存則の保持を保証する。
• カシミール（ポテンシャル・エンストロフィー）の離散的汎関数微分を定義し、運動量方程式における標的的散逸を可能にする。
• f-平面および球面浅い水方程式に本手法を適用し、エッジベースの有限要素風離散微分作用素を用いる。
• 得られた半離散系を、元のスケームと同一の次数でエネルギーを保存する時間積分法で統合する。

実験結果

リサーチクエスチョン

• RQ1 浅い水方程式においてエネルギーを保存しつつポテンシャル・エンストロフィーを散逸させる離散的選択的減衰枠組みを構築できるか？
• RQ2 カシミール散逸を施した本手法の変分積分法は、標準的なバイハーモニック散乱と比較してエネルギー保存性および解の品質においてどのように異なるか？
• RQ3 本手法は、人工的なエネルギー損失を導入せずに、平面上および球面上の長期的シミュレーションをどの程度安定化できるか？
• RQ4 離散的交換子構造を維持でき、連続理論と整合性を持つと期待できるか？
• RQ5 複雑なテストケース（例：山地上を流れる流れ）において、従来の散乱法と比較して小スケール構造をよりよく解像できるか？

主な発見

• 提案手法は、散乱なしの変分積分法と同等のオーダーでエネルギーを保存するが、バイハーモニック散乱スキームでは顕著なエネルギー損失を示す。
• カシミール散逸を施したシミュレーションでは、わずかで一貫したエネルギー増加が観測され、エネルギー保存構造への最小限の乱れであることが示された。
• エンストロフィー散逸スキームは、不自然な小スケールノイズを効果的に抑制し、球面上で100日間の統合においても解の品質を維持した。
• エネルギーおよびエンストロフィースペクトルの結果から、カシミール散逸スキームはバイハーモニック散乱法よりも小スケール構造をよりよく解像していることが明らかになった。
• 山地上を流れる流れのテストケースでは、カシミール散逸を施したシミュレーションが高分解能で計算された基準解とよく一致したが、散乱なしのスケームはノイズが多くなっていた。
• 離散的交換子構造は連続理論と整合的であり、離散的選択的減衰枠組みの理論的基盤が妥当であることが検証された。