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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Global gyrokinetic simulations of the impact of magnetic island on ion temperature gradient driven turbulence

The effect of island width on the multi-scale interactions between magnetic island (MI) and ion temperature gradient (ITG) turbulence has been investigated based on the global gyrokinetic approach. It is found that the coupling between the island and turbulence is enhanced when the MI width (w) becomes larger. A vortex flow that is highly sensitive to the width of the MI can be triggered, ultimately resulting in a potent E × B shear flow and a consequent reduction in turbulent transport. The shearing rate induced by the vortex flow is minimum at the O-point while it is maximum at both of the two reconnection points of the island, i.e. the X-points, regardless of the island width. There exists a nonmonotonic relationship between zonal flow (ZF) amplitude and island width, showing that the ZF is partially suppressed by medium-sized MIs whereas enhanced in the case of large island. A larger MI can tremendously damage the ITG mode structure, resulting in higher turbulent transport at the X-point whereas a lower one at the O-point, respectively. Such phenomenon will be less distinct at very small island widths below w/a ~ 8% ≈ 12ρ i (a is the minor radius and ρ i the ion gyroradius), where it shows that turbulence near the X-point is hardly affected although it is still suppressed inside the island. Furthermore, the influence of different island sizes on turbulence transport level is also discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Global gyrokinetic simulations of electrostatic microturbulent transport in LHD stellarator with boron impurity

Abstract Global gyrokinetic simulations of electrostatic microturbulent transport for discharge # 166256 of the Large Helical Device stellarator in the presence of boron impurity show the co-existence of the ion temperature gradient (ITG) turbulence and trapped electron mode (TEM) turbulence before and during boron powder injection. ITG turbulence dominates in the core, whereas TEM dominates near the edge, consistent with the experimental observations. Linear TEM frequency increases from ∼ 80 kHz to ∼ 100 kHz during boron injection, and ITG frequency decreases from ∼ 20 kHz to ∼ 13 kHz, consistent with the experiments. The poloidal wave number spectrum is broad for both ITG (0–0.5 mm −1 ) and TEM (0–0.25 mm −1 ). The nonlinear simulations with boron impurity show a reduction in the heat conductivity compared to the case without boron. The comparison of the nonlinear transport before and during boron injection shows that the ion heat transport is substantially reduced in the region where the TEM is dominant. However, the average electron heat transport throughout the radial domain and the average ion heat transport in the region where the ITG is dominant are similar. The simulations with boron show the effective heat conductivity values qualitatively agree with the estimate obtained from the experiment.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Stability and transport of gyrokinetic critical pedestals

Abstract A gyrokinetic threshold model for pedestal width–height scaling prediction is applied to multiple devices. A shaping and aspect ratio scan is performed on National Spherical Torus Experiment (NSTX) equilibria, finding Δ ped = 0.92 A 1.04 κ − 1.24 0.38 δ β θ , ped 1.05 for the wide-pedestal branch with pedestal width Δ ped , aspect ratio A , elongation κ , triangularity δ , and normalized pedestal height β θ , ped . The width–transport scaling is found to vary significantly if the pedestal height is varied either with a fixed density or fixed temperature, showing how fueling and heating sources affect the pedestal density and temperature profiles for the kinetic-ballooning-mode (KBM) limited profiles. For an NSTX equilibrium, at fixed density, the wide branch is Δ ped = 0.028 ( q e / Γ e − 1.7 ) 1.5 ∼ η e 1.5 and at fixed temperature Δ ped = 0.31 ( q e / Γ e − 4.7 ) 0.85 ∼ η e 0.85 , where q e and Γ e are turbulent electron heat and particle fluxes and η e = ∇ ln ⁡ T e / ∇ ln ⁡ n e for an electron temperature T e and density n e . Pedestals close to the KBM limit are shown to have modified turbulent transport coefficients compared to the strongly driven KBMs. The role of flow shear is studied as a width–height scaling constraint and pedestal saturation mechanism for a standard and lithiated wide pedestal discharge. Finally, the stability, transport, and flow shear constraints are combined and examined for an NSTX experiment.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Gyrokinetic prediction of core tungsten peaking in a WEST plasma with nitrogen impurities

Tungsten peaking is predicted in the core of a WEST plasma with total-f gyrokinetic simulations, including both collisional and turbulent transport. This prediction is validated with a synthetic diagnostic of the bolometry. Although nitrogen impurities are shown to reduce the neoclassical peaking of tungsten on-axis, the overall tungsten peaking increases when nitrogen impurities are present, as they reduce the turbulence screening off-axis. This finding is important for the plasma current ramp-up phase of ITER, where light impurities seeding will be desirable to achieve low temperatures at the plasma-facing components and reduce tungsten sputtering. It provides further argument for applying early ECRH heating to maintain margins on the core power balance. The neoclassical peaking factor is cross-verified between XGC and FACIT. The heat flux at separatrix and the heat load width are modeled by XGC and compared to WEST data.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Gyrokinetic profile prediction and validation of a negative triangularity plasma in ASDEX Upgrade

