No Arabic abstract
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.
The description of the local turbulent energy transfer, and the high-resolution ion distributions measured by the Magnetospheric Multiscale mission, together provide a formidable tool to explore the cross-scale connection between the fluid-scale energy cascade and plasma processes at sub-ion scales. When the small-scale energy transfer is dominated by Alfvenic, correlated velocity and magnetic field fluctuations, beams of accelerated particles are more likely observed. Here, for the first time we report observations suggesting the nonlinear wave-particle interaction as one possible mechanism for the energy dissipation in space plasmas.
We provide numerical evidence that a Kelvin-Helmholtz instability occurs in the Dirac fluid of electrons in graphene and can be detected in current experiments. This instability appears for electrons in the viscous regime passing though a micrometer-scale obstacle and affects measurements on the time scale of nanoseconds. A possible realization with a needle-shaped obstacle is proposed to produce and detect this instability by measuring the electric potential difference between contact points located before and after the obstacle. We also show that, for our setup, the Kelvin-Helmholtz instability leads to the formation of whirlpools similar to the ones reported in Bandurin et al. [Science 351, 1055 (2016)]. To perform the simulations, we develop a lattice Boltzmann method able to recover the full dissipation in a fluid of massless particles.
Observations show various jets in the solar atmosphere with significant rotational motions, which may undergo instabilities leading to heat ambient plasma. We study the Kelvin-Helmholtz (KH) instability of twisted and rotating jets caused by the velocity jumps near the jet surface. We derive a dispersion equation with appropriate boundary condition for total pressure (including centrifugal force of tube rotation), which governs the dynamics of incompressible jets. Then, we obtain analytical instability criteria of Kelvin-Helmholtz instability in various cases, which were verified by numerical solutions to the dispersion equation. We find that twisted and rotating jets are unstable to KH instability when the kinetic energy of rotation is more than the magnetic energy of the twist. Our analysis shows that the azimuthal magnetic field of 1-5 G can stabilize observed rotations in spicule/macrospicules and X-ray/EUV jets. On the other hand, non-twisted jets are always unstable to KH instability. In this case, the instability growth time is several seconds for spicule/macrospicules and few minutes (or less) for EUV/X-ray jets. We also find that standing kink and torsional Alfven waves are always unstable near the antinodes due to the jump of azimuthal velocity at the surface, while the propagating waves are generally stable. KH vortices may lead to enhanced turbulence development and heating of surrounding plasma, therefore rotating jets may provide energy for chromospheric and coronal heating.
There has been interest in recent years to assess the ability of astrophysical hydrodynamics codes to correctly model the Kelvin-Helmholtz instability. Smoothed particle hydrodynamics (SPH), in particular, has received significant attention, though there has yet to be a clear demonstration that SPH yields converged solutions that are in agreement with other methods. We have performed SPH simulations of the Kelvin-Helmholtz instability using the test problem put forward by Lecoanet et al (2016). We demonstrate that the SPH solutions converge to the reference solution in both the linear and non-linear regimes. Quantitative convergence in the strongly non-linear regime is achieved by using a physical Navier-Stokes viscosity and thermal conductivity. We conclude that standard SPH with an artificial viscosity can correctly capture the Kelvin-Helmholtz instability.
We perform simulations of the Kelvin-Helmholtz instability using smoothed particle hydrodynamics (SPH). The instability is studied both in the linear and strongly non-linear regimes. The smooth, well-posed initial conditions of Lecoanet et al. (2016) are used, along with an explicit Navier-Stokes viscosity and thermal conductivity to enforce the evolution in the non-linear regime. We demonstrate convergence to the reference solution using SPH. The evolution of the vortex structures and the degree of mixing, as measured by a passive scalar `colour field, match the reference solution. Tests with an initial density contrast produce the correct qualitative behaviour. The L2 error of the SPH calculations decreases as the resolution is increased. The primary source of error is numerical dissipation arising from artificial viscosity, and tests with reduced artificial viscosity have reduced L2 error. A high-order smoothing kernel is needed in order to resolve the initial velocity amplitude of the seeded mode and inhibit excitation of spurious modes. We find that standard SPH with an artificial viscosity has no difficulty in correctly modelling the Kelvin-Helmholtz instability and yields convergent solutions.