Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ELMs”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Simultaneous ELM suppression and divertor detachment via synergistic boron powder and neon injection in EAST

A novel approach for simultaneous power exhaust and edge-localized mode (ELM) control is presented in the Experimental Advanced Superconducting Tokamak discharges, which utilize an ITER-like tungsten divertor. Real-time injection of boron (B) powder and neon (Ne) gas overcomes their limitations encountered when used separately. Pure Ne seeding leads to a narrow operational window constrained by core impurity accumulation and H-mode to L-mode back transitions, while pure solid B injection (SBI) is insufficient for effective divertor cooling. In comparison, their combined use achieves a stable, stationary, ELM-suppressed H-mode with adequate power exhaust. This synergistic scenario features partial energy detachment at the outer divertor while maintaining good plasma confinement (H 98 ∼ 1) with minimal degradation. Two key features of this scenario are: (1) the SBI triggers a persistent Edge Harmonic Mode (EHM), which provides a crucial continuous particle transport channel, preventing Ne and tungsten/molybdenum accumulation without flushing out by ELM, and (2) the B + Ne mixture allows for active optimization of the radiated power profile. Core radiation can be reduced by substituting a portion of the Ne with B, leveraging their complementary non-coronal equilibrium radiation efficiencies. This combined B + Ne injection scheme presents a promising pathway toward integrated core-edge scenarios, offering the potential to minimize total impurity throughput while leveraging an actuator (powder injection) already being considered for ITER.

Boron↗
DOE OSTI2025

Assessing Edge Localized Mode (ELM) Mitigation by Pellet Triggering in Smallest-Possible Affordable Robust Compact (SPARC) Low-Collisionality Discharges

Simulations of pellet edge-localized mode (ELM) triggering with the M3D-C1 code confirmed experimentally observed pellet mass thresholds in low-collisionality DIII-D discharges and predicted extremely large thresholds in ITER [A. Wingen et al, Nucl. Fusion 64, 066015(2024)]. The linear (static) and nonlinear (dynamic) workflows developed in this previous work are now applied to the SPARC tokamak. Multiple synthetic equilibria for the high-performance H-mode scenario are tested for MHD stability without pellets using the linear workflow and ELITE simulations; a stable one is identified. Pellets of various mass are injected along an inboard, near-midplane trajectory in the M3D-C1 simulation. No pellet mass threshold for ELM triggering can be found for realistic pellet masses, concluding that fueling from this injection location would not trigger ELMs. A poloidal scan of possible injection locations shows very high thresholds on the high-field side and reasonable thresholds on the low-field-side (LFS). Based on results of this work the location for a LFS injection was finalized in the SPARC design. For this LFS injection the pellet mass threshold is determined using the linear workflow, aligning well with the poloidal scan.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Developing an ELM Ecosystem Dynamics Model on GPU with OpenACC

Porting a complex scientific code, such as the E3SM land model (ELM), onto a new computing architecture is challenging. The paper presents design strategies and technical approaches to develop an ELM ecosystem dynamics model with compiler directives (OpenACC) on NVIDIA GPUs. The code has been refactored with advanced OpenACC features (such as deepcopy and routine directives) to reduce memory consumption and to increase the levels of parallelism through parallel loop reconstruction and new data structures. As a result, the optimized parallel implementation achieved more than a 140-time speedup (50 ms vs 7600 ms), compared to a naive implementation that uses OpenACC routine directive and parallelizes the code across existing loops on a single NVIDIA V100. On a fully loaded computing node with 44 CPUs and 6 GPUs, the code achieved over a 3.0-times speedup, compared to the original code on the CPU. Furthermore, the memory footprint of the optimized parallel implementation is 300 MB, which is around 15% of the 2.15 GB of memory consumed by a naive implementation. This study is the first effort to develop the ELM component on GPUs efficiently to support ultra-high-resolution land simulations at continental scales.

Schwartz, Peter↗
DOE OSTI2026

