Том 42, № 5 (2016)

Magnetic reconnection processes that can be excited in collisionless plasma regimes are of interest to space and astrophysics to the extent that the layers in which reconnection takes place are not rendered unrealistically small by their unfavorable dependence on relevant macroscopic distances. The equations describing new modes producing magnetic reconnection over relatively small but significant distances, unlike tearing types of mode, even when dealing with large macroscopic scale lengths, are given. The considered modes are associated with a finite electron temperature gradient and have a phase velocity in the direction of the electron diamagnetic velocity that can reverse to the opposite direction as relevant parameters are varied over a relatively wide range. The electron temperature perturbation has a primary role in the relevant theory. In particular, when referring to regimes in which the longitudinal (to the magnetic field) electron thermal conductivity is relatively large, the electron temperature perturbation becomes singular if the ratio of the transverse to the longitudinal electron thermal conductivity becomes negligible.

The nonlinear equation describing the Farley–Buneman (FB) waves in multispecies collisional plasmas is derived by employing the multiple-scale reduction analysis. It is shown that the presence of several ion species with different collisionalities and different ion masses removes the degeneracy of the nonlinear equation and generates the nonlinear terms resulting in wave steepening and wave breaking. This effect may be responsible for formation of one-dimensional coherent FB waves of a finite amplitude.

The problem of dispersion of geodesic acoustic modes is revisited with two different methods for the solution of the kinetic equation. The dispersive corrections to the mode frequency are calculated by including the m = 2 poloidal harmonics. Our obtained results agree with some earlier results but differ in various ways with other previous works. Limitations and advantages of different approaches are discussed.

A kinetic treatment of geodesic acoustic modes (GAMs), taking into account ion parallel dynamics, drift and the second poloidal harmonic effects is presented. It is shown that first and second harmonics of the ion sound modes, which have respectively positive and negative radial dispersion, can be coupled due to the geodesic and drift effects. This coupling results in the drift geodesic ion sound eigenmode with a frequency below the standard GAM continuum frequency. Such eigenmode may be able to explain the split modes observed in some experiments.

An analytic approach combining the effect of equilibrium diamagnetic flows and the finite ionsound gyroradius associated with electron−ion decoupling and kinetic Alfvén wave dispersion is derived to study resistive drift instabilities in a plasma slab. Linear numerical computations using the NIMROD code are performed with cold ions and hot electrons in a plasma slab with a doubly periodic box bounded by two perfectly conducting walls. A linearly unstable resistive drift mode is observed in computations with a growth rate that is consistent with the analytic dispersion relation. The resistive drift mode is expected to be suppressed by magnetic shear in unbounded domains, but the mode is observed in numerical computations with and without magnetic shear. In the slab model, the finite slab thickness and the perfectly conducting boundary conditions are likely to account for the lack of suppression.

The role of the bootstrap and polarization currents for the stability of neoclassical tearing modes is investigated employing both a drift kinetic and a gyrokinetic approach. The adiabatic response of the ions around the island separatrix implies, for island widths below or around the ion thermal banana width, density flattening for islands rotating at the ion diamagnetic frequency, while for islands rotating at the electron diamagnetic frequency the density is unperturbed and the only contribution to the neoclassical drive arises from electron temperature flattening. As for the polarization current, the full inclusion of finite orbit width effects in the calculation of the potential developing in a rotating island leads to a smoothing of the discontinuous derivatives exhibited by the analytic potential on which the polarization term used in the modeling is based. This leads to a reduction of the polarization-current contribution with respect to the analytic estimate, in line with other studies. Other contributions to the perpendicular ion current, related to the response of the particles around the island separatrix, are found to compete or even dominate the polarization-current term for realistic island rotation frequencies.

The influence of current density and pressure gradient profiles in the pedestal on the access to the regimes free from edge localized modes (ELMs) like quiescent H-mode in ITER is investigated. Using the simulator of MHD modes localized near plasma boundary based on the KINX code, calculations of the ELM stability were performed for the ITER plasma in scenarios 2 and 4 under variations of density and temperature profiles with the self-consistent bootstrap current in the pedestal. Low pressure gradient values at the separatrix, the same position of the density and temperature pedestals and high poloidal beta values facilitate reaching high current density in the pedestal and a potential transition into the regime with saturated large scale kink modes. New version of the localized MHD mode simulator allows one to compute the growth rates of ideal peeling-ballooning modes with different toroidal mode numbers and to determine the stability region taking into account diamagnetic stabilization. The edge stability diagrams computations and sensitivity studies of the stability limits to the value of diamagnetic frequency show that diamagnetic stabilization of the modes with high toroidal mode numbers can help to access the quiescent H-mode even with high plasma density but only with low pressure gradient values at the separatrix. The limiting pressure at the top of the pedestal increases for higher plasma density. With flat density profile the access to the quiescent H-mode is closed even with diamagnetic stabilization taken into account, while toroidal mode numbers of the most unstable peeling-ballooning mode decrease from n = 10−40 to n = 3−20.

