Search NASA⌕ Search

SEARCH · Search NASA

Results for “Gyrokinetics”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

ContourNet: Salient Local Contour Identification for Blob Detection in Fusion Plasma Simulation Data

We present ContourNet, a deep learning approach to identify salient local isocontours as blobs in large-scale 5D gyrokinetic tokamak simulation data. Blobs-regions of high turbulence that run along the edge wall down toward the diverter and can damage the tokamak-are non-well-defined features but have been empirically localized by isocontours in 2D normalized fluctuating density fields. The key of our study is to train ContourNet to follow the empirical rules to detect blobs over the time-varying simulation data. The architecture of ContourNet is a convolutional neural segmentation network: the inputs are the density field and a rasterized isocontour; the output is a set of isocontour encircling blobs. At the training stage, we feed the network with manually identified isocontours and propagated labels. At the inference stage, we extract isocontours from the segmented blob regions. Results show that our approach can achieve both high accuracy and performance, which enables scientists to understand the blob dynamics influencing the confinement of the plasma.

Blob detection↗
DOE OSTI2023

Context-aware learning of hierarchies of low-fidelity models for multi-fidelity uncertainty quantification

Multi-fidelity Monte Carlo methods leverage low-fidelity and surrogate models for variance reduction to make tractable uncertainty quantification even when numerically simulating the physical systems of interest with high-fidelity models is computationally expensive. This work proposes a context-aware multi-fidelity Monte Carlo method that optimally balances the costs of training low-fidelity models with the costs of Monte Carlo sampling. It generalizes the previously developed context-aware bi-fidelity Monte Carlo method to hierarchies of multiple models and to more general types of low-fidelity models. When training low-fidelity models, the proposed approach takes into account the context in which the learned low-fidelity models will be used, namely for variance reduction in Monte Carlo estimation, which allows it to find optimal trade-offs between training and sampling to minimize upper bounds of the mean-squared errors of the estimators for given computational budgets. This is in stark contrast to traditional surrogate modeling and model reduction techniques that construct low-fidelity models with the primary goal of approximating well the high-fidelity model outputs and typically ignore the context in which the learned models will be used in upstream tasks. Further, the proposed context-aware multi-fidelity Monte Carlo method applies to hierarchies of a wide range of types of low-fidelity models such as sparse-grid and deep-network models. Numerical experiments with the gyrokinetic simulation code Gene show speedups of up to two orders of magnitude compared to standard estimators when quantifying uncertainties in small-scale fluctuations in confined plasma in fusion reactors. This corresponds to a runtime reduction from 72 days to four hours on one node of the Lonestar6 supercomputer at the Texas Advanced Computing Center.

42 ENGINEERING↗
DOE OSTI2020

Moment preserving constrained resampling with applications to particle-in-cell methods

The Moment Preserving Constrained Resampling (MPCR) algorithm for particle resampling is introduced and applied to particle-in-cell (PIC) methods to increase simulation accuracy, reduce compute cost, and/or avoid numerical instabilities. The general algorithm partitions the system space into smaller subsets and resamples the distribution within each subset. Further, the algorithm is designed to conserve any number of particle and grid moments with a high degree of accuracy (i.e. machine accuracy). The effectiveness of MPCR is demonstrated with several numerical tests, including a use-case study in gyrokinetic fusion plasma simulations. Finally, the computational cost of MPCR is negligible compared to the cost of particle evolution in PIC methods, and the tests demonstrate that periodic particle resampling yields a significant improvement in the accuracy and stability of the results.

63 RADIATION, THERMAL, AND OTHER ENVIRON. POLLUTAN↗
DOE OSTI2020

How eigenmode self-interaction affects zonal flows and convergence of tokamak core turbulence with toroidal system size

Self-interaction is the process by which a microinstability eigenmode that is extended along the direction parallel to the magnetic field interacts non-linearly with itself. This effect is particularly significant in gyrokinetic simulations accounting for kinetic passing electron dynamics and is known to generate stationary $E\times B$ zonal flow shear layers at radial locations near low-order mode rational surfaces (Weikl et al. Phys. Plasmas, vol. 25, 2018, 072305). Here, we find that self-interaction, in fact, plays a very significant role in also generating fluctuating zonal flows, which is critical to regulating turbulent transport throughout the radial extent. Unlike the usual picture of zonal flow drive in which microinstability eigenmodes coherently amplify the flow via modulational instabilities, the self-interaction drive of zonal flows from these eigenmodes are uncorrelated with each other. It is shown that the associated shearing rate of the fluctuating zonal flows therefore reduces as more toroidal modes are resolved in the simulation. In simulations accounting for the full toroidal domain, such an increase in the density of toroidal modes corresponds to an increase in the toroidal system size, leading to a finite system size effect that is distinct from the well-known profile shearing effect.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Kinetic-ballooning-mode turbulence in low-average-magnetic-shear equilibria

