Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “gyrokinetic simulation”

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

NSTX-U research advancing the physics of spherical tokamaks

Abstract The objectives of NSTX-U research are to reinforce the advantages of STs while addressing the challenges. To extend confinement physics of low- A , high beta plasmas to lower collisionality levels, understanding of the transport mechanisms that set confinement performance and pedestal profiles is being advanced through gyrokinetic simulations, reduced model development, and comparison to NSTX experiment, as well as improved simulation of RF heating. To develop stable non-inductive scenarios needed for steady-state operation, various performance-limiting modes of instability were studied, including MHD, tearing modes, and energetic particle instabilities. Predictive tools were developed, covering disruptions, runaway electrons, equilibrium reconstruction, and control tools. To develop power and particle handling techniques to optimize plasma exhaust in high performance scenarios, innovative lithium-based solutions are being developed to handle the very high heat flux levels that the increased heating power and compact geometry of NSTX-U will produce, and will be seen in future STs. Predictive capabilities accounting for plasma phenomena, like edge harmonic oscillations, ELMs, and blobs, are being tested and improved. In these ways, NSTX-U researchers are advancing the physics understanding of ST plasmas to maximize the benefit that will be gained from further NSTX-U experiments and to increase confidence in projections to future devices.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Cross-scale interaction between microturbulence and meso-scale reversed shear Alfvén eigenmodes in DIII-D plasmas

Abstract This paper reports global nonlinear gyrokinetic simulations that couple meso-scale reversed shear Alfvén eigenmodes (RSAEs) driven by energetic particles (EPs) and ion temperature gradient (ITG) microturbulence driven by thermal plasma, using equilibrium and profiles from DIII-D discharge #159243. In simulations focusing only on the ITG, electrostatic ITG drives a huge thermal ion heat transport, which is reduced by a factor of 10 to a level close to the experimental value in electromagnetic simulation due to finite β effect. In the simulations coupling the RSAE and ITG, ITG can scatter the resonant EP nonlinearly trapped by the RSAE and damp the zonal flows generated by the RSAE. The regulation of the RSAE by the ITG greatly reduces the initial saturation amplitude of the RSAE but increases the RSAE amplitude and associated EP transport to experimental levels in the quasi-steady state. The RSAE effects on the ITG, specifically the stronger zonal flows generated by the RSAE and the RSAE frequency modulation of the ITG-induced thermal ion heat transport, in turn, leads to a reduction of the thermal ion heat transport by more than a factor of 2 . For a stronger background ITG, the regulation of the RSAE by the ITG is stronger, while the RSAE effects on the ITG are weaker. This work highlights the importance of cross-scale coupling in the dynamics of the AE turbulence and EP transport.

Physics↗
DOE OSTI2024

Enhancing predictive capabilities in fusion burning plasmas through surrogate-based optimization in core transport solvers

Abstract This work presents the PORTALS framework (Rodriguez-Fernandez et al 2022 Nucl. Fusion 62 076036), which leverages surrogate modeling and optimization techniques to enable the prediction of core plasma profiles and performance with nonlinear gyrokinetic simulations at significantly reduced cost, with no loss of accuracy. The efficiency of PORTALS is benchmarked against standard methods, and its full potential is demonstrated on a unique, simultaneous 5-channel (electron temperature, ion temperature, electron density, impurity density and angular rotation) prediction of steady-state profiles in a DIII-D ITER Similar Shape plasma with GPU-accelerated, nonlinear CGYRO (Candy et al 2016 J. Comput. Phys. 324 73–93). This paper also provides general guidelines for accurate performance predictions in burning plasmas and the impact of transport modeling in fusion pilot plants studies.

Physics↗
DOE OSTI2024

Density pedestal prediction model for tokamak plasmas

