No Arabic abstract
In the linear collisionless limit, a zonal potential perturbation in a toroidal plasma relaxes, in general, to a non-zero residual value. Expressions for the residual value in tokamak and stellarator geometries, and for arbitrary wavelengths, are derived. These expressions involve averages over the lowest order particle trajectories, that typically cannot be evaluated analytically. In this work, an efficient numerical method for the evaluation of such expressions is reported. It is shown that this method is faster than direct gyrokinetic simulations performed with the GENE and EUTERPE codes. Calculations of the residual value in stellarators are provided for much shorter wavelengths than previously available in the literature. Electrons must be treated kinetically in stellarators because, unlike in tokamaks, kinetic electrons modify the residual value even at long wavelengths. This effect, that had already been predicted theoretically, is confirmed by gyrokinetic simulations.
In stellarators, zonal flow activity depends sensitively on geometry of the three dimensional magnetic field, via an interplay of mechanisms that is not fully understood. In this work, we investigate this by studying three magnetic configurations of the Wendelstein 7-X stellarator. We find that variation in linear zonal flow damping is accompanied by variation in nonlinear drive, and identify key geometric features that control these effects. Understanding the resulting balance is important for the development of reduced models of turbulent transport.
The collisionless axisymmetric zonal flow residual calculation for a tokamak plasma is generalized to include electromagnetic perturbations. We formulate and solve the complete initial value zonal flow problem by retaining the fully self-consistent axisymmetric spatial perturbations in the electric and magnetic fields. Simple expressions for the electrostatic, shear and compressional magnetic residual responses are derived that provide a fully electromagnetic test of the zonal flow residual in gyrokinetic codes. Unlike the electrostatic potential, the parallel vector potential and the parallel magnetic field perturbations need not relax to flux functions for all possible initial conditions.
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 $k_x$, but the oscillation and damping time scales are longer in local representations, particularly for small $k_x$ approaching the system size.
In tokamak plasmas, the interaction among the micro-turbulence, zonal flows (ZFs) and energetic particles (EPs) can affect the turbulence saturation level and the consequent confinement quality and thus, is important for future burning plasmas. In this work, the EP anisotropy effects on the ZF residual level are studied by using anisotropic EP distributions with dependence on pitch. Significant effects on the long wavelength ZFs have been found when small to moderate width around the dominant pitch in the EP distribution function is assumed. In addition, it is found that ZF residual level is enhanced by barely passing/trapped and/or deeply trapped EPs, but it is suppressed by well passing and/or intermediate trapped EPs. Numerical calculation shows that for ASDEX Upgrade plasmas, typical EP distribution functions can bring in -3%~+5.5% mitigation/enhancement in ZF residual level, depending on the EP distribution functions.
We present a theory of the nonlinear growth of zonal flows in magnetized plasma turbulence, by the mechanism of secondary instability. The theory is derived for general magnetic geometry, and is thus applicable to both tokamaks and stellarators. The predicted growth rate is shown to compare favorably with nonlinear gyrokinetic simulations, with the error scaling as expected with the small parameter of the theory.