MS-FLUKSS and Its Application to Modeling Flows of Partially Ionized Plasma in the Heliosphere

Flows of partially ionized plasma are frequently characterized by the presence of both thermal and nonthermal populations of ions. This occurs, e. g., in the outer heliosphere -- the part of interstellar space beyond the solar system whose properties are determined by the solar wind (SW) interaction with the local interstellar medium (LISM). Understanding the behavior of such flows requires us to investigate a variety of physical phenomena occurring throughout the solar system. These include charge exchange processes between neutral and charged particles, the birth of pick-up ions (PUIs), the origin of energetic neutral atoms (ENAs), SW turbulence, etc. Collisions between atoms and ions in the heliospheric plasma are so rare that they should be modeled kinetically. PUIs born when LISM neutral atoms charge-exchange with SW ions represent a hot, non-equilibrium component and also require a kinetic treatment. The behavior of PUIs at the SW termination shock (TS) is of major importance for the interpretation of the puzzling data from the Voyager 1 and 2 spacecraft, which are now the only in situ space mission intended to investigate the boundary of the solar system. We have recently proposed an explanation of the sky-spanning \"ribbon\" of unexpectedly intense emissions of ENAs detected by the Interstellar Boundary Explorer (IBEX) mission. Numerical solution of these problems with the realistic boundary conditions provided by remote and in situ observations of the SW properties, requires the application of adaptive mesh refinement (AMR) technologies and petascale supercomputers. Supported by the NSF ITR program and various NASA projects, we have implemented these in our Multi-Scale FLUid-Kinetic Simulation Suite, which is a collection of problem-oriented routines incorporated into the Chombo AMR framework. For the next 5--10 years, heliophysics research is faced with an extraordinary opportunity that cannot be soon repeated. This is to make in situ measurements of the SW from the Sun to the heliospheric boundaries and, at the same time, extract information about the global behavior of the evolving heliosphere through ENA observations by IBEX. In this paper, we describe the application of new possibilities provided within our Extreme Science and Engineering Discovery Environment (XSEDE) project to model challenging space physics and astrophysics problems. We used XSEDE supercomputers to analyze flows of magnetized, rarefied, partially-ionized plasma, where neutral atoms experience resonant charge exchange and collisions with ions. We modeled the SW flows in the inner and outer heliosphere and compared our results with in situ measurements performed by the ACE, IBEX, and Voyager spacecraft.