Search NASA⌕ Search

SEARCH · Search NASA

Results for “gyrokinetic theory”

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

Magnetic geometry effects on turbulent density fluctuations in Wendelstein 7-X

Line-integrated density fluctuation measurements from the core of the optimised stellarator Wendelstein 7-X (W7-X) are studied in the context of magnetic field geometry changes. The amplitude as well as the spectral distribution of fluctuations is affected by varying the magnetic field configuration. Weaker turbulent fluctuations are observed in configurations with higher rotational transform and lower magnetic mirror depth. Simulations with the gyrokinetic codes stella and GENE-3D are used to investigate the role of magnetic geometry for the experimental observations. The stabilising effect of the high rational transform configuration is reproduced by simulations and can be associated with stronger shaping of the magnetic field. This leads to increased finite Larmor radius and local shear effects, which limit the toroidal extent of ion temperature gradient (ITG) modes. For W7-X configurations with higher magnetic mirror depth, gyrokinetic theory and simulations predict weaker fluctuations, since they are closer to having the maximum-J property, which has a stabilising effect on trapped-electron mode and ITG. The opposite trend is observed experimentally, which can therefore not be explained by geometry effects on turbulence alone.

Wendelstein-7X↗
DOE OSTI2024

First Measurement of Drift-Alfvén Wave Polarization in Magnetically Confined Fusion Plasmas

Polarization of drift-Alfvén waves, defined as the ratio of electrostatic to electromagnetic fluctuations, has remained unmeasurable in fusion plasmas for decades, despite its pivotal role in understanding wave dynamics and their impact on plasmas. We report the first measurements of drift-Alfvén wave polarization in a hot, magnetically-confined plasma. Here, the breakthrough is enabled by a novel methodology developed from gyrokinetic theory, utilizing fluctuations of electron temperature and density. Analysis of data from the DIII-D tokamak reveals that the waves above the geodesic acoustic mode frequency exhibit dominant electromagnetic polarization, whereas lower frequency waves show a mix of electromagnetic and electrostatic polarization, indicating a strong coupling between shear Alfvén waves and drift-acoustic waves.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Turbulence Transition in Magnetically Confined Hydrogen and Deuterium Plasmas

In this study, we discovered a turbulence transition in a large helical device. The turbulence level and turbulence-driven energy transport decrease to a specific transition density and increase above it. The ruling turbulences below and above the transition density were ion-temperature gradient (ITG) and resistive-interchange (RI) turbulences, consistent with the predictions of gyrokinetic theory and two-fluid MHD model, respectively. Isotope experiments on hydrogen (H) and deuterium (D) clarified the role of transitions. In the ITG regime, turbulence levels and energy transport were comparable in the H and D plasmas. In contrast, in the RI regime, they were clearly suppressed in the D plasma. The results provide crucial knowledge for understanding isotope effects and future optimization of stellarator and heliotron devices. Published by the American Physical Society 2024

Physics↗
DOE OSTI2025

Gyrokinetic Simulation of Tokamak Fusion Plasmas

The final technical report for project DE‐FG02‐08ER54954 “Gyrokinetic Simulation of Tokamak Fusion Plasmas”. This report covers the entire period of the sponsored project. This grant supported in part the salary, travel, and research supplies for the principal investigator Scott E. Parker, graduate research associate Stefan Tirkas, who is a Ph.D. graduate student in the Physics Department. Postdoctoral Scholar Dr. Neeraj Kumar is also collaborating on this project. The main goals of the project are to investigate the nonlinear saturation of electron scale turbulence using gyrokinetic theory and simulation and to investigate impurity transport in the edge pedestal region.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Quasilinear theory and modelling of gyrokinetic turbulent transport in tokamaks

The theory, development, and validation of reduced quasilinear models of gyrokinetic turbulent transport in the closed flux surface core of tokamaks is reviewed. In combination with neoclassical collisional transport, these models are successful in accurately predicting core tokamak plasma temperature, density, rotation, and impurity profiles in a variety of confinement regimes. Refined experimental tests have been performed to validate the predictions of the quasilinear models, probing changes in the dominant gyrokinetic instabilities, as reflected in fluctuation measurements, cross-phases, and transport properties. These tests continue to produce a deeper understanding of the complex mix of instabilities at both electron and ion gyroradius scales.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Dynamic mode decomposition for gyrokinetic eigenmode analysis

