[论文解读] 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

一句话总结

本论文提出了一种完全解析的五维数值方法，用于在离子 Vlasov-Fokker-Planck 模拟中计算初级中子能谱，无需对离子速度分布函数或微分截面进行近似。该方法能够精确模拟非麦克斯韦分布离子等非平衡等离子体中的中子能谱，揭示了由于非平衡离子分布导致的显著谱形偏差，特别是在冲击驱动的惯性约束聚变内爆中，从而提升了模拟与实验之间的基准对比精度。

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.

研究动机与目标

• 解决现有流体动力学方法和粒子-网格法在离子分布偏离麦克斯韦分布时模拟中子能谱的局限性。
• 开发一种高保真度方法，用于计算初级中子能谱，同时保持离子速度分布函数和微分截面的完整分辨率。
• 实现高能量密度物理实验中更精确的合成诊断，特别是在观察到非麦克斯韦离子分布的惯性约束聚变中。
• 通过消除中子能谱计算中的近似，提升模拟与实验之间对比的保真度。

提出的方法

• 该方法通过在质心系中对反应离子的速度分布函数和微分截面进行五维数值积分来计算中子能谱。
• 利用递推关系推导出的解析解计算质心系中散射角余弦的积分，即 $ 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 与标准数值积分相比，所提出的五维积分方法在计算成本和数值精度方面表现如何？
• RQ4 局部麦克斯韦分布的偏离在多大程度上导致了用于模拟与实验基准对比的合成中子能谱的可测量变化？

主要发现

• 该方法可在不近似离子速度分布函数或微分截面的情况下计算初级中子能谱，确保了谱形计算的完整保真度。
• 研究发现，离子分布函数偏离麦克斯韦分布对冲击驱动 ICF 内爆中的中子能谱有显著影响，尤其体现在高能尾部。
• 对角度积分的解析解显著降低了计算成本，相较于数值积分尤其在高能反应（如 $ K \sim 10\,\text{MeV} $，其中 $ N \sim 10 $）中优势明显。
• 该方法适用于氘-氘和氘-氚聚变反应，并支持任意离子分布函数。
• 对 $ v_{cm} = 0 $ 和 $ \sin\alpha = 0 $ 等极限情况的包含，使得在特殊构型下能够实现高效且精确的计算，增强了数值稳定性。
• 该方法通过提供更真实的中子能谱合成诊断，促进了模拟与实验数据之间更精确的基准对比。