Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userEnergy conversion in turbulent weakly-collisional plasmas: Eulerian Hybrid Vlasov-Maxwell simulations
111
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
Kinetic simulations based on the Eulerian Hybrid Vlasov-Maxwell (HVM) formalism permit the examination of plasma turbulence with useful resolution of the proton velocity distribution function (VDF). The HVM model is employed here to study the balance of energy, focusing on channels of conversion that lead to proton kinetic effects, including growth of internal energy and temperature anisotropies. We show that this Eulerian simulation approach, which is almost noise-free, is able to provide an accurate energy balance for protons. The results demonstrate explicitly that the recovered temperature growth is directly related to the role of the pressure-strain interaction. Furthermore, analysis of local spatial correlations indicates that the pressure-strain interaction is qualitatively associated with strong-current, high-vorticity structures, although other local terms -- such as the heat flux -- weaken the correlation. These numerical capabilities based on the Eulerian approach will enable deeper study of transfer and conversion channels in weakly collisional Vlasov plasmas.
rate research
Read More
The physical foundations of the dissipation of energy and the associated heating in weakly collisional plasmas are poorly understood. Here, we compare and contrast several measures that have been used to characterize energy dissipation and kinetic-scale conversion in plasmas by means of a suite of kinetic numerical simulations describing both magnetic reconnection and decaying plasma turbulence. We adopt three different numerical codes that can also include interparticle collisions: the fully kinetic particle-in-cell VPIC, the fully kinetic continuum Gkeyll, and the Eulerian Hybrid Vlasov-Maxwell (HVM) code. We differentiate between (i) four energy-based parameters, whose definition is related to energy transfer in a fluid description of a plasma, and (ii) four distribution function-based parameters, requiring knowledge of the particle velocity distribution function. There is an overall agreement between the dissipation measures obtained in the PIC and continuum reconnection simulations, with slight differences due to the presence/absence of secondary islands in the two simulations. There are also many qualitative similarities between the signatures in the reconnection simulations and the self-consistent current sheets that form in turbulence, although the latter exhibits significant variations compared to the reconnection results. All the parameters confirm that dissipation occurs close to regions of intense magnetic stresses, thus exhibiting local correlation. The distribution function-based measures show a broader width compared to energy-based proxies, suggesting that energy transfer is co-localized at coherent structures, but can affect the particle distribution function in wider regions. The effect of interparticle collisions on these parameters is finally discussed.
The mechanism of heating for hot, dilute, and turbulent plasmas represents a long-standing problem in space physics, whose implications concern both near-Earth environments and astrophysical systems. In order to explore the possible role of interparticle collisions, simulations of plasma turbulence -- in both collisionless and weakly collisional regimes -- have been compared by adopting Eulerian Hybrid Boltzmann-Maxwell simulations, being proton-proton collisions explicitly introduced through the nonlinear Dougherty operator. Although collisions do not significantly influence the statistical characteristics of the turbulence, they dissipate nonthermal features in the proton distribution function and suppress the enstrophy/entropy cascade in the velocity space, damping the spectral transfer toward large Hermite modes. This enstrophy dissipation is particularly effective in regions where the plasma distribution function is strongly distorted, suggesting that collisional effects are enhanced by fine velocity-space structures. A qualitative connection between the turbulent energy cascade in fluids and the enstrophy cascade in plasmas has been established, opening a new path to the understanding of astrophysical plasma turbulence
This work describes a new 1D hybrid approach for modeling atmospheric pressure discharges featuring complex chemistry. In this approach electrons are described fully kinetically using Particle-In-Cell/Monte-Carlo (PIC/MCC) scheme, whereas the heavy species are modeled within a fluid description. Validity of the popular drift-diffusion approximation is verified against a full fluid model accounting for the ion inertia and a fully kinetic PIC/MCC code for ions as well as electrons. The fluid models require knowledge of the momentum exchange frequency and dependence of the ion mobilities on the electric field when the ions are in equilibrium with the latter. To this end an auxiliary Monte-Carlo scheme is constructed. It is demonstrated that the drift-diffusion approximation can overestimate ion transport in simulations of RF-driven discharges with heavy ion species operated in the $gamma$ mode at the atmospheric pressure or in all discharge simulations for lower pressures. This can lead to exaggerated plasma densities and incorrect profiles provided by the drift-diffusion models. Therefore, the hybrid code version featuring the full ion fluid model should be favored against the more popular drfit-diffusion model, noting that the suggested numerical scheme for the former model implies only a small additional computational cost.
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.
Vlasiator is a new hybrid-Vlasov plasma simulation code aimed at simulating the entire magnetosphere of the Earth. The code treats ions (protons) kinetically through Vlasovs equation in the six-dimensional phase space while electrons are a massless charge-neutralizing fluid [M. Palmroth et al., Journal of Atmospheric and Solar-Terrestrial Physics 99, 41 (2013); A. Sandroos et al., Parallel Computing 39, 306 (2013)]. For first global simulations of the magnetosphere, it is critical to verify and validate the model by established methods. Here, as part of the verification of Vlasiator, we characterize the low-$beta$ plasma wave modes described by this model and compare with the solution computed by the Waves in Homogeneous, Anisotropic Multicomponent Plasmas (WHAMP) code [K. Ronnmark, Kiruna Geophysical Institute Reports 179 (1982)], using dispersion curves and surfaces produced with both programs. The match between the two fundamentally different approaches is excellent in the low-frequency, long wavelength range which is of interest in global magnetospheric simulations. The left-hand and right-hand polarized wave modes as well as the Bernstein modes in the Vlasiator simulations agree well with the WHAMP solutions. Vlasiator allows a direct investigation of the importance of the Hall term by including it in or excluding it from Ohms law in simulations. This is illustrated showing examples of waves obtained using the ideal Ohms law and Ohms law including the Hall term. Our analysis emphasizes the role of the Hall term in Ohms law in obtaining wave modes departing from ideal magnetohydrodynamics in the hybrid-Vlasov model.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
