Semianalytical calculation of the zonal-flow oscillation frequency in stellarators

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.

Comments: 28 pages, 13 figures

Similar Publications

This work examines the formation of wake fields caused by ions streaming around a charged dust particle, using three-dimensional particle-in-cell (PIC) simulations with charge-neutral collisions included. The influence of an external driving electric field, which leads to a non-Maxwellian distribution of ions, is investigated in detail. The wake features formed for non-Maxwellian ions exhibit significant deviations from those observed within the model of a shifted Maxwellian distribution. Read More


We demonstrate the application of the Dynamic Mode Decomposition (DMD) for the diagnostic analysis of the nonlinear dynamics of a magnetized plasma in resistive magnetohydrodynamics. The DMD method is an ideal spatio-temporal matrix decomposition that correlates spatial features of computational or experimental data while simultaneously associating the spatial activity with periodic temporal behavior. DMD can produce low-rank, reduced order surrogate models that can be used to reconstruct the state of the system and produce high-fidelity future state predictions. Read More


The effects that energetic trapped ions have on linear resistive magnetohydrodynamic (MHD) instabilities are studied in a reduced model that captures the essential physics driving or damping the modes through variations in the magnetic shear. The drift-kinetic orbital interaction of a slowing down distribution of trapped energetic ions with a resistive MHD instability is integrated to a scalar contribution to the perturbed pressure, and entered into an asymptotic matching formalism for the resistive MHD dispersion relation. Toroidal magnetic field line curvature is included to model trapping in the particle distribution, in an otherwise cylindrical model. Read More


Microplasma generation using microwaves in an electromagnetically induced transparency (EIT)-like metasurface composed of two types of radiatively coupled cut-wire resonators with slightly different resonance frequencies is investigated. Microplasma is generated in either of the gaps of the cut-wire resonators as a result of strong enhancement of the local electric field associated with resonance and slow microwave effect. The threshold microwave power for plasma ignition is found to reach a minimum at the EIT-like transmission peak frequency, where the group index is maximized. Read More


Electron dynamics surrounding the X-line in magnetopause-type asymmetric reconnection is investigated using a two-dimensional particle-in-cell simulation. We study electron properties of three characteristic regions in the vicinity of the X-line. The fluid properties, velocity distribution functions (VDFs), and orbits are studied and cross-compared. Read More


In this paper, the structural property of complex plasma and self-diffusion coefficient of dust particles in presence of such wake potential have been investigated using Langevin dynamics simulation in the subsonic regime of ion flow. The study reveals that the self-diffusion of dust grains is strongly affected by ion flow and it changes its character in the wake potential dominant to the Debye-H\"uckel potential dominant regimes. The dependence of self-diffusion coefficient on parameters such as neutral pressure, dust size, ion flow velocity, Coulomb coupling parameter have been calculated for subsonic regime using Green-Kubo expression which is based on integrated velocity autocorrelation function (VACF). Read More


The excitation of THz radiation by the interaction of an ultra short laser pulse with the modes of a miniature corrugated plasma waveguide is considered. The axially corrugated waveguide supports the electromagnetic (EM) modes with appropriate polarization and subluminal phase velocities that can be phase matched to the ponderomotive potential associated with laser pulse, making significant THz generation possible. This process is studied via full format Particle-in-Cell (PIC) simulations that, for the first time, model the nonlinear dynamics of the plasma and the self-consistent evolution of the laser pulse in the case where the laser pulse energy is entirely depleted. Read More


Whether Turbulence-induced anomalous resistivity (AR) can facilitate a fast magnetic reconnection in collisionless plasma is a subject of active debate for decades. A particularly difficult problem in experimental and numerical simulation studies of the problem is how to distinguish the effects of AR from those originating from Hall-effect and other non-turbulent processes in the generalized Ohm's. In this paper, using particle-in-cell simulations, we present a case study of how AR produced by Buneman Instability accelerates magnetic reconnection. Read More


The Principle of Maximum Entropy, a powerful and general method for inferring the distribution function given a set of constraints, is applied to deduce the overall distribution of plasmoids (flux ropes/tubes). The analysis is undertaken for the general 3D case, with mass, total flux and (3D) velocity serving as the variables of interest, on account of their physical and observational relevance. The distribution functions for the mass, width, total flux and helicity exhibit a power-law behavior with exponents of $-4/3$, $-2$, $-3$ and $-2$ respectively for small values, whilst all of them display an exponential falloff for large values. Read More


A novel capsule target design to improve the hot-spot pressure in the high-adiabat implosion for inertial confinement fusion is proposed, where a layer of comparatively high-density material is used as a pusher between the fuel and the ablator. This design is based on our theoretical finding of the stagnation scaling laws, which indicates that the hot spot pressure can be improved by increasing the kinetic energy density $\rho_d V_{imp}^2/2$ ($\rho_d$ is the shell density when the maximum shell velocity is reached, $V_{imp}$ is the implosion velocity.) of the shell. Read More