Stellar feedback drives the circulation of matter from the disk to the halo of galaxies. We perform three-dimensional magnetohydrodynamic simulations of a vertical column of the interstellar medium with initial conditions typical of the solar circle in which supernovae drive turbulence and determine the vertical stratification of the medium. The simulations were run using a stable, positivity-preserving scheme for ideal MHD implemented in the FLASH code. We find that the majority (\approx 90 %) of the mass is contained in thermally-stable temperature regimes of cold molecular and atomic gas at T < 200 K or warm atomic and ionized gas at 5000 K < T < 10^{4.2} K, with strong peaks in probability distribution functions of temperature in both the cold and warm regimes. The 200 - 10^{4.2} K gas fills 50-60 % of the volume near the plane, with hotter gas associated with supernova remnants (30-40 %) and cold clouds (< 10 %) embedded within. At |z| ~ 1-2 kpc, transition-temperature (10^5 K) gas accounts for most of the mass and volume, while hot gas dominates at |z| > 3 kpc. The magnetic field in our models has no significant impact on the scale heights of gas in each temperature regime; the magnetic tension force is approximately equal to and opposite the magnetic pressure, so the addition of the field does not significantly affect the vertical support of the gas. The addition of a magnetic field does reduce the fraction of gas in the cold (< 200 K) regime with a corresponding increase in the fraction of warm (~ 10^4 K) gas. However, our models lack rotational shear and thus have no large-scale dynamo, which reduces the role of the field in the models compared to reality. The supernovae drive oscillations in the vertical distribution of halo gas, with the period of the oscillations ranging from ~ 30 Myr in the T < 200 K gas to ~ 100 Myr in the 10^6 K gas, in line with predictions by Walters & Cox.

We present new 3D magnetohydrodynamic (MHD) simulations of a supernova-driven, stratified interstellar medium. These simulations were run using the Waagan (2009) positivity preserving scheme for ideal MHD implemented in the Flash code. The scheme is stable even for the Mach numbers approaching 100 found in this problem. We have previously shown that the density distribution arising from hydrodynamical versions of these simulations creates low-density pathways through which Lyman continuum photons can travel to heights |z| > 1 kpc. This naturally produces the warm ionized medium through photoionization due primarily to O stars near the plane. However, our earlier models reproduce the peak but not the width of the observed emission measure distribution. Here, we examine whether inclusion of magnetic fields and a greater vertical extent to the simulation domain produce a gas distribution that better matches the observations. We further study the change of magnetic energy over time in our models, showing that it appears to reach a steady state after a few hundred megayears, presumably supported by a turbulent dynamo driven by the supernova explosions.

The ideal MHD equations are a central model in astrophysics, and their solution relies upon stable numerical schemes. We present an implementation of a new method, which possesses excellent stability properties. Numerical tests demonstrate that the theoretical stability properties are valid in practice with negligible compromises to accuracy. The result is a highly robust scheme with state-of-the-art efficiency. The scheme's robustness is due to entropy stability, positivity and properly discretised Powell terms. The implementation takes the form of a modification of the MHD module in the FLASH code, an adaptive mesh refinement code. We compare the new scheme with the standard FLASH implementation for MHD. Results show comparable accuracy to standard FLASH with the Roe solver, but highly improved efficiency and stability, particularly for high Mach number flows and low plasma beta. The tests include 1D shock tubes, 2D instabilities and highly supersonic, 3D turbulence. We consider turbulent flows with RMS sonic Mach numbers up to 10, typical of gas flows in the interstellar medium. We investigate both strong initial magnetic fields and magnetic field amplification by the turbulent dynamo from extremely high plasma beta. The energy spectra show a reasonable decrease in dissipation with grid refinement, and at a resolution of 512^3 grid cells we identify a narrow inertial range with the expected power-law scaling. The turbulent dynamo exhibits exponential growth of magnetic pressure, with the growth rate twice as high from solenoidal forcing than from compressive forcing. Two versions of the new scheme are presented, using relaxation-based 3-wave and 5-wave approximate Riemann solvers, respectively. The 5-wave solver is more accurate in some cases, and its computational cost is close to the 3-wave solver.

The idea of this work is to compare a new positive and entropy stable approximate Riemann solver by Francois Bouchut with a state-of the-art algorithm for astrophysical fluid dynamics. We implemented the new Riemann solver into an astrophysical PPM-code, the Prometheus code, and also made a version with a different, more theoretically grounded higher order algorithm than PPM. We present shock tube tests, two-dimensional instability tests and forced turbulence simulations in three dimensions. We find subtle differences between the codes in the shock tube tests, and in the statistics of the turbulence simulations. The new Riemann solver increases the computational speed without significant loss of accuracy.