We explore the effect of magnetic fields on the vertical distribution and multiphase structure of the supernova-driven interstellar medium (ISM) in simulations that admit dynamo action. As the magnetic field is amplified to become dynamically significant, gas becomes cooler and its distribution in the disc becomes more homogeneous. We attribute this to magnetic quenching of vertical velocity, which leads to a decrease in the cooling length of hot gas. A non-monotonic vertical distribution of the large-scale magnetic field strength, with the maximum at |z| $\approx$ 300 pc causes a downward pressure gradient below the maximum which acts against outflow driven by SN explosions, while it provides pressure support above the maximum.

We examine the effect of uniform ambient magnetic fields on the evolution of supernova-driven blast waves into a homogeneous ambient ISM in thermal equilibrium. Using the Pencil Code we simulate high resolution nonideal magnetohydrodynamic simulations in 3D. We find that supernova blast waves are sensitive to plane-parallel magnetic fields of strength in excess of 1 $\mu$G for ambient gas number density 1 cm$^{-3}$ . Perpendicular to the field, the inward magnetic pressure gradient induces retrograde mass accretion in the wake of the primary shock front. Subsequently, we find that the primary shockwave expands faster perpendicular to the field, but with reduced momentum, while the remnant core is subject to magnetic confinement. This leads to a decrease in fractional volume of hot gas but also an increase in the density and temperature of hot gas in the magnetically confined remnant. The magnetic pressure gradient behind the shock front generates enhanced regions favourable to UV- heating and thus reduces net radiative losses. Although the presence of a strong uniform magnetic field can reduce momentum early on, and hence residual kinetic energy, it increases the efficiency of residual total energy injection by the SN into the ISM by up to 40% within 1 Myr.

The magnetic network extending from the photosphere (solar radius $\simeq R_\odot$) to lower corona ($R_\odot + 10$ Mm) plays an important role in the heating mechanisms of the solar atmosphere. Here we further develop the models with realistic open magnetic flux tubes of the authors in order to model more complicated configurations. Closed magnetic loops, and combinations of closed and open magnetic flux tubes are modelled. These are embedded within a stratified atmosphere, derived from observationally motivated semi-empirical and data-driven models subject to solar gravity and capable of spanning from the photosphere up into the chromosphere and lower corona. Constructing a magnetic field comprising self-similar magnetic flux tubes, an analytic solution for the kinetic pressure and plasma density is derived. Combining flux tubes of opposite polarity it is possible to create a steady background magnetic field configuration modelling a solar atmosphere exhibiting realistic stratification. The result can be applied to SOHO/MDI and SDO/HMI and other magnetograms from the solar surface, upon which photospheric motions can be simulated to explore the mechanism of energy transport. We demonstrate this powerful and versatile method with an application to Helioseismic and Magnetic Imager data.

High Mach number shocks are ubiquitous in interstellar turbulence. The Pencil Code is particularly well suited to the study of magnetohydrodynamics in weakly compressible turbulence and the numerical investigation of dynamos because of its high-order advection and time evolution algorithms. However, the high-order algorithms and lack of Riemann solver to follow shocks make it less well suited to handling high Mach number shocks, such as those produced by supernovae (SNe). Here, we outline methods required to enable the code to efficiently and accurately model SNe, using parameters that allow stable simulation of SN-driven turbulence, in order to construct a physically realistic galactic dynamo model. These include the resolution of shocks with artificial viscosity, thermal conductivity, and mass diffusion; the correction of the mass diffusion terms; and a novel generalization of the Courant condition to include all source terms in the momentum and energy equations. We test our methods with the numerical solution of the one-dimensional (1D) Riemann shock tube (Sod, J. Comput. Phys. 1978, 27), also extended to a 1D adiabatic shock with parameters and Mach number relevant to SN shock evolution, including shocks with radiative losses. We extend our test with the three-dimensional (3D) numerical simulation of individual SN remnant evolution for a range of ambient gas densities typical of the interstellar medium and compare these to the analytical solutions of Sedov-Taylor (adiabatic) and the snowplough and Cioffi, McKee and Bertschinger (Astrophys. J. 1988, 334) results incorporating cooling and heating processes. We show that our new timestep algorithm leads to linear rather than quadratic resolution dependence as the strength of the artificial viscosity varies, because of the corresponding change in the strength of interzone gradients.

