Damping Of Zonal Flows Research Articles

Moment-based approach to the flux-tube linear gyrokinetic model

This work reports on the development and numerical implementation of the linear electromagnetic gyrokinetic (GK) model in a tokamak flux-tube geometry using a moment approach based on the expansion of the perturbed distribution function on a velocity-space Hermite–Laguerre polynomials basis. A hierarchy of equations of the expansion coefficients, referred to as the gyro-moments (GMs), is derived. We verify the numerical implementation of the GM hierarchy in the collisionless limit by performing a comparison with the continuum GK code GENE, recovering the linear properties of the ion temperature gradient, trapped electron, kinetic ballooning and microtearing modes, as well as the collisionless damping of zonal flows. An analysis of the distribution functions and ballooning eigenmode structures is performed. The present investigation reveals the ability of the GM approach to describe fine velocity-space-scale structures appearing near the trapped and passing boundary and kinetic effects associated with parallel and perpendicular particle drifts. In addition, the effects of collisions are studied using advanced collision operators, including the GK Coulomb collision operator. The main findings are that the number of GMs necessary for convergence decreases with plasma collisionality and is lower for pressure gradient-driven modes, such as in H-mode pedestal regions, compared with instabilities driven by trapped particles and magnetic gradient drifts often found in the core. The accuracy of approximations often used to model collisions (relative to the GK Coulomb operator) is studied in the case of trapped electron modes, showing differences between collision operator models that increase with collisionality and electron temperature gradient, consistent with the results of Pan et al. (Phys. Rev. E, vol. 103, 2021, L051202). Such differences are not observed in other edge microinstabilities, such as microtearing modes. The importance of a proper collision operator model is also confirmed by analysing the collisional damping of geodesic acoustic modes and zonal flows.

Nonlinear Fokker-Planck collision operator in Rosenbluth form for gyrokinetic simulations using discontinuous Galerkin method

A gyroaveraged nonlinear collision operator is formulated based on the Fokker-Planck operator in the Rosenbluth-MacDonald-Judd (RMJ) potential form and implemented for the gyrokinetic simulations with the discontinuous Galerkin scheme. The divergence structure of the original RMJ form is carefully preserved throughout the formulation to guarantee the density conservation while neglecting the finite Larmor radius effect. The B-spline finite element method is used to calculate the Rosenbluth potentials for the nonlinear collision operator. In addition to the nonlinear collision operator, linear and Dougherty collision models are also implemented to assess the benefits and drawbacks of each model. For the conservation of the parallel momentum and energy, we adopt a simple advection-diffusion model which numerically enforces the conservation of physical quantities. From bump-on-tail relaxation tests, the monotonically increasing entropy in time and conservation properties are demonstrated for the developed collision operator. Also, a few theoretical predictions for the neoclassical physics such as the neoclassical heat flux, poloidal flow and collisional damping of zonal flow are successfully reproduced by numerical simulations.

Effect of triangularity on ion-temperature-gradient-driven turbulence

The linear and nonlinear properties of ion-temperature-gradient-driven turbulence with adiabatic electrons are modeled for axisymmetric configurations for a broad range of triangularities δ, both negative and positive. Peak linear growth rates decrease with negative δ but increase and shift toward a finite radial wavenumber kx with positive δ. The growth-rate spectrum broadens as a function of kx with negative δ and significantly narrows with positive δ. The effect of triangularity on linear instability properties can be explained through its impact on magnetic polarization and curvature. Nonlinear heat flux is weakly dependent on triangularity for |δ|≤0.5, decreasing significantly with extreme δ, regardless of sign. Zonal modes play an important role in nonlinear saturation in the configurations studied, and artificially suppressing zonal modes increased nonlinear heat flux by a factor of about four for negative δ, increasing with positive δ by almost a factor of 20. Proxies for zonal-flow damping and drive suggest that zonal flows are enhanced with increasing positive δ.

Global gyrokinetic simulation with kinetic electron for collisionless damping of zonal flow in stellarators

