The development of smooth particle magnetohydrodynamic (SPMHD) has significantly improved the simulation of complex astrophysical processes. However, the preservation the solenoidality of the magnetic field is still a severe problem for the MHD. A fo
rmulation of the induction equation with a vector potential would solve the problem. Unfortunately all previous attempts suffered from instabilities. In the present work, we evolve the vector potential in the Coulomb gauge and smooth the derived magnetic field for usage in the momentum equation. With this implementation we could reproduce classical test cases in a stable way. A simple test case demonstrates the possible failure of widely used direct integration of the magnetic field, even with the usage of a divergence cleaning method.
We aim to estimate the contribution of the radial component of the Lorentz force to the gas rotation in several types of galaxies. Using typical parameters for the exponential scale of synchrotron emission and the scale length of HI gas, under the as
sumption of equipartition between the energies of cosmic rays and total magnetic fields, we derive the Lorentz force and compare it to the gravitational force in the radial component of the momentum equation. We distinguish the different contributions between the large-scale and the small-scale turbulent fields by Reynolds averaging. We compare these findings with a dynamical dynamo model. We find a possible reduction of circular gas velocity in the very outer parts and an increase inside a radius of four times the synchrotron scale length. Sufficiently localized radial reversals of the magnetic field may cause characteristic modulations in the gas rotation curve with typical amplitudes of 10-20 km/s. It is unlikely that the magnetic field contributes to the flat rotation in the outer parts of galaxies. If anything, it will emph{impede} the gravitationally supported rotation, demanding for an even higher halo mass to explain the observed rotation profile. We speculate that this may have consequences for ram pressure stripping and the truncation of the stellar disc.
(Abridged) Based on the rapidly increasing all-sky data of Faraday rotation measures and polarised synchrotron radiation, the Milky Ways magnetic field is now modelled with an unprecedented level of detail and complexity. We aim to complement this he
uristic approach with a physically motivated, quantitative Galactic dynamo model -- a model that moreover allows for the evolution of the system as a whole, instead of just solving the induction equation for a fixed static disc. Building on the framework of mean-field magnetohydrodynamics and extending it to the realm of a hybrid evolution, we perform three-dimensional global simulations of the Galactic disc. Closure coefficients embodying the mean-field dynamo are calibrated against resolved box simulations of supernova-driven interstellar turbulence. The emerging dynamo solutions comprise a mixture of the dominant axisymmetric S0 mode, with even parity, and a subdominant A0 mode, with odd parity. Notably, such a superposition of modes creates a strong localised vertical field on one side of the Galactic disc. We moreover find significant radial pitch angles, which decay with radius -- explained by flaring of the disc. In accordance with previous work, magnetic instabilities appear to be restricted to the less-stirred outer Galactic disc. Their main effect is to create strong fields at large radii such that the radial scale length of the magnetic field increases from 4 kpc (for the case of a mean-field dynamo alone) to about 10 kpc in the hybrid models. There remain aspects (e.g., spiral arms, X-shaped halo fields, fluctuating fields) that are not captured by the current model and that will require further development towards a fully dynamical evolution. Nevertheless, the work presented demonstrates that a hybrid modelling of the Galactic dynamo is feasible and can serve as a foundation for future efforts.
The emergence of large-scale magnetic fields observed in the diffuse interstellar medium is explained by a turbulent dynamo. The underlying transport coefficients have previously been extracted from numerical simulations. So far, this was restricted
to the kinematic regime, but we aim to extend our analysis into the realm of dynamically important fields. This marks an important step on which derived mean-field models rely to explain observed equipartition strength fields. As in previous work, we diagnose turbulent transport coefficients by means of the test-field method. We derive quenching functions for the dynamo {alpha} effect, diamagnetic pumping, and turbulent diffusivity, which are compared with theoretical expectations. At late times, we observe the suppression of the vertical wind. Because this potentially affects the removal of small-scale magnetic helicity, new concerns arise about circumventing constraints imposed by the conservation of magnetic helicity at high magnetic Reynolds numbers. While present results cannot safely rule out this possibility, the issue only becomes important at late stages and is absent when the dynamo is quenched by the wind itself.
Magnetic field amplification by a fast dynamo is seen in local box simulations of SN-driven ISM turbulence, where the self-consistent emergence of large-scale fields agrees very well with its mean-field description. We accordingly derive scaling laws
of the turbulent transport coef- ficients in dependence of the SN rate, density and rotation. These provide the input for global simulations of regular magnetic fields in galaxies within a mean-field MHD framework. Using a Kennicutt-Schmidt relation between the star formation (SF) rate and midplane density, we can reduce the number of free parameters in our global models. We consequently present dynamo models for different rotation curves and radial density distributions.
The fractal shape and multi-component nature of the interstellar medium together with its vast range of dynamical scales provides one of the great challenges in theoretical and numerical astrophysics. Here we will review recent progress in the direct
modelling of interstellar hydromagnetic turbulence, focusing on the role of energy injection by supernova explosions. The implications for dynamo theory will be discussed in the context of the mean-field approach. Results obtained with the test field-method are confronted with analytical predictions and estimates from quasilinear theory. The simulation results enforce the classical understanding of a turbulent Galactic dynamo and, more importantly, yield new quantitative insights. The derived scaling relations enable confident global mean-field modelling.
Supernovae are the dominant energy source for driving turbulence within the interstellar plasma. Until recently, their effects on magnetic field amplification in disk galaxies remained a matter of speculation. By means of self-consistent simulations
of supernova-driven turbulence, we find an exponential amplification of the mean magnetic field on timescales of a few hundred million years. The robustness of the observed fast dynamo is checked at different magnetic Reynolds numbers, and we find sustained dynamo action at moderate Rm. This indicates that the mechanism might indeed be of relevance for the real ISM. Sensing the flow via passive tracer fields, we infer that SNe produce a turbulent alpha effect which is consistent with the predictions of quasilinear theory. To lay a foundation for global mean-field models, we aim to explore the scaling of the dynamo tensors with respect to the key parameters of our simulations. Here we give a first account on the variation with the supernova rate.
Within the interstellar medium, supernovae are thought to be the prevailing agents in driving turbulence. Until recently, their effects on magnetic field amplification in disk galaxies remained uncertain. Analytical models based on the uncorrelated-e
nsemble approach predicted that any created field would be expelled from the disk before it could be amplified significantly. By means of direct simulations of supernova-driven turbulence, we demonstrate that this is not the case. Accounting for galactic differential rotation and vertical stratification, we find an exponential amplification of the mean field on timescales of several hundred million years. We especially highlight the importance of rotation in the generation of helicity by showing that a similar mechanism based on Cartesian shear does not lead to a sustained amplification of the mean magnetic field.
Supernovae are known to be the dominant energy source for driving turbulence in the interstellar medium. Yet, their effect on magnetic field amplification in spiral galaxies is still poorly understood. Previous analytical models, based on the evoluti
on of isolated, non-interacting supernova remnants, predicted a dominant vertical pumping that would render dynamo action improbable. In the present work, we address the issue of vertical transport, which is thought to be the key process that inhibits dynamo action in the galactic context. We aim to demonstrate that supernova driving is a powerful mechanism to amplify galactic magnetic fields. We conduct direct numerical simulations in the framework of resistive magnetohydrodynamics. Our local box model of the interstellar medium comprises optically-thin radiative cooling, an external gravitational potential, and background shear. Dynamo coefficients for mean-field models are measured by means of passive test fields. Our simulations show that supernova-driven turbulence in conjunction with shear leads to an exponential amplification of the mean magnetic field. We found turbulent pumping to be directed inward and approximately balanced by a galactic wind.
