Search NASA⌕ Search

SEARCH · Search NASA

Results for “gyrokinetic”

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

Simultaneous reproduction of experimental profiles, fluxes, transport coefficients, and turbulence characteristics via nonlinear gyrokinetic profile predictions in a DIII-D ITER similar shape plasma

Experimental conditions obtained on the DIII-D tokamak in the ITER Similar Shape (ISS) have been compared extensively with nonlinear gyrokinetic simulation using the CGYRO code [Candy et al., J. Comput. Phys. 324, 73–93 (2016)] with comparisons spanning ion and electron heat fluxes, electron and impurity particle transport, and turbulent fluctuation levels and characteristics. Bayesian optimization techniques [Rodriguez-Fernandez et al., Nucl. Fusion 62(7), 076036 (2022)], combined with nonlinear gyrokinetics, have been used to obtain simultaneously Qi, Qe, and Γe flux-matched profiles that are found to be in good agreement with experimental profile measurements. Synthetic diagnostics were used to compare measured beam emission spectroscopy and correlation electron cyclotron emission turbulent fluctuations with nonlinear simulation. Although some disagreements exist, nonlinear simulations are found to be in generally good agreement with measured fluctuation levels, spectral shapes, and measured radial trends in low-k δne/ne and δTe/Te. Low (Li and C) and mid-Z (Ca) impurity transport was also compared with these flux-matched simulations. Fully stripped, low-Z impurities are well reproduced by the gyrokinetic modeling while clear disagreement exists in comparisons with mid-Z impurities. Nonlinear gyrokinetic investigation into the Z dependence of impurity transport in the ISS conditions is also performed, demonstrating clear trends of impurity diffusion with impurity Z (both D∝Z and D∝1/Z) that vary with the radial location studied. These trends are shown to result from the local dominance of ion temperature gradient or ∇n driven trapped electron mode turbulence and may contribute to the disagreement between simulation and experiment in mid-Z impurity transport. The results of this work represent one of the most complete validation studies of the gyrokinetic model performed to date and provide an example of new capabilities for predicting performance in future fusion devices.

Physics↗
DOE OSTI2021

Discovering exact, gauge-invariant, local energy-momentum conservation laws for the electromagnetic gyrokinetic system by high-order field theory on heterogeneous manifolds

Gyrokinetic theory is arguably the most important tool for numerical studies of transport physics in magnetized plasmas. However, exact local energy-momentum conservation laws for the electromagnetic gyrokinetic system have not been found despite the continuous effort. Without such local conservation laws, energy and momentum can be instantaneously transported across spacetime, which is unphysical and casts doubt on the validity of numerical simulations based on the gyrokinetic theory. The standard Noether procedure for deriving conservation laws from corresponding symmetries does not apply to gyrokinetic systems because the gyrocenters and electromagnetic fields reside on different manifolds. To overcome this difficulty, we develop a high-order field theory on heterogeneous manifolds for classical particle-field systems and apply it to derive exact, local conservation laws, in particular the energy-momentum conservation laws, for the electromagnetic gyrokinetic system. A weak Euler-Lagrange equation is established to replace the standard Euler-Lagrange equation for the particles. It is discovered that an induced weak Euler-Lagrange current enters the local conservation laws. And it is the new physics captured by the high-order field theory on heterogeneous manifolds. A recently developed gauge-symmetrization method for high-order electromagnetic field theories using the electromagnetic displacement-potential tensor is applied to render the derived energy-momentum conservation laws electromagnetic gauge-invariant.

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 OSTI2022

Toward the core-edge coupling of delta-f and total-f gyrokinetic models

The coupling of core reduced delta-f and edge total-f gyrokinetic models could enable a significant speed-up of the whole device modeling simulations. It is thus interesting to compare these models in the core where they will be coupled. The main difference between these models is the presence of the zeroth order term on the right-hand side of the total-f gyrokinetic equation. This term is responsible for driving a radial electric field that is also associated with Geodesic acoustic mode-like oscillations. To investigate the coupling between these models, the subsequent large-scale gyrokinetic 3D turbulence simulations have been performed starting from an already saturated gyrokinetic axisymmetric equilibria (collisions are not included). This work has also been extended to couple different axisymmetric and turbulent models available in XGC. It is shown that the radial electric field and its drive have to be taken care of consistently while coupling different models together.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Toward continuum gyrokinetic study of high-field mirrors

