Fluid theory and simulations of instabilities, turbulent transport and coherent structures in partially-magnetized plasmas of discharges To cite this article: A I Smolyakov et al 2017 Plasma Phys. Control. Fusion 59 014041 View the article online for updates and enhancements. Related content Anomalous transport in high-temperature plasmas with applications to solenoidal fusion systems R.C. Davidson and N.A. Krall-Modelling electron transport in magnetized low-temperature discharge plasmas G J M Hagelaar-Physics, simulation and diagnostics of Hall effect thrusters J C Adam, J P Boeuf, N Dubuit et al.-Recent citations Nonlinear structures and anomalous transport in partially magnetized E×B plasmas Salomon Janhunen et al-Centrifugal instability in the regime of fast rotation R.

Nonlinear dynamics of the electron-cyclotron instability driven by the electron E × B current in crossed electric and magnetic field is studied. In nonlinear regime the instability proceeds by developing a large amplitude coherent wave driven by the energy input from the fundamental cyclotron resonance. Further evolution shows the formation of the long wavelength envelope akin to the modulational instability. Simultaneously, the ion density shows the development of high-k content responsible for wave focusing and sharp peaks on the periodic cnoidal wave structure. It is shown that the anomalous electron transport (along the direction of the applied electric field) is dominated by the long wavelength part of the turbulent spectrum.

The current flow in two-fluid plasma is inherently unstable if plasma components (e.g., electrons and ions) are in different collisionality regimes. A typical example is a partially magnetized E×B plasma discharge supported by the energy released from the dissipation of the current in the direction of the applied electric field (perpendicular to the magnetic field). Ions are not magnetized so they respond to the fluctuations of the electric field ballistically on the inertial time scale. In contrast, the electron current in the direction of the applied electric field is dissipatively supported either by classical collisions or anomalous processes. The instability occurs due to a positive feedback between the electron and ion current coupled by the quasi-neutrality condition. The theory of this instability is further developed taking into account the electron inertia, finite Larmor radius and nonlinear effects. It is shown that this instability results in highly nonlinear quasi-coherent structures resembling breathing mode oscillations in Hall thrusters.

No abstract

Low-frequency axial oscillations in the range of 5–50 kHz stand out as a pervasive feature observed in many types of Hall thrusters. While it is widely recognized that the ionization effects play the central role in this mode, as manifested via the large-scale oscillations of neutral and plasma density, the exact mechanism(s) of the instabilities remain unclear. To gain further insight into the physics of the breathing mode and evaluate the role of kinetic effects, a one-dimensional time-dependent full nonlinear low-frequency model describing neutral atoms, ions, and electrons is developed in full fluid formulation and compared to the hybrid model in which the ions and neutrals are kinetic. Both models are quasi-neutral and share the same electron fluid equations that include the electron diffusion, mobility across the magnetic field, and the electron energy evolution. The ionization models are also similar in both approaches. The predictions of fluid and hybrid simulations are compared for different test cases. Two main regimes are identified in both models: one with pure low-frequency behavior and the other one, where the low-frequency oscillations coexist with high-frequency oscillations in the range of 100–200 kHz, with the characteristic time scale of the ion channel fly-by time, 100–200 kHz. The other test case demonstrates the effect of a finite temperature of injected neutral atoms, which has a substantial suppression effect on the oscillation amplitude.

The structure and various components of the electron drift velocity are discussed in application to plasma discharges with the E×B drift. In high density plasmas, the contribution of the diamagnetic drift can be of the same order magnitude as the E×B drift. It is pointed out that curvature and gradient drifts associated with magnetic field inhomogeneities manifest themselves via the electron pressure anisotropy. Estimates show that the components of the diamagnetic drift related to the electron pressure anisotropy and magnetic field gradients can be important for the parameters of modern magnetrons and Hall thrusters. Similar additional terms appear in the momentum balance as mirror forces which may affect the distribution of the electrostatic potential in Hall devices.

Low-frequency ionization oscillations involving plasma and neutral density (breathing modes) are the most violent perturbations in Hall thrusters for electric propulsion. Because of its simplicity, the zero-dimensional (0D) predator–prey model of two nonlinearly coupled ordinary differential equations for plasma and neutral density has often been used for the characterization of such oscillations and scaling estimates. We investigate the properties of its continuum analog, the one-dimensional (1D) system of two nonlinearly coupled equations in partial derivatives (PDEs) for plasma and neutral density. This is a more general model, of which the standard 0D predator–prey model is a special limit case. We show that the 1D model is stable and does not show any oscillations for the boundary conditions relevant to Hall thrusters and the uniform ion velocity. We then propose a reduced 1D model based on two coupled PDEs for plasma and neutral densities that is unstable and exhibit oscillations if the ion velocity profile with the near-the-anode back-flow (toward the anode) region is used. Comparisons of the reduced model with the predictions of the full model that takes into account the self-consistent plasma response show that the main properties of the breathing mode are well captured. In particular, it is shown that the frequency of the breathing mode oscillations is weakly dependent on the final ion velocity but shows a strong correlation with the width of the ion back-flow region.

The effects of noise in particle-in-cell (PIC) and Vlasov simulations of the Buneman instability in unmagnetized plasmas are studied. It is found that, in the regime of low drift velocity, the linear stage of the instability in PIC simulations differs significantly from the theoretical predictions, whereas in the Vlasov simulations it does not. A series of highly resolved PIC simulations with increasingly large numbers of macroparticles per cell is performed using a number of different PIC codes. All the simulations predict highly similar growth rates that are several times larger than those calculated from the linear theory. As a result, we find that the true convergence of the PIC simulations in the linear regime is elusive to achieve in practice and can easily be misidentified. The discrepancy between the theoretical and the observed growth rates is attributed to the initial noise inherently present in PIC simulations, but not in Vlasov simulations, that causes particle trapping even though the fraction of trapped particles is low. We show analytically that even weak distortions of the electron velocity distribution function (such as flattening due to particle trapping) result in significant modifications of the growth rates. It is also found that the common quiet-start method for PIC simulations leads to more accurate growth rates but only if the maximum growth rate mode is perturbed initially. We demonstrate that the quiet-start method does not completely remedy the noise problem because the simulations generally exhibit inconsistencies with the linear theory.