Extension of MIDAS-1D2V model: fusion reactions and neutral beam capture

MIDAS-1D2V is a plasma simulation code for axially symmetric open traps. It solves the non-stationary kinetic equation for the ion distribution function, accounting for Coulomb collisions. This paper presents two key modifications to the model. The first modification enables the calculation of fusion reactions in the plasma, with charged reaction products treated as an additional source term in the kinetic equation. The second modification incorporates beam capture dynamics, including ionisation and charge-exchange processes. Both modifications have been implemented for hydrogen isotopes.