High-temperature superconducting (HTS) magnetic mirrors under development exploit strong fields with high mirror ratio to compress loss cones and enhance confinement and may offer cheaper, more compact fusion power plant candidates. This new class of devices could exhibit largely unexplored interchange and gradient-driven modes. Such instabilities, and methods to stabilize them, can be studied with gyrokinetics, given the strong magnetization and prevalence of kinetic effects. Our focus here is to (a) determine if oft-used gyrokinetic models for open field lines produce the electron-confining (Pastukhov) electrostatic potential and (b) examine and address challenges faced by gyrokinetic codes in studying HTS mirrors. Here, we show that a one-dimensional limit of said models self-consistently develops a potential qualitatively approaching the analytical Pastukhov level. Additionally, we describe the computational challenges of studying high mirror ratios with open field line gyrokinetic solvers and offer a force softening method to mitigate small time steps needed for time integration in colossal magnetic field gradients produced by HTS coils, providing a 19X speedup.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

A new synthetic correlation electron cyclotron emission diagnostic for validating nonlinear gyrokinetic simulations of electron temperature turbulence

To validate nonlinear gyrokinetic simulations of electron temperature turbulence, the experimental correlation electron cyclotron emission (CECE) measurements are to be compared using a synthetic CECE diagnostic, which generates modeled CECE measurement quantities by implementing realistic measurement parameters (e.g., spatial and wavenumber resolutions, radial location, etc.) to nonlinear gyrokinetic simulations. In this work, we calculate the radial and vertical spatial and wavenumber transfer functions, which are defined by the electron cyclotron emission emissivity radial profile and vertical probing antenna pattern, respectively. These transfer functions are applied to nonlinear gyrokinetic simulations of electron temperature turbulence using the continuum gyrokinetic code. A simultaneous comparison of the experimental electron temperature turbulence power spectrum and root-mean-square (RMS) level, as well as the radial correlation length with the new synthetic CECE diagnostic at a core location ρ ∼ 0.75 in an L-mode DIII-D tokamak plasma, is presented. Here, the preliminary result shows that the synthetic CECE output underestimates the RMS level by ∼42% and overestimates the radial correlation length by ∼40%.

Cyclotron radiation↗
DOE OSTI2024

Prediction of performance and turbulence in ITER burning plasmas via nonlinear gyrokinetic profile prediction

Burning plasma performance, transport, and the effect of hydrogen isotope (H, D, D-T fuel mix) on confinement has been predicted for ITER baseline scenario (IBS) conditions using nonlinear gyrokinetic profile predictions. Accelerated by surrogate modeling (Rodriguez-Fernandez et al 2022 Nucl. Fusion 62 076036), high fidelity, nonlinear gyrokinetic simulations performed with the CGYRO code (Candy et al 2016 J. Comput. Phys. 324 73), were used to predict profiles of T i , T e , and n e while including the effects of alpha heating, auxiliary power (NBI + ECH), collisional energy exchange, and radiation losses inside of $r/a$ = 0.9. Predicted profiles and resulting energy confinement are found to produce fusion power and gain that are approximately consistent with mission goals ($P_\textrm{fusion} = 500$ MW at Q = 10) for the baseline scenario and exhibit energy confinement that is within 1σ of the H-mode energy confinement scaling. The power of the surrogate modeling technique is demonstrated through the prediction of alternative ITER scenarios with reduced computational cost. These scenarios include conditions with maximized fusion gain and an investigation of potential resonant magnetic perturbation (RMP) effects on performance with a minimal number of gyrokinetic profile iterations required (3–6). These predictions highlight the stiff ITG nature of the core turbulence predicted in the ITER baseline and demonstrate that $Q \gt$ 17 conditions may be accessible by reducing auxiliary input power while operating in IBS conditions. Prediction of full kinetic profiles allowed for the projection of hydrogen isotope effects around ITER baseline conditions. The gyrokinetic fuel ion species was varied from H, D, and 50/50 D-T and kinetic profiles were predicted. Results indicate that a weak or negligible isotope effect will be observed to arise from core turbulence in IBS conditions. The resulting energy confinement, turbulence, and density peaking, and the implications for ITER operations will be discussed.