We test the sensitivity of hydrodynamic and magnetohydrodynamic turbulent convection simulations with respect to Mach number, thermal and magnetic boundary conditions, and the centrifugal force. We find that varying the luminosity, which also controls the Mach number, has only a minor effect on the large-scale dynamics. A similar conclusion can also be drawn from the comparison of two formulations of the lower magnetic boundary condition with either vanishing electric field or current density. The centrifugal force has an effect on the solutions, but only if its magnitude with respect to acceleration due to gravity is by two orders of magnitude greater than in the Sun. Finally, we find that the parameterisation of the photospheric physics, either by an explicit cooling term or enhanced radiative diffusion, is more important than the thermal boundary condition. In particular, runs with cooling tend to lead to more anisotropic convection and stronger deviations from the Taylor-Proudman state. In summary, the fully compressible approach taken here with the Pencil Code is found to be valid, while still allowing the disparate timescales to be taken into account.

We present high resolution ($1024^3$) simulations of super-/hyper-sonic isothermal hydrodynamic turbulence inside an interstellar molecular cloud (resolving scales of typically 20 -- 100 AU), including a multi-disperse population of dust grains, i.e., a range of grain sizes is considered. Due to inertia, large grains (typical radius $a \gtrsim 1.0\,\mu$m) will decouple from the gas flow, while small grains ($a\lesssim 0.1\,\mu$m) will tend to better trace the motions of the gas. We note that simulations with purely solenoidal forcing show somewhat more pronounced decoupling and less clustering compared to simulations with purely compressive forcing. Overall, small and large grains tend to cluster, while intermediate-size grains show essentially a random isotropic distribution. As a consequence of increased clustering, the grain-grain interaction rate is locally elevated; but since small and large grains are often not spatially correlated, it is unclear what effect this clustering would have on the coagulation rate. Due to spatial separation of dust and gas, a diffuse upper limit to the grain sizes obtained by condensational growth is also expected, since large (decoupled) grains are not necessarily located where the growth species in the molecular gas is.

The role of magnetic fields in the multi-phase interstellar medium (ISM) is explored using magnetohydrodynamic (MHD) simulations that include energy injection by supernova (SN) explosions and allow for dynamo action. Apart from providing additional pressure support of the gas layer, magnetic fields reduce the density contrast between the warm and hot gas phases and quench galactic outflows. A dynamo-generated, self-consistent large-scale magnetic field affects the ISM differently from an artificially imposed, unidirectional magnetic field.

We apply Gaussian smoothing to obtain mean density, velocity, magnetic and energy density fields in simulations of the interstellar medium based on three-dimensional magnetohydrodynamic equations in a shearing box $1\times1\times2 \, \rm{kpc}$ in size. Unlike alternative averaging procedures, such as horizontal averaging, Gaussian smoothing retains the three-dimensional structure of the mean fields. Although Gaussian smoothing does not obey the Reynolds rules of averaging, physically meaningful central statistical moments are defined as suggested by Germano (1992). We discuss methods to identify an optimal smoothing scale $\ell$ and the effects of this choice on the results. From spectral analysis of the magnetic, density and velocity fields, we find a suitable smoothing length for all three fields, of $\ell \approx 75 \, \rm{pc}$. We discuss the properties of third-order statistical moments in fluctuations of kinetic energy density in compressible flows and suggest their physical interpretation. The mean magnetic field, amplified by a mean-field dynamo, significantly alters the distribution of kinetic energy in space and between scales, reducing the magnitude of kinetic energy at intermediate scales. This intermediate-scale kinetic energy is a useful diagnostic of the importance of SN-driven outflows.

