Three-dimensional Modeling of the Solar Wind and Local Interstellar Medium Interaction with Pickup Ions in the Presence of Heliospheric Current Sheet

Our three-dimensional, time-dependent, multi-fluid model has been used to investigate the solar wind (SW)-local interstellar medium (LISM) interaction with pickup ions (PUIs) treated as a separate fluid. A non-zero, but fixed, angle between the Sun's magnetic and rotation axis is adopted. The flow of the plasma mixture (thermal SW protons, PUIs, and electrons), is described by the system of ideal magnetohydrodynamic equations with the source terms responsible for charge exchange between ions and neutral atoms. Different populations of neutral atoms are governed by the individual sets of the Euler equations. As the standard Rankine-Hugoniot relations are not appropriate to describe the anisotropic behavior of PUIs at the termination shock, we use a kinetically-derived set of boundary conditions at it. We extend our previous work [1] and perform these new simulations on a Cartesian grid. This approach allows us to maintain a uniform grid resolution in all directions, without compromising resolution, at large distances from the Sun. The possibility of transition of the SW flow to a stochastic regime in the region between the termination shock and heliopause is further investigated.