Spatiotemporal evolution of runaway electrons from synchrotron images in Alcator C-Mod

In the Alcator C-Mod tokamak, relativistic runaway electron (RE) generation can occur during the flattop current phase of low density, diverted plasma discharges. Due to the high toroidal magnetic field (B = 5.4 T), RE synchrotron radiation is measured by a wide-view camera in the visible wavelength range (~400-900 nm). In this paper, a statistical analysis of over one thousand camera images is performed to investigate the plasma conditions under which synchrotron emission is observed in C-Mod. In addition, the spatiotemporal evolution of REs during one particular discharge is explored in detail via a thorough analysis of the distortion-corrected synchrotron images. To accurately predict RE energies, the kinetic solver CODE [Landreman et al 2014 Comput. Phys. Commun. 185 847-855] is used to evolve the electron momentum-space distribution at six locations throughout the plasma: the magnetic axis and flux surfaces q = 1, 4/3, 3/2, 2, and 3. These results, along with the experimentally-measured magnetic topology and camera geometry, are input into the synthetic diagnostic SOFT [Hoppe et al 2018 Nucl. Fusion 58 026032] to simulate synchrotron emission and detection. Interesting spatial structure near the surface q = 2 is found to coincide with the onset of a locked mode and increased MHD activity. Furthermore, the RE density profile evolution is fit by comparing experimental to synthetic images, providing important insight into RE spatiotemporal dynamics

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.Comment: 23 pages, 12 figure

SOFT: A synthetic synchrotron diagnostic for runaway electrons

Improved understanding of the dynamics of runaway electrons can be obtained by measurement and interpretation of their synchrotron radiation emission. Models for synchrotron radiation emitted by relativistic electrons are well established, but the question of how various geometric effects -- such as magnetic field inhomogeneity and camera placement -- influence the synchrotron measurements and their interpretation remains open. In this paper we address this issue by simulating synchrotron images and spectra using the new synthetic synchrotron diagnostic tool SOFT (Synchrotron-detecting Orbit Following Toolkit). We identify the key parameters influencing the synchrotron radiation spot and present scans in those parameters. Using a runaway electron distribution function obtained by Fokker-Planck simulations for parameters from an Alcator C-Mod discharge, we demonstrate that the corresponding synchrotron image is well-reproduced by SOFT simulations, and we explain how it can be understood in terms of the parameter scans. Geometric effects are shown to significantly influence the synchrotron spectrum, and we show that inherent inconsistencies in a simple emission model (i.e. not modeling detection) can lead to incorrect interpretation of the images.Comment: 24 pages, 12 figure

Transport Phenomena in Alcator C-Mod H-Modes

Abstract. Several interesting new results have come from studies of ICRF-heated, H-mode plasmas in Alcator C-Mod. Dimensionless scaling studies have found gyro-Bohm-like transport similar to that reported on other devices; however, the dependence on collisionality was surprisingly strong, with BE −1. Despite high edge temperatures and strong edge pressure gradients, type I edge-localized modes (ELMs) are not observed in C-Mod. Instead we obtain a regime that we have dubbed enhanced D (EDA) which is accompanied by high-frequency density fluctuations. For all H-modes, core gradients were found to increase linearly with edge temperature, suggesting the importance of critical gradient/marginal stability behaviour. Comparisons with the IFS-PPPL model have begun, showing quantitative agreement in some cases. Impurity particle transport was studied via the laser blow-off technique with impurity confinement found to be effectively infinite for ELM-free discharges but reduced into the range 0.1–0.2 s for the EDA plasmas. 1

On the minimum transport required to passively suppress runaway electrons in SPARC disruptions

In [V.A. Izzo et al 2022 Nucl. Fusion 62 096029], state-of-the-art modeling of thermal and current quench (CQ) MHD coupled with a self-consistent evolution of runaway electron (RE) generation and transport showed that a non-axisymmetric (n = 1) in-vessel coil could passively prevent RE beam formation during disruptions in SPARC, a compact high-field tokamak projected to achieve a fusion gain Q > 2 in DT plasmas. However, such suppression requires finite transport of REs within magnetic islands and re-healed flux surfaces; conservatively assuming zero transport in these regions leads to an upper bound of RE current ~1 MA compared to ~8.7 MA of pre-disruption plasma current. Further investigation finds that core-localized electrons, within r/a < 0.3 and with kinetic energies 0.2-15 MeV, contribute most to the RE plateau formation. Yet only a relatively small amount of transport, i.e. a diffusion coefficient ~18 $\mathrm{m^2/s}$, is needed in the core to fully mitigate these REs. Properly accounting for (i) the CQ electric field's effect on RE transport in islands and (ii) the contribution of significant RE currents to disruption MHD may help achieve this