Global gyrokinetic simulations with kinetic electrons for collisionless damping of zonal flows in LHD and W7-X stellarators show that the helical components of the equilibrium magnetic field responsible for helically trapped particles have significant impacts on zonal flow. Kinetic electrons reduce zonal flow residue and increase the frequency of low frequency oscillation (LFO). The LFO is induced by dominant helical harmonics of magnetic field strength. Furthermore, linear toroidal coupling of multiple toroidal n-harmonics barely affects the zonal flows, but can generate long wavelength toroidal harmonics with the same toroidal number as the helical magnetic field.

Development of advanced linearized gyrokinetic collision operators using a moment approach

The derivation and numerical implementation of a linearized version of the gyrokinetic (GK) Coulomb collision operator (Jorge et al., J. Plasma Phys., vol. 85, 2019, 905850604) and of the widely used linearized GK Sugama collision operator (Sugama et al., Phys. Plasmas, vol. 16, 2009, 112503) is reported. An approach based on a Hermite–Laguerre moment expansion of the perturbed gyrocentre distribution function is used, referred to as gyromoment expansion. This approach allows the considering of arbitrary perpendicular wavenumber and expressing the two linearized GK operators as a linear combination of gyromoments where the expansion coefficients are given by closed analytical expressions that depend on the perpendicular wavenumber and on the temperature and mass ratios of the colliding species. The drift-kinetic (DK) limits of the GK linearized Coulomb and Sugama operators are also obtained. Comparisons between the gyromoment approach and the DK Coulomb and GK Sugama operators in the continuum GK code GENE are reported, focusing on the ion-temperature-gradient instability and zonal flow damping, finding an excellent agreement. It is confirmed that stronger collisional damping of the zonal flow residual by the Sugama GK model compared with the GK linearized Coulomb (Pan et al., Phys. Plasmas, vol. 27, 2020, 042307) persists at higher collisionality. Finally, we show that the numerical efficiency of the gyromoment approach increases with collisionality, a desired property for boundary plasma applications.

Effects of resonant magnetic perturbations on radial electric fields in DIII-D tokamak

Gyrokinetic simulations of DIII-D tokamak equilibrium find that resonant magnetic perturbation (RMP) drives a neoclassical non-ambipolar electron particle flux, which causes a rapid change of equilibrium radial electric fields consistent with experimental observations during the suppression of the edge localized mode (ELM). The simulation results provide a support for the conjecture that RMP-induced changes of radial electric fields lead to the enhanced turbulent transport at the pedestal top during the ELM suppression (Taimourzadeh et al 2019 Nucl. Fusion 59 046005). Furthermore, gyrokinetic simulations of collisionless damping of zonal flows show that resonant responses to the RMP decrease the residual level of the zonal flows and damp the geodesic acoustic mode.

Importance of gyrokinetic exact Fokker-Planck collisions in fusion plasma turbulence.

Gyrokinetic simulations of turbulence are fundamental to understanding and predicting particle and energy loss in magnetic fusion devices. Previous works have used model collision operators with approximate field-particle terms of unknown accuracy and/or have neglected collisional finite Larmor radius effects. This work moves beyond models to demonstrate important corrections using a gyrokinetic Fokker-Planck collision operator with the exact field-particle terms, in realistic simulations of turbulence in magnetically confined fusion plasmas. The exact operator shows significant corrections for temperature-gradient-driven trapped electron mode turbulence and zonal flow damping, and for microtearing modes in a Joint European Torus pedestal under ITER-like wall conditions. Analysis of the corrections using parameter scans motivates an accurate model which closely reproduces the exact results while reducing computational demands.

Comparison of local and global gyrokinetic calculations of collisionless zonal flow damping in quasi-symmetric stellarators

