Search NASA⌕ Search

SEARCH · Search NASA

Results for “Runaway electrons”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Nonlinear saturation of whistler modes driven by runaway electrons

Runaway electrons exhibit kinetic instabilities with potentially beneficial consequences. The anomalous Doppler resonance between the electrons and whistler modes is a primary underlying mechanism. These instabilities require a first-principle nonlinear theoretical analysis, which would ultimately enable assessment of their impact on disruption mitigation and diagnostics, especially in ITER-relevant conditions. This paper presents recent progress in developing the required theoretical framework for isolated nonlinear resonances. By employing the generic bump-on-tail model, we predict the wave saturation levels in both near-threshold and strongly driven regimes. Remarkably, for the waves of interest, the parallel component of the wave vector dictates the parallel momentum of the resonant particles. This feature, together with the expected wave saturation level, provides a complete description of the resonance impact on the runaway electron distribution function. We show that the parametric dependence of the nonlinearly saturated mode provides a way to recover certain features of the runaway momentum distribution function.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Wall heating by subcritical energetic electrons generated by the runaway electron avalanche source *

Abstract Subcritical energetic electrons (SEEs) produced by the runaway electron (RE) avalanche source at energies below the runaway threshold are found to be the primary contributor to surface heating of plasma-facing components (PFCs) during final loss events. This finding is supported by theoretical analysis, computational modeling with the Kinetic Orbit Runaway electrons Code (KORC), and qualitative agreement with DIII-D experimental observations. The avalanche source generates significantly more secondary electrons below the runaway threshold, which thermalize rapidly when well-confined. However, during a final loss event, the RE beam impacts the first wall, and SEEs are deconfined before they can thermalize. Additionally, because the energy deposition length decreases faster than energy, the deposited energy density, and thus the maximum PFC surface temperature change, is larger for SEEs than REs. KORC simulations employ an analytic first wall to model particle deconfinement onto a non-axisymmetric wall composed of individual tiles. PFC surface heating is calculated using a 1D model extended to include an energy-dependent deposition length scale. Simulations of DIII-D qualitatively agree with infrared (IR) imaging only when SEEs from the avalanche source are included. These results demonstrate that SEEs are the dominant contributor to PFC surface heating and indicate that the avalanche source plays a critical role in the PFC damage caused during final loss events. The prominence of SEEs also has important implications for interpreting IR imaging, one of the primary diagnostics for RE-wall interaction diagnosis, despite REs dominating the energy and current density. This result improves predictions of wall damage due to post-disruption REs to estimate material lifetime and design RE mitigation systems for ITER and future reactors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1978

Spontaneous emission near the electron plasma frequency in a plasma with a runaway electron tail

Spontaneous emission of radiation with frequencies near the electron plasma frequency is studied for a plasma which consists of both thermal and runaway electrons. It is found that a substantial enhancement of the spontaneous radiation intensity can occur in this frequency regime via a Cherenkov resonance with the runaway electrons. Numerical analysis indicates that, for reasonable estimates of densities and energies, the plasma-frequency radiation can attain levels greater than the peak thermal emission at the second gyroharmonic.

Freund, H. P.↗
DOE OSTI2024

Resistive hose modes in tokamak runaway electron beams

Beams of energetic runaway electrons are generated during disruptions in tokamaks, and fluid models are used to study their effects on macroscale dynamics. Linear computations of a massless, runaway electron beam coupled to MHD plasma show that resistive hose instabilities grow faster than tearing modes at large resistivity. Eigenvalue results with reduced models of the resistive hose instability are compared with results from the full MHD and beam system, showing that the resistive hose decouples from any plasma response. An estimate of plasma temperature at which growth of the resistive hose dominates tearing for post-disruption DIII-D plasma parameters is in a physically relevant regime.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Simulation Center for Runaway Electron Avoidance and Mitigation (SCREAM SciDAC) (Technical Final Report)