Machine Learning Calibration of Groundwater Table Depth in ELM: Impact on Land Surface Hydrology and Land‐Atmosphere Fluxes

Accurate representation of groundwater table depth (GWTD) is crucial for simulating hydrological cycling in Earth system models (ESM). Nevertheless, there is a notable gap in the literature regarding the validation of GWTD simulations in ESMs and their subsequent impact on downstream hydrological components. This study explores the calibration of parameterization of global GWTD using machine learning within the Energy Exascale Earth System Model (E3SM) Land Model (ELM). Despite achieving significant gains in simulating GWTD through calibration, offline ELM simulations unexpectedly show that these improvements do not translate to substantial enhancements in model performance for other key hydrological variables, including soil moisture (SM), runoff, groundwater contribution to runoff or base flow index (BFI), and evapotranspiration and its partitioning. The performance in SM and runoff was even degraded in some regions, while BFI was mostly overestimated. Although there is significant improvement in GWTD within the critical range of 1–5 m, where groundwater traditionally influences land surface energy fluxes, these improvements occurred mostly in humid areas where the impact of GWTD on surface processes is minimal. Although the impacts of model calibration are generally small in offline ELM simulations, coupled land-atmosphere simulations exhibit much stronger responses to GWTD calibration, highlighting the role of land-atmosphere feedbacks in Earth system modeling. These findings underscore the need for integrated calibration strategies that simultaneously optimize multiple hydrological variables. However, if a single-variable approach is necessary, it is crucial to establish clear priorities for calibration, identifying the most critical variables that have the greatest impact on overall model performance.

