Search NASA⌕ Search

SEARCH · Search NASA

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Runaway electron deconfinement in SPARC and DIII-D by a passive 3D coil

The operation of a 3D coil—passively driven by the current quench (CQ) loop voltage—for the deconfinement of runaway electrons (REs) is modeled for disruption scenarios in the SPARC and DIII-D tokamaks. Nonlinear magnetohydrodynamic (MHD) modeling is carried out with the NIMROD code including time-dependent magnetic field boundary conditions to simulate the effect of the coil. Further modeling in some cases uses the ASCOT5 code to calculate advection and diffusion coefficients for REs based on the NIMROD-calculated fields, and the DREAM code to compute the runaway evolution in the presence of these transport coefficients. Compared with similar modeling in Tinguely et al (2021 Nucl. Fusion 61 124003), considerably more conservative assumptions are made with the ASCOT5 results, zeroing low levels of transport, particularly in regions in which closed flux surfaces have reformed. Of three coil geometries considered in SPARC, only the n = 1 coil is found to have sufficient resonant components to suppress the runaway current growth. Without the new conservative transport assumptions, full suppression of the RE current is maintained when the thermal quench MHD is included in the simulation or when the RE current is limited to 250kA, but when transport in closed flux regions is fully suppressed, these scenarios allow RE beams on the order of 1–2 MA to appear. Additional modeling is performed to consider the effects of the close ideal wall. In DIII-D, the CQ is modeled for both limited and diverted equilibrium shapes. In the limited shape, the onset of stochasticity is found to be insensitive to the coil current amplitude and governed largely by the evolution of the safety-factor profile. In both devices, prediction of the q-profile evolution is seen to be critical to predicting the later time effects of the coil.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

The SN Ia runaway LP 398-9: detection of circumstellar material and surface rotation

ABSTRACT A promising progenitor scenario for Type Ia supernovae (SNeIa) is the thermonuclear detonation of a white dwarf in a close binary system with another white dwarf. After the primary star explodes, the surviving donor can be spontaneously released as a hypervelocity runaway. One such runaway donor candidate is LP 398-9, whose orbital trajectory traces back ≈105 yr to a known supernova remnant. Here, we report the discovery of carbon-rich circumstellar material around LP 398-9, revealed by a strong infrared excess and analysed with follow-up spectroscopy. The circumstellar material is most plausibly composed of inflated layers from the star itself, mechanically and radioactively heated by the past companion’s supernova. We also detect a 15.4 h periodic signal in the UV and optical light curves of LP 398-9, which we interpret as surface rotation. The rotation rate is consistent with theoretical predictions from this supernova mechanism, and the brightness variations could originate from surface inhomogeneity deposited by the supernova itself. Our observations strengthen the case for this double-degenerate SNIa progenitor channel, and motivate the search for more runaway SNIa donors.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Utilizing the deuterium-tritium fusion resonance to diagnose thermal runaway in igniting plasmas

For high-efficiency inertial confinement fusion implosions, it is predicted that a burning hot spot will successfully encompass all surrounding fuel and then transition into a thermal runaway where the internal energy increase from fusion occurs on a timescale faster than the expansion of the fuel is able to quench the fusion chain reaction after ignition occurs. Observation of this dynamic phase transition would indicate distinct burn properties and indicate an implosion's robustness. A technique for diagnosing the presence of thermal runaway from measurements of nuclear reaction history is presented. The technique is based on taking the logarithmic derivative of the nuclear reaction history, called the 𝛼 curve, and allowing a mathematical decoupling of the mass, volume, and thermal reactivity in the fusion reaction rate equation. During thermal runaway, where the thermal temperature dominates the burn dynamics, a maximum in the 𝛼 curve is found where there is a maximum in the first derivative of the thermal fusion reactivity, an effect to the deuterium-tritium (DT) fusion cross-section resonance. This provides a distinct signature related to the fundamental nature of the DT fusion nuclear resonance and signifies the transition into the fusion thermal instability. Impacts of charged particle transport on the effect are also assessed and the analytical formulas are compared and found to be in agreement with radiation hydrodynamic codes.

