No Arabic abstract
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.
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.
Zero frequency zonal flow (ZFZF) excitation by trapped energetic electron driven beta-induced Alfven eigenmode (eBAE) is investigated using nonlinear gyrokinetic theory. It is found that, during the linear growth stage of eBAE, resonant energetic electrons (EEs) not only effectively drive eBAE unstable, but also contribute to the nonlinear coupling, leading to ZFZF excitation. The trapped EE contribution to ZFZF generation is dominated by EE responses to eBAE in the ideal region, and is comparable to thermal plasma contribution to Reynolds and Maxwell stresses.
Finite Larmor radius (FLR) effects on non-diffusive transport in a prototypical zonal flow with drift waves are studied in the context of a simplified chaotic transport model. The model consists of a superposition of drift waves of the linearized Hasegawa-Mima equation and a zonal shear flow perpendicular to the density gradient. High frequency FLR effects are incorporated by gyroaveraging the ExB velocity. Transport in the direction of the density gradient is negligible and we therefore focus on transport parallel to the zonal flows. A prescribed asymmetry produces strongly asymmetric non- Gaussian PDFs of particle displacements, with Levy flights in one direction but not the other. For zero Larmor radius, a transition is observed in the scaling of the second moment of particle displacements. However, FLR effects seem to eliminate this transition. The PDFs of trapping and flight events show clear evidence of algebraic scaling with decay exponents depending on the value of the Larmor radii. The shape and spatio-temporal self-similar anomalous scaling of the PDFs of particle displacements are reproduced accurately with a neutral, asymmetric effective fractional diffusion model.
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.
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.