Self-consistent solar wind-LISM interaction including neutrals: A 4-fluid model

We investigate the interaction between the heliosphere and local interstellar medium, including charge-exchange self-consistently between the neutral and ionized components. The system is modelled as consisting of neutral hydrogen, protons and electrons. The electrons are assumed to contribute to the dynamics via a pressure that is everywhere equal to the proton pressure. The fluid equations for the protons and hydrogen are solved self-consistently, assuming the two species interact only via charge exchange. The protons are treated as a single isotropic fluid. The hydrogen is decomposed into three isotropic fluids: hydrogen of interstellar origin with characteristic temperature of approximately 104 K, hydrogen produced via charge exchange in the post-shock solar wind, with characteristic temperature of approximately 105 K, and hydrogen produced via charge exchange in the supersonic solar wind, relatively cool but expanding radially outward with the solar wind speed. Two-Dimensional time-dependent numerical simulations of this system are performed to calculate the number density, temperature and bulk velocity of neutral hydrogen and plasma everywhere in the heliosphere. These quantities are compared to relevant solar wind observations. Implications for the pick-up process in the supersonic solar wind, as well as heating and deceleration of the wind, are addressed.