1.School of National Defense Science and Technology, Southwest University of Science and Technology, Mianyang 621010, China 2.College of Liberal Arts and Science, National University of Defense Technology, Changsha 410073, China 3.University of Science and Technology of China, Hefei 230027, China
Fund Project:Project supported by the Doctoral Fund of Southwest University of Science and Technology, China (Grant No. 13zx7138), the Natural Science Foundation of Hunan Province, China (Grant No. 2018JJ2473), the National Natural Science Foundation of China (Grant Nos. 11475260, 11475150), the Opening Foundation of Key Laboratory of Neutronics and Radiation Safety of CAS (Grant No. NEUTRON201707), the Fund of Fundamental Science on Nuclear Safety and Simulation Technology Laboratory, Harbin Engineering University, China (Grant No. HEUNSS18SF04), and the Fund of Robot Technology Used for Special Environment Key Laboratory of Sichuan Province, China (Grant No. 17kftk01)
Received Date:28 March 2019
Accepted Date:13 August 2019
Available Online:01 November 2019
Published Online:05 November 2019
Abstract:Inertial confinement fusion (ICF) is one of the possible ways to realize controlled thermonuclear fusion. The fusion neutron source term is one of the important parameters in the physical design and analysis of laser plasma. The accuracy of the fusion neutron source term directly affects the reliability of the analysis results. At present, the neutron source term of deuterium-tritium fusion reaction in ICF is mainly based on formula method. It has limited applications in temperature and reaction type. Because of a large quantity of data, it is impossible to simulate the fusion reaction of each particle. In this paper, the concept of particle cloud is introduced, that is, the collection of the like particles with the same position and speed, and it is considered that the action of particle cloud is the same reaction. Because the particles should satisfy the Maxwell velocity distribution at a certain temperature and the direction is all around the circumference angle, the collision cross sections between the incident particle and different target particles are different. Therefore, the design program takes all the possible velocities, polar angles and direction angles, reads the collision cross sections between deuterium and tritium and makes corrections, and obtains the multi-temperature differential correction cross sections of deuterium and tritium fusion with Doppler energy broadening. On these bases, Monte Carlo method and discrete ordinate method method are used. A numerical simulation program for the fusion rate of D-T particles in laser plasma is developed in this paper. It is found that there are significant differences between the DT, DD, TD cross sections and the original cross sections after Doppler broadening. In a range of plasma temperature between 20 keV and 100 keV, the simulation results are more consistent with the cross section data of ENDF/B-VI and ENDF/B-VII databases of deuterium-tritium fusion reaction than those from the analytical formula method. There is a large error between the numerical simulation results and the analytical formula method in the low energy region. It may be caused by the difference of calculation methods and too big difference among the used fusion cross sections at low temperature. Keywords:laser plasma/ deuterium-tritium fusion reaction/ differential cross-section temperature correction/ Monte Carlo method/ discrete ordinate method