[論文レビュー] A Proof of Convergence for Numerical Approximations Generated by the Locally Inertial Godunov Method in General Relativity

Zeke Vogler, Blake Temple|arXiv (Cornell University)|Dec 21, 2011

Differential Equations and Boundary Problems参考文献 9被引用数 4

ひとこと要約

本稿は、一般相対性理論におけるアインシュタイン方程式を解くための局所的慣性ガドゥノフ法の収束性を証明しており、有界で振動を示さない数値近似が弱解に収束することを示している。この手法は動的時間刻みと境界条件を伴うリーマン問題の解を用い、フラックスおよびソース項の滑らかさと全変動の有界性の仮定の下で収束が確立される。

ABSTRACT

We provide details for the proof of convergence of the locally inertial Godunov method with dynamical time dilation.

研究の動機と目的

• 局所的慣性ガドゥノフ法の数値相対性理論における収束性を厳密に証明すること。
• 有限シミュレーションにおける可変時間刻みと空間的境界を伴うガドゥノフスキームの収束性を証明する課題に対処すること。
• グロハとテンプルによるグリムスキームの議論をガドゥノフフレームワークに拡張し、ランダムサンプリングではなく平均化の手法に適応すること。
• 数値的収束性と全変動の有界性が、アインシュタイン方程式の弱解への収束を示すことを検証すること。
• 局所的慣性ガドゥノフ法を用いた相対性理論における衝撃波相互作用の数値シミュレーションの理論的基盤を提供すること。

提案手法

• 標準的なシュバルツシルト座標系におけるアインシュタイン方程式の近似解を計算するために、動的時間歪度を用いた局所的慣性ガドゥノフ法を用いる。
• 各空間セル内でリーマン問題の解を適用して系を進める。フラックスおよびソース項は、弱く同等なアインシュタイン＝ユーオール系から導出される。
• 滑らかな関数 G, f, h を用いた時間依存のODEステップを組み込み、計量 A をリプシッツ連続とみなす。
• 有限なシミュレーション領域を反映する境界条件を課し、左および右の境界データを含む。
• 各時間ステップにおける解の全変動に有界性を課すことにより、振動を制御し、安定性を確保する。
• グロハとテンプルの議論の修正版を用い、ランダムサンプリングではなくガドゥノフスキームの平均化アプローチに適応する。

実験結果

リサーチクエスチョン

• RQ1 局所的慣性ガドゥノフ法が、アインシュタイン方程式の弱解に収束する条件は何か？
• RQ2 全変動の有界性が導出されるのではなく仮定される場合、数値近似の収束性をどのように証明できるか？
• RQ3 可変時間刻みと空間的境界は、一般相対性理論におけるガドゥノフ型スキームの収束解析において果たす役割は何か？
• RQ4 項 A′·∇A f(Aij,û,x) の含め方が、リーマンセル境界における計量の不連続性をどのように扱うか？
• RQ5 数値的収束性と有界な全変動が、全変動の事前証明がなくても弱解への収束を示すことができるか？

主な発見

• 定理により、局所的慣性ガドゥノフ法で生成された近似解が有界かつ振動を示さない場合、それらはアインシュタイン方程式の弱解に収束することが保証される。
• 証明により、数値スキームの残差誤差が O(Δx) であることが示され、与えられた仮定の下で一貫性と収束性が確認される。
• セル平均値とリーマン問題解との差は、セル内での関数の全変動によって有界であることが、補題2.2で形式化されている。
• G の滑らかさと全変動の有界性を用いてODEステップの誤差を有界化し、全体のスキームで O(Δx) の誤差が得られる。
• 計量の不連続性は、ODE更新に項 A′·∇A f(Aij,û,x) を含めることで、セル境界で適切に扱われる。
• この結果は、数値的収束性と全変動の有界性の検証が容易な衝撃波相互作用のシミュレーションに特に適している。