We consider the use of “extended fluid models” as a viable alternative to computationally demanding kinetic simulations in order to manage the global large scale evolution of a collisionless plasma while accounting for the main effects that come into play when spatial micro-scales of the order of the ion inertial scale di and of the thermal ion Larmor radius .i are formed. We present an extended two-fluid model that retains finite Larmor radius (FLR) corrections to the ion pressure tensor while electron inertia terms and heat fluxes are neglected. Within this model we calculate analytic FLR plasma equilibria in the presence of a shear flow and elucidate the role of the magnetic field asymmetry. Using a Hybrid Vlasov code, we show that these analytic equilibria offer a significant improvement with respect to conventional magnetohydrodynamic shear-flow equilibria when initializing kinetic simulation

Two-dimensional simulations of the Kelvin-Helmholtz instability in an inhomogeneous compressible plasma with a density gradient show that, in a transverse magnetic field configuration, the vortex pairing process and the Rayleigh-Taylor secondary instability compete during the nonlinear evolution of the vortices. Two different regimes exist depending on the value of the density jump across the velocity shear layer. These regimes have different physical signatures that can be crucial for the interpretation of satellite data of the interaction of the solar wind with the magnetospheric plasma.

Three-dimensional simulations of the Kelvin-Helmholtz (KH) instability in a magnetic configuration reproducing typical conditions at the flank Earth's magnetosphere during northward periods show the system's ability to generate favorable conditions for magnetic reconnection to occur at mid-latitude. Once these conditions are established, magnetic reconnection proceeds spontaneously in both hemispheres generating field lines that close on Earth but are connected to the solar wind at low-latitude, allowing direct entrance of solar wind plasma into the magnetosphere. These results are consistent with recent observations of KH vortices showing the signature of reconnection events occurring well outside the equatorial plane (Bavassano et al., 2010 Ann. Geophys. 28, 893 (2010)).

The nonlinear evolution of collisionless plasmas is typically a multi-scale process where the energy is injected at large, fluid scales and dissipated at small, kinetic scales. Accurately modelling the global evolution requires to take into account the main micro-scale physical processes of interest. This is why comparison of different plasma models is today an imperative task aiming at understanding cross-scale processes in plasmas. We report here the first comparative study of the evolution of a magnetized shear flow, through a variety of different plasma models by using magnetohydrodynamic, Hall-MHD, two-fluid, hybrid kinetic and full kinetic codes. Kinetic relaxation effects are discussed to emphasize the need for kinetic equilibriums to study the dynamics of collisionless plasmas in non trivial configurations. Discrepancies between models are studied both in the linear and in the nonlinear regime of the magnetized Kelvin-Helmholtz instability, to highlight the effects of small scale processes on the nonlinear evolution of collisionless plasmas. We illustrate how the evolution of a magnetized shear flow depends on the relative orientation of the fluid vorticity with respect to the magnetic field direction during the linear evolution when kinetic effects are taken into account. Even if we found that small scale processes differ between the different models, we show that the feedback from small, kinetic scales to large, fluid scales is negligable in the nonlinear regime. This study show that the kinetic modeling validates the use of a fluid approach at large scales, which encourages the development and use of fluid codes to study the nonlinear evolution of magnetized fluid flows, even in the colisionless regime.