[論文レビュー] Radiation hydrodynamics with neutrinos: Variable Eddington factor method for core-collapse supernova simulations

Markus Rampp, Hans‐Thomas Janka|ArXiv.org|Mar 7, 2002

Neutrino Physics Research参考文献 86被引用数 213

ひとこと要約

本論文は、時間的・エネルギー的依存のボルツマン方程式を解くために可変エディングトン因子法を用いた、核崩壊超新星におけるニュートリノ輸送をシミュレートする新しい数値コードを提示する。時間陰性積分を用いたレイグリッド上でのモーメント方程式とモデルボルツマン方程式の繰り返し解法により、ニュートリノエネルギーおよびレプトン数の交換を高精度にモデル化可能であり、モジュラー設計により多次元および相対論的応用が可能である。

ABSTRACT

Neutrino transport and neutrino interactions in dense matter play a crucial role in stellar core collapse, supernova explosions and neutron star formation. Here we present a detailed description of a new numerical code for treating the time and energy dependent neutrino transport in hydrodynamical simulations of such events. The code is based on a variable Eddington factor method to deal with the integro-differential character of the Boltzmann equation. The moments of the neutrino distribution function and the energy and lepton number exchange with the stellar medium are determined by iteratively solving the zeroth and first order moment equations in combination with a model Boltzmann equation. The latter is discretized on a grid of tangent rays. The integration of the transport equations and the neutrino source terms is performed in a time-implicit way. In the present version of the program, the transport part is coupled to an explicit hydrodynamics code which follows the evolution of the stellar plasma by a finite-volume method with piecewise parabolic interpolation, using a Riemann solver for calculating the hydrodynamic states. The neutrino source terms are implemented in an operator-split step. Neutrino transport and hydrodynamics can be calculated with different spatial grids and different time steppings. The structure of the described code is modular and offers a high degree of flexibility for an application to relativistic and multi-dimensional problems at different levels of refinement and accuracy. We critically evaluate results for a number of test cases, including neutrino transport in rapidly moving stellar media and approximate relativistic core collapse, and suggest a path for generalizing the code to be used in multi-dimensional simulations of convection in neutron stars and supernovae.

研究の動機と目的

• 核崩壊超新星におけるニュートリノ輸送を高精度でシミュレートするための堅牢で柔軟な数値フレームワークの開発。
• 多次元で時間依存の流体力学的環境における積分微分ボルツマン方程式を解く計算上の課題に対処すること。
• コア崩壊およびシャクルリバーブの過程で、高密度の星間物質とのニュートリノエネルギーおよびレプトン数の交換を高精度にモデル化すること。
• 将来的な相対論的および多次元超新星シミュレーションへの応用に適したモジュラーで拡張可能なコード構造を提供すること。
• 急速な流れや近似的相対論的崩壊を含むテストケースを通じて、本手法の信頼性を検証すること。

提案手法

• ニュートリノ分布関数の方向依存性を近似するために可変エディングトン因子法を用い、完全なボルツマン方程式をモーメント方程式の集合に簡略化する。
• 零次および一次モーメント方程式を繰り返し解き、接線レイのグリッド上に離散化されたモデルボルツマン方程式と結合する。
• 剛性の高い、急激に変化する系においても安定性を確保するため、輸送および源項方程式に時間陰性積分を適用する。
• 流体力学状態のためのピecewiseパラボリック再構築とリーマンソルバを用いた明示的有限体積流体力学コードとニュートリノ輸送を結合する。
• 演算子分割を用いてニュートリノ源項を扱い、流体力学とニュートリノ輸送のための独立した空間グリッドおよび時間刻みを可能にする。
• 核種組成の遷移を3つの異なる密度・温度領域でモデル化し、解離／再結合要因 $f_{\mathrm{I}}, f_{\mathrm{II}}, f_{\mathrm{III}}$ によって連続的遷移を制御する。

実験結果

リサーチクエスチョン

• RQ1 時間的・エネルギー的依存の核崩壊超新星シミュレーションにおいて、ニュートリノ輸送の積分微分ボルツマン方程式を効率的かつ高精度に解く方法は何か？
• RQ2 繰り返しモーメント解法とレイベース離散化を併用した可変エディングトン因子法は、ニュートリノ加熱およびエネルギー沈着のモデル化に十分な精度を達成できるか？
• RQ3 コードのモジュラー設計は、将来的な多次元および相対論的シミュレーションへの拡張をどのように支援するか？
• RQ4 ニッケルの生成（$^{56}$Ni）、アルファ粒子の再結合、重元素核の解離などの核種組成遷移が、崩壊過程におけるエネルギーおよびレプトン輸送に与える影響は何か？
• RQ5 数値結果は、急速に移動する媒体内でのニュートリノ輸送や近似的相対論的崩壊といった既知のテストケースとどのように一致するか？

主な発見

• 可変エディングトン因子法は、星間物質とのエネルギーおよびレプトン数の交換を含め、核崩壊超新星におけるニュートリノ輸送の本質的物理を的確に捉えている。
• 時間陰性積分スキームにより、コアバンスやシャクル形成時に生じる極めて動的で剛性の高い環境でも数値的安定性が保証されている。
• コードのモジュラーなアーキテクチャにより、流体力学とニュートリノ輸送のための独立した空間的・時間的離散化が可能となり、計算の柔軟性と効率性が向上している。
• 急速に移動する媒体内でのニュートリノ輸送および近似的相対論的崩壊のシミュレーションは、期待される物理的挙動と良好に一致しており、本手法の正確性が検証された。
• $f_{\mathrm{I}}, f_{\mathrm{II}}, f_{\mathrm{III}}$ 要因による核種組成遷移の取り扱いは、プラズマ内での相転移に伴うエネルギー保存および熱力学的一致性を確実にしている。
• 本手法は、中性子星および超新星における対流およびニュートリノ駆動爆発機構の多次元シミュレーションへの将来的な応用に適している。