In this work, kinetic-ballooning-mode (KBM) turbulence is studied via gyrokinetic flux-tube simulations in three magnetic equilibria that exhibit small average magnetic shear: the Helically Symmetric eXperiment (HSX), the helical-axis Heliotron-J and a circular tokamak geometry. For HSX, the onset of KBM being the dominant instability at low wavenumber occurs at a critical value of normalized plasma pressure $\beta ^{\rm KBM}_{\rm crit}$ that is an order of magnitude smaller than the magnetohydrodynamic (MHD) ballooning limit $\beta ^{\rm MHD}_{\rm crit}$ when a strong ion temperature gradient (ITG) is present. However, $\beta ^{\rm KBM}_{\rm crit}$ increases and approaches the MHD ballooning limit as the ITG tends to zero. For these configurations, $\beta ^{\rm KBM}_{\rm crit}$ also increases as the magnitude of the average magnetic shear increases, regardless of the sign of the normalized magnetic shear. Simulations of Heliotron-J and a circular axisymmetric geometry display behaviour similar to HSX with respect to $\beta ^{\rm KBM}_{\rm crit}$ . Despite large KBM growth rates at long wavelengths in HSX, saturation of KBM turbulence with $\beta > \beta _{\rm crit}^{\rm KBM}$ is achievable in HSX and results in lower heat transport relative to the electrostatic limit by a factor of roughly five. Nonlinear simulations also show that KBM transport dominates the dynamics when KBMs are destabilized linearly, even if KBM growth rates are subdominant to ITG growth rates.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Microstability of β ~ 1 tokamak equilibria

High-power-density tokamaks offer a potential solution to design cost-effective fusion devices. One way to achieve high power density is to operate at a high β value (the ratio of thermal to magnetic pressure), i.e. β~1 . However, a β~1 state may be unstable to various pressure- and current-driven instabilities or have unfavourable microstability properties. To explore these possibilities, we generate β~1 equilibria and investigate their stability. First, we demonstrate the generation of high-β equilibria with the computer code VMEC. We then analyse these equilibria to determine their stability against the infinite-n ideal-ballooning mode. We follow that by engaging in a detailed microstability study using the GS2 code, beginning with assessments of electrostatic ion-temperature-gradient and trapped election mode instabilities. We observe interesting behaviour for the high-β equilibria – stabilization of these modes through two distinct mechanisms – large negative local shear and reversal of electron precession drift. Finally, we perform electromagnetic gyrokinetic simulations and observe enhanced stability in the outer core of high-β equilibria and absence of kinetic ballooning modes in the negative-triangularity, high-β equilibria. The enhanced outer-core stability of high-β equilibria is different from their lower-β counterparts and offers an alternative, potentially favourable regime of tokamak operation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Quasilinear modelling of collisional trapped electron modes

The simplification of collision operators is necessary for quasilinear turbulence modelling used with integrated modelling frameworks, such as the gyrokinetic code QuaLiKiz. The treatment of collisions greatly impacts the accuracy of trapped electron mode (TEM) modelling, which is necessary to predict the electron heat flux and the balance between inward and outward particle fluxes. In particular, accurate particle flux predictions are necessary to successfully model density peaking in the tokamak core. We explored two ways of improving collisional TEM model reduction for tokamak core plasmas. First, we carried out linear GENE simulations to study the complex interplay between collisions and trapped electrons. We then used these simulations to define an effective trapped fraction to characterize the collisional TEM based on two key parameters, the local inverse-aspect ratio $\epsilon$ and the collisionality $\nu ^\ast$. One aspect missing from analytical TEM research is that the collisional frequency and the bounce-transit frequency are both velocity dependent; this effective trapped fraction takes both into account. In doing so, we determined that two parameters are not enough to model the collisional TEM, as an additional third free parameter was necessary. We determined that this model, as currently formulated, is not suitable for integrated modelling purposes. Second, we directly improved QuaLiKiz's Krook operator, which relies on two free parameters. We determined that these parameters required adjustments against higher-fidelity collisional models. In order to improve density profile predictions when paired with integrated models, we refined the Krook operator by using GENE simulations as a higher-fidelity point of comparison. We then demonstrate strong improvement of density peaking predictions of QuaLiKiz within the integrated modelling framework JETTO.

