Search NASA⌕ Search

SEARCH · Search NASA

Results for “Separatrix density”

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

Effect of the divertor configuration on the JET edge radial electric field

The role of the divertor configuration and divertor plasma physics on the L–H transition is poorly understood, leading to large uncertainties in predicting the L–H power threshold in future devices. This contribution reports on edge perpendicular plasma flow measurements by Doppler backscattering in JET L–H transition experiments with the outer divertor strike-point at different positions: horizontal target (HT), vertical target (VT), and in the corner configuration (between the HT and VTs). The edge perpendicular flow was found to be significantly affected by changes in the divertor configuration in the region inside the separatrix. Our results do not show evidence for the existence of a critical edge flow shear needed to achieve H-mode for different divertor configurations, with a larger shear observed for the VT configuration. No significant change in the shear flow and in the density fluctuation level is measured preceding the L–H transition in the region just inside the separatrix. The dynamics of the L–H transitions are also influenced by the divertor configuration with divertor oscillations (DOs) observed only in the HT configuration. Interestingly, DOs are associated with marked changes in the edge perpendicular flow around the separatrix.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

2D analysis of tokamak divertor-plasma detachment-bifurcation with operational parameters and geometries

UEDGE simulations with density scans for various input power, transport coefficients and outer poloidal leg length are performed to study the conditions for the existence of a bifurcation-like drop of T e at the outer strike point, commonly referred to as a detachment cliff, when transitioning to a detached plasma from an attached plasma in the outer divertor as the upstream density increases (McLean et al., 2015). The simulation results show that a detachment cliff tends to occur with a higher power input regardless of diffusivities and leg length. Further analysis of change of plasma profiles at a cliff indicate that, in addition to the sharp reduction of the E x B drift fluxes in the outer divertor studied in Jaervinen et al., (2018), the substantial change of the Mach number in the outer divertor and the decrease of the outer mid-plane T e due to the radiation front moving across the separatrix into the confinement region above the X-point consistently occur for all UEDGE density scans that have a detachment cliff. UEDGE time-dependent simulation of the evolution of a detachment cliff shows that the rapid increase of radiation above the X-point occurs in a time scale of ~0.3–0.5, which could possibly be the trigger for the formation of a detachment cliff, quicker than the Mach number change in a time scale of ~1 ms and the drop of T e in a time scale of ~2–3 ms in the outer divertor.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Detachment in Fusion Plasmas with Symmetry Breaking Magnetic Perturbation Fields

Power exhaust from the bulk plasma is significantly altered by symmetry breaking magnetic perturbation fields, because these create direct connections (perturbed field lines) from the confined high temperature plasma to solid surfaces. Additionally, the same amount of power is distributed among those new exhaust channels as for a symmetric magnetic configuration, which reduces the local upstream heat flux flowing down the perturbed field lines, thereby making access to detachment easier (i.e., at lower upstream density) for the divertor plasma near the location corresponding to the symmetric magnetic separatrix. However, the divertor plasma regions with connection to the bulk plasma are extended nonaxisymmetrically further outside, where significant heat loads occur, unlike in the symmetric configuration. The temperature remains high at those locations, which reduces the divertor plasma dissipation capacity, making the mitigation of heat loads more difficult to achieve.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Prediction of DIII-D Pedestal Structure From Externally Controllable Parameters

The sharp increase of pressure at the edge of a high confinement mode (H-mode) plasma, the pedestal, strongly impacts overall plasma performance. Predicting the pedestal is a necessity to control and optimize tokamak operations. Here, an experimental data-driven machine learning (ML) approach is presented that predicts the pedestal heights and widths of electron density (n e ) and electron temperature (T e ) profiles as well as the separatrix ne from externally controllable parameters such as the plasma shape, heating method and power, and gas puff rate and integrated gas puff. The OMFIT framework was used with DIII-D data to efficiently, robustly, and automatically build a database of pedestal parameters to train machine learning models. Database creation was enabled by the search engine tool for DIII-D data, TokSearch, which parallelizes data fetching, enabling fast searches through basic signals of thousands of DIII-D shots and selection of relevant time intervals. Principal Component Analysis (PCA) separated the database into three clusters that represent classes of plasma shapes that are regularly used in DIII-D. The most important parameters for setting the pedestal structure were plasma current (I p ), toroidal magnetic field (B Φ ), neutral beam heating power (P NBI ) and shaping quantities. The Deep Jointly Informed Neural Networks (DJINN) algorithm was applied to identify suitable neural network (NN) architectures that appropriately capture the features of the pedestal database. Separate NNs were implemented for each pedestal parameter, and ensembling methods were used to improve the prediction accuracy and allowed estimation of the prediction uncertainty. The pedestal predictions of the test dataset lie within the measurement uncertainties of the pedestal parameters. The NN outperformed simple Linear Regression (LR) analysis, indicating non-linear dependencies in the pedestal structure. The presented achievements illustrate a promising path for future research, using feature extraction to infer experimental trends and thereby improve pedestal models as well as deploying NN for a fast pedestal prediction in DIII-D scenario development.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Optimization of the N ∥ Upshift in the DIII-D high field side lower hybrid current drive experiment

