[論文レビュー] Primary Neutron Spectra in Ion Vlasov-Fokker-Planck Simulations

B. D. Appelbe, W. T. Taitano+4 more|arXiv (Cornell University)|May 3, 2023

Laser-induced spectroscopy and plasma被引用数 4

ひとこと要約

本稿では、イオン速度分布関数や微分断面積を近似せずに、完全に解像度を保った5次元数値手法を提示し、イオン・ブレージェル＝フォッカー・プランク（Vlasov-Fokker-Planck）シミュレーションにおいて一次中性子スペクトルを計算する。この手法により、非マクスウェル分布を示す非平衡プラズマにおける中性子スペクトルの正確なシミュレーションが可能となり、衝撃駆動型慣性核融合（ICF）インプロージョンにおいて特に顕著なスペクトルのずれが明らかになった。これにより、シミュレーションと実験のベンチマーク手法が向上する。

ABSTRACT

The energy spectra of unscattered neutrons produced by deuterium-deuterium and deuterium-tritium fusion reactions are an important diagnostic in High Energy Density Physics experiments as the spectra are sensitive to the velocities of reacting ions. Methods exist for calculating these spectra in radiation-hydrodynamic ("hydro") and Particle-in-Cell ("PiC") simulations. The spectra are particularly sensitive to the high energy tail of ion velocity distribution functions since reaction cross sections increase rapidly with the kinetic energy of a reacting ion pair at the energies achieved in laboratory experiments. This means both the hydro and PiC method may not be suitable in certain plasma regimes. The hydro method assumes that the ion velocity distribution is locally Maxwellian, while the PiC method is subject to statistical noise that makes it challenging to accurately simulate finer details of the spectra. In this work, we present a complementary approach: a method for calculating the neutron spectra in ion Vlasov-Fokker-Planck simulations in which the velocity distribution function is fully-resolved. The method is implemented in the spherically-symmetric code iFP which is used to simulate laser-driven Inertial Confinement Fusion experiments. The method is computationally intensive as it requires a five-dimensional numerical integral, but no approximations of the distribution functions or differential cross sections are required. Results show that deviations of the ion distribution functions from Maxwellian can have a noticeable effect on neutron spectra in shock-driven ICF implosions. The method should facilitate more accurate benchmarking of simulations and experiments.

研究の動機と目的

• イオン分布関数がマクスウェル分布から逸脱する状況において、従来の流体力学的および粒子法（PIC）手法が中性子スペクトルをシミュレートする際に抱える限界を解消すること。
• イオン速度分布関数および微分断面積の完全な解像度を保った高精度な一次中性子スペクトル計算手法の開発。
• 非マクスウェル分布が観察される慣性核融合（ICF）を含む高エネルギー密度物理学実験における、正確な合成診断の実現。
• 中性子エネルギースペクトルの計算において近似を排除することで、シミュレーションと実験の間の忠実度の高いベンチマークを向上させること。

提案手法

• 反応するイオンの速度分布関数および中心系フレームにおける微分断面積を5次元数値積分により中性子スペクトルを計算する。
• 中心系フレームにおける散乱角の余弦に関する積分に対して、再帰的関係式 $ F_j(x) = \int \frac{x^j}{\sqrt{1-x^2}} dx $ を用いて解析的解を得る。
• 粒子法に内在する統計的ノイズを回避し、放射流体力学コードと同様に局所的マクスウェル分布を仮定しない。
• 球対称な iFP コードに実装されており、レーザー駆動型慣性核融合実験をシミュレートする。
• 主要な方程式には、$ \gamma_a, \gamma_b, \gamma_c $ および $ p $ を含む核関数の変換が含まれ、$ x_{U,L} $ を用いて境界が定義され、角度積分の正確な評価が可能になる。
• $ v_{cm} = 0 $ や $ \sin\alpha = 0 $ といった極限ケースは、デルタ関数や簡略化された式を用いて解析的に取り扱い、数値積分に代わる。

実験結果

リサーチクエスチョン

• RQ1 衝撃駆動型慣性核融合インプロージョンにおいて、非マクスウェル分布のイオン速度分布関数が一次中性子スペクトルの形状にどのように影響を与えるか？
• RQ2 完全に解像度を保ったVlasov-Fokker-Planck手法は、イオン分布関数や断面積を近似せずに中性子スペクトルを計算可能か？また、流体力学的または粒子法との比較において、その性能はいかにか？
• RQ3 提案された5次元積分法は、標準的な数値積分と比較して、計算コストおよび数値的精度においてどの程度優れているか？
• RQ4 局所的マクスウェル分布からの逸脱が、シミュレーションと実験のベンチマークに用いられる合成中性子スペクトルに、どの程度の測定可能な変化をもたらすか？

主な発見

• 本手法により、イオン速度分布関数や微分断面積を近似せずに、スペクトル計算の完全な忠実度を確保した一次中性子スペクトルの計算が可能になった。
• 衝撃駆動型ICFインプロージョンにおいて、イオン分布関数がマクスウェル分布から逸脱していることが、特に高エネルギー尾部に顕著なスペクトルのずれを引き起こすことが判明した。
• 角度積分の解析的解法により、特に高エネルギー反応（例：$ K \sim 10\,\text{MeV} $、ここで $ N \sim 10 $）において、数値積分に比べて計算コストを顕著に低減できた。
• 本手法は重水素-重水素反応および重水素-トリチウム反応の両方に対応でき、任意のイオン分布関数を扱える。
• $ v_{cm} = 0 $ や $ \sin\alpha = 0 $ といった極限ケースの取り扱いにより、特殊な構成でも効率的かつ正確に評価可能となり、数値的安定性が向上した。
• より現実的な合成診断としての中性子エネルギースペクトルを提供するため、本手法により、シミュレーションと実験データの間のより正確なベンチマークが可能になった。