In this work, gyrokinetic simulations are performed with the CGYRO code (Candy et al 2016 J. Comput. Phys. 324 73–93) for a negative triangularity H-mode plasma in ASDEX Upgrade, and compared with experimental measurements. The PORTALS framework (Rodriguez-Fernandez et al 2024 Nucl. Fusion 64 076034) is used to accelerate the prediction of kinetic profiles for this plasma, using surrogate modeling and Bayesian optimization. Ion heat flux, electron heat flux, and electron particle flux are simultaneously matched across the simulated radial regime of the plasma (normalized radius $r/a = 0.35-0.90$), and the resulting ion temperature, electron temperature, and electron density profiles match well with the experimental profile data within this radial range. A synthetic Correlation Electron Cyclotron Emission diagnostic is applied to find well-matched electron temperature fluctuation properties between simulation and experiment. The flux-matched profiles provide a basis for investigation of the turbulence nature across the plasma radius, revealing the dominance of Trapped Electron Mode turbulence at $r/a = 0.35$, the dominance of Ion Temperature Gradient turbulence at $r/a = 0.55$, 0.75, and 0.83, and an instability boundary at $r/a = 0.90$.

gyrokinetic simulation↗
DOE OSTI2025

Gyrokinetic Particle Simulations of Core Turbulence in Fusion Plasmas

In this project, we have further developed GTC gyrokinetic particle simulations to study plasma confinement properties in the stellarators and tokamaks with 3D magnetic fields, where both neoclassical and turbulent transport are self-consistently treated on an equal footing.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

A learned closure method applied to phase mixing in a turbulent gradient-driven gyrokinetic system in simple geometry

We present a new method for formulating closures that learn from kinetic simulation data. We apply this method to phase mixing in a simple gyrokinetic turbulent system – temperature-gradient-driven turbulence in an unsheared slab. The closure, called the learned multi-mode (LMM) closure, is constructed by, first, extracting an optimal basis from a nonlinear kinetic simulation using singular value decomposition. Subsequent nonlinear fluid simulations are projected onto this basis and the results are used to formulate the closure. We compare the closure with other closures schemes over a broad range of the relevant two-dimensional parameter space (collisionality and gradient drive). We find that the turbulent kinetic system produces phase-mixing rates much lower than the linear expectations, which the LMM closure is capable of capturing. We also compare radial heat fluxes. A Hammett–Perkins closure, generalized to include collisional effects, is quite successful throughout the parameter space, producing ${\sim }14\,\%$ root-mean-square (r.m.s.) error. The LMM closure is also very effective: when trained at three (two) points (in a 35 point parameter grid), the LMM closure produces $8\,\%$ ( $12\,\%$ ) r.m.s. errors. The LMM procedure can be readily generalized to other closure problems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Constructing field-aligned coordinate systems for gyrokinetic simulations of tokamaks in X-point geometries

Structures in tokamak plasmas are elongated along the direction of the magnetic field and short in the directions perpendicular to the magnetic field. Many tokamak simulation codes take advantage of this by using a field-aligned coordinate system. However, field-aligned coordinate systems have a coordinate singularity at magnetic X-points where the poloidal magnetic field vanishes, which makes it difficult to use field-aligned coordinate systems when simulating the core and scrape-off layer simultaneously. Here, we present an algorithm for grid generation and computing geometric quantities in a standard field-aligned coordinate system that avoids the singularity and allows one to conduct two-dimensional gyrokinetic axisymmetric simulations in X-point geometries. Convergence tests of advection, boundary value problems and geometric quantities all show greater than first-order convergence even in the vicinity of the X-point. We also demonstrate the geometric consistency of our algorithm with an example simulation of the spherical tokamak for energy production, which shows machine-precision particle conservation.

fusion plasma↗
DOE OSTI2020

Nonlinear global gyrokinetic delta- f turbulence simulations in a quasi-axisymmetric stellarator

