Engineering topics

Grierson, B. A.

Publications and source records attributed to Grierson, B. A..

32 records · Page 2

DOE OSTI2021

Radially resolved active charge exchange measurements of the hydrogenic isotope fraction on DIII-D

Radially resolved hydrogenic isotope fraction measurement capabilities have been developed for DIII-D using the main-ion charge exchange recombination (MICER) spectroscopy system in preparation for mixed hydrogen and deuterium experiments. Constraints on the hydrogenic ion temperatures and velocities based on measurements of the impurity ion properties are required to accurately fit the spectrum. Corrections for cross sectional distortions, spatial smearing due to the halo, and a neoclassical offset between the impurity and hydrogenic toroidal rotation are applied to the constraints prior to fitting the MICER spectrum. Extensive atomic physics calculations have been performed using the FIDASIM code, which has recently been improved to allow simulations using mixtures of hydrogenic species. Overall, these results demonstrate that for the same plasma parameters, the Dα emission is 20%–30% brighter than Hα due to differences in rate coefficients associated with the different ion thermal velocities for the same temperature and therefore must be taken into consideration when calculating absolute densities. However, despite these differences, the absolute error when estimating the hydrogen isotope fraction [n H /(n H + n D )] by using the Hα radiance fraction [L Hα /(L Hα + L Dα )] is typically less than 5% due to the way the fraction is formed, making the radiance fraction a reasonably accurate estimate of the isotope fraction for most cases.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Testing predictions of electron scale turbulent pedestal transport in two DIII-D ELMy H-modes

In this study, we present linear and nonlinear gyrokinetic analyses in the pedestal region of two DIII-D ELMy H-mode discharges using the CGYRO code. The otherwise matched discharges employ different divertor configurations to investigate the impact of varying recycling and particle source on pedestal profiles. Linear gyrokinetic simulations find electrostatic ion-scale instabilities (ion temperature gradient and trapped electron modes, ITG-TEM) are present just inside the top of the pedestal with growth rates that are enhanced significantly by parallel velocity shear. In the sharp gradient region, E x B shearing rates are comparable or larger than ion scale growth rates, suggesting the suppression of ITG-TEM modes in this region. Instead, the electron temperature profiles are found to be correlated with and just above the electron temperature gradient (ETG) instability thresholds. Using gradients varied within experimental uncertainties, nonlinear electron-scale gyrokinetic simulations predict electron heat fluxes from ETG turbulence, that when added to neoclassical ion thermal transport simulated by NEO, account for 30-60% of the total experimental heat flux. In addition, the neoclassical electron particle flux is found to contribute significantly to the experimental fluxes inferred from SOLPS-ITER analysis. Additional nonlinear gyrokinetic simulations are run varying input gradients to develop a threshold-based reduced model for ETG transport, finding a relatively simple dependence on η e = L ne /L Te . Predictive transport simulations are used to validate this pedestal-specific ETG model, in conjunction with a model for neoclassical particle transport. In both discharges, the predicted electron temperatures are always overpredicted, indicative of the insufficient stiffness in the ETG pedestal model to account for all of the experimental electron thermal transport. In the case of the closed divertor discharge with lower particle source, the predicted electron density is close to the experiment, consistent with the magnitude of neoclassical particle transport in that discharge. However, the density profiles are overpredicted in the open divertor discharge (larger particle source), due to insufficient model transport. In conclusion, the implications for other mechanisms accounting for the remainder of transport in the sharp gradient region in the two discharges are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Charge exchange recombination spectroscopy measurements of DIII-D poloidal rotation with poloidal asymmetry in angular rotation

Sixteen new tangential views for the charge exchange recombination (CER) spectroscopy diagnostic at DIII-D were installed in 2019 on the high-field side (HFS) of the tokamak with the main goal being the measurement of main-ion (deuterium) poloidal rotation. Eight of the new views are connected to spectrometers, which view the main-ion spectrum, adding main-ion measurements where there were previously none, and another eight new views increased the spatial resolution of existing impurity (carbon) measurements on the HFS. When combined with the existing low-field side measurements, measurements at two locations on flux surfaces out to a normalized minor radius of ≈0.6 are possible. The new tangential views have been used to measure the deuterium poloidal rotation directly for the first time using the Poloidal Asymmetry in Angular Rotation (PAAR) method. Furthermore, these new measurements enable further testing of the validity of neoclassical poloidal rotation predictions. Separate measurements of the radial electric field can be made for an impurity ion and the main-ion by combining the PAAR measurements with additional CER measurements of toroidal rotation, temperature, and density. These independent measurements of the radial electric field agree reasonably well.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Pedestal collapse by resonant magnetic perturbations

