Dynamo saturation in direct simulations of the multi‐phase turbulent interstellar medium

The ordered magnetic field observed via polarised synchrotron emission in nearby disc galaxies can be explained by a mean‐field dynamo operating in the diffuse interstellar medium (ISM). Additionally, vertical‐flux initial conditions are potentially able to influence this dynamo via the occurrence of the magnetorotational instability (MRI). We aim to study the influence of various initial field configurations on the saturated state of the mean‐field dynamo. This is motivated by the observation that different saturation behaviour was previously obtained for different supernova rates. We perform direct numerical simulations (DNS) of three‐dimensional local boxes of the vertically stratified, turbulent interstellar medium, employing shearing‐periodic boundary conditions horizontally. Unlike in our previous work, we also impose a vertical seed magnetic field. We run the simulations until the growth of the magnetic energy becomes negligible. We furthermore perform simulations of equivalent 1D dynamo models, with an algebraic quenching mechanism for the dynamo coefficients. We compare the saturation of the magnetic field in the DNS with the algebraic quenching of a mean‐field dynamo. The final magnetic field strength found in the direct simulation is in excellent agreement with a quenched α) dynamo. For supernova rates representative of the Milky Way, field losses via a Galactic wind are likely responsible for saturation. We conclude that the relative strength of the turbulent and regular magnetic fields in spiral galaxies may depend on the galaxy's star formation rate. We propose that a mean field approach with algebraic quenching may serve as a simple sub‐grid scale model for galaxy evolution simulations including a prescribed feedback from magnetic fields. (© 2015 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Dimensionless measures of turbulent magnetohydrodynamic dissipation rates

The magnetic Reynolds number, R _M, is defined as the product of a characteristic scale and associated flow speed divided by the microphysical magnetic diffusivity. For laminar flows, R _M also approximates the ratio of advective to dissipative terms in the total magnetic energy equation, but for turbulent flows this latter ratio depends on the energy spectra and approaches unity in a steady state. To generalize for flows of arbitrary spectra we define an effective magnetic dissipation number,   R _M,e  , as the ratio of the advection to microphysical dissipation terms in the total magnetic energy equation, incorporating the full spectrum of scales, arbitrary magnetic Prandtl numbers, and distinct pairs of inner and outer scales for magnetic and kinetic spectra. As expected, for a substantial parameter range   R _M,e∼ O (1) ≪ R _M  . We also distinguish   R _M,e  from     where the latter is an effective magnetic Reynolds number for the mean magnetic field equation when a turbulent diffusivity is explicitly imposed as a closure. That   R _M,e  and     approach unity even if   R _M≫ 1  highlights that, just as in hydrodynamic turbulence, energy dissipation of large-scale structures in turbulent flows via a cascade can be much faster than the dissipation of large-scale structures in laminar flows. This illustrates that the rate of energy dissipation by magnetic reconnection is much faster in turbulent flows, and much less sensitive to microphysical reconnection rates compared to laminar flows.     

LOLAS: an optical turbulence profiler in the atmospheric boundary layer with extreme altitude resolution

We report the development and first results of an instrument called Low Layer SCIDAR (Scintillation Detection and Ranging) (LOLAS) which is aimed at the measurement of optical-turbulence profiles in the atmospheric boundary layer with high altitude resolution. The method is based on the Generalized SCIDAR (GS) concept, but unlike the GS instruments which need a 1-m or larger telescope, LOLAS is implemented on a dedicated 40-cm telescope, making it an independent instrument. The system is designed for widely separated double-star targets, which enables the high altitude resolution. Using a 200-arcsec-separation double star, we have obtained turbulence profiles with unprecedented 12-m resolution. The system incorporates necessary novel algorithms for autoguiding, autofocus and image stabilization. The results presented here were obtained at Mauna Kea Observatory. They show LOLAS capabilities but cannot be considered as representative of the site. A forthcoming paper will be devoted to the site characterization. The instrument was built as part of the Ground Layer Turbulence Monitoring Campaign on Mauna Kea for Gemini Observatory.     

Dynamo coefficients from local simulations of the turbulent ISM

