[論文レビュー] On Euler's equation and `EPDiff'

David Mumford, Peter W. Michor|arXiv (Cornell University)|Sep 28, 2012

Nonlinear Waves and Solitons参考文献 6被引用数 3

ひとこと要約

本稿では、パラメータ $ \varepsilon, \eta \geq 0 $ を持つ正則化作用素 $ L_{\varepsilon,\eta} $ を用いて、非圧縮性で非粘性流体の Euler 方程式の解が、微分同相群上の geodesic 方程式族である EPDiff（Euler-Poincaré 方程式）の解の極限として得られることを確立している。主な結果は、流れによる運動量保存が安定性とコンパクトな台の伝播を保証するため、EPDiff 解と Euler 方程式の間の局所時間内収束バウンドが得られている。

ABSTRACT

We study a family of approximations to Euler's equation depending on two parameters $\varepsilon,η\ge 0$. When $\varepsilon=η=0$ we have Euler's equation and when both are positive we have instances of the class of integro-differential equations called EPDiff in imaging science. These are all geodesic equations on either the full diffeomorphism group $\operatorname{Diff}_{H^\infty}(\mathbb R^n)$ or, if $\varepsilon = 0$, its volume preserving subgroup. They are defined by the right invariant metric induced by the norm on vector fields given by $$ \|v\|_{\varepsilon,η} = \int_{\mathbb R^n} dx $$ where $L_{\varepsilon,η} = (I- frac{η^2}{p} riangle)^p \circ (I- frac1{\varepsilon^2} abla \circ ÷)$. All geodesic equations are locally well-posed, and the $L_{\varepsilon,η}$-equation admits solutions for all time if $η>0$ and $p\ge (n+3)/2$. We tie together solutions of all these equations by estimates which, however, are only local in time. This approach leads to a new notion of momentum which is transported by the flow and serves as a generalization of vorticity. We also discuss how delta distribution momenta lead to "vortex-solitons", also called "landmarks" in imaging science, and to new numeric approximations to fluids.

研究の動機と目的

• Euler 方程式と微分同相群上での geodesic 方程式族である EPDiff の間の厳密な関係を確立すること。
• 作用素 $ L_{\varepsilon,\eta} = (I - \frac{\eta^2}{p}\Delta)^p \circ (I - \frac{1}{\varepsilon^2}\nabla\circ\text{div}) $ が流体方程式に与える正則化効果を分析すること。
• EPDiff 方程式が $ \eta > 0 $ かつ $ p \geq (n+3)/2 $ のとき、すべての時間にわたって一意に解が存在すること（グローバルに適切に定義されていること）を示し、Euler 方程式の近似が可能になるようにすること。
• 流れに伴って輸送される一般化運動量の概念を導入し、渦度の概念を拡張すること。
• デルタ関数分布の運動量から生じるヴォーテックスソリトン（ランドマーク）の出現と、それらが数値的流体近似において果たす役割を調査すること。

提案手法

• 微分同相群 $ \text{Diff}_{H^\infty}(\mathbb{R}^n) $ 上に右不変な計量を定義し、ノルム $ \|v\|_{\varepsilon,\eta}^2 = \int \langle L_{\varepsilon,\eta}v, v \rangle \, dx $ を用いる。ここで $ L_{\varepsilon,\eta} $ は滑らかさと発散補正を組み合わせたものである。
• EPDiff 方程式 $ \partial_t m = -(v\cdot\nabla)m - \text{div}(v)m - m\cdot(Dv)^T $（ただし $ m = L_{\varepsilon,\eta}v $）を、微分同相群上の geodesic 方程式として導出する。
• 流れ $ \varphi $ に対して、運動量 1-形式 $ \widetilde{m} $ が不変であることを証明する。すなわち $ \widetilde{m}(\cdot,t) = \varphi(\cdot,t)_* \widetilde{m}(\cdot,0) $ が成り立つ。これにより、台の保存と減衰性の保存が保証される。
• EPDiff 解が $ \varepsilon, \eta > 0 $ のとき、Euler 方程式の解へ $ \varepsilon, \eta \to 0 $ の極限で局所時間内に収束することを示す推定式を確立する。
• 初期運動量を用いて流れを表現するため、Lagrangian 形式 $ \partial_t \varphi = K^{\varphi} \ast (\varphi_* \widetilde{m}(\cdot,0)) $ を用いる。ここで $ K $ は $ L_{\varepsilon,\eta} $ のグリーン関数である。
• $ \delta $-関数運動量を用いたランドマーク geodesic の数値的探索により、角運動量とエネルギーに応じて反発または吸引のダイナミクスを示す。

実験結果

リサーチクエスチョン

• RQ1 正則化パラメータ $ \varepsilon, \eta > 0 $ を持つ EPDiff 方程式の解は、$ \varepsilon, \eta \to 0 $ のとき、Euler 方程式の解へどのように収束するか？
• RQ2 一般化運動量 $ m = L_{\varepsilon,\eta}v $ は、流体の流れ下で構造をどのように保ち、渦度の概念をどのように一般化するか？
• RQ3 EPDiff 方程式がグローバル解をもつ条件は何か？ そしてこれは Euler 方程式の適切な定義（well-posedness）とどのように関係するか？
• RQ4 デルタ関数分布の運動量からヴォーテックスソリトン（ランドマーク）はどのように生じるか？ また、それらはどのような力学的挙動を示すか？
• RQ5 EPDiff フレームワークは、特に $ \eta \to 0 $ の極限において、Euler 方程式の安定な数値的近似として機能できるか？

主な発見

• Euler 方程式の解は、$ L_{\varepsilon,\eta} $ を用いた EPDiff 方程式の解の極限として得られ、$ \varepsilon, \eta \to 0 $ のとき、局所時間内に収束バウンドが確立された。
• $ \eta > 0 $ かつ $ p \geq (n+3)/2 $ のとき、EPDiff 方程式はグローバルに適切に定義されており、すべての時間にわたって一意な解が存在する。これにより、正則化された近似が可能になる。
• 運動量 1-形式 $ \widetilde{m} $ は流れ $ \varphi $ に対して不変であるため、コンパクトな台と急速な減衰性が保たれ、構造的安定性が示唆される。
• デルタ関数分布の運動量から生じるヴォーテックスソリトン（ランドマーク）は、エネルギーと角運動量に応じて散乱または捕捉のダイナミクスを示し、臨界エネルギー閾値 $ E = \frac{8}{5}|\omega|^2 $ が存在する。
• 核の微分 $ K_{0,1} $ の原点近傍におけるベクトル場は、$ (x_1,x_2) $-平面上でらせん状の流れ線を示し、$ x_3 $-軸に沿って流出する。これは局所的ヴォーテックスリングに類似した挙動を示している。
• ランドマーク geodesic の数値的統合により、背景運動量 $ \overline{m} $ を変化させることで、反発的および吸引的相互作用が生じることが示された。これは、流体不安定性のモデル化に有効である可能性を示唆している。