Pedestal collapse (i.e., the complete loss of the edge transport barrier (ETB)) in DIII-D H-mode plasmas occurs when resonant magnetic perturbations (RMPs) penetrate the steep gradient region at the plasma edge. Normally, RMP driven magnetic islands can occur at the top and bottom of the H-mode pedestal and these islands generate conditions consistent with edge-localized-mode (ELM) suppression and density pump-out, respectively, based on nonlinear two-fluid MHD simulations. In contrast, MHD simulations show that the steep pressure gradient region between the top and bottom of the DIII-D pedestal is generally immune to resonant field penetration due to large local E × B and diamagnetic flows. By this fortuitous circumstance, the edge-transport-barrier and H-mode confinement can be maintained while achieving ELM suppression. However, pedestal collapse can occur in DIII-D when the screening flows are inadequate to prevent field penetration in the steep gradient region of the pedestal. Non-linear two-fluid MHD simulations support the role of resonant field penetration in pedestal collapse for DIII-D H-mode plasmas with weak edge E × B and diamagnetic screening flows. ITER will likely have weaker edge screening flows than present experiments due to its much larger size, making it more susceptible to resonant field penetration in the steep gradient region of the pedestal. Analysis of model ITER equilibria demonstrates that resonant field penetration in the steep pressure gradient region is possible for RMP levels of the order required for ELM suppression. The effect of such penetration on the ITER pedestal will depend sensitively on the resulting degree of island overlap.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2021

Response of thermal and fast-ion transport to beam ion population, rotation and T e / T i in the DIII-D steady state hybrid scenario

The thermal and fast-ion transport properties of DIII-D steady-state hybrid discharges with normalized beta β N ≳ 3 are studied at low injected torque and an increased electron to ion temperature ratio T e / T i . Linear stability analysis performed with the TGLF turbulent code indicates that a high-k mode is usually dominant at smaller radii, whereas a low-k mode is usually dominant at larger radii in these plasmas. A reduction in the net injected torque from 8.6 to 4.3 N-m leads to reduced E × B shear and hence, an enhanced turbulence that was observed on the Doppler backscattering diagnostic and was also computed with TGLF. As T e / T i in the core was increased from 0.57 to 0.66 by adding electron cyclotron current drive (ECCD) to these plasmas, higher levels of transport are observed with increased high-k modes indicated by TGLF. The fast-ion transport level varied over an order of magnitude in these discharges depending on whether Alfvén eigenmodes, fishbones, or no instabilities were observed. Hybrid plasmas with fishbones have decreased fast-ion transport, compared to plasmas with Alfvén eigenmodes, since they are resonant with a smaller portion of phase space and their resonance is farther from the wall. This reduction in fast-ion transport with ECCD mitigates the increase in turbulent transport, resulting in higher performance than expected during strong electron heating. Similarly, the lowest fast-ion transport was observed in the low torque plasma, which also led to better than expected performance at this torque value. Finally, the thermal and fast-ion transport changes observed as the torque/rotation and T e / T i are varied indicate possible methods for transferring this scenario to a reactor.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

Magnetic shear effect on plasma transport at T e / T i ~ 1 through electron cyclotron heating in DIII-D plasmas

The effect of magnetic shear on plasma transport for an electron to ion temperature ratio ( T e / T i ) near unity has been explored in DIII-D utilizing electron cyclotron heating (ECH). Previous reports showed that significant confinement degradation occurred at T e / T i ~ 1 in positive shear (PS) plasmas in DIII-D, whereas reduced confinement degradation was observed in negative central shear (NCS) plasmas. In this study, plasma transport in weak magnetic shear (WS) plasmas with ECH is investigated and compared with that in NCS and PS plasmas. Here the magnetic shears ($\hat s$) are $\hat s$ > 0.5, ~0 and <-0.1 in the core region (ρ~ 0.3–0.4) of PS, WS and NCS plasmas, respectively, and flat or negative inside ρ~ 0.4 in the WS and NCS plasmas. Weak magnetic shear is found to be effective in minimizing degradation of ion thermal confinement as T e / T i increases through ECH application, and an improved confinement factor of H 98y2 ~ 1.2 is maintained, similar to NCS plasmas. At T e / T i ~ 1, the ion thermal diffusivity around an internal transport barrier decreases when changing the magnetic shear from positive to weak or negative shear. Also, reduced local particle and momentum transport was indicated by steeper density and toroidal rotation profiles in the weak and negative shear regimes. Linear gyrokinetic simulations predict little change in growth rates of low-k turbulence with ECH application in the WS and NCS plasmas, which is consistent with the transport and profile analyses.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

Wide Operational Windows of Edge-Localized Mode Suppression by Resonant Magnetic Perturbations in the DIII-D Tokamak

