No Arabic abstract
Incorporation of kinetic effects such as Landau damping into a fluid framework was pioneered by Hammett and Perkins PRL 1990, by obtaining closures of the fluid hierarchy, where the gyrotropic heat flux fluctuations or the deviation of the 4th-order gyrotropic fluid moment, are expressed through lower-order fluid moments. To obtain a closure of a fluid model expanded around a bi-Maxwellian distribution function, the usual plasma dispersion function $Z(zeta)$ that appears in kinetic theory or the associated plasma response function $R(zeta)=1 + zeta Z(zeta)$, have to be approximated with a suitable Pade approximant in such a way, that the closure is valid for all $zeta$ values. Such closures are rare, and the original closures of Hammett and Perkins are often employed. Here we present a complete mapping of all plausible Landau fluid closures that can be constructed at the level of 4th-order moments in the gyrotropic limit and we identify the most precise closures. Furthermore, by considering 1D closures at higher-order moments, we show that it is possible to reproduce linear Landau damping in the fluid framework to any desired precision, thus showing convergence of the fluid and collisionless kinetic descriptions.
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.
We discuss the self-consistent dynamics of plasmas by means of hamiltonian formalism for a system of $N$ near-resonant electrons interacting with a single Langmuir wave. The connection with the Vlasov description is revisited through the numerical calculation of the van Kampen-like eigenfrequencies of the linearized dynamics for many degrees of freedom. Both the exponential-like growth as well as damping of the Langmuir wave are shown to emerge from a phase mixing effect among beam modes, revealing unexpected similarities between the stable and unstable regimes.
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.
Kinetic treatments of drift-tearing modes that match an inner resonant layer solution to an external MHD region solution, characterised by $Delta^{prime}$, fail to properly match the ideal MHD boundary condition on the parallel electric field, $E_{parallel}.$ In this paper we demonstrate how consideration of ion sound and ion Landau damping effects achieves this and place the theory on a firm footing. As a consequence, these effects contribute quite significantly to the critical value of $Delta^{prime}$ for instability of drift-tearing modes and play a key role in determining the minimum value for this threshold.
This paper brings further insight into the recently published N-body description of Debye shielding and Landau damping [Escande D F, Elskens Y and Doveil F 2014 Plasma Phys. Control. Fusion 57 025017]. Its fundamental equation for the electrostatic potential is derived in a simpler and more rigorous way. Various physical consequences of the new approach are discussed, and this approach is compared with the seminal one by Pines and Bohm [Pines D and Bohm D 1952 Phys. Rev. 85 338--353].