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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

TEM turbulence in simulation and experiment in quasisymmetric and degraded-symmetry configurations of the HSX stellarator

The Helically Symmetric Experiment (HSX) has demonstrated reduced neoclassical transport and flow damping with quasisymmetry, but the difference of turbulent transport between quasisymmetric and degraded-symmetry configurations has not been definitively assessed. Profile analysis presented here shows that anomalous transport at the mid-radius is larger with degraded quasisymmetry when temperature and density profiles are matched. An initial comparison (${T_\mathrm{i}/T_\mathrm{e}} = 1$) between the quasisymmetric and degraded symmetry configurations finds that trapped-electron mode (TEM) turbulence can support larger heat fluxes in the degraded symmetry configuration despite smaller linear growth rates. However, further simulations at the experimental temperature ratio ($T_\mathrm{i}/T_\mathrm{e} = 0.2$) do not reproduce this effect. Experimental measurements of heat transport and density fluctuation amplitude are compared to nonlinear gyrokinetic simulations at the experimental temperature ratio. The heat flux and density fluctuation amplitude in simulation reproduce a substantial dependence on the density gradient, and the simulated heat flux matches measurements within experimental uncertainties, indicating that TEM turbulence is responsible for the anomalous transport at the steep density gradient in HSX.

HSX↗
DOE OSTI2025

Separability of microtearing mode and electron temperature gradient turbulence regimes

The separability of microtearing mode (MTM)-dominated and electron temperature gradient (ETG)-driven turbulence regimes is studied with multiscale nonlinear gyrokinetic simulations. The simulations are based on National Spherical Torus Experiment-like, high-confinement mode pedestal parameters, where electromagnetic perturbations are large. Linear analysis indicates a wide scale-separation between the MTM and ETG modes in binormal wavenumber space (perpendicular to the magnetic field line), with no unstable modes at intermediate scales. Likewise, single-scale nonlinear analyses, retaining ion-only or electron-only spatio-temporal scales, produce seemingly well-converged transport states. Surprisingly, the multiscale simulation, which contains both the ion and electron scales, closely follows the transport from the electron-scale simulation. This trend is robust over a wide range of electron temperature gradient. Remarkably, compared to ion-scale simulations, MTM turbulence is significantly reduced at multiscale resolution even when ETG turbulence is low. In this case, traditional ion-scale resolution overestimates the electron energy flux, and it is not possible to accurately simulate the MTM turbulence with separable ion-scale simulations. Here, while the analysis confirms the validity of electron-scale simulations for predicting the electron transport, it also indicates that multiscale simulation may be required for reproducing the turbulence spectrum for systems with coupled MTM-ETG turbulence.

H-mode↗
DOE OSTI2021

The role of ion and electron-scale turbulence in setting heat and particle transport in the DIII-D ITER baseline scenario

In this work, the heat and particle transport in a DIII-D ITER Baseline Scenario (IBS) discharge has been investigated using both linear and nonlinear gyrokinetic simulations performed with the CGYRO code [J. Candy et al: Journal Comp. Phys. (2016)].These simulations were used to investigate the role that ion-scale ($k_θp_s$ < 1:0) and electron-scale ($k_θp_s$ > 1:0) turbulence play in determining heat and particle transport in the core of conditions believed to directly extrapolate to ITER operation. This investigation spans over nearly half of the plasma minor radius, from ρ = 0:45 - 0:85. To probe the nature of the transport and turbulence in these conditions and to validate the model against experimental fluxes, scans of $a/L_{T_i}$, $a/L_{T_e}$, and $a/L_n$ were performed at all radial locations. Long wavelength turbulence is found to be dominated by ITG modes with strong non-adiabatic electron effects and appears unable to reproduce ion and electron heat fluxes and electron particle fluxes simultaneously at most radial positions when single parameter scans are performed. To investigate the potential role of short wavelength turbulence, a series of electron-scale simulations are performed that indicate that experimentally relevant levels of electron heat flux could arise at sub-ion scales. Quantitative comparison of the simulated fluxes obtained from ion and electron-scale simulations with experimental levels along with the sensitivity of simulations results to changes in inputs is presented. Through extensive sensitivity scans and comparison with multiple transport channels, this work demonstrates a clear need for self-consistent modification of multiple inputs and suggests multi-scale interactions play a significant role at many of the radial locations studied. The results of this analysis help shape our understanding of the model fidelity needed to predict turbulence and transport reactor-relevant conditions and have important implications for the prediction of future fusion devices.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

