Search NASA⌕ Search

Engineering topics

Zhu, Hongxuan

Publications and source records attributed to Zhu, Hongxuan.

13 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Intrinsic Toroidal Rotation Driven by Turbulent and Neoclassical Processes in Tokamak Plasmas from Global Gyrokinetic Simulations

Gyrokinetic tokamak plasmas can exhibit intrinsic toroidal rotation driven by the residual stress. While most studies have attributed the residual stress to the parallel-momentum flux from the turbulent E × B motion, the parallel-momentum flux from the drift-orbit motion (denoted $Π^D_\parallel$) and the E × B-momentum flux from the E × B motion (denoted $Π_{E×B}$) are often neglected. Here, we use the global total-f gyrokinetic code XGC to study the residual stress in the core and the edge of a DIII-D H-mode plasma. Numerical results show that both $Π^D_\parallel$ and $Π_{E×B}$ make up a significant portion of the residual stress. In particular, $Π^D_\parallel$ in the core is higher than the collisional neoclassical level in the presence of turbulence, while in the edge it represents an outflux of countercurrent momentum even without turbulence. Using a recently developed “orbit-flux” formulation, we show that the higher-than-neoclassical-level $Π^D_\parallel$ in the core is driven by turbulence, while the outflux of countercurrent momentum from the edge is mainly due to collisional ion orbit loss. In conclusion, these results suggest that $Π^D_\parallel$ and $Π_{E×B}$ can be important for the study of intrinsic toroidal rotation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Effects of collisional ion orbit loss on neoclassical tokamak radial electric fields

We report that ion orbit loss is considered important for generating the radially inward electric field E r in a tokamak edge plasma. In particular, this effect is emphasized in diverted tokamaks with a magnetic X point. In neoclassical equilibria, Coulomb collisions can scatter ions onto loss orbits and generate a radially outward current, which in steady state is balanced by the radially inward current from viscosity. To quantitatively measure this loss-orbit current in an edge pedestal, an ion-orbit-flux diagnostic has been implemented in the axisymmetric version of the gyrokinetic particle-in-cell code XGC. As the first application of this diagnostic, a neoclassical DIII-D H-mode plasma is studied using gyrokinetic ions and adiabatic electrons. The validity of the diagnostic is demonstrated by studying the collisional relaxation of E r in the core. After this demonstration, the loss-orbit current is numerically measured in the edge pedestal in quasisteady state. In this plasma, it is found that the radial electric force on ions from E r approximately balances the ion radial pressure gradient in the edge pedestal, with the radial force from the plasma flow term being a minor component. The effect of orbit loss on E r is found to be only mild.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Effects of collisional ion orbit loss on neoclassical tokamak radial electric fields

Ion orbit loss is considered important for generating the radially inward electric field Er in a tokamak edge plasma. In particular, this effect is emphasized in diverted tokamaks with a magnetic X point. In neoclassical equilibria, Coulomb collisions can scatter ions onto loss orbits and generate a radially outward current, which in steady state is balanced by the radially inward current from viscosity. To quantitatively measure this loss-orbit current in an edge pedestal, an ion-orbit-flux diagnostic has been implemented in the axisymmetric version of the gyrokinetic particle-in-cell code XGC. As the first application of this diagnostic, a neoclassical DIII-D H-mode plasma is studied using gyrokinetic ions and adiabatic electrons. The validity of the diagnostic is demonstrated by studying the collisional relaxation of Er in the core. After this demonstration, the loss-orbit current is numerically measured in the edge pedestal in quasisteady state. In this plasma, it is found that the radial electric force on ions from Er approximately balances the ion radial pressure gradient in the edge pedestal, with the radial force from the plasma flow term being a minor component. The effect of orbit loss on Er is found to be only mild.

gyrokinetic simulations↗
DOE OSTI2021

Orbit-modulated transport and sources in time-dependent plasmas

The orbit-modulation framework for steady-state ion-orbit boundary fluxes is reformulated to apply to time-dependent plasmas. The new formulation can quantify the extent to which orbit fluxes are affected by the time dependence of the distribution function, and may be applied both to time-dependent physical scenarios and to computationally expensive simulations that are unable to reach a true steady state. As the system tends to a time-independent state, the orbit-loss fluxes approach their steady-state values, which are determined entirely by upstream transport (turbulent and collisional) and upstream sources and sinks. This result holds for the fluxes of any moments that are expressible in terms of the invariants that determine the orbits, including parallel toroidal angular momentum.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Wave-kinetic approach to zonal-flow dynamics: Recent advances

The basic physics of drift-wave turbulence and zonal flows has long been studied within the framework of the wave-kinetic theory. Recently, this framework has been reexamined from first principles, which has led to more accurate yet still tractable “improved” wave-kinetic equations. In particular, these equations reveal an important effect of the zonal-flow “curvature” (the second radial derivative of the flow velocity) on the dynamics and stability of drift waves and zonal flows. Here, we overview these recent findings and present a consolidated high-level picture of (mostly quasilinear) zonal-flow physics within reduced models of drift-wave turbulence.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Wave-kinetic approach to zonal-flow dynamics: recent advances