high-energy-density plasmas↗
DOE OSTI2022

Thermal Runaway of Li-Ion Cells: How Internal Dynamics, Mass Ejection, and Heat Vary with Cell Geometry and Abuse Type

Thermal runaway of lithium-ion batteries can involve various types of failure mechanisms each with their own unique characteristics. Using fractional thermal runaway calorimetry and high-speed radiography, the response of three different geometries of cylindrical cell (18650, 21700, and D-cell) to different abuse mechanisms (thermal, internal short circuiting, and nail penetration) are quantified and statistically examined. Correlations between the geometry of cells and their thermal behavior are identified, such as increasing heat output per amp-hour (kJ Ah −1 ) of cells with increasing cell diameter during nail penetration. High-speed radiography reveals that the rate of thermal runaway propagation within cells is generally highest for nail penetration where there is a relative increase in rate of propagation with increasing diameter, compared to thermal or internal short-circuiting abuse. For a given cell model tested under the same conditions, a distribution of heat output is observed with a trend of increasing heat output with increased mass ejection. Finally, internal temperature measurements using thermocouples embedded in the penetrating nail are shown to be unreliable thus demonstrating the need for care when using thermocouples where the temperature is rapidly changing. All data used in this manuscript are open access through the NREL and NASA Battery Failure Databank.

25 ENERGY STORAGE↗
DOE OSTI2020

Gas Generation Models for Thermal Runaway in Lithium Ion Batteries

A numerical model is implemented to study cell venting, internal pressure, and gas-phase dynamics behavior of Li-ion cells undergoing thermal runaway. A k-ɛ Reynolds-Averaged Navier-Stokes (RANS) model is adopted to describe the turbulent flow out of the cells, while the fluid dynamics inside the cells is described by Darcy-Forchheimer’s equation. Thermal abuse reactions and gas generation kinetics are described by a single-step lumped reaction model. A series of computational fluid dynamics (CFD) simulations are conducted on a single cell at various states-of-charge (100%, 50%, 25%) to demonstrate detailed flow and thermal behavior as a function of quantity of gas generated during cell venting. Venting events are categorized into two stages: i) breaching of the cell container and ii) thermal runaway reactions. It is shown that the cell response is dominated by the second stage since most of the gases are generated during thermal runaway. Also, the propensity for propagation is highly affected by state-of-charge (SOC). Cells at higher SOCs produce more heat and gas during the venting event, owing to higher mass and concentrations of reacting gases, and consequently reach higher internal cell pressures which increase the risk of side-wall breaching.

Kim, Jinyong↗
DOE OSTI2021

A laser inverse Compton scattering diagnostic to study runaway electron dynamics during tokamak disruptions

The goal is to design a new diagnostic based on laser inverse Compton scattering (LICS) to study the dynamics of runaway electron formation and their subsequent loss during triggered disruptions in DIII-D. It is known that better measurements of the runaway electron population, throughout the evolution of a disruption, is critical for mitigating the impact of energetic electron beams especially in high plasma current discharges such as those in reactor relevant burning plasmas. The goal of the LICS approach being designed for DIII-D is to measure the timeresolved spatial profile and energy distribution of the runaway electrons while they are in the core of the tokamak plasma. During the last year, a code was developed to simulate the signals from a new laser scattering diagnostic for DIII-D, which will measure the relativistic electron population expected in DIII-D after disruptions triggered by argon pellet injection.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

LIM1TR: Lithium-Ion Modeling with 1-D Thermal Runaway (V.1.0)

LIM1TR (Lithium-Ion Modeling with 1-D Thermal Runaway) is an open-source code that uses the finite volume method to simulate heat transfer and chemical kinetics on a quasi 1-D domain. The target application of this software is to simulate thermal runaway in systems of lithium-ion batteries. The source code for LIM1TR can be found at https://github.com/ajkur/lim1tr. This user guide details the steps required to create and run simulations with LIM1TR starting with setting up the Python environment, generating an input file, and running a simulation. Additional details are provided on the output of LIM1TR as well as extending the code with custom reaction models. This user guide concludes with simple example analyses of common battery thermal runaway scenarios. The corresponding input files and processing scripts can be found in the “Examples” folder in the on-line repository, with select input files included in the appendix of this document.