The role of ETG modes in JET–ILW pedestals with varying levels of power and fuelling

Abstract We present the results of GENE gyrokinetic calculations based on a series of JET–ITER-like-wall (ILW) type I ELMy H-mode discharges operating with similar experimental inputs but at different levels of power and gas fuelling. We show that turbulence due to electron-temperature-gradient (ETGs) modes produces a significant amount of heat flux in four JET–ILW discharges, and, when combined with neoclassical simulations, is able to reproduce the experimental heat flux for the two low gas pulses. The simulations plausibly reproduce the high-gas heat fluxes as well, although power balance analysis is complicated by short ELM cycles. By independently varying the normalised temperature gradients ( ω T e ) and normalised density gradients ( ω n e ) around their experimental values, we demonstrate that it is the ratio of these two quantities η e = ω T e / ω n e that determines the location of the peak in the ETG growth rate and heat flux spectra. The heat flux increases rapidly as η e increases above the experimental point, suggesting that ETGs limit the temperature gradient in these pulses. When quantities are normalised using the minor radius, only increases in ω T e produce appreciable increases in the ETG growth rates, as well as the largest increases in turbulent heat flux which follow scalings similar to that of critical balance theory. However, when the heat flux is normalised to the electron gyro-Bohm heat flux using the temperature gradient scale length L T e , it follows a linear trend in correspondence with previous work by different authors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Isotope effects under the influence of global radial electric fields in a helical configuration

Isotope effects under the influence of a radial electric field are examined in a helical magnetic field configuration. We perform global gyrokinetic simulations with additional poloidal rotations to estimate quasi-linear heat flux due to ion temperature gradient mode under the mixing length model. In single-ion-species plasmas, the mass number dependency of heat flux agrees with gyro-Bohm scaling in the absence of a radial electric field. Favorable mass number dependencies violating gyro-Bohm scaling are observed in the presence of a global radial electric field or a heavy hydrogen component in multi-ion-species plasmas. The radial electric field and the heavy hydrogen component affect the heat flux through an increase of wavelength as well as mode stabilization. Poloidal Mach number characterizes the transition from unfavorable to favorable mass number dependency under radial electric fields. While the heat flux is independent of mass number for a given poloidal Mach number, the heat flux decreases for higher mass numbers in a given radial electric field. The heat flux is also independent of average mass number in multi-ion-species plasmas because the heavy hydrogen component effectively enhances the light hydrogen heat flux. The present results are potentially relevant to the violation of gyro-Bohm scaling observed in the recent deuterium experiments in the Large Helical Device.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Nonlinear saturation of reverse shear Alfvén eigenmodes induced by kinetic electrons

A strong nonlinear saturation mechanism of the reverse shear Alfvén eigenmode (RSAE) induced by thermal electrons is observed in gyrokinetic $\delta$$f$ particle-in-cell simulation. This nonlinear effect occurs at moderate mode amplitude, $\delta$$B$ ⟂ /$B$ < 10 -3 , and is associated with the electron parallel streaming along the perturbed magnetic field lines, the term v $\|$ $\frac{δΒ_{⟂}}{B}$ ∙ ∇$\delta$$f$, in the electron drift-kinetic equation. In the case of an n = 4 RSAE, where n is the toroidal mode number, this magnetic fluttering nonlinearity leads to strong mode coupling and broadening of the saturation spectrum. The nonlinear components, which are highly damped through the interaction between thermal electrons and the parallel electric field, modifies the RSAE mode structure and suppresses energetic particle drive significantly. The split-weight scheme (Chen and Parker 2007 J. Comput. Phys. 220 839–55) is used to simulate kinetic electrons.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

First-principles based plasma profile predictions for optimized stellarators

