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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Passive deconfinement of runaway electrons using an in-vessel helical coil

A helical coil designed to passively generate non-axisymmetric fields during a plasma disruption is shown (via electromagnetic analysis, linear MHD modeling, and relativistic drift orbit tracing) to be effective at deconfining runaway electrons (REs) on a time scale significantly faster than the plasma current quench. Magnetic equilibria from DIII-D RE-producing scenarios are used to calculate the toroidal electric field generated during the current quench phase of a disruption, which in turn drives current in the proposed n = 1 in-vessel helical coil, without the need for any external power supplies or disruption detection or prediction techniques. Simulations of the plasma evolution using the TokSys GS Evolve code predict the inductive coupling of coil currents up to 12% of the pre-disruption plasma current into the helical coil. The coil geometry is parametrically varied to maximize both the non-resonant and resonant components of the 3D magnetic perturbation, resulting in δB/B ≈ 10 –2 and a vacuumisland overlapwidth of up to 0.7ψ N . The REORBIT module of the MARS-F code is used to model the full non-axisymmetric magnetic field and trace RE drift orbits to determine the effect on RE deconfinement, with up to 70% of the RE orbits lost after 0.2 ms. A two-stage evolution of the RE orbit loss fraction is observed to be caused by resonant trapping between multiple magnetic island chains. Finally, electromagnetic and thermal stresses on the coil are calculated to be within operational limits for installation in DIII-D, and scale favorably to a reactor-size device. Furthermore, these findings motivate future experimental study of the helical coil concept in DIII-D or other tokamaks.

3D coil↗
DOE OSTI2021

Modeling the complete prevention of disruption-generated runaway electron beam formation with a passive 3D coil in SPARC

Here, the potential formation of multi-mega-ampere beams of relativistic 'runaway' electrons (REs) during sudden terminations of tokamak plasmas poses a significant challenge to the tokamak's development as a fusion energy source. Here, we use state-of-the-art modeling of disruption magnetohydrodynamics coupled with a self-consistent evolution of RE generation and transport to show that a non-axisymmetric in-vessel coil will passively prevent RE beam formation during disruptions in the SPARC tokamak, a compact, high-field, high-current device capable of achieving a fusion gain Q > 2 in deuterium–tritium plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Toroidal modeling of runaway electron loss due to 3D fields in ITER

Mitigation of runaway electrons (REs) by three-dimensional (3D) magnetic field perturbations is numerically investigated for the ITER 15 MA baseline D–T scenario, utilizing the MARS-F code (Liu et al Phys. Plasmas 7 3681) with a drift orbit test particle tracing module. Considered are two types of 3D fields: the n = 3 (n is the toroidal mode number) resonant magnetic perturbation (RMP) utilized for the purpose of controlling the edge localized modes in ITER, and perturbations generated by the n = 1 magneto-hydrodynamic (MHD) instabilities in a post-disruption plasma. The RMP field, applied to a pre-disruption plasma, is found to be moderately effective in mitigating the RE seeds in ITER when vacuum field model is assumed. Up to ~40% loss fraction is possible at 90 kA-turn coil current. The mitigation efficiency is however substantially reduced, down to less than 5%, when the plasma response is taken into account. This is due to strong screening of the resonant magnetic field components by the plasma response resulting in much less field line stochasticity. On the other hand, the MARS-F modeling, based on the DINA-simulated post-disruption equilibria, shows that the n = 1 resistive kink instabilities develop in these plasmas, as the edge safety factor qa evolves and drops below integer numbers. RE mitigation by these MHD instabilities is sensitive to the eigenmode structure. The best mitigation is achieved as q a drops below 3, when a global kink instability occurs that encompasses both internal and external components. This global instability is found to be capable of mitigating over 80% MeV-level passing RE orbits at a field perturbation |δB|/B 0 that is comparable to that observed in DIII-D experiments, and full mitigation if the perturbation amplitude is doubled. The 'wetted' area on the ITER limiting surface, due to MHD instability induced RE loss, generally increases with the perturbation amplitude (together with increasing loss fraction). At the highest perturbation level assumed in this study, the wetted area reaches ~60% of the total limiting surface area. Lastly, the lost RE orbits mainly strike the outer divertor region of the limiting surface, with some fraction also hitting a wide area along the inboard side of the surface.

