Kinetic modelling of runaway electron avalanches in tokamak plasmas
Runaway electrons (REs) can be generated in tokamak plasmas if the accelerating force from the toroidal electric field exceeds the collisional drag force due to Coulomb collisions with the background plasma. In ITER, disruptions are expected to generate REs mainly through knock-on collisions, where enough momentum can be transferred from existing runaways to slow electrons to transport the latter beyond a critical momentum, setting off an avalanche of REs. Since knock-on runaways are usually scattered off with a significant perpendicular component of the momentum with respect to the local magnetic field direction, these particles are highly magnetized. Consequently, the momentum dynamics require a full 3-D kinetic description, since these electrons are highly sensitive to the magnetic non-uniformity of a toroidal configuration. A bounce-averaged knock-on source term is derived. The generation of REs from the combined effect of Dreicer mechanism and knock-on collision process is studied with the code LUKE, a solver of the 3-D linearized bounce-averaged relativistic electron Fokker-Planck equation, through the calculation of the response of the electron distribution function to a constant parallel electric field. This work shows that the avalanche effect can be important even in non-disruptive scenarios. RE formation through knock-on collisions is found to be strongly reduced when taking place off the magnetic axis, since trapped electrons cannot contribute to the RE population. The relative importance of the avalanche mechanism is investigated as a function of the key parameters for RE formation; the plasma temperature and the electric field strength. In agreement with theoretical predictions, the simulations show that in low temperature and E-field knock-on collisions are the dominant source of REs and can play a significant role for RE generation, including in non-disruptive scenarios.