In the present Letter, first-of-its-kind computer simulations predicting plasma profiles for modern optimized stellarators—while self-consistently retaining neoclassical transport, turbulent transport with 3D effects, and external physical sources—are presented. These simulations exploit a newly developed coupling framework involving the global gyrokinetic turbulence code GENE-3D, the neoclassical transport code KNOSOS, and the 1D transport solver TANGO. This framework is used to analyze the recently observed degradation of energy confinement in electron-heated plasmas in the Wendelstein 7-X stellarator, where the central ion temperature was 'clamped' to $T_\textrm i \approx 1.5$ keV regardless of the external heating power. By performing first-principles based simulations, we provide key evidence to understand this effect, namely the inefficient thermal coupling between electrons and ions in a turbulence-dominated regime, which is exacerbated by the large $T_\textrm e/T_\textrm i$ ratios, and show that a more efficient ion heat source, such as direct ion heating, will increase the on-axis ion temperature. This work paves the way towards the use of high-fidelity models for the development of the next generation of stellarators, in which neoclassical and turbulent transport are optimized simultaneously.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Transport barrier and spinning blob dynamics in the tokamak edge

Abstract In this work, we investigate the dynamics of plasma blobs in the edge of magnetic confinement devices using a full-f gyrokinetic particle-in-cell code with X-point geometry. In simulations, the evolution of a seeded blob is followed as it approaches a naturally-forming zonal shear layer near the separatrix, where the blob is stabilized by a large spin induced by the self-consistent adiabatic electron response, and blob bifurcation and trapping are observed during the cross-field propagation of blobs. A new theoretical explanation in both the zonal free and zonal shear layer is constructed, where the dominant E × B spin motion is included. A theoretical condition for a transport barrier induced by the interaction between spinning blobs and the zonal shear layer is obtained, and its scaling is verified with simulations. The new theoretical framework, especially the transport barrier, is applicable to explain and predict various experimental phenomena. In particular, the transport barrier condition calculated with experimental parameters demonstrates that the blob radial transport for H mode is smaller than L mode in experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Reconstruction and interpretation of ionization asymmetry in magnetic confinement via synthetic diagnostics

Abstract Strong poloidal refueling asymmetry in the DIII-D tokamak is inferred from line radiation measurements. Synthetic diagnostics in neutral transport modeling coupled to gyrokinetic simulations illuminate implications for the plasma flow profile in the scrape-off layer of single-null beam-driven discharges. Recycling occurs primarily either on the inner or outer divertor legs, depending on the toroidal magnetic field direction. By reversing the toroidal magnetic field, the observed line radiation asymmetry is nearly eliminated or reversed. It is determined that, while relatively simple physics can describe the observed ionization asymmetry, predicting the overall brightness of the hydrogenic Lyman- α signal requires detailed simulation of the plasma and resulting turbulence. To this end, kinetic plasma simulations fully coupled to comprehensive neutral transport calculations—a novel capability—provide first-principles reproduction of Lyman- α observations on DIII-D.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Assessing the impact of alpha particles on thermal confinement in JET D-T plasmas through global GENE-Tango simulations

The capability of the global, electromagnetic gyrokinetic GENE code interfaced with the transport Tango solver is exploited to address the impact of fusion alpha particles (in their dual role of fast particles and heating source) on plasma profiles and performance at JET in the discharges with the highest quasi-stationary peak fusion power during the DTE2 experimental campaigns. Employing radially global nonlinear electromagnetic GENE-Tango simulations, we compare results with/without alpha particles and alpha heating. Our findings reveal that alpha particles have a negligible impact on turbulent transport, with GENE-Tango converging to similar plasma profiles regardless of their inclusion as a kinetic species in GENE. On the other hand, alpha heating is found to contribute to the peaking of the electron temperature profiles, leading to a 1 keV drop on the on-axis electron temperature when alpha heating is neglected in Tango. The minimal impact of alpha particles on turbulent transport in this JET discharge–despite this being the shot with the highest fusion output–is attributed to the low content of fusion alpha in this discharge. To assess the potential impact of alpha particles on turbulent transport in regimes with higher alpha particle density, as expected in ITER and fusion reactors, we artificially increased the alpha particle concentration to levels expected for ITER. By performing global nonlinear GENE standalone simulations, we found that increasing the alpha particle density beyond five times the nominal value lead to significant overall turbulence destabilization. These results demonstrate that an increased alpha particle concentration can significantly impact transport properties under simulated JET experimental conditions. However, these findings cannot be directly extrapolated to ITER due to the substantial differences in parameters such as plasma size, magnetic field, plasma current, and thermal pressure.