3D fields↗
DOE OSTI2022

Simulation of runaway electron production with CQL3D coupled to NIMROD

Abstract A coupling between two distinctly different codes—one magnetohydrodynamic (MHD) and another kinetic—is achieved and applied for simulation of runaway electron (RE) production. The 3D initial value MHD code NIMROD simulates a DIII-D pure neon shattered pellet injection plasma quench including the propagation and ablation of the fragments, ionization and recombination of the impurities, and the radiated and transported energies. The field data from NIMROD is then used by the bounce-averaged Fokker–Planck Collisional QuasiLinear 3D (CQL3D) kinetic code to simulate the production of REs and their radial transport. The coupling procedure involves mapping of data between different grids and adjustment of the NIMROD toroidal electric field when REs appear. It is shown that without the radial transport, a large RE current is generated, up to 30% of the pre-pellet ohmic current. However, when the radial transport is included in CQL3D, the RE current is reduced to undetectable level, consistent with experiment. Various forms of the radial diffusion are surveyed to determine conditions when the fast electrons would not have time to be accelerated to relativistic energies before they are lost to chamber wall.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

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

Abstract Experimental trends in thermal plasma partial recombination resulting from massive D 2 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 ( n e ≈ 10 18 m − 3 ) to increase RE final loss MHD levels. In both DIII-D and JET, thermal electron density n e is found to drop by ∼100 × when the thermal plasma partially recombines, with a minimum at a vacuum vessel-averaged D 2 density in the range 10 20 − 10 21 m − 3 . RE effective resistivity also drops after partial recombination, indicating expulsion of the Ar content. The n e 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 n e values on time scales faster than expected RE plateau vertical drift timescales (of order 100 ms), provided sufficient D 2 or H 2 is injected. In SPARC, it is predicted that achieving significant n e recombination will be challenging, due to the very high RE current density. In both ITER and SPARC, it is predicted that achieving low n e will be easier with Ar as a background impurity (rather than Ne).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Parametric study of Alfvénic instabilities driven by runaway electrons during the current quench in DIII-D

Abstract To avoid or mitigate runaway electron (RE) beams in ITER, RE-driven instabilities are actively studied as a complimentary technique to massive material injection. In this work we report experimental dependencies of Alfvénic instabilities driven by REs during the current quench in DIII-D on plasma and material injection parameters. These instabilities, observed in the frequency range of 0.1–3 MHz, correlate with increased RE loss and thus may play a role in non-sustained RE beams. It was found that as the toroidal magnetic field ( B T ) decreases, the RE population becomes more energetic, the energy of instabilities increases, and no RE beam is observed when the maximum energy of REs exceeds 15 MeV (or when B T is below 1.8 T). Analysis of disruptions at plasma core temperature ( T e ) of 1 keV and 8 keV shows that the RE population is much less energetic (with the maximum energy of only about 3 MeV) when T e is high, and no instabilities are observed in this case. Besides disruptions above caused by Ar injection, cases with Ne and D 2 injections were also studied. Both Ne and D 2 injections cause no sustained RE beams, however, for different reasons. Measurements of the instability polarization indicate that it is of predominantly compressional nature at the edge, which is consistent with modeling suggesting excitation of compressional Alfvén eigenmodes. However, drive of global Alfvén eigenmodes is also possible at low frequencies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Cross-machine comparison of runaway electron generation during tokamak start-up for extrapolation to ITER