Abstract A model for the pedestal density prediction based on neutral penetration combined with pedestal transport is presented. The model is tested against a pedestal database of JET-ILW Type I ELMy H-modes showing good agreement over a wide range of parameters both in standalone modelling (using the experimental temperature profile) and in full Europed modelling that predicts both density and temperature pedestals simultaneously. The model is further tested for ASDEX Upgrade and MAST-U Type I ELMy H-modes and both are found to agree with the same model parameters as for JET-ILW. The JET-ILW experiment where the isotope of the main ion is varied in a D/T scan at constant gas rate and constant β N is successfully modelled as long as the separatrix density ( n e,sep ) and pedestal transport coefficient ratio ( D / χ ) are varied in accordance with the experimentally observed variation of n e,sep and the isotope dependence of D / χ found in gyrokinetic simulations. The predictions are found to be sensitive to n e,sep which is why the model is combined with an n e,sep model to predict the pedestal for the STEP fusion reactor.

Physics↗
DOE OSTI2024

Transport Barriers in magnetized plasmas- general theory with dynamical constraints

Abstract A fundamental dynamical constraint—that fluctuation induced charge-weighted particle flux must vanish- can prevent instabilities from accessing the free energy in the strong gradients characteristic of Transport Barriers (TBs). Density gradients, when large enough, lead to a violation of the constraint and hence preclude unstable modes and turbulent transport. This mechanism, then, broadens the class of configurations (in magnetized plasmas) where these high confinement states can be formed and sustained. The need for velocity shear, the conventional agent for TB formation, is obviated. The most important ramifications of the constraint is to permit a charting out of the domains conducive to TB formation and hence to optimally confined fusion worthy states; the detailed investigation is conducted through new analytic methods and extensive gyrokinetic simulations.

Physics↗
DOE OSTI2024

Inter-ELM pedestal turbulence dynamics dependence on q 95 and temperature gradient

A series of dedicated experiments from DIII-D tokamak provide spatially and temporally resolved measurements of electron density and temperature, and multiscale and multichannel fluctuations over a wide range of conditions. Measurements of long wavelength density fluctuations in the type-I ELMing H-mode pedestal routinely revealed a coexistence of multiple instabilities that exhibit dramatic different dynamic behaviors as q 95 and temperature gradient are varied, apparently responsible for limiting pedestal temperature profiles. Two distinct frequency bands of density fluctuations are modulated with ELM cycle with frequency above 200kHz propagating in the electron diamagnetic direction in the lab frame (electron mode) and below 200kHz propagating in the ion diamagnetic direction (ion mode). The electron mode amplitude peaks near the electron temperature gradient region and increases with q 95 which seems to be correlated with the increased χ e at higher q 95 , similar to the characteristics expected for Micro-tearing Mode (MTM). At higher q 95 , during the inter-ELM period, the ion mode decays at later phase of the ELM cycle. Consistently, the poloidal correlation length of the ion mode is also found to reduce which suggests the possible E×B flow shear suppression of ion mode at later phase of the ELM cycle as the Er well recovers. In contrast, the electron mode grows during the ELM cycle and reaches saturation at around 50-60% of ELM period. Linear gyrokinetic simulations find the MTMs as the most unstable mode in the pedestal electron temperature gradient region. The higher q 95 and lower magnetic shear destabilize MTMs. These observations provide key insights of the underlying physics of multifield properties and rich dataset of experimental 'fingerprints' that enable new tests of the theoretical pedestal models and lead to developing a predictive model for pedestal formation on ITER and future burning plasma experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Successful prediction of tokamak transport in the L-mode regime

A long standing shortfall in the predicted L-mode edge energy transport by reduced quasi-linear models of gyrokinetic turbulent transport has been resolved. The improved model TGLF-SAT2 has higher fidelity to gyrokinetic simulations of the electron-scale contribution to the electron energy transport and the ion-scale flux surface shape dependence of energy transport. The success of TGLF-SAT2 in predicting the L-mode and Ohmic edge profiles is critical to whole pulse simulation and opens the door to prediction of the H-mode power threshold.

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

Understanding the L-H isotope effect at the DIII-D tokamak and advancements in synthetic turbulence diagnostics