gyrokinetics↗
DOE OSTI2021

Gyrokinetic simulation of turbulence and transport in the SPARC tokamak

The turbulence and transport expected in the SPARC tokamak Primary Reference Discharge (PRD) [P. Rodriguez-Fernandez et al., J. Plasma Phys. 86, 865860503 (2020)] have been investigated with the gyrokinetic code CGYRO [J. Candy et al., J. Comput. Phys. 324, 73–93 (2016)]. Linear and nonlinear simulations that focus on ion (kθρs<1.0) and electron-scale (kθρs>1.0) turbulence were used to probe the nature of the turbulence and the resulting transport in the fusion core. It is found that in the SPARC PRD, ion temperature gradient (ITG) turbulence is expected to dominate transport over most of the profile with some potential trapped electron mode impact in the near edge. Stiff turbulence is observed over a part of the plasma core such that SPARC's ion temperature profile will likely be pinned to just above the critical gradient for ITG. The role of electromagnetic turbulence, rotation, and electron-scale turbulence was investigated to provide some insight into the physics required to accurately predict SPARC performance via gyrokinetics. Additionally, predictions of impurity peaking for potential low- and high-Z SPARC first-wall materials are probed using ion-scale simulation. The dominance of low-k turbulence in SPARC provides a potential opportunity for more tractable prediction of plasma profiles using nonlinear gyrokinetics. This work is the first step toward full gyrokinetic profile prediction of SPARC kinetic profiles and the resulting fusion power and plasma gain.

Howard, N. T. (ORCID:0000000287876309)↗
DOE OSTI2020

Applying the Gyrokinetic Formulation of Magnetized Plasma Dynamics to Z-pinch plasmas

This project adapts and applies the gyrokinetic Eulerian code COGENT to an emerging area of highly-magnetized Z-pinch plasma systems. Simulation efforts for these plasmas currently use either magneto-hydrodynamic (MHD) fluid codes, which suffer from fidelity being restricted to narrow plasma conditions, or fully kinetic particle-in-cell (PIC) methods, which suffer from the run-time needed to resolve the very small spatial and temporal scales required for the plasma conditions. The progress with achieving detailed theoretical understanding of Z-pinch plasma systems can be greatly expedited using some combination of gyrokinetic and drift-ordered MHD models, developed in the tokamak community and implemented in the COGENT code. These formulations are applicable to a much larger parameter space than standard MHD models, while providing a substantial speed-up in run-time over fully kinetic codes. In particular, the gyrokinetic formalism is successfully applied to a flow-stabilized Z pinch (FSZP) system – a promising fusion concept where a sheared axial flow is used to provide stability against large-scale ideal-MHD modes. Results of gyrokinetic simulations with COGENT agree well with those from fully kinetic simulations, while requiring ~500x less run-time.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Fluid and gyrokinetic turbulence in open field-line, helical plasmas

Two-fluid Braginskii codes have simulated open-field line turbulence for over a decade, and only recently, it has become possible to study these systems with continuum gyrokinetic codes. This work presents a first-of-its-kind comparison between fluid and (long-wavelength) gyrokinetic models in open field-lines, using the GDB and Gkeyll codes to simulate interchange turbulence in the Helimak device at the University of Texas [T. N. Bernard et al., Phys. Plasmas 26, 042301 (2019)]. Partial agreement is attained in a number of diagnostic channels when the GDB sources and sheath boundary conditions (BCs) are selected carefully, especially the heat-flux BCs that can drastically alter the temperature. The radial profile of the fluctuation levels is qualitatively similar and quantitatively comparable on the low-field side, although statistics such as moments of the probability density function and the high-frequency spectrum show greater differences. This comparison indicates areas for future improvement in both simulations, such as sheath BCs, and improvements in GDB like particle conservation and spatially varying thermal conductivity, in order to achieve better fluid-gyrokinetic agreement and increase fidelity when simulating experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Validation of gyrokinetic simulations in NSTX and projections for high-k turbulence measurements in NSTX-U