High field side lower hybrid current drive (LHCD) is one potential candidate for efficient non-inductive current drive in tokamak power plants, and the first test of this technology will occur on the DIII-D tokamak during the 2024 campaign. Previous LFS launch experiments operated in the multi-pass regime and relied on scrape-off layer interactions to close the spectral gap. In the DIII-D experiment, single-pass damping is achievable via an upshift in the parallel refractive index N ∥ caused by mode converting twice (slow → fast → slow). This mode conversion affects the ray trajectories and can lead to enhanced N ∥ upshift depending on where mode conversion occurs. Compared to multi-pass absorption experiments, the optimization of launched N ∥ and plasma parameters can be counter-intuitive: increased density may increase efficiency and smaller N ∥,launch tend to damp closer to the separatrix. A hard x-ray camera installed to measure the bremsstrahlung (50–250 keV) radiation from LHCD-generated fast electrons is capable of verifying the trends reporting in this paper through comparison to the ray-tracing/Fokker–Planck codes GENRAY/CQL3D.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Separating divertor closure effects on divertor detachment and pedestal shape in DIII-D

Comparison between an open divertor and a more-closed divertor in DIII-D demonstrates detachment at an up to 40% lower pedestal density (n e,ped ) in the closed divertor, due to a combination of decreased fueling of the pedestal and increased dissipation in the scrape off layer (SOL) in the closed divertor, both resulting from increased neutral trapping in the divertor. Predicting whether the relationship between divertor closure and detachment will hold for an opaque SOL, in which the contribution of ionizing neutrals to fueling the pedestal is lessened, requires separating out the different mechanisms contributing to the density difference at detachment. A series of experiments on DIII-D characterizes matched discharges using various divertor configurations to isolate the effects of divertor closure. These experiments show detachment at up to 25% lower n e,sep in the closed divertor than the open, supported by simulations showing increased neutral trapping, and hence increased dissipation, in the closed divertor. A difference in n e,ped /n e,sep is also seen: for matched n e,sep , the closed divertor has up to 20% lower ne,ped, consistent with modeling showing a smaller ionization fraction inside the separatrix in this case. Understanding how these pieces fit together will help in develop of predictive models of pedestal density and detached divertors compatible with a high-performance core.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Electron transport barrier and high confinement in configurations with internal islands close to the plasma edge of W7-X

The low magnetic shear in the Wendelstein 7-X (W7-X) stellarator makes it feasible to shape the separatrix by the large islands constituting an island-divertor, and this can be exploited to access various magnetic configurations, including samples of different internal island sizes and locations. To investigate the configuration effects on the plasma confinement, a configuration scan was performed by changing the coil currents to vary the rotational transform between values 5/4 and 5/6 at the plasma boundary with different power levels (2, 4, 6 MW) of electron cyclotron resonance heating (ECRH) at a maximum plasma density of 8 × 10 19 m -2 . neutral beam injection (NBI) heating was also applied during some configurations of the scan to create a density ramp and access high densities beyond the X2 ECRH cutoff. For the magnetic configurations, where the 5/5 and 5/6 island chains were moved closer to separatrix but remaining inside the last closed flux surface, the electron cyclotron emission shows that an electron temperature, T e , pedestal develops already during ECRH heated plasma buildup phase indicating a transport barrier, and the barrier sustains irrespective of changed plasma heating conditions such as NBI in the later part of discharge. The transport barrier is broken by subsequent fast crashes, observed with multiple plasma diagnostics with characteristics such as tokamak edge localized modes, and the corresponding crash amplitude and frequency vary with plasma pressure. The impact of the transport barrier on plasma confinement can be seen through the increased core T e profile, which could be responsible for the overall increase in the stored diamagnetic energy by approximately 10% for these configurations. After the plasma heating is terminated, a backwards transition to a degraded confinement state is also observed. These observations indicate a configuration triggered high confinement mode in low shear W7-X. This work focuses on the occurrence of this transport barrier for different magnetic configurations and its relation to internal magnetic islands.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Enhanced pedestal transport driven by edge collisionality on Alcator C-Mod and its role in regulating H-mode pedestal gradients

