No Arabic abstract
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.
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.
A massively parallel simulation code, called textit{dHybrid}, has been developed to perform global scale studies of space plasma interactions. This code is based on an explicit hybrid model; the numerical stability and parallel scalability of the code are studied. A stabilization method for the explicit algorithm, for regions of near zero density, is proposed. Three-dimensional hybrid simulations of the interaction of the solar wind with unmagnetized artificial objects are presented, with a focus on the expansion of a plasma cloud into the solar wind, which creates a diamagnetic cavity and drives the Interplanetary Magnetic Field out of the expansion region. The dynamics of this system can provide insights into other similar scenarios, such as the interaction of the solar wind with unmagnetized planets.
An axisymmetric fully ionized plasma rotates around its axis when a charge separation between magnetic surfaces is produced from DC fields or RF waves. On each magnetic surface both electrons and ions obey the isorotation law and perform an azimuthal E cross B rotation at the same angular velocity. When Coulomb collisions are taken into account such a flow displays no Ohmic current short circuiting of the charge separation and thus no linear dissipation. A nonlinear Ohmic response appears when inertial effects are considered, providing a dissipative relaxation of the charge separation between the magnetic surfaces. This nonlinear conductivity results from an interplay between Coriolis, centrifugal and electron-ion collisional friction forces. This phenomena is identified, described and analyzed. In addition, both the quality factor of angular momentum storage as well as the efficiency of wave driven angular momentum generation are calculated and shown to be independent of the details of the charge separation processes.
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.
We propose a new approach to high-intensity laser-driven electron acceleration in a plasma. Here, we demonstrate that a plasma wave generated by a stimulated forward-scattering of an incident laser pulse can be in a longest acceleration phase with an incident laser wave. This is why the plasma wave has the maximum amplification coefficient which is determined by the breakdown (overturn) electric field in which the acceleration of injected relativistic beam electrons occurs. We estimate qualitatively the acceleration parameters of relativistic electrons in the field of a plasma wave generated at the stimulated forward scattering of a high-intensity laser pulse in a plasma.