[論文レビュー] bhlight: General Relativistic Radiation Magnetohydrodynamics with Monte Carlo Transport

Benjamin R. Ryan, Joshua C. Dolence+1 more|arXiv (Cornell University)|May 19, 2015

Astrophysical Phenomena and Observations被引用数 6

ひとこと要約

bhlight は、中程度の降着レートにおけるブラックホール降着流における放射移動を周波数依存モンテカルロ輸送でモデル化する、一般相対論的放射磁気流体力学（GRRMHD）の新しい数値スキームを提示する。閉じ込め近似を用いず、放射移動方程式を直接解くことで、光学的に薄いから中程度の厚さの範囲でコンプトン冷却およびシンクロtron冷却を正確に捉えることができ、放射が力学に支配的でないが無視できない寄与を果たす相対論的で磁化されたプラズマの高精度なシミュレーションを可能にする。

ABSTRACT

We present bhlight, a numerical scheme for solving the equations of general relativistic radiation magnetohydrodynamics (GRRMHD) using a direct Monte Carlo solution of the frequency-dependent radiative transport equation. bhlight is designed to evolve black hole accretion flows at intermediate accretion rate, in the regime between the classical radiatively efficient disk and the radiatively inefficient accretion flow (RIAF), in which global radiative effects play a sub-dominant but non-negligible role in disk dynamics. We describe the governing equations, numerical method, idiosyncrasies of our implementation, and a suite of test and convergence results. We also describe example applications to radiative Bondi accretion and to a slowly accreting Kerr black hole in axisymmetry.

研究の動機と目的

• 一般相対論的放射磁気流体力学（GRRMHD）のための数値スキームを開発し、薄いディスクと RIAF の間の中程度の降着レート領域における放射移動を正確にモデル化すること。
• Eddington やモーメント閉じ込め近似などの近似スキームに依存する既存の GRRMHD コードの限界を克服すること。これらの近似は中程度の光学的厚さで破綻する。
• 放射が支配的でないが冷却およびエネルギー輸送において重要な役割を果たすブラックホール降着流の高精度なシミュレーションを可能にすること。
• 特に M87* のような系に特に関連する、強い重力場における相対論的で磁化されたプラズマと放射の連成された力学をモデル化するフレームワークを提供すること。
• 厳密なテスト問題と収束性評価を通じて、複雑な天体物理学的環境でも数値的安定性と正確性を保証する方法の妥当性を検証すること。

提案手法

• 周波数依存放射移動方程式を曲がった時空内で直接モンテカルロ法で解く。閉じ込め近似を避ける。
• 放射は、局所的なプラズマ状態に基づいて確率的に発生・吸収・散乱されるスーパーフォトンとしてモデル化され、光線に沿って進む。
• 放射エネルギーおよび運動量は、光子分布関数から計算される放射のストレステンソルを通じて流体と結合される。
• Kerr-Schild または Kerr 度様を用いた 3+1 形式の一般相対論的フレームワークで動作し、強い重力場でのシミュレーションを可能にする。
• 流体のダイナミクスは、高分解能有限体積 MHD スキームで進化させ、放射力および加熱項はモンテカルロ光子輸送から計算される。
• 軸対称および 3D シミュレーションをサポートし、極端な状態でも精度を保つためにアダプティブメッシュ細分化と相対論的流体力学ソルバーを組み込む。

実験結果

リサーチクエスチョン

• RQ1 中程度の光学的厚さを有する一般相対論的磁気流体力学的流れにおいて、周波数依存放射移動を正確かつ効率的にシミュレートする方法は何か？
• RQ2 直接モンテカルロアプローチは、放射非効率降着流におけるコンプトン冷却およびシンクロtron冷却のモデル化において、閉じ込めに基づく手法を上回る性能を示せるか？
• RQ3 中程度の降着レート（$\dot{m} \sim 10^{-2}$）におけるブラックホール降着流の構造およびダイナミクスに、放射が及える主要な物理的効果は何か？
• RQ4 完全な放射移動を組み込むことで、標準的な RIAF モデルと比較して降着流の熱的および力学的構造にどのような影響が生じるか？
• RQ5 閉じ込め近似に依存せずに、シュバルツシルト時空における放射大気の既知の解析解をどの程度正確に再現できるか？

主な発見

• シュバルツシルト時空における放射大気の正確な解析解が再現され、流体温度が $\sqrt{-g_{00}}T(r) = T_\infty$ を満たす場合に $I_\nu / \nu^3 = B_\nu / \nu^3$ が常に成り立つことから、スキームの正確性が裏付けられた。
• 線形モード解析により、放射媒長波動モード（放射支持音波）およびコンプトン冷却による減衰が理論的期待と整合して正しく捉えられていることが示された。
• 放射的ボーンディ降着において、コードは期待される冷却 timescale および半径方向温度プロファイルを再現し、放射と物質間のエネルギー交換が正確に再現されていることを示した。
• ゆっくり降着する Kerr ブラックホールのシミュレーションでは、放射力および冷却が純粋 MHD と比較して、特にハネル領域およびコロナ領域の流れ構造に顕著な変化をもたらすことが明らかになった。
• 放射エネルギー密度および運動量フラックスといった主要な診断量について、空間的および時間的分解能で2次収束を達成した。
• 光学的厚さ領域および薄さ領域の両方で、人工的な減衰や不安定性が生じず、安定で強固な性能を示した。特に $\tau \sim 10^{-5}$ から $10^{-3}$ の領域でも、閉じ込め近似が破綻する領域で同様の性能を発揮した。