[論文レビュー] A new three-dimensional general-relativistic hydrodynamics code

Luca Baiotti, Ian Hawk+2 more|arXiv (Cornell University)|Apr 22, 2010

Cosmology and Gravitation Theories被引用数 15

ひとこと要約

本稿では、孤立した相対論的星の長期的かつ安定なシミュレーションを目的とした、新しい3次元一般相対論的流体力学コードWhiskyを提示する。このコードは、線形振動および非線形ダイナミクス（ブラックホール形成や安定配置への移行を含む）を的確にモデル化でき、Cactusフレームワーク内に組み込まれた高度な数値解法（PPM、ENO、HLLEリーマンソルバーなど）を用いて、複雑な時空環境下でも高い精度、収束性、耐障害性を示している。

ABSTRACT

We present a new three-dimensional general relativistic hydrodynamics code, the Whisky code. This code incorporates the expertise developed over the past years in the numerical solution of Einstein equations and of the hydrodynamics equations in a curved spacetime, and is the result of a collaboration of several European Institutes. We here discuss the ability of the code to carry out long-term accurate evolutions of the linear and nonlinear dynamics of isolated relativistic stars.

研究の動機と目的

• 強い重力場下の複雑な天体物理学的系をシミュレートできる、堅牢で精度が高くスケーラブルな3次元一般相対論的流体力学コードの開発。
• コードが球形および回転する相対論的星の線形および非線形ダイナミクス（振動や重力的崩壊を含む）を処理できるかの検証。
• 特にブラックホール形成やショック加熱を伴う時空曲率下での、長期間にわたる安定な星のモデル進化の実現。
• 今後のコア崩壊超新星や中性子星合体に伴う重力波放射の研究のための信頼できる計算ツールの提供。
• 将来的な拡張（固定メッシュ細分化など）の基盤を築くこと、特に波領域での分解能向上と波動抽出の改善を目的とする。

提案手法

• 一般相対論的流体力学方程式を、一般相対論の共形的横断的トレースレス（CTT）形式に基づいて3次元直交座標グリッド上で解く。
• 流体の不連続性を高解像度で捉えるために、パiecewise Parabolic Method（PPM）および本質的に振動しない（ENO）再構築法を用いる。
• 流体力学ソルバーにおける正確なフラックス計算のために、HLLEおよびMarquinaフラックス式のリーマンソルバーを採用。
• ロウのリーマンソルバーは、解析的左固有ベクトルとコンactなフラックス式を用いて、計算効率と安定性を向上させる。
• コードはCactus計算ツールキット内に構築されており、並列処理、モジュラー設計、アインシュタイン方程式ソルバーとの統合を可能にしている。
• 初期データは、Tolman-Oppenheimer-Volkoff（TOV）解を用い、崩壊や振動を引き起こすために微小な摂動を導入する。

実験結果

リサーチクエスチョン

• RQ1 Whiskyコードは、3次元の時空曲率下で、球対称および高速回転する相対論的星の線形的径方向および準径方向振動を正確にシミュレートできるか？
• RQ2 摂動論的手法、Cowling近似、または低次元コードと比較して、コードが振動の固有周波数をどの程度正確に再現できるか？
• RQ3 径方向モード不安定な星の非線形ダイナミクス（特にブラックホールへの重力的崩壊を含む）を安定して進化させることができるか？
• RQ4 ブラックホール形成過程における有効視界の挙動はいかなるものか？また、視界質量はどの程度正確に追跡できるか？
• RQ5 非摂動的・非線形領域において、大振幅振動とショック加熱をコードが捉えることができるか？

主な発見

• コードは、球対称および高速回転星の径方向および準径方向振動周波数を的確に再現しており、摂動論的手法および低次元数値コードと優れた一致を示している。
• 高速回転星では、摂動論的解が存在しない状況においても、固有周波数を正確に計算でき、顕著な前進を示している。
• 径方向モード不安定なモデルでは、ブラックホールへの崩壊または中心密度が低下した安定平衡配置への移行という2つの結果が得られた。
• ブラックホール形成は的確に追跡され、t = 0.246 msに有効視界が形成され、t = 0.306 msまで質量が増加したが、その後数値誤差が支配的になる。
• 理想気体状態方程式を用いたモデルでは、非線形振動がショック加熱により減衰し、エネルギー損失により最終的な静止質量が初期値よりわずかに小さくなった。
• コードは2次収束性を示し、長時間のシミュレーションにおいても高い安定性を示しており、数値的実装の堅牢性が確認された。