The linear collisionless damping of zonal flows is calculated for quasi-symmetric stellarator equilibria in flux-tube, flux-surface, and full-volume geometry. Equilibria are studied from the quasi-helical symmetry configuration of the Helically Symmetric eXperiment (HSX), a broken symmetry configuration of HSX, and the quasi-axial symmetry geometry of the National Compact Stellarator eXperiment (NCSX). Zonal flow oscillations and long-time damping affect the zonal flow evolution, and the zonal flow residual goes to zero for small radial wavenumber. The oscillation frequency and damping rate depend on the bounce-averaged radial particle drift in accordance with theory. While each flux tube on a flux surface is unique, several different flux tubes in HSX or NCSX can reproduce the zonal flow damping from a flux-surface calculation given an adequate parallel extent. The flux-surface or flux-tube calculations can accurately reproduce the full-volume long-time residual for moderate kx, but the oscillation and damping time scales are longer in local representations, particularly for small kx approaching the system size.

Logarithmically discretized model of bounce averaged gyrokinetics and its implications on tokamak turbulence

A logarithmically discretized model, which consists of writing the system in log polar coordinates in wave-number domain and reducing the nonlinear interactions to a sum over neighboring scales that satisfy the triad conditions, is proposed for bounce averaged gyrokinetics, where the energy dependence is kept over a semi-regular grid that allows quadrature calculations in order to guarantee quasi-neutrality. The resulting model is a cheaper implementation of nonlinear multi-scale physics involving trapped electron modes, trapped ion modes, and zonal flows, which can handle anisotropy. The resulting wave-number spectrum is anisotropic at large scales, where the energy injection is clearly anisotropic, but is isotropised rapidly, leading generally towards an isotropic k−4 spectrum for spectral potential energy density for fully kinetic system and a k−5 spectrum for the system with one adiabatic species. Zonal flow damping, which is necessary for reaching a steady state in this model, plays an important role along with electron adiabaticity. Interesting dynamics akin to predator-prey evolution is observed among zonal flows and similarly large scale but radially elongated structures.

Neoclassical transport benchmark of global full-f gyrokinetic simulation in stellarator configurations

A global full-f gyrokinetic simulation code, GT5D, is extended to incorporate stellarator equilibria, provided by a widely used three-dimensional equilibrium code, VMEC [S. P. Hirshman and O. Betancourt, J. Comput. Phys. 96, 99 (1991)]. An interface code is newly developed, in which the so-called VMEC coordinates are transformed into the coordinates of GT5D. The following benchmark calculations in a typical LHD equilibrium are carried out for the numerical verification of GT5D in stellarator configurations. First, zonal flow damping tests are performed, and good agreement in the residual levels with an analytical theory is confirmed. Second, the neoclassical transport and the ambipolar radial electric field computed by a global δf neoclassical transport code, FORTEC-3D [S. Satake et al., Plasma Fusion Res. 3, S1062 (2008)], are successfully reproduced.

Observation of Oscillatory Radial Electric Field Relaxation in a Helical Plasma

Measurements of the relaxation of a zonal electrostatic potential perturbation in a nonaxisymmetric magnetically confined plasma are presented. A sudden perturbation of the plasma equilibrium is induced by the injection of a cryogenic hydrogen pellet in the TJ-II stellarator, which is observed to be followed by a damped oscillation in the electrostatic potential. The waveform of the relaxation is consistent with theoretical calculations of zonal potential relaxation in a nonaxisymmetric magnetic geometry. The turbulent transport properties of a magnetic confinement configuration are expected to depend on the features of the collisionless damping of zonal flows, of which the present Letter is the first direct observation.

Zonostrophic instability driven by discrete particle noise

The consequences of discrete particle noise for a system possessing a possibly unstable collective mode are discussed. It is argued that a zonostrophic instability (of homogeneous turbulence to the formation of zonal flows) occurs just below the threshold for linear instability. The scenario provides a new interpretation of the random forcing that is ubiquitously invoked in stochastic models such as the second-order cumulant expansion or stochastic structural instability theory; neither intrinsic turbulence nor coupling to extrinsic turbulence is required. A representative calculation of the zonostrophic neutral curve is made for a simple two-field model of toroidal ion-temperature-gradient-driven modes. To the extent that the damping of zonal flows is controlled by the ion–ion collision rate, the point of zonostrophic instability is independent of that rate.