An extensive validation effort performed for a modest-beta NSTX NBI-heated H-mode discharge predicts that electron thermal transport can be entirely explained by electron-scale turbulence fluctuations driven by the electron temperature gradient mode (ETG), both in conditions of strong and weak ETG turbulence drive. Thermal power-balance estimates computed by TRANSP as well as the shape of the high-k density fluctuation wavenumber spectrum and the fluctuation level ratio between strongly driven and weakly driven ETG-turbulence conditions can be matched by nonlinear gyrokinetic simulations and a synthetic diagnostic for high-k scattering. Linear gyrokinetic simulations suggest that the ion-scale instability in the weak ETG condition is close to the critical threshold for the kinetic ballooning mode instability, and nonlinear ion-scale gyrokinetic simulations show that turbulence might be in a state reminiscent of a Dimits' shift regime, opening speculation on the role that ion-scale turbulence might play for the weak ETG condition. A simulation that matched all experimental constraints is chosen to project high-k turbulence spectra in NSTX-U, revealing that the new high-k system [R. Barchfeld et al., Rev. Sci. Instrum. 89, 10C114 (2018)] should be sensitive to density fluctuations from radially elongated streamer structures. In this study, two schemes are designed to characterize the radial and poloidal wavenumber dependence of the density fluctuation wavenumber power spectrum around the streamer peak, suggesting future high-k fluctuation measurements could be sensitive to an asymmetry in the k r spectrum introduced due to the presence of strong background flow shear.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Tokamak ITG-KBM transition benchmarking with the mixed variables/pullback transformation electromagnetic gyrokinetic scheme

Electromagnetic gyrokinetic simulation of high temperature plasma is required to predict confinement in magnetic fusion devices and has posed challenges for existing codes. In this paper, we demonstrate successful global gyrokinetic simulation of the ion temperature gradient-driven mode-kinetic ballooning mode transition in a toroidal fusion plasma test case using the mixed variables/pullback transformation (MV/PT) scheme with the particle-in-cell codes XGC and ORB5, and compare to results from a conventional continuum code from the literature. Furthermore, the MV/PT scheme combines explicit time integration with mitigation of the well-known electromagnetic gyrokinetic “cancelation problem.” We calculate eigenmodes in the electrostatic and parallel vector potentials, and find good agreement in growth rate, real frequency, and the normalized plasma pressure of mode transition.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Gyrokinetic and extended-MHD simulations of a flow shear stabilized Z-pinch experiment

Axisymmetric (m = 0) gyrokinetic and extended-MHD simulations of a sheared-flow Z-pinch plasma are performed with the high-order finite volume code COGENT. The present gyrokinetic model solves the long-wavelength limit of the gyrokinetic equation for both ion and electron species coupled to the electrostatic gyro-Poisson equation for the electrostatic potential. The extended-MHD model is electromagnetic and includes the effects of the gyro-viscous pressure tensor, diamagnetic electron and ion heat fluxes, and generalized Ohm's law. Additionally, a prominent feature of this work is that the radial profiles for the plasma density and temperature are taken from the fusion Z-pinch experiment (FuZE), and the magnetic field profile is obtained as a solution of the MHD force balance equation. Such an approach allows to address realistic plasma parameters and provide insights into the current and planned experiments. In particular, it is demonstrated that the radial profiles play an important role in stabilization, as the embedded guiding center (E×B ) drift has a strong radial shear, which can contribute to the Z-pinch stabilization even in the absence of the fluid flow shear. The results of simulations for the FuZE plasma parameters show a decrease in the linear growth rate with an increase in the flow shear; however, full stabilization in the linear regime is not observed even for large (comparable to the Alfvén velocity) radial variations of the axial flow. Nonlinear stability properties of the FuZE plasmas are also studied, and it is found that profile broadening can have a pronounced stabilizing effect in the nonlinear regime.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Hybrid-gyrokinetic simulations of low- n toroidal Alfvén eigenmodes using gKPSP

