The heliosphere is formed due to interaction between the solar wind (SW) and local interstellar medium (LISM). The shape and position of the heliospheric boundary, the heliopause, in space depend on the parameters of interacting plasma flows. The interplay between the asymmetrizing effect of the interstellar magnetic field and charge exchange between ions and neutral atoms plays an important role in the SW-LISM interaction. By performing three-dimensional, MHD plasma / kinetic neutral atom simulations, we determine the width of the outer heliosheath - the LISM plasma region affected by the presence of the heliosphere - and analyze quantitatively the distributions in front of the heliopause. It is shown that charge exchange modifies the LISM plasma to such extent that the contribution of a shock transition to the total variation of plasma parameters becomes small even if the LISM velocity exceeds the fast magnetosonic speed in the unperturbed medium. By performing adaptive mesh refinement simulations, we show that a distinct boundary layer of decreased plasma density and enhanced magnetic field should be observed on the interstellar side of the heliopause. We show that this behavior is in agreement with the plasma oscillations of increasing frequency observed by the plasma wave instrument onboard Voyager 1. We also demonstrate that Voyager observations in the inner heliosheath between the heliospheric termination shock and the heliopause are consistent with dissipation of the heliospheric magnetic field. The choice of LISM parameters in this analysis is based on the simulations that fit observations of energetic neutral atoms performed by IBEX.

The magnetometer (MAG) on Voyager 1 (V1) has been sampling the interstellar magnetic field (ISMF) since August 2012. The V1 MAG observations have shown draped ISMF in the very local interstellar medium disturbed occasionally by significant enhancements in magnetic field strength. Using a three-dimensional, data driven, multi-fluid model, we investigated these magnetic field enhancements beyond the heliopause that are supposedly associated with solar transients. To introduce time-dependent effects at the inner boundary at 1 astronomical unit, we used daily averages of the solar wind parameters from the OMNI data set. The model ISMF strength, direction, and proton number density are compared with V1 data beyond the heliopause. The model reproduced the large-scale fluctuations between 2012.652 and 2016.652, including major events around 2012.9 and 2014.6. The model also predicts shocks arriving at V1 around 2017.395 and 2019.502. Another model driven by OMNI data with interplanetary coronal mass ejections (ICMEs) removed at the inner boundary suggests that ICMEs may play a significant role in the propagation of shocks into the interstellar medium.