Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “runaway”

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 OSTI2024

Modeling delayed thermal runaway in nitric acid-soaked cat litter mixed with radioactive waste

Thermal ignition of radioactive waste within a 55-gallon drum was simulated by using a pressure-dependent waste decomposition model (Hobbs et al. in Process Saf Environ Prot https://doi.org/10.1016/j.psep.2022.09.047, 2022) calibrated with data from full-scale drum experiments (Parker et al. in The thermolytic response of a surrogate RNS waste mixture at the drum scale. Los Alamos National Laboratory Report LA-UR-16-21760, 2016) and validated with experiments from multiple laboratories (Hobbs et al. in Thermal analysis of aged nitric acid-soaked kitty litter in TRU waste drums-23370.WM2023 Conference, Phoenix, AZ, 2023). The acceleration of nitric acid chemistry reacting with an organic cat litter leading to thermal ignition was likely triggered by a restricted vent in the drum. Here, we address whether the form of the rate equation in (Hobbs et al. in Process Saf Environ Prot https://doi.org/10.1016/j.psep.2022.09.047, 2022) is sufficient to extrapolate thermal ignition within aged drums of similar content that have been stored in Texas for over nine years by investigating four different reaction rate forms for waste decomposition. A critical reaction rate reduction analysis is performed on each of these models to determine if delayed thermal runaway within vented aged waste is possible after nine years. We found that a pressure-dependent first-order rate expression not only predicted the accidental ignition of the waste drum, but the form also matches multiple experiments from different laboratories. Even though the waste composition decreases over time, the model predicts that acceleration leading to thermal runaway is possible if the waste is confined, even after 9 years. In conclusion, waste containing oxidizers such as nitric acid should not be mixed with organic adsorbents, especially if the waste is confined.

cookoff↗
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 OSTI2024

Simulations of stand-off runaway electron beam termination by tungsten particulates for tokamak disruption mitigation

Abstract Stand-off runaway electron termination by injected tungsten particulates offers a plausible option in the toolbox of disruption mitigation. Tungsten is an attractive material choice for this application due to large electron stopping power and high melting point. To assess the feasibility of this scheme, we simulate runaway collisions with tungsten particulates using the MCNP program for incident runaway energies ranging from 1 to 10 MeV. We assess runaway termination from energetics and collisional kinematics perspectives. Energetically, the simulations show that 99% of runaway beam energy is removed by tungsten particulates on a timescale of 4–9 µ s. Kinematically, the simulations show that 99% of runaways are terminated by absorption or backscattering on a timescale of 3–4 µ s. By either metric, the runaway beam is effectively terminated before the onset of particulate melting. Furthermore, the simulations show that secondary radiation emission by tungsten particulates does not significantly impact the runaway termination efficacy of this scheme. Secondary radiation is emitted at lower particle energies than the incident runaways and with a broad angular distribution such that the majority of secondary electrons emitted will not experience efficient runaway re-acceleration. Overall, the stand-off runaway termination scheme is a promising concept as a last line of defense against runaway damage in ITER, SPARC, and other future burning-plasma tokamaks.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Corrigendum: Simulations of stand-off runaway electron beam termination by tungsten particulates for tokamak disruption mitigation (2024 Nucl. Fusion 64 056019)

In the original article, we inadvertently omitted a reference which previously introduced the concept of tungsten injection for 'Disruption Mitigation in Tokamak Reactor via Reducing the Seed Electrons of Avalanche.' In that reference, the authors have proposed injection of a tungsten cylinder of 80 mm by 8 mm dimensions, by rail guns, to absorb the runaway seed population. This can be contrasted with another interesting idea of using small tungsten pellets coated with a low Z material for depleting the runaway seeds as an option for upgrade of the ITER Disruption Mitigation System. Interested readers are referred to those two papers on the tungsten injection for runaway seed removal/reduction, and the related issue of not shortening the current quench excessively. The application focus of was not on reducing the runaway seeds in the initial Ohmic-to-runaway current conversion phase, but on safe termination of a fully formed and likely decaying runaway beam that is about to scrape off against the first wall, which can be accelerated by vertical displacement events. This scheme thus serves as a last line of defense against potential wall damage. The choice here is locally released tungsten particulates, similar to the previous idea of using tungsten particulates as a dust shield for the divertor. The most significant finding of is that strong pitch angle scattering, in addition to energy attenuation and absorption, of high-energy runaways by tungsten particulates, provides a transport mechanism by which runaways damage of the first wall can be mitigated by reducing and spreading the runaway wall load. We regret the omission in the original article, and thank the authors of the cited article for bringing this matter to our attention. We also find the possibility of using the injected tungsten rod from to terminate the runaways, as opposed to the cloud of tungsten particulates in intriguing. In that context, it is of interest to mention another possibility of a retractable tungsten metal arm that can be swung out for deployment in runaway termination, which would remove the need for post-mitigation recovery of tungsten debris in injection schemes.

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

Mechanically induced thermal runaway severity analysis of Li-ion batteries and continuous energy release monitoring

The large-scale deployment of Li-ion batteries in stationary energy storage and electrical vehicle applications demands a strong focus on safety, particularly on the thermal runaway risk and severity evaluation. A standardized single-side mechanical indentation test protocol was developed to induce an internal short-circuit (ISC) and evaluate cells' thermal runaway severity at different state of charge (SOC). The observed hazard severity (OHS in five categories) and evaluated scores in this work have a comprehensive consideration of each cell's capacity, initial voltage, SOC, temperature and voltage change, allowing a better evaluation of the cells' thermal runaway potential. This method was applied to about 200 Li-ion batteries in order to build an extensive thermal runaway database covering various SOCs, capacities and chemistries. In this study, we monitored the transitions of stored electrochemical energy and applied mechanical energy into both thermal energy and acoustic emissions (AE). The surface temperature and mechanical failures were monitored by infrared imaging and AE to capture critical events within battery cells throughout the mechanical indentation tests. Furthermore, the initial temperature maps can predict two types of follow-up events: thermal runaway or gradual heat release via conduction. Analyzing each cell's severity, AEs, and leveraging the evolving database offer insights into predicting occurrences of thermal runaway. The test method, thermal runaway severity evaluation and prediction, and the corresponding database provide battery designers, manufacturers, and end-users a clear overview of Li-ion batteries' thermal runaway potential under mechanical abuse, advancing the safety design of Li-ion batteries.

Acoustic emission↗
DOE OSTI2026

Microwave-Assisted Hydrogen Generation from Hydrocarbon-Bearing Reservoir Rocks: Stage-Dependent Thermal Runaway and In-Situ Carbonate Engineering

Microwave-assisted hydrogen generation from hydrocarbon-bearing reservoir rocks is strongly influenced by mineralogy, methane activation, carbonate reactions, and thermal runaway behavior. This study investigates a new approach in which carbonate phases are generated in-situ through the reaction of internally produced CO2 with Ca(OH)2 under microwave heating conditions. The objective is to evaluate how rock mineralogy, methane injection, and Ca(OH)2 addition influence hydrogen generation, carbon redistribution, and stage-dependent reaction pathways during microwave exposure. Microwave heating experiments were conducted using Permian Basin reservoir rocks under three experimental conditions: rock-only experiments under Ar atmosphere, CH4–Ar experiments without additive, and CH4–Ar experiments containing 5 wt% Ca(OH)2. Methane-assisted experiments were performed under continuous injection of 30 standard cubic centimeters per minute (sccm) CH4 and 30 sccm Ar. Based on thermal runaway behavior, each experiment was divided into three operational stages: Before Thermal Runaway (BR), After Thermal Runaway–Decrease in Microwave Power (ARD), and After Thermal Runaway– Increase in Microwave Power (ARI). Temperature and gas composition were continuously monitored throughout the experiments. The rock-only experiments demonstrated that hydrogen generation can occur intrinsically from hydrocarbon-bearing rocks under microwave heating, even without externally injected methane. However, hydrogen production did not correlate solely with kerogen content, indicating that mineralogy strongly influences hydrogen-generation pathways. Correlation analyses suggested that kerogen decomposition initially generated CH4, CO, and CO2, followed by secondary hydrocarbon reactions associated with H2 and C2 hydrocarbon formation. Methane-assisted experiments substantially increased hydrogen production; however, identical methane injection rates produced significantly different hydrogen yields among the rock samples, confirming that mineralogical composition controls methaneconversion behavior under microwave heating conditions. The addition of Ca(OH)2 significantly altered carbon evolution behavior in a stage-dependent manner. During the BR stage, Ca(OH)2 reduced gas-phase CO2 production, particularly in carbonate-rich rocks, indicating favorable conditions for in-situ carbonation and carbonate deposition prior to extensive thermal decomposition. The suppression of CO2 during BR became more pronounced with increasing carbonate content of the rock system. After thermal runaway, carbonate-containing systems exhibited enhanced hydrogen generation behavior, suggesting that carbonatederived mineral transformations and carbonate-mediated reactions contribute to high-temperature hydrogen-generation pathways. In carbonate-poor rocks, Ca(OH)2 addition enabled simultaneous URTeC 4493775 2 enhancement of hydrogen production and partial suppression of CO2 release during the post-runaway stages. Overall, the results demonstrate that microwave-assisted hydrogen generation is governed by dynamically evolving interactions among kerogen decomposition, methane activation, mineral transformations, carbonate formation/decomposition, and thermal runaway behavior. This work introduces in-situ carbonate engineering with Ca(OH)2 as a strategy for coupling hydrogen generation with partial insitu carbon management under microwave heating conditions.

03 NATURAL GAS↗
DOE OSTI2023

Runaway electron current reconstitution after a nonaxisymmetric magnetohydrodynamic flush

Benign termination of mega-ampere (MA) level runaway current has been convincingly demonstrated in recent JET and DIII-D experiments, establishing it as a leading candidate for runaway mitigation on ITER. This comes in the form of a runaway flush by parallel streaming loss along stochastic magnetic field lines formed by global magnetohydrodynamic instabilities, which are found to correlate with a low-Z injection that purges the high-Z impurities from a post-thermal-quench plasma. Here, in this study, we show the competing physics that govern the postflush reconstitution of the runaway current in an ITER-like reactor where significantly higher current is expected. The trapped “runaways” are found to dominate the seeding for runaway reconstitution, and the incomplete purge of high-Z impurities helps drain the seed but produces a more efficient avalanche, two of which compete to produce a 2–3 MA step in current drop before runaway reconstitution of the plasma current.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
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 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↗
DOE OSTI2025

Predictions of Cell-to-Cell Propagation and Vent Gas Production in the Thermal Runaway of Lithium-Ion Battery Stacks

This work presents the thermal runaway propagation model LIM1TR (Lithium-ion Modeling with 1-D Thermal Runaway) as an efficient tool to predict different cell-to-cell thermal runaway propagation scenarios. Here, we explored the vent gas volume production and reaction duration highlighting the relationship between these parameters and thermal runaway propagation due to convection by the vented gases. Two metrics based on gas production rate and heating rate are utilized as good indicators of the start and end of thermal runaway. LIM1TR results are compared with and validated by experiments from the literature for single-cell and multicell array experiments of 5 Ah and 10 Ah cells. By accounting for intraparticle diffusion of reacting species in the electrodes, we were able to capture the general dynamics of thermal runaway propagation and estimate acceptable reaction durations compared with the experimental values. Simulation results further demonstrated that varying heating modes lead to distinct reaction durations, consistent with experimental observations. Vent gas volume predictions indicate the need to consider both full and partial oxidation of the electrolyte. The outcomes of this work are building blocks for further investigations of module-to-module propagation by vented gases through convective heat transfer.

25 ENERGY STORAGE↗
DOE OSTI2025

A round bobbin critical current measurement and thermal runaway simulation of REBCO coated conductors

A round bobbin test has been widely used to measure critical current Ic of practical superconductors including Nb–Ti, Nb3Sn, and Bi-2212 round strands at 4.2 K and in high magnetic fields but rarely used for rare earth-barium-copper-oxide (REBCO) coated conductor tapes. Here, we applies this method to REBCO tapes and test their Ic at 77 K, self-field, and 4.2 K with a magnetic field of 14 T applied parallel to the tape. The sample carries very high current densities; it is 1.3 m long and has a dense array of voltage taps that allows probing localized thermal runaways and correlation with Ic variations along the length. At 77 K, the Ic values of all sections were measured and found to be rather uniform. However, at 4.2 K, we found that a 3 cm section has 90% of the overall voltage drop and it went into thermal runaways upon further increasing current. Thus, at 4.2 K, we were prevented by localized thermal runaways and unable to determine Ic of most of the REBCO conductor sections; this behavior is different from Bi-2212 and Bi-2223 for which Ic of all turns of a round barrel sample can be determined. We simulate REBCO tapes during current ramping to a thermal runaway and analyze impacts of various factors including cooling, conductor n-values, and thickness of stabilizer on the thermal runaway current Iq. The combined experiments and simulations identify n-value as a crucial factor that leads to the observed differences between Bi-2212 and REBCO. Together with cooling and the standard deviation of Ic variations along the length, the n-value determines whether Iq is above or below the average Ic of a long piece of conductor and their ratio Iq/Ic,ave, greatly affecting quench detection and thermal runaway prevention strategies of Bi-2212 and REBCO magnets.

Mosat, Marek↗
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

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 OSTI2024

Runaway Gravitational Production of Dark Photons

We demonstrate that gravitational particle production of a massive, Abelian, vector (Proca) field during inflation in the presence of nonminimal coupling to gravity may suffer from an instability which leads to runaway production of high-momentum modes. This is untenable unless there is some mechanism to regulate the runaway. We discuss the parameter space of the particle mass and nonminimal couplings where such a runaway occurs and possible ways to tame the runaway. We find that there is no obvious way to resolve the runaway in a UV completion or with kinetic mixing to the standard model.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