fusion plasma↗
DOE OSTI2025

Collisionless zonal-flow dynamics in quasisymmetric stellarators

The linear collisionless plasma response to a zonal-density perturbation in quasisymmetric stellarators is studied, including the geodesic-acoustic-mode oscillations and the Rosenbluth–Hinton residual flow. While the geodesic-acoustic-mode oscillations in quasiaxisymmetric configurations are similar to tokamaks, they become non-existent in quasi-helically symmetric configurations when the effective safety factor in helical-angle coordinates is small. Compared with concentric-circular tokamaks, the Rosenbluth–Hinton residual is also found to be multiplied by a geometric factor $\mathcal {C}$ that arises from the flux-surface-averaged classical polarization. Using the near-axis-expansion framework, we derive an analytic expression for $\mathcal {C}$ , which varies significantly among different configurations. These analytic results are compared with numerical simulation results from the global gyrokinetic particle-in-cell code GTC, and good agreement with the theoretical Rosenbluth–Hinton residual level is achieved when the quasisymmetry error is small enough.

Zhu, Hongxuan (ORCID:0000000198446972)↗
DOE OSTI2025

Enhanced Collisional Losses from a Magnetic Mirror Using the Lenard-Bernstein Collision Operator

Collisions are crucial in governing particle and energy transport in plasmas confined in a magnetic mirror trap. Modern gyrokinetic codes model transport in magnetic mirrors, but some utilize approximate model collision operators. This study focuses on a Pastukhov-style method of images calculation of particle and energy confinement times using a Lenard-Bernstein model collision operator. Prior work on parallel particle and energy balances used a different Fokker-Planck plasma collision operator. The method must be extended in non-trivial ways to study the Lenard-Bernstein operator. To assess the effectiveness of our approach, we compare our results with a modern finite element solver. Our findings reveal that the particle confinement time scales like a exp( a 2 ) using the Lenard-Bernstein operator, in contrast to the more accurate scaling that the Coulomb collision operator would yield a 2 exp( a 2 ), where a 2 is approximately proportional to the ambipolar potential. We propose that codes solving for collisional losses in magnetic mirrors utilizing the Lenard-Bernstein or Dougherty collision operator scale their collision frequency of any electrostatically confined species. This study illuminates the collision operator’s intricate role in the Pastukhov-style method of images calculation of collisional confinement.

fusion plasma↗
DOE OSTI2020

Dependence of the impurity transport on the dominant turbulent regime in ELM-y H-mode discharges on the DIII-D tokamak

Laser blow-off injections of aluminum and tungsten have been performed on the DIII-D tokamak to investigate the variation of impurity transport in a set of dedicated ion and electron heating scans with a fixed value of the external torque. The particle transport is quantified via the Bayesian inference method, which, constrained by a combination of a charge exchange recombination spectroscopy, soft X-ray measurements, and VUV spectroscopy provides a detailed uncertainty quantification of transport coefficients. Contrasting discharge phases with a dominant electron and ion heating reveal 3-fold drop in the impurity confinement time and order of magnitude increase in midradius impurity diffusion, when additional electron heating is applied. Further, the calculated stationary aluminum density profiles reverse from peaked in electron heated to hollow in the ion heated case, following a similar trend as electron and carbon density. Comparable values of a core diffusion have been observed for W and Al ions, while differences in the propagation dynamics of these impurities are attributed to pedestal and edge transport. Modeling of the core transport with non-linear gyrokinetics code CGYRO, significantly underpredicts the magnitude of the variation in Al transport. Diffusion increases 3-times steeper with additional electron heat flux and 10-times lower diffusion is observed in ion heated case than predicted by the modeling. Here, the CGYRO model quantitatively matches the increase of the Al diffusion when approaching the linear threshold for the transition from the ion temperature gradient (ITG) to trapped electron mode (TEM).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Constructing a new predictive scaling formula for ITER's divertor heat-load width informed by a simulation-anchored machine learning

