A two-time-layer scheme has been developed for solving the simplified kinetic Fokker–Planck equation related to the transport of charged products of thermonuclear reactions, which includes an interpolation procedure in four-dimensional grid space. Instabilities in the scheme were detected at low particle velocities and for a specific choice of particle deceleration in the ion field, which enters the kinetic equation as a parameter. It was shown that the thermalization condition, which prohibits solving the kinetic equation for particles with energy lower than the average ion energy, significantly limits the number of thermonuclear reactions where instability can manifest. The scheme was tested on the problem of relaxation to a stationary state and on a problem with a prescribed time-dependent thermonuclear reaction rate, for which an exact solution to the kinetic equation can be found.
Read full abstract