The efforts for comparing dust polarization measurements with synthetic observations from MHD models have previously concentrated on the scale of molecular clouds. Here we extend the model comparisons to kiloparsec scales, taking into account hot shocked gas generated by supernovae, and a non-uniform dynamo-generated magnetic field at both large and small scales down to 4 pc spatial resolution. Radiative transfer calculations are used to model dust emission and polarization on the top of MHD simulations. We compute synthetic maps of column density $N_H$, polarization fraction $p$, and polarization angle dispersion $S$, and study their dependencies on the important properties of the MHD simulations. These include the large-scale magnetic field and its orientation, the small-scale magnetic field, and supernova-driven shocks. Similar filament-like structures of $S$ as seen in the Planck all-sky maps are visible in our synthetic results, although the smallest scale structures are absent from our maps. Supernova-driven shock fronts and $S$ do not show significant correlation. Instead, $S$ can clearly be attributed to the distribution of the small-scale magnetic field. We also find that the large-scale magnetic field influences the polarization properties, such that, for a given strength of magnetic fluctuation, a strong plane-of-the-sky mean field weakens the observed $S$, while strengthening $p$. The anticorrelation of $p$ and $S$, and decreasing $p$ as a function of $N_H$ are consistent across all synthetic observations. The magnetic fluctuations follow an exponential distribution, rather than Gaussian, characteristic of flows with intermittent repetitive shocks. The observed polarization properties and column densities are sensitive to the line of sight distance over which the emission is integrated.

Vortex motions are frequently observed on the solar photosphere. These motions may play a key role in the transport of energy and momentum from the lower atmosphere into the upper solar atmosphere, contributing to coronal heating. The lower solar atmosphere also consists of complex networks of flux tubes that expand and merge throughout the chromosphere and upper atmosphere. We perform numerical simulations to investigate the behaviour of vortex driven waves propagating in a pair of such flux tubes in a non-force-free equilibrium with a realistically modelled solar atmosphere. The two flux tubes are independently perturbed at their footpoints by counter-rotating vortex motions. When the flux tubes merge, the vortex motions interact both linearly and nonlinearly. The linear interactions generate many small-scale transient magnetic substructures due to the magnetic stress imposed by the vortex motions. Thus, an initially monolithic tube is separated into a complex multi-threaded tube due to the photospheric vortex motions. The wave interactions also drive a superposition that increases in amplitude until it exceeds the local Mach number and produces shocks that propagate upwards with speeds of approximately $ 50$ km s$^{-1}$. The shocks act as conduits transporting momentum and energy upwards, and heating the local plasma by more than an order of magnitude, with peak temperature approximately $60,000$ K. Therefore, we present a new mechanism for the generation of magnetic waveguides from the lower solar atmosphere to the solar corona. This wave guide appears as the result of interacting perturbations in neighbouring flux tubes. Thus, the interactions of photospheric vortex motions is a potentially significant mechanism for energy transfer from the lower to upper solar atmosphere.

