[論文レビュー] Minimal geodesics along volume preserving maps, through semi-discrete optimal transport

Quentin Mérigot, Jean‐Marie Mirebeau|arXiv (Cornell University)|May 13, 2015

Geometric Analysis and Curvature Flows参考文献 2被引用数 4

ひとこと要約

本稿は、半離散的最適輸送を用いて体積保存写像群における最小測地線を計算する数値的手法を導入し、非古典的で多値な解が、非圧縮性流体のオイラー方程式に対して初めて数値的に再構成可能となった。この手法はブレニエの一般化極分解を活用し、一般化流れの緩和された設定において解に収束し、フラクタル的な補集合構造と古典的流れ領域を超えた非決定的粒子軌道を明らかにした。

ABSTRACT

We introduce a numerical method for extracting minimal geodesics along the group of volume preserving maps, equipped with the L2 metric, which as observed by Arnold solve Euler's equations of inviscid incompressible fluids. The method relies on the generalized polar decomposition of Brenier, numerically implemented through semi-discrete optimal transport. It is robust enough to extract non-classical, multi-valued solutions of Euler's equations, for which the flow dimension is higher than the domain dimension, a striking and unavoidable consequence of this model. Our convergence results encompass this generalized model, and our numerical experiments illustrate it for the first time in two space dimensions.

研究の動機と目的

• 非圧縮性流体のオイラー方程式の解に対応する、$L^2$ 距離における体積保存微分同相写像群における最小測地線を数値的に計算する堅牢な数値的手法の開発。
• 次元 $d \geq 3$ における古典的最小化子の失敗を、粒子軌道が分岐・交差する可能性がある一般化流れの緩和された設定で扱う。
• 半離散的最適輸送を用いてブレニエの一般化極分解を数値実装し、古典的解が存在しない場合でも最小測地線を抽出する。
• 2次元空間においてフラクタル的補集合構造を有する非古典的で多値な流体流れの初の数値的可視化を提供する。

提案手法

• 測地線問題を、連続的経路 $\Omega = C^0([0,1], X)$ 上の確率測度上の凸最適化問題として定式化し、非圧縮性および境界条件を満たす条件下で作用 $\int_0^1 |\dot{\omega}(t)|^2 dt$ を最小化する。
• ブレニエの一般化極分解を用い、体積保存写像を勾配写像と測度保存写像の合成に分解することで、半離散的最適輸送を介した数値的計算を可能にする。
• 半離散的最適輸送フレームワークは、ターゲット測度（最終流体配置）を離散化し、ソース（初期配置）を連続的に保つことで、輸送計画の効率的計算を可能にする。
• 圧力のラグランジュ乗数としての圧力の役割を含む双対問題を解くことで一般化流れを計算し、緩和された一般化流れ設定において収束性を確立する。
• 粒子軌道は最適輸送計画からのサンプリングによって再構成され、非決定的挙動はクラスタリングおよびボックスカウント次元推定によって可視化される。
• 本手法は、https://github.com/mrgt/EulerSemidiscrete で公開されたオープンソースソフトウェアとして実装されており、再現性とさらなる実験を可能としている。

実験結果

リサーチクエスチョン

• RQ1 高次元における振動的挙動により古典的解が存在しない場合でも、体積保存写像群における最小測地線を数値的に計算可能か？
• RQ2 2次元空間において、粒子が分岐・交差する非古典的で多値な流体流れはどのように振る舞い、数値的に捉えることができるか？
• RQ3 このような一般化流れの補集合の幾何的構造はどのようなものか？理論的分析が示唆するようにフラクタル的性質を示すか？
• RQ4 圧力場の正則性は、再構成された流体軌道の非決定的性質とどのように関係するか？
• RQ5 グローバルに非決定的性質が顕在化する内部領域とは対照的に、境界付近では数値解がどの程度決定的のままで保たれるか？

主な発見

• $t_{\text{max}} < 1$ の場合、計算された解は単位正方形における古典的ベルトラミー流に一致し、決定的粒子軌道と滑らかな圧力場を示した。
• $t_{\text{max}} \geq 1$ の場合、解は非古典的かつ非決定的となり、粒子軌道が交差・分岐するようになり、一般化流れの挙動の兆候が現れた。
• $t_{\text{max}}=0.9$ の場合の解の補集合のボックスカウント次元は $D=2$ であったが、$t_{\text{max}} \in \{1.1, 1.3, 1.5\}$ の場合、2より大きく3より小さい値に増加し、フラクタル構造を示唆した。
• ドメインの境界付近では粒子運動は概ね決定的で時計回りであったが、内部では極めて非決定的かつ反時計回りの混合が顕著に観察された。
• 圧力勾配場は初期および最終時刻付近で正則性の低下の兆候を示し、理論的結果 $L^2_{\text{loc}}(]0,T[, \text{BV}(X))$ と整合的であった。
• アルゴリズムはディスク反転およびベルトラミー流の両方の一般化流れを成功裏に再構成した。後者は $t_{\text{max}} = 1.5$ でより顕著な非決定的挙動を示した。