In order to study chromospheric magnetosonic wave propagation including, for the first time, the effects of ionneutral interactions in the partially ionized solar chromosphere, we have developed a new multi-fluid computational modelaccounting 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.
We present a reduced kinetic mechanism for the modeling of the behavior of the electronic states of the atomic species in air mixtures. The model is built by lumping the electronically excited states of the atomic species and by performing Maxwell-Boltzmann averages of the rate constants describing the elementary kinetic processes of the individual states within each group. The necessary reaction rate coefficients are taken from the model compiled by Bultel et al. [“Collisional-radiative model in air for earth re-entry problems,” Phys. Plasmas 13, 043502 (2006)10.1063/1.2194827]. The reduced number of pseudo-states considered leads to a significant reduction of the computational cost, thus enabling the application of the state of the art collisional radiative models to bi-dimensional and three-dimensional problems. The internal states of the molecular species are assumed to be in equilibrium. The rotational energy mode is assumed to quickly equilibrate with the translational energy mode at the kinetic temperature of the heavy species as opposed to the electronic and the vibrational modes, assumed to be in Maxwell-Boltzmann equilibrium at a common temperature TV. In a first step we validate the model by using simple zero- and one-dimensional test cases for which the full kinetic mechanism can be run efficiently. Finally, the reduced kinetic model is used to analyze the strong non-equilibrium flow surrounding the FIRE II flight experiment during the early part of its re-entry trajectory. It is found that the reduced kinetic mechanism is capable of reproducing the ionizational non-equilibrium phenomena, responsible for the drastic reduction of the radiative heat loads on the space capsules during the re-entry phase.
We present a novel global 3D coronal MHD model called COCONUT, polytropic in its first stage and based on a time-implicit backward Euler scheme. Our model boosts run-time performance in comparison with contemporary MHD-solvers based on explicit schemes, which is particularly important when later employed in an operational setting for space-weather forecasting. It is data-driven in the sense that we use synoptic maps as inner boundary inputs for our potential-field initialization as well as an inner boundary condition in the further MHD time evolution. The coronal model is developed as part of the EUropean Heliospheric FORecasting Information Asset (EUHFORIA) and will replace the currently employed, more simplistic, empirical Wang–Sheeley–Arge (WSA) model. At 21.5 R
⊙ where the solar wind is already supersonic, it is coupled to EUHFORIA’s heliospheric model. We validate and benchmark our coronal simulation results with the explicit-scheme Wind-Predict model and find good agreement for idealized limit cases as well as real magnetograms, while obtaining a computational time reduction of up to a factor 3 for simple idealized cases, and up to 35 for realistic configurations, and we demonstrate that the time gained increases with the spatial resolution of the input synoptic map. We also use observations to constrain the model and show that it recovers relevant features such as the position and shape of the streamers (by comparison with eclipse white-light images), the coronal holes (by comparison with EUV images), and the current sheet (by comparison with WSA model at 0.1 au).
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.