[論文レビュー] A General Approach for Producing Hamiltonian Numerical Schemes for Fluid Equations

Colin J. Cotter|ArXiv.org|Jan 26, 2005

Numerical methods for differential equations参考文献 34被引用数 6

ひとこと要約

本稿では、ラグランジュ粒子動力学とEuler座標系に基づくハミルトニアン評価を組み合わせることで、流体方程式のための半離散ハミルトニアンスキームを構築する一般化された粒子-メッシュ法を提示する。この手法は、正則構造を保証し、長時間にわたりO(Δt)の誤差内でハミルトニアンを保存するためのシンプレクティック時刻積分を可能にする。2次元EP-Diff方程式に対する数値実験では、優れた保存性特性が確認された。

ABSTRACT

Given a fluid equation with reduced Lagrangian $l$ which is a functional of velocity $\MM{u}$ and advected density $D$ given in Eulerian coordinates, we give a general method for semidiscretising the equations to give a canonical Hamiltonian system; this system may then be integrated in time using a symplectic integrator. The method is Lagrangian, with the variables being a set of Lagrangian particle positions and their associated momenta. The canonical equations obtained yield a discrete form of Euler-Poincaré equations for $l$ when projected onto the grid, with a new form of discrete calculus to represent the gradient and divergence operators. Practical symplectic time integrators are suggested for a large family of equations which include the shallow-water equations, the EP-Diff equations and the 3D compressible Euler equations, and we illustrate the technique by showing results from a numerical experiment for the EP-Diff equations.

研究の動機と目的

• Euler座標系における縮約ラグランジアンを用いた、流体方程式の正則ハミルトニアン半離散化を体系的に構築するためのアプローチを開発すること。
• Euler座標系における非正則なLie-Poisson構造の離散化の困難さを、ラグランジュ粒子を用いて正則形式に回復させることで克服すること。
• 粒子ベースの動力学とグリッドベースのハミルトニアン評価を組み合わせることで、流体シミュレーションにおける長時間にわたるハミルトニアンの保存を実現すること。
• 浅水域方程式、EP-Diff方程式、圧縮性Euler方程式を含む広範な流体モデルに適用可能なフレームワークを提供すること。
• 2次元EP-Diff方程式に対する数値実験を通じて、計算の実行可能性と安定性を示すこと。

提案手法

• 正則シンプレクティック構造を離散系に保証するため、ラグランジュ粒子の位置と運動量を正則変数として用いる。
• 粒子ベースの運動量と密度をEulerグリッドに補間し、グリッド点上の和としてハミルトニアンを評価する。
• 離散的微分積分法の枠組みを用いて、離散Euler-Poincaré方程式における勾配および発散作用素を表現する。
• 縮約ラグランジアンのLegendre変換を用いて離散ハミルトニアン系を導出し、粒子-メッシュ補間を介してグリッド上に射影する。
• 半離散系に対してシンプレクティック時刻積分法（例：シンプレクティックEuler法）を適用し、長時間にわたりハミルトニアンを保存する。
• 離散方程式が離散Euler-Poincaré方程式の形を満たすようにし、変分的構造を維持する。

実験結果

リサーチクエスチョン

• RQ1 Euler座標系における縮約ラグランジアンを用いた、流体方程式の正則ハミルトニアン半離散化を生成する一般化手法を開発可能か？
• RQ2 Euler座標系における流体方程式の非正則構造を、シンプレクティック時刻積分に適した正則形式に変換するにはどうすればよいか？
• RQ3 粒子-メッシュ補間は、グリッドベースのハミルトニアン評価を可能にしつつ、ハミルトニアン構造をどのように保存するのか？
• RQ4 シンプレクティック時刻積分法は、半離散系において長時間にわたりハミルトニアンをどの程度正確に保存するか？
• RQ5 本フレームワークは、浅水域-α方程式や3次元圧縮性Euler方程式のような多様な流体モデルに適用可能か？

主な発見

• 長時間にわたり離散ハミルトニアンが初期値からO(Δt)の誤差内に保たれ、シンプレクティック時刻積分の有効性が裏付けられた。
• 2次元EP-Diff方程式に対する数値実験では、発散や物理的に不自然な振動を示さない良好な解が得られた。
• H¹空間における弱解を扱えるため、δ関数のような特異的運動量分布の取り扱いが可能になった。
• 粒子-メッシュ形式は、速度場に有限要素法を用いるのを支援し、頑健な空間離散化を可能にした。
• ラグランジュ変数が複雑な微分作用素に絡まない形で、Euler形式における複雑な微分作用素の組み込みが可能になった。
• スカラー、ベクトル、テンソルを含む多様な運搬量を伴う半直接積系への応用が可能であり、フレームワークは拡張可能である。