Dynamic mode decomposition (DMD) is a post-processing approach to decompose a complex time series into a set of modes via spectral analysis. DMD provides a new and powerful method to recover gyrokinetic drift-wave eigenfrequencies and eigenfunctions based only on the solution of the gyrokinetic-Maxwell initial value problem with almost no added cost to the initial value solver. In the present paper, DMD is applied to the CGYRO gyrokinetic code using a newly-developed CGYRO-DMD post-processor. CGYRO-DMD is numerically efficient, even on a single CPU. It does not set any restrictions on the plasma shape, beta (ratio of the plasma pressure to the magnetic field pressure), collisionality or number of species, and allows one to resolve numerous eigenmodes, even of comparable growth rates. In addition, DMD is not limited to unstable modes, but rather can capture stable and unstable branches simultaneously. In this work, we illustrate the accuracy of DMD through gyrokinetic analysis of mode transition for electromagnetic drift wave instabilities.

drift-wave eigenmodes↗
DOE OSTI2025

Omnigenous stellarators with improved ideal and kinetic ballooning stability

Omnigenity is a property of a magnetic field which ensures confinement of trapped particles. It is a necessary requirement for any high-performance stellarator. After creating an omnigenous equilibrium, one must also ensure reduced transport resulting from kinetic and magnetohydrodynamic (MHD) instabilities. To this end, we leverage the GPU-accelerated DESC optimization suite, which is used to design stable, finite-β omnigenous equilibria with poloidal, toroidal, and helical symmetry, achieving Mercier, ideal ballooning, and as a consequence, improved kinetic ballooning stability. We discover stellarators with second stability, a regime of large pressure gradient where an equilibrium becomes ideal ballooning stable, and demonstrate and explore both using theory and gyrokinetic simulations the connection between ideal and kinetic ballooning stability.

optimization↗
DOE OSTI2026

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)↗
DOE OSTI2024

Turbulence link to L-mode, I-mode, and H-mode confinement in the DIII-D tokamak

Understanding the physics of low-confinement (L-), improved-confinement (I-), and high-confinement (H-) modes is critical for fusion reactors. The finding herein reports observations of two types of turbulence coexisting near the L-mode edge, one magnetohydrodynamic (MHD)-like and another micro-tearing mode (MTM)-like, linked to the H-mode and I-mode confinement in the DIII-D tokamak. Ion-scale magnetic and density turbulence is measured using a Faraday-effect radial-interferometer-polarimeter and beam-emission-spectroscopy (BES). Broadband turbulence spectra of up to ~600 kHz are observed in two discharges where transitions between L-mode, I-mode, and H-mode occurs. Turbulence is found to be inversely correlated with confinement, meaning lower turbulence power at higher confinement. Distinctively, the high-frequency (HF, >~100 kHz) magnetic turbulence power changes by the most (55%) during transitions primarily involving energy confinement change, whereas the low-frequency (LF, <~100 kHz) magnetic and density turbulence power changes by the most (80%) during transitions primarily involving particle confinement change. The LF turbulence amplitude oscillates with and leads to deuterium-alpha emission oscillations before an H-mode. These results imply that HF turbulence mainly affects energy confinement whereas LF turbulence can affect particle confinement. The magnetic and density turbulence exhibits coherence up to 0.6 and cross-phase magnitude close to π/2 in most cases, suggesting they have a common origin in both the LF and HF ranges. BES suggests that LF turbulence resides at the edge (ρ = 0.95) and HF turbulence can be at the outer core (ρ = 0.8) or edge (ρ = 0.95). Comparisons of measurements, theory, and gyrokinetic simulations suggest that HF turbulence is MTM-like in all cases, whereas LF turbulence is more consistent with MHD-like modes and the exact instability might change during transitions—except that a drift-wave origin is possible in a low collisionality H-mode. These results suggest that the H-mode involves suppressed MHD-like turbulence, whereas the I-mode mitigates MTM-like turbulence along with largely unchanged MHD-like turbulence.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Development of a data-driven neural network model for electron thermal transport in NSTX