Abstract It is determined that while heat flux differences between hydrogen and deuterium isotope experiments result from natural differences in carbon impurity content at DIII-D, it is not the origin of the low to high confinement mode (L-H) transition isotope effect. More specifically, a two times larger edge radial electric field in hydrogen compared to deuterium is uncovered and believed to play an important role. The origin of this radial electric field difference is determined to have two possible origins: differences in poloidal rotation and turbulent Reynolds stress in the closed field line region, and increased outer strike point temperatures and space potentials on open field lines. Experimental observations from both profile and turbulence diagnostics are supported by nonlinear gyrokinetic simulations using the code CGYRO. Simulations illustrated heat transport isotope effects in the plasma edge and shear layer resulting from differences in impurity content, electron non-adiabaticity, and main ion mass dependent E × B shear stabilization. Turbulence prediction comparisons from flux-matched CGYRO simulations to experimental measurements including electron temperature, density and velocity fluctuations are found to be in good agreement with available data. A dedicated DIII-D experiment in hydrogen was performed to seed more carbon than naturally occurring, to match deuterium experiments, and possibly reduce the L-H power threshold based on gyro-kinetic predictions. To our surprise, while ion temperature gradient (ITG) turbulence was stabilized, nodiscernible change in L-H power threshold were observed in these special hydrogen experiments. In particular, it is noticed that the edge radial electric field and Reynolds stress were observed as nearly unchanging in the presence of ITG stabilization. These experimental data have enabled a more comprehensive picture of the multitude of isotope effects at play in fusion experiments, and the important potential connection between the confined and unconfined plasma regions in regulating L-H transition dynamics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Plasma self-driven current in tokamaks with magnetic islands

Magnetic island perturbations may cause a reduction in plasma self-driven current that is needed for tokamak operation. A novel effect on tokamak self-driven current revealed by global gyrokinetic simulations is due to magnetic-island-induced 3D electric potential structures, which have the same dominant mode numbers as that of the magnetic island, whereas centered at both the inner and outer edge of the island. The non-resonant potential islands are shown to drive a current through an efficient nonlinear parallel acceleration of electrons. In large aspect ratio (large-A) tokamak devices, this new effect can result in a significant global reduction of the electron bootstrap current when the island size is sufficiently large, in addition to the local current loss across the island region due to the pressure profile flattening. It is shown that there exists a critical magnetic island width for large-A tokamaks beyond which the electron bootstrap current loss is global and increases rapidly with the island size. As such, this process may introduce a size limit for tolerable magnetic islands in large-A tokamak devices in the context of steady state operation. On the other hand, the current loss caused by magnetic islands in low-A tokamaks such as spherical tokamak (ST) NSTX/U is minor. The reduction of the axisymmetric current by magnetic islands scales with the square of island width. However, the loss of the current is mainly local to the island region, and the pace of current loss as the island size increases is substantially slower compared to large-A tokamaks. In particular, the bootstrap current reduction in STs is even smaller in the reactor-relevant high-β p regime where neoclassical tearing modes are more likely to develop.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Understanding the formation of a low-pressure pedestal in the presence of a strong internal transport barrier in DIII-D high β p plasmas

As a promising scenario for fusion reactors, the high poloidal-beta (β P ) scenario is characterized by a strong large radius internal transport barrier (ITB), which significantly enhances the overall confinement quality and the bootstrap current fraction for fully non-inductive operation. It is frequently observed that in the presence of a strong ITB, the pedestal height is lower and is accompanied by small edge localized modes (ELMs), which further improves the compatibility of a high performance core with an edge solution. A mechanism for the formation of the low pedestal is proposed in this paper. It is found that the strong ITB creates an off-axis bootstrap current to clamp the local safety factor q, and thus the magnetic shear in the outer core/pedestal region is increased. Gyrokinetic simulations with the CGYRO code show that the higher magnetic shear brings the experimental profiles into the range where the growth rate of drift-wave instabilities and thus transport is higher, and therefore a lower pedestal gradient is expected. Here, the combination of low pedestal and high magnetic shear further enhances the turbulent transport across the whole pedestal, consistent with power balance analysis. Such a positive feedback mechanism ultimately results in a lower pressure pedestal as observed in experiments. Under such a low pedestal, linear simulations with BOUT++ predict the growth rates of peeling–ballooning modes to be lower across the whole toroidal mode number spectra, and the nonlinear BOUT++ simulation exhibits lower saturated fluctuation intensity as well, consistent with the experimentally observed lower ELM size.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Suppressing trapped-electron-mode-driven turbulence via optimization of three-dimensional shaping

