Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ELM”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

The effectiveness of D 2 pellet injection in reducing intra-ELM and inter-ELM tungsten divertor erosion rates in DIII-D during the Metal Rings Campaign

Abstract Edge localized modes (ELMs) in H-mode plasmas erode plasma-facing components (PFCs) and lead to impurities in the core, reducing confinement. This study analyzes D 2 pellet injection on the DIII-D fusion experiment used as an ELM mitigation technique applied during the 2016 tungsten Metal Rings Campaign to reduce W erosion during ELMs. The 400.9 nm photon wavelength line emission intensity of tungsten atoms (WI) filterscope channels and Langmuir probes were used to infer the gross erosion rate of tungsten-coated tiles installed in the divertor of DIII-D. D 2 mass injection rates ranging from 34 to 41 arbitrary units (A.U.) and no D 2 injection resulted in a similar total W erosion rate during ELMs (intra-ELM). On average, results show a 29% increase in the total gross W erosion rate with intermediate mass injection rates (∼13–23 A.U.) compared to the no pellets and the highest injection rate cases. On average, the fast D 2 mass injection rate cases had 15% less erosion in the inter-ELM phase than the case with no pellets. Generally, higher D 2 mass injection rates increased the ELM frequency, and the highest injection rates reduced the average erosion per ELM and fractional carbon impurities at the top of the pedestal by nearly 40% when compared to the no-pellet case. As expected, a higher D 2 pellet injection rate led to a higher plasma density and lower plasma temperature in the divertor. Additionally, an increasing divertor inter-ELM plasma electron density directly correlated to more frequent pellet injection and a decrease in both the average gross intra-ELM W erosion and the total gross intra-ELM W erosion rate. Simulations of intra-ELM erosion using the ‘free-streaming plus recycling model’ (FSRM) underestimate W erosion during pellet injection by about 30% on average. The discrepancies between the experimental measurements and the FSRM intra-ELM W erosion predictions are postulated to be due to C/W material mixing. A simple analytic mixed-material model is presented and results in better agreement with the experimental data. These results highlight the importance of incorporating the effects of a mixed-material layer in the analysis of PFC erosion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Design and experimental demonstration of feedback adaptive RMP ELM controller toward complete long pulse ELM suppression on KSTAR

Operation of a fusion power plant requires robust edge localized mode (ELM) suppression simultaneously with high plasma performance. In this study, we describe a novel feedback adaptive resonant magnetic perturbation (RMP) ELM controller designed to address this problem by achieving optimized ELM suppression through the advanced application of 3D RMPs. From real-time $D_α$ data, the controller is achievable robust ELM suppression while simultaneously minimizing the applied RMP in order to enhance plasma performance. In real time, the instantaneous ELM-frequency is analyzed with an adaptive feedback algorithm to determine amplitudes and phases of RMP coil currents that will maximize plasma performance while maintaining ELM suppression. When applied through the KSTAR plasma control system in several experiments using $textiti{n}$ = 1 RMPs, robust ELM suppression is achieved and sustained in feedback while reducing the RMP strength to ~65% of its initial value. Minimization of the RMP strength in this manner not only allows for operation of longer discharges due to a decrease in flux consumption but also allows for a strong recovery of up to ~60% of $β_N$ throughout the ELM-free period.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Evolution of ELMs, pedestal profiles and fluctuations in the inter-ELM period in NBI and ECH dominated discharges in DIII-D

