Search NASA⌕ Search

SEARCH · Search NASA

Results for “gyrofluid & 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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Comparing particle transport in JET and DIII-D plasmas: gyrokinetic and gyrofluid modelling

The paper models the impact of fueling versus particle transport in the DIII-D and JET tokamaks. Experimentally, it was found that in JET an increase in the density peaking at lower collisionality could be attributed to an increase in core fueling, however on DIII-D, the opposite was observed. In DIII-D, the increase in density peaking was the result of an increase in inward transport, not an increase in core fueling. The paper shows that the difference in temperature gradients and thus turbulence dynamics between the two experiments can explain the experimentally observed discrepancy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Radiative and continuum dampings of reversed shear Alfvén eigenmodes and perturbative analysis limitations for tokamaks

A careful theoretical analysis of the excitation of Alfvén eigenmodes (AEs), such as TAE (toroidicity-induced AE) and RSAE (reversed shear AE), by superalfvenic energetic particles is required for reliable predictions of energetic ion relaxation in present day fusion experiments. This includes the evaluation of different AE damping mechanisms including radiative and continuum dampings which are the focus of this study. A recent comprehensive benchmark of different eigenmode solvers including gyrokinetic, gyrofluid and hybrid magenetohydrodynamics (MHD) has shown that employed models may have deficiencies when addressing some of them (Taimourzadeh et al., Nucl. Fusion, vol. 59, 2019, 066006). Here, in this paper, we are studying the radiative and continuum dampings of RSAEs in details which were missing in hybrid NOVA/NOVA-C calculations to prepare a NOVA-C package with a substantial upgrade. Both dampings require the finite Larmor radius (FLR) corrections to AE mode structures to be accounted for. Accurately calculating different damping rates and understanding their parametric dependencies, we resolve the limitation coming out of the perturbative approach. In particular, here, the radiative damping is included perturbatively, whereas the continuum damping is computed non-perturbatively. Our comparison leads to the conclusion that the non-perturbative treatment of the unstable RSAE modes is needed to find the agreement with the gyrokinetic calculations. We expect that the RSAE mode structure modification plays a dominant role in determining the RSAE stability.

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 OSTI2025

Paradigm for global gyrokinetic turbulence

Using a novel wave-number-advection algorithm, we describe a method to add global profile curvature terms to the local gyrokinetic equations. This approach enables a high-precision global simulation capability without sacrificing the efficiency and spectral accuracy of local simulations. Preliminary numerical results suggest that there are two distinct types of global correction to the local result: (1) a trivial shift of the local result in radius, (2) a non-trivial decorrelation that is roughly equivalent to the well-known E×B shear effect. This interpretation provides a new, intuitive framework for understanding global gyrokinetic simulations.

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 OSTI2024

Saturation of Fishbone Instability by Self-Generated Zonal Flows in Tokamak Plasmas

Gyrokinetic simulations of the fishbone instability in DIII-D tokamak plasmas find that self-generated zonal flows can dominate the nonlinear saturation by preventing coherent structures from persisting or drifting in the energetic particle phase space when the mode frequency down-chirps. Results from the simulation with zonal flows agree quantitatively, for the first time, with experimental measurements of the fishbone saturation amplitude and energetic particle transport. Moreover, the fishbone-induced zonal flows are likely responsible for the formation of an internal transport barrier that was observed after fishbone bursts in this DIII-D experiment. Finally, gyrokinetic simulations of a related ITER baseline scenario show that the fishbone induces insignificant energetic particle redistribution and may enable high performance scenarios in ITER burning plasma experiments. Published by the American Physical Society 2024

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 OSTI2024

Intrinsic Toroidal Rotation Driven by Turbulent and Neoclassical Processes in Tokamak Plasmas from Global Gyrokinetic Simulations

Gyrokinetic tokamak plasmas can exhibit intrinsic toroidal rotation driven by the residual stress. While most studies have attributed the residual stress to the parallel-momentum flux from the turbulent E × B motion, the parallel-momentum flux from the drift-orbit motion (denoted $Π^D_\parallel$) and the E × B-momentum flux from the E × B motion (denoted $Π_{E×B}$) are often neglected. Here, we use the global total-f gyrokinetic code XGC to study the residual stress in the core and the edge of a DIII-D H-mode plasma. Numerical results show that both $Π^D_\parallel$ and $Π_{E×B}$ make up a significant portion of the residual stress. In particular, $Π^D_\parallel$ in the core is higher than the collisional neoclassical level in the presence of turbulence, while in the edge it represents an outflux of countercurrent momentum even without turbulence. Using a recently developed “orbit-flux” formulation, we show that the higher-than-neoclassical-level $Π^D_\parallel$ in the core is driven by turbulence, while the outflux of countercurrent momentum from the edge is mainly due to collisional ion orbit loss. In conclusion, these results suggest that $Π^D_\parallel$ and $Π_{E×B}$ can be important for the study of intrinsic toroidal rotation.

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

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↗
DOE OSTI2021

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