A data-driven electron thermal transport neural network (ETT-NN) model, trained on TRANSP interpretative analysis results of National Spherical Torus Experiment (NSTX), was developed to enable faster and more accurate ETT computation for spherical tokamaks (STs). The model incorporates both convolutional NNs and recurrent NNs, allowing it to simultaneously account for the spatial and temporal non-localities and multi-scale features of turbulent transport, which have been considered only in a limited manner in conventional models. The model was validated through interpretative analysis and predictive simulations using Tokamak Reactor Integrated Automated Suite for Simulation and Computation, demonstrating relatively high accuracy. Additionally, parameter scans were performed on test discharges known to exhibit specific turbulent modes, such as microtearing mode, trapped electron mode, kinetic ballooning mode, and electron temperature gradient mode. The scanning results revealed that the ETT-NN model exhibits the same trends as those observed in conventional gyrokinetic simulations or theories, while also capturing the global nature of turbulent transport, indicating that the data-driven model accurately reflects the underlying physical characteristics. Furthermore, due to the dimensionless nature of the model, we can feasibly expand its applicability by incorporating data from other devices and uncovering the characteristics of ETT in STs in the future.

NSTX↗
DOE OSTI2024

Potential vorticity conservation for plasma turbulence in an inhomogeneous magnetic field: Theory and implications

The concept and theory of potential vorticity in drift wave turbulence are extended to the case of an inhomogeneous magnetic field. A one-field magnetic potential vorticity conserving equation is derived via the use of conservative gyrokinetics. The similarity between the corresponding systems for drift wave turbulence and shallow water theory is discussed in detail. Zonal flow physics in an inhomogeneous magnetic field is discussed. In particular, a Charney–Drazin type nonacceleration theorem is derived from the novel system, which conserves magnetic potential vorticity. Extensions of the turbulent equipartition theory to the transport of magnetic potential vorticity are proposed.

Physics↗
DOE OSTI2025

Turbulent Heating in Collisionless Low-beta Plasmas: Imbalance, Landau Damping, and Electron–Ion Energy Partition

An understanding of how turbulent energy is partitioned between ions and electrons in weakly collisional plasmas is crucial for modeling many astrophysical systems. Using theory and simulations of a four-dimensional reduced model of low-beta gyrokinetics (the “Kinetic Reduced Electron Heating Model”), we investigate the dependence of collisionless heating processes on plasma beta and imbalance (normalized cross-helicity). These parameters are important because they control the helicity barrier, the formation of which divides the parameter space into two distinct regimes with remarkably different properties. In the first, at lower beta and/or imbalance, the absence of a helicity barrier allows the cascade of injected power to proceed to small (perpendicular) scales, but its slow cascade rate makes it susceptible to significant electron Landau damping, in some cases leading to a marked steepening of the magnetic spectra on scales above the ion Larmor radius. In the second, at higher beta and/or imbalance, the helicity barrier halts the cascade, confining electron Landau damping to scales above the steep “transition-range” spectral break, resulting in dominant ion heating. We formulate quantitative models of these processes that compare well to simulations in each regime, and combine them with results of previous studies to construct a simple formula for the electron–ion heating ratio as a function of beta and imbalance. This model predicts a “winner takes all” picture of low-beta plasma heating, where a small change in the fluctuations' properties at large scales (the imbalance) can cause a sudden switch between electron and ion heating.

Interplanetary turbulence↗
DOE OSTI2024

Gyrokinetic simulation of pedestal degradation correlated with enhanced magnetic turbulence in a DIII-D ELMy H-mode discharge

Gyrokinetic simulation of a dedicated pedestal density ramping-up discharge on DIII-D can reproduce the enhancement of magnetic turbulence in the pedestal, which is identified to be caused by micro-tearing modes (MTMs). Here, an increase of MTM amplitude results in higher electron thermal diffusivity, consistent with experimentally observed lower electron temperature gradient and degraded pedestal height. Gyrokinetic simulation identifies the major cause of MTM enhancement to be the increase of collisionality, which has a significant impact on the MTM intensity and is beyond the description of any (quasi-)linear theory.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Understanding plasma turbulence through exact coherent structures

Plasma turbulence is a key challenge in understanding transport phenomena in magnetically confined plasmas. This work presents a generalized framework to analyze plasma turbulence that utilizes periodic orbit theory. In periodic orbit theory, doubly periodic solutions (coherent structures) of the governing equation(s) serve as building blocks of the considered turbulent dynamics. To illustrate the concept and method, the particularly simple Kuramoto–Sivashinsky (referred to here as LMRT for the original authors: LaQuey, Mahajan, Rutherford, and Tang) trapped-ion mode toy model is used. By applying numerical optimization techniques to the LMRT equation, we extract coherent spacetime patterns that represent the library of allowable fundamental structures of the equation. These structures provide a framework to systematically describe turbulence as a composition of recurrent solutions, revealing an underlying order within chaotic plasma motion. Although illustrated here using the simplified LMRT model for clarity, this framework provides a general strategy that can be extended to more complex and realistic models of plasma turbulence, including gyrokinetic systems. This offers a new method for predicting and potentially controlling transport processes in fusion plasmas by providing a bridge between nonlinear dynamical systems theory and plasma physics in the form of a generalized framework with which to analyze and understand spatially extended nonlinear partial differential equations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Saturation of fishbone instability through zonal flows driven by energetic particle transport in tokamak plasmas