We apply correlation analysis to random fields in numerical simulations of the supernova-driven interstellar medium (ISM) with the magnetic field produced by dynamo action. We solve the thermo-magnetohydrodynamic (MHD) equations in a shearing, Cartesian box representing a local region of the ISM, subject to thermal and kinetic energy injection by supernova explosions, and parameterized optically-thin radiative cooling. We consider the cold, warm\, and hot phases of the ISM separately; the analysis mostly considers the warm gas, which occupies the bulk of the domain. Various physical variables have different correlation lengths in the warm phase: $40 pc$, $50 pc$, and $60 pc$ for random magnetic field, density, and velocity, respectively, in the midplane. The correlation time of the random velocity is comparable to the eddy turnover time, about $10^7 yr$, although it may be shorter in regions with higher star formation rate. The random magnetic field is anisotropic, with the standard deviations of its components $b_x/b_y/b_z$ having approximate ratios $0.5/0.6/0.6$ in the midplane. The anisotropy is attributed to the global velocity shear from galactic differential rotation, and locally inhomogeneous outflow to the galactic halo. The correlation length of Faraday depth along the $z$-axis, $120 pc$, is greater than for electron density, $60 - 90 pc$, and vertical magnetic field, $60 pc$. Such comparisons may be sensitive to the orientation of the line of sight. Uncertainties of the structure functions of synchrotron intensity rapidly increase with the scale. This feature is hidden in power spectrum analysis, which can undermine the usefulness of power spectra for detailed studies of interstellar turbulence.

The forcing of interstellar turbulence, driven mainly by supernova explosions, is irrotational in nature, but the development of significant amounts of vorticity and helicity, accompanied by large-scale dynamo action, has been reported. Several earlier investigations examined vorticity production in simpler systems; here all the relevant processes can be considered simultaneously. We also investigate the mechanisms for the generation of net helicity and large-scale flow in the system. We use a three-dimensional, stratified, rotating and shearing local simulation domain of the size 1x1x2 kpc$^3$, forced with SN explosions occurring at the rate typical of the solar neighbourhood in the Milky Way. In addition to the nominal simulation run with realistic Milky Way parameters, we vary the rotation and shear rates, but keep the absolute value of their ratio fixed. Reversing the sign of shear vs. rotation allows us to separate the rotation- and shear-generated contributions. As in earlier studies, we find the generation of significant amounts of vorticity, with on average 65% of the kinetic energy being in the rotational modes. The vorticity production can be related to the baroclinicity of the flow, especially in the regions of hot, dilute clustered supernova bubbles. In these regions, the vortex stretching acts as a sink of vorticity. The net helicities produced by rotation and shear are of opposite signs for physically motivated rotation laws, with the solar neighbourhood parameters resulting in the near cancellation of the total net helicity. We also find the excitation of oscillatory mean flows, the strength and oscillation period of which depend on the Coriolis and shear parameters; we interpret these as signatures of the anisotropic kinetic (AKA) effect. We use the method of moments to fit for the turbulent transport coeffcients, and find $\alpha_{\rm AKA}$ values of the order 3-5 km/s.

The Sun, aside from its eleven year sunspot cycle is additionally subject to long term variation in its activity. In this work we analyse a solar-like convective dynamo simulation, containing approximately 60 magnetic cycles, exhibiting equatorward propagation of the magnetic field, multiple frequencies, and irregular variability, including a missed cycle and complex parity transitions between dipolar and quadrupolar modes. We compute the turbulent transport coefficients, describing the effects of the turbulent velocity field on the mean magnetic field, using the test-field method. The test-field analysis provides a plausible explanation of the missing cycle in terms of the reduction of $\alpha_{\phi\phi}$ in advance of the reduced surface activity, and enhanced downward turbulent pumping during the event to confine some of the magnetic field at the bottom of the convection zone, where local maximum of magnetic energy is observed during the event. At the same time, however, a quenching of the turbulent magnetic diffusivities is observed, albeit differently distributed in depth compared to the other transport coefficients. Therefore, dedicated mean-field modelling is required for verification.