Experimental measurements of plasma and neutral profiles across the pedestal are used in conjunction with 2D edge modeling to examine pedestal stiffness in Alcator C-Mod H-mode plasmas. Enhanced D α experiments on Alcator C-Mod observed pedestal degradation and loss in confinement below a critical value of net power crossing the separatrix, P net = $P^{crit}_{net}$ ≈ 2.3 MW, in the absence of any external fueling. New analysis of ionization and particle flux profiles reveal saturation of the pedestal electron density, $n^{ped}_{e}$, despite continuous increases in ionization throughout the pedestal, inversely related to P net . A limi to the pedestal $\nabla$n e emerges as the particle flux, Γ D , continues to grow, implying increases in the effective particle diffusivity, D eff . This is well-correlated with the separatrix collisionality, $v^*_{sep}$ and a turbulence control parameter, α t , implying a possible transition in type of turbulence. The transition is well correlated with the experimentally observed value of $P^{crit}_{net}$. SOLPS-ITER modeling is performed for select discharges from the power scan, constrained with experimental electron and neutral densities, measured at the outer midpane. The modeling confirms general growth in D eff , consistent with experimental findings, and additionally suggests even larger growth in Χ e at the same $P^{crit}_{net}$.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Isotope dependence of the type I ELMy H-mode pedestal in JET-ILW hydrogen and deuterium plasmas

The pedestal structure, edge transport and linear MHD stability have been analyzed in a series of JET with the ITER-like wall hydrogen (H) and deuterium (D) type I ELMy H-mode plasmas. The pedestal pressure is typically higher in D than in H at the same input power and gas rate, with the difference mainly due to lower density in H than in D (Maggi et al (JET Contributors) 2018 Plasma Phys. Control. Fusion 60 014045). A power balance analysis of the pedestal has shown that higher inter-ELM separatrix loss power is required in H than in D to maintain a similar pedestal top pressure. This is qualitatively consistent with a set of interpretative EDGE2D-EIRENE simulations for H and D plasmas, showing that higher edge particle and heat transport coefficients are needed in H than in D to match the experimental profiles. It has also been concluded that the difference in neutral penetration between H and D leads only to minor changes in the upstream density profiles and with trends opposite to experimental observations. This implies that neutral penetration has a minor role in setting the difference between H and D pedestals, but higher ELM and/or inter-ELM transport are likely to be the main players. The interpretative EDGE2D-EIRENE simulations, with simultaneous upstream and outer divertor target profile constraints, have indicated higher separatrix electron temperature in H than in D for a pair of discharges at low fueling gas rate and similar stored energy (which required higher input power in H than in D at the same gas rate). The isotope dependence of linear MHD pedestal stability has been found to be small, but if a higher separatrix temperature is considered in H than in D, this could lead to destabilization of peeling-ballooning modes and shrinking of the stability boundary, qualitatively consistent with the reduced pedestal confinement in H.

08 HYDROGEN↗

Overview of density pedestal structure: Role fueling versus transport

