Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userBridging hybrid- and full-kinetic models with Landau-fluid electrons: I. 2D magnetic reconnection
374
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
Magnetic reconnection (MR) plays a fundamental role in plasma dynamics under many different conditions, from space and astrophysical environments to laboratory devices. High-resolution in-situ measurements from space missions allow to study naturally occurring MR processes in great detail. Alongside direct measurements, numerical simulations play a key role in investigating the fundamental physics underlying MR. The choice of an adequate plasma model to be employed in numerical simulations, while also compromising with their computational cost, is crucial to efficiently address the problem. We consider a new plasma model that includes a refined electron response within the hybrid-kinetic framework (kinetic ions, fluid electrons). The extent to which this new model can reproduce a full-kinetic description of 2D MR, with particular focus on its robustness during the non-linear stage, is evaluated. We perform 2D simulations of MR with moderate guide field by means of three different plasma models: a hybrid-Vlasov-Maxwell model with isotropic, isothermal electrons, a hybrid-Vlasov-Landau-fluid (HVLF) model where an anisotropic electron fluid is equipped with a Landau-fluid closure, and a full-kinetic one. When compared to the full-kinetic case, the HVLF model effectively reproduces the main features of MR, as well as several aspects of the associated electron micro-physics and its feedback onto proton dynamics. This includes the global evolution of MR and the local physics occurring within the so-called electron-diffusion region, as well as the evolution of species pressure anisotropy. In particular, anisotropy driven instabilities (such as firehose, mirror, and cyclotron instabilities) play a relevant role in regulating electrons anisotropy during the non-linear stage of MR. As expected, the HVLF model captures all these features, except for the electron-cyclotron instability.
rate research
Read More
The space-filter approach has proved a fundamental tool in studying turbulence in neutral fluids, providing the ability to analyze scale-to-scale energy transfer in configuration space. It is well known that turbulence in plasma presents challenges different from neutral fluids, especially when the scale of interests include kinetic effects. The space-filter approach is still largely unexplored for kinetic plasma. Here we derive the space-filtered (or, equivalently coarse-grained) equations in configuration space for a quasi-neutral hybrid-kinetic plasma model, in which ions are fully kinetic and electrons are a neutralizing fluid. Different models and closures for the electron fluid are considered, including finite electron-inertia effects and full electrons pressure-tensor dynamics. Implications for the cascade of turbulent fluctuations in real space depending on different approximations are discussed.
As modeling of collisionless magnetic reconnection in most space plasmas with realistic parameters is beyond the capability of todays simulations, due to the separation between global and kinetic length scales, it is important to establish scaling relations in model problems so as to extrapolate to realistic scales. Recently, large scale particle-in-cell (PIC) simulations of island coalescence have shown that the time averaged reconnection rate decreases with system size, while fluid systems at such large scales in the Hall regime have not been studied. Here we perform the complementary resistive MHD, Hall MHD and two fluid simulations using a ten-moment model with the same geometry. In contrast to the standard Harris sheet reconnection problem, Hall MHD is insufficient to capture the physics of the reconnection region. Additionally, motivated by the results of a recent set of hybrid simulations which show the importance of ion kinetics in this geometry, we evaluate the efficacy of the ten-moment model in reproducing such results.
In Part 2 of our guide to collisionless fluid models, we concentrate on Landau fluid closures. These closures were pioneered by Hammett and Perkins and allow for the rigorous incorporation of collisionless Landau damping into a fluid framework. It is Landau damping that sharply separates traditional fluid models and collisionless kinetic theory, and is the main reason why the usual fluid models do not converge to the kinetic description, even in the long-wavelength low-frequency limit. We start with a brief introduction to kinetic theory, where we discuss in detail the plasma dispersion function $Z(zeta)$, and the associated plasma response function $R(zeta)=1+zeta Z(zeta)=-Z(zeta)/2$. We then consider a 1D (electrostatic) geometry and make a significant effort to map all possible Landau fluid closures that can be constructed at the 4th-order moment level. These closures for parallel moments have general validity from the largest astrophysical scales down to the Debye length, and we verify their validity by considering examples of the (proton and electron) Landau damping of the ion-acoustic mode, and the electron Landau damping of the Langmuir mode. We proceed by considering 1D closures at higher-order moments than the 4th-order, and as was concluded in Part 1, this is not possible without Landau fluid closures. We show that it is possible to reproduce linear Landau damping in the fluid framework to any desired precision, thus showing the convergence of the fluid and collisionless kinetic descriptions. We then consider a 3D (electromagnetic) geometry in the gyrotropic (long-wavelength low-frequency) limit and map all closures that are available at the 4th-order moment level. In the Appendix A, we provide comprehensive tables with Pade approximants of $R(zeta)$ up to the 8th-pole order, with many given in an analytic form.
A Landau fluid model for a collisionless electron-proton magnetized plasma, that accurately reproduces the dispersion relation and the Landau damping rate of all the magnetohydrodynamic waves, is presented. It is obtained by an accurate closure of the hydrodynamic hierarchy at the level of the fourth order moments, based on linear kinetic theory. It retains non-gyrotropic corrections to the pressure and heat flux tensors up to the second order in the ratio between the considered frequencies and the ion cyclotron frequency.
Reduced fluid models for collisionless plasmas including electron inertia and finite Larmor radius corrections are derived for scales ranging from the ion to the electron gyroradii. Based either on pressure balance or on the incompressibility of the electron fluid, they respectively capture kinetic Alfven waves (KAWs) or whistler waves (WWs), and can provide suitable tools for reconnection and turbulence studies. Both isothermal regimes and Landau fluid closures permitting anisotropic pressure fluctuations are considered. For small values of the electron beta parameter $beta_e$, a perturbative computation of the gyroviscous force valid at scales comparable to the electron inertial length is performed at order $O(beta_e)$, which requires second-order contributions in a scale expansion. Comparisons with kinetic theory are performed in the linear regime. The spectrum of transverse magnetic fluctuations for strong and weak turbulence energy cascades is also phenomenologically predicted for both types of waves. In the case of moderate ion to electron temperature ratio, a new regime of KAW turbulence at scales smaller than the electron inertial length is obtained, where the magnetic energy spectrum decays like $k_perp^{-13/3}$, thus faster than the $k_perp^{-11/3}$ spectrum of WW turbulence.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