In DIII-D, it has been observed that ELM frequency decreases by 40% and ELM spacing becomes more regular in time when heating is changed from pure neutral beam injection (NBI) to predominantly electron cyclotron heating (ECH) in ITER Similar Shape plasmas. In comparison with pure the NBI discharges, pedestal fluctuations in magnetics and density increase in the ECH dominated discharges. Recovery of the pedestal profiles like electron density (ne), temperature (Te) and pressure (pe) shows marked differences for these two heating schemes. Average profiles in the last 30% of the ELM cycle shows higher Te, lower ne, and similar pe at the pedestal top for the ECH discharge when compared to the NBI discharge. Here, the gradient of Te (∇Te) is also steeper at the pedestal in the ECH discharge. Magnetic fluctuations show 3 distinct modes in 13~116 kHz in the ECH discharges only. ne fluctuations show two modes evolving in the inter-ELM period of the ECH discharge, a low frequency (400 kHz) quasi-coherent mode (LFQC) and high frequency (~2 MHz) broadband (HFB) fluctuations. Evolution of these modes has marked correspondence with the inter-ELM ∇Te recovery. A sharp decrease in the Dα baseline is observed whenever the LFQC weakens and the HFB grows, prior to each large ELM. Transport coefficients obtained from TRANSP shows that MTM and/or TEM are plausible candidates for the observed fluctuations. Linear gyrofluid simulation (TGLF) corroborates this characterization. TGLF shows that the linear growth rate of the most dominant mode peaks at ion-scale (kθρs ~ 0.4) at the pedestal steep gradient and the frequency is in the electron diamagnetic direction. It is proposed that increased fluctuations in the ECH dominated case, due to increased ∇Te, caused an increase in fluctuation-driven transport in the pedestal and slowed the pedestal recovery between ELMs, leading to a reduction in the ELM frequency.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Active tungsten expulsion in ELM-absent H-mode plasmas via on-demand ELM triggering with lithium granule injection

Lithium granules gravitationally injected into the upper X-point region demonstrated on-demand edge-localized modes (ELM) triggering in otherwise ELM-suppressed H-mode plasmas on the Experimental Advanced Superconducting Tokamak. Sub-mm lithium granules dropped into enhanced D-alpha H-mode plasmas achieved a high triggering efficiency, while enabling ELM frequencies from several to hundreds of hertz. Core radiation from heavy impurities, dominated by W, was reduced by up to 60%, and the normalized energy confinement increased by up to 30%. At low injection frequencies, ELMs of substantially reduced size compared to spontaneous type-I ELMs were observed. At high injection frequencies, a transition to a mixed ELM phase occurred, characterized by intermittent larger ELMs and suppression of the quasi-coherent mode, achieving the most significant W reduction and energy confinement improvement. These results highlight a promising pathway for active W control via controlled, small ELMs in long-pulse, high-performance scenarios.

EDA H-mode↗
DOE OSTI2021

ELM and inter-ELM heat and particle flux to a secondary divertor in the DIII-D tokamak

DIII-D ELMing discharges with a secondary X-point (XPT) inside the vessel were used to evaluate the power distribution between the primary and secondary divertors. The magnetic balance, defined by the distance between the primary and secondary separatrices at the outer midplane (dRsep), was varied from -5 mm (lower-single-null, LSN) to +16 mm (upper-single-null, USN). In this work, we find that the secondary divertor receives up to 1/3 of the total heat flux (primary + secondary divertors), and it is dominated by the ELM-driven flux. Both the radially-integrated and the peak heat flux to the secondary divertor decay below ~50% of the maximum as dRsep is varied from -5 to +6mm, but the integrated heat flux decay flattens out and it would take dRsep above ~25mm to reduce the ELM heat flux to the secondary divertor below 10% of the one deposited to a well-defined SN. Both the secondary inner and outer strike points receive heat flux during ELMs. The peak heat flux is comparable in both strike points for dRsep < 10mm, but the inner strike point receives less than 15% of the total flux. Furthermore, this is experimental evidence that the secondary inner divertor region receives significant ELM flux which should be taken into account for the design of future generation tokamaks.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

ELM-Wet: Inclusion of a wetland landunit with sub-grid representation of eco-hydrological patches and hydrological forcing in E3SM Land Model (ELM)