Understanding and predicting divertor heat-load width λq is a critically important problem for an easier and more robust operation of ITER with high fusion gain. Previous predictive simulation data for λ q using the extreme-scale edge gyrokinetic code XGC1 [S. Ku et al., Phys. Plasmas 25, 056107 (2018)] in the electrostatic limit under attached divertor plasma conditions in three major US tokamaks [C. S. Chang et al., Nucl. Fusion 57, 116023 (2017)] reproduced the Eich and Goldston attached-divertor formula results [formula #14 in T. Eich et al., Nucl. Fusion 53, 093031 (2013) and R. J. Goldston, Nucl. Fusion 52, 013009 (2012)] and furthermore predicted over six times wider λ q than the maximal Eich and Goldston formula predictions on a full-power (Q = 10) scenario ITER plasma. After adding data from further predictive simulations on a highest current JET and highest-current Alcator C-Mod, a machine learning program is used to identify a new scaling formula for λ q as a simple modification to the Eich formula #14, which reproduces the Eich scaling formula for the present tokamaks and which embraces the wide λ q XGC for the full-current Q = 10 ITER plasma. Additionally, the new formula is then successfully tested on three more ITER plasmas: two corresponding to long burning scenarios with Q = 5 and one at low plasma current to be explored in the initial phases of ITER operation. The new physics that gives rise to the wider λ q XGC is identified to be the weakly collisional, trapped-electron-mode turbulence across the magnetic separatrix, which is known to be an efficient transporter of the electron heat and mass. Electromagnetic turbulence and high-collisionality effects on the new formula are the next study topics for XGC1.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Saturation physics of threshold heat-flux reduction

The saturation physics of ion-temperature-gradient-driven turbulence is examined in relation to the temperature-gradient variation of the heat flux, which can exhibit an upshift of the critical gradient for significant flux relative to the linear instability threshold. Gyrokinetic measurements of saturation properties and spectral energy transfer, which will be defined in Sec. II, are presented, indicating that the physics of saturation is fundamentally unchanged on either side of the upshifted gradient. To analyze heat transport below and above the upshifted critical gradient, a fluid model for toroidal ion-temperature-gradient turbulence is modified to include the kinetic instability threshold. The model and the heat flux are rendered in the eigenmode decomposition to track the dominant mode-coupling channel of zonal-flow-catalyzed transfer to a conjugate stable mode. Given linear and nonlinear symmetries, the stable mode level and the cross-correlation of the unstable and stable mode amplitudes are related to the unstable mode level via linear physics. Here, the heat flux can then be written in terms of the unstable-mode level, which through a nonlinear balance depends on the eigenmode-dependent coupling coefficients and the triplet correlation time of the dominant coupled modes. Resonance in these quantities leads to suppressed heat flux above the linear threshold, with a nonlinear upshift of the critical gradient set by the resonance broadening of a finite perpendicular wavenumber and collisionality.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Edge turbulence measurements in L-mode and I-mode at ASDEX Upgrade

The I-mode confinement regime is promising for future reactor operation due to high energy confinement without high particle confinement. However, the role of edge turbulence in creating I-mode's beneficial transport properties is still unknown. New measurements of edge turbulence (ρpol=0.9−1.0) in L-modes and I-modes at low and high densities at ASDEX Upgrade are presented in this paper. A high radial resolution correlation electron cyclotron emission radiometer measures the broadband turbulence throughout the L-mode and I-mode edge and pedestal. The weakly coherent mode (WCM) is measured in both L-mode and I-mode near the last closed flux surface with Te fluctuation levels of 2.3%–4.2%, with a frequency shift between the two phases related to a deeper Er well in I-mode. An neTe phase diagnostic captures a change of the WCM neTe phase between L-mode and I-mode from −171° to −143°. The thermal He beam diagnostic measures a WCM wavenumber range of −0.5 to −1.0 cm−1. A low-frequency edge oscillation (LFEO) appears in the I-mode phase of these discharges and displays coupling to the WCM, but the LFEO does not appear in the L-mode phase. Linear gyrokinetic simulations of the outer core and pedestal top turbulence indicate that while the dominant turbulent modes in the outer core are ion directed and electrostatic, the turbulence becomes increasingly electron directed and electromagnetic with increasing radius. Collisionality is not found to impact characteristics of the L-mode and I-mode edge turbulence with respect to the presence of the WCM; however, the quality of global confinement decreases with collisionality.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Verification of a fully kinetic ion model for electromagnetic simulations of high-frequency waves in toroidal geometry

For the research of high-frequency electromagnetic waves in tokamaks, an electromagnetic simulation model, in which the ion dynamics is described by a six-dimensional Vlasov equation and the electron dynamics is described by a drift kinetic equation, is formulated and implemented in the global gyrokinetic toroidal code (GTC). Analytic dispersion relations are derived in reduced systems and compared with various theories to verify the model. Linear simulations of a generalized ion Bernstein wave and ion cyclotron emission are verified by comparing the GTC simulation results with analytic dispersion relation theory and magnetoacoustic cyclotron instability theory, respectively, in cylindrical geometry.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Electron temperature gradient driven transport model for tokamak plasmas

A new model for electron temperature gradient (ETG) modes is developed as a component of the multi-mode anomalous transport module [Rafiq et al., Phys Plasmas 20, 032506 (2013)] to predict a time-dependent electron temperature profile in conventional and low aspect ratio tokamaks. This model is based on two-fluid equations that govern the dynamics of low-frequency short- and long-wavelength electromagnetic toroidal ETG driven drift modes. A low collisionality NSTX discharge is used to scan the plasma parameter dependence on the ETG real frequency, growth rate, and electron thermal diffusivity. Electron thermal transport is discovered in the deep core region where modes are more electromagnetic in nature. Several previously reported gyrokinetic trends are reproduced, including the dependencies of density gradients, magnetic shear, β and gradient of β (β'), collisionality, safety factor, and toroidicity, where β is the ratio of the plasma pressure to the magnetic pressure. Furthermore, the electron heat diffusivity associated with the ETG mode is discovered to be on a scale consistent with the experimental diffusivity determined by power balance analysis.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Experimental observation of a field-aligned ion beam produced by magnetic reconnection of two flux ropes

An ion beam field-aligned to the background guide field ($B$ 0 = $330$ G) was observed in a reconnection experiment on the Large Plasma Device and, to the authors' knowledge, is the first experimental observation of its kind. Two kink-unstable flux ropes ( L = 11 m, d = 7.6 cm) were made to collide, which allows magnetic reconnection to occur. Sub-Alfvénic ion beams with energies of up to 15 eV were then observed from measurements of the local ion energy distribution function. The beam ions do not appear to be heated. Further, they were correlated with the collision of the ropes and appear to be energized by magnetic reconnection. The results and interpretation of the measurements are supported by three-dimensional gyrokinetic particle simulations of the merging flux ropes and electric field measurements from previous experiments [W. Gekelman et al. , Astrophys. J. 853, 33 (2018)]. The mechanism behind the acceleration appears to be non-local.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Nonlinear dynamics of small-scale Alfvén waves

Herein we study the nonlinear evolution of very oblique small-scale Alfvén waves with k ⊥ d i ≳1. At these scales, the waves become significantly compressive, unlike in magnetohydrodynamics, due to the Hall term in the equations. We demonstrate that when frequencies are small compared to the ion gyrofrequency and amplitudes are small compared to unity, no new nonlinear interaction appears due to the Hall term alone at the lowest non-trivial order, even when k ⊥ d i ~1. However, at the second non-trivial order, we discover that the Hall physics leads to a slow but resonant nonlinear interaction between co-propagating Alfvén waves, an inherently three-dimensional effect. Including the effects of finite temperature, finite frequency, and electron inertia, the two-fluid Alfvén wave also becomes dispersive once one or more of k ⊥ ρ s , k ⊥ d e , or k ∥ d i becomes significant: for oblique waves at low β as studied here, this can be at a much smaller scale than d i . We show that the timescale for one-dimensional steepening of two-fluid Alfvén waves is only significant at these smaller dispersive scales, and also derive an expression for the amplitude of driven harmonics of a primary wave. Importantly, both new effects are absent in gyrokinetics and other commonly used reduced two-fluid models. Our calculations have relevance for the interpretation of laboratory Alfvén wave experiments, as well as shedding light on the physics of turbulence in the solar corona and inner solar wind, where the dominant nonlinear interaction between counter-propagating waves is suppressed, allowing these new effects to become important.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Ballooning theory for micro-tearing mode in tokamak

This paper aims to investigate the impact of magnetic drift on the linear micro-tearing mode by using a kinetic approach to derive a reduced two-field eigen system in real space. Here, since the magnetic drift in real space has derivatives, it is more efficient to solve the mode equations in a Fourier-ballooning representation using the two-dimensional (2D) ballooning transform. The lowest-order eigen system in the Fourier-ballooning representation consists of two integral equations, which are numerically solved using the finite difference method for both eigenvalues and wave functions. The main results will be presented through graphical eigenvalue scans for each parameter. Furthermore, we present a graphical comparison between the predictions of the ballooning theory and GENE gyrokinetic code simulation in the pedestal region.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