Large-scale Kinetic Simulations of Colliding Plasmas within a Hohlraum of Indirect Drive Inertial Confinement Fusions

Tianyi Liang,Dong Wu,Xiaochuan Ning,Lianqiang Shan,Zongqiang Yuan,Hongbo Cai,Zhengmao Sheng,Xiantu He
2023-09-20
Abstract:The National Ignition Facility has recently achieved successful burning plasma and ignition using the inertial confinement fusion (ICF) approach. However, there are still many fundamental physics phenomena that are not well understood, including the kinetic processes in the hohlraum. Shan et al. [Phys. Rev. Lett, 120, 195001, 2018] utilized the energy spectra of neutrons to investigate the kinetic colliding plasma in a hohlraum of indirect drive ICF. However, due to the typical large spatial-temporal scales, this experiment could not be well simulated by using available codes at that time. Utilizing our advanced high-order implicit PIC code, LAPINS, we were able to successfully reproduce the experiment on a large scale of both spatial and temporal dimensions, in which the original computational scale was increased by approximately 7 to 8 orders of magnitude. When gold plasmas expand into deuterium plasmas, a kinetic shock is generated and propagates within deuterium plasmas. Simulations allow us to observe the entire progression of a strong shock wave, including its initial formation and steady propagation. Although both electrons and gold ions are collisional (on a small scale compared to the shock wave), deuterium ions seem to be collisionless. This is because a quasi-monoenergetic spectrum of deuterium ions can be generated by reflecting ions from the shock front, which then leads to the production of neutrons with unusual broadening due to beam-target nuclear reactions. This work displays an unprecedented kinetic analysis of an existing experiment, shedding light on the mechanisms behind shock wave formation. It also serves as a reference for benchmark simulations of upcoming new simulation codes and may be relevant for future research on mixtures and entropy increments at plasma interfaces.
Plasma Physics
What problem does this paper attempt to address?
This paper attempts to solve some basic physical phenomena in the inertial confinement fusion (ICF) process, especially the dynamic processes generated by the collision between gold plasma and deuterium plasma in the indirect - drive ICF hohlraum. Specifically, the research aims to: 1. **Understand and simulate large - scale dynamic processes**: Due to the large spatiotemporal scale in the experiment, traditional particle simulation codes cannot accurately simulate these processes. The authors used their developed high - order implicit PIC (particle - in - cell) code LAPINS and successfully reproduced the experiment on a large spatiotemporal scale. 2. **Explain the abnormal broadening of the neutron energy spectrum**: Previous experiments observed an abnormally large width of the neutron energy spectrum, which may be related to the dynamic shock wave. Through simulation, the authors found that when the gold plasma expands into the deuterium plasma, a dynamic shock wave will be generated, and this shock wave can reflect the upstream deuterium ions, forming quasi - mono - energetic deuterium ions. These deuterium ions undergo nuclear reactions with the target material, resulting in the broadening of the neutron energy spectrum. 3. **Study the formation mechanism of the shock wave**: Through simulation, the authors were able to observe the entire evolution process of the strong shock wave, including its initial formation and stable propagation stages. Although electrons and gold ions are collisional on a small scale, deuterium ions behave almost collisionlessly. This is because the deuterium ions reflected at the shock front can generate a quasi - mono - energetic spectrum, and then produce neutrons with abnormal broadening through beam - target nuclear reactions. 4. **Provide references for future research**: This work not only provides an unprecedented dynamic analysis for existing experiments, but also provides references for future entropy increase research on mixtures and plasma interfaces. ### Main Conclusions - Through the high - order implicit PIC code LAPINS, the authors successfully simulated the collision process between gold plasma and deuterium plasma on a large spatiotemporal scale. - The simulation results show that when the gold plasma expands into the deuterium plasma, a dynamic shock wave will be generated, which reflects the upstream deuterium ions, forming quasi - mono - energetic deuterium ions. - These deuterium ions undergo nuclear reactions with the target material, resulting in the broadening of the neutron energy spectrum, explaining the abnormally wide neutron energy spectrum observed in the experiment. - The research results are helpful for understanding the formation mechanism of the shock wave and provide important references for future related research. ### Formula Summary - The velocity of the shock wave \( V_{sh} = 662 \, \text{km/s} \) - Mach number \( M=\frac{V_{sh}}{c_s} = 5.25 \), where \( c_s=\sqrt{\frac{k T_{e0}}{m_D}} = 126.4 \, \text{km/s} \) - The relationship between neutron energy broadening and temperature \( \Delta E_n = 82.5 T_D^{1/2} \) - The energy of the emitted neutrons \( E_n=\frac{E_D m_n}{4 m_D}\left[\left(2-\frac{m_n}{m_D}\right)\left(\frac{2Q}{E_D}+1\right)+2\sqrt{\left(2-\frac{m_n}{m_D}\right)\left(\frac{2Q}{E_D}+1\right)}\cos\theta + 1\right] \) These formulas help to explain the experimental phenomena and verify the simulation results.