In order to study chromospheric magnetosonic wave propagation including, for the first time, the effects of ion-neutral interactions in the partially ionized solar chromosphere, we have developed a ...new multi-fluid computational model accounting for ionization and recombination reactions in gravitationally stratified magnetized collisional media. The two-fluid model used in our 2D numerical simulations treats neutrals as a separate fluid and considers charged species (electrons and ions) within the resistive MHD approach with Coulomb collisions and anisotropic heat flux determined by Braginskiis transport coefficients. The electromagnetic fields are evolved according to the full Maxwell equations and the solenoidality of the magnetic field is enforced with a hyperbolic divergence-cleaning scheme. The initial density and temperature profiles are similar to VAL III chromospheric model in which dynamical, thermal, and chemical equilibrium are considered to ensure comparison to existing MHD models and avoid artificial numerical heating. In this initial setup we include simple homogeneous flux tube magnetic field configuration and an external photospheric velocity driver to simulate the propagation of MHD waves in the partially ionized reactive chromosphere. In particular, we investigate the loss of chemical equilibrium and the plasma heating related to the steepening of fast magnetosonic wave fronts in the gravitationally stratified medium.
This paper describes the main features of a pioneering unsteady solver for simulating ideal two-fluid plasmas on unstructured grids, taking profit of GPGPU (General-purpose computing on graphics ...processing units). The code, which has been implemented within the open source COOLFluiD platform, is implicit, second-order in time and space, relying upon a Finite Volume method for the spatial discretization and a three-point backward Euler for the time integration. In particular, the convective fluxes are computed by a multi-fluid version of the AUSM+up scheme for the plasma equations, in combination with a modified Rusanov scheme with tunable dissipation for the Maxwell equations. Source terms are integrated with a one-point rule, using the cell-centered value. Some critical aspects of the porting to GPU’s are discussed, as well as the performance of two open source linear system solvers (i.e. PETSc, PARALUTION). The code design allows for computing both flux and source terms on the GPU along with their Jacobian, giving a noticeable decrease in the computational time in comparison with the original CPU-based solver. The code has been tested in a wide range of mesh sizes and in three different systems, each one with a different GPU. The increased performance (up to 14x) is demonstrated in two representative 2D benchmarks: propagation of circularly polarized waves and the more challenging Geospace Environmental Modeling (GEM) magnetic reconnection challenge.
Arc jets are unique facilities used to evaluate the performance of Thermal Protection Systems (TPS) for hypersonic vehicles. They produce high pressure and high enthalpy plasma flow to simulate the ...extreme heat encountered during atmospheric entry. The constricted arc heater part of an arc jet increases the test gas temperature by Joule heating. This study details the development of the three-dimensional unsteady plasma flow analysis tool, ARCHeS (ARC Heater Simulator). Coupled Navier-Stokes, radiative transfer and Maxwell equations yield current density, magnetism, radiation, and flow field solutions. The present work constitutes the first demonstration of an unsteady three-dimensional plasma flow simulation of high pressure and high enthalpy arc heater that captures kink instabilities of the electric arc. It is found that the arc attachment is mainly driven by upstream arc instabilities.
We study a 1D geometry of a plasma confined between two conducting floating walls with applications to laboratory plasmas. These plasmas are characterized by a quasi-neutral bulk that is joined to ...the wall by a thin boundary layer called sheath that is positively charged. Although analytical solutions are available in the sheath and the pre-sheath, joining the two areas by one analytical solution is still an open problem which requires the numerical resolution of the fluid equations coupled to Poisson equation. Current numerical schemes use high-order discretizations to correctly capture the electron current in the sheath, presenting unsatisfactory results in the boundary layer and they are not adapted to all the possible collisional regimes. In this work, we identify the main numerical challenges that arise when attempting the simulations of such configuration and we propose explanations for the observed phenomena via numerical analysis. We propose a numerical scheme with controlled diffusion as well as new discrete boundary conditions that address the identified issues.
We present a novel numerical model that simulates ideal two-fluid plasmas coupled to the full set of Maxwell’s equations with application to space and laboratory plasmas. We use a fully-implicit ...finite volume method for unstructured meshes, that uses an advection upstream splitting method (i.e., AUSM+-up) for all speeds to discretize the numerical fluxes of the fluids. In addition, we discretize Maxwell’s equations with a modified-Rusanov scheme. The electromagnetic numerical dissipation is scaled using the scales of the fluid-electromagnetics coupled problem that are found to be very different from those of the uncoupled problem. Our numerical scheme guarantees that the elliptical constraints of Maxwell’s equations are satisfied by using hyperbolic divergence cleaning. We validate the performance and accuracy of our model by simulating the following conventional cases: a circularly polarized wave, a Brio–Wu type shock tube, and a two-fluid plasma reconnection with the GEM challenge set up. Our model reveals the complexity of the two-fluid model compared to magnetohydrodynamics (MHD) models, as the inclusion of charge separation, the displacement current and the electron dynamics present are ignored by the MHD simplifications. The two-fluid model shows the presence of electromagnetic and plasma waves and the effect that they have in even the simplest cases. We also compare our model to other available two-fluid models and find our results to be in good agreement.