Wetlands emit the most biogenic methane (CH4) and present the greatest uncertainty in the global CH4 budget. Modeling these emissions is challenging due to the temporal and spatial variability in wetland structure and CH4 flux rates, along with complex interactions among hydrological, ecological, meteorological, and microbial processes that govern CH4 dynamics. To address these issues, we aim to enhance the accuracy of wetland representation in the U.S. Department of Energy’s Exascale Earth System Model (E3SM) Land Model, ELM. This effort led to the development of ELM-Wet, which incorporates a dedicated wetland landunit with subgrid representation of eco-hydrological patch types. We implement wetland-specific hydrology by imposing site-specific constraints on surface water levels, thereby allowing different patches to sustain varying depths of inundation. Additionally, we refined the calculation of aerenchyma transport diffusivity based on observed conductance across different vegetation types. We validated these enhancements through site-specific simulations of a coastal freshwater wetland, the Salvador WMA Freshwater Marsh (Ameriflux, site ID US-LA2), located in the coast of Louisiana (29.85N,90.29W) at an elevation of 0 m. The site was simulated with ELM-Wet and the default version ELMv1. We use Bayesian Optimization to parameterize CO2 and CH4 fluxes. Eddy covariance observations of CO2 and CH4 fluxes from 2012-2013 were used to train the model and data from 2021 were used for validation. In this repository, we include all the output of all the simulations performed using different versions of ELMv1 and ELM-Wet, all input required to the model, and the Matlab scripts we used to processed the output data. ELM-Wet_flmd.csv includes a detailed description of the datasets files.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Integrated RMP-based ELM-crash-control process for plasma performance enhancement during ELM crash suppression in KSTAR

The integrated Resonant Magnetic Perturbation (RMP)-based Edge-Localized Mode (ELM)-crash-control process aims to enhance the plasma performance during the RMP-driven ELM crash suppression, where the RMP induces an unwanted confinement degradation. In this study, the normalized beta ($\beta_\textrm{N}$) is introduced as a metric for plasma performance. The integrated process incorporates the latest achievements in the RMP technique to enhance $\beta_\textrm{N}$ efficiently. The integrated process triggers the n = 1 Edge-localized RMP (ERMP) at the L–H transition timing using the real-time Machine Learning (ML) classifier. The pre-emptive RMP onset can reduce the required external heating power for achieving the same $\beta_\textrm{N}$ by over 10% compared to the conventional onset. During the RMP phase, the adaptive feedback RMP ELM controller, demonstrating its performance in previous experiments, plays a crucial role in maximizing $\beta_\textrm{N}$ during the suppression phase and sustaining the $\beta_\textrm{N}$-enhanced suppression state by optimizing the RMP strength. The integrated process achieves $\beta_\textrm{N}$ up to ~2.65 during the suppression phase, which is ~10% higher than the previous KSTAR record but ~6% lower than the target of the K-DEMO first phase ($\beta_\textrm{N}$ = 2.8), and maintains the suppression phase above the lower limit of target $\beta_\textrm{N}$ (= 2.4) for ~4 s (~60$\tau_\textrm{E}$). In addition to $\beta_\textrm{N}$ enhancement, the integrated process demonstrates quicker restoration of the suppression phase and recovery of $\beta_\textrm{N}$ compared to the adaptive control with the n = 1 Conventional RMP (CRMP). The post-analysis of the experiment shows the localized effect of the ERMP spectrum in radial and the close relationship between the evolution of $\beta_\textrm{N}$ and the electron temperature.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Modeling, analysis, and code/data validation of DIII-D tokamak divertor experiments on ELM and non-ELM plasma tungsten sputtering erosion

Abstract We analyzed recent DIII-D tokamak tungsten divertor probe experiments using advanced, coupled, sputter erosion/redeposition, plasma, and surface response code packages. Modeling is done for ELMing H-mode, and L-mode plasmas, impinging on various size tungsten deposits on Divertor Material Evaluation System (DiMES) carbon probes. The simulations compute 3D, full kinetic, sub-gyromotion, impurity sputtering and transport, including changes in tungsten surface composition and response due to mixed deuterium and carbon ions irradiation. Per our analysis, ELM (edge localized mode) plasma sputtering in DIII-D mostly involves free-streaming high energy (∼500–1000 eV) D + and C +6 ions, with high near-surface plasma density. L-Mode sputtering is due to impurity sputtering (C, W) only, with lower density. All cases show complete redeposition of tungsten on the divertor, with significant redeposition on the tungsten spots themselves, and low self-sputtering. Comparison of ELM plasma gross tungsten erosion simulation results with in-situ spectroscopic data is good, as are code/data comparisons of net erosion using post-exposure Rutherford backscattering (RBS) data for the L-mode probes. The analysis, extrapolated to a full tungsten divertor, implies low net erosion and negligible plasma contamination from sputtering. These results support the use of high-Z plasma facing surfaces in ITER and beyond.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