I model the multi-phase interstellar medium (ISM) randomly heated and shocked by supernovae, with gravity, differential rotation and other parameters we understand to be typical of the solar neighbourhood. The simulations are 3D extending horizontally 1 x 1 kpc squared and vertically 2 kpc, symmetric about the galactic mid-plane. They routinely span gas number densities 1/10000 to 100 per cubic cm, temperatures 100 to 100 MK, speeds up to 10000 km/s and Mach number up to 25. Radiative cooling is applied from two widely adopted parameterizations, and compared directly to assess the sensitivity of the results to cooling. There is strong evidence to describe the ISM as comprising well defined cold, warm and hot regions, which are statistically close to thermal and total pressure equilibrium. This result is not sensitive to the choice of parameters considered here. The distribution of the gas density within each can be robustly modelled as lognormal. Appropriate distinction is required between the properties of the gases in the supernova active mid-plane and the more homogeneous phases outside this region. The connection between the fractional volume of a phase and its various proxies is clarified. An exact relation is then derived between the fractional volume and the filling factors defined in terms of the volume and probabilistic averages. These results are discussed in both observational and computational contexts. The correlation scale of the random flows is calculated from the velocity autocorrelation function; it is of order 100 pc and tends to grow with distance from the mid-plane. The origin and structure of the magnetic fields in the ISM is also investigated in non-ideal MHD simulations. A seed magnetic field, with volume average of roughly 4 nG, grows exponentially to reach a statistically steady state within 1.6 Gyr.

We explore the effects of the multi-phase structure of the interstellar medium (ISM) on galactic magnetic fields. Basing our analysis on compressible magnetohydrodynamic (MHD) simulations of supernova-driven turbulence in the ISM, we investigate the properties of both the mean and fluctuating components of the magnetic field. We find that the mean magnetic field preferentially resides in the warm phase and is generally absent from the hot phase. The fluctuating magnetic field does not show such pronounced sensitivity to the multi-phase structure.

A system of multiple open magnetic flux tubes spanning the solar photosphere and lower corona is modelled analytically, within a realistic stratified atmosphere subject to solar gravity. This extends results for a single magnetic flux tube in magnetohydrostatic equilibrium, described in Gent et al. (MNRAS, 435, 689, 2013). Self-similar magnetic flux tubes are combined to form magnetic structures, which are consistent with high-resolution observations. The observational evidence supports the existence of strands of open flux tubes and loops persisting in a relatively steady state. Self-similar magnetic flux tubes, for which an analytic solution to the plasma density and pressure distribution is possible, are combined. We calculate the appropriate balancing forces, applying to the equations of momentum and energy conservation to preserve equilibrium. Multiplex flux tube configurations are observed to remain relatively stable for up to a day or more, and it is our aim to apply our model as the background condition for numerical studies of energy transport mechanisms from the solar surface to the corona. We apply magnetic field strength, plasma density, pressure and temperature distributions consistent with observational and theoretical estimates for the lower solar atmosphere. Although each flux tube is identical in construction apart from the location of the radial axis, combinations can be applied to generate a non-axisymmetric magnetic field with multiple non-uniform flux tubes. This is a considerable step forward in modelling the realistic magnetized three-dimensional equilibria of the solar atmosphere.

The hydrodynamic state of the interstellar medium (ISM) heated and randomly stirred by supernovae (SNe) is investigated. We use a three-dimensional non-ideal hydrodynamic ISM model in a domain extending 0.5 x 0.5 kpc horizontally and 2 kpc vertically to explore the relative importance of various physical and numerical effects on the multi-phase, turbulent ISM. We include both Type I and II SNe, the latter occurring only in dense regions. First we investigate the role of the thermal instability in the temperature range 300-6100 K, comparing results obtained for two different cooling functions, one susceptible to the instability, the other stable. The presence of thermal instability in the system is mainly visible as the tendency of the gas to avoid the relevant temperature range, as it quickly evolves towards either colder or warmer phases. Nevertheless, the formation of dense structures for both cooling functions appears to be dominated by expanding and colliding supernova remnants, rather than by the thermal instability. As we need to include a finite thermal conduction coefficient to resolve the thermal instability with our uniform grid, we also explore the effects of changing Prandtl number on the system. The purely divergent SN forcing is found to produce significant amounts of vorticity. The relative vorticity is around 0.6-0.7 for the highest Prandtl numbers explored, Pr=40, and is observed to diminish almost by a factor of two for the lowest Prandtl numbers studied, Pr=1. Rotation laws with angular velocity decreasing or increasing outwards are investigated, enabling us to separate the contributions to kinetic helicity due to rotation and shear. When angular velocity decreases outwards, these two contributions partly cancel each other, resulting in a smaller net helicity.