Abstract Turbulent transport driven by trapped electron modes (TEMs) is believed to drive significant heat and particle transport in quasihelically symmetric stellarators. Two three-dimensionally-shaped magnetic configurations with suppressed TEM-driven turbulence were generated through optimization that targeted quasihelical symmetry and the available energy of trapped electrons. Initial equilibria have flux surface shapes with a helically rotating negative triangularity (NT) and positive triangularity (PT). In gyrokinetic simulations, TEMs are suppressed in the reduced-TEM NT and PT configurations, showing that NT does not have the same beneficial turbulence properties over PT as seen in tokamaks. Heat fluxes from TEMs are also suppressed. Without temperature gradients and with a strong density gradient, the most unstable modes at low k y were consistent with toroidal universal instabilities (UIs) in the NT case and slab UIs in the PT case. Nonlinear simulations show that UIs drive substantial heat flux in both the NT and PT configurations. A moderate increase in β halves the heat flux in the NT configuration, while suppressing the heat flux in the PT geometry. Based on the present work, future optimizations aimed at reducing electrostatic drift wave-driven turbulent transport will need to consider UIs if β is sufficiently small.

Duff, J. M. (ORCID:0000000153086694)↗
DOE OSTI2025

Shafranov-shift destabilization of ballooning-type micro-instabilities

Contrary to common wisdom, we argue that the Shafranov shift is destabilizing for the ambient ballooning-type instabilities, which account for most of the transport flux in tokamak core plasmas. Higher Shafranov shift indeed reduces the magnetic drift frequency ω d , i.e. improves the bad curvature mildly, around the outboard midplane. However, this improvement is limited to a very narrow spatial region while ω d is increased over the remaining poloidal space. The eigenfunction averaged $\langle$ω d $\rangle$ is effectively enhanced due to the finite mode width of the eigenfunction and is thus destabilizing to the ballooning-type mode, as demonstrated by gyrokinetic simulations with the CGYRO code using local Miller equilibrium geometry. The predicted nonlinear flux also increases with Shafranov shift, consistent with linear simulations. The reduced transport model TGLF can capture the physics reasonably well.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Progress and innovations in the TCV tokamak research programme

Research on the Tokamak à Configuration Variable addresses a wide range of key questions relevant to ITER and future fusion power plants. Over the past two years, highly productive experimental campaigns have led to major advances across several areas: the ITER baseline scenario; pedestal properties in low-collisionality, peeling-limited conditions; and the development of high- β N , non-inductive regimes. Alternative high-confinement scenarios have likewise received significant attention, with remarkable progress in quasi-continuous exhaust operation, X-point radiator plasmas, and negative triangularity configurations. Substantial achievements were also made in the mitigation or benign termination of runaway electron beams, in elucidating fast-ion loss mechanisms, and in improving exhaust behaviour in both conventional and alternative divertor geometries. These experimental results have been strongly supported by advances in modelling and their direct application to the experiment, ranging from gyrokinetic simulations of core and pedestal turbulence to fluid-based studies of scrape-off layer and divertor physics in diverse geometries. Plasma control has taken on an increasingly important role, with model-based and data-driven approaches now closely intertwined with physics studies. This article provides a overview of these recent activities, together with a brief outlook on forthcoming upgrades and next steps.

EPFL↗
DOE OSTI2026