energetic particles↗
DOE OSTI2026

Turbulence-Driven Edge-Localized-Mode-Free High-Confinement Mode with Divertor Detachment in a Metal-Wall Tokamak

We report the first demonstration of a minute-scale, edge-localized-mode-free high-confinement plasma regime compatible with divertor partial detachment and enhanced pedestal performance in a metal-wall tokamak, the Experimental Advanced Superconducting Tokamak. This regime is enabled by a newly identified mechanism: during divertor partial detachment, reduced ionization and enhanced pumping in a closed divertor lead to less cooling of the pedestal by recycling neutrals and seeding impurities, resulting in an increased pedestal temperature gradient, which excites high-frequency broadband turbulence. Gyrokinetic simulations identify the high-frequency broadband turbulence as a temperature-gradient-driven trapped electron mode (𝜂 𝑒 -TEM), which drives outward transport of particles and heat, thereby maintaining the edge-localized-mode-free state. This pedestal regime is particularly promising for the International Thermonuclear Experimental Reactor, where the anticipated lower density gradient, reduced E×B shear, and lower collisionality in the pedestal will further facilitate 𝜂 𝑒 -TEM excitation. The achieved integrated scenario with a detached divertor and turbulence-dominated pedestal thus offers a compelling solution for managing heat loads and metal impurity sources for long-pulse high-performance operation in future fusion reactors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Direct microstability optimization of stellarator devices

Turbulent transport is regarded as one of the key issues in magnetic confinement nuclear fusion, both for tokamaks and stellarators. Here, in this work, we show that a significant decrease in a microstability-based proxy, as opposed to a geometric one, for the turbulent heat flux, namely the quasilinear heat flux, can be obtained in an efficient manner by coupling stellarator optimization with linear gyrokinetic simulations. This is accomplished by computing the quasilinear heat flux at each step of the optimization process, as well as the deviation from quasisymmetry, and minimizing their sum, leading to a balance between neoclassical and the turbulent transport proxy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Impact of Magnetic Field Configuration on Heat Transport in Stellarators and Heliotrons

We assess the magnetic field configuration in modern fusion devices by comparing experiments with the same heating power, between a stellarator and a heliotron. The key role of turbulence is evident in the optimized stellarator, while neoclassical processes largely determine the transport in the heliotron device. Gyrokinetic simulations elucidate the underlying mechanisms promoting stronger ion scale turbulence in the stellarator. Similar plasma performances in these experiments suggests that neoclassical and turbulent transport should both be optimized in next step reactor designs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Partnership Center for High-fidelity Boundary Plasma Simulation

The University of Colorado Boulder contributed to the DOE Partnership Center for High‑fidelity Boundary Plasma Simulation (award DE‑SC0018271) by advancing gyrokinetic modeling of core–edge plasma interactions in tokamaks. Using the GEM and XGC codes, the project investigated turbulence spreading from the plasma edge to the core, finding that while edge‑driven trapped electron mode turbulence enhances core ion heat flux in DIII‑D L‑mode plasmas, it does not fully explain experimentally observed transport shortfalls, indicating a need for flux‑driven models. The team also demonstrated that edge micro‑instabilities are highly sensitive to small variations in magnetic equilibrium, particularly near the pedestal and edge, highlighting the importance of accurately characterizing equilibrium uncertainties. Significant progress was made in spatially coupling the core GEM and edge XGC simulations, with successful validation against standalone XGC results in realistic tokamak geometries, laying the groundwork for future extensions to kinetic electrons and electromagnetic effects.

