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.

Abstract. Due to their capability to reduce turbulent transport in magnetized plasmas, understanding the dynamics of zonal flows is an important problem in the fusion programme. Since the pioneering work by Rosenbluth and Hinton in axisymmetric tokamaks, it is known that studying the linear and collisionless relaxation of zonal flow perturbations gives valuable information and physical insight. Recently, the problem has been investigated in stellarators and it has been found that in these devices the relaxation process exhibits a characteristic feature: a damped oscillation. The frequency of this oscillation might be a relevant parameter in the regulation of turbulent transport, and therefore its efficient and accurate calculation is important. Although an analytical expression can be derived for the frequency, its numerical evaluation is not simple and has not been exploited systematically so far. Here, a numerical method for its evaluation is considered, and the results are compared with those obtained by calculating the frequency from gyrokinetic simulations. This "semianalytical" approach for the determination of the zonal-flow frequency reveals accurate and faster than the one based on gyrokinetic simulations.

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.

Abstract. The results of global linear gyrokinetic simulations of residual flows carried out with the code EUTERPE in the TJ-II three dimensional geometry are reported. The linear response of the plasma to potential perturbations homogeneous in a magnetic surface shows several oscillation frequencies: a GAMlike frequency, in qualitative agreement with the formula given by Sugama et al [1], and a much lower frequency oscillation in agreement with the predictions of Mishchenko et al [2,3] for stellarators. The dependency of both oscillations with ion and electron temperatures and the magnetic configuration is studied. The low frequency oscillations are in the frequency range supporting the longrange correlations between potential signals experimentally observed in TJ-II.

The low frequency oscillatory relaxation of zonal potential perturbations is studied numerically in the TJ-II stellarator (where it was experimentally detected for the first time). It is studied in full global gyrokinetic simulations of multi-species plasmas. The oscillation frequency obtained is compared with predictions based on single-species simulations using simplified analytical relations. It is shown that the frequency of this oscillation for a multi-species plasma can be accurately obtained from single-species calculations using extrapolation formulas. The damping of the oscillation and the influence of the different inter-species collisions is studied in detail. It is concluded that taking into account multiple kinetic ions and electrons with impurity concentrations realistic for TJ-II plasmas allows to account for the values of frequency and damping rate in zonal flows relaxations observed experimentally.