ELM–Wet: Inclusion of a Wet–Landunit With Sub–Grid Representation of Eco–Hydrological Patches and Hydrological Forcing Improves Methane Emission Estimations in the E3SM Land Model (ELM)

Wetlands are the largest emitters of biogenic methane (CH 4 ) and represent the highest source of uncertainty in global CH 4 budgets. Here, we aim to improve the realism of wetland representation in the U.S. Department of Energy's Exascale Earth System Model land surface model, ELM, thereby reducing uncertainty of CH 4 flux predictions. We develop an updated version, ELM-Wet, where we activate a separate landunit for wetlands that handles multiple wetland-specific eco-hydrological patch functional types. We introduce more realistic hydrological forcing through prescribing site-level constraints on surface water elevation, which allows resolving different sustained inundation depth for different patches, and if data exists, prescribing inundation depth. We modified the calculation of aerenchyma transport diffusivity based on observed conductance per leaf area for different vegetation types. We use Bayesian Optimization to parameterize CO 2 and CH 4 fluxes in the developed wet-landunit. Site-level simulations of a coastal non-tidal freshwater wetland in Louisiana were performed with the updated model. Eddy covariance observations of CO 2 and CH 4 fluxes from 2012 to 2013 were used to train the model and data from 2021 were used for validation. Patch-specific chamber flux observations and observations of CH 4 concentration profiles in the soil porewater from 2021 were used for evaluation of the model performance. Our results show that ELM-Wet reduces the model's CH 4 emission root mean squared error by up to 33% and is able to represent inter-daily CO 2 and CH 4 flux variability across the wetland's eco-hydrological patches, including during periods of extreme dry or wet conditions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Understanding the negative triangularity ELM trigger and ELM free state on DIII-D with ECE-imaging

Here, the Electron Cyclotron Emission Imaging (ECEI) diagnostic was used to observe a finite-n interchange mode structure in the edge of negative triangularity shaped plasmas on DIII-D. At a small negative triangularity (δ u = -0.2), the plasma is in the H-mode with ELMs that are triggered by a low-n interchange mode. At a larger negative triangularity (δ u = -0.4) and low NBI power (2 MW), a dithering oscillation is observed that is triggered by a low-n interchange mode, whereas at higher NBI power (>2 MW), the edge reverts to L-mode and the low-n interchange mode is present continuously. In all cases, the edge pressure gradient is clamped when the interchange mode is present. It is concluded that the low-n interchange mode prevents the plasma from transitioning to H-mode at a large negative triangularity. This agrees with linear BOUT++ simulations which suggest that the interchange-type MHD can be a resistive ballooning mode whereby resistivity can significantly increase the finite-n ballooning mode growth rate. The absence of H-mode at large negative triangularity can, thus, be explained by the excitation of low-n pressure driven resistive ballooning modes in the plasma edge.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Experimental evidence of enhanced radial transport in small ELM regimes at DIII-D

Small/type-II edge-localized-modes (ELMs), carrying 1% of the plasma stored energy, are found to deposit only 45 ± 5% of the ELM power near the strike point, and the remaining 55 ± 5% to the far scrape-off-layer (SOL). Small ELMs spread their power over a larger area compared to type-I ELMs, where such a ratio is about 60% and 40% for near- and far-SOL regions, respectively. The larger spread is reflected in the heat flux width (λq) in the SOL for the intra-small ELMs profile of 6.0 mm, almost a factor 2 larger than that of type-I ELMs of 3.15 mm, for similar plasma conditions and magnetic configuration. At the ELM peak, the small ELMs λq is found to be up to 4 times larger than for the type-I ELMs, going from 2 to 7.9 mm, indicating enhanced radial transport in the neon-seeded small ELM scenario. Inter-ELM λqs have been also calculated at the secondary outer divertor in quasi-double-null (QDN) discharges. It is found that, on average, λq is 2.2 times larger in the high-separatrix-density small ELM regime, compared to a reference type-I ELM one. These findings are supported by small ELMs radial velocity profiles, measured at the outer midplane with a fast reciprocating probe, showing a decay length (λvr) in the SOL of 12.8 cm, which is 3.3 times larger than that for the type-I ELMs of 3.9 cm. This analysis shows that small ELMs, although attractive for future machines due to low peak heat flux and large λq, might be of concern for the larger flux to the outer wall.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