Runaway electrons can severely damage the plasma facing components on ITER during a major disruption and pose a major risk for tokamak fusion. It has been recognized that an adequate disruption mitigation system (DMS) is essential for the safe operation of ITER. The United States is responsible for the design and implementation of the disruption mitigation system on ITER, and in July 2016 the Simulation Center for Runaway Electron Avoidance and Mitigation (SCREAM) was launched by DOE, in a joint Fusion Energy Sciences (FES) and Advanced Scientific Computing Research (ASCR) collaboration. SCREAM was a comprehensive theory and simulation SciDAC center that provided physics guidance in the avoidance and mitigation of runaway electrons, and in tandem with domestic and international experiments, helped establish the qualitative and quantitative bases for safe operational scenarios and viable mitigation techniques. The SCREAM center assembled a national team of experts in runaway electron physics, tokamak disruptions, magnetohydrodynamic (MHD) simulation, and advanced algorithms and computing. The team combined advanced simulation and analysis capability facilitated by direct participation of ASCR SciDAC institutes with theoretical models and code development by FES scientists to focus on the runaway risk for ITER and tokamaks in general. The research scope was focussed on integrated simulations of kinetic runaway electrons, including MHD and fluid models of impurity transport, within a research plan guided by theory. The specific research tasks were (1) establish the fundamental physics of runaway generation, saturation, and dynamical evolution in a tokamak; (2) examine the critical path toward runaway avoidance; and (3) investigate the viability and effectiveness of the leading candidate schemes for runaway mitigation. In all three areas, members of the team carried out scoping studies that established the readiness for rapid and critical advances, especially in the deployment and further development of large-to extreme-scale simulation tools. Our multi-pronged computational approach included (1) relativistic Fokker-Planck solvers with discretization in phase space, (2) self-consistent particle-in-cell techniques, (3) particle-based Monte-Carlo, and (4) MHD-particle hybrid simulations. Cross-check between these different methods provided an additional means for verification and further bolstered the fidelity of our physics prediction. Validation against experimental results brings confidence to the predictive capability for ITER and frequently leads to new ideas for understanding and mitigating the thermal quench driven runaway electron phenomenon.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Orbit-averaging and deposition accuracy for runaway electron beams in hybrid kinetic-MHD simulations of the runaway plateau

We develop a new procedure that combines the kinetic orbit runaway electron code (KORC) and the NIMROD extended-magnetohydrodynamic code to simulate runaway electrons (REs) in the post-disruption plateau. KORC integrates guiding-center orbits, with a barycentric-based binary search strategy providing initial guesses for the Newton–Raphson logical-to-physical coordinate inversion, ensuring reliable particle-to-mesh mapping in NIMROD, whose fields remain static for the present study. Samples are drawn in accord with experimental parallel current profiles of RE beams during the plateau phase. Deposition in NIMROD is verified through comparison with a Python-based finite-element code that ensures periodicity in the poloidal direction and continuity at the magnetic axis. Accurate representation of near-axis fields requires finer mesh resolution to prevent under- and overshoots in current density from orbit inaccuracies. Yet, at a fixed particle count, increasing mesh resolution amplifies statistical noise in the deposited fields. An orbit-averaging method accumulates partial current deposits over multiple kinetic steps and reduces the statistical noise with little added computational cost. By coupling kinetic routines from KORC directly into the NIMROD codebase, these developments lay essential groundwork for future self-consistent KORC–NIMROD coupling.

Algorithms and data structure↗
DOE OSTI2025

An upper pressure limit for low- Z benign termination of runaway electron beams in TCV

Abstract We present a model for the particle balance in the post-disruption runaway electron plateau phase of a tokamak discharge. The model is constructed with the help of, and applied to, experimental data from TCV discharges investigating the so-called ‘low- Z benign termination’ runaway electron mitigation scheme. In the benign termination scheme, the free electron density is first reduced in order for a subsequently induced MHD instability to grow rapidly and spread the runaway electrons widely across the wall. We show that the observed non-monotonic dependence of the free electron density with the measured neutral pressure is due to plasma re-ionization induced by runaway electron impact ionization. At higher neutral pressures, more target particles are present in the plasma for runaway electrons to collide with and ionize. Parameter scans are conducted to clarify the role of the runaway electron density and energy on the free electron density, and it is found that only the runaway electron density has a noticeable impact. While the free electron density is shown to be related to the spread of heat fluxes at termination, the exact cause for the upper neutral pressure limit remains undetermined and an object for further study.