Edge-Localized-Mode (ELM) suppression by resonant magnetic perturbations (RMPs) generally occurs over very narrow ranges of the plasma current (or magnetic safety factor $q_{95}$) in the DIII-D tokamak. However, wide $q_{95}$ ranges of ELM suppression are needed for the safety and operational flexibility of ITER and future reactors. In DIII-D ITER Similar Shape (ISS) plasmas with $n$=3 RMPs, the range of $q_{95}$ for ELM suppression is found to increase with decreasing electron density. Nonlinear two-fluid MHD simulations reproduce the observed $q_{95}$ windows of ELM suppression and the dependence on plasma density, based on the conditions for resonant field penetration at the top of the pedestal. When the RMP amplitude is close to the threshold for resonant field penetration, only narrow isolated magnetic islands form near the top of the pedestal, leading to narrow $q_{95}$ windows of ELM suppression. However, as the threshold for field penetration decreases with decreasing density, resonant field penetration can take place over a wider range of $q_{95}$. For sufficiently low density (penetration threshold) multiple magnetic islands form near the top of the pedestal giving rise to continuous $q_{95}$ windows of ELM suppression. The model predicts that wide $q_{95}$ windows of ELM suppression can be achieved at substantially higher pedestal pressure in DIII-D by shifting to higher toroidal mode number ($n$=4) RMPs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

Creation and sustainment of wide pedestal quiescent H-mode with zero net neutral beam torque

Recent experiments on DIII-D have shown it is possible to create and sustain wide pedestal quiescent H-mode (QH-mode) plasmas with zero net torque from neutral beam injection (NBI) for the full discharge duration. Wide pedestal QH-mode has many of the features of the previously investigated QH-mode while having the advantage of increased edge pedestal pressure and excellent energy confinement time. Both QH-mode variants operate without edge localized modes. Accordingly, these discharges demonstrate that significant input torque is not essential to the exploitation of wide pedestal QH-mode in future devices that are expected to have small or nonexistent NBI torque. In conclusion, developing operating conditions that allowed net zero torque access to wide pedestal QH-mode required implementing several techniques to avoid locked modes including minimizing intrinsic error fields, avoiding large sawteeth, and driving toroidal rotation via neoclassical toroidal viscosity.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

Improved core-edge compatibility using impurity seeding in the small angle slot (SAS) divertor at DIII-D

Impurity seeding studies in the slot divertor (SAS) at DIII-D have revealed a strong relationship between the detachment onset and pedestal characteristics with both target geometry and impurity species. N 2 seeding in the slot has led to the first simultaneous observation of detachment on the entire suite of boundary diagnostics viewing the SAS without degradation of core confinement. SOLPS-ITER simulations with D+C+N, full cross-field drifts and n-n collisions activated are performed for the first time in DIII-D to interpret the behavior. This illustrates a strong effect of divertor configuration and plasma drifts on the recycling source distribution with significant consequences on plasma flows. Flow reversal is found for both main ions and impurities affecting strongly the impurity transport and providing an explanation for the observed dependence on strike point location of the detachment onset and impurity leakage found in the experiments. Matched discharges with either nitrogen or neon injection show that while nitrogen does not significantly affect the pedestal, neon leads to increased pedestal pressure gradient and improved pedestal stability. Little nitrogen penetrates in the core, but a great amount of neon is found in the pedestal consistent with the different ionization potentials of the two impurities. This research demonstrates that neutral and impurity distributions in the divertor can be controlled through variations in strike point locations in a fixed baffle structure. Divertor geometry combined with impurity seeding enables mitigated divertor heat flux balancing core contamination and thus leading to enhanced divertor dissipation and improved core-edge compatibility which are essential for ITER and for future fusion reactors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

The role of edge resonant magnetic perturbations in edge-localized-mode suppression and density pump-out in low-collisionality DIII-D plasmas

Two-fluid nonlinear MHD simulations using the TM1 code demonstrate that the formation of magnetic islands at the top and bottom of the H-mode pedestal, together with the strong screening of resonant fields in the gradient region of the pedestal, can account for ELM suppression and density pump-out by n = 2 resonant magnetic perturbations (RMPs) in low-collisionality DIII-D ITER similar shape (ISS) plasmas. Using experimentally relevant transport coefficients, neoclassical resistivity, electron collisionality, and RMP amplitudes, nonlinear MHD simulations reproduce the observed level of density reduction (density pump-out) in DIII-D due to the formation of narrow magnetic islands in the resistive foot of pedestal. For large amplitude RMPs ( B r / B t > 1 × 10 -4 ) simulations predict resonant field penetration and resulting in pressure reduction at the top of the pedestal, consistent with experimental observations at the onset of ELM suppression. The predicted reduction in the height and width of the pedestal by magnetic island enhanced collisional transport provides a quantitative mechanism for the stabilization of the peeling-ballooning mode (PBM). Importantly, these simulations predict strong screening of resonant fields in the steep gradient region of the pedestal due to strong E × B rotation and diamagnetic flows. Yet, if the plasma resistivity is made artificially larger (~ 10X) than neoclassical, the simulations predict magnetic stochasticity throughout the plasma edge and the collapse of the pedestal due to the reduction in the penetration threshold with increasing resistivity. A scaling relation for resonant field penetration at the pedestal top, using several hundred nonlinear simulations, reproduces the density and E × B rotation dependence of the ELM suppression threshold observed in DIII-D. Simulations using ITER model equilibria indicate that the penetration threshold at the top of the ITER pedestal will be lower than in DIII-D due to the anticipated lower perpendicular flow velocities in ITER, resulting in the weaker screening of resonant fields.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

