[論文レビュー] General Relativistic Hydrodynamics with a Roe solver

Frits Eulderink, Garrelt Mellema|arXiv (Cornell University)|Nov 15, 1994

Pulsars and Gravitational Waves Research参考文献 1被引用数 15

ひとこと要約

本稿では、外部重力場下での任意座標系における一般相対論的流体動力学へのロウ型近似リーマンスolverを提案する。この手法は、古典的流体動力学から相対論的流体動力学へとロウ法を拡張したものであり、超相対論的衝撃波（0.99999cに達する）を的確に捉えることができ、従来の明示的スキームを上回る性能を示す。また、保存則に基づく特徴的数値フレームワークを提供し、特殊相対性理論および非デカルト座標系に応用可能である。

ABSTRACT

We present a numerical method to solve the equations of general relativistic hydrodynamics in a given external gravitational field. The method is based on a generalization of Roe's approximate Riemann solver for the non relativistic Euler equations in Cartesian coordinates. The new method is applied to a set of standard test problems for general relativistic hydrodynamics, and is shown to perform well in comparison to existing numerical schemes. In contrast to existing explicit methods the present method can cope with strong relativistic shocks. By-products are: the characteristic form of the general relativistic Euler equations, a numerical method for special relativity that can deal with strong discontinuities, a numerical scheme for the integration of the Euler equations in an arbitrary coordinate system, possibly under the influence of (external) gravity, and a novel method to incorporate source terms in numerical schemes.

研究の動機と目的

• 外部重力場下での任意座標系における一般相対論的オイラー方程式を解くための堅牢な数値法の開発。
• 光速に近い超相対論的衝撃波を扱えない既存の明示的スキームの限界を克服すること。
• 古典的流体動力学で効果的であることが実証済みのロウの近似リーマンスローバーを相対論的領域に一般化し、任意の調整パrameterを用いずに物理的一致性を保つこと。
• 特殊相対性理論、任意のフレームにおける古典的流体動力学、およびソース項を含む一般相対論的流れに適用可能な統一的フレームワークの提供。
• 重力や曲がった時空内での流れに対して、物理的に整合的かつ正確にソース項を数値スキームに組み込む方法の確立。

提案手法

• 任意座標系における一般相対論的オイラー方程式の保存形に、ロウ平均線形化されたフラックスヤコビアンを適用し、ロウの近似リーマンスローバーを適応化する。
• 一次元における相対論的オイラー方程式の特徴的形を導出し、左・右状態の差が小さい極限での正確な波動伝播解析を可能にする。
• 時間分割法を用いたソース項の新規処理を導入し、ソース項は2次精度のストラング分割法で統合する。
• 保存則に基づく有限体積離散化を用い、数値フラックスはロウスローバーで計算することで、強い衝撃領域における保存性と安定性を確保する。
• 特殊相対性理論および一般相対性理論の両極限に適用：ミンコフスキー計量を用いることで特殊相対論的流体動力学用スキームが得られ、非相対論的極限を取ることで任意フレームにおける古典的流体動力学が回復される。
• 標準的なテスト問題（相対論的衝撃波、希薄波、ステーショナリーショックチューブなど）を用いた検証により、高い精度と頑健性を示す。

実験結果

リサーチクエスチョン

• RQ1 ロウの近似リーマンスローバーは、一般相対論的流体動力学方程式へと成功裏に一般化可能であり、強力な衝撃波に対してもその頑健性を保つのか？
• RQ2 本スキームは、0.99999cに達する超相対論的流れを扱う際、既存の手法と比較してどの程度優れているのか？
• RQ3 本手法は一般相対性理論を越えて、特殊相対論的流体動力学や非デカルト座標系における古典的流体動力学にどの程度応用可能なのか？
• RQ4 重力や時空曲率効果といったソース項は、どのように保存則に基づく数値スキームに正確に組み込むことができるのか？
• RQ5 左・右状態の差が小さい極限において、相対論的方程式の特徴的形がロウスローバーから自然に導かれるのか？これは物理的一致性の裏付けとなる。

主な発見

• 提案されたロウスローバーは、ローレンツ因子が約224（0.99999cに対応）に達する超相対論的衝撃波を的確に捉えることができ、これまでは失敗していた領域をカバーする。
• 空間的・時間的2次精度を維持し、任意の人工的粘性を用いずに保存則に基づく近似リーマンスローバーで数値フラックスを計算する。
• 一般相対論的オイラー方程式の特徴的形が導出され、左・右状態の差が小さい極限で回復されることを示し、スキームの物理的一致性を裏付ける。
• 計量をミンコフスキーに設定することで特殊相対論的流体動力学に適用可能であり、非相対論的極限を取ることで任意フレームにおける古典的流体動力学に回復される。広範な一般性を示す。
• ソース項は2次精度のストラング分割法で組み込まれ、曲がった時空内での重力などの外部力の正確な取り扱いが可能である。
• ショックチューブや希薄波といった標準的な相対論的流体動力学ベンチマーク問題に対する数値的テストにより、正確な解と優れた一致を示し、スキームの頑健性と精度を確認した。