This paper presents comparisons of predictions from two models of turbulent transport in tokamaks. The focus is on assessing how close a reduced but significantly less computationally demanding model (TGLF) does in matching the more detailed model (CGYRO) in predicting scalings with driving gradient ("stiffness"). The impacts of other physics mechanisms such as flow shear, fast ions, and electromagnetic fluctuations on these predictions is also studied. Overall, while there are clear differences in the scaling, these do not appear to significantly impact the fidelity of the reduced model in prediction the experimental conditions studied relative to the more advanced model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

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. 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.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Investigation of core impurity transport in DIII-D diverted negative triangularity plasmas

Abstract Tokamak operation at negative triangularity has been shown to offer high energy confinement without the typical disadvantages of edge pedestals (Marinoni et al 2021 Nucl. Fusion 61 116010). In this paper, we examine impurity transport in DIII-D diverted negative triangularity experiments. Analysis of charge exchange recombination spectroscopy reveals flat or hollow carbon density profiles in the core, and impurity confinement times consistently shorter than energy confinement times. Bayesian inferences of impurity transport coefficients based on laser blow-off injections and forward modeling via the Aurora package (Sciortino et al 2021 Plasma Phys. Control. Fusion 63 112001) show core cross-field diffusion to be higher in L-mode than in H-mode. Impurity profile shapes remain flat or hollow in all cases. Inferred radial profiles of diffusion and convection are compared to neoclassical, quasilinear gyrofluid, and nonlinear gyrokinetic simulations. Heat transport is observed to be better captured by reduced turbulence models with respect to particle transport. State-of-the-art gyrokinetic modeling compares favorably with measurements across multiple transport channels. Overall, these results suggest that diverted negative triangularity discharges may offer a path to a highly-radiative L-mode scenario with high core performance.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

A tale of two polarization paradoxes: The diamagnetic polarization paradox

An accurate calculation of the total polarization charge density in a plasma is essential for a self-consistent determination of the electric field. Yet, the “diamagnetic polarization paradox” refers to the fact that there is a paradoxical factor of 1/2 difference between the pressure-driven “diamagnetic polarization” density calculated using real space drift theory vs action-angle space guiding center and gyrokinetic theory that has not been explained before. Here, we show that both results can be made consistent with one another. Half of the diamagnetic polarization is due to the transformation from the guiding center density to the real space density. The other half is due to the fact that, within the drift kinetic ordering assumptions, the guiding center density should be expressed as the gyroaverage of the density in the limit of vanishing Larmor radius. A comprehensive review of polarization is presented to complete the derivation, and then we derive results that are required for the polarization calculations to agree. Expressions for the diamagnetic polarization density are given that are accurate to first order in amplitude and all orders in gyroradius within the gyrokinetic theory for a constant magnetic field. Applications to Maxwell–Boltzmann particle distribution functions (PDFs), including anisotropic temperature, are presented. Local invariants, like total energy and toroidal momentum, do not generate net polarization effects; the electric and thermodynamic polarizations must precisely cancel. In contrast, anisotropic dependence on the magnetic moment generates a net polarization proportional to the temperature anisotropy. Finally, when sources are present, the equilibrium PDF is approximately the ratio of two orbit averages.

adiabatic theory↗
DOE OSTI2021

Importance of gyrokinetic exact Fokker-Planck collisions in fusion plasma turbulence

Gyrokinetic simulations of turbulence are fundamental to understanding and predicting particle and energy loss in magnetic fusion devices. Previous works have used model collision operators with approximate field-particle terms of unknown accuracy and/or have neglected collisional finite Larmor radius effects. This research moves beyond models to demonstrate important corrections using a gyrokinetic Fokker-Planck collision operator with the exact field-particle terms, in realistic simulations of turbulence in magnetically confined fusion plasmas. The exact operator shows significant corrections for temperature-gradient-driven trapped electron mode turbulence and zonal flow damping, and for microtearing modes in a Joint European Torus pedestal under ITER-like wall conditions. Analysis of the corrections using parameter scans motivates an accurate model which closely reproduces the exact results while reducing computational demands.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