25 ENERGY STORAGE↗
DOE OSTI2020

A fully implicit, scalable, conservative nonlinear relativistic Fokker–Planck 0D-2P solver for runaway electrons

Upon application of a sufficiently strong electric field, electrons break away from thermal equilibrium and approach relativistic speeds. These highly energetic ‘runaway’ electrons (~ MeV) play a significant role in tokamak disruption physics, and therefore their accurate understanding is essential to develop reliable mitigation strategies. As such, we have developed a fully implicit solver for the 0D-2P (i.e., including two momenta coordinates) relativistic nonlinear Fokker–Planck equation (rFP). As in earlier implicit rFP studies (NORSE, CQL3D), electron–ion interactions are modeled using the Lorentz operator, and synchrotron damping using the Abraham–Lorentz–Dirac reaction term. However, our implementation improves on these earlier studies by (1) ensuring exact conservation properties for electron collisions, (2) strictly preserving positivity, and (3) being scalable algorithmically and in parallel. Key to our proposed approach is an efficient multigrid preconditioner for the linearized rFP equation, a multigrid elliptic solver for the Braams–Karney potentials, and a novel adaptive technique to determine the associated boundary values. We verify the accuracy and efficiency of the proposed scheme with numerical results ranging from small electric-field electrical conductivity measurements to the accurate reproduction of runaway tail dynamics when strong electric fields are applied.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Mechanically induced thermal runaway severity analysis for Li-ion batteries

Thermal runaway is one of the most important safety concerns in the deployment of Li-ion batteries. For this study, we developed a standardized single-side indentation test protocol to induce an internal short-circuit. Cell voltage, temperature, and applied compressive force are monitored as a function of time. Each cell is given an observed hazard severity (OHS in five categories) modified from the EUCAR table. Meanwhile a calculated hazard severity (CHS) is obtained solely based on temperature and voltage curves. The calculation formula considers the cell temperature, rate of temperature increase, cell capacity, state of charge (SOC), voltage drop, and voltage drop rate. Each term is assigned with an appropriate weight to scale the calculated score from 0 to100. This method was applied to >100 Li-ion batteries with various SOCs and chemistries, and their OHS and CHS were displayed against SOC. The purpose is to provide battery designers, manufacturers, and end-users a clear comparison of thermal runaway severity of different batteries.

25 ENERGY STORAGE↗
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 OSTI2026

Collision operator for electron runaway in cold weakly ionised plasmas

In cold weakly ionised plasmas, the Dreicer generation mechanism can be non-diffusive as demonstrated in a recent study. By expanding the previous letter, we present the detailed description of a proper collision operator to precisely account for the non-diffusive electron kinetics. The operator appropriately combines the Fokker–Planck operator and Boltzmann operator where free-bound collision cross-sections are valid in the low-energy region. The proposed operator is envisaged to predict runaway electron generation in cold weakly ionised plasmas, particularly to design a runaway-free reactor tokamak start-up.

fusion plasma↗
DOE OSTI2022

Quantitative spatiotemporal mapping of thermal runaway propagation rates in lithium-ion cells using cross-correlated Gabor filtering