A cross-machine comparison of global parameters that determine the runaway electron (RE) generation and loss process during tokamak start-up was carried out with the aim to extrapolate these to ITER. The study found that all considered discharges, also those that do not show signs of RE, are non-thermal at the start, i.e. have a streaming parameter larger than 0.1. During the current ramp-up the electric field, E, remains above the critical value, E c , that allows RE in the plasma. The distinction to be made is not if RE can form but, if sufficient RE can form fast enough such that they are detected or start to dominate the dynamics of the tokamak discharge. The dynamics of the value of E, density and temperature during tokamak are key to the formation of RE. It was found that larger devices operate with E closer to E c , due to their higher temperatures, hence the RE generation is relatively slower. The slower time scales for the formation of RE, estimated to be of the order of 100s of ms in ITER simplifies the development of avoidance schemes. The RE confinement time is also an important determinant of the entire process and is found to increase with the device size. The study also revealed that drift orbit losses, a mechanism often attributed as the main RE loss mechanism during the early tokamak discharge, are actually more difficult to achieve. RE losses might be more likely attributed to RE diffusion due to magnetic turbulence.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Runaway electron plateau current profile reconstruction from synchrotron imaging and Ar-II line polarization angle measurements in DIII-D

Abstract Current profile reconstructions are obtained for high current ( I p ≃ 550 kA) post-disruption runaway electron (RE) plateau plasmas in DIII-D. Two novel methods of measuring the RE current profile in high-current RE plateaus are introduced and compared: localization of the q = 2 rational surface using visible synchrotron emission (SE) imaging and the measurement of the polarization angle of line-integrated Ar-II line emission. The two methods are found to be consistent with each other within the data uncertainties. Different simulations of the RE current profile are compared with the measurements: the toroidal fluid RE model is found to best fit the data, within the measurement uncertainties. In addition to introducing two novel methods to measure the RE current profile and validating present simulation capabilities, this work demonstrates that instabilities can grow at q = 2 and q = 1 surfaces without necessarily causing a RE final loss instability. Numerical simulations are also presented to elucidate the role of these instabilities on synchrotron emission.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Modelling the brittle failure of graphite induced by the controlled impact of runaway electrons in DIII-D

The thermo-mechanical response of an ATJ graphite sample to controlled runaway electron (RE) dissipation, realized in DIII-D, is modelled with a novel work-flow that features the RE orbit code KORC, the Monte Carlo particle transport code Geant4 and the finite element multiphysics software COMSOL. KORC provides the RE striking positions and momenta, Geant4 calculates the volumetric energy deposition and COMSOL simulates the thermoelastic response. Brittle failure is predicted according to the maximum normal stress criterion, which is suitable for ATJ graphite owing to its linear elastic behavior up to fracture and its isotropic mechanical properties. Measurements of the conducted energy, damage topology, explosion timing and blown-off material volume, impose a number of empirical constraints that suffice to distinguish between different RE impact scenarios and to identify RE parameters which provide the best match to the observations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Toroidal modeling of Alfvén eigenmodes excited by runaway electrons in DIII-D and ITER

The non-perturbative MHD-kinetic hybrid code MARS-K (Liu et al 2014 Phys. Plasmas 21 056105) is updated to include relativistic effects for kinetic fast particles, enabling the code to model excitation of Alfvén eigenmodes (AEs) by runaway electrons (REs) in post-disruption tokamak plasmas. Applying the updated code to RE beams in both DIII-D and ITER, a zoo of AE modes triggered by trapped REs due to precessional drift-kinetic resonances is computed while scanning the RE energy. At fixed RE energy, multiple unstable roots are also excited. These AE modes possess radially different eigenmode structures, ranging from global modes to core-localized ones. The computed mode frequency is in the Alfvén frequency range, increasing with the assumed RE energy in a staircase fashion and quantitatively matching the experimental measurement (in DIII-D). At the (more relevant) high-frequency range (above 1 MHz), the modeled eigenmodes are identified as compressional AEs (CAEs) in DIII-D and a mixture of CAE and shear Alfvén waves in ITER.

Alfvén eigenmodes↗
DOE OSTI2026

Variability of MHD instabilities in benign termination of high-current runaway electron beams in the JET and DIII-D tokamaks

