No Arabic abstract
In a cross-field (ExB) setup, the electron ExB flow relative to the unmagnetized ions can cause the Electron Cyclotron Drift Instability (ECDI) due to resonances of the ion-acoustic mode and the electron cyclotron harmonics. This occurs in collisionless shocks in magnetospheres and in ExB discharge devices such as Hall thrusters. ECDI induces an electron flow parallel to the background E field at a speed greatly exceeding predictions by classical collision theory. Such anomalous transport might cause unfavorable plasma flows towards the walls of ExB devices. Prediction of ECDI and anomalous transport is often thought to require a fully kinetic treatment. In this work, however, we demonstrate that a reduced variant of this instability, and more importantly, the anomalous transport, can be treated self-consistently in a collisionless two-fluid framework without any adjustable collision parameter, by treating both electron and ion species on an equal footing. We will first present linear analyses of the instability in the two-fluid 5- and 10-moment models, and compare them against the fully kinetic theory. At low temperatures, the two-fluid models predict the fastest-growing mode comparable to the kinetic results. Also, by including more moments, secondary (and possibly higher) unstable branches can be recovered. The dependence of the instability on ion-to-electron mass ratio, plasma temperature, and the background field strength is also thoroughly explored. We then carry out 5-moment simulations of the cross-field setup. The development of the instability and the anomalous transport are confirmed and in excellent agreement with theoretical predictions. The force balance properties are also studied. This work casts new insights into the nature of ECDI and the induced anomalous transport and demonstrates the potential of the two-fluid moment model in the efficient modeling of ExB plasmas.
Using fully kinetic 3D simulations, the reconnection dynamics of asymmetric current sheets are examined at the Earths magnetopause. The plasma parameters are selected to model MMS magnetopause diffusion region crossings with guide fields of 0.1, 0.4, and 1 of the reconnecting magnetosheath field. In each case, strong drift-wave fluctuations are observed in the lower-hybrid frequency range at the steep density gradient across the magnetospheric separatrix. These fluctuations give rise to cross-field electron particle transport. In addition, this turbulent mixing leads to significantly enhanced electron parallel heating in comparison to 2D simulations. We study three different methods of quantifying the anomalous dissipation produced by the drift fluctuations, based on spatial averaging, temporal averaging, and temporal averaging followed by integrating along magnetic field lines. Comparison of the different methods reveals complications in identifying and measuring the anomalous dissipation. Nevertheless, the anomalous dissipation from short wavelength drift fluctuations appears weak for each case, and the reconnection rates observed in 3D are nearly the same as in 2D models. The 3D simulations feature a number of interesting new features that are consistent with recent MMS observations, including cold beams of magnetosheath electrons that penetrate into the hotter magnetospheric inflow, the related observation of decreasing temperature in regions of increasing total density, and an effective turbulent diffusion coefficient that agrees with predictions from quasi-linear theory.
We studied the role of electron physics in 3D two-fluid 10-moment simulation of the Ganymedes magnetosphere. The model captures non-ideal physics like the Hall effect, the electron inertia, and anisotropic, non-gyrotropic pressure effects. A series of analyses were carried out: 1) The resulting magnetic field topology and electron and ion convection patterns were investigated. The magnetic fields were shown to agree reasonably well with in-situ measurements by the Galileo satellite. 2) The physics of collisionless magnetic reconnection were carefully examined in terms of the current sheet formation and decomposition of generalized Ohms law. The importance of pressure anisotropy and non-gyrotropy in supporting the reconnection electric field is confirmed. 3) We compared surface brightness morphology, represented by surface electron and ion pressure contours, with oxygen emission observed by the Hubble Space Telescope (HST). The correlation between the observed emission morphology and spatial variability in electron/ion pressure was demonstrated. Potential extension to multi-ion species in the context of Ganymede and other magnetospheric systems is also discussed.
A numerically efficient framework that takes into account the effect of the Coulomb collision operator at arbitrary collisionalities is introduced. Such model is based on the expansion of the distribution function on a Hermite-Laguerre polynomial basis, to study the effects of collisions on magnetized plasma instabilities at arbitrary mean-free path. Focusing on the drift-wave instability, we show that our framework allows retrieving established collisional and collisionless limits. At the intermediate collisionalities relevant for present and future magnetic nuclear fusion devices, deviations with respect to collision operators used in state-of-the-art turbulence simulation codes show the need for retaining the full Coulomb operator in order to obtain both the correct instability growth rate and eigenmode spectrum, which, for example, may significantly impact quantitative predictions of transport. The exponential convergence of the spectral representation that we propose makes the representation of the velocity space dependence, including the full collision operator, more efficient than standard finite difference methods.
Fluid models that approximate kinetic effects have received attention recently in the modelling of large scale plasmas such as planetary magnetospheres. In three-dimensional reconnection, both reconnection itself and current sheet instabilities need to be represented appropriately. We show that a heat flux closure based on pressure gradients enables a ten moment fluid model to capture key properties of the lower-hybrid drift instability (LHDI) within a reconnection simulation. Characteristics of the instability are examined with kinetic and fluid continuum models, and its role in the three-dimensional reconnection simulation is analysed. The saturation level of the electromagnetic LHDI is higher than expected which leads to strong kinking of the current sheet. Therefore, the magnitude of the initial perturbation has significant impact on the resulting turbulence.
After reviewing the variational approach to splitting mean flow and fluctuation kinetics in the standard Vlasov theory, the same method is applied to the drift-kinetic equation from Littlejohns theory of guiding-center motion. This process sheds a new light on drift-ordered fluid (drift-fluid) models, whose anisotropic pressure tensor is then considered in detail. In addition, current drift-fluid models are completed by the insertion of magnetization terms ensuring momentum conservation. Magnetization currents are also shown to lead to challenging aspects when drift-fluid models are coupled to Maxwells equations for the evolution of the electromagnetic field. In order to overcome these difficulties, a simplified guiding-center theory is proposed along with its possible applications to hybrid kinetic-fluid models.