Monte Carlo Simulation of D-D Neutron Emission Rate in NBI-Heated Deuterium Plasma of LHD

Abstract

As the first step to analyze the orbit of the triton produced by the fusion reaction between the beam-injected deuteron and the deuteron in the bulk plasma, we developed the code to calculate the velocity and spatial profiles of the fusion products in the Large Helical Device (LHD) that properly treats the particle loss boundary and the velocity anisotropy of the triton at birth. The neutron emission and neutron counting rates calculated by the code were compared with experimental results. The experimental results are reproduced by the developed code when the bulk plasma density, including the impurity ion, was used.