Abuse testing of lithium-ion batteries is widely performed in order to develop new safety standards and strategies. However, testing methodologies are not standardised across the research community, especially with failure mechanisms being inherently difficult to reproduce. High-speed X-ray radiography is proven to be a valuable tool to capture events occurring during cell failure, but the observations made remain largely qualitative. We have therefore developed a robust image processing toolbox that can quantify, for the first time, the rate of propagation of battery failure mechanisms revealed by high-speed X-ray radiography. Using Gabor filter, the toolbox selectively tracks the electrode structure at the onset of failure. This facilitated the estimation of the displacement of electrodes undergoing abuse via nail penetration, and also the tracking of objects, such as the nail, as it propagates through a cell. Further, by cross-correlating the Gabor signals, we have produced practical, illustrative spatiotemporal maps of the failure events. From these, we can quantify the propagation rates of electrode displacement prior to the onset of thermal runaway. The highest recorded acceleration (≈514 mm s -2 ) was when a nail penetrated a cell radially (perpendicular to the electrodes) as opposed to axially (parallel to the electrodes). The initiation of thermal runaway was also resolved in combination with electrode displacement, which occurred at a lower acceleration (≈108 mm s -2 ). Our assistive toolbox can also be used to study other types of failure mechanisms, extracting otherwise unattainable kinetic data. Ultimately, this tool can be used to not only validate existing theoretical mechanical models, but also standardise battery failure testing procedures.

25 ENERGY STORAGE↗
DOE OSTI2020

Runaway electron transport in stochastic toroidal magnetic fields

Here, we study the transport and confinement properties of runaway electrons (RE) in the presence of magnetic fields with perturbations producing different levels of stochasticity. We use Kinetic Orbit Runaway Electron Code (KORC) [Carbajal et al. , Phys. Plasmas 24, 042512 (2017) and del-Castillo-Negrete et al. , Phys. Plasmas 25, 056104 (2018)] for simulating the full-orbit (FO) and guiding-center (GC) dynamics of RE in perturbed magnetic fields that exhibit magnetic islands. We extend previous works on this problem [Wingen et al. , Nucl. Fusion 46, 941 (2006); Izzo et al. , Nucl. Fusion 51, 063032 (2011); Papp et al. , Nucl. Fusion 51, 043004 (2011); V. Izzo and P. Parks, Phys. Plasmas 24, 060705 (2017); and Sommariva et al. , Nucl. Fusion 58, 016043 (2018)] by studying in detail full-orbit effects on the RE dynamics. We quantify FO effects on RE transport by performing one-to-one comparisons between FO and GC simulations. It is found that, for the magnetic field configurations considered, GC simulations predict twice the RE losses of FO simulations for 1 MeV and four times the RE losses of FO simulations for 25 MeV. Similarly, we show how different GC and FO dynamics of RE moving around magnetic islands can be, especially in the scenario where the RE Larmor radius is on the order of the size of the magnetic island. We also study the role of rotation of the magnetic islands on RE confinement, and we find that low-frequency toroidal rotation has no observable effect on RE transport in the cases considered. These results shed some light into the potential of avoidance or mitigation mechanisms based on magnetic perturbations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Toroidal modeling of runaway electron loss due to 3-D fields in DIII-D and COMPASS

The 3-D field induced relativistic runaway electron (RE) loss has been simulated for DIII-D, COMPASS and ITER plasmas, utilizing the MARS-F code incorporated with the recently developed and updated RE orbit module (REORBIT). Modeling shows effectively 100% loss of a post-disruption, high-current runaway beam in DIII-D, due to the 1 kG level of magnetic field perturbation produced by a fast growing n = 1 resistive kink instability. The RE loss is shown to be independent of the particle energy or the initial location of particles in the configuration space. Applied resonant magnetic perturbation (RMP) fields from in-vessel coils are not effective for RE beam mitigation in DIII-D, but do produce finite (>10%) RE loss in COMPASS, consistent with experimental observations in above two devices. The major reasons for this difference in RE control by RMP between these two devices are (i) the coil proximity to the RE beam and (ii) the effective coil current scaling versus the machine size and the toroidal magnetic field. About 10% RE loss fraction is also predicted for an ITER 15 MA scenario with pre-disruption plasma, highlighting the role of the plasma response. Up to 30% loss is computed, however, by artificially scaling the equilibrium pressure to zero. This is due to the more resistive plasma response and stronger resulting field line stochasticity. Distributions of the lost REs to the limiting surface show poloidally peaked profile near the high-field-side in both DIII-D and COMPASS, covering about 100o poloidal angle. Higher perturbation field level and/or higher particle energy also result in REs being lost to the low-field-side of the limiting surface of these two devices, increasing the effective wetted area. Finally, for the modeled ITER plasmas, the applied RMP field, with optimal poloidal spectrum that maximizes the edge localized mode control, leads to REs being lost to the lower divertor region of the limiting surface within a narrow poloidal band.

