Search NASA⌕ Search

Engineering topics

Burrell, K. H.

Publications and source records attributed to Burrell, K. H..

16 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

The role of ion-scale micro-turbulence in pedestal width of the DIII-D wide-pedestal QH mode

The low-edge rotation, intrinsically ELM-free, and improved confinement wide-pedestal quiescent H-mode (QH-mode), discovered in DIII-D tokamak, has pedestal widths exceeding the EPED-kinetic-ballooning mode (KBM) model scaling typically by at least 25%. Ion-scale ($k_yp_s$ < 1) microturbulence and its role in setting the pedestal structure is investigated using the radially local δ$f$ gyrokinetic code CGYRO. The electromagnetic trapped electron mode (TEM) is unstable at the pedestal top, while plasma beta (β$_e$) is ∼60% below the KBM onset threshold and the electron temperature gradient mode is found to be unstable in the peak gradient region. Nonlinear simulation reveals that the ion-scale turbulence could produce electron energy flux consistent with the flux inferred from power balance at the pedestal top, with a reasonable variation of the local shearing rate; and the local neoclassical transport from NEO is dominant over the simulated turbulent transport in the ion energy flux channel. The simulated ion-scale turbulence produces much lower electron energy flux than inferred from experiment in the pedestal peak gradient region. A correction to the EPED-KBM pedestal width scaling is obtained based on the two-dimensional scan of pedestal top plasma beta (β$_e$) and normalized electron density and temperature scale lengths,${a}$/${L_n}_e$, ${a}$/${L_T}_n$using CGYRO linear simulations. Mode transitions among TEM, micro-tearing mode, ion-temperature gradient mode and KBM, are observed in the 2D scan at the pedestal top. A fixed normalized growth rate for these drift-type modes is taken to determine the pedestal width scaling, which shows good consistency with the QH experimental database on pedestal heights and widths. The onset of KBM instabilities and the local E × B shear suppression criterion set the lower and upper limit for the pedestal width of standard QH-mode, wide-pedestal QH-mode and type-I ELMy H mode. A potentially higher and wider pedestal is expected from the new scaling of pedestal width. This work presents an improved understanding of the ion-scale micro-turbulence of wide-pedestal QH-mode and sheds light on a promising scenario for future reactors, including ITER and beyond.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Comparison of MHD stability properties between QH-mode and ELMy H-mode plasmas by considering plasma rotation and ion diamagnetic drift effects

Magnetohydrodynamic (MHD) stability at tokamak edge pedestal in a quiescent H-mode (QH-mode) and type-I ELMy H-mode plasmas in DIII-D experiment was analyzed by considering plasma rotation and ion diamagnetic drift effects. QH-mode plasma is marginally stable to kink/peeling mode (K/PM), but ELMy H-mode one is almost unstable to peeling-ballooning mode (PBM). It was identified that there are three physics features responsible for the difference in the MHD stability properties between QH-mode plasma and ELMy H-mode one. These are the distance of pedestal foot from the last closed flux surface (LCFS), the amount of the ion diamagnetic drift frequency at pedestal, and impact of coupled rotation and ion diamagnetic drift effects. These features were confirmed through the numerical experiments that the stability properties of the QH-mode plasma can be changed to that of the ELMy H-mode one by shifting the plasma profiles inward in the radial direction and halving the ion diamagnetic drift frequency. The reasons of the change in the stability properties are thought as that K/PM is stabilized due to the inward shift of the bootstrap current profile, and PBM is destabilized due to the reduction of the coupled rotation and ion diamagnetic drift stabilizing effect. Importance of these features was validated through numerical experiments with experimental data of other QH-mode plasmas in DIII-D. All the results show that MHD stability properties of QH-mode plasma can be obtained in case that pedestal foot is close to LCFS, ion diamagnetic drift frequency is large due to high ion temperature, and strong rotation shear exists near pedestal.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Regulation of the central safety factor and normalized beta under low NBI torque in DIII-D

An algorithm has been designed to simultaneously control the central safety factor (q 0 ) and normalized beta (β N ) while ensuring near-zero torque from the neutral beam injection in DIII-D. Feedback control of q 0 and β N in tokamaks can be beneficial due to the close relationship that these variables have with plasma performance and magneto-hydrodynamic stability. In addition, low neutral-beam-torque conditions are of special interest in present devices because future burning-plasma tokamaks such as ITER will most likely operate at very low plasma rotation. The control synthesis of the algorithm presented in this work is based on a linearized, one-dimensional (1D) model of the current-profile dynamics coupled with a zero-dimensional (0D) plasma-energy balance. The actuators considered are neutral beam injection and electron-cyclotron heating and current drive, and discrete logic determines the neutral-beam injection powers that deliver near-zero torque. Here, the algorithm has been tested in nonlinear, 1D simulations using COTSIM (Control-Oriented Transport SIMulator) and in DIII-D experiments, demonstrating satisfactory performance.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
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 OSTI2022