How much of the density pedestal structure profile is determined by fueling versus transport is still a valid question. In most current tokamaks, the Scrape-Off Layer is not opaque enough to screen the neutrals and limit their penetration past the separatrix. As such recycled neutrals as well as additional gas puffing can directly affect the pedestal structure through ionization. Yet, in ITER and other future burning plasma magnetic confinement devices, this ionization source will be pushed further out into the Scrape-Off Layer. This poses the question, how much of the pedestal density structure is governed by this edge ionization source versus plasma transport effects. Theoretically, several turbulent modes have been identified which could provide a 'pinch' like up-gradient transport mechanism. Up-gradient transport is necessary in a source-free region to obtain a peaked density profile, which is often observed in the core of tokamaks. In an opaque Scrape-Off Layer, without the existence of such a pinch, the pedestal density structure would eventually disappear. Not only would this limit our ability to fuel a fusion reactor, it would also affect the pedestal stability. In this paper we will give an overview of the state of the art on what sets the pedestal density structure based on experimental observations as well as theoretical models. We will complement these studies with new results from Alcator C-Mod and DIII-D where the opaqueness is increased to reach values similar to those in ITER.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The effect of plasma shaping on high density H-mode SOL profiles and fluctuations in TCV

Abstract The impact of plasma shaping on the properties of high density H-mode scrape-off layer (SOL) profiles and transport at the outer midplane has been investigated on Tokamakà configuration variable. The experimental dataset has been acquired by evolving the upper triangularity while keeping the other parameters constant. The scan comprises δ u p values between 0.0 and 0.6, excluding negative triangularity scenarios. Within this study, a transition from type-I edge localised modes to the quasi-continuous exhaust regime takes place from low to high δ u p . The modification of the upstream SOL profiles has been assessed, in terms of separatrix quantities, within the α t turbulence control parameter theoretical framework (Eich et al 2020 Nucl. Fusion 60 056016). The target parallel heat load and the upstream near-SOL density profiles have been shown to broaden significantly for increasing α t . Correspondingly, in the far SOL a density shoulder formation is observed when moving from low to high δ u p . These behaviours have been correlated with an enhancement of the SOL fluctuation level, as registered by wall-mounted Langmuir probes as well as the thermal helium beam diagnostic. Specifically, both the background and the filamentary-induced fluctuating parts of the first wall ion saturation current signal are larger at higher δ u p , with filaments being ejected more frequently into the SOL. Comparison of two pulses at the extremes of the δ u p scan range, but with otherwise same input parameters, shows that the midplane neutral pressure does not change much during the H-mode phase of the discharge. This indicates that indirect effects of the change in geometry, linked to first wall recycling sources, should not play a significant role. The total core radiation increases at high δ u p , on account of a stronger plasma–wall interaction and resulting larger carbon impurity intake from the first wall. This is likely associated to the enhanced first wall fluctuations, as well as a smaller outer gap and the close-to-double-null magnetic topology at high shaping.

Stagni, A. (ORCID:0000000180841544)↗

Plasma and Energetic Particle Behaviors During Asymmetric Magnetic Reconnection at the Magnetopause

The factors controlling asymmetric reconnection and the role of the cold plasma population in the reconnection process are two outstanding questions. We present a case study of multipoint Cluster observations demonstrating that the separatrix and flow boundary angles are greater on the magnetosheath than on the magnetospheric side of the magnetopause, probably due to the stronger density than magnetic field asymmetry at this boundary. The motion of cold plasmaspheric ions entering the reconnection region differs from that of warmer magnetosheath and magnetospheric ions. In contrast to the warmer ions, which are probably accelerated by reconnection in the diffusion region near the subsolar magnetopause, the colder ions are simply entrained by 𝐄×𝐁 drifts at high latitudes on the recently reconnected magnetic field lines. This indicates that plasmaspheric ions can sometimes play only a very limited role in asymmetric reconnection, in contrast to previous simulation studies. Three cold ion populations (probably H+, He+, and O+) appear in the energy spectrum, consistent with ion acceleration to a common velocity.

Cluster↗

Electron Densities in Solar Flare Loops, Chromospheric Evaporation Upflows, and Acceleration Sites