WRF-ELM v1.0: a regional climate model to study land–atmosphere interactions over heterogeneous land use regions

Abstract. The Energy Exascale Earth System Model (E3SM) Land Model (ELM) is a state-of-the-art land surface model that simulates the intricate interactions between the terrestrial land surface and other components of the Earth system. Originating from the Community Land Model (CLM) version 4.5, ELM has been under active development, with added new features and functionality, including plant hydraulics, radiation–topography interaction, subsurface multiphase flow, and more explicit land use and management practices. This study integrates ELM v2.1 with the Weather Research and Forecasting (WRF; WRF-ELM) model through a modified Lightweight Infrastructure for Land Atmosphere Coupling (LILAC) framework, enabling affordable high-resolution regional modeling by leveraging ELM's innovative features alongside WRF's diverse atmospheric parameterization options. This framework includes a top-level driver for variable communication between WRF and ELM and Earth System Modeling Framework (ESMF) caps for the WRF atmospheric component and ELM workflow control, encompassing initialization, execution, and finalization. Importantly, this LILAC–ESMF framework demonstrates a more modular approach compared to previous coupling efforts between WRF and land surface models. It maintains the integrity of ELM's source code structure and facilitates the transfer of future developments in ELM to WRF-ELM. To test the ability of the coupled model to capture land–atmosphere interactions over regions with a variety of land uses and land covers, we conducted high-resolution (4 km) WRF-ELM ensemble simulations over the Great Lakes region (GLR) in the summer of 2018 and systematically compared the results against observations, reanalysis data, and WRF-CTSM (WRF coupled with the Community Terrestrial Systems Model). In general, the coupled WRF-ELM model has reasonably captured the spatial distribution of surface state variables and fluxes across the GLR, particularly over the natural vegetation areas. The evaluation results provide a baseline reference for further improvements in ELM in the regional application of high-resolution weather and climate predictions. Our work serves as an example to the model development community for expanding an advanced land surface model's capability to represent fully-coupled land–atmosphere interactions at fine spatial scales. The development and release of WRF-ELM marks a significant advancement for the ELM user community, providing opportunities for fine-scale regional representation, parameter calibration in coupled mode, and examination of new schemes with atmospheric feedback.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Toward holistic understanding of the ITER-like resonant magnetic perturbation (RMP) ELM control on KSTAR

Here, KSTAR has clarified a set of unresolved 3D physics issues utilizing the ITER-like in-vessel, three-row, resonant magnetic perturbation (RMP) configurations. Since RMP-driven, edge-localized-modes (ELMs)-crash control elevates the divertor heat flux peak through its impact on edge plasma parameters and transport, a series of intentionally misaligned RMP configurations (IMCs) have been explored to investigate the relationship between RMP ELM control and divertor heat fluxes, while searching for an ideal IMC that could be favorable in both aspects. First of all, the contrasting influence of kink vs anti-kink phasing on the ELM-crash suppression has been articulated, demonstrating the synergistic benefit of 'kink' phasing on ELM-crash-suppression. On the other hand, the three-row IMC in the anti-kink phasing becomes more insensitive to the ELM-crashes at the sub-marginal level of RMP, consistent with theory. Meanwhile, the divertor 'wetted' area of ELM-crash-suppression gets narrower than that of ELM-crash-mitigation, suggesting that ELM-crash-mitigation remains advantageous over ELM-crash-suppression in terms of time-averaged divertor thermal loading. In comparison, based on a set of two-row IMCs, no evidence of divertor heat flux broadening was found during ELM-crash-suppression, supporting a hypothesis that the dispersal of the divertor heat flux in three-row IMCs cannot be driven by helically structured two-row RMPs alone. Among ITER-like three-rows, lower two-row RMPs have been found to be much more effective in suppressing the ELM-crashes than upper two-row RMPs. Although it is quite preliminary, the up/down asymmetric dependence of RMP coupling may be generically attributed to lower-single-null plasmas. Such a holistic understanding of RMP-driven, ELM-crash-control in KSTAR is expected not only to elucidate various subtle points in the vicinity of ELM-crash-suppression, but also to clarify the relevant divertor thermal loading issues for ITER and beyond.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Investigation of pedestal parameters and divertor heat fluxes in small ELM regimes in DIII-D

