[論文レビュー] A Lagrangian scheme for the incompressible Euler equation using optimal transport

Thomas Gallouët, Quentin Mérigot|HAL (Le Centre pour la Communication Scientifique Directe)|May 2, 2016

Navier-Stokes equation solutions参考文献 15被引用数 6

ひとこと要約

本稿では、測度保持微分同相群の空間におけるアーノルドの測地線的解釈と、最適輸送を用いたブレニエの外挿的アプローチを組み合わせることで、非圧縮性Euler方程式のラグランジュ数値スキームを提案する。この手法は、非圧縮性を罰則化によって強制するために半離散的最適輸送ソルバーを用い、正則解への収束を達成し、ベルトラミー流、ケルビン＝ヘルムホルツ不安定性、レイノルズ＝ターナー不安定性を含む2次元テストケースにおいても頑健な性能を示す。

ABSTRACT

We approximate the regular solutions of the incompressible Euler equation by the solution of ODEs on finite-dimensional spaces. Our approach combines Arnold's interpretation of the solution of Euler's equation for incompressible and inviscid fluids as geodesics in the space of measure-preserving diffeomorphisms, and an extrinsic approximation of the equations of geodesics due to Brenier. Using recently developed semi-discrete optimal transport solvers, this approach yields numerical scheme able to handle problems of realistic size in 2D. Our purpose in this article is to establish the convergence of these scheme towards regular solutions of the incompressible Euler equation, and to provide numerical experiments on a few simple testcases in 2D.

研究の動機と目的

• 非圧縮性Euler方程式のラグランジュ数値スキームを、最適輸送を用いて非圧縮性を保存する方法を考案すること。
• このスキームがEuler方程式の正則解に向かって理論的に収束することを確立すること。
• 現実的な粒子数を想定した2次元ベンチマーク問題において、このスキームを数値的に検証すること。
• 重力などの外部力が作用する多相流にこの手法を拡張すること。
• ケルビン＝ヘルムホルツ不安定性やレイノルズ＝ターナー不安定性のような複雑な流体不安定性を、このスキームがいかに捉えることができるかを示すこと。

提案手法

• 非圧縮性Euler方程式を測度保持微分同相群の空間における測地線として解釈する。
• 非圧縮性制約を測度保持写像の集合への距離の二乗によって罰則化する定式化を用いる。
• 最適輸送を用いて、ブレニエの極分解定理により、測度保持写像の集合への距離を計算する。
• 半離散的最適輸送ソルバーを用いてラゲール細胞とその重心を計算し、粒子の位置と速度を更新する。
• 時間積分にはシンプレクティックEulerスキームを用い、速度更新にはラゲール細胞の変位の重心を含む。
• 速度更新式に力項を追加することで、外部力（例：重力）を含める拡張が可能となる。

実験結果

リサーチクエスチョン

• RQ1 最適輸送に基づくラグランジュスキームは、非圧縮性Euler方程式の正則解に向かって収束するか？
• RQ2 ベルトラミー流や不安定性を含む標準的な2次元流体力学テストケースにおいて、このスキームはどの程度の性能を示すか？
• RQ3 重力のような外部力が導入された場合、このスキームは安定的かつ正確に保たれるか？
• RQ4 粒子数が大規模（例：300,000粒子）であっても、幾何学的・物理的整合性を維持できるか？
• RQ5 粒子の分割方法（規則的 vs. 最適化済み）の選択が、シミュレーションの精度および力学的挙動にどのように影響するか？

主な発見

• 適切な仮定の下で、数値スキームが非圧縮性Euler方程式の正則解に収束することを理論的に確立した。
• N = 900粒子を用いた単位正方形上のベルトラミー流では、安定的で構造的なラゲール細胞を示し、正しく再現された。
• ケルビン＝ヘルムホルツ不安定性のシナリオでは、最適化済みパーティションを用いたスキームが対称性を破り、物理的期待と整合する渦を形成した。N = 300,000粒子を用いた。
• 重力が作用する50,000粒子を用いたレイノルズ＝ターナー不安定性実験では、界面不安定性が形成され、多相流におけるスキームの頑健性が確認された。
• レイノルズ＝ターナー系の計算は、1つのラップトップコアで6時間未満で完了し、計算効率の高さが示された。