We compare electron densities measured at three different locations in solar flares: (1) in Soft X-Ray (SXR) loops, determined from SXR emission measures and loop diameters from Yohkoh Soft X-Ray Telescope maps (n(sub e, sup SXR) = (0.2-2.5) x 10(exp 11)/ cu cm); (2) in chromospheric evaporation upflows, inferred from plasma frequency cutoffs of decimetric radio bursts detected with the 0.1-3 GHz spectrometer Phoenix of ETH Zuerich (n(sub e, sup upflow) = (0.3-11) x 10(exp 10)/cu cm; and (3) in acceleration sites, inferred from the plasma frequency at the separatrix between upward-accelerated (type III bursts) and downward-accelerated (reverse-drift bursts) electron beams [n(sub e, sup acc) = (0.6-10) x 10(exp 9)/cu cm]. The comparison of these density measurements, obtained from 44 flare episodes (during 14 different flares), demonstrates the compatibility of flare plasma density diagnostics with SXR and radio methods. The density in the upflowing plasma is found to be somewhat lower than in the filled loops, having ratios in a range n(sub e, sup upflow)/n(sub e, sup SXR) = 0.02-1.3, and a factor of 3.6 higher behind the upflow front. The acceleration sites are found to have a much lower density than the SXR-bright flare loops, i.e., n(sub e, sup acc)/n(sub e, sup SXR) = 0.005- 0.13, and thus must be physically displaced from the SXR-bright flare loops. The scaling law between electron time-of-flight distances l' and loop half-lengths s, l'/s = 1.4 +/- 0.3, recently established by Aschwanden et al. suggests that the centroid of the acceleration region is located above the SXR-bright flare loop, as envisioned in cusp geometries (e.g., in magnetic reconnection models).

Aschwanden, Markus J.↗

Impurity leakage and radiative cooling in the first nitrogen and neon seeding study in the closed DIII-D SAS configuration

A comparative study of nitrogen versus neon has been carried out to analyze the impact of the two radiative species on power dissipation, SOL impurity distribution, divertor and pedestal characteristics. The experimental results show that N remains compressed in the divertor, thereby providing high radiative losses without affecting the pedestal profiles and displacing carbon as dominant radiator. Neon, instead, radiates more upstream than N thus reducing the power flux through the separatrix leading to a reduced ELM frequency and compression in the divertor. A significant amount of neon is measured in the plasma core leading to a steeper density gradient. The different behavior between the two impurities is confirmed by SOLPS-ITER modeling which for the first time at DIII-D includes multiple impurity species and a treatment of full drifts, currents and neutral–neutral collisions. The impurity transport in the SOL is studied in terms of the parallel momentum balance showing that N is mostly retained in the divertor whereas Ne leaks out consistent with its higher ionization potential and longer mean free path. This is also in agreement with the enrichment factor calculations which indicate lower divertor enrichment for neon. The strong ionization source characterizing the SAS divertor causes a reversal of the main ions and impurity flows. The flow reversal together with plasma drifts and the effect of the thermal force contribute significantly in the shift of the impurity stagnation point affecting impurity leakage. This work provides a demonstration of the impurity leakage mechanism in a closed divertor structure and the consequent impact on pedestal. Since carbon is an intrinsic radiator at DIII-D, in this paper we have also demonstrated the different role of carbon in the N vs Ne seeded cases both in the experiments and in the numerical modeling. Here, carbon contributes more when neon seeding is injected compared to when nitrogen is used. Finally, the results highlight the importance of accompanying experimental studies with numerical modeling of plasma flows, drifts and ionization profile to determine the details of the SOL impurity transport as the latter may vary with changes in divertor regime and geometry. In the cases presented here, plasma drifts and flow reversal caused by high level of closure in the slot upper divertor at DIII-D play an important role in the underlined mechanism.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Near- and far scrape-off layer transport studies in detached, small-ELM ASDEX Upgrade discharges by means of EMC3-EIRENE

We compare various diagnostics characterizing an ASDEX Upgrade upper single-null discharge to EMC3-EIRENE simulations now including volumetric recombination and main chamber plasma-wall interaction but not yet drifts. The discharge is in a small-ELM regime and is approaching detachment due to a density ramp and nitrogen seeding. Time/ELM averaged upstream density, electron and ion temperature measurements as well as downstream ion saturation current, electron temperature and power flux measurements match the simulation within the error bars before and after reaching detachment, qualifying the code for this regime. Although the Ion Cyclotron Resonance Heating (ICRH) antenna limiters are 35 mm away from the separatrix, i.e. ten times the typical near-SOL power fall-off length $\lambda^{near}_q\ =\ $3.8 mm, they receive power fluxes of several hundreds of kW exceeding the ones to the outer target in the detached phase and strongly contribute to the particle fueling. This is explained by a substantially larger far-SOL power fall-off length $\lambda^{far}_q\ =\ $25 mm possibly caused by the small ELMs in the detached phase and the occurrence of a density shoulder.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Dependence of scrape-off layer profiles and turbulence on gas fuelling in high density H-mode regimes in TCV

A set of high density, highly shaped H-mode discharges has been performed in the TCV tokamak with the aim of assessing the effects of increasing divertor neutral recycling on the properties of upstream inter-ELM scrape-off layer (SOL) profiles and transport. An increase of divertor neutral pressure has been correlated with the evolution of separatrix properties and turbulence level. The latter has been quantified by means of the α t parameter introduced in (Eich 2020 Nucl. Fusion 60 056016), describing the contribution of resistive-interchange turbulence in the SOL relative to drift wave transport. The analysis reveals a general broadening of the upstream SOL profiles as α t increases, with the SOL power width measured by the vertical IR thermography system increasing significantly. In a similar way, the upstream density profile widens in the near SOL, whereas in the far SOL a density shoulder is observed to progressively form and increase in amplitude. This behaviour is associated with an enhancement of far SOL turbulent transport in the form of blob-filaments travelling radially faster across the far SOL and becoming bigger at higher α t . Finally, the detected filaments, evaluated from the fast reciprocating probe at the outer midplane, are determined to mostly belong to the resistive ballooning and resistive X-point regimes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Diagnosing the pedestal magnetic field and magnetohydrodynamics radial structure with pedestal–scrape of layer electron cyclotron emission radiation inversion in H-mode plasma (invited)

Forward modeling is used to interpret inversion patterns of the pedestal–Scrape of Layer (SOL) Electron Cyclotron Emission (ECE) in DIII-D H-mode experiments. The modeling not only significantly improves the ECE data interpretation quality but also leads to the potential measurements of (1) the magnetic field strength |B| at the separatrix, (2) the pedestal |B| evolution during an inter-Edge Localized Mode (ELM) period, and (3) the pedestal Magnetohydrodynamics (MHD) radial structure. The ECE shine-through effect leads to three types of pedestal–SOL radiation inversions that are discussed in this paper. The first type of inversion is the non-monotonic T e,rad profile with respect to the major radius. Using the ECE frequency at the minimum T e,rad , the inversion can be applied to measure the magnetic field |B| at the separatrix and calibrate the mapping of the ECE channels with respect to the separatrix. The second type of inversion refers to the opposite phase between the radiation fluctuations δT e,rad at the pedestal and SOL. This δT e,rad phase inversion is sensitive to density and temperature fluctuations at the pedestal foot and, thus, can be used to qualitatively measure the MHD radial structure. The third type of inversion appears when the pedestal and SOL T e,rad evolve in an opposite trend, which can be used to infer the pedestal |B| field change during an inter-ELM period. Finally, the bandwidth effect on measuring δT e,rad due to pedestal MHD is also investigated in the radiation modeling.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Type-I ELM mitigation by continuous lithium granule gravitational injection into the upper tungsten divertor in EAST

Large edge-localized modes (ELMs) were mitigated by gravitational injection of lithium granules into the upper X-point region of the experimental advanced superconducting tokamak (EAST) device with tungsten plasma-facing components. The maximum ELM size was reduced by ~70% in high β N H-mode plasmas. Large ELM stabilization was sustained for up to about 40 energy confinement times, with constant core radiated power and no evidence of high-Z or low-Z impurity accumulation. The lithium granules injection reduced the edge plasma pedestal density and temperature and their gradients, due to increased edge radiation and reduced recycling from the plasma-facing components. Ideal stability calculations using the ELITE code indicate that the stabilization of large ELMs correlates with improved stability of intermediate-n peeling-ballooning modes, due to reduced edge current resulting from the profile changes. The pedestal pressure reduction was partially offset by a core density increase, which resulted in a modest ~7% drop in core stored energy and normalized energy confinement time. We surmise that the remnant small ELMs are triggered by the penetration of multiple Li granules just past the separatrix, similar to small ELMs triggered by deuterium pellet Futatani et al (2014 Nucl. Fusion 54 073008). This study extends previous ELM elimination with Li powder injection Maingi et al (2018 Nucl. Fusion 58 024003) in EAST because (1) use of small, dust-like powder and the related potential health hazards were eliminated, and (2) use of macroscopic granules should be more applicable to future devices, due to deeper penetration than dust particles, e.g. inside the separatrix with velocities ~10 m s -1 in EAST.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