Multicomponent diffusion and energy characteristics of partially ionized plasma in the ionosphere of a planet

The problem of energy and multicomponent ambipolar diffusion of plasma in the lower ionosphere of a planet with a weak magnetic field is considered

Two-Fluid MHD Simulations of Converging HI Flows in the Interstellar Medium. I: Methodology and Basic Results

We develop an unconditionally stable numerical method for solving the coupling between two fluids (frictional forces/heatings, ionization, and recombination), and investigate the dynamical condensation process of thermally unstable gas that is provided by the shock waves in a weakly ionized and magnetized interstellar medium by using two-dimensional two-fluid magnetohydrodynamical simulations. If we neglect the effect of magnetic field, it is known that condensation driven by thermal instability can generate high density clouds whose physical condition corresponds to molecular clouds (precursor of molecular clouds). In this paper, we study the effect of magnetic field on the evolution of supersonic converging HI flows and focus on the case in which the orientation of magnetic field to converging flows is orthogonal. We show that the magnetic pressure gradient parallel to the flows prevents the formation of high density and high column density clouds, but instead generates fragmented, filamentary HI clouds. With this restricted geometry, magnetic field drastically diminishes the opportunity of fast molecular cloud formation directly from the warm neutral medium, in contrast to the case without magnetic field.Comment: ApJ accepte

A three-dimensional numerical method for modelling weakly ionized plasmas

Astrophysical fluids under the influence of magnetic fields are often subjected to single-fluid or two-fluid approximations. In the case of weakly ionized plasmas however, this can be inappropriate due to distinct responses from the multiple constituent species to both collisional and non-collisional forces. As a result, in dense molecular clouds and proto-stellar accretion discs for instance, the conductivity of the plasma may be highly anisotropic leading to phenomena such as Hall and ambipolar diffusion strongly influencing the dynamics. Diffusive processes are known to restrict the stability of conventional numerical schemes which are not implicit in nature. Furthermore, recent work establishes that a large Hall term can impose an additional severe stability limit on standard explicit schemes. Following a previous paper which presented the one-dimensional case, we describe a fully three-dimensional method which relaxes the normal restrictions on explicit schemes for multifluid processes. This is achieved by applying the little known Super TimeStepping technique to the symmetric (ambipolar) component of the evolution operator for the magnetic field in the local plasma rest-frame, and the new Hall Diffusion Scheme to the skew-symmetric (Hall) component.Comment: 13 pages, 9 figures, accepted for publication in MNRA

Thermodynamic evolution of cosmological baryonic gas: I. Influence of non-equipartition processes

Using N-body/hydrodynamic simulations, the influence of non-equipartition processes on the thermal and dynamical properties of cosmological baryonic gas is investigated. We focus on a possible departure from equilibrium between electrons, ions and neutral atoms in low temperature (10^4-10^6 K) and weakly ionized regions of the intergalactic medium. The simulations compute the energy exchanges between ions, neutrals and electrons, without assuming thermal equilibrium. They include gravitation, shock heating and cooling processes, and follow self-consistently the chemical evolution of a primordial composition hydrogen-helium plasma without assuming collisional ionization equilibrium. At high redshift, a significant fraction of the intergalactic medium is found to be warmer and weakly ionized in simulations with non-equipartition processes than in simulations in which the cosmological plasma is considered to be in thermodynamic equilibrium. With a semi-analytical study of the out of equilibrium regions we show that, during the formation of cosmic structures, departure from equilibrium in accreted plasma results from the competition between the atomic cooling processes and the elastic processes between heavy particles and electrons. Our numerical results are in agreement with this semi-analytical model. Therefore, since baryonic matter with temperatures around 10^4 K is a reservoir for galaxy formation, non-equipartition processes are expected to modify the properties of the objects formed.Comment: 15 pages, 16 figures. Accepted for publication in A&A. For a version with high-resolution figures, see http://www.raunvis.hi.is/~courty/series.htm

Solar wind collisional heating

To properly describe heating in weakly collisional turbulent plasmas such as the solar wind, inter-particle collisions should be taken into account. Collisions can convert ordered energy into heat by means of irreversible relaxation towards the thermal equilibrium. Recently, Pezzi et al. (Phys. Rev. Lett., vol. 116, 2016, p. 145001) showed that the plasma collisionality is enhanced by the presence of fine structures in velocity space. Here, the analysis is extended by directly comparing the effects of the fully nonlinear Landau operator and a linearized Landau operator. By focusing on the relaxation towards the equilibrium of an out of equilibrium distribution function in a homogeneous force-free plasma, here it is pointed out that it is significant to retain nonlinearities in the collisional operator to quantify the importance of collisional effects. Although the presence of several characteristic times associated with the dissipation of different phase space structures is recovered in both the cases of the nonlinear and the linearized operators, the influence of these times is different in the two cases. In the linearized operator case, the recovered characteristic times are systematically larger than in the fully nonlinear operator case, this suggesting that fine velocity structures are dissipated slower if nonlinearities are neglected in the collisional operator