[論文レビュー] A Mathematical Description of the IDSA for Supernova Neutrino transport, its discretization and a comparison with a finite volume scheme for Boltzmann's Equation

Heiko Berninger, Emmanuel Frénod+4 more|Archive ouverte UNIGE (University of Geneva)|Nov 29, 2012

Numerical methods for differential equations参考文献 7被引用数 4

ひとこと要約

本稿は、核崩壊超新星におけるニュートリノ輸送のための等方的拡散源近似（IDSA）の厳密な数学的定式化を提供し、時間分割と有限体積法を用いて解かれた簡略化されたボルツマン方程式と比較している。IDSAは、低分散誤差にもかかわらず、吸収・放出と自由拡散の間の遷移領域に主に局在する、相対$L^2$ノルムで約20％の顕著なモデル化誤差を生じることを示している。

ABSTRACT

In this paper we give an introduction to the Boltzmann equation for neutrino transport used in core collapse supernova models as well as a detailed mathematical description of the \emph{Isotropic Diffusion Source Approximation} (IDSA). Furthermore, we present a numerical treatment of a reduced Boltzmann model problem based on time splitting and finite volumes and revise the discretization of the IDSA for this problem. Discretization error studies carried out on the reduced Boltzmann model problem and on the IDSA show that the errors are of order one in both cases. By a numerical example, a detailed comparison of the reduced model and the IDSA is carried out and interpreted. For this example the IDSA modeling error with respect to the reduced Boltzmann model is numerically determined and localized.

研究の動機と目的

• 核崩壊超新星におけるニュートリノ輸送のための等方的拡散源近似（IDSA）の数学的に厳密な導出および記述を提供すること。
• 静的背景物質を仮定した簡略化されたボルツマンモデルに基づく、時間分割と有限体積法を用いた数値スキームの開発。
• IDSAの離散化を、一貫性のある空間的および角度的グリッドを保ったまま、簡略化モデルに再表現・適応すること。
• 数値実験を用いて、簡略化されたボルツマンモデルに対するIDSAのモデル化誤差を定量化および局在化すること。
• 特に吸収・放出と自由拡散の挙動の遷移領域において、IDSAが重要な物理的状態を捉える有効性と限界を評価すること。

提案手法

• 球対称性を仮定したニュートリノ輸送の$Ø(v/c)$ボルツマン方程式を導出する。エネルギー、半径、角度依存性を含む。
• 分布関数を吸収・放出成分と自由拡散成分に分解し、拡散に基づく源項で結合するIDSAを導入する。
• 反応項と輸送項を時間分割で分離し、反応項は解析的に、輸送項は保存的な有限体積法で解く。
• 先行研究におけるIDSAの離散化を、簡略化モデルに適応させ、空間的および角度的グリッドの一貫性を保つ。
• 分散誤差が1％未満になるよう、高解像度グリッド（$N_r=257$, $N_\mu=129$, $\Delta t=1.5625\cdot10^{-3}$）を用いる。
• IDSAとボルツマン解との間の$L^2$相対差をモデル化誤差として計算し、空間的に絶対差をマップすることで誤差の局在化を図る。

Figure 1: Flow map of the flux of neutrinos due to transport and without reactions. We can see that incoming neutrinos ( $\mu<0$ ) turn into outgoing neutrinos ( $\mu>0$ ) and, therefore, eventually leave the computational domain again.

実験結果

リサーチクエスチョン

• RQ1 IDSAモデルは、簡略化されたボルツマン方程式と比較して、解の精度においてどのように定量的に異なるか？
• RQ2 IDSAのモデル化誤差は、空間的にどこに顕著に現れ、物理的状態とどのように関連しているか？
• RQ3 IDSAのモデル化誤差の大きさは、簡略化されたボルツマンモデルに対してどの程度で、遷移領域に依存するか？
• RQ4 両モデルにおける離散化誤差は、比較の信頼性にどのように影響するか？
• RQ5 なぜIDSAは反応支配領域においてボルツマンモデルよりも速やかに平衡状態に収束するのか？

主な発見

• 簡略化されたボルツマンモデルに対するIDSAのモデル化誤差は、$L^2$ノルムで約20％であり、誤差範囲は±2％である。
• モデル化誤差は、主に反応支配領域と自由拡散領域の間の遷移領域に局在しており、特に空間的区間$(3,6)$に顕著に現れる。
• 簡略化されたボルツマンモデルおよびIDSAの両方における離散化誤差は1次オーダーであるが、高解像度グリッドを用いることで1％未満に制御されている。
• IDSAは反応領域および輸送領域の両方でボルツマン解を低く見積もっており、最大の低減は遷移領域で発生する。
• IDSAは、定常状態近似による吸収・放出成分と拡散成分の強い結合のおかげで、反応支配領域において平衡状態への収束が速い。
• IDSAの拡散成分は、期待どおりに反応支配領域で消えるべきであるが、実際には残存しており、これが吸収・放出成分の減少を引き起こし、モデル化誤差の要因となっている。

(a) Reduced Boltzmann model with $T=0.0005$ : The dependency of the error in the infinity norm on $\Omega$ on the spatial mesh size $h=\Delta r$ is displayed. The curve with round markers represents the relative error, the curve with the plus markers displays the absolute error. The straight line ha