Examination of stiff ion temperature gradient mode physics in simulations of DIII-D H-mode transport

A systematic evaluation of gyrokinetic and gyrofluid model predictions of ion temperature gradient (ITG) stability and transport using parameters from DIII-D high confinement mode (H-mode) plasmas has been performed. The nonlinear CGYRO code is used to make the gyrokinetic predictions, and the quasilinear TGLF model for the corresponding gyrofluid predictions. The assessments are made at three radii (normalized toroidal flux ρtor = 0.4, 0.55, and 0.7) in three different plasma scenarios with varying levels of neutral beam heating and torque. For each of the nine cases (3 radii × 3 scenarios) considered, ITG turbulence is found to be the dominant long-wavelength instability and transport mechanism. The inclusions of both transverse magnetic fluctuations and dynamic fast beam ions are stabilizing for all cases considered, with strongest effects seen at ρor = 0.4 where the fast ion population and normalized plasma pressure β = 2μ0nT/B2 are highest. The further inclusion of parallel magnetic fluctuations does not have a meaningful impact on the ITG turbulence in these scenarios, but does destabilize (in combination with fast ions) new high-frequency instabilities at ρtor = 0.4 in the high power scenarios. In each case the linear and nonlinear ITG critical gradients are predicted to be lower than the measured ITG scale lengths and their associated uncertainties. Inclusion of equilibrium flow shear in the transport predictions generally leads to an upshift in effective critical gradient rather than a qualitative change in the predicted stiffness, with stronger responses typically seen in the gyrokinetic predictions than in the gyrofluid results. However, in most cases these upshifted gradients still remain below the measured values and their uncertainties. Although the predicted critical gradients are below the measured gradients, both models predicted flux-matching gradients consistent with measured values in six of the nine cases considered, with no clear systematic over- or underprediction. Thus, while the experimental ion temperature profiles do not appear to be closely pinned to the ITG critical gradient, both gyrokinetic and gyrofluid models are able to accurately match the measured gradients reasonably well in most cases.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Numerical modeling of pedestal stability and broadband turbulence of wide-pedestal QH-mode plasmas on DIII-D

The wide-pedestal quiescent high confinement mode discovered on DIII-D in recent years is a stationary and quiescent H-mode (QH-mode) with the pedestal width exceeding EPED prediction by at least 25%. Its characteristics, such as low rotation, high energy confinement and edge localized mode-free operation, make it an attractive operation mode for future reactors. Linear and nonlinear simulations using BOUT++ reduced two fluid MHD models and awere carried out to investigate the bursty broadband turbulence often observed in the edge of wide-pedestal QH-mode plasmas. Two kinds of MHD-scale instabilities in different spatial locations within the pedestal were found in the simulations: one mild peeling–ballooning (PB) mode γ PB < 0.04ω A ) located near the minimum in E r</:sub> well propagating in ion diamagnetic drift direction; and one drift-Alfvén wave locates at smaller radius compared to E r</:sub> well propagating in the electron diamagnetic drift direction and unstable only when the parallel electron dynamics is included in the simulation. The coupling between drift wave and shear Alfvén wave provides a possible cause of the experimentally observed local profile flattening in the upper-pedestal. The rotation direction, mode location, as well as the wavenumber of these two modes from BOUT++ simulations agree reasonably well with the experimental measurements, while the lack of quantitative agreement is likely due to the lack of trapped electron physics in current fluid model. This work presents improved physics understanding of the pedestal stability and turbulence dynamics for wide-pedestal QH-mode.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Explaining the lack of power degradation of energy confinement in wide pedestal quiescent H-modes via transport modeling