Basic physics of drift-wave turbulence and zonal flows has long been studied within the framework of wave-kinetic theory. Recently, this framework has been re-examined from first principles, which has led to more accurate yet still tractable "improved" wave-kinetic equations. In particular, these equations reveal an important effect of the zonal-flow "curvature" (the second radial derivative of the flow velocity) on dynamics and stability of drift waves and zonal flows. We overview these recent findings and present a consolidated high-level picture of (mostly quasilinear) zonal-flow physics within reduced models of drift-wave turbulence.

fusion plasmas↗
DOE OSTI2020

Theory of the tertiary instability and the Dimits shift within a scalar model

The Dimits shift is the shift between the threshold of the drift-wave primary instability and the actual onset of turbulent transport in a magnetized plasma. It is generally attributed to the suppression of turbulence by zonal flows, but developing a more detailed understanding calls for consideration of specific reduced models. The modified Terry–Horton system has been proposed by St-Onge ( J. Plasma Phys. , vol. 83, 2017, 905830504) as a minimal model capturing the Dimits shift. Here, we use this model to develop an analytic theory of the Dimits shift and a related theory of the tertiary instability of zonal flows. We show that tertiary modes are localized near extrema of the zonal velocity $U(x)$ , where $x$ is the radial coordinate. By approximating $U(x)$ with a parabola, we derive the tertiary-instability growth rate using two different methods and show that the tertiary instability is essentially the primary drift-wave instability modified by the local $U'' \doteq {\rm d}^2 U/{\rm d} x^2 $ . Then, depending on $U''$ , the tertiary instability can be suppressed or unleashed. The former corresponds to the case when zonal flows are strong enough to suppress turbulence (Dimits regime), while the latter corresponds to the case when zonal flows are unstable and turbulence develops. This understanding is different from the traditional paradigm that turbulence is controlled by the flow shear $| {\rm d} U / {\rm d} x |$ . Our analytic predictions are in agreement with direct numerical simulations of the modified Terry–Horton system.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Theory of the tertiary instability and the Dimits shift within a scalar model

The Dimits shift is the shift between the threshold of the drift-wave primary instability and the actual onset of turbulent transport in magnetized plasma. It is generally attributed to the suppression of turbulence by zonal flows, but developing a more detailed understanding calls for consideration of specific reduced models. The modified Terry--Horton system has been proposed by St-Onge [J. Plasma Phys. {\bf 83}, 905830504 (2017)] as a minimal model capturing the Dimits shift. Here, we use this model to develop an analytic theory of the Dimits shift and a related theory of the tertiary instability of zonal flows. We show that tertiary modes are localized near extrema of the zonal velocity $U(x)$, where $x$ is the radial coordinate. By approximating $U(x)$ with a parabola, we derive the tertiary-instability growth rate using two different methods and show that the tertiary instability is essentially the primary drift-wave instability modified by the local $U''$. Then, depending on $U''$, the tertiary instability can be suppressed or unleashed. The former corresponds to the case when zonal flows are strong enough to suppress turbulence (Dimits regime), while the latter corresponds to the case when zonal flows are unstable and turbulence develops. This understanding is different from the traditional paradigm that turbulence is controlled by the flow shear $U'$. Our analytic predictions are in agreement with direct numerical simulations of the modified Terry--Horton system.

fusion plasmas↗
DOE OSTI2020

Solitary zonal structures in subcritical drift waves: a minimum model

Solitary zonal structures have recently been identified in gyrokinetic simulations of subcritical drift-wave (DW) turbulence with background shear flows. However, the nature of these structures has not been fully understood yet. Here, we show that similar structures can be obtained within a reduced model, which complements the modified Hasegawa–Mima equation with a generic primary instability and a background shear flow. We also find that these structures can be qualitatively reproduced in the modified Hasegawa–Wakatani equation, which subsumes the reduced model as a limit. In particular, we illustrate that in both cases, the solitary zonal structures approximately satisfy the same 'equation of state', which is a local relation connecting the DW envelope with the zonal-flow velocity. Due to this generality, our reduced model can be considered as a minimum model for solitary zonal structures in subcritical DWs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Theory of the Tertiary Instability and the Dimits Shift from Reduced Drift-Wave Models

Tertiary modes in electrostatic drift-wave turbulence are localized near extrema of the zonal velocity U(x) with respect to the radial coordinate x. We argue that these modes can be described as quantum harmonic oscillators with complex frequencies, so their spectrum can be readily calculated. The corresponding growth rate $γ$ TI is derived within the modified Hasegawa-Wakatani model. We show that $γ$TI equals the primary-instability growth rate plus a term that depends on the local U''; hence, the instability threshold is shifted compared to that in homogeneous turbulence. This provides a generic explanation of the well-known yet elusive Dimits shift, which we find explicitly in the Terry-Horton limit. Linearly unstable tertiary modes either saturate due to the evolution of the zonal density or generate radially propagating structures when the shear |U'| is sufficiently weakened by viscosity. The Dimits regime ends when such structures are generated continuously.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