The key challenge of our project is to perform multi-scale simulations of the solar wind (SW) from the solar surface to the boundary of the heliosphere, and its interaction with the local interstellar medium (LISM). Our simulations have both fundamental and practical applications: (1) the physics of SW-LISM interaction is common to many astrophysical objects and fusion plasma, (2) description of the SW flow and magnetic field at Earth and planets allows us to predict hazardous conditions affecting humans and electronics in space. Our accomplishments: (1) SW simulations with pickup ions (PUIs) and turbulence included, (2) development of a data-driven model for coronal mass ejections; and (3) modeling the SW-LISM interaction to interpret observations from Voyagers and IBEX spacecraft. MS-FLUKSS solves self-consistently the system of MHD equations for the ion-electron mixture, kinetic Boltzmann equation for neutral atoms, and two separate systems for PUIs and turbulence generated by them.