Spatial variation of the supersonic thermal plasma flow downstream of the termination shock

In this paper we start from the most recently observed fact that the solar wind plasma after passage over the termination shock is still supersonic with a Mach number of about 2. To explain this unexpected phenomenon and to predict the evolution of properties of the downstream plasma flow we here consider a two-fluid proton plasma with pick-up protons as a separate suprathermal, second proton fluid. We then formulate a self-consistent system of hydrodynamical conservation equations coupling the two fluids by dynamical and thermodynamical coupling terms and taking into account the effects of newly incorporated protons due to charge exchange with the H-atoms in the heliosheath. This then allows us to predict that in the most probable case the solar wind protons will become subsonic over a distance of about 30 AU downstream of the shock. As we can also show, it may, however, happen that the plasma mixture later again reconverts to a supersonic signature and has to undergo a second shock before meeting the heliopause.