We use the global gyrokinetic stellarator code XGC-S to look at turbulence physics in geometry of the proposed quasi-axisymmetric stellarator, QUASAR. In this article, the XGC-S code is first verified for nonlinear simulations of ion temperature gradient-driven turbulence by a successful benchmark with the established nonlinear tokamak code XGC1 in a CYCLONE-like tokamak case. We then apply the XGC-S code to QUASAR for nonlinear turbulence physics. We find that the highly poloidally localized linear mode gives way to a more evenly distributed nonlinear turbulent structure. Furthermore, the radial extent of the mode increases and occupies a large fraction of the volume, despite the initially localized temperature profile, as the profiles relax.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Examination of synthetic gas puff imaging diagnostic data from a gyrokinetic turbulence code

A synthetic gas puff imaging (GPI) diagnostic has been developed for the purpose of validating the three-dimensional gyrokinetic turbulence code XGC. The synthetic diagnostic is described and applied to XGC simulations of two Alcator C-Mod discharges. The turbulence characteristics deduced from the resulting simulated GPI frames, using analysis techniques similar to those applied to experimental data, are compared with locally derived characteristics extracted directly from the XGC output. The comparison of the two is shown to be potentially impacted by misalignment between the GPI view and the magnetic field, the dependence of the light emission on the electron density and temperature, and spatial and temporal variations in the neutral gas cloud induced by the turbulent plasma fluctuations. We conclude that quantitative and, in some cases even qualitative, validation of turbulence simulations need to account for these effects. While we cannot directly compare our results with experimental data due to the absence of high quality GPI data from the shots and times simulated by XGC, we do relate the overall characteristics obtained from the synthetic GPI analysis to published Alcator C-Mod GPI data.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Investigating shear flow through continuum gyrokinetic simulations of limiter biasing in the Texas Helimak

Previous limiter-biasing experiments on the Texas Helimak, a simple magnetized torus, have been inconclusive on the effect of flow shear on turbulence levels. To explore this, the first gyrokinetic simulations of limiter biasing in the Helimak using the plasma physics code Gkeyll have been carried out, and results are presented here. For the scenarios considered, turbulence is mostly driven by the interchange instability, which depends on gradients of steady state density profiles. An analysis of both experimental and simulation data demonstrates that shear rates are mostly less than than local linear growth rates, and not all requirements for shear stabilization are met. Yet, the mostly vertical shear flow has an important effect on bulk transport and experimental steady-state density profiles, and changes to the gradients correspond to changes in turbulence levels.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Spatial core-edge coupling of the particle-in-cell gyrokinetic codes GEM and XGC

Two existing particle-in-cell gyrokinetic codes, GEM for the core region and XGC for the edge region, have been successfully coupled with a spatial coupling scheme at the interface in a toroidal geometry. Additionally, a mapping technique is developed for transferring data between GEM's structured and XGC's unstructured meshes. Two examples of coupled simulations are presented to demonstrate the coupling scheme. The optimization of GEM for graphics processing unit is also presented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Property of neoclassical GAMs induced by pellet generated plasma perturbations in the gyrokinetic code XGC

A pellet source model has been implemented into the gyrokinetic code XGC and applied to neoclassical geodesic acoustic mode (GAM) study as the first step to investigate potentially rapid kinetic spread of pellet-born particles following their injection into plasma. Here in this study, GAM oscillations of the radial electric field, as well as of their effect on the radial particle and energy fluxes, are studied and the observed frequencies agree with theoretical expectations. It is verified that (i) GAM oscillations are driven both within the pellet-fuelled region and for locations at a larger radius than this region and not for locations radially inside the pellet-fuelled region, and that, (ii) even though the poloidal spread along the magnetic field lines is fast, the time-averaged radial plasma transport under GAM oscillations is kept to roughly the same level as the neoclassical transport in the absence of plasma turbulence.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Gyrokinetic benchmark of the electron temperature-gradient instability in the pedestal region

Transport from turbulence driven by the electron temperature-gradient (ETG) instability is likely a major source of electron heat losses through the pedestal. Due to extreme gradients and strong shaping, ETG instabilities in the pedestal are distinct from those in the core, having, for example, multiple branches (toroidal and slab) in different wavenumber ranges. Due to its importance for pedestal transport, and its rather exotic character, a rigorous multi-code benchmarking exercise is imperative. In this work, we describe such an exercise, wherein we have carried out a detailed comparison of local linear pedestal ETG simulations using three gyrokinetic codes, CGYRO, GEM, and GENE and testing different geometric parameters (such as circular, Miller, and equilibrium EFIT geometry). The resulting linear frequencies, growth rates, and eigenfunctions show very good agreement between the codes in the three types of employed geometries. A nonlinear benchmark between CGYRO and GENE is also described, exhibiting good agreement (a maximum of 20% difference in the heat fluxes computed) at two locations in the pedestal. This lays the foundation for confidently modeling ETG turbulence in the pedestal

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Toward exascale whole-device modeling of fusion devices: Porting the GENE gyrokinetic microturbulence code to GPU