A single open magnetic flux tube spanning the solar photosphere (solar radius R) and the lower corona (R + 10 Mm) is modelled in magnetohydrostatic equilibrium within a realistic stratified atmosphere subject to solar gravity. Such flux tubes are observed to remain relatively stable for up to a day or more, and it is our aim to apply the model as the background condition for numerical studies of energy transport mechanisms from the surface to the corona. We solve analytically an axially symmetric 3D structure for the model, with magnetic field strength, plasma density, pressure and temperature all consistent with observational and theoretical estimates. The self similar construction ensures the magnetic field is divergence free. The equation of pressure balance for this particular set of flux tubes can be integrated analytically to find the pressure and density corrections required to preserve the magnetohydrostatic equilibrium. The model includes a number of free parameters, which makes the solution applicable to a variety of other physical problems and it may therefore be of more general interest.

We simulate the multi-phase interstellar medium randomly heated and stirred by supernovae, with gravity, differential rotation and other parameters of the solar neighbourhood. Here we describe in detail both numerical and physical aspects of the model, including injection of thermal and kinetic energy by SN explosions, radiative cooling, photoelectric heating and various transport processes. With 3D domain extending 1 kpc^2 horizontally and 2 kpc vertically, the model routinely spans gas number densities 10^-5 - 10^2 cm^-3, temperatures 10-10^8 K, local velocities up to 10^3 km s^-1 (with Mach number up to 25). The thermal structure of the modelled ISM is classified by inspection of the joint probability density of the gas number density and temperature. We confirm that most of the complexity can be captured in terms of just three phases, separated by temperature borderlines at about 10^3 K and 5x10^5 K. The probability distribution of gas density within each phase is approximately lognormal. We clarify the connection between the fractional volume of a phase and its various proxies, and derive an exact relation between the fractional volume and the filling factors defined in terms of the volume and probabilistic averages. These results are discussed in both observational and computational contexts. The correlation scale of the random flows is calculated from the velocity autocorrelation function; it is of order 100 pc and tends to grow with distance from the mid-plane. We use two distinct parameterizations of radiative cooling to show that the multi-phase structure of the gas is robust, as it does not depend significantly on this choice.

The origin and structure of the magnetic fields in the interstellar medium of spiral galaxies is investigated with 3D, non-ideal, compressible MHD simulations, including stratification in the galactic gravity field, differential rotation and radiative cooling. A rectangular domain, 1x1x2 kpc^{3} in size, spans both sides of the galactic mid-plane. Supernova explosions drive transonic turbulence. A seed magnetic field grows exponentially to reach a statistically steady state within 1.6 Gyr. Following Germano (1992) we use volume averaging with a Gaussian kernel to separate magnetic field into a mean field and fluctuations. Such averaging does not satisfy all Reynolds rules, yet allows a formulation of mean-field theory. The mean field thus obtained varies in both space and time. Growth rates differ for the mean-field and fluctuating field and there is clear scale separation between the two elements, whose integral scales are about 0.7 kpc and 0.3 kpc, respectively.

We review progress in the development of physically realistic three dimensional simulated models of the galaxy.We consider the scales from star forming molecular clouds to the full spiral disc. Models are computed using hydrodynamic (HD) or magnetohydrodynamic (MHD) equations and may include cosmic ray or tracer particles. The range of dynamical scales between the full galaxy structure and the turbulent scales of supernova (SN) explosions and even cloud collapse to form stars, make it impossible with current computing tools and resources to resolve all of these in one model. We therefore consider a hierarchy of models and how they can be related to enhance our understanding of the complete galaxy.