High Runaway Electrons Research Articles

Trends in runaway electron plateau partial recombination by massive H2 or D2 injection in DIII-D and JET and first extrapolations to ITER and SPARC

Experimental trends in thermal plasma partial recombination resulting from massive injection into high-Z (Ar) containing runaway electron (RE) plateaus in DIII-D and JET are studied for the purpose of achieving sufficiently low electron density () to increase RE final loss MHD levels. In both DIII-D and JET, thermal electron density is found to drop by ∼100 when the thermal plasma partially recombines, with a minimum at a vacuum vessel-averaged density in the range . RE effective resistivity also drops after partial recombination, indicating expulsion of the Ar content. The level after partial recombination is found to increase as RE current is increased. The amount of initial Ar in the RE plateau is not observed to have a strong effect on partial recombination. Partial recombination timescales of order 5 ms in DIII-D and 15 ms in JET are observed. These basic trends and timescales are matched with a 1D diffusion model, which is then used to extrapolate to ITER and SPARC tokamaks. Within the approximations of this model, it is predicted that ITER will be able to achieve sufficiently low values on time scales faster than expected RE plateau vertical drift timescales (of order 100 ms), provided sufficient or is injected. In SPARC, it is predicted that achieving significant recombination will be challenging, due to the very high RE current density. In both ITER and SPARC, it is predicted that achieving low will be easier with Ar as a background impurity (rather than Ne).

MARS-F modeling of post-disruption runaway beam loss by magnetohydrodynamic instabilities in DIII-D

A drift orbit model for relativistic test electrons has been incorporated into the MARS-F code (Liu et al 2000 Phys. Plasmas 7 3681), in order to study the runaway electron (RE) behavior in the presence of magneto-hydrodynamic perturbations computed by MARS-F. By implementing the model directly into the MARS-F curve-linear magnetic coordinates, maximal accuracy in representing the full field perturbation is preserved. The updated code is utilized to study the high current RE beam loss in a post-disruption DIII-D plasma, revealing that a fast growing, n = 1 (n is the toroidal mode number) resistive kink instability, at ∼100 Gauss level, can induce significant fraction of RE loss, largely by perturbing drift orbits of REs. A ∼1000 Gauss perturbation fully terminates the RE beam, as found in both experiment and modeling. The 3D field induced loss increases with the perturbation amplitude but decreases with the particle energy. The loss fraction is generally not sensitive to the initial particle pitch angle. The particle velocity change, due to electric field acceleration/deceleration, small pitch angle scattering, synchrotron radiation and Bremsstrahlung, further perturbs the RE trajectory but plays a minor role in prompt RE loss within microseconds time scale. Therefore, the dominant dependencies are simply the RE energy and instability strength. For comparison, a resonant magnetic perturbation field, generated by 4 kAt n = 3 even parity I-coil currents in DIII-D and with the plasma response field included, is found to induce almost no loss for the same RE beam.

Synchrotron and collisional damping effects on runaway electron distributions

Validated models of runaway electron (RE) dissipation are required to confidently approve safe ITER Q=10 operation. DIII-D experiments using quiescent REs are exploring the importance of synchrotron and collisional damping terms to RE dissipation. New time and energy-resolved measurements of RE bremsstrahlung hard X-ray (HXR) emission reveal stark differences between high and low energy REs as damping terms are varied. Previously reported anomalously high RE dissipation only applies to low energy REs. At high energy (where synchrotron effects are strongest) low synchrotron damping cases reach higher peak RE energy despite weaker particle confinement. Low-energy RE decay is observed concurrently with high-energy RE growth. RE dissipation models predict bump-on-tail distributions whose properties depend on the damping terms. Measured HXR spectra are very broad, as expected for bump-on-tail distributions.

Measurement of synchrotron radiation from runaway electrons during the TEXTOR tokamak disruptions

Investigations of runaway electrons (REs) during induced disruptions are performed at the TEXTOR tokamak. The synchrotron radiation generated by REs in the plasma core is detected using an infrared camera. The measurements enable the observation of the structure and dynamics of the runaway beam. In particular, the runaway beam is investigated as a function of the vertical and horizontal control fields. From the plasma current a runaway number of 4.33 × 1016 was calculated. The number of REs with energies exceeding 25 MeV is 6.4 × 1015 according to the synchrotron measurement at the end of the current quench phase. The mean pitch angle of the >25 MeV runaways is found to be 52 mrad which is smaller than in the low density discharge scenario. In addition, the synchrotron measurements show for the first time that a significant number of the high energy REs can survive after the end of the current plateau phase. In order to understand this, runaway orbit calculations have been performed which are in good agreement with the measurements.

Runaway Electron Suppression by ECRH and RMP in KSTAR

We have observed reduction of the runaway electron synchrotron radiation, hard X-ray (HXR) intensity, and HXR energy after applying 110 GHz 2nd harmonic electron cyclotron resonant heating (ECRH) during runaway electron (RE) discharges at low density with startup runaway electrons. However, we did not see a significant reduction of X-rays from 170 GHz 2nd harmonic ECRH at a higher field. A recently installed IR TV camera was used to observe the forward cone of synchrotron radiation from high energy REs in KSTAR. We have observed changes to the synchrotron images and reduction of the HXR by application of resonant magnetic perturbations (RMP) from in-vessel control coils (IVCC) installed inside KSTAR in the n = 1 configuration.

Novel rapid shutdown strategies for runaway electron suppression in DIII-D

New rapid shutdown strategies have been recently tested in the DIII-D tokamak to mitigate runaway electrons (REs). Disruptions in ITER are predicted to generate multi-MeV REs that could damage the machine. The RE population in large tokamaks is expected to be dominated by avalanche amplification which can be mitigated at high density levels by collisional drag. Particle injection schemes for collisional suppression of RE have been developed and tested in ITER-relevant scenarios: massive gas injection, shattered pellet injection (SPI) and shell pellet injection. The results show an improved penetration of particles injected with the SPI. Another strategy has been developed to harmlessly deconfine REs by applying a non-axisymmetric magnetic perturbation to worsen their confinement. This technique appeared to deconfine seed RE before the avalanche process could amplify the RE beam. The last method tested was to use the plasma position control system on the RE beam to prevent it from contacting the wall. This proved effective in preventing high current RE beam from touching the wall and providing more time to mitigate an existing RE beam but a successful ‘soft landing’ (without fast final losses) of the RE has not been observed yet.

Calculation of runaway electrons stopping power in ITER

The energy loss rate of runaway electrons (RE) was analysed for ITER plasma facing components materials (Be and W). The stopping power, the energy deposition profiles, and the material erosion are estimated by using the codes MEMOS and ENDEP. The latter has been updated by including the effect of the target’s polarizability. Our calculations show that this effect is significant for high RE energies and low Z materials such as Be. We also find that the conversion of the RE’s magnetic energy into heat can explain the temperature rise on dump plate in JET. In the case of ITER, the calculated heat deposition due to RE is almost two times the melting threshold energy of Be but well below that of W.