Benign termination, in which magnetohydrodynamic (MHD) instabilities deconfine runaway electrons (REs) following hydrogenic injections, is a promising strategy for mitigating dangerous RE loads after disruptions. Recent experiments on the Joint European Torus (JET) have explored this scenario at higher pre-disruptive plasma currents than are achievable on other devices, revealing challenges in obtaining benign terminations at I p ≥2.5 MA. This work analyzes the evolution of these high-current RE beams and their terminating MHD events using fast magnetic sensor measurements and EFIT equilibrium reconstructions for approximately 40 JET and 20 DIII-D tokamak discharges. On JET, unsuccessful non-benign terminations occur at low edge safety factor (q edge ≈ 2), and are preceded by intermittent, non-terminating MHD events at higher rational qedge. Trends in the internal inductance I i indicate more peaked RE current profiles in the high-I p non-benign population, which may hinder successful recombination through re-ionization of the companion plasma. In contrast, benign terminations on JET typically occur at higher q edge ≥3 and exhibit less peaked RE current profiles. DIII-D displays a broader range of terminating edge safety factors, again correlated with the measured values. Across both tokamaks, the RE current peaking is therefore found to determine which MHD instability boundary is encountered, a result confirmed by linear resistive MHD modeling with the CASTOR3D code. Measured growth rates are similar for benign and non-benign cases, indicating that ideal MHD timescales at low density after hydrogenic injection do not alone explain efficient RE deconfinement. Instead, non-benign cases are most readily characterized by their comparably lower overall MHD perturbation amplitudes δB. These observations suggest that the interplay between ideal and resistive dynamics governs the termination process, with implications for extrapolating benign RE termination to high-I p reactor scenarios.

MHD instabilities↗
DOE OSTI2025

Ultrasonic Characterization of Lithium Ion Thermal Runaway Conditions for Real-time Ultrasonic Enabled BMS Integration

The demand for energy storage is growing, and lithium-ion batteries are a promising technology to meet this need due to their high power/energy density, high round-trip efficiency, rapid response time, and portability. However, recent catastrophic events caused by thermal runaway have slowed their adoption, highlighting the need for an early warning system for battery failure. In this work, ultrasound is used to detect physical changes in 950 mAh batteries by identifying material property changes independent of voltage and current. [Ultrasound signal features (e.g., time of flight, maximum frequency component) were extracted as the batteries were cycled and subjected to both constant current and constant voltage overcharge and were used to develop two metrics identifying failure: a warning to detect the start of overcharge and an emergency stop (E-stop) to immediately take the battery out of service. The identification method involved locating magnitude differences of several ultrasound features compared to baseline operation considering different currents and temperatures, and the warning/stopping metrics were consistent across all experiments. For an average overcharge time of 140 minutes, the average warning was issued 124 minutes before the failure and the average E-stop was triggered 94 minutes before failure. As a test of using ultrasound for early warning detection, a battery was forced into overcharge and returned to normal cycling conditions based on the previously determined warning metric. Both the voltage profile and the ultrasound measurements returned to their baseline behavior, indicating that ultrasonic detection can not only identify battery failure before a catastrophic event, but can also provide early enough warnings such that overcharges can be detected and corrected quickly enough so the battery does not need to be decommissioned.]

25 ENERGY STORAGE↗
NASA NTRS2017

A Thermal Runaway Failure Model for Low-Voltage BME Ceramic Capacitors with Defects

Reliability of base metal electrode (BME) multilayer ceramic capacitors (MLCCs) that until recently were used mostly in commercial applications, have been improved substantially by using new materials and processes. Currently, the inception of intrinsic wear-out failures in high quality capacitors became much greater than the mission duration in most high-reliability applications. However, in capacitors with defects degradation processes might accelerate substantially and cause infant mortality failures. In this work, a physical model that relates the presence of defects to reduction of breakdown voltages and decreasing times to failure has been suggested. The effect of the defect size has been analyzed using a thermal runaway model of failures. Adequacy of highly accelerated life testing (HALT) to predict reliability at normal operating conditions and limitations of voltage acceleration are considered. The applicability of the model to BME capacitors with cracks is discussed and validated experimentally.

reliability↗
NASA NTRS2017

Statistical Characterization of 18650-Format Lithium-Ion Cell Thermal Runaway Energy Distributions