Observations in polarized emission reveal the existence of large‐scale coherent magnetic fields in a wide range of spiral galaxies. Radio‐polarization data show that these fields are strongly inclined towards the radial direction, with pitch angles up to 35° and thus cannot be explained by differential rotation alone. Global dynamo models describe the generation of the radial magnetic field from the underlying turbulence via the so called α ‐effect. However, these global models still rely on crude assumptions about the small‐scale turbulence. To overcome these restrictions we perform fully dynamical MHD simulations of interstellar turbulence driven by supernova explosions. From our simulations we extract profiles of the contributing diagonal elements of the dynamo α ‐tensor as functions of galactic height. We also measure the coefficients describing vertical pumping and find that the ratio between these two effects has been overestimated in earlier analytical work, where dynamo action seemed impossible. In contradiction to these models based on isolated remnants we always find the pumping to be directed inward. In addition we observe that depends on whether clustering in terms of superbubbles is taken into account. Finally, we apply a test field method to derive a quantitative measure of the turbulent magnetic diffusivity which we determine to be ∼2 kpckms^–1. (© 2008 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Hopf bifurcations with fluctuating gain

A Hopf bifurcation describes a transition from a steady state to an oscillatory one. This report describes how the presence of white Gaussian noise in the term driving the instability (the gain) affects the character of such a transition. We find that for small noise intensities the effective critical value of the bifurcation parameter increases linearly with the noise intensity. This result is used to quantitatively describe the onset of the magneto‐rotational instability in a laboratory experiment using liquid metal with turbulent background flow. (© 2008 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Nonhelical mean‐field dynamos in a sheared turbulence

Mechanisms of nonhelical large‐scale dynamos (shear‐current dynamo and effect of homogeneous kinetic helicity fluctuations with zero mean) in a homogeneous turbulence with large‐scale shear are discussed. We have found that the shearcurrent dynamo can act even in random flows with small Reynolds numbers. However, in this case mean‐field dynamo requires small magnetic Prandtl numbers (i.e., when Pm < Pm^cr < 1). The threshold in the magnetic Prandtl number, Pm^cr = 0.24, is determined using second order correlation approximation (or first‐order smoothing approximation) for a background random flow with a scale‐dependent viscous correlation time τ_c = (νk 2)^–1 (where ν is the kinematic viscosity of the fluid and k is the wave number). For turbulent flows with large Reynolds numbers shear‐current dynamo occurs for arbitrary magnetic Prandtl numbers. This dynamo effect represents a very generic mechanism for generating large‐scale magnetic fields in a broad class of astrophysical turbulent systems with large‐scale shear. On the other hand, mean‐field dynamo due to homogeneous kinetic helicity fluctuations alone in a sheared turbulence is not realistic for a broad class of astrophysical systems because it requires a very specific random forcing of kinetic helicity fluctuations that contains, e.g., low‐frequency oscillations. (© 2008 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

The inner scale of the plasma turbulence towards PSR J1644−4559

By measuring the decaying shape of the scatter-broadened pulse from the bright distant pulsar PSR J1644−4559, we probe waves scattered at relatively high angles by very small spatial scales in the interstellar plasma, which allows us to test for a wavenumber cutoff in the plasma density spectrum. Under the hypothesis that the density spectrum is due to plasma turbulence, we can thus investigate the (inner) scale at which the turbulence is dissipated. We report observations carried out with the Parkes radio telescope at 660 MHz from which we find strong evidence for an inner scale in the range 70–100 km, assuming an isotropic Kolmogorov spectrum. By identifying the inner scale with the ion inertial scale, we can also estimate the mean electron density of the scattering region to be 5–10 cm⁻³. This is comparable with the electron density of H  ii region G339.1−0.4, which lies in front of the pulsar, and so confirms that this region dominates the scattering. We conclude that the plasma inside the region is characterized by fully developed turbulence with an outer scale in the range 1–20 pc and an inner scale of 70–100 km. The shape of the rising edge of the pulse constrains the distribution of the strongly scattering plasma to be spread over about 20 per cent of the 4.6 kpc path from the pulsar, but with similarly high electron densities in two or more thin layers, their thicknesses can only be 10–20 pc.     

Large-scale galactic turbulence: can self-gravity drive the observed H i velocity dispersions?

Observations of turbulent velocity dispersions in the H  i component of galactic discs show a characteristic floor in galaxies with low star formation rates and within individual galaxies the dispersion profiles decline with radius. We carry out several high-resolution adaptive mesh simulations of gaseous discs embedded within dark matter haloes to explore the roles of cooling, star formation, feedback, shearing motions and baryon fraction in driving turbulent motions. In all simulations the disc slowly cools until gravitational and thermal instabilities give rise to a multiphase medium in which a large population of dense self-gravitating cold clouds are embedded within a warm gaseous phase that forms through shock heating. The diffuse gas is highly turbulent and is an outcome of large-scale driving of global non-axisymmetric modes as well as cloud–cloud tidal interactions and merging. At low star formation rates these processes alone can explain the observed H  i velocity dispersion profiles and the characteristic value of  ∼10 km s⁻¹  observed within a wide range of disc galaxies. Supernovae feedback creates a significant hot gaseous phase and is an important driver of turbulence in galaxies with a star formation rate per unit area  ≳10⁻³ M_⊙ yr⁻¹ kpc⁻²  .     

Decameter Type III-Like Bursts

We report the first observations of Type III-like bursts at frequencies 10 – 30 MHz. More than 1000 such bursts during 2002 – 2004 have been analyzed. The frequency drift of these bursts is several times that of decameter Type III bursts. A typical duration of the Type III-like bursts is 1 – 2 s. These bursts are mainly observed when the source active region is located within a few days from the central meridian. The drift rate of the Type III-like bursts can take a large value by considering the velocity of Type III electrons and the group velocity of generated electromagnetic waves.