Fang, Yilin [Pacific Northwest National Laboratory↗
DOE OSTI2024

Two-stage crash process in resistive drift ballooning mode driven ELM crash

We report a two-stage crash process in edge localized mode (ELM) driven by resistive drift-ballooning modes (RDBMs) numerically simulated in a full annular torus domain with a scale-separated four-field reduced MHD (RMHD) model using the BOUT++ framework. In the early nonlinear phase, the small first crash is triggered by linearly unstable RDBMs, and m/n=2/1 magnetic islands are nonlinearly excited by nonlinear coupling of RDBMs as well as their higher harmonics. Here, m is the poloidal mode number, n is the toroidal mode number, the q = 2 rational surface exists near the pressure gradient peak, and q is the safety factor. Simultaneously, middle-n RDBM turbulence develops but is poloidally localized around X-points of the magnetic islands, leading to the small energy loss. The second large crash occurs in the late nonlinear phase. Higher harmonics of m/n=2/1 magnetic islands well develop around the q = 2 surface via nonlinear coupling and make the magnetic field stochastic by magnetic island overlapping. Turbulence heat transport develops at X-points of higher harmonics of m/n=2/1 magnetic islands, resulting in the turbulence spreading in the poloidal direction. The large second crash is triggered when the turbulence covers the whole poloidal region so that the magnetic island generation and magnetic field stochastization before the large crash can be interpreted as ELM precursors. It is concluded that the ELM trigger is attributed to the turbulent spreading in the poloidal direction in synchronization with the magnetic field stochastization and the crash is driven by E × B convection rather than the conventional Rechester–Rosenbluth anomalous electron heat transport.

Seto, H. (ORCID:0000000259278643)↗
DOE OSTI2024

Synergistic effects of W–Ta alloys exposed to various ELM-like fusion events

Abstract Recently tungsten alloys have shown capabilities in suppressing morphological changes when exposed to harsh fusion conditions. Morphological changes due to He + and D + ions at high temperatures (1000–2000 K) can lead to easier erosion and sputtering of tungten (W) particles towards quenching the fusion plasma. This study investigated W–3Ta alloy as potential plasma facing materials (PFMs) exposed to different edge-localized mode (ELM)-like energies and ion loadings. Performance of W–3Ta has been tested as a function of ELM-like energies (0.6–1.5 GW m −2 ) and varying He + and D + ion concentrations (100% He + , 50% He + :50% D + ,10% He + :90% D + (similar ELM-like reactor condition), and 100% D + ions). Early-stage fuzz formation and surface pores were observed during 100% He + ion irradiation and mixed ion loadings but not observed during 100% D + ion loadings. Pores decreased in quantity and increased in size with increasing transient heat loading. Grain boundaries acted as trapping sites of helium and deuterium atoms during dual beam ion irradiations, as observed through increased pore size along the grain boundary. Qualitative erosion estimates were measured using witness plates. Increased erosion was observed with increasing D + ion concentration, which is attributed to D supersaturation in the near surface region affecting thermal and mechanical properties. More erosion was observed for W–3Ta samples compared to pure W, evident through large (>10 µ m) cracks on the surface, and an increase in W metal x-ray photoelectron spectroscopy peaks under similar loading conditions. The slight enhancement in erosion is attributed to the slightly decreased thermal properties, ∼10% thermal diffusivity drop, of W–3Ta at high temperatures. The damage threshold was impacted by ion loading, where distinct cracking was observed for 100% He + at a lower heat loading. This matrix study using reactor-like loading supports the need for extensive studies on future W alloy PFMs before selecting one as an alternative.

Physics↗
DOE OSTI2024

Exploring the effect of ELM and code-coupling frequencies on plasma and material modeling of dynamic recycling in divertors

Abstract Integrated modeling of plasma-surface interactions provides a comprehensive and self-consistent description of the system, moving the field closer to developing predictive and design capabilities for plasma facing components. One such workflow, including descriptions for the scrape-off-layer plasma, ion-surface interactions and the sub-surface evolution, was previously used to address steady-state scenarios and has recently been extended to incorporate time-dependence and two-way information flow. The new model can address dynamic recycling in transient scenarios, such as the application presented in this paper: the evolution of W samples pre-damaged by helium and exposed to ELMy H-mode plasmas in the DIII-D DiMES. A first set of simulations explored the effect of ELM frequency. This study was discussed in detail in this conference’s proceedings and is summarized here. The 2nd set of simulations, which is the focus of this paper, explores the effect of code-coupling frequency. These simulations include initial SOLPS solutions converged to the inter-ELM state, ion impact energy ( E in ) and angles ( A in ) calculated by hPIC2, and an improved heat transfer description in Xolotl. The model predicts increases in particle fluxes and decreases in heat fluxes by 10%–20% with the coupling time-step. Compared with the first set of simulations, the less shallow impact angle leads to smaller reflection rates and significant D implantation. The higher fraction of implanted flux (and deeper), in particular during ELMs, increases the accumulated D content in the W near-surface region. Future expansion of the workflow includes coupling to hPIC2 and GITR to ensure accurate descriptions of E in and A in , and W impurity transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Wall conditioning and ELM mitigation with boron nitride powder injection in KSTAR

Results from KSTAR powder injection experiments, in which tens of milligrams of boron nitride (BN) were dropped into low-power H-mode plasmas, show an improvement in wall conditions in subsequent discharges and, in some cases, a reduction or elimination of edge-localized modes (ELMs). Injected powder is distributed by the plasma flow and is deposited on the wall and, over the course of several discharges, was observed to gradually reduce recycling by 33%, and decrease both the ELM amplitude and frequency. This is the first demonstration of the use of BN for ELM mitigation. In all of these experiments, an Impurity Powder Dropper (IPD) was used to introduce precise, controllable amounts of the materials into ELMy H-mode KSTAR discharges. The plasma duration was between 10 s and 15 s, 𝐼𝑝 = 500 kA, 𝐵𝑇 = 1.8 T, 𝑃NBI = 1.6 MW, and 𝑃ECH = 0.6 MW. Plasma densities were between 2 and 3 × 1019 m−3. In all cases, the pre-fill and startup gas-fueling was kept constant, suggesting that the decrease in baseline D𝛼 emission is in fact due to a reduction in recycling. The results presented herein highlight the viability of powder injection for intra-shot and between-shot wall conditioning.

Gilson, Erik↗
DOE OSTI2026

elm-diagnostics

elm-diagnostics is a Python package for computing diagnostic analyses and visualizations for the E3SM Land Model (ELM) component and is meant to support new feature development in ELM. The tool reads model history files and performs quantitative analyses including budget-closure checking, variable transformations, temporal aggregations, and statistical summaries to support model evaluation, validation, and scientific interpretation. The framework is designed for extensibility, with modular architecture enabling straightforward addition of new diagnostic methods, derived variables, analysis types, visualization approaches, and model-specific adaptations

Hoffman, Matt [Los Alamos National Laboratory]↗
DOE OSTI2025

Assessing Edge Localized Mode (ELM) Mitigation by Pellet Triggering in Smallest-Possible Affordable Robust Compact (SPARC) Low-Collisionality Discharges

Simulations of pellet edge-localized mode (ELM) triggering with the M3D-C1 code confirmed experimentally observed pellet mass thresholds in low-collisionality DIII-D discharges and predicted extremely large thresholds in ITER. The linear (static) and nonlinear (dynamic) workflows developed in this previous work are now applied to the SPARC tokamak. Multiple synthetic equilibria for the high-performance H-mode scenario are tested for MHD stability without pellets using the linear workflow and ELITE simulations; a stable one is identified. Pellets of various mass are injected along an inboard, near-midplane trajectory in the M3D-C1 simulation. No pellet mass threshold for ELM triggering can be found for realistic pellet masses, concluding that fueling from this injection location would not trigger ELMs. A poloidal scan of possible injection locations shows very high thresholds on the high-field side and reasonable thresholds on the low-field-side (LFS). Based on results of this work the location for a LFS injection was finalized in the SPARC design. For this LFS injection the pellet mass threshold is determined using the linear workflow, aligning well with the poloidal scan.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Active wall conditioning through boron powder injection compatible $\mathrm{ELM}$ control in $\mathrm{EAST}$

Real-time wall conditioning and edge-localized mode (ELM) suppression through real-time boron powder injection into high-confinement mode discharges of the Experimental Advanced Superconducting Tokamak (EAST) have been investigated. The boron powder injection effectively conditioned the plasma-facing components, as shown by the reductions in low-Z and high-Z impurities and the global recycling level in the whole inner vessel, and suppressed the strong magnetohydrodynamic (MHD) activity in the plasma core. The improved wall conditions also resulted in an earlier transition from L- to H-mode. Boron powder injection also effectively suppressed ELM in the lower single null configuration with BT in the favorable direction while leaving the plasma density and stored energy largely unaffected.

36 MATERIALS SCIENCE↗
DOE OSTI2024

ELM-free H-mode phase and decoupling of peeling–ballooning stability boundary in the MAST Upgrade tokamak

A linear magnetohydrodynamic (MHD) peeling–ballooning stability analysis of the edge-localized mode (ELM)-free phase of a MAST Upgrade (MAST-U) H-mode plasma is presented. In contrast to other similar discharges, #47018 is found to have a significantly higher and wider pedestal during its ELM-free H-mode phase that lasts for approximately 80 ms; this is made possible by the reduced core MHD mode activity on the q = 2 surface. During this period, there is sustained decoupling of peeling and ballooning branches of the stability boundary on J–α space, opening an access channel to the second stability regime with higher peaks in pedestal current density $J_{\mathrm{N,ped}}$ and pressure gradient (α). Such decoupling of the stability boundary has not previously been observed in MAST-U H-modes, and if such a condition can be readily reproduced, it opens a wide range of opportunities for MAST-U to explore low-collisionality peeling-limited pedestal regimes as well as advanced scenarios such as quiescent H-modes that are relevant to future reactors such as STEP and ITER.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Time-dependent experimental identification of inter-ELM microtearing modes in the tokamak edge on DIII-D

In a series of discharges on the DIII-D tokamak, fast vertical plasma jogs are used to induce current perturbations in the steep gradient region of the H-mode edge. These current perturbations directly influence the edge q profile, decoupling the resonant location and instability drive of pedestal-localized microtearing modes (MTMs.) By exploiting this effect, we develop and apply a new experimental technique to track the dynamical frequency evolution of MTMs in the pedestal region, providing a compelling validation of the MTM model. The frequency of potential MTMs is calculated as the Doppler-shifted electron diamagnetic frequency at rational q=m/n surfaces, showing remarkable agreement with chirped frequency behavior of n=3, 4 and 5 modes detected with fast magnetics. Data is collected throughout multiple ELM cycles in order to build robust statistics describing the time-dependent frequency evolution of MTMs, which can be explained by examining the recovery of pedestal gradients after an ELM event. MTMs have a dominant transport contribution in the electron thermal channel, so the presented results indicate that reduced models of pedestal transport must be electromagnetic in nature and constructed with accurate calculations of MTM stability; inclusion of this physics is essential for accurate predictions of the electron temperature pedestal profile. Supporting measurements of mode saturation, propagation direction and transport fingerprints are made to support the dynamic frequency determination, unambiguously and experimentally identifying MTMs in the pedestal region of DIII-D.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Permanent magnets for ELM suppression in tokamaks: feasibility and operational compatibility

Permanent magnets can provide static magnetic fields without power supplies or feed lines, offering a simpler alternative to conventional electromagnets in tokamaks. This work examines permanent magnet arrays (PMAs) for edge localized mode (ELM) control in double-null (DN) configurations, where traditional RMP coils exhibit limited effectiveness. Linear plasma response calculations using IPEC assess high-field-side (HFS) permanent magnet placement, benchmarking performance against DIII-D low-field-side (LFS) internal coils (I-coils) while evaluating engineering constraints. Modeling results demonstrate that HFS permanent magnet configurations generate HFS plasma response amplitudes more than 5 times larger than conventional I-coil systems. While other resonant response metrics show comparable or reduced response relative to I-coils, the combination of conventional RMP systems and permanent magnets is expected to provide enhanced performance. Operational impact assessments of persistent magnetic fields, including sideband field effects and startup/ramp-up compatibility, reveal acceptable performance within established operational boundaries. This analysis establishes PMAs as a technically viable approach for DN ELM control with reductions in system complexity, motivating further experimental validation on existing tokamak facilities.

3D magnetic perturbation↗
DOE OSTI2026

Spatiotemporal structure of edge harmonic oscillation and its role in ELM-free QH-mode at KSTAR

In Quiescent H-mode (QH-mode), edge-localized modes (ELMs) are naturally replaced by a low-n edge harmonic oscillation (EHO), yet the self-regulating transport mechanism driven by the EHO remains insufficiently understood. Using high-spatiotemporal-resolution imaging diagnostics on KSTAR—electron cyclotron emission imaging and broadband electron cyclotron emission—we resolve the eigenmode structure of the EHO and elucidate its regulatory role in edge transport. The EHO is localized within the pedestal near the maximum pressure gradient, propagates in the ion-diamagnetic direction, and its radial envelope expands with increasing shear, suggesting that rotational shear is closely associated with the structural evolution of the EHO. Information-theoretic Transfer Entropy analysis identifies a distinct ‘dual-stabilization’ regulatory interaction pattern: the EHO is associated with enhanced outward energy transport to limit the pressure gradient, while the background shear flow is associated with preferential regulatory influence on the EHO saturation amplitude and energy flux. These results demonstrate that the QH-mode pedestal is sustained by a shear-associated regulation mechanism, where the coupling between rotational shear and mode structure appears closely linked to the sustained ELM-free state.

ECEI↗
DOE OSTI2022

Integration of high confinement, high poloidal beta plasma with dual radiated power and detachment controls for divertor protection and ELM suppression

Divertor detachment without serious core confinement quality loss in DIII-D’s high poloidal β scenario has been combined with impurity-induced ELM mitigation without disruption. Use of Ne previously granted access to a detached, non-ELMing regime that retained high confinement quality due to stimulation of Internal Transport Barrier (ITB) growth, but suffered from on-going core fuel dilution and high disruptivity. Excess Ne accumulation in the core plasma has now been avoided by feeding back core radiated power ($P_{rad}$) measurements to control Ne seeding, rather than using attachment fraction ($A_{frac}$) control with Ne; this also reduces disruptivity. At the same time, N 2 seeding is used in a feedback loop with $A_{frac}$ measurements, which previously posed low disruption risk. In this way, the effect of Ne in the core is managed while avoiding excess seeding, and N 2 acts to correct for any excess heat exhaust that might interfere with detachment. The average Ne flow rate was 38% of what was used in pure Ne $A_{frac}$ control, plus average N 2 flow that was 43% of pure N 2 $A_{frac}$ control, even while meeting an even deeper detachment target. Meanwhile, a steep pressure gradient in the core plasma reduces sensitivity to impurity-induced degradation of the pedestal and reasonable confinement quality was maintained despite operational challenges that blocked formation of an ITB in these experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

New proposed ITER divertor design using carbon insert on tungsten to mitigate ELMs and secondary radiation effects on nearby components

Abstract Building a successful device for the magnetic fusion energy production is a great challenge. ITER is an international project of the tokamak based magnetic fusion design being developed for the demonstration of the feasibility of thermonuclear technologies for future realization of successful commercial fusion energy. A key obstacle to a successful magnetic fusion energy production is however, the performance during abnormal events including plasma disruptions and edge-localized modes (ELMs). A credible reactor design must tolerate at least a few of these transient events without serious consequences such as melting of the structure. This paper investigates and compares the performance of the current ITER tokamak design during two types of transient events, i.e., ELMs occurring at normal operation and disruptions during abnormal operation. We simulated the divertor components response using our integrated 3D HEIGHTS package. The simulations include self-consistent modeling of the interaction of the released core plasma particles with the initial solid divertor material, energy deposition processes, vaporization of divertor material, secondary plasma formation and MHD evolution, incident core particles collisions and scattering from this dense secondary plasma, photon radiation of secondary plasma, and the resulting heat loads on nearby components. Our simulations showed that using a small carbon insert around the strike point can significantly reduce the overall expected damage on the tungsten dome structure, reflector plates, and prevent tungsten vaporization and its potential core plasma contamination.

47 OTHER INSTRUMENTATION↗
DOE OSTI2024

ELM-free enhanced D α H-mode with near zero NBI torque injection in DIII-D tokamak

Enhanced D α H-mode (EDA H-mode), an ELM-free H-mode regime, and the concomitant quasi-coherent mode (QCM) are explored in neutral beam heated, lower single null plasmas with near zero external torque injection. This regime exhibits good energy confinement (H 98y2 ~ 1) with β N ~ 2, high density, regime access at low input power, and no ELMs. New time–resolved measurements of electron and ion density, temperature, plasma rotation, and radial electric field as well as turbulence, flows, thermal and particle transport, and linear instability calculations examining the dynamics of the QCM are presented. Measurements using multiple fluctuation diagnostics reveal the QCM to be a separatrix spanning mode, peaking just inside the separatrix, existing in a wide range of k ⊥ ρ s ~ 0.1–1.2 with multiple harmonics, and propagating with a very small phase velocity in the plasma frame, where k ⊥ is the binormal wavenumber and ρ s is the ion sound radius. Linear gyrokinetic simulations of this EDA H-mode discharge with CGYRO indicate that the trapped electron mode (TEM) and electron temperature gradient mode are the dominant instabilities in the radial region where QCM is also unstable. Qualitative analysis indicates that the properties of TEM are consistent with the experimentally observed characteristics of the QCM. These similarities suggest that the QCM is consistent with a TEM instability existing in the edge region of the EDA H-mode plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