Hoppe, M. (ORCID:0000000339948977)↗
DOE OSTI2025

Excitation of whistler and slow-X waves by runaway electrons in a collisional plasma

Runaway electrons are known to provide robust ideal or collisionless kinetic drive for plasma wave instabilities in both the whistler and slow-X branches, via the anomalous Doppler-shifted cyclotron resonances. In a cold and dense post-thermal-quench plasma, collisional damping of the plasma waves can compete with the collisionless drive. Previous studies have found that, due to their higher wavelength and frequency, slow-X waves suffer stronger collisional damping than the whistlers, while the ideal growth rate of slow-X modes is higher. Here, we study runaway avalanche distributions that maintain the same eigen distribution and increase only in magnitude over time. The distributions are computed from the relativistic Fokker–Planck–Boltzmann solver, upon which a linear dispersion analysis is performed to search for the most unstable or least damped slow-X and whistler modes. Taking into account the effect of plasma density, plasma temperature, and effective charge number, we find that the slow-X modes tend to be excited before the whistlers in a runaway current ramp-up. Furthermore, even when the runaway current density is sufficiently high that both branches are excited, the most unstable slow-X mode has a much higher growth rate than the most unstable whistler mode. The qualitative and quantitative trends uncovered in the current study indicate that even though past experiments and modeling efforts have concentrated on whistler modes, there is a compelling case that slow-X modes should also be a key area of focus in the runaway self-mediation through wave instabilities.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Radiation modeling and experimental validation of sub-cyclotron frequency microwave synchrotron radiation as a diagnostic tool for runaway electrons in present and future tokamaks

The microwave synchrotron radiation in the sub-cyclotron frequency region has been systematically validated as a diagnostic method for runaway electrons (>1 MeV) with radiation modeling and experiments on the DIII-D tokamak. Experimental observations reveal that sub-cyclotron frequency microwave radiation emerges only after runaway electrons are excited during plasma disruptions. Electron gyro-magnetic radiation modeling confirms that only runaway electrons with energies exceeding 1 MeV can emit radiation in the observed sub-cyclotron frequency band on DIII-D. The synchrotron Razin and Beaming effects are modeled in both the sub-cyclotron and standard Electron Cyclotron Emission (ECE) frequency region. Due to the pronounced beaming effect, the modeling supports that a standard ECE system is insensitive to synchrotron radiation from beamed runaway electrons (p ∥ >> p ⟂ ), but very sensitive to ECE radiation from isotropic suprathermal electrons (100 keV or less). DIII-D experiments further show that combining ECE and sub-cyclotron receivers reveals new insights into transient runaway dynamics during the plateau phase. Certain transient crash events appear only in the sub-cyclotron band or lead to opposite changes on the sub-cyclotron and standard ECE signals, supporting models where transient MHD events can exert a different effect on electrons of different energy levels. This novel diagnostic technique offers a powerful tool for studying and detecting runaway electron dynamics in both present-day and high-field future tokamaks. The diagnostic approach has been extrapolated to the SPARC tokamak, where sub-cyclotron radiation is predicted to respond sensitively and solely to runaway electrons during SPARC’s L-mode startup phase.

