Calculating the neutron yield of the Globus-M2 tokamak with allowance for the anisotropy of the ion velocity distribution function during neutral injection of high-energy atoms

Two algorithms for calculating the neutron yield of the Globus-M2 tokamak are considered for the case of anisotropic ion velocity distribution functions. Such calculations are necessary for estimating the neutron flux produced as a result of the interaction between fast particles (deuterons) arising from the injection of a high-energy atomic beam. The first algorithm comprises the calculation of the six-dimensional integral in the velocity space, and the second uses the expansion of the ion distribution function in Legendre polynomials, which allows one to reduce the dimension of integration. The results obtained using these algorithms are compared, and the advantages and disadvantages of both algorithms are discussed.