Wide pedestal quiescent H (WPQH)-mode is an attractive scenario for future burning plasmas as they operate without ELMs. WPQH is characterized by formation of a wider and higher pedestal (than quiescent H-mode), and broadband fluctuations in the pedestal. Unlike conventional H-modes, where the energy confinement time reduces with increasing heating power, the WPQH plasmas reported in this paper do not show power degradation of the energy confinement. As the injected neutral beam power was increased, reduced core (ρ ≤ 0.45) transport calculated by transp, as well as increased core temperatures, pressure gradient and diamagnetic E × B shear rate were observed. The reduction in the heat transport and rapid decrease in the ion temperature gradient scale length suggest the formation of an ion internal transport barrier (ITB) that was accompanied by increased stored energy in the core. Quasilinear turbulent transport modeling using the trapped gyro Landau fluid (tglf) code was used to predict the ITB and its turbulence stability properties. By using profiles and equilibria produced by matching the transp transport fluxes with the tglf transport model within the tgyro transport solver, the energy confinement time captures the experimentally observed insensitivity to the increased PNBI. Linear stability analysis reveals that drift-wave instabilities in the core are stabilized by E × B shear, T i /T e ratio and Shafranov shift; the latter was found to have the strongest effect on the turbulence suppression at the highest heating level.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Stabilization of kink/peeling modes by coupled rotation and ion diamagnetic drift effects in quiescent H-mode plasmas in DIII-D and JT-60U

Magnetohydrodynamic stability at the edge pedestal in several quiescent H-mode (QH-mode) plasmas in DIII-D and JT-60U experiments was analyzed by considering plasma rotation and ion diamagnetic drift effects. It was identified that a kink/peeling mode, which is a prime candidate for a trigger of edge harmonic oscillation in QH-mode, is stabilized by plasma rotation when considering the ion diamagnetic drift simultaneously in both experiments. Furthermore, the stabilizing effect by rotation becomes more effective in case using the rotation profile of the main ion species evaluated by assuming radial force balance. In addition, when inverting the rotation direction, it was found that the kink/peeling mode is more stabilized when considering the rotation of the main ion species, though the mode is less stabilized by the rotation of impurity ion species. The result implies that the kink/peeling mode stability in QH-mode plasmas is sensitive to how the rotation profile is evaluated, but it is shown that a qualitative trend stabilizing the kink/peeling mode by rotation can be reproduced with the poloidal rotation profile of an impurity predicted numerically based on the neoclassical theory.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Turbulence-driven flow dynamics in general axisymmetric toroidal geometry

The present work gives the equations governing the generation of toroidally axisymmetric flows by turbulent Reynolds and Maxwell stresses in finite aspect ratio, general cross section tokamak plasmas. Inclusion of the divergence-free flow constraint in lowest order changes the time scale for evolution of the poloidal flows driven by turbulence by substantial factors. In the pedestal region for present day machines, comparing to earlier cylindrical models, the time scale evaluated using a large aspect ratio circular cross section model can be two orders of magnitude longer while the present, general geometry result can be about one order of magnitude longer. Inclusion of gyroviscosity in the calculation shows that the only lowest order radial velocity fluctuations that enter the problem are those due to fluctuating E×B flows. Toroidal and poloidal flow effects on the toroidally axisymmetric flows are inextricably coupled due to the neoclassical poloidal viscosity. Accordingly, the physics is inherently three dimensional and measurements of all three velocity components are required to obtain the information needed to quantitatively test the theory. The parallel and angular momentum equations for the lowest order, toroidally axisymmetric flows look like radial transport equations when the turbulence is included. The turbulence terms provide the radial transport fluxes. In standard neoclassical theory, the parallel flow equation is local on each flux surface; there is no radial derivative term. However, adding turbulence gives a way, in principle, for radial transport to lead to poloidal flows that deviate from the neoclassical prediction. As a result, the inclusion of the Maxwell stress provides a mechanism for MHD fluctuations to alter the toroidally axisymmetric flows.

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

Impact of shape on pedestal characteristics in the wide pedestal quiescent H-mode in the DIII-D tokamak

Experiments on the DIII-D tokamak have advanced the operational limits of wide pedestal quiescent H-mode (WPQH) plasmas towards increased ITER relevance by simultaneously demonstrating well-matched plasma shape and net zero injected torque. Wide pedestal QHmodes are a compelling candidate regime for a future power producing device because they maintain a stationary pedestal without ELMs via additional edge transport. The pedestal is wider than what would be predicted from kinetic ballooning mode physics due to enhanced edge transport generated by broadband turbulence, a limit cycle oscillation, or some combination thereof. Here, compared to the double null shape, the lower single null shape is observed to have a lower density, narrower pedestal width, larger density fluctuations over a broad range of wavenumber, and operates closer to the peeling-ballooning instability boundary calculated from the simple pedestal scaling, Δ ψΝ , which is still observed to be wider than the EPED prediction.

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↗