Two-fluid and kinetic transport physics of Kelvin–Helmholtz instabilities in nonuniform low-beta plasmas

• Published in: Physics of plasmas 2020-10, Vol.27 (10)
• Main Authors: Vogman, G. V., Hammer, J. H., Shumlak, U., Farmer, W. A.
• Format: Article
• Language: English
• Subjects: 70 PLASMA PHYSICS AND FUSION TECHNOLOGY; Collisionless plasmas; Computational fluid dynamics; Density gradients; Diamagnetism; Diffusion layers; Distribution functions; Drift; Electric fields; Fluid flow; Kelvin-Helmholtz instability; Magnetic properties; Magnetization; Magnetohydrodynamics; MATHEMATICS AND COMPUTING; Nonuniform plasmas; Parameter modification; Plasma physics; Shear layers; Simulation; Space charge; Stability analysis; Transport properties
• ISSN: 1070-664X; 1089-7674
• DOI: 10.1063/5.0014489

Hall-magnetohydrodynamic (Hall-MHD) theory, two-fluid simulations, and kinetic simulations are used to investigate the cross-field transport properties of Kelvin–Helmholtz instabilities in nonuniform low-beta collisionless plasmas. Hall-MHD analysis shows how the linear properties of the instability are modified by density gradients and magnetization. High-order accurate two-fluid and kinetic simulations, with complete dynamics of finite-mass electrons and ions, are applied to a suite of parameter cases to systematically assess the effects of diamagnetic drift, magnetization, charge separation, and finite Larmor motion. Initialization of exact two-species kinetic equilibria facilitates the study of isolated physical effects and enables detailed cross-comparisons between two-fluid and kinetic simulations, including for cases where ion gyroradii are comparable to gradient scale lengths. For nonuniform plasmas with significant space charge, the results of two-fluid and kinetic simulations are found to disagree with Hall-MHD predictions. Kelvin–Helmholtz instability growth rates, per unit shear, are shown to be smaller when ion diamagnetic drift and E × B drift are parallel and larger when the two drifts are antiparallel. The effect is attributed to polarization drift in the shear layer, which leads to redistribution of charge, alters the electric field that drives plasma advection, and consequently modifies growth rates. Instability-induced mass transport for different parameters is characterized in terms of the flux across the shear layer and a simplified diffusion model. Distribution functions from kinetic simulations are shown to deviate substantially from Maxwellian reconstructions, indicating the importance of kinetic physics during the nonlinear phase of the instability.