Search NASA⌕ Search

SEARCH · Search NASA

Results for “Gyrokinetics”

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 307 records · Page 17

Theory of zonal flow growth and propagation in toroidal geometry

The toroidal geometry of tokamaks and stellarators is known to play a crucial role in the linear physics of zonal flows (ZFs), leading to e.g. the Rosenbluth–Hinton residual and geodesic acoustic modes. However, descriptions of the nonlinear ZF growth from a turbulent background typically resort to simplified models of the geometry. We present a generalised theory of the secondary instability to model the ZF growth from turbulent fluctuations in toroidal geometry, demonstrating that the radial magnetic drift substantially affects the nonlinear ZF dynamics. In particular, the toroidicity gives rise to a new branch of propagating ZFs, the toroidal secondary mode, which is nonlinearly supported by the turbulence. We present a theory of this mode and compare the theory against gyrokinetic simulations of the secondary mode. The connection with other secondary modes—the ion-temperature-gradient and Rogers–Dorland–Kotschenreuther secondary modes—is also examined.

Nies, Richard [Princeton Univ., NJ (United States)↗

Investigation of mode activity in NBI-heated experiments of Wendelstein 7-X

The 2018 operation phase (OP 1.2b) of the stellarator Wendelstein 7-X (W7-X) included, for the first time, neutral beam injection (NBI) to heat the plasma. Since the injection geometry at W7-X is not parallel, this generates both passing and trapped fast particles. During longer phases of NBI injection, with the primary purpose to study the heating efficiency of this system, Alfvén eigenmodes (AEs) were observed by a number of diagnostics, including the phase contrast imaging (PCI) system, the magnetic pick-up coils (Mirnov coils), and the soft x-ray multi-camera tomography system (XMCTS). Alfvén eigenmodes are of great interest for future fusion reactors as it has been shown that the resonant interaction of fast ions with self-excited AEs can lead to enhanced transport of fast ions and potentially to energy losses. This is especially true for so-called gap-modes, Alfvén eigenmodes with frequencies in gaps of the continuous spectrum, since they lack continuum damping. These modes are commonly known to be excited by fast ions, but other destabilizing mechanisms, e.g. the electron-pressure gradient are also possible. In this article we present a first analysis of the experimentally observed frequencies from the theoretical side. The calculation of shear Alfvén wave continua for selected cases and the assignment of observed frequencies to the gaps of the continuous spectra are presented. Using the ideal-MHD code CKA (Könies A. 2007 10th IAEA TM on Energetic Particles in Magnetic Confinement System), we find gap modes that match the experimental measurements in terms of the observed frequencies. We emphasize the crucial roles played by the coupling of sound and Alfvén waves as well as of the Doppler shift arising as a consequence of the radial electric field in W7-X. We employ the perturbative gyrokinetic code CKA-EUTERPE (Feh´er 2013 Simulation of the interaction between Alfv´en waves and fast particles), using a slowing-down distribution function for the fast ions as calculated by the Monte-Carlo particle following code ASCOT (Hirvijoki et al 2014 Comput. Phys. Commun. 185 1310–21) to assess the fast-ion drive. Here we find that the fast-ion drive is insufficient to overcome the background-plasma damping. The fact that unstable modes were observed experimentally may point to problems with the modelling or indicate the existence of other destabilizing mechanisms, e.g. associated with the electron-pressure gradient (Windisch et al 2017 Plasma Phys. Control. Fusion 59 105002) that sensitively depend on the profiles of the background plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Impact of resonant magnetic perturbations on zonal flows and microturbulence

Results from a joint experimental and computational effort studying the effect of resonant magnetic perturbations (RMPs) on microturbulence levels and their connection to zonal flows in the DIII-D tokamak L-mode are presented. Beam emission spectroscopy measurements show a direct increase in density fluctuations at microturbulent scales with increasing RMP amplitude, suggesting that magnetic activity introduced by the RMP affects the regulation of microturbulence on DIII-D. This is analogous to how MHD-scale magnetic fluctuations arising from tearing modes have been observed in simulations to increase microturbulence levels in the reversed-field pinch (RFP). In the RFP, this is attributed to magnetic fluctuations eroding turbulence-limiting zonal flows; this work examines if a similar mechanism is present for DIII-D microturbulence. Gyrokinetic simulations find that the application of an RMP corresponds directly to a decrease in zonal flow levels, producing a similar increase of turbulent fluctuation levels over a range of RMP amplitudes as observed in the experiment.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Inference of experimental radial impurity transport on Alcator C-Mod: Bayesian parameter estimation and model selection

We present a fully Bayesian approach for the inference of radial profiles of impurity transport coefficients and compare its results to neoclassical, gyrofluid and gyrokinetic modeling. Using nested sampling, the Bayesian Impurity Transport InferencE (BITE) framework can handle complex parameter spaces with multiple possible solutions, offering great advantages in interpretative power and reliability with respect to previously demonstrated methods. BITE employs a forward model based on the pySTRAHL package, built on the success of the well-known STRAHL code [Dux, IPP Report, 2004], to simulate impurity transport in magnetically-confined plasmas. In this paper, we focus on calcium (Ca, Z=20) Laser Blow-Off injections into Alcator C-Mod plasmas. Multiple Ca atomic lines are diagnosed via high-resolution X-ray Imaging Crystal Spectroscopy and Vacuum Ultra-Violet measurements. We analyze a sawtoothing I-mode discharge for which neoclassical and turbulent (quasilinear and nonlinear) predictions are also obtained. Here, we find good agreement in diffusion across the entire radial extent, while turbulent convection and density profile peaking are estimated to be larger in experiment than suggested by theory. Efforts and challenges associated with the inference of experimental pedestal impurity transport are discussed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

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↗

Radial electric field and density fluctuations measured by Doppler reflectometry during the post-pellet enhanced confinement phase in W7-X

Radial profiles of density fluctuations and the radial electric field, Er, have been measured using Doppler reflectometry during the post-pellet enhanced confinement phase achieved, under different heating power levels and magnetic configurations, during the 2018 W7-X experimental campaign. A pronounced Er-well is measured with local values as high as -40 kV m -1 in the radial range ρ ~ 0.7–0.8 during the post-pellet enhanced confinement phase. The maximum Er intensity scales with both the plasma density and electron cyclotron heating power level, following a similar trend to the plasma energy content. A good agreement is found when the experimental Er profiles are compared to simulations carried out using the neoclassical codes, the drift kinetic equation solver (DKES) and kinetic orbit-averaging solver for stellarators (KNOSOS). The density fluctuation level decreases from the plasma edge toward the plasma core and the drop is more pronounced in the post-pellet enhanced confinement phase than in reference gas-fuelled plasmas. Besides, in the post-pellet phase, the density fluctuation level is lower in the high iota magnetic configuration than in the standard one. Therefore, to determine whether this difference is related to the differences in the plasma profiles or to the stability properties of the two configurations, gyrokinetic simulations have been carried out using the codes stella and EUTERPE. The simulation results point to the plasma profile evolution after the pellet injection and the stabilization effect of the radial electric field profile as the dominant players in the stabilization of the plasma turbulence.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Alpha particle driven Alfvénic instabilities in ITER post-disruption plasmas

Fusion-born alpha particles in ITER disruption simulations are investigated as a possible drive of Alfvénic instabilities. The ability of these waves to expel runaway electron (RE) seed particles is explored in the pursuit of a passive, inherent RE mitigation scenario. The spatiotemporal evolution of the alpha particle distribution during the disruption is calculated using the linearized Fokker–Planck solver CODION coupled to a fluid disruption simulation. These simulations are done in the limit of no alpha particle transport during the thermal quench, which can be seen as a most pessimistic situation where there is also no RE seed transport. Under these assumptions, the radial anisotropy of the resulting alpha population provides free energy to drive Alfvénic modes during the quench phase of the disruption. We use the linear gyrokinetic magnetohydrodynamic code LIGKA to calculate the Alfvén spectrum and find that the equilibrium is capable of sustaining a wide range of modes. The self-consistent evolution of the mode amplitudes and the alpha distribution is calculated utilizing the wave-particle interaction tool HAGIS. Intermediate mode number (n = 7–15, 22–26) toroidal Alfvén eigenmodes are shown to saturate at an amplitude of up to δB/B ≈ 0.1% in the spatial regimes crucial for RE seed formation. We find that the mode amplitudes are predicted to be sufficiently large to permit the possibility of significant radial transport of REs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Effects of zonal flows on ion temperature gradient instability in the scrape-off layer of a field-reversed configuration

Gyrokinetic simulations of long wavelength ion temperature gradient (ITG) turbulence in the scrape-off layer (SOL) of a field-reversed configuration (FRC) find that zonal flows are nonlinearly generated and are the dominant mechanism for the nonlinear saturation of the ITG instability. After the ITG saturation, zonal flows remain undamped and gradually suppress the turbulent transport to a very low level. In the simulations with collisions, collisional damping gradually reduces zonal flow amplitude to a lower level, which allows finite ITG turbulence intensity and ion heat transport in the SOL. The steady state turbulence intensity and ion heat transport are found to be proportional to the collision frequency. Furthermore, this favorable scaling suggests that minimizing collisions (e.g. increasing temperature, reducing impurity content, etc) and preserving toroidal symmetry could improve plasma confinement in the FRC.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Experimental inference of neutral and impurity transport in Alcator C-Mod using high-resolution x-ray and ultra-violet spectra

We present experimental inferences of cross-field impurity transport coefficients for AlcatorC-Mod plasmas without edge-localized modes, leveraging a novel forward model for the entire Ca Kα spectrum, including satellite lines within the spectral range, to compare to high-resolution x-ray imaging crystal spectroscopy (XICS). These measurements are complemented by extreme ultra-violet (EUV) spectroscopy that constrains transport closer to the edge. Using new atomic data sets for both XICS and EUV analysis has enabled consideration of line ratios across both spectral ranges and has increased the accuracy of inferred transport coefficients. Inclusion of charge exchange between edge thermal neutrals and impurities is shown to be extremely important in C-Mod pedestals. Atomic D neutral densities from experimental D Lyα measurements at the midplane are compared to SOLPS-ITER simulations, finding good agreement. Bayesian inferences of impurity transport coefficients are presented for L-, EDA H-, and I-mode discharges, making use of the Aurora package for forward modeling and combining our spectroscopic constraints. Experimentally inferred diffusion profiles are found to match turbulent transport models at midradius within uncertainties, using both quasilinear gyro-fluid TGLF SAT-1 and nonlinear ion-scale gyrokinetic CGYRO simulations. Furthermore, significant discrepancies in convection are observed in some cases, suggesting difficulties in predictions of flat or hollow impurity profiles.

neutrals↗

Generation of E × B flow shear by finite orbit width effects from heat sources in tokamaks

Here we present a possible mechanism for the generation of strong E × B flow shear relevant to internal transport barrier formation in tokamak plasmas. From gyrokinetic calculations, we show that strong E × B flow shear can be generated by finite orbit width (FOW) effects associated with a non-uniform heat source and is sufficient to lead to transport barrier formation in the core region with a moderate power level. Two FOW effects inducing neoclassical polarization are shown to be responsible for this: (1) the radial drift of particle orbit center due to the variation of the heat source within orbit width and (2) the non-uniformly evolved orbit width by the non-uniform heating.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Deconvolving the roles of E × B shear and pedestal structure in the energy confinement quality of super H-mode experiements

Analysis of “super H-mode” experiments on DIII-D has put forward that high plasma toroidal rotation, not high pedestal, plays the essential role in achieving energy confinement quality H 98y2 >>1. Recently, super H-mode experiments with variable input torque have confirmed that high rotation shear discharges have very high levels of H 98y2 (>1.5), independent of the pedestal height, and that high pedestal discharges with low rotation shear have levels of H 98y2 only slightly above 1 (≤1.2). Although some increase in stored energy with higher pedestal occurs, the energy confinement quality mainly depends on the toroidal rotation shear, which varies according to different levels of injected neutral beam torque per particle. Here, quasi-linear gyrofluid modeling achieves a good match of the experiment when including the E×B shear; without including plasma rotation, the modeling predicts a confinement quality consistent with the empirical observation of H 98y2 ~1.2 at low rotation. Nonlinear gyrokinetic transport modeling shows that the effect of E×B turbulence stabilization is far larger than other mechanisms, such as the so-called hot-ion stabilization (T i /T e ) effect. Consistent with these experimental and modeling results are previous simulations of the ITER Baseline Scenario using a super H-mode pedestal solution, which showed the potential to exceed the Q=10 target if the pedestal density could be increased above the Greenwald limit. A close look at these simulations reveals that the predicted energy confinement quality is below 1 even at the highest pedestal pressure. The improvement in Q at higher pedestal density is due to the improved fusion power generation at the higher core density associated with higher pedestal density, not to an improved energy confinement quality.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Three-dimensional inhomogeneity of electron-temperature-gradient turbulence in the edge of tokamak plasmas

Nonlinear multiscale gyrokinetic simulations of a Joint European Torus edge pedestal are used to show that electron-temperature-gradient (ETG) turbulence has a rich three-dimensional structure, varying strongly according to the local magnetic-field configuration. In the plane normal to the magnetic field, the steep pedestal electron temperature gradient gives rise to anisotropic turbulence with a radial (normal) wavelength much shorter than in the binormal direction. In the parallel direction, the location and parallel extent of the turbulence are determined by the variation in the magnetic drifts and finite-Larmor-radius (FLR) effects. The magnetic drift and FLR topographies have a perpendicular-wavelength dependence, which permits turbulence intensity maxima near the flux-surface top and bottom at longer binormal scales, but constrains turbulence to the outboard midplane at shorter electron-gyroradius binormal scales. Here our simulations show that long-wavelength ETG turbulence does not transport heat efficiently, and significantly decreases overall ETG transport—in our case by ~40%—through multiscale interactions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Zonal flow screening in negative triangularity tokamaks

Here, this paper presents a comparative study of zonal flow screening in positive and negative triangularity tokamaks. Neoclassical screening sets the strength of zonal flow shear for fixed drive and damping. Orbit calculations show that the radial excursions of trapped particles are reduced in negative triangularity tokamaks, as compared to positive triangularity. Yet surprisingly, the neoclassical dielectric susceptibility actually increases with decreasing triangularity, such that the susceptibility is higher for negative triangularity than for positive triangularity. This is because the reduction in neoclassical polarization by shrinking the banana width is offset by the increase in neoclassical polarization by the enhancement of trapped fraction for negative triangularity. As a result, the zonal flow screening length is actually enhanced for negative triangularity, as compared to positive triangularity. Hence, the zonal flow residual is smaller for negative triangularity than for positive triangularity. Results from gyrokinetic simulations support the analytic calculations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

On the role of mode resonances in regulating zonal-flow-moderated plasma microturbulence

Abstract The onset of turbulent heat transport at a higher temperature gradient than the critical gradient of linear instability, known as the Dimits shift, is a recurring feature of nonlinear simulations for magnetically confined fusion plasmas. Resonance in the nonlinear coupling between the modes that dominate energy transfer can lead to suppression of turbulence and transport above the linear critical gradient. As an expression of this resonance, gyrokinetic simulations show a quasi-coherent interaction between streamers and sidebands coupled through the zonal flow within the Dimits regime. This mechanism is further confirmed by use of artificial complex frequencies which break the resonance. By incorporating corresponding saturation physics, the standard quasilinear model for rapid head flux prediction is improved, which can now predict reduced heat flux in the Dimits regime. In particular, the triplet correlation time, the lifetime of the nonlinear interaction, is shown to be well-approximated by combinations of linear eigenvalues, and yields good representations of the heat flux variation both in and above the Dimits regime. Thus, a reduced but predictive model for transport near the critical gradient of zonal-flow saturated turbulence now exists.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Electrostatic turbulence in EAST plasmas with internal transport barrier

Based on first-principles nonlinear gyrokinetic simulations, the electrostatic turbulence properties in the internal transport barrier (ITB) region of an Experimental Advanced Superconducting Tokamak discharge (#93890) are investigated. Specifically, ITBs with steep density and temperature gradients are located in the weakly negative magnetic shear region at the plasma center. In the linear stage, the growth rate and frequency of the ion temperature gradient (ITG) mode increase significantly due to resonant excitation by trapped electrons. That is, the resonance between trapped electrons and the ITG becomes strong due to the precession drift reversal of trapped electrons by the negative magnetic shear and Shafranov shift. Meanwhile, the trapped electron mode is stable in the ITB region due to only a very small fraction of electrons precessing in the direction of the electron diamagnetic drift. Nonlinear simulations show that, after considering the non-adiabatic effect of trapped electrons, the heat conductivity of ions and the turbulence intensity increase by at least a factor of 7 compared with the results only considering the adiabatic effect of electrons. The zonal charge density of trapped electrons can partially cancel that of ions, which weakens the intensity of the zonal flow, and consequently reduces the zonal flow regulation and enhances the turbulent transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Micro-tearing mode dominated electron heat transport in DIII-D H-mode pedestal

Abstract A new, comprehensive set of evidence reveals that Micro-Tearing Modes (MTMs) dominate pedestal electron heat transport in an H-mode experiment in the DIII-D tokamak. The experiment investigates the role of MTMs by scanning pedestal collisionality, a main drive of MTM instability, from 0.43 to 0.84 on the pedestal top. Broadband (150–800 kHz) magnetic and density fluctuations originating from the pedestal gradient region and highly consistent with MTMs are observed, with amplitude increasing during the scan. The higher magnetic fluctuation amplitude correlates with a lower pedestal electron temperature gradient, implying MTMs may regulate the pedestal electron heat transport. The collisionality scan results in profile and transport changes consistent with predicted transport capability of MTMs: (1) experimentally-determined electron heat diffusivity increases ∼40% at the location where the broadband density fluctuations peak; (2) ion heat diffusivity has less increase (<20%); and (3) a locally flattened region in the electron temperature pedestal is observed at high collisionality. A local, linear gyrokinetic simulation finds MTMs as the most unstable mode in the pedestal gradient region. In addition, local, nonlinear simulations suggest MTMs can dominate and drive experimentally-relevant, megawatt-level electron heat flux. This result establishes MTMs as an effective transport mechanism in the H-mode pedestal, in particular at high collisionality.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Direct measurement of the electron turbulence-broadening edge transport barrier to facilitate core–edge integration in tokamak fusion plasmas

Abstract The integration of a high-performance core and a dissipative divertor, or the so-called ‘core–edge integration,’ has been widely identified as a critical gap in the design of future fusion reactors. In this letter, we report, for the first time, direct experimental evidence of electron turbulence at the DIII-D H-mode pedestal that correlates with the broadening of the pedestal and thus facilitates core–edge integration. In agreement with gyrokinetic simulations, this electron turbulence is enhanced by high η e ( η e = L n / L T e , where L n is the density scale length and L T e is the electron temperature scale length), which is due to a strong shift between the density and temperature pedestal profiles associated with a closed divertor. The modeled turbulence drives significant heat transport with a lower pressure gradient that may broaden the pedestal to a greater degree than the empirical and theoretically predicted pedestal width scalings. Such a wide pedestal, coupled with a closed divertor, enables us to achieve a good core–edge scenario that integrates a high-temperature low-collisionality pedestal (pedestal top temperature T e,ped > 0.8 keV and a pedestal top collisionality ν * ped < 1) under detached divertor conditions. This paves a new path toward solving the core–edge integration issue in future fusion reactors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