GENE solves the five-dimensional gyrokinetic equations to simulate the development and evolution of plasma microturbulence in magnetic fusion devices. The plasma model used is close to first principles and computationally very expensive to solve in the relevant physical regimes. In order to use the emerging computational capabilities to gain new physics insights, several new numerical and computational developments are required. Here, we focus on the fact that it is crucial to efficiently utilize GPUs (graphics processing units) that provide the vast majority of the computational power on such systems. In this paper, we describe the various porting approaches considered and given the constraints of the GENE code and its development model, justify the decisions made, and describe the path taken in porting GENE to GPUs. We introduce a novel library called gtensor that was developed along the way to support the process. Performance results are presented for the ported code, which in a single node of the Summit supercomputer achieves a speed-up of almost 15× compared to running on central processing unit (CPU) only. Typical GPU kernels are memory-bound, achieving about 90% of peak. Our analysis shows that there is still room for improvement if we can refactor/fuse kernels to achieve higher arithmetic intensity. We also performed a weak parallel scalability study, which shows that the code runs well on a massively parallel system, but communication costs start becoming a significant bottleneck.

Germaschewski, K. (ORCID:0000000284956354)↗
DOE OSTI2022

A survey of pedestal magnetic fluctuations using gyrokinetics and a global reduced model for microtearing stability

This article presents a global reduced model for slab-like microtearing modes (MTMs) in the H-mode pedestal, which reproduces distinctive features of experimentally observed magnetic fluctuations, such as chirping and discrete frequency bands at noncontiguous mode numbers. Our model, importantly, includes the global variation of the diamagnetic frequencies, which is necessary to reproduce the experimental observations. The key insight underlying this model is that MTM instability is enabled by the alignment of a rational surface with the peak in the profile of the diamagnetic frequency. Conversely, MTMs are strongly stabilized for toroidal mode numbers for which these quantities are misaligned. This property explains the discrete fluctuation bands in several DIII-D and JET discharges, which we survey using our reduced model in conjunction with global gyrokinetic simulations. A fast yet accurate reduced model for MTMs enables rapid interpretation of magnetic fluctuation data from a wide range of experimental conditions to help assess the role of MTM in the pedestal.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Gyrokinetic analysis of inter-edge localized mode transport mechanisms in a DIII-D pedestal

In this study, gyrokinetic simulations are used to study pedestal fluctuations for DIII-D discharge 174082 using the GENE code. Nonlinear local simulations indicate that electron heat flux has contributions from electron temperature gradient-driven transport but at levels insufficient to satisfy power balance. We show that microtearing modes (MTM) and neoclassical transport are likely to account for the remaining observed energy losses in the electron and ion channels, respectively. The MTM instabilities found in the simulations are consistent with the high-frequency fluctuations identified in the magnetic fluctuation data from Mirnov coils. The fluctuation data in this discharge also exhibit a low-frequency band of fluctuations. By modifying the equilibrium profiles and plasma β, simulations produce MHD modes, which may be responsible for these observed low-frequency fluctuations. We compare several metrics involving ratios of fluctuation amplitudes and transport quantities for both MTMs and MHD modes. This analysis suggests that the available data are consistent with the simultaneous activity of both MHD modes and MTMs provided that the former is limited largely to the particle transport channel.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

A gyrokinetic simulation model for 2D equilibrium potential in the scrape-off layer of a field-reversed configuration

The equilibrium potential structure in the scrape-off layer (SOL) of the field-reversed configuration (FRC) can be affected by the penetration of edge biasing applied at the divertor ends. The primary focus of the paper is to establish a formulation that accurately captures both parallel and radial variations of the two-dimensional (2D) potential in SOL. The formulation mainly describes a quasi-neutral plasma with a logical sheath boundary. A full-f gyrokinetic ion model and a massless electron model are implemented in the GTC-X code to solve for the self-consistent equilibrium potential, given fixed radial potential profiles at the boundaries. The first essential point of this 2D model lies in its ability to couple radial and parallel dynamics stemming from resistive currents and drag force on ions. The model successfully recovers the fluid force balance and continuity equations. These collisional effects on 2D potential mainly appear through the density profile changes, modifying the potential through electron pressure gradient. This means an accurate prescription of electron density and temperature profiles is important in predicting the potential structure in the FRC SOL. The Debye sheath potential and the potential profiles applied at the boundaries can be additional factors contributing to the 2D variations in SOL. This comprehensive full-f scheme holds promise for future investigations into turbulent transport in the presence of the self-consistent 2D potential together with the non-Maxwellian distributions and open boundary conditions in the FRC SOL.

Physics↗