Yu, Guanying [Hefei Comprehensive National Science↗
NASA NTRS1982

Electrons in the ionospheric source cone - Evidence for runaway electrons as carriers of downward Birkeland currents

Extremely intense field-aligned fluxes of low energy electrons have been observed with ISIS-2 streaming out of the ionosphere at auroral latitudes. Fluxes in excess of 10 billion per sq cm sec ster at energies below 500 eV with peak fluxes from 10 to 100 eV were detected at 1400 km. The electrons are very strongly field-aligned, having pitch angles confined within 10 deg of the magnetic field. Since they are so intense and so highly collimated they cannot be produced by atmospheric backscattering of a primary auroral beam. These electrons are sometimes associated with ionospheric ions that have been accelerated transverse to the magnetic field. They occur in regions of downward field-aligned current, and may thus be carriers of the current, thus far unidentified. It is suggested that they are runaway electrons from the ionosphere produced by a downward field-aligned component of the electric field.

Klumpar, D. M.↗
DOE OSTI2023

Marginal stability constraint on runaway electron distribution

High-frequency kinetic instabilities of the strongly anisotropic runaway electrons (RE) can enhance the pitch-angle scattering of the runaways significantly. This wave-induced scattering can easily prevail over runaway scattering on high-Z impurities. In a steady state, collisional damping balances the kinetic drive of the unstable waves, keeping the RE distribution function at marginal stability. The marginal stability constraint limits the achievable RE densities and the shape of the RE distribution function. In this study, we consider whistler and compressional Alfvén waves as the primary source of enhanced elastic scattering of the runaways. By balancing the anomalous Doppler resonance drive with the collisional wave damping, we find the RE distribution function in the ultra-relativistic range of the phase space. We also derive an expression for the spectral energy density of the waves. We show that the power needed to compensate for the wave dissipation is negligible compared to the work of the electric field. The latter is in balance with the synchrotron losses of the runaway electrons.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Coupled 2D MHD and runaway electron fluid simulations of SPARC disruptions

Runaway electrons (REs) generated during disruption events in tokamaks can carry mega-Ampère level currents, potentially causing damage to plasma-facing components. Understanding RE evolution during disruption events is important for evaluating strategies to mitigate RE damage. Using two-dimensional toroidally symmetric magnetohydrodynamic (MHD) simulations in M3D-C1, which incorporates a fluid RE model evolved self-consistently with the bulk MHD fluid, we examine the seeding and avalanching of REs during disruptions in the SPARC tokamak – a compact, high-field, high-current device designed to achieve a fusion gain Q > 2 in deuterium–tritium plasmas. The M3D-C1 simulations of unmitigated disruptions demonstrate RE plateau formation and peaking of the final current density, which agree well with the results of lower-fidelity reduced RE fluid models. This work provides the first systematic comparison and benchmarking of different primary sources, including activated tritium beta decay and Compton scattering, in SPARC disruption simulations with self-consistent MHD and RE coupling.

Compton scattering↗
DOE OSTI2025

Resistive hose modes in tokamak runaway electron beams II

Resistive hose instabilities of runaway electron (RE) beams immersed in resistive background plasmas are examined with analytic and numerical calculations. The RE beam-plasma equilibria considered are characteristic of the situation observed post-thermal quench in a tokamak disruption. An analytic linear dispersion relation is presented for the case of a uniform RE current density profile with a sharp boundary in cylindrical geometry. Initial value linear calculations for a more general profile in toroidal geometry find that reducing aspect ratio increases the resistive hose mode growth rates with fixed safety factor profile. Nonlinear calculations in cylindrical and toroidal geometry find that the resistive hose instability-driven fluctuations relax the gradient of the current density profile. In toroidal geometry, changes to the magnetic topology are observed as a result of the resistive hose activity.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Simulation of DIII-D disruption with argon pellet injection and runaway electron beam

Abstract The next generation of large tokamaks, including ITER, will be equipped with a disruption mitigation system (DMS) that can be activated if a disruption is deemed to be imminent. Introducing impurities by pellet (large or shattered) or massive gas injection has been shown to be an effective mitigation mechanism on many tokamaks. The goal of the mitigation is to lessen the thermal and electromagnetic loads from the disruption without generating enough high-energy (runaway) electrons to damage the device. Variations of this mitigation process with impurity injection are presently being tested on many experiments. We have modeled one such impurity injection experiment on DIII-D using the M3D-C1 nonlinear 3D extended MHD code (Jardin et al 2012 Comput. Sci. Discovery 6 014002), The model includes an argon large pellet injection and ablation model, impurity ionization, recombination, and radiation, and runaway electron formation and subsequent evolution, including both Dreicer and avalanche sources. We obtain reasonable agreement with the experimental results for the timescale of the thermal and current quench and for the magnitude of the runaway electron plateau formed during the mitigation. This is the first 3D full MHD simulation with pellets and REs to simulate the disruption process and it also provides a partial validation of the M3D-C1 DMS model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Analysis of runaway electron driven whistler wave instability experiments

Data acquired on the DIII-D tokamak were analyzed. The data are from experiments that were conducted to study an instability that is driven unstable by intense populations of electrons with MeV energies that are known as runaway electrons. The instability is a type of plasma wave called a whistler wave that occurs at frequencies above the ion cyclotron frequency but well below the electron cyclotron frequency. The waves were measured by magnetic fluctuation coils that are embedded in the DIII-D vacuum vessel wall. After upgrades to this diagnostic were completed, new experiments were conducted on July 13, 2020 in order to measure the toroidal mode number of the whistler waves and to extend the frequency of the detected waves. Through the use of mixers, instability between 600-700 MHz was detected. (The initial experiments only measured up to 200 MHz.) Analysis of the data was led by Hari Choudhury, a PhD student at Columbia University. Mr. Choudhury has submitted two papers for publication that include contributions by UC Irvine (UCI) Professor Heidbrink and his graduated PhD student Genevieve DeGrandchamp: “Detailed Characterization of Runaway Electron Driven Whistler Waves in Low-Density DIII-D Discharges” and “First Demonstration of Resonant Pitch-Angle Scattering of Relativistic Electrons by Externally-Launched Helicon Waves.” The first paper, which has been submitted to Physics of Plasmas, has significant contributions to both the data and the interpretation by UCI scientists. In contrast, UCI contributions to the second paper, which has been submitted to Physical Review Letters, are relatively minor.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

A physics-constrained deep learning surrogate model of the runaway electron avalanche growth rate

A surrogate model of the runaway electron avalanche growth rate in a magnetic fusion plasma is developed. This is accomplished by employing a physics-informed neural network (PINN) to learn the parametric solution of the adjoint to the relativistic Fokker–Planck equation. The resulting PINN is able to evaluate the runaway probability function across a broad range of parameters in the absence of any synthetic or experimental data. This surrogate of the adjoint relativistic Fokker–Planck equation is then used to infer the avalanche growth rate as a function of the electric field, synchrotron radiation and effective charge. Predictions of the avalanche PINN are compared against first principle calculations of the avalanche growth rate with excellent agreement observed across a broad range of parameters.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

The impact of collisionality on the runaway electron avalanche during a tokamak disruption

The exponential growth (avalanching) of runaway electrons (REs) during a tokamak disruption continues to be a large uncertainty in RE modeling. The present work investigates the impact of tokamak geometry on the efficiency of the avalanche mechanism across a broad range of disruption scenarios. It is found that the parameter ν *,crit ⁠, describing the collisionality at the critical energy to run away, delineates how toroidal geometry impacts RE formation. In particular, utilizing a reduced but self-consistent description of plasma power balance, it is shown that for a high-density deuterium-dominated plasma, ν *,crit is robustly less than one, resulting in a substantial decrease in the efficiency of the RE avalanche compared to predictions from slab geometry. In contrast, for plasmas containing a substantial quantity of neon or argon, ν *,crit ≳1, no reduction of the avalanche is observed due to toroidal geometry. This sharp contrast in the impact of low- vs high-Z material results primarily from the relatively strong radiative cooling from high-Z impurities, enabling the plasma to be radiatively pinned at low temperatures and thus large electric fields, even for modest quantities of high-Z material.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