42 ENGINEERING↗
DOE OSTI2024

The Information Length Concept Applied to Plasma Turbulence

A methodology to study statistical properties of anomalous transport in fusion plasma is investigated. Three time traces generated by the full-f gyrokinetic code GKNET are analyzed for this purpose. The time traces consist of heat flux as a function of the radial position, which is studied in a novel manner using statistical methods. The simulation data exhibit transport processes with both medium and long correlation length along the radius. A typical example of a phenomenon with long correlation length is avalanches. In order to investigate the evolution of the turbulent state, two basic configurations are studied, one flux-driven and one gradient-driven with decaying turbulence. The information length concept in tandem with Boltzmann–Gibbs and Tsallis entropy is used in the investigation. It is found that the dynamical states in both flux-driven and gradient-driven cases are surprisingly similar, but the Tsallis entropy reveals differences between them. This indicates that the types of probability distribution function are nevertheless quite different since the higher moments are significantly different.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Data for "Verification of the kinetic electron role in the microinstabilities in a negative triangularity model equilibrium"

This contains the dataset used in "Verification of the kinetic electron role in the microinstabilities in a negative triangularity model equilibrium" published in Physics of Plasmas in Oct 2024. It consists of plaintext files as well as .bp files (which may be read by ADIOS open-source software) used for creating the plots that appear in the paper. These are derived from simulation outputs from XGC (X-point Gyrokinetic Code). The data mainly consists of growth rate and frequency measurements as well as poloidal cross-section data.

gyrokinetic↗
DOE OSTI2026

Modeling Temperature Profiles in the Pedestal of NSTX with Reduced Models

This paper describes new modeling capabilities for predicting H-mode pedestal profiles in spherical tokamaks. Temperature profiles for NSTX discharges 132543 and 132588 are modeled by coupling the \textsc{astra} transport solver with neoclassical transport and gyrokinetic-based reduced models for electron temperature gradient (ETG) and kinetic ballooning mode (KBM) instabilities. A quasi-linear surrogate model for ion-scale transport is developed using linear \textsc{gene} simulations, requiring only a single free parameter calibrated to one discharge. Time-evolving the temperatures with fixed density yields good agreement with experiments for both discharges. Systematic analysis of the transport mechanisms reveals that neoclassical transport is huge across the entire pedestal region for the ion channel. ETG turbulence is large in the plasma edge and low density gradient region, contributing substantially to the electron channel. However, KBM/MHD-like modes also drive significant transport in both the ion and electron thermal channels, making them essential for accurate pedestal modeling. Further refinements, including explicit E×B shear suppression and scaled ETG transport, produce quantitative but not qualitative improvements. This work lays the foundation for predictive modeling of future devices. This paper is on arxiv and has been submitted to Nuclear Fusion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Effect of Magnetic Islands on Neoclassical Heat Diffusivity in a Global Tokamak Simulation

Resonant magnetic perturbations (RMPs) alter magnetic field topology via island formation and can modify core transport in toroidal fusion devices. Using the global gyrokinetic particle‐in‐cell code XGC‐S, originally developed for stellarator geometries, we quantify how island topology affects neoclassical heat transport in a circular tokamak with ($m,n$) = (5,2) RMP‐induced islands. Electron radial heat diffusivity exhibits a two‐peak structure with a dominant enhancement near the O/X‐point region and a secondary peak at the outer island boundary, whereas ion transport remains close to the neoclassical prediction. The observed electron diffusivity is strongly enhanced over the neoclassical baseline but remains well below the Rechester–Rosenbluth estimate, consistent with the non‐ergodic character of the island topology. The density dependence of diffusivity is weaker inside the islands than outside, indicating the presence of topology‐driven transport channels that do not depend strongly on collisionality. Two‐dimensional heat flux maps reveal an up–down antisymmetric structure at the island separatrices consistent with magnetic gradient and curvature drift effects, in contrast to the symmetric temperature flattening inside the islands. These results identify topology‐dependent avenues for electron heat transport in RMP‐perturbed cores and suggest regimes in which island geometry, rather than collisions alone, governs radial transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