Quantitative analysis of zonal flow influence on turbulent plasmas driven by trapped electrons

The role of zonal flow (ZF) in the turbulence saturation of trapped electron modes (TEM) in magnetically confined plasmas is revisited. Here, in this study, we examine ZF excitation and saturation mechanisms in TEM turbulence using detailed free-energy transfer diagnostics from nonlinear gyrokinetic simulations. When the ion channel becomes subdominant, electrons shift to transferring energy to the zonal component, which results in zonal flows always playing a significant though subdominant role, in the case of the temperature gradient driven TEM. The velocity-space structure of energy transfer reveals that trapped electrons directly couple to ZF. Regarding the saturation physics, the zonal flow advection, drift-wave–drift-wave interactions, and stable modes at the pump waves are compared. When ZF saturation is weak, drift-wave–drift-wave interactions play a larger role in maintaining turbulence saturation. The findings elucidate the roles of nonzonal to zonal flow coupling, stable modes, and nonzonal wave-wave interactions, and provide a basis for improving reduced models of turbulent transport.

Song, Jiheon [Hanyang University, Seoul (Korea, Re↗
DOE OSTI2024

Reduction or Enhancement of Stellarator Turbulence by Impurities

A systematic study of the impact of impurities on the turbulent heat fluxes is presented for the stellarator Wendelstein 7-X (W7-X) and, for comparison, the Large Helical Device and ITER. By means of nonlinear multispecies gyrokinetic simulations, it is shown that impurities, depending on the sign of their density gradient, can significantly enhance or reduce turbulent ion heat losses. For the relevant scenario of turbulence reduction, an optimal impurity concentration that minimizes the ion heat diffusivity emerges as a universal feature. This result demonstrates the potential of impurities for controlling turbulence and accessing enhanced confinement regimes in fusion plasmas and, in particular, in W7-X. Published by the American Physical Society 2024

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Ion versus Electron Heating in Compressively Driven Astrophysical Gyrokinetic Turbulence

The partition of irreversible heating between ions and electrons in compressively driven (but subsonic) collisionless turbulence is investigated by means of nonlinear hybrid gyrokinetic simulations. We derive a prescription for the ion-to-electron heating ratio Q i /Q e as a function of the compressive-to-Alfvénic driving power ratio P compr /P AW , of the ratio of ion thermal pressure to magnetic pressure β i , and of the ratio of ion-to-electron background temperatures T i /T e . It is shown that Q i /Q e is an increasing function of P compr /P AW . When the compressive driving is sufficiently large, Q i /Q e approaches ≃P compr /P AW . This indicates that, in turbulence with large compressive fluctuations, the partition of heating is decided at the injection scales, rather than at kinetic scales. Analysis of phase-space spectra shows that the energy transfer from inertial-range compressive fluctuations to sub-Larmor-scale kinetic Alfvén waves is absent for both low and high β i , meaning that the compressive driving is directly connected to the ion-entropy fluctuations, which are converted into ion thermal energy. This result suggests that preferential electron heating is a very special case requiring low βi and no, or weak, compressive driving. Our heating prescription has wide-ranging applications, including to the solar wind and to hot accretion disks such as M87 and Sgr A*.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Saturation of magnetized plasma turbulence by propagating zonal flows

Strongly driven ion-scale turbulence in tokamak plasmas is shown to be regulated by a new propagating zonal flow mode, the toroidal secondary mode, which is nonlinearly supported by the turbulence. The mode grows and propagates due to the combined effects of zonal flow shearing and advection by the magnetic drift. Above a threshold in the turbulence level, small-scale toroidal secondary modes become unstable and shear apart turbulent eddies, forcing the turbulence level to remain near the threshold. This threshold condition is used to derive scaling laws for the turbulent heat flux, fluctuation spectra, and zonal flow amplitude, which are validated in nonlinear gyrokinetic simulations and explain previous experimental observations.

Nies, Richard [Princeton Univ., NJ (United States)↗