In conditions of ideal axisymmetry, for a magnetized plasma in a generic bounded domain, necessarily toroidal, the uniform absorption of external energy (e.g., RF or any isotropic auxiliary heating) cannot give rise to net forces or torques. Experimental evidence on contemporary tokamaks shows that the near central absorption of RF heating power (ICH and ECH) and current drive in presence of MHD activity drives a bulk plasma rotation in the co-I_p direction, opposite to the initial one. Also the appearance of classical or neoclassical tearing modes provides a nonlinear magnetic braking that tends to clamp the rotation profile at the q-rational surfaces. The physical origin of the torque associated with P_RF absorption could be due the effects of asymmetry in the equilibrium configuration or in power deposition, but here we point out also an effect of the response of the so-called neoclassical offset velocity to the power dependent heat flow increment. The neoclassical toroidal viscosity due to internal magnetic kink or tearing modes tends to relax the plasma rotation to this asymptotic speed, which in absence of auxiliary heating is of the order of the ion diamagnetic velocity. It can be shown by kinetic and fluid calculations, that the absorption of auxiliary power by ions modifies this offset proportionally to the injected power thereby forcing the plasma rotation in a direction opposite to the initial, to large values. The problem is discussed in the frame of the theoretical models of neoclassical toroidal viscosity.

Bifurcation of solutions to the Grad–Shafranov-type equation for helically symmetric plasma near the threshold for tearing instability are analyzed. Quadratic and cubic nonlinearities were added to the linear dependence of the current density on the helical flux. Depending on the character of nonlinearity, two types of bifurcation can be observed, the “small” and the “large” ones. The small bifurcation is typical of cubic nonlinearity and reveals itself in the growth of the magnetic island from zero as the profile parameter increases above the instability threshold. The large bifurcation is typical of quadratic nonlinearity and causes jumplike formation of a large-scale magnetic island upon exceeding the instability threshold. As the profile parameter decreases below the instability threshold, the large-scale island continues to persist for some time (the hysteresis effect) and then suddenly disappears.

The properties of two-dimensional turbulence in a circular domain are investigated within the framework of the punctuated point-vortex model. Vortex dynamics is governed by Hamiltonian equations, and it is interrupted by instantaneous events resulting in vortex merging. Motion of about 100 point vortices is simulated using an accurate, symplectic integration method. Ensembles of like-sign vortices relax to a quasi-lattice state. Vortices with zero total vorticity tend to be randomized. Their motion still does not become fully chaotic. We observe emergence of long lived large dipoles (co-propagating pairs of vortices with opposite signs), which affect the evolution of the whole vortex ensemble. The presence of such dipoles accelerate the vortex decay rate. The decay exponent has been estimated as ξ ≃ 1.7, which is much larger than ξ ≃ 0.7, reported in previous studies of decaying turbulence. Since dipole dynamics depends on specific properties of the point vortex system, our findings suggest that a universal decay exponent in such systems does not exist.

We study the modulational instability of an intense photon beam in a relativistic pair plasma. We use the wave-kinetic description of the photon field and relativistic fluid equations for electrons and positrons. This allows us to consider the influence of the photon spectral distribution and photon recoil effects on the instability threshold and growth rates. The case of very low frequencies modulations, well below plasma frequency, is compared to that of high-frequency modulations corresponding to the plasmon decay instability.

Dusty plasma in the region of the lunar terminator is considered. It is shown that, in this region, a structure resembling a plasma sheath forms near the lunar surface. This sheath creates a potential barrier, due to which electrons over the illuminated part of the Moon are confined by electrostatic forces. The width of the sheath-like structure is on the order of the ion Debye length. In this structure, significant (about several hundred V/m) electric fields arise, which lift charged micron-size dust grains to heights of several tens of centimeters. The suggested effect may be used to explain the glow observed by the Surveyor spacecraft over the lunar terminator.