H-mode pedestal improvements with neon injection in DIII-D

Improvements in the H-mode pedestal pressure and global energy confinement were determined on DIII-D with both 3 and 6 MW of neutral beam (NBI) heating with neon injection. The pedestal pressure increased, primarily associated with a significant increase in the pedestal density with a relatively smaller decrease in the pedestal temperature. At the same neon injection rate, a rapid decrease in edge localized mode (ELM) frequency and ~35% increase in plasma stored energy were observed in the 3 MW discharge, while in the 6 MW discharge, the ELM frequency was gradually decreased and the plasma stored energy gradually increased by up to 20%. The measured pedestal widths from both 3 and 6 MW discharges matched the EPED1.0 width scaling, ${\Delta _{{\psi _N}}} = 0.089\beta _{p,ped}^{1/2}$ ,data-, within 20% deviations before and after neon injection. The peeling-ballooning mode (PBM) stability analysis showed that after neon injection, the ballooning stability boundary was increased, while the peeling stability boundary only dropped once the neon level was relatively high. The ballooning stability boundary increase with neon injection is consistent with increase of the ion diamagnetic stabilization (i.e. including the contribution from carbon, neon as well as deuterium ions) which is used in the criterion to determine when the calculated linear growth rates of PBM are unstable. Consequently, the new ELITE computed PBM boundaries were more consistent with the operating points than when only deuterium was included. This in turn suggests that the impurity may help to improve pedestal pressure through affecting ion diamagnetic frequency in the pedestal region, which could positively lead to the ballooning stability boundary extension.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

Enhanced helium exhaust during edge-localized mode suppression by resonant magnetic perturbations at DIII-D

Effective helium confinement time, τ p * ,He , and its ratio with energy confinement time, τ p * ,He / τ E , are key metrics quantifying the suitability of fusion plasmas for a continuous burn. Comparisons in the DIII-D tokamak of discharges with suppression of edge localized modes (ELMs) by resonant magnetic field perturbations (RMPs) to the corresponding unperturbed ELMy discharges found that both of these metrics were strongly reduced, by a factor of approximately 2, after application of the RMPs. This reduction in τ p * ,He during RMP ELM suppression was observed in the plasma core, edge, and pumping plenum, where higher neutral helium concentration during RMPs was also measured. In conclusion, these findings provide evidence that future devices employing RMP ELM suppression may meet or even exceed the helium exhaust provided by the ELMs themselves, reducing helium ash, and thus maintaining high fusion gain.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

On the very high energy confinement observed in super H-mode DIII-D experiments

Analysis of recent super H-mode experiments on DIII-D shows that high rotation, not high pedestal, plays the essential role in achieving very high confinement H 98y2 >1.5. Very high confinement is reached early on in the H-mode phase of these discharges, when the pedestal is still very low, but after the toroidal rotation has already built-up to very high levels in the core. As the discharge evolves, the rotation drops, and so does the energy confinement, despite a sustained very high-pressure pedestal. During this evolution, the confinement quality is linearly correlated with the core toroidal rotation, which varies according to different levels of injected neutral beam torque per particle. Core transport modeling shows that the contribution from rotation in the E×B shear is responsible for confinement quality significantly in excess of standard H-mode (H 98y2 ~1).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

Predicting the rotation profile in ITER

Determining the toroidal rotation for future tokamaks like ITER is a challenging and important problem. By combining empirical scalings for the intrinsic rotation at the top of the pedestal with the expected neutral beam torque and modeling of momentum transport, the toroidal rotation profile for ITER is predicted with TGYRO using TGLF (SAT0 and SAT1). On axis rotation exceeds 20 krad s –1 and the E x B shear is significant enough to reduce turbulent transport and significantly increase confinement and fusion power when comparing to cases that ignore the effect of rotation. Here, the prediction of the rotation at the top of the pedestal is made with increased confidence due to experiments and modeling in DIII-D that have determined the importance of fast-ion and neutral particle transport effects on intrinsic rotation at this location. In particular, the effect of neutral particles on momentum transport in the pedestal region is found to be insignificant.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