A. I. Smolyakov

Zonal flow generation by parametric instability in magnetized plasmas and geostrophic fluids

Two-dimensional magnetized plasmas and geostrophic fluids exhibit a common nonlinearity due to the advection of vorticity. It is shown here that due to this nonlinearity, the propagation of small scale wave packets is accompanied by instability of a low frequency, long wavelength component. This instability is the coherent hydrodynamic generalization of the resonant type mean flow instability identified recently [P. H. Diamond, M. N. Rosenbluth, F. L. Hinton, M. Malkov, J. Fleischer, and A. Smolyakov, 17th IAEA Fusion Energy Conference, IAEA-CN-69/TH3/1, Yokohama, 1998 (to be published, International Atomic Energy Agency, Vienna)]. The mechanism discussed here, along with the resonant type, constitutes the “hydrodynamic” and “kinetic” regimes of the same process, similar to the case of plasma-beam instabilities. It is suggested that this generic mechanism is responsible for the generation of mean flow in atmospheres of rotating planets and magnetized plasmas.

Rotating nonlinear magnetic islands in a tokamak plasma

The nonlinear dynamics of rotating low m (poloidal mode number) tearing modes in a tokamak with external resonant magnetic perturbations is examined. Nonlinear evolution equations for the island width and the toroidal rotation frequency are derived within the two‐fluid magnetohydrodynamic model, taking into account the plasma rotation and neoclassical parallel viscosity. The nonlinear stability of magnetic islands interacting with a static external magnetic perturbation is considered, and the critical magnetic field for the appearance of a locked mode is determined. It is shown that the coupling of the perpendicular and longitudinal plasma flow due to the neoclassical plasma viscosity enhances the amplitude of the critical magnetic field compared to the value obtained in a slab approximation. The perpendicular plasma viscosity causes a finite phase shift between the applied external field and the magnetic island, and further increases the value of the critical magnetic field required to induce a magnetic ...

Nonlinear evolution of tearing modes in inhomogeneous plasmas

The nonlinear behaviour of tearing modes in a plasma with density and temperature gradients is reviewed. The effects of inhomogeneities can essentially modify the evolution of small scale islands from that predicted by Rutherford theory. Plasma gradient effects provide the mechanism for island excitation even in situations when the linear tearing mode stability parameter Delta is negative. The magnetic islands are sustained by the differential response of electron and ion components of a plasma in a fluctuating electric field. Such nonlinear magnetic islands are not related to the linear instability of drift-tearing modes. The nonlinear equations describing the evolution of the width and frequency of the rotating islands are derived. In the framework of one-fluid MHD, the general equation for a neighbouring equilibrium in a finite pressure plasma is considered. The dynamics of unstable m=2 and m=1 magnetic islands based on this equation is described. The quasilinear saturation of island growth in a finite pressure plasma leads to the bifurcation of the island type equilibrium into states without islands. A new evolution equation of m=1 islands is derived. For monotonic safety factor and temperature profiles this equation predicts saturation of the m=1 island growth.

Secondary instability in drift wave turbulence as a mechanism for zonal flow and avalanche formation

The article reports on recent developments in the theory of secondary instability in drift-ion temperature gradient turbulence. Specifically, the article explores secondary instability as a mechanism for zonal flow generation, transport barrier dynamics and avalanche formation. These in turn are related to the space-time statistics of the drift wave induced flux, the scaling of transport with collisionality and β, and the spatio-temporal evolution of transport barriers.

Kinetic effects in a Hall thruster discharge

Recent analytical studies and particle-in-cell simulations suggested that the electron velocity distribution function in E×B discharge of annular geometry Hall thrusters is non-Maxwellian and anisotropic. The average kinetic energy of electron motion in the direction parallel to the thruster channel walls (across the magnetic field) is several times larger than that in the direction normal to the walls. Electrons are stratified into several groups depending on their origin (e.g., plasma or channel walls) and confinement (e.g., lost on the walls or trapped in the plasma). Practical analytical formulas are derived for the plasma flux to the wall, secondary electron fluxes, plasma potential, and electron cross-field conductivity. Calculations based on these formulas fairly agree with the results of numerical simulations. The self-consistent analysis demonstrates that the elastic electron scattering in collisions with atoms and ions plays a key role in formation of the electron velocity distribution function ...

