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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 73 records · Page 4

All sources NASA NTRS DOE OSTI

DOE OSTI2026

Impact of toroidal magnetic field direction on integrated ELM-stable operation and divertor power exhaust via boron powder injection in EAST

We report the first in-depth comparison of the impact of toroidal magnetic field direction on solid boron injection used for Edge-Localized Mode (ELM) control, power exhaust, and core high-Z impurity control in the Experimental Advanced Superconducting Tokamak. With favorable ion ∇B drift towards the upper X-point in an upper-single-null configuration, boron injection effectively suppresses ELMs, produces a detachment of the inner divertor target, and leads to improved energy confinement. ELM suppression in this configuration is accompanied by the excitation of an Edge Harmonic Mode. In contrast, with unfavorable ion ∇B drift away from the upper X-point, boron injection also suppresses ELMs but leads to a more symmetric detachment state of both the inner and outer divertor targets, while plasma energy confinement is slightly degraded despite similar boron injection levels; a different low-frequency coherent mode without multiple harmonics is observed. Measurements from toroidally separated views show that the divertor response to boron injection is essentially toroidally symmetric, supporting the use of two-dimensional SOLPS-ITER modeling with a toroidally uniform impurity source. These experimental observations are qualitatively consistent with SOLPS-ITER simulations, which highlight the critical role of E × B drift effects in setting the Bt-dependent in–out asymmetry of detachment and in asymmetrically transporting particles and injected impurities within the scrape-off layer and private-flux region. These findings underscore the importance of drift physics and real-time wall conditioning in controlling low-Z impurity transport and optimizing edge solutions for integrated, ELM-stable, high-performance tokamak operation.

E × B drift↗

DOE OSTI2025

Representing Soil Microbial Dynamics and Organo‐Mineral Interactions in the E3SM Land Model (ELM‐ReSOM)

Explicit representation of soil microbial processes and interactions with biotic and abiotic processes in Earth System Models (ESMs) remains limited, despite their importance in biogeochemical cycles. To address this gap, which hinders prediction of global biogeochemial cycling and responses to atmospheric conditions, we integrated a microbe- and mineral-surface-explicit model, the Reaction-network-based model of soil organic matter and Microbes (ReSOM), into the Energy Exascale ESM (E3SM) land model (ELM). Here, we describe ELM-ReSOM and show a case study at a conifer forest in California. ELM-ReSOM accurately simulated surface CO 2 fluxes and SOM stocks, demonstrating improved representations of microbial and mineral interactions compared to the default ELM. We examined ELM-ReSOM sensitivity to microbial traits, enzyme properties, and organo-mineral interactions. Microbial traits such as the maximum mortality rate, transporter-density scaling factor, and maximum monomer assimilation rate were strong controllers of heterotrophic respiration, while these microbial traits and enzyme-related properties collectively influenced SOM stocks. Mineral surfaces primarily affected SOM stocks by adsorbing enzymes, thereby limiting depolymerization. Synergies among processes led to stronger impacts of parameters when evaluated together versus separately (i.e., most parameters had greater indirect than direct effects). For example, due to interactions of microbial necromass with mineral surface adsorption, the indirect effect of the maximum microbial mortality rate was 33% larger than its direct effect on SOM stock. Thus, microbial and enzyme dynamics and their interactions with mineral surfaces play critical roles in SOM cycling. Tackling the challenges of microbe-explicit models will advance understanding and modeling of SOM dynamics.

Tao, Jing [Lawrence Berkeley National Laboratory (↗

DOE OSTI2022

Local measurements of the pedestal magnetic field profile throughout the ELM cycle on DIII-D

We report new high speed localized measurements of the pedestal magnetic field during the edge localized mode (ELM) cycle of a DIII-D High confinement mode (H-mode) discharge indicate a temporally and spatial complex redistribution of the edge current density profile, $j_{edge}$. The measurement technique extracts the magnetic field magnitude, $\textit{B}$, via the spectral separation of Stark-split neutral beam radiation in the pedestal. Single spatial channel measurements from a novel spatial heterodyne spectrometer are validated in discharges with core current profile changes. The technique measures Stark-splitting changes that imply $\textit{B}$ changes as small as 1 mT with high time resolution (50 μs). At normalized poloidal flux $ψ_n$ = 1.0, $\textit{B}$ appears saturated in the inter-ELM period and then rapidly decreases in <200 μs by ~1%, before edge recycling emission begins to increase. Radially inboard of $j_{edge}$, $\textit{B}$ increases at the ELM crash. The behavior is consistent with a rapid collapse of $j_{edge}$ at the ELM crash and subsequent pedestal recovery. In some discharges, at $ψ_n$ < 0.96, changes in $\textit{B}$ are observed throughout the ELM cycle. In others, $\textit{B}$ recovers and is relatively stable until a few ms leading up to the next crash. Measurements of $\textit{B}$ during the H-mode transition show a large increase at $ψ_n$ = 1 with little change at $ψ_n$ = 0.9, consistent with the formation of the edge bootstrap current density peak. The $ψ_n$ = 0.9 spectrum is complicated by predicted changes to the Stark component intensities with density at the L–H transition.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

New understanding of inter-ELM pedestal turbulence, transport, and gradient behavior in the DIII-D tokamak

New observations of pedestal localized turbulence in the inter-ELM period of H-mode plasmas in DIII-D show that ion temperature gradient mode scale (ITG-scale) density fluctuation (ñ) increases immediately after each ELM crash and is quickly suppressed during the increase in local E × B shear. This excitation and subsequent suppression of ITG-scale ñ can explain the previously reported anomalous ion heat flux, Qi during the ELM (Viezzer et al 2017 Nucl. Fusion 57 022020). Shorter wavelength trapped electron mode scale (TEM-scale) ñ starts to increase at a critical pedestal temperature gradient (∇Te,ped) and saturates as local E × B shear, ${T}_{\mathrm{i}}^{\mathrm{C}\mathrm{6}+}/{T}_{\mathrm{e}}$ ratio, and ∇Te,ped saturate. This TEM-scale ñ, which has the potential to cause electron (and also ion) heat transport, is consistent with driving an anomalous electron heat flux Qe, where Qe is estimated between ELMs using experimental profiles and power balance. Both ITG- and TEM-scale ñ amplitude variations with background Ti/Te and ∇ne,ped are found to be consistent with theoretical predictions of these measured density fluctuations being ITG and TEM instabilities respectively. Furthermore, these new and unique observations on the nature of turbulence and their potential contributions to electron and ion heat fluxes at different ELM periods (i.e. collapse, recovery, and saturation phases of pedestal gradients) can significantly test and improve our pedestal predictive capabilities for ITER and other future fusion devices.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Simulation of pellet ELM triggering in low-collisionality, ITER-like discharges

3D nonlinear, as well as 2D linear M3D-C1 simulations are used to model ELM triggering by small pellets in DIII-D discharges in the ITER relevant, peeling-limited pedestal stability regime. A critical pellet size threshold is found in both experiment and modeling depending on pedestal conditions, pellet velocity and injection direction. Using radial injection at the outboard midplane, the threshold is determined by M3D-C1 for multiple time slices of a DIII-D low-collisionality discharge that has pellet ELM triggering. Experimental observations show that a larger pellet size than the standard 1.3 mm diameter is necessary for ELM triggering; 1.8 mm pellets triggered several ELMs in cases where a smaller pellet failed. The M3D-C1 simulations are in good agreement with these observations. While the 2D linear simulations give insight into the change of growth rates for various toroidal modes with pellet size, the 3D nonlinear simulations apply a pellet ablation model that mimics the actual injection with good match to the experiment. The 3D nonlinear simulation confirms the pellet ELM triggering for a pellet size larger than the threshold found by the linear simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Nonlinear MHD modeling of n = 1 RMP -induced pedestal transport and mode coupling effects on ELM suppression in KSTAR

We report fully suppressing edge-localized modes (ELMs), e.g., with resonant magnetic perturbations (RMPs), is essential to reach and sustain high-performance steady-state H-mode plasmas because large ELMs can significantly reduce the lifetime of divertor components in future tokamak reactors. RMP-driven ELM suppression in KSTAR has been modeled by coupling the neoclassical transport code PENTRC to the nonlinear 3D MHD code JOREK. We have found that the radial transport from the combined effects of the kink-peeling, tearing response, and neoclassical toroidal viscosity can explain the pedestal degradation observed in experiments. In addition, it has been found that the RMP response can increase the inter-ELM heat flux on the lower outer divertor by redistributing the heat transport between the divertor plates. In addition to the degraded pedestal, ELM suppression is also attributable to the RMP-induced mode interactions. While the linear stability of peeling-ballooning mode (PBMs) improves owing to the degraded pedestal, the PBM and RMP interaction increases the spectral transfer between edge harmonics, preventing catastrophic growth and the crash of unstable modes. Here, it turns out that the magnetic islands near the pedestal top can play a vital role in mediating the mode interactions.

mode coupling↗

DOE OSTI2024

Nonlinear modeling of ELM mitigation with RMP on HL-2A

Abstract Nonlinear modeling of mitigation of the edge localized mode (ELM) with resonant magnetic perturbation (RMP) is performed for the HL-2A tokamak, utilizing the three-dimensional (3D) magnetohydrodynamic code JOREK. Based on the 3D equilibrium established after application of the n = 1 ( n is the toroidal mode number) RMP at 4.9 kAt coil current with odd parity, ELM mitigation is successfully simulated consistent with the experimental result. Nonlinear simulations show strong mode coupling among toroidal Fourier harmonics, allowing redistribution of the magnetic energy such that the most unstable toroidal mode saturates at a lower level. This magnetic energy cascade offers an explanation of the RMP-induced ELM mitigation achieved in HL-2A. Detailed examination of the simulation results shows persistent resonant field screening even during the ELM mitigation phase. Finite plasma resistivity however does enable partial penetration of the resonant field thus modifying the edge magnetic topology and characteristics of the edge transport. Plasma radial profiles undergo pronounced changes around the pedestal region, when the magnetic energy of the most unstable toroidal mode reaches the maximum value. Systematic scans of the applied RMP coil current with the JOREK simulations find a threshold value of around 4.5 kAt required for achieving the ELM mitigation on HL-2A.

Physics↗

DOE OSTI2021

Soil Water Potentials (1990-2018) from a calibrated ELM-FATES, and rooting depth analyses scripts, PA-BCI, Panama

This dataset has ELM-FATES simulation outputs and scripts that accompany and reproduce the manuscript, Chitra-Tarak et al. (2021) Hydraulically-vulnerable trees survive on deep-water access during droughts in a tropical forest. New Phytologist. https: //doi.org/10.1111/nph.17464”. This dataset is organized into four parts: (1) source code for the model, Energy Exascale Earth System Land Model-Functionally Assembled Advanced Terrestrial simulator (ELM-FATES) (2) R scripts associated with ELM-FATES simulations: to prepare files for parameter ensemble simulation runs, calibrate the model with observed hydrological fluxes observed at BCI, re-run the model with tuned parameters and obtain long-term soil matric potentials by depth (3) ELM-FATES hydrological simulation outputs at PA-BCI, Panama (1990-2018, daily), and (4) R scripts for inverse rooting depth modeling and post-processing: prepare datasets and inversely model tree rooting depths at BCI, test hypotheses about the relationship of rooting depths with plant hydraulic traits and drought-induced mortality, and reproduce the manuscript itself. n.b. 1. Data needed to run R scripts in (2) and (4) and thus to reproduce the manuscript are sourced from various public datasets. Refer to the data availability section of the accompanying manuscript. 2. This dataset version (prepared on May 12, 2021) matches the revised and published version of the accompanying manuscript (published May 16, 2021) replacing a previous version. In the publication, this dataset is cited as follows: Chitra-Tarak R, Xu C, Aguilar S, Anderson-Teixeira K, Chambers J, Detto M, Faybishenko B, Fisher R, Knox R, Koven C et al. 2020. Soil water potentials (1990–2018) from a calibrated ELM-FATES, and rooting depth analyses scripts, PA-BCI, Panama. 2.0. NGEE Tropics Data Collection. doi: 10.15486/ngt/1696806

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

Landsat near-infrared (NIR) band and ELM-FATES sensitivity to forest disturbances and regrowth in the Central Amazon

Forest disturbance and regrowth are key processes in forest dynamics, but detailed information on these processes is difficult to obtain in remote forests such as the Amazon. We used chronosequences of Landsat satellite imagery (Landsat 5 Thematic Mapper and Landsat 7 Enhanced Thematic Mapper Plus) to determine the sensitivity of surface reflectance from all spectral bands to windthrow, clear-cut, and clear-cut and burned (cut + burn) and their successional pathways of forest regrowth in the Central Amazon. We also assessed whether the forest demography model Functionally Assembled Terrestrial Ecosystem Simulator (FATES) implemented in the Energy Exascale Earth System Model (E3SM) Land Model (ELM), ELM-FATES, accurately represents the changes for windthrow and clear-cut. The results show that all spectral bands from the Landsat satellites were sensitive to the disturbances but after 3 to 6 years only the near-infrared (NIR) band had significant changes associated with the successional pathways of forest regrowth for all the disturbances considered. In general, the NIR values decreased immediately after disturbance, increased to maximum values with the establishment of pioneers and early successional tree species, and then decreased slowly and almost linearly to pre-disturbance conditions with the dynamics of forest succession. Statistical methods predict that NIR values will return to pre-disturbance values in about 39, 36, and 56 years for windthrow, clear-cut, and cut + burn disturbances, respectively. The NIR band captured the observed, and different, successional pathways of forest regrowth after windthrow, clear-cut, and cut + burn. Consistent with inferences from the NIR observations, ELM-FATES predicted higher peaks of biomass and stem density after clear-cuts than after windthrows. ELM-FATES also predicted recovery of forest structure and canopy coverage back to pre-disturbance conditions in 38 years after windthrows and 41 years after clear-cut. The similarity of ELM-FATES predictions of regrowth patterns after windthrow and clear-cut to those of the NIR results suggests the NIR band can be used to benchmark forest regrowth in ecosystem models. Our results show the potential of Landsat imagery data for mapping forest regrowth from different types of disturbances, benchmarking, and the improvement of forest regrowth models.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Optimization of 3D controlled ELM-free state with recovered global confinement for KSTAR with n = 1 resonant magnetic field perturbation

Mitigation of deleterious heat flux from edge-localized modes (ELMs) on fusion reactors is often attempted with 3D perturbations of the confining magnetic fields. However, the established technique of resonant magnetic perturbations (RMPs) also degrades plasma performance, complicating implementation on future fusion reactors. Here, we introduce an adaptive real-time control scheme on the KSTAR tokamak as a viable approach to achieve an ELM-free state and simultaneously recover high-confinement (βN ~ 1.91, β p ~ 1.53, and H 98 ~ 0.9), demonstrating successful handling of a volatile complex system through adaptive measures. We show that, by exploiting a salient hysteresis process to adaptively minimize the RMP strength, stable ELM suppression can be achieved while actively encouraging confinement recovery. This is made possible by a self-organized transport response in the plasma edge which reinforces the confinement improvement through a widening of the ion temperature pedestal and promotes control stability, in contrast to the deteriorating effect on performance observed in standard RMP experiments. These results establish the real-time approach as an up-and-coming solution toward an optimized ELM-free state, which is an important step for the operation of ITER and reactor-grade tokamak plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Critical role of current-driven instabilities for ELMs in NSTX

The impact of different extended-magnetohydrodynamic (MHD) contributions on the stability thresholds of peeling-ballooning modes in ELMing and ELM-free plasmas in the spherical tokamak NSTX is investigated with the initial value code M3D-C1. We show that ELMing discharges in NSTX are limited by resistive current-driven peeling modes, whereas non-ELMing wide-pedestal H-mode discharges are located near the ideal pressure-driven ballooning threshold. It is demonstrated that extended-MHD can lead to more reliable edge stability predictions than existing ideal-MHD models. Resistive peeling-ballooning modes are found to exist well before the ideal stability threshold is met, and kink-peeling modes exhibit considerable sensitivity to plasma resistivity. Other effects not considered in ideal-MHD models affect PB modes in NSTX in a weaker way. Gyroviscous stress appears stabilizing such that the stability boundary lies closer to the experimental point. Equilibrium rotation can suppress ideal core modes and thus isolate edge modes. These results are important for the development of a predictive pedestal model for low-aspect ratio tokamaks.

pedestal↗

DOE OSTI2022

Characteristics of grassy ELMs and their impact on the divertor heat flux width

We report BOUT++ turbulence simulations are conducted for a 60 s steady-state long pulse high β p EAST grassy ELM discharge. BOUT++ linear simulations show that the unstable mode spectrum covers a range of toroidal mode numbers from low-n (n = 10–15) peeling–ballooning modes (P–B) to high-n (n = 40–80) drift-Alfvén instabilities. Nonlinear simulations show that the ELM crash is triggered by low-n peeling modes and fluctuation is generated at the peak pressure gradient position and radially spread outward into the scrape-off-layer, even though the drift-Alfvén instabilities dominate the linear growth phase. However, drift-Alfvén turbulence delays the onset of the grassy ELM and enhances the energy loss with the fluctuation extending to pedestal top region. Simulations further show that if the peeling drive is removed, the fluctuation amplitude drops by an order of magnitude and the ELM crashes disappear. The divertor heat flux width is ~2 times larger than the estimates based on the HD model and the Eich's ITPA multi-tokamak scaling (or empirical Eich scaling) due to the strong radial turbulence transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Metrics and extrapolation of resonant magnetic perturbation thresholds for ELM suppression

This large database study of resonant magnetic perturbation (RMP) edge localized mode (ELM) suppression thresholds in the AUG, DIII-D, EAST, and KSTAR tokamaks details the key strengths and weaknesses of RMP metrics. The RMP ELM suppression database used for this work contains plasma information at the time of transition from ELMing to ELM suppressed states where a clear experimental threshold is identified. The experimental threshold distributions are compared for five metrics: (1) the island overlap width, (2) pedestal top Chirikov overlap, (3) peeling edge displacement, (4) pedestal top resonant drive, and (5) edge dominant mode overlap. The distributions, the regularity of the dependence on RMP coil currents, and the sensitivities of a given metric to equilibrium reconstruction details are compared. The overlap metric proves to be a good compromise between including the appropriate plasma response physics and maintaining a numerical robustness. This quantity does not exhibit clear power-law scalings for projection, but machine learning can assist in predicting thresholds within the existing parameter ranges and providing uncertainty quantification of those predictions. Two new first-principles models, one utilizing a threshold from the non-linear Modified Rutherford equation evaluated at the pedestal top and one utilizing the SLAYER code to calculate the linear tearing threshold from torque balance, offer possible paths to extrapolation beyond the existing database parameter space.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

HL-2A's ELM cycle simulations by integrating BOUT++'s drift MHD and transport code

A new integrating model has been developed to couple tokamak edge multiscale magnetohydrodynamic (MHD) events and transport simulations, such as edge-localized mode (ELM) cycles. As a proof of principle, we first start from a set of three-field two-fluid model equations, which includes the pressure, current, and vorticity. Here, the equations are separated into the slowly evolving part of the axisymmetric component by taking a time average of the axisymmetric component. The time-averaged fluxes, which are quadratic in fluctuating quantities, act as driven terms for the time-averaged axisymmetric quantities that determine the plasma transport, and therefore the large-scale evolution of the plasma profiles. Then the HL-2A's ELM cycles are simulated using the model. Good agreements of ELM size and pedestal recovery time have been achieved for the solutions obtained from the coupled simulation compared with experiment. For one ELM cycle simulation, the coupled code can achieve a speedup of a factor of up to 30 over standalone code.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2026

ELM‐MOSART‐DOC: A Large‐Scale Riverine Dissolved Organic Carbon Model and Its Application Over the United States

Riverine dissolved organic carbon (DOC), primarily sourced from soil organic carbon (SOC), plays a crucial role in regional and global carbon cycles. However, the complexities of the underlying mechanisms and limited observations present significant challenges for predictive understanding of DOC at regional or larger scales. Recently, we developed a machine learning‐based (ML) map of DOC transformation rates, bridging the gap between SOC and DOC leaching flux and simplifying terrestrial DOC representation. Building on this advancement, we introduce ELM‐MOSART‐DOC, a DOC module integrated into the riverine component of the Energy Exascale Earth System Model (E3SM)—the Model for Scale Adaptive River Transport (MOSART). ELM‐MOSART‐DOC simulates DOC transport and transformation across both headwater streams and river networks, including those managed. Model validation demonstrates the ability of ELM‐MOSART‐DOC to accurately capture long‐term average DOC concentrations, with Kling‐Gupta Efficiency (KGE) scores of 0.58 and 0.76 at large and local stations, respectively. We further assess the impact of reservoirs through different simulation schemes, revealing that reservoirs significantly alter DOC fluxes by regulating streamflow patterns and promoting DOC mineralization. Model simulations indicate that reservoirs reduce total DOC flux from the Mississippi River into the ocean by 7.5%, with the long‐term average annual export decreasing from 3.34 to 3.14 teragrams (Tg) per year. ELM‐MOSART‐DOC integrates process‐based modeling with ML parameterization to enhance the predictive understanding of riverine biogeochemical processes. This approach reduces uncertainties in modeling regional and global carbon cycle ESMs and provides new insights into carbon cycling and its implications for global environmental change.

Li, Lingbo [Univ. of Houston, TX (United States); ↗

DOE OSTI2021

Enhancement of edge turbulence concomitant with ELM suppression during boron powder injection in EAST

A reproducible, quasi-stationary edge localized mode (ELM)-suppressed scenario was obtained over a wide range of plasma parameters by continuous injection of boron (B) powder into an upper-single null discharge in the experimental advanced superconducting tokamak [Sun et al., Nucl. Fusion 61, 014002 (2021)]. This powder-induced ELM-absent regime is associated with an edge harmonic mode (EHM) that provides continuous particle exhaust to maintain constant density without confinement degradation and impurity accumulation, the latter of which is often observed in ELM-free regimes. A flow rate threshold of B powder injection, leading to a threshold intensity of the EHM, is necessary for full ELM suppression. The fundamental harmonic of the EHM exhibits a toroidal mode number n = 1. The mode is observable in the entire poloidal cross section with a peak near the upper X-point in an upper-single null configuration. Here, the EHM spans radially across the pedestal and scrape-off layer, peaking inside the separatrix. The EHM appears to be insensitive to q 95 , heating power, plasma toroidal rotation, and pedestal collisionality.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Plasma performance and operational space without ELMs in DIII-D

A database of DIII-D plasmas without ELMs compares the operating space and plasma performance of stationary no-ELM regimes: RMP-ELM suppression, QH-mode, I-mode, EDA H-mode, L-mode, and negative triangularity L-mode (Neg-D). Here, operational space is documented in terms of engineering and physics parameters, revealing divergent constraints. Some operational space discriminants (like pedestal collisionality) are well known, while others, such as low torque & safety factor, or high power & density, are less commonly emphasized. Normalized performance (H-factor and betan) also discriminate the no-ELM regimes and favor regimes tolerant to power in DIII-D: RMP, QH, and Neg-D.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

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, I p = 500 kA, B T = 1 . 8 T, P NBI = 1 . 6 MW, and P ECH = 0 . 6 MW. Plasma densities were between 2 and 3 × 1 0 19 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

36 MATERIALS SCIENCE↗