Implications of advanced collision operators for gyrokinetic simulation

In this work, we explore both the potential improvements and pitfalls that arise when using advanced collision models in gyrokinetic simulations of plasma microinstabilities. Comparisons are made between the simple-but-standard electron Lorentz operator and specific variations of the advanced Sugama operator. The Sugama operator describes multi-species collisions including energy diffusion, momentum and energy conservation terms, and is valid for arbitrary wavelength. We report scans over collision frequency for both low and high modes, with relevance for multiscale simulations that couple ion and electron scale physics. The influence of the ion–ion collision terms—not retained in the electron Lorentz model—on the damping of zonal flows is also explored. Collision frequency scans for linear and nonlinear simulations of ion-temperature-gradient instabilities including impurity ions are presented. Finally, implications for modeling turbulence in the highly collisional edge are discussed.

Effect of resonant magnetic perturbations on microturbulence in DIII-D pedestal

Vacuum resonant magnetic perturbations (RMP) applied to otherwise axisymmetric tokamak plasmas produce in general a combination of non-resonant effects that preserve closed flux surfaces (kink response) and resonant effects that introduce magnetic islands and/or stochasticity (tearing response). The effect of the plasma kink response on the linear stability and nonlinear transport of edge turbulence is studied using the gyrokinetic toroidal code GTC for a DIII-D plasma with applied n = 2 vacuum RMP. GTC simulations use the 3D equilibrium of DIII-D discharge 158103 (Nazikian et al 2015 Phys. Rev. Lett. 114 105002), which is provided by nonlinear ideal MHD VMEC equilibrium solver in order to include the effect of the plasma kink response to the external field but to exclude island formation at rational surfaces. Analysis using the GTC simulation results reveal no increase of growth rates for the electrostatic drift wave instability and for the electromagnetic kinetic-ballooning mode in the presence of the plasma kink response to the RMP. Furthermore, nonlinear electrostatic simulations show that the effect of the 3D equilibrium on zonal flow damping is very weak and found to be insufficient to modify turbulent transport in the electrostatic turbulence.

Nonlinear damping of zonal flows

The modulatonal instability theory for the generation of large-scale (zonal) modes by drift modes has been extended to the second order including the effects of finite amplitude zonal flows, ϕq. The nonlinear (second-order) sidebands are included in the perturbative expansion to derive the nonlinear equation for the evolution of ϕq. It is shown that effects of finite ϕq reduce the growth rate of zonal flow with a possibility of oscillatory regimes at a later stage.

Direct identification of predator-prey dynamics in gyrokinetic simulations

The interaction between spontaneously formed zonal flows and small-scale turbulence in nonlinear gyrokinetic simulations is explored in a shearless closed field line geometry. It is found that when clear limit cycle oscillations prevail, the observed turbulent dynamics can be quantitatively captured by a simple Lotka-Volterra type predator-prey model. Fitting the time traces of full gyrokinetic simulations by such a reduced model allows extraction of the model coefficients. Scanning physical plasma parameters, such as collisionality and density gradient, it was observed that the effective growth rates of turbulence (i.e., the prey) remain roughly constant, in spite of the higher and varying level of primary mode linear growth rates. The effective growth rate that was extracted corresponds roughly to the zonal-flow-modified primary mode growth rate. It was also observed that the effective damping of zonal flows (i.e., the predator) in the parameter range, where clear predator-prey dynamics is observed, (i.e., near marginal stability) agrees with the collisional damping expected in these simulations. This implies that the Kelvin-Helmholtz-like instability may be negligible in this range. The results imply that when the tertiary instability plays a role, the dynamics becomes more complex than a simple Lotka-Volterra predator prey.

Collisionless kinetic-fluid simulation of zonal flows in non-circular tokamaks