Kinetic simulation of secondary electron emission effects in Hall thrusters

The particle-in-cell code has been developed for kinetic simulations of Hall thrusters with a focus on plasma-wall interaction. It is shown that the effect of secondary electron emission on wall losses is different from predictions of previous fluid and kinetic studies. In simulations, the electron velocity distribution function is strongly anisotropic, depleted at high energy, and nonmonotonic. Secondary electrons form two beams propagating between the walls of a thruster channel in opposite radial directions. The beams produce secondary electron emission themselves depending on their energy at the moment of impact with the wall, which is defined by the electric and magnetic fields in the thruster as well as by the electron transit time between the walls. The condition for the space-charge-limited secondary electron emission depends not only on the energy of bulk plasma electrons but also on the energy of beam electrons. The contribution of the beams to the particles and energy wall losses may be much larger than that of the plasma bulk electrons. Recent experimental studies may indirectly support the results of these simulations, in particular, with respect to the electron temperature saturation and the channel width effect on the thruster discharge.

Screening of resonant magnetic perturbations by flows in tokamaks

The non-linear reduced four-field RMHD model in cylindrical geometry was extended to include plasma rotation, neoclassical poloidal viscosity and two fluid diamagnetic effects. Interaction of the static resonant magnetic perturbations (RMPs) with the rotating plasmas in tokamaks was studied. The self-consistent evolution of equilibrium electric field due to RMP penetration is taken into account in the model. It is demonstrated that in the pedestal region with steep pressure gradients, mean flows perpendicular to the magnetic field, which includes and electron diamagnetic components plays an essential role in RMP screening by plasma. Generally, the screening effect increases for lower resistivity, stronger rotation and smaller RMP amplitude. Strong screening of central islands was observed limiting RMP penetration to the narrow region near the separatrix. However, at certain plasma parameters and due to the non-linear evolution of the radial electric field produced by RMPs, the rotation can be compensated by electron diamagnetic rotation locally. In this case, RMPs can penetrate and form magnetic islands. Typical plasma parameters and RMPs spectra on DIII-D, JET and ITER were used in modelling examples presented in the paper.

Surface waves on a quantum plasma half-space

Surface modes are coupled electromagnetic/electrostatic excitations of free electrons near the vacuum-plasma interface and can be excited on a sufficiently dense plasma half-space. They propagate along the surface plane and decay in either sides of the boundary. In such dense plasma models, which are of interest in electronic signal transmission or in some astrophysical applications, the dynamics of the electrons is certainly affected by the quantum effects. Thus, the dispersion relation for the surface wave on a quantum electron plasma half-space is derived by employing the quantum hydrodynamical (QHD) and Maxwell-Poison equations. The QHD include quantum forces involving the Fermi electron temperature and the quantum Bohm potential. It is found that, at room temperature, the quantum effects are mainly relevant for the electrostatic surface plasma waves in a dense gold metallic plasma.

Generalized action invariants for drift waves-zonal flow systems

Generalized action invariants are identified for various models of drift wave turbulence in the presence of the mean shear flow. It is shown that the wave kinetic equation describing the interaction of the small scale turbulence and large scale shear flow can be naturally writen in terms of these invariants. Unlike the wave energy, which is conserved as a sum of small- and large-scale components, the generalized action invariant is shown to correspond to a quantity which is conserved for the small scale component alone. This invariant can be used to construct canonical variables leading to a different definition of the wave action (as compared to the case without shear flow). It is suggested that these new canonical action variables form a natural basis for the description of the drift wave turbulence with a mean shear flow.

Plateau regime dynamics of the relaxation of poloidal rotation in tokamak plasmas

A theory of the relaxation dynamics of the radial electric field toward its neoclassical value in the regime of subsonic poloidal rotation is presented. It is shown that the relaxation occurs via damped oscillations on time scales proportional to the ion transit time.