3-D field↗
DOE OSTI2021

Simulation of decelerating streamers in inhomogeneous atmosphere with implications for runaway electron generation

The dynamics of positive and negative streamers is numerically simulated in atmospheric pressure air in the range of parameters corresponding to the streamer deceleration and termination in the middle of the discharge gap. A detailed comparison with experiments in air at constant and variable density demonstrates good agreement between the 2D simulation results and the observations. It is shown that positive and negative streamers behave in radically different ways when decelerating and stopping. When the head potential drops, the negative streamer transits to the mode in which the propagation is due to the forward electron drift. In this case, the radius of the ionization wave front increases, whereas the electric field at the streamer head decreases further and the streamer stops. Its head diameter continues to increase due to the slow drift of free electrons in the residual under-breakdown field. On the contrary, the only advancement mechanism for a positive streamer with a decreasing head potential is a decrease in the effective radius of the ionization wave, leading to a local increase in the electric field. This mechanism makes it possible to compensate for the reduction in the efficiency of gas photoionization at small head diameters. A qualitative 1D model is suggested to describe streamer deceleration and stopping for different discharge polarities. Estimates show that, during positive streamer stopping, the local electric field at the streamer head can exceed the threshold corresponding to the transition of electrons to the runaway mode when the head potential (relative to the surrounding space) decreases to~1.2 kV in atmospheric pressure air. In this case, pulsed generation of a beam of runaway electrons directed into the channel of a stopping positive streamer can occur. The energy of the formed pulsed electron beam depends on the intensity of photoionization in front of the streamer head. This energy can vary from 700 V (when increasing the photoionization rate by a factor of 10 with respect to the value in atmospheric pressure air) to 2.6 kV (when decreasing the photoionization rate by a factor of 1000). It is possible that this behavior of decelerating positive streamers can explain the observed bursts of x-ray radiation during the streamer propagation in long air gaps.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Polarized imaging of visible synchrotron emission from runaway electron plateaus in DIII-D

A new approach to estimate the pitch angle of highly energetic postdisruption runaway electrons is presented by imaging the polarized runaway electron synchrotron emission. The vertical to horizontal polarization ratio is found to range between 3 and 14 in image-average studies looking at repeated shots as well as the polarization in a single discharge. Numerical studies with guiding center and full orbit simulations performed with the SOFT and KORC codes are given and compared to the experiment. The pitch angle evolution analyzed using a 1D impurity diffusion model in combination with kinetic simulations is also presented. Furthermore, these studies find agreement in the magnitude of the pitch angle, but not in the temporal rate of change of pitch angle. The measured pitch angle evolution timescale of order 50 ms falls between a modeled timescale of 20 ms in an equilibrium model and >>100 ms in a nonequilibrium model.

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

Green Additive Limits Runaway Crystallinity in PM6:Y6 Organic Solar Cells but Causes Field-Independent Geminate Recombination

Solvent additives are well-established for improving device performance of organic solar cells by enhancing nanodomain phase separation and crystallization. However, optimum additive concentrations are typically at trace amounts with slight fluctuations causing runaway crystallinity that short the device and limit scale-up potential. Furthermore, we show that a nonhalogenated solvent additive 1-phenylnaphthalene (PN) can be used to optimize high-performance PM6:Y6 devices, but without runaway crystallinity even at concentrations 10× above the optimum. Our holistic analysis of device losses shows that PN increases photoabsorption while decreasing bimolecular recombination. However, above the optimum, charge generation is reduced due to field-independent geminate recombination of charge transfer states that overcomes the other benefits. Rotated Y6 crystallite orientation is a likely culprit, where the device electric field is ineffective at separating the charges. Despite this, device performance remains higher than that of traditional additives, indicating a route to avoiding the high sensitivity of performance to processing conditions.

14 SOLAR ENERGY↗