Here, we report a benchmark study of toroidal Alfvén eigenmode (TAE) simulation using the hybrid-gyrokinetic code GyroKinetic Plasma Simulation Program (gKPSP). A simulation capability for energetic particles based on the gyrokinetic δf method has been newly implemented in the gKPSP code. Benchmark simulations have been performed in both circular and realistic tokamak geometries. Good agreement has been found with previously reported results, demonstrating the new capability of the gKPSP code. We have investigated the effects of the distribution function on TAE stability by examining both isotropic and anisotropic slowing-down distributions of energetic particles. The slowing-down distribution produces a higher linear growth rate than a Maxwellian distribution, while a growth rate scan with the anisotropy shows an opposite trend. This can be attributed to competition between Landau damping and the linear drive, which are correlated with the fraction of resonant passing particles and their distribution in phase space.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Core-edge modeling of gyrokinetic turbulence by coupling the delta-f and total-f models in the XGC code

Whole device gyrokinetic turbulence is modeled using a delta-f model in the core coupled to a total-f model in the edge. This new model is implemented in the gyrokinetic code X-Point Included Gyrokinetic Code (XGC) and works within a single simulation. It employs a canonical Maxwellian background in the core and the usual local Maxwellian in the edge. These different core and edge backgrounds are coupled over an overlap region. A thermal-bath operator is used to keep δf small in the core.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Verification of electromagnetic simulation capabilities in global gyrokinetic particle-in-cell code GTS

Recently, the numerical scheme presented by Mishchenko et al. enabled explicit gyrokinetic simulations of low-frequency electromagnetic instabilities in tokamaks at experimentally relevant values of plasma β⁠. This scheme resolved the long-standing cancellation problem that previously hindered gyrokinetic particle-in-cell code simulations of magnetohydrodynamic phenomena with inherently small parallel electric fields. Moreover, the scheme did not employ approximations that eliminate critical tearing-type instabilities. Here, we report on the implementation of this numerical scheme in the global gyrokinetic particle-in-cell code GTS. This implementation allows for a more complete and accurate picture of interaction between small scale turbulence and MHD modes in tokamaks. Additionally, we present a comprehensive set of verification simulations of numerous electromagnetic instabilities relevant to present-day tokamaks. These simulations encompass the kinetic ballooning mode, the internal kink mode, the tearing mode, the micro-tearing mode, and the toroidal Alfven eigenmode destabilized by energetic ions, which are all instrumental in understanding tokamak physics. We will also showcase the preliminary nonlinear simulations of kinetic ballooning instabilities and (2,1) island formation due to tearing mode instability. These simulations validate the accuracy of the scheme implementation and pave the way for studying how these instabilities affect plasma confinement and performance.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Linearised Fokker–Planck collision model for gyrokinetic simulations

We introduce a gyrokinetic, linearised Fokker–Planck collision model that satisfies conservation laws and is accurate at arbitrary collisionalities. The differential test-particle component of the operator is exact; the integral field-particle component is approximated using a spherical harmonic and a modified Laguerre polynomial expansion developed by Hirshman and Sigmar (1976 Phys. Fluids 19 1532). The numerical methods of the implementation in the δf-gyrokinetic code stella (Barnes et al 2019 J. Comput. Phys. 391 365–80) are discussed, and conservation properties of the operator are demonstrated. The collision model is then benchmarked against the collision model of the gyrokinetic solver GS2 in the limiting cases of a reduced test-particle collision operator and energy- and momentum-conserving operator. The accuracy of the full collision model is investigated by solving the parallel Spitzer-Härm problem for the transport coefficients. It is shown that retaining collisional energy flux and higher-order terms in the field-particle operator reduces errors in the transport coefficients from 10%–25% for a simple momentum- and energy-conserving model to under 1%.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