Gyrokinetic and kinetic-MHD simulations are performed for the fishbone instability in the DIII-D discharge #178631, chosen for validation of first-principles simulations to predict the energetic particle (EP) transport in an ITER prefusion baseline scenario. Fishbone modes are found to generate zonal flows, which dominate the fishbone saturation. The underlying mechanisms of the two-way fishbone-zonal flows nonlinear interplay are discussed in details. Numerical and analytical analyses identify the fishbone-induced EP redistribution as the dominant generation mechanism for zonal flows. The zonal flows modify the nonlinear dynamics of phase space zonal structures, which reduces the amount of EPs able to resonate with the mode, leading to an early fishbone saturation. Simulation results including zonal flows agree quantitatively with DIII-D experimental measurements of the fishbone saturation amplitude and EP transport, supporting this novel saturation mechanism by self- generated zonal flows. Moreover, the wave-particle mode-locking mechanism is shown to determine quantitatively the fishbone frequency down-chirping, as evident in GTC simulation results in agreement with predictions from analytical theory. Finally, the&#xD;fishbone-induced zonal flows are possibly responsible for the formation of an ion-ITB in the DIII-D discharge. Based on the low EP transport and the large zonal flow shearing rates associated with the fishbone instability in gyrokinetic simulations of the ITER scenario, it is conjectured that high performance scenarios could be designed in ITER burning plasmas through fishbone-induced ITBs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Generalized kinetic equation for tokamak plasma equilibrium distribution function

A generalised kinetic equation for the equilibrium distribution function in a finite beta, arbitrary tokamak plasma is derived. The equation is correct to second order in ρ/L (ρ is the particle Larmor radius and L is the system size). Resolving finite Larmor radius length scales with no restriction on the ratio of poloidal to total equilibrium magnetic field, B $\vartheta$ /B, it generalises the drift kinetic theory of [Hazeltine Phys. Plasmas 15 (1973) 77] to the limit of B $\vartheta$ /B ~ 1 (e.g. to ensure validity for spherical tokamaks). Two cases are considered. The first provides the equilibrium distribution function, consistent with the generalised gyrokinetic formalism of [Dudkovskaia et al. Plasma Phys. Control. Fusion 65 (2023) 045010], derived specifically to capture neoclassical equilibrium currents in gyrokinetic stability analyses in strong gradient regions. The second assumes short length scales in the direction perpendicular to the magnetic field, which can occur as a result of small coherent magnetic structures in the plasma, such as neoclassical tearing mode magnetic islands close to threshold. This then extends the drift island equations of [Dudkovskaia et al. Nucl. Fusion 63 (2023) 016020] for the plasma response to magnetic islands to a spherical tokamak plasma configuration. In conclusion, resolving ρ ~ ρ $\vartheta$ (or B $\vartheta$ ~ B), where ρ $\vartheta$ is the particle poloidal Larmor radius, is also expected to influence calculations of the magnetic island propagation frequency and the associated contributions to the island onset conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Global linear and nonlinear gyrokinetic simulations of tearing modes

Abstract To better understand multi-scale interactions between global tearing modes and microturbulence in the Madison Symmetric Torus (MST) reversed-field pinch (RFP), the global gyrokinetic code Gene is modified to describe global tearing mode instability via a shifted Maxwellian (SM) distribution consistent with experimental equilibria. The implementation of the SM is tested and benchmarked by comparisons with different codes and models. Good agreement is obtained in code-code and code-theory comparisons. Linear stability of tearing modes of a non-reversed MST discharge is studied. A collisionality scan is performed to the lowest order unstable modes ( n = 5, n = 6) and shown to behave consistently with theoretical scaling. The nonlinear evolution is simulated, and saturation is found to arise from mode coupling and transfer of energy from the most unstable tearing mode to small-scale stable modes mediated by the m = 2 tearing mode. The work described herein lays the foundation for nonlinear simulation and analysis of the interaction of tearing modes and gyroradius-scale instabilities in RFP plasmas.

Physics↗