Numerical modeling of transient phenomena in the distant solar wind and in the heliosheath

In this paper we summarize our group's efforts to develop a self-consistent, data-driven model of the solar wind interaction with the local interstellar medium. In this model, we describe the motion of plasma with the MHD approach, while the transport of neutral atoms is described by either kinetic or multi-fluid equations. The model and its implementation in the Multi-Scale Fluid-Kinetic Simulation Suite (MS-FLUKSS) is thoroughly tested and validated by comparing our results with other models and spacecraft measurements. The model appears to be rather accurate in reproducing the Interstellar Boundary Explorer ribbon and the general features of the plasma flow available from the Voyager 1 and 2 in-situ data. In particular, the timing of the Voyager crossing of the termination shock is reasonably accurately reproduced without involving exotically strong interstellar magnetic field values. We added a turbulence model to our governing equations and started treating pickup ions as a separate fluid, which made it possible to greatly improve our ability to reproduce the proton temperature. We discuss the results of our modeling of corotating interaction regions, global merged interaction regions, and solar cycle effects. A possible explanation is provided for the low radial velocity measurements from Voyager 1.