[論文レビュー] Discrete momentum maps for lattice EPDiff

Colin J. Cotter, Darryl D. Holm|arXiv (Cornell University)|Feb 14, 2006

Nonlinear Waves and Solitons参考文献 13被引用数 12

ひとこと要約

本稿は、微分同相のためのEuler-Poincaré方程式（EPDiff）の変分的粒子メッシュ（VPM）離散化を発展させ、ラグランジュ粒子上の速度場の左および右作用に対する離散運動量写像を導出する。この手法は、離散ケルビンの定理により循環を正確に保存し、数値的に検証されており、運動量保存を伴う非線形内部波の相互作用を幾何学的に一貫した枠組みでシミュレート可能である。

ABSTRACT

We focus on the spatial discretization produced by the Variational Particle-Mesh (VPM) method for a prototype fluid equation the known as the EPDiff equation}, which is short for Euler-Poincaré equation associated with the diffeomorphism group (of $\mathbb{R}^d$, or of a $d$-dimensional manifold $Ω$). The EPDiff equation admits measure valued solutions, whose dynamics are determined by the momentum maps for the left and right actions of the diffeomorphisms on embedded subspaces of $\mathbb{R}^d$. The discrete VPM analogs of those dynamics are studied here. Our main results are: (i) a variational formulation for the VPM method, expressed in terms of a constrained variational principle principle for the Lagrangian particles, whose velocities are restricted to a distribution $D_{\VPM}$ which is a finite-dimensional subspace of the Lie algebra of vector fields on $Ω$; (ii) a corresponding constrained variational principle on the fixed Eulerian grid which gives a discrete version of the Euler-Poincaré equation; and (iii) discrete versions of the momentum maps for the left and right actions of diffeomorphisms on the space of solutions.

研究の動機と目的

• 変分的粒子メッシュ（VPM）法を用いてEPDiff方程式の幾何学的に一貫した空間離散化を開発すること。
• 半離散EPDiffに対する制約付き変分原理を導出し、ラグランジュ粒子の運動とエーランジアングリッドの運動を結びつけること。
• 微分同相の粒子配置への左および右作用に対する離散運動量写像を構築すること。
• VPM離散化下でケルビンの循環則を保存する離散循環定理を確立すること。
• 今後の研究において質量密度などの運搬量を含めるようにフレームワークを拡張し、分散効果のモデル化を可能にする。

提案手法

• 領域上のベクトル場の有限次元部分空間に制限された粒子速度に制約を課したVPMのための制約付き変分原理を定式化する。
• 変分原理からエーランジアングリッド上での離散Euler-Poincaré方程式を導出し、連続EPDiff構造と整合性を保つようにする。
• 左作用運動量写像を、粒子運動量と形状関数からの標準的なVPMグリッド運動量再構成として定義する。
• グリッド速度を保存しながら運動量を再割り当てする粒子再ラベル化に対応する右作用運動量写像を導入する。
• 流れマップのヤコビアンを用いて循環ループの離散線素を構築し、ループ積分の数値近似を可能にする。
• 離散変形行列を用いて線素を発展させ、循環積分を計算し、丸め誤差の範囲での正確な保存を示す。

実験結果

リサーチクエスチョン

• RQ1 VPM法は、EPDiffの幾何的構造を保存する制約付き変分原理としてどのように定式化できるか？
• RQ2 微分同相の粒子配置への左および右作用に対する離散アナログとしての運動量写像は何か？
• RQ3 VPMフレームワーク下でケルビンの循環則を保存する離散循環定理を導出可能か？
• RQ4 離散運動量写像は、数値スキームにおける運動量保存および粒子再ラベル化の物理的解釈とどのように関係するか？
• RQ5 ヤコビアンは、離散設定における循環積分のための無限小線素の発展において果たす役割は何か？

主な発見

• VPM法はエーランジアングリッド上に離散Euler-Poincaré方程式を満たし、連続EPDiff力学と整合性を持つ。
• 離散状態における左作用運動量写像は、標準的なVPMグリッド運動量再構成式と一致する。
• 右作用運動量写像は、グリッド速度を保存しながら粒子運動量を再割り当てる離散再ラベル変換に対応する。
• 離散循環積分はシミュレーション中に正確に保存され、丸め誤差の範囲での変動のみが生じる。数値結果で確認された。
• 本手法は、流れ場における線素の伸びを的確に捉えており、時間経過に伴うループの輸送を示す図版で可視化されている。
• 本フレームワークは、質量密度やその他の運搬量を含めるように拡張可能であり、分散効果のモデル化に向けた基盤を提供する。