We performed three-dimensional magnetohydrodynamic simulations to study the evolution of a supernova remnant (SNR) in a turbulent neutral atomic interstellar medium. The media used as background shares characteristics with the Solar neighbourhood and the SNR has mass and energy similar to those of a Type Ia object. Our initial conditions consist of dense clouds in a diluted medium, with the main difference between simulations being the average magnitude of the magnetic field. We measured amplifications of the magnetic energy of up to 34$\%$ and we generated synthetic maps that illustrate how the same object can show different apparent geometries and physical properties when observed through different lines of sight.

We present numerical experiments aimed to study the correlation between the magnetic field strength, $B$, and the density, $n$, in the cold atomic interstellar medium (CNM). We analyze 24 magneto-hydrodynamic models with different initial magnetic field intensities ($B_0=$0.4, 2.1, 4.2, and 8.3 $\mu$G) and/or mean densities (2, 3, and 4 cm$^{-3}$), in the presence of driven and decaying turbulence, with and without self-gravity, in a cubic computational domain with 100 pc by side. Our main findings are: i) For forced simulations, which reproduce the main observed physical conditions of the CNM in the Solar neighborhood, a positive correlation between $B$ and $n$ develops for all the $B_0$ values. ii) The density at which this correlation becomes significant ($\lesssim 30$ cm$^{-3}$) depends on $B_0$ but is not sensitive to the presence of self-gravity. iii) The effect of self-gravity, when noticeable, consists of producing a shallower correlation at high densities, suggesting that, in the studied regime, self-gravity induces motions along the field lines. iv) Self-gravitating decaying models where the CNM is subsonic and sub-Alfv\'enic with $\beta\lesssim 1$ develop a high density positive correlation whose slopes are consistent with a constant $\beta(n)$. v) Decaying models where the low density CNM is subsonic and sub-Alfv\'enic with $\beta>1$ show a negative correlation at intermediate densities, followed by a high density positive correlation.

To contribute to the understanding of the magnetic field's influence on the segregation of CNM in the solar neighbourhood we analyse MHD simulations which include the main physical characteristics of the local neutral atomic ISM. The simulations have a continuous solenoidal Fourier forcing in a periodic box of 100 pc per side and an initial uniform magnetic field ($\vec{B_0}$) with intensities ranging between $\sim 0.4$ $\mu$G and $\sim 8$ $\mu$G. Our main results are: i) the CNM mass fraction diminishes with the increase in magnetic field intensity. ii) There is a preferred alignment between CNM structures and $\vec{B}$ in all our $B_0$ range but the preference weakens as $B_{0}$ increases. It is worth noticing that this preference is also present in two-dimensional projections making an extreme angle ($0$ or $\pi / 2$) with respect to $\vec{B_0}$ and it is only lost for the strongest magnetic field when the angle of projection is perpendicular to $\vec{B_0}$. iii) The aforementioned results are prevalent despite the inclusion of self-gravity in our continuously forced simulations with a mean density similar to the average value of the solar neighbourhood. iv) Given a fixed $B_0$ and slightly higher mean densities, up to double, the effects of self-gravity are still not qualitatively significant.