Effective thermal management systems, designed to handle the impacts of thermal runaway (TR) and to prevent cell-to-cell propagation, are key to safe operation of lithium-ion (Li-ion) battery assemblies. Critical factors for optimizing these systems include the total energy released during a single cell TR event and the fraction of the total energy that is released through the cell casing vs. through the ejecta material. A unique calorimeter was utilized to examine the TR behavior of a statistically significant number of 18650-format Li-ion cells with varying manufacturers, chemistries, and capacities. The calorimeter was designed to contain the TR energy in a format conducive to discerning the fractions of energy released through the cell casing vs. through the ejecta material. Other benefits of this calorimeter included the ability to rapidly test of large quantities of cells and the intentional minimization of secondary combustion effects. High energy (270 Wh/kg) and moderate energy (200 Wh/kg) 18650 cells were tested. Some of the cells had an imbedded short circuit (ISC) device installed to aid in the examination of TR mechanisms under more realistic conditions. Other variations included cells with bottom vent (BV) features and cells with thin casings (0.22 1/4m). After combining the data gathered with the calorimeter, a statistical approach was used to examine the probability of certain TR behavior, and the associated energy distributions, as a function of capacity, venting features, cell casing thickness and temperature.

Thermal Runaway Energy Distribution↗
NASA NTRS1995

Disk Tides and Accretion Runaway

It is suggested that tidal interaction of an accreting planetary embryo with the gaseous pre-planetary disk may provide a mechanism to breach the so-called runaway limit during the formation of the giant planet cores.

accretion↗
DOE OSTI2021

An asymptotic-preserving 2D-2P relativistic Drift-Kinetic-Equation solver for runaway electron simulations in axisymmetric tokamaks

We propose an asymptotic-preserving (AP), uniformly convergent numerical scheme for the relativistic collisional Drift-Kinetic Equation (rDKE) to simulate runaway electrons in axisymmetric toroidal magnetic field geometries typical of tokamak devices. The approach is derived from an exact Green's function solution with numerical approximations of quantifiable impact, and results in a simple, two-step operator-split algorithm, consisting of a collisional Eulerian step, and a Lagrangian orbit-integration step with analytically prescribed kernels. The AP character of the approach is demonstrated by analysis of the dominant numerical errors, as well as by numerical experiments. We demonstrate the ability of the algorithm to provide accurate answers regardless of plasma collisionality on a circular axisymmetric tokamak geometry.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Overcharge-to-thermal-runaway behavior and safety assessment of commercial lithium-ion cells with different cathode materials: A comparison study

In this paper, overcharge behaviors and thermal runaway (TR) features of large format lithium-ion (Li-ion) cells with different cathode materials (LiFePO 4 (LFP), Li[Ni 1/3 Co 1/3 Mn 1/3 ]O 2 (NCM111), Li[Ni 0.6 Co 0.2 Mn 0.2 ]O 2 (NCM622) and Li[Ni 0.8 Co 0.1 Mn 0.1 ]O 2 (NCM811)) were investigated. The results showed that, under the same overcharge condition, the TR of LFP Li-ion cell occurred earlier compared with the NCM Li-ion cells, indicating its poor overcharge tolerance and high TR risk. However, when TR occurred, LFP Li-ion cell exhibited lower maximum temperature and mild TR response. All NCM Li-ion cells caught fire or exploded during TR, while the LFP Li-ion cell only released a large amount of smoke without fire. According to the overcharge behaviors and TR features, a safety assessment score system was proposed to evaluate the safety of the cells. In short, NCM Li-ion cells have better performance in energy density and overcharge tolerance (or low TR risk), while LFP Li-ion cell showed less severe response to overcharging (or less TR hazards). For NCM Li-ion cells, as the ratio of nickel in cathode material increases, the thermal stability of the cathode materials becomes poorer, and the TR hazards increase. Such a comparison study on large format Li-ion cells with different cathode materials can provide deeper insights into the overcharge behaviors and TR features, and provide guidance for engineers to reasonably choose battery materials in automotive applications.

25 ENERGY STORAGE↗