Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userViDA: a Vlasov-DArwin solver for plasma physics at electron scales
248
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
We present a Vlasov-DArwin numerical code (ViDA) specifically designed to address plasma physics problems, where small-scale high accuracy is requested even during the non linear regime to guarantee a clean description of the plasma dynamics at fine spatial scales. The algorithm provides a low-noise description of proton and electron kinetic dynamics, by splitting in time the multi-advection Vlasov equation in phase space. Maxwell equations for the electric and magnetic fields are reorganized according to Darwin approximation to remove light waves. Several numerical tests show that ViDA successfully reproduces the propagation of linear and nonlinear waves and captures the physics of magnetic reconnection. We also discuss preliminary tests of the parallelization algorithm efficiency, performed at CINECA on the Marconi-KNL cluster. ViDA will allow to run Eulerian simulations of a non-relativistic fully-kinetic collisionless plasma and it is expected to provide relevant insights on important problems of plasma astrophysics such as, for instance, the development of the turbulent cascade at electron scales and the structure and dynamics of electron-scale magnetic reconnection, such as the electron diffusion region.
rate research
Read More
Kinetic-range turbulence in magnetized plasmas and, in particular, in the context of solar-wind turbulence has been extensively investigated over the past decades via numerical simulations. Among others, one of the widely adopted reduced plasma model is the so-called hybrid-kinetic model, where the ions are fully kinetic and the electrons are treated as a neutralizing (inertial or massless) fluid. Within the same model, different numerical methods and/or approaches to turbulence development have been employed. In the present work, we present a comparison between two-dimensional hybrid-kinetic simulations of plasma turbulence obtained with two complementary approaches spanning about two decades in wavenumber - from MHD inertial range to scales well below the ion gyroradius - with a state-of-the-art accuracy. One approach employs hybrid particle-in-cell (HPIC) simulations of freely-decaying Alfvenic turbulence, whereas the other consists of Eulerian hybrid Vlasov-Maxwell (HVM) simulations of turbulence continuously driven with partially-compressible large-scale fluctuations. Despite the completely different initialization and injection/drive at large scales, the same properties of turbulent fluctuations at $k_perprho_igtrsim1$ are observed. The system indeed self-consistently reprocesses the turbulent fluctuations while they are cascading towards smaller and smaller scales, in a way which actually depends on the plasma beta parameter. Small-scale turbulence has been found to be mainly populated by kinetic Alfven wave (KAW) fluctuations for $betageq1$, whereas KAW fluctuations are only sub-dominant for low-$beta$.
Comparisons are presented between a hybrid Vlasov-Maxwell (HVM) simulation of turbulence in a collisionless plasma and fluid reductions. These include Hall-magnetohydrodynamics (HMHD) and Landau fluid (LF) or FLR-Landau fluid (FLR-LF) models that retain pressure anisotropy and low-frequency kinetic effects such as Landau damping and, for the last model, finite Larmor radius (FLR) corrections. The problem is considered in two space dimensions, when initial conditions involve moderate-amplitude perturbations of a homogeneous equilibrium plasma subject to an out-of-plane magnetic field. LF turns out to provide an accurate description of the velocity field up to the ion Larmor radius scale, and even to smaller scales for the magnetic field. Compressibility nevertheless appears significantly larger at the sub-ion scales in the fluid models than in the HVM simulation. High frequency kinetic effects, such as cyclotron resonances, not retained by fluid descriptions, could be at the origin of this discrepancy. A significant temperature anisotropy is generated, with a bias towards the perpendicular component, the more intense fluctuations being rather spread out and located in a broad vicinity of current sheets. Non-gyrotropic pressure tensor components are measured and their fluctuations are shown to reach a significant fraction of the total pressure fluctuation, with intense regions closely correlated with current sheets.
To explain energy dissipation via turbulence in collisionless, magnetized plasmas, the existence of a dual real- and velocity-space cascade of ion-entropy fluctuations below the ion gyroradius has been proposed. Such a dual cascade, predicted by the gyrokinetic theory, has previously been observed in gyrokinetic simulations of two-dimensional, electrostatic turbulence. For the first time we show evidence for a dual phase-space cascade of ion-entropy fluctuations in a three-dimensional simulation of hybrid-kinetic, electromagnetic turbulence. Some of the scalings observed in the energy spectra are consistent with a generalized theory for the cascade that accounts for the spectral anisotropy of critically balanced, intermittent, sub-ion-Larmor-scale fluctuations. The observed velocity-space cascade is also anisotropic with respect to the magnetic-field direction, with linear phase mixing along magnetic-field lines proceeding mainly at spatial scales above the ion gyroradius and nonlinear phase mixing across magnetic-field lines proceeding at perpendicular scales below the ion gyroradius. Such phase-space anisotropy could be sought in heliospheric and magnetospheric data of solar-wind turbulence and has far-reaching implications for the dissipation of turbulence in weakly collisional astrophysical plasmas.
Vlasov solvers that operate on a phase-space grid are highly accurate but also numerically demanding. Coarse velocity space resolutions, which are unproblematic in particle-in-cell (PIC) simulations, lead to strong numerical heating or oscillations in standard continuum Vlasov methods. We present a new dual Vlasov solver which is based on an established positivity preserving advection scheme for the update of the distribution function and an energy conserving partial differential equation solver for the kinetic update of mean velocity and temperature. The solvers work together via moment fitting during which the maximum entropy part of the distribution function is replaced by the solution from the partial differential equation solver. This numerical scheme makes continuum Vlasov methods competitive with PIC methods concerning computational cost and enables us to model large scale reconnection in Earths magnetosphere with a fully kinetic continuum method. The simulation results agree well with measurements by the MMS spacecraft.
Wave properties and instabilities in a magnetized, anisotropic, collisionless, rarefied hot plasma in fluid approximation are studied, using the 16-moments set of the transport equations obtained from the Vlasov equations. These equations differ from the CGL-MHD fluid model (single fluid equations by Chew, Goldberger, and Low, 1956) by including two anisotropic heat flux evolution equations, where the fluxes invalidate the double polytropic CGL laws. We derived the general dispersion relation for linear compressible wave modes. Besides the classic incompressible fire hose modes there appear four types of compressible wave modes: two fast and slow mirror modes - strongly modified compared to the CGL model - and two thermal modes. In the presence of initial heat fluxes along the magnetic field the wave properties become different for the waves running forward and backward with respect to the magnetic field. The well known discrepancies between the results of the CGL-MHD fluid model and the kinetic theory are now removed: i) The mirror slow mode instability criterion is now the same as that in the kinetic theory. ii) Similarly, in kinetic studies there appear two kinds of fire hose instabilities - incompressible and compressible ones. These two instabilities can arise for the same plasma parameters, and the instability of the new compressible oblique fire hose modes can become dominant. The compressible fire hose instability is the result of the resonance coupling of three retrograde modes - two thermal modes and a fast mirror mode. The results can be applied to the theory of solar and stellar coronal and wind models.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
