SEARCH · Engineering Papers
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.
This article explores the "Emergency Operations Center of the Future", including outcomes from an in-depth landscape assessment, stakeholder outreach, and tabletop exercises that explored how the rapid pace of advancements in science and technology can both accelerate and be highly disruptive to emergency services.
Critical minerals (CM), such as rare earth elements (REE), cobalt, nickel, and lithium, have important uses in modern energy and technologies, yet are vulnerable to potential supply chain disruptions. One potential domestic CM source is fossil energy wastes, such as coal combustion ash, acid mine drainage (AMD) and treatment solids (AMD solids), and Oil and Gas (O&G) drilling wastes (drill cuttings and produced waters). CM recovery from these feedstocks is promising due to their abundant quantity and fast availability as waste products. To develop informed and effective CM recovery, DOE’s National Energy Technology Laboratory (NETL) have collected and analyzed CM data for aforementioned fossil energy wastes, and utilized advanced geochemical characterization (e.g., synchrotron microprobe, sequential extraction and geochemical modeling) to identify the CM speciation and binding environments. Novel methods that recover multiple CMs while co-producing other valuable byproducts from these feedstocks have been developed. Successful examples include: (1) the discovery of easily mobile REE phases in Ca-rich coal combustion ash resulted in a patented REE recovery process from the ash feedstock while producing zeolite sorbents from the extraction wastes; (2) the successful identification of REE/Co/Ni/Zn hosting phases in acid mine drainage treatment solids (AMD solids) has informed the sequential CM recovery from AMD solids and has inspired lithium sorbent development from the extraction wastes; (3) the recovery potential of Li and other CMs in O&G produced waters and drill cuttings has been explored while the extraction residuals have been demonstrated to support plant growth as soil supplements. The innovations driven by characterization information have the potential to maximize CM recovery revenue, offset the cost of waste management and wastewater treatments while reducing the cost and environmental footprint of CM extraction.
Rare earth elements (REE) and other critical minerals (CM, e.g., Co, Ni, Li) have important uses in green energy and modern technologies, yet are vulnerable to potential supply chain disruptions. One potential domestic CM source is fossil energy wastes, such as acid mine drainage (AMD) and treatment solids (AMD solids), coal combustion ash, and Oil and Gas (O&G) drilling wastes (drill cuttings and produced waters). While they can contain lower CM concentrations then traditional ore, the quantity and fast availability as waste feedstock makes them a promising CM resource. To explore their promise, researchers at DOE’s National Energy Technology Laboratory (NETL) have collected and analyzed CM data for aforementioned fossil energy wastes, and utilized advanced geochemical characterization (e.g., synchrotron microprobe and sequential extraction) to identify the CM speciation and binding environments, and developed sustainable and targeted CM recovery. Successful examples include: (1) the discovery of easily mobile REE phases in Ca-rich coal combustion ash and developing a patented REE recovery process from Ca-rich Powder River Basin coal ash; (2) the successful identification of REE/Co/Ni/Zn hosting phases in acid mine drainage treatment solids (AMD solids) with diverse chemical composition (Al, Mn, or Fe-rich) informing the sequential recovery of different REE/CMs from AMD solids; (3) the recovery potential of Li and other CMs in O&G produced waters and drill cuttings. The innovations driven by characterization have the potential to offset the cost of waste management and wastewater treatments while reducing the cost and environmental footprint of CM extraction.
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.
Rare earth elements (REE) and other critical minerals (CM, e.g., Co, Li) have important uses in green energy and modern technologies, yet are vulnerable to potential supply chain disruptions. One potential domestic source of CM is fossil energy wastes, such as acid mine drainage (AMD) and treatment solids (AMD solids), coal combustion ash, and Oil and Gas (O&G) drilling wastes (drill cuttings and produced waters). CM recovery from these waste feedstocks is promising due to their abundances and fast availability as waste products. On the other hand, the CM occurrence in these wastes is in small quantities compared to traditional ore bodies. Thus, novel and strategic separation and extraction processes are under exploration. Researchers at DOE’s National Energy Technology Laboratory (NETL) are collecting and analyzing CM data for aforementioned fossil energy wastes, and utilizing advanced geochemical characterization (e.g., synchrotron microprobe and sequential extraction) to understand the CM speciation and binding environments in these materials to better inform sustainable and effective recovery. Successful examples discussed in this talk include: (1) the discovery of easily mobile REE phases in Ca-rich coal combustion ash and developing a patented REE recovery process from Ca-rich Powder River Basin coal ash; (2) the successful identification of REE/Co/Ni/Zn hosting phases in acid mine drainage treatment solids (AMD solids) with diverse chemical composition (Al, Mn, or Fe-rich) informing the sequential recovery of different REE/CMs from AMD solids; (3) the recovery potential of Li and other CMs in O&G produced waters and drill cuttings. The innovations driven by characterization have the potential to offset the cost of waste management and wastewater treatments while reducing the cost and environmental footprint of CM extraction.
Novel disruption prevention solutions spanning a range of control regimes are being developed and tested on DIII-D to enable ITER success. First, a new real-time control algorithm has been developed and tested for regulating nearness to stability limits and maintaining safety-margins. Its first application has been for reliable prevention of vertical displacement events (VDEs) by adjusting plasma elongation (κ) and the inner-gap between the plasma and inner-wall in response to real- time open-loop VDE growth rate (γ) estimators. VDEs were robustly prevented up to average open-loop growth rates of 800 rad/s with initial tunings, with only applying shape modification when near safety limits. Second, the disruption risk during fast, emergency shutdown after large tearing and locked modes can be significantly improved by transitioning to a limited topology during shutdown. More than 50% of emergency limited shutdowns after locked modes reach a final normalized current I_N < 0.3 before terminating, scaling to the 3 MA ITER requirement. Furthermore, this is in contrast to diverted shutdowns, the majority of which disrupt at I_N > 0.8. Despite improvements, these results highlight the critical importance of early prevention. Third, a novel emergency shut down method has been developed which excites instabilities to form a warm, helical core post-thermal quench. The current quench extends to ~100ms and avoids VDEs and runaway electron generation. Novel real-time machine learning disruption prediction has been integrated with the DIII-D proximity controller, and a real- time compatible multi-mode MHD spectroscopy technique has been developed. Results presented here were enabled by a focused effort, the Disruption Free Protocol, in DIII-D’s 2019-20 campaign to complement disruption prevention experiments with a large piggy-back program. In addition to testing novel techniques, it is estimated to have helped avoid 32 potential disruptions in piggyback operations with rapid, early shutdowns after large rotating n=1 or locked modes.
We report the acceleration of the magnetic island rotation by the modulated resonant magnetic perturbation (MRMP) has been studied in J-TEXT tokamak experiments. After applying the MRMP, the phase difference between the tearing mode (TM) and MRMP, Δξ, oscillated near the effective phase difference, Δξ eff , which was defined as the time averaged value of Δξ. When the Δξ eff was closed to the—π/2, the MRMP only contributed an accelerating torque on the TM. As the result, the TM rotation frequency was increased by a few kilohertz for the optimized relative phase by small RMPs of the order of 10 -5 of the toroidal field and the locked mode induced disruption was avoided. It is found that the TM rotation could be increased to a higher frequency by applying a stronger MRMP. There is a negative sinusoidal relationship between TM frequency and Δξ eff .
Tokamak disruptions are associated with breaking magnetic surfaces, which makes magnetic field lines chaotic in large regions of the plasma. The enforcement of quasi-neutrality in a region of chaotic field lines requires an electric potential that has both short and long correlation distances across the magnetic field lines. The short correlation distances produce a Bohm-like diffusion coefficient ∼Te/eB and the long correlation distances aT produce a large scale flow ∼Te/eBaT. This cross-field diffusion and flow are important for sweeping impurities into the core of a disrupting tokamak. The analysis separates the electric field in a plasma into the sum of a divergence-free, E→B, and a curl-free, E→q, part, a Helmholtz decomposition. The divergence-free part of E→ determines the evolution of the magnetic field. The curl-free part enforces quasi-neutrality, E→q=−∇→Φq. Magnetic helicity evolution gives the required boundary condition for a unique Helmholtz decomposition and an unfortunate constraint on steady-state tokamak maintenance.
The FESAC long range plan calls out disruption avoidance and mitigation as key remaining technical gaps. In discussing the roles of DIII-D and NSTX-U, the FESAC long range plan says “Additional research on these facilities, in combination with private and international collaborations, continuing support of existing university tokamak programs, and utilization of US expertise in theory and simulation, is needed to find solutions to remaining technical gaps. These gaps include disruption prediction, avoidance, and mitigation …”. Disruptions pose an existential threat to ITER and to FPPs. For a fusion reactor, unplanned shutdowns caused by disruptions will be a significant barrier to connecting such a reactor to the electric grid, even if disruption mitigation is successful. Disruption studies for ITER in recent years have largely focused on disruption mitigation (e.g., pellet injection), motivated by near-term deadlines for finalizing the design of the mitigation hardware. It is recognized, however, that mitigation alone will not suffice. The 2022 U.S. ITER Research Needs Workshop Report states that ”[d]isruptions are considered the largest threat to the ITER Research Program”, and that “[m]itigation should be a last resort”. As we discuss below, there are unresolved foundational issues that play a critical role in avoidance, and DIII-D is an ideal device for generating the data needed to address these issues.
Both predicting and controlling disruptions are critical and urgent issues for ITER as some disruptions with a short warning time may be unavoidable. For these cases, a rapid response disruption mitigation system (DMS) referred to as the electromagnetic particle injector (EPI) is being developed. The primary advantages of the EPI are its fast response time and high velocity, which have been demonstrated in offline experiments (Raman et al., 2021). The EPI is capable of accelerating a metallic sabot electromagnetically using a rail gun to the required velocities ( > 2 km/s) within 2 ms. Two high-field racetrack magnets able to generate fields over 2 T are positioned above and below the rails to permit high velocity at low rail currents, a requirement to minimize electrode erosion. At the end of the acceleration phase, a sabot capture mechanism retains the spent sabot inside the vacuum chamber that houses the EPI. At this point, it releases well-defined microspheres, or a shell pellet, of a radiative payload into the disrupting plasma. A remotely operated sabot loading system positioned behind the injector contains several pre-equipped sabots that can be loaded by an operator from the tokamak control room. The injector is interfaced to the tokamak through a guide tube attached to the front of the EPI vacuum chamber. Finally, the advantages of the EPI system over other DMSs under consideration are described in conjunction with the details of an EPI system designed for near-term test on an existing large tokamak.
Fusion-born alpha particles in ITER disruption simulations are investigated as a possible drive of Alfvénic instabilities. The ability of these waves to expel runaway electron (RE) seed particles is explored in the pursuit of a passive, inherent RE mitigation scenario. The spatiotemporal evolution of the alpha particle distribution during the disruption is calculated using the linearized Fokker–Planck solver CODION coupled to a fluid disruption simulation. These simulations are done in the limit of no alpha particle transport during the thermal quench, which can be seen as a most pessimistic situation where there is also no RE seed transport. Under these assumptions, the radial anisotropy of the resulting alpha population provides free energy to drive Alfvénic modes during the quench phase of the disruption. We use the linear gyrokinetic magnetohydrodynamic code LIGKA to calculate the Alfvén spectrum and find that the equilibrium is capable of sustaining a wide range of modes. The self-consistent evolution of the mode amplitudes and the alpha distribution is calculated utilizing the wave-particle interaction tool HAGIS. Intermediate mode number (n = 7–15, 22–26) toroidal Alfvén eigenmodes are shown to saturate at an amplitude of up to δB/B ≈ 0.1% in the spatial regimes crucial for RE seed formation. We find that the mode amplitudes are predicted to be sufficiently large to permit the possibility of significant radial transport of REs.
Here in this paper we lay the groundwork for a robust cross-device comparison of data-driven disruption prediction algorithms on DIII-D and JET tokamaks. In order to consistently carry on a comparative analysis, we define physics-based indicators of disruption precursors based on temperature, density, and radiation profiles that are currently not used in many other machine learning predictors for DIII-D data. These profile-based indicators are shown to well-describe impurity accumulation events in both DIII-D and JET discharges that eventually disrupt. The univariate analysis of the features used as input signals in the data-driven algorithms applied on the data of both tokamaks statistically highlights the differences in the dominant disruption precursors. JET with its ITER-like wall is more prone to impurity accumulation events, while DIII-D is more subject to edge-cooling mechanisms that destabilize dangerous magnetohydrodynamic modes. Even though the analyzed data sets are characterized by such intrinsic differences, we show through a few examples that the inclusion of physics-based disruption markers in data-driven algorithms is a promising path toward the realization of a uniform framework to predict and interpret disruptive scenarios across different tokamaks. As long as the destabilizing precursors are diagnosed in a device-independent way, the knowledge that data-driven algorithms learn on one device can be re-used to explain a disruptive behavior on another device.
Structural analyses of disruptions and the electromagnetic (EM) analyses to support them typically use symmetry to produce a manageable model size. Nonaxisymmetric halo loads require larger analysis models. Here, in this paper, the results of transient EM analyses of halo strikes during a vertical displacement event are presented. The plasma motions and halo characteristics are prescribed based on analyses performed by the Institute of Plasma Physics of the Czech Republic. The time transient EM solution provides loads to a full three-dimensional transient dynamic analysis of the vacuum vessel. The responses to large lateral halo loads are altered and mitigated by the mode responses of the vessel. Dynamic load factors are computed for specific locations in the vessel and supports.
Abstract Precise values for radiated energy in tokamak disruption experiments are needed to validate disruption mitigation techniques for burning plasma tokamaks like ITER and SPARC. Control room analysis of radiated power ( P rad ) on JET assumes axisymmetry, since fitting 3D radiation structures with limited bolometry coverage is an under-determined problem. In mitigated disruptions, radiation is toroidally asymmetric and 3D, due to fast-growing 3D MHD modes and localized impurity sources. To address this problem, Emis3D adopts a physics motivated forward modeling (‘guess and check’) approach, comparing experimental bolometry data to synthetic data from user-defined radiation structures. Synthetic structures are observed with the Cherab modeling framework and a best fit chosen using a reduced χ 2 statistic. 2D tomographic inversion models are tested, as well as helical flux tubes and 3D MHD simulated structures from JOREK. Two nominally identical pure neon shattered pellet injection (SPI) mitigated discharges in JET are analyzed. 2D tomographic inversions with added toroidal freedom are the best fits in the thermal quench (TQ) and current quench (CQ). In the pre-TQ, 2D reconstructions are statistically the best fits, but are likely over-optimized and do not capture the 3D radiation structure seen in fast camera images. The next-best pre-TQ fits are helical structures that extend towards the high-field side, consistent with an impurity flow under the magnetic nozzle effect also observed in JOREK simulations. Whole-disruption radiated fractions of 0.98 + 0.03 / − 0.29 and 1.01 + 0.02 / − 0.17 are found, suggesting that the stored energy may have been fully mitigated by each SPI, although mitigation efficiencies well below ITER and SPARC requirements for high energy pulses are still within the large uncertainties. Emis3D is also used to validate JOREK SPI simulations, and confirms improvements in matching experiment from changes to impurity modeling. Time-dependent toroidal peaking factors are calculated and discussed.
Mitigation of the multiple risks associated with disruptions and runaway electrons in tokamaks involves competing demands. Success requires that each risk be understood sufficiently that appropriate compromises can be made. Here the focus is on the interaction of short timescale magnetic-perturbations with the structure in ITER that is closest to the plasma, blanket modules covered by separated beryllium tiles. The effect of this tiled surface on the perturbations and on the forces on structures is subtle. Indeterminacy can be introduced by tile-to-tile shorting. A determinate subtlety is introduced because electrically separated tiles can act as a conducting surface for magnetic perturbations that have a normal component to the surface. A practical method for including this determinate subtlety into plasma simulations is developed. The shorter the timescales and the greater the localization, particularly in the toroidal direction, the more important the magnetic effects of the tiles become.
Shattered pellet injector systems have been installed on DIII-D, JET, and KSTAR and used to experimentally determine the effectiveness of the shattered pellet injection (SPI) process in mitigating the deleterious effects of a tokamak plasma disruption. Pellets are fired, and before entering the plasma, strike a bent tube known as a shatter tube causing the pellet to shatter. The process of pellet fragmentation is a chaotic process that can be described in terms of fragment size distribution through a statistical model that incorporates the effects of the pellet material and impact characteristics. In addition to the fragment size distribution, the shatter plume has other characteristics of interest, such as a fragment velocity distribution and temporal mass evolution. The fragment velocity distribution is important because it is needed to accurately model the spread and location of the ablation and the deposition of impurities in the plasma over time. The temporal mass evolution is necessary to determine the time-resolved delivery of mass to the plasma.Due to installation constraints, the shatter tube currently installed on JET has a unique geometry with a modest S-bend followed by a 20-deg bend at the end of the tube. The DIII-D and KSTAR shatter tube design is a simple tube bent through an angle of 20 deg followed by a straight section. The resulting shatter sprays from the JET shatter tube and a 20-deg miter bend shatter tube were experimentally characterized for various pellet materials and speeds. Laboratory testing of these shatter tubes allows for the use of fast cameras to capture the fragment spray traveling through a large vacuum chamber. These high-speed videos of the shatter plumes allow the fragment size distribution, temporal mass evolution, and velocity distribution of the fragments within the plume to be determined. Overall, this paper presents a comparison of the unique geometry of the JET shatter tube to the miter bend geometries used for shattering and some insight into the variables that may be adjusted to produce the optimal shatter spray. The impact of entrained propellant gas on the resulting shatter spray was examined during testing.
Rare earth elements (REEs), as critical minerals, have important uses in modern energy and technologies, yet are vulnerable to potential supply chain disruptions. To establish domestic REE supply chain, efficient REE detection methods for resource characterization and mineral processing will be needed to accelerate innovations for domestic REE recovery. This study developed a rapid, novel, and cost-effective for REE detection method using ion chromatography (IC) for aqueous samples. Various REE-targeted eluent gradients and post-column agent compositions were tested on the chelation ion chromatography (CIC) with UV-vis detector for optimal separation and quantification of REEs within approximately 20 min. The single-channel pump to deliver the post-column solution to UV-vis detector was replaced with a 4-channel gradient pump, to increase operation and maintenance efficiencies. After method optimization, resulting calibration curves for more than ten REEs achieved high coefficients of determination (R2>0.999) and low relatively standard deviations (below 3.24%), demonstrating sub-ppm level detection limits (0.0897 to 0.1149 mg/L). The reliability of the CIC method was validated through comparison with inductively coupled plasma mass spectrometry (ICP-MS), showing strong agreement in REE recovery from certified standards. The impact of metal ions and salts on REE recovery using CIC was also systematically investigated. CIC consistently exhibited reliable performance in the presence of salt solutions such as NaCl and Na₂SO₄ (up to 10,000 mg/L). Our study also found the presence of high concentrations of Al ions (at 10,000 mg/L) significantly influenced REE determination, and elevated concentrations of Ca ions affected the recovery of specific REEs, including La, Ce, and Pr. The CIC method was further tested on REE-containing eluents from solvent extraction tests out of fly ash leachates. REE detection from these real processing fluids were reported to achieve 90% to 100% recovery rate from our IC method, compared to ICP-MS results. This study underscores the potential of CIC as a reliable and efficient alternative for REE determination in complex matrices. It also highlights the importance of minimizing select interfering metal ions in solutions to ensure accurate results. The REE CIC method presents a promising, low-maintenance, salt-tolerant, and cost-effective alternative to traditional analytical methods for REE analysis.
A numerical study of magnetohydrodynamics (MHD) and tracer-particle evolution investigates the effects of resonant magnetic perturbations (RMPs) on the confinement of runaway electrons (REs) in tokamak discharges conducted in the Madison Symmetric Torus. In computational results of applying RMPs having a broad toroidal spectrum but a single poloidal harmonic, m = 1 RMP does not suppress REs, whereas m = 3 RMP achieves significant deconfinement but not the complete suppression obtained in the experiment. MHD simulations with the NIMROD code produce sawtooth oscillations, and the associated magnetic reconnection can affect the trajectory of REs starting in the core region. Simulations with m = 3 RMP produce chaotic magnetic topology over the outer region, but the m = 1 RMP produces negligible changes in field topology, relative to applying no RMP. Using snapshots of the MHD simulation fields, full-orbit relativistic electron test particle computations with KORC show ≈50% loss from the m = 3 RMP compared to the 10%–15% loss from the m = 1 RMP. Here, test particle computations of the m = 3 RMP in the time-evolving MHD simulation fields show correlation between MHD activity and late-time particle losses, but total electron confinement is similar to computations using magnetic-field snapshots.