GX: a GPU-native gyrokinetic turbulence code for tokamak and stellarator design

GX is a code designed to solve the nonlinear gyrokinetic system for low-frequency turbulence in magnetized plasmas, particularly tokamaks and stellarators. In GX, our primary motivation and target is a fast gyrokinetic solver that can be used for fusion reactor design and optimization along with wide-ranging physics exploration. Here, this has led to several code and algorithm design decisions, specifically chosen to prioritize time to solution. First, we have used a discretization algorithm that is pseudospectral in the entire phase space, including a Laguerre–Hermite pseudospectral formulation of velocity space, which allows for smooth interpolation between coarse gyrofluid-like resolutions and finer conventional gyrokinetic resolutions and efficient evaluation of a model collision operator. Additionally, we have built GX to natively target graphics processors (GPUs), which are among the fastest computational platforms available today. Finally, we have taken advantage of the reactor-relevant limit of small $\rho _*$ by using the radially local flux-tube approach. In this paper we present details about the gyrokinetic system and the numerical algorithms used in GX to solve the system. We then present several numerical benchmarks against established gyrokinetic codes in both tokamak and stellarator magnetic geometries to verify that GX correctly simulates gyrokinetic turbulence in the small $\rho _*$. Moreover, we show that the convergence properties of the Laguerre–Hermite spectral velocity formulation are quite favourable for nonlinear problems of interest. Coupled with GPU acceleration, which we also investigate with scaling studies, this enables GX to be able to produce useful turbulence simulations in minutes on one (or a few) GPUs and higher fidelity results in a few hours using several GPUs. GX is open-source software that is ready for fusion reactor design studies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Energetic particle physics: Chapter 7 of the special issue: on the path to tokamak burning plasma operation

We review the physics of energetic particles (EPs) in magnetically confined burning fusion plasmas with focus on advances since the last update of the ITER Physics Basis (Fasoli et al 2007 Nucl. Fusion 47 S264). Topics include basic EP physics, EP generation, diagnostics of EPs and instabilities, the interaction of EPs and thermal plasma instabilities, EP-driven instabilities, energetic particle modes (EPMs), and turbulence, linear and nonlinear stability and simulation of EP-driven instabilities and EPMs, 3D effects, scenario optimization strategies based on EP phase-space control, EPs in reduced field scenarios in ITER before DT, and the physics of runaway electrons. We describe the simulation and modeling of EPs in fusion plasmas, including instability drive and damping as well as EP transport, with a range of approaches from first-principles to reduced models, including gyrokinetic simulations, kinetic-MHD models, gyrofluid models, reduced models, and semi-analytical approaches.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Isotope effects on energy transport in the core of ASDEX-Upgrade tokamak plasmas: Turbulence measurements and model validation

Design and operation of future tokamak fusion reactors using a deuterium–tritium 50:50 mix requires a solid understanding of how energy confinement properties change with ion mass. This study looks at how turbulence and energy transport change in L-mode plasmas in the ASDEX Upgrade tokamak when changing ion species between hydrogen and deuterium. For this purpose, both experimental turbulence measurements and modeling are employed. Local measurements of ion-scale (with wavevector of fluctuations perpendicular to the B-field k⊥< 2 cm−1, k⊥ρs< 0.2, where ρs is the ion sound Larmor radius using the deuterium ion mass) electron temperature fluctuations have been performed in the outer core (normalized toroidal flux ρTor=0.65−0.8) using a multi-channel correlation electron cyclotron emission diagnostic. Lower root mean square perpendicular fluctuation amplitudes and radial correlation lengths have been measured in hydrogen vs deuterium. Measurements of the cross-phase angle between a normal-incidence reflectometer and an ECE signal were made to infer the cross-phase angle between density and temperature fluctuations. The magnitude of the cross-phase angle was found larger (more out-of-phase) in hydrogen than in deuterium. TRANSP power balance simulations show a larger ion heat flux in hydrogen where the electron-ion heat exchange term is found to play an important role. These experimental observations were used as the basis of a validation study of both quasilinear gyrofluid trapped gyro-Landau fluid-SAT2 and nonlinear gyrokinetic GENE codes. Linear solvers indicate that, at long wavelengths (k⊥ρs<1), energy transport in the deuterium discharge is dominated by a mixed ion-temperature-gradient (ITG) and trapped-electron mode turbulence while in hydrogen transport is exclusively and more strongly driven by ITG turbulence. The Ricci validation metric has been used to quantify the agreement between experiments and simulations taking into account both experimental and simulation uncertainties as well as four different observables across different levels of the primacy hierarchy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