Fluid simulation of linear zonal flow damping is done with a closure model based on the collisionless gyrokinetics [Sugama et al., Phys. Plasmas 14, 022502 (2007)]. Simulation results of residual zonal flow for low radial wavenumbers are compared with theoretical formulas for circular and non-circular tokamaks. The effects of the elongation and the triangularity are shown to be properly treated in the closure model. Effects of initial parallel flows on zonal flow evolution are also clarified. An appropriate choice of the initial parallel flow gives a much higher residual level than the conventional result with no initial parallel flow. Besides, the zonal flow simulations are done with the E × B nonlinearity as initial sources, which is evaluated from linear gyrokinetic microinstabilities such as ion temperature gradient modes, trapped electron modes, and electron temperature gradient modes, in order to estimate efficiency of zonal flow generation by the source instabilities.

Drift hole structure and dynamics with turbulence driven flows

The role of turbulence driven flows in describing drift hole structure and dynamics is discussed. Turbulence driven flows enter the plasma medium response and alter drift hole structures by changing the screening length of the drift hole potential. Specifically, turbulence driven flows shift the drift hole potential radially, and absorb drift hole energy via the hole-flow resonance. It is shown that the absorption shifts the phase of a momentum flux, and so enables the irreversible coupling of drift holes to turbulence driven flows. We show that drift holes and turbulence driven flows are dynamically coupled, and self-regulate each other, so that a stationary state can be achieved with non-zero turbulence driven flows. As an application, a bound on the fluctuation amplitude in the coupled system is derived. The bound is obtained by requiring that the resultant zonal flow velocity should be smaller than the critical flow velocity for the drift hole potential to be self-bound (i.e., the velocity that the screening length be positive). The result predicts |φ̃|max2∼(νd/ωci)(k∥/ky), where zonal flow damping appears as a control parameter. The implications of this result for the problem of edge-core coupling (i.e., explaining turbulence and transport in “No Man’s Land”) are discussed.

Role of stable modes in zonal flow regulated turbulence

Zonal flows are shown to regulate ion temperature gradient turbulence by enabling efficient energy transfer from the instability to a damped eigenmode in the unstable wavenumber range. The damped mode also saturates turbulence when zonal flows are not active in saturation dynamics, for example, in electron temperature gradient turbulence, but the transfer from unstable to stable mode is less efficient and requires a larger amplitude to balance the instability drive. From numerical solutions of a fluid model with a single damped eigenmode, an eigenmode decomposition of the nonlinear evolution shows that the dominant energy transfer involves the triplet correlation of the unstable mode, the zonal flow, and the stable mode at three wavenumbers satisfying k=k′+k″. In this triplet, nearly all of the energy from the instability goes to the damped mode. The very small fraction going to the zonal flow is balanced by small zonal flow damping. This combination of unstable mode, zonal flow, and stable mode minimizes the nonlinear frequency mismatch and avails itself of large coupling strengths associated with the zonal flow.

Impact of resonant magnetic perturbations on nonlinearly driven modes in drift-wave turbulence

In this work, we study the effects of resonant magnetic perturbations (RMPs) on turbulence, flows, and confinement in the framework of resistive drift wave turbulence. We extend the Hasegawa-Wakatani model to include RMP fields. The effect of the RMPs is to induce a linear coupling between the zonal electric field and the zonal density gradient, which drives the system to a state of electron radial force balance for large δBrB0. Both the vorticity flux (Reynolds stress) and particle flux are modulated. We derive an extended predator prey model which couples zonal potential and density dynamics to the evolution of turbulence intensity. This model has both turbulence drive and RMP amplitude as control parameters and predicts a novel type of transport bifurcation in the presence of RMPs. We find states that are similar to the ZF-dominated state of the standard predator-prey model, but for which the power threshold is now a function of the RMP strength. For small RMP amplitude, the energy of zonal flows decreases and the turbulence energy increases with δBrB0, corresponding to a damping of zonal flows.