Runaway electron generation during plasma shutdown by killer pellet injection
Tokamak discharges are sometimes terminated by disruptions that may cause large mechanical and thermal loads on the vessel. To mitigate disruption-induced problems it has been proposed that ‘killer’ pellets could be injected into the plasma in order to safely terminate the discharge. Killer pellets enhance radiative energy loss and thereby lead to rapid cooling and shutdown of the discharge. But pellets may also cause runaway electron generation, as has been observed in experiments in several tokamaks. In this work, runaway dynamics in connection with deuterium or carbon pellet-induced fast plasma shutdown is considered. A pellet code, which calculates the material deposition and initial cooling caused by the pellet is coupled to a runaway code, which determines the subsequent temperature evolution and runaway generation. In this way, a tool has been created to test the suitability of different pellet injection scenarios for disruption mitigation. If runaway generation is avoided, the resulting current quench times are too long to safely avoid large forces on the vessel due to halo currents.
Numerical calculation of the runaway electron distribution function and associated synchrotron emission
Synchrotron emission from runaway electrons may be used to diagnose plasma conditions during a tokamak disruption, but solving this inverse problem requires rapid simulation of the electron distribution function and associated synchrotron emission as a function of plasma parameters. Here we detail a framework for this forward calculation, beginning with an efficient numerical method for solving the Fokker-Planck equation in the presence of an electric field of arbitrary strength. The approach is continuum (Eulerian), and we employ a relativistic collision operator, valid for arbitrary energies. Both primary and secondary runaway electron generation are included. For cases in which primary generation dominates, a time-independent formulation of the problem is described, requiring only the solution of a single sparse linear system. In the limit of dominant secondary generation, we present the first numerical verification of an analytic model for the distribution function. The numerical electron distribution function in the presence of both primary and secondary generation is then used for calculating the synchrotron emission spectrum of the runaways. It is found that the average synchrotron spectra emitted from realistic distribution functions are not well approximated by the emission of a single electron at the maximum energy. © 2013 Elsevier B.V.
Runaway electrons in an impure plasma
The electron runaway theories of Kruskal and Bernstein and Lebedev are extended to apply to multiply ionized, multiple species plasmas. The resulting expressions for electron runaway rates are compared with results obtained by Kulsrud, Sun, Winsor, and Fallon from numerical solutions of the Fokker–Planck equation.