The invention belongs to the field of magnetic confinement controlled
nuclear fusion, and discloses a numerical calculation method and
system for
phase space distribution function evolution of
helium ions in magnetic confinement
deuterium-
tritium fusion, and the method comprises the steps: carrying out the
numerical stability verification based on a self-adaptive CFL condition, and determining a proper step
size value, carrying out
discretization processing on the input data by utilizing a multi-dimensional grid; establishing a Fokker-Planck equation containing friction, energy
diffusion, position space transport and
helium ion generation source items, and converting the Fokker-Planck equation into a difference format capable of being solved in a recursive manner through a
finite difference method; in the calculation process, a corresponding boundary condition is set in combination with
physical reality, a difference format of a boundary area is corrected, self-consistent updating is carried out on the main
ion density based on a quasi-neutral condition, and information of evolution of a
helium ion phase space distribution function along with time is obtained; and introducing an energy boundary judgment method related to the
equivalent temperature of helium ash, distinguishing the helium ash from high-energy particles, and outputting a result of evolution of density profiles of the helium ash and the high-energy particles along with time.