Abstract Divertor heat flux and its correlation with pedestal parameters within various small edge localized mode (ELM) regimes, including high beta poloidal, type-II and ELMs with negative triangularity H-modes were investigated in DIII-D. The parallel energy fluences of type-II and high beta poloidal small ELM regimes fall below the linear scaling with pedestal electron pressure for type-I ELMs put forward in Eich et al 2017 ( Nucl. Mater. Energy 12 84–90). The negative triangularity of H-mode ELMs follow the Eich scaling for type-I ELMs. The parallel heat flux and total heat loads to the divertor were determined using high-time resolution infrared thermography, while pedestal parameters were obtained through self-consistent kinetic equilibrium reconstructions. Linear regressions for the type-II and high beta poloidal regimes demonstrate that an equivalent 7.5 MA small ELM scenario in ITER would fall below the ~5 MJ m − 2 leading edge melting limit for tungsten (Gunn et al 2017 Nucl. Fusion 57 046025). Utilizing fast thermography, the scrape-off layer power fall-off length for both inter-ELM and intra-ELM was determined and compared to the Eich scaling with poloidal magnetic field in Eich et al (ASDEX Upgrade Team and JET EFDA Contributors 2013 Nucl. Fusion 53 093031). Except for the high beta poloidal scenario, all the small ELM regimes during both inter- and intra-ELM periods had power fall-off lengths ( λ q ) larger then would be expected from the B pol , MP − 1 scaling associated with type-I ELMs, signifying their potential in managing heat loads and offering a solution for core–edge integration.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Exploring the transition from continuous turbulence fluctuations to bursting ELMs in high SOL density regimes

BOUT++ turbulence simulations of the DIII-D reveal that the density profile between the separatrix and pedestal plays a crucial role in the dynamics of edge localized modes (ELMs) and edge plasma turbulent transport. Nonlinear simulations demonstrate that small ELMs in the DIII-D hybrid scenario under high SOL density conditions are predominantly driven by local ballooning modes near the separatrix, stabilizing global instabilities while enhancing localized pressure fluctuations. A key control parameters for ELM dynamics is the separatrix-to-pedestal density ratio, n e,sep /n e,ped . A high ratio indicates a shallow gradient, favoring small ELMs, while a lower ratio signals a steep gradient, which increases the likelihood of large ELMs. Comprehensive parameter scans, including n e,sep /n e,ped , density gradient profiles near the separatrix, and resistivity, reveal the critical role of these parameters in shaping transitions between turbulence-driven transport and ELM bursting. The scans demonstrate that in high SOL density regimes, small ELMs can result from either global resistive MHD instabilities or local ballooning modes near the separatrix, depending on the steepness of the separatrix density gradient. These findings also highlight the transition from continuous turbulence to small ELMs. The post-crash peak in pressure fluctuations, δP rms serves as a critical metric for identifying transition from continuous turbulence fluctuations to ELM bursting. Larger δP rms values correlate with ELM bursts driven by local or global instabilities, whereas smaller values indicate turbulence-dominated transport. Drift-Alfvén and resistive ballooning turbulence enhance the entrainment of fluctuations from the pedestal to the SOL, contributing to the complex interplay of dynamics in this regime. These findings emphasize the importance of separatrix density shaping and pedestal gradient control for optimizing ELM behavior in ITER and future fusion devices.

Li, Nami [Lawrence Livermore National Laboratory (↗
DOE OSTI2020

Grassy ELM regime at low pedestal collisionality in high-power tokamak plasma

Spontaneous mitigation of giant ELMs and appearance of grassy ELMs have been observed repeatedly at low pedestal collisionality (${\nu }_{\mathrm{e},\mathrm{p}\mathrm{e}\mathrm{d}}^{{\ast}}\sim 0.15$) in the high-power (P inj > 13 MW) hybrid scenario in DIII-D. Higher β p and higher q 95 appear to be beneficial to achieving the grassy ELM regime. The grassy ELM H-mode plasma shows high energy confinement performance (H 98y2 up to over 1.5) under the conditions of high neutral beam torque and high core rotation. The pedestal width appears to exceed the EPED1.0 model prediction by more than 50%. Pedestal stability analysis performed with the ELITE code indicates that the stability against low-n kink/peeling modes is improved with increased auxiliary heating power and the operational point in the grassy ELM regime is located near the ballooning boundary. The pedestal stability characteristics during the grassy ELM crashes have been investigated in comparison with the giant ELM crashes based on plasma profiles experimentally measured with high resolution and accuracy. It has been found that the underlying mechanism for the observed small-amplitude ELM crashes is mainly the expansion of the ballooning stability boundary induced by an initial radially localized collapse in the pedestal, which helps to stop the growth of instabilities and further collapse of the pedestal. The effect of electron density pedestal on mitigating edge localized instabilities has been analyzed by numerical simulation, suggesting that the electron density pedestal characterized by high n e,sep /n e,ped and low density gradient helps to stabilize peeling-ballooning modes because of a low pressure pedestal gradient and to lower the ballooning boundary mainly because of a low ion diamagnetic frequency in the pedestal region, thus triggering ballooning instabilities and producing the intrinsic grassy ELMs. Numerical simulation of the Chinese fusion engineering test reactor (CFETR) with the SOLPS code indicates that the separatrix density might be insensitive to the electron diffusivity in the pedestal region and increase with the power flowing from the core region to the edge region. Furthermore, pedestal stability analysis suggests that the flat density pedestal with high separatrix density obtained in the high-power plasma in CFETR would make the operational point close to the ballooning boundary, which is considered to help destabilize ballooning instabilities and facilitate the access to the grassy ELM regime.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Prediction of pellet mass thresholds for ELM triggering in low-collisionality, ITER-like discharges

Abstract In ITER, pellets are calculated to require more than 8 times the mass than currently planned to reliably trigger edge-localized modes (ELMs). Unmitigated heat flux impulses from ELMs are intolerable in ITER at full power and current. Therefore, ITER operation relies on multiple approaches to control ELM heat fluxes. One method is pellet ELM pacing to instigate small rapid ELMs with low heat flux. Predicting the performance of pellet pacing is critical for ITER, which is expected to operate in a regime with a low-collisionality, peeling-limited pedestal. However, to trigger ELMs the local pressure increase in the expanding pellet cloud pushes the equilibrium over the ballooning stability limit. In this work, linear and nonlinear M3D-C1 simulations are used to predict pellet mass thresholds in DIII-D discharges and ITER scenarios with peeling-limited pedestals. It is found that the distance of the equilibrium’s operational point from the ballooning branch of the pedestal stability boundary strongly changes thresholds. Linear M3D-C1 simulations find a strong dependence of the pellet mass threshold on the poloidal injection location for ITER’s 15 MA, Q = 10 scenario. The required pellet mass at the planned injection locations is 8 to 17 times larger than currently considered. However, such linear simulations do not include pellet ablation physics or time evolution of density and temperature. A new scheme of 2D nonlinear simulations, coupled with linear stability analysis at various steps throughout the nonlinear time evolution, was developed to include such physics and improve on the linear results. These new nonlinear-to-linear simulations confirm previous findings. This result suggests that pellet ELM triggering in ITER could require pellets much larger than those currently planned, which makes ELM-pacing operationally challenging. On the other hand, fueling pellets injected from the high-field side will likely not unintentionally trigger ELMs in an otherwise ELM-stable plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