Search

Small-scale dynamos play important roles in modern astrophysics, especially on galactic and extragalactic scales. Owing to dynamo action, purely hydrodynamic Kolmogorov turbulence hardly exists and is often replaced by hydromagnetic turbulence. Understanding the size of dissipative magnetic structures is important in estimating the time-scale of galactic scintillation and other observational and theoretical aspects of interstellar and intergalactic small-scale dynamos. Here we show that, during the kinematic phase of the small-scale dynamo, the cutoff wavenumber of the magnetic energy spectra scales as expected for large magnetic Prandtl numbers, but continues in the same way also for moderately small values – contrary to what is expected. For a critical magnetic Prandtl number of about 0.3, the dissipative and resistive cutoffs are found to occur at the same wavenumber. In the non-linearly saturated regime, the critical magnetic Prandtl number becomes unity. The cutoff scale now has a shallower scaling with magnetic Prandtl number below a value of about three, and a steeper one otherwise compared to the kinematic regime. [ABSTRACT FROM AUTHOR]

The magnetic fields in galaxy clusters and probably also in the interstellar medium are believed to be generated by a small-scale dynamo. Theoretically, during its kinematic stage, it is characterized by a Kazantsev spectrum, which peaks at the resistive scale. It is only slightly shallower than the Saffman spectrum that is expected for random and causally connected magnetic fields. Causally disconnected fields have the even steeper Batchelor spectrum. Here, we show that all three spectra are present in the small-scale dynamo. During the kinematic stage, the Batchelor spectrum occurs on scales larger than the energy-carrying scale of the turbulence, and the Kazantsev spectrum on smaller scales within the inertial range of the turbulence – even for a magnetic Prandtl number of unity. In the saturated state, the dynamo develops a Saffman spectrum on large scales, suggestive of the build-up of long-range correlations. At large magnetic Prandtl numbers, elongated structures are seen in synthetic synchrotron emission maps showing the parity-even E polarization. We also observe a significant excess in the E polarization over the parity-odd B polarization at subresistive scales, and a deficiency at larger scales. This finding is at odds with the observed excess in the Galactic microwave foreground emission, which is believed to be associated with larger scales. The E and B polarizations may be highly non-Gaussian and skewed in the kinematic regime of the dynamo. For dust emission, however, the polarized emission is always nearly Gaussian, and the excess in the E polarization is much weaker. [ABSTRACT FROM AUTHOR]

In many astrophysical environments, self-gravity can generate kinetic energy, which, in principle, is available for driving dynamo action. Using direct numerical simulations, we show that in unstirred self-gravitating subsonic turbulence with helicity and a magnetic Prandtl number of unity, there is a critical magnetic Reynolds number of about 25 above which the work done against the Lorentz force exceeds the Ohmic dissipation. The collapse itself drives predominantly irrotational motions that cannot be responsible for dynamo action. We find that, with a weak magnetic field, one-third of the work done by the gravitational force goes into compressional heating and the remaining two-thirds go first into kinetic energy of the turbulence before a fraction of it is converted further into magnetic and finally thermal energies. Close to the collapse, however, these fractions change toward 1/4 and 3/4 for compressional heating and kinetic energy, respectively. When the magnetic field is strong, the compressional heating fraction is unchanged. Out of the remaining kinetic energy, one quarter goes directly into magnetic energy via work against the Lorentz force. The fraction of vortical motions diminishes in favour of compressive motions that are almost exclusively driven by the Jeans instability. For an initially uniform magnetic field, field amplification at scales larger than those of the initial turbulence are driven by tangling. [ABSTRACT FROM AUTHOR]

We study the prospects of detecting magnetic helicity in galaxies by observing the dust polarization of the edge-on galaxy NGC 891. Our numerical results of mean-field dynamo calculations show that there should be a large-scale component of the rotationally invariant parity-odd B polarization that we predict to be negative in the first and third quadrants, and positive in the second and fourth quadrants. The large-scale parity-even E polarization is predicted to be negative near the axis and positive further away in the outskirts. These properties are shown to be mostly a consequence of the magnetic field being azimuthal and the polarized intensity being maximum at the centre of the galaxy and are not a signature of magnetic helicity. [ABSTRACT FROM AUTHOR]

We study the effects of ambipolar diffusion (AD) on hydromagnetic turbulence. We consider the regime of large magnetic Prandtl number, relevant to the interstellar medium. In most of the cases, we use the single-fluid approximation where the drift velocity between charged and neutral particles is proportional to the Lorentz force. In two cases we also compare with the corresponding two-fluid model, where ionization and recombination are included in the continuity and momentum equations for the neutral and charged species. The magnetic field properties are found to be well represented by the single-fluid approximation. We quantify the effects of AD on total and spectral kinetic and magnetic energies, the ohmic and AD dissipation rates, the statistics of the magnetic field, the current density, and the linear polarization as measured by the rotationally invariant E and B mode polarizations. We show that the kurtosis of the magnetic field decreases with increasing AD. The E mode polarization changes its skewness from positive values for small AD to negative ones for large AD. Even when AD is weak, changes in AD have a marked effect on the skewness and kurtosis of E , and only a weak effect on those of B. These results open the possibility of employing E and B mode polarizations as diagnostic tools for characterizing turbulent properties of the interstellar medium. [ABSTRACT FROM AUTHOR]

Recent direct numerical simulations (DNS) of large-scale turbulent dynamos in strongly stratified layers have resulted in surprisingly sharp bipolar structures at the surface. Here, we present new DNS of helically and non-helically forced turbulence with and without rotation and compare with corresponding mean-field simulations (MFS) to show that these structures are a generic outcome of a broader class of dynamos in density-stratified layers. The MFS agree qualitatively with the DNS, but the period of oscillations tends to be longer in the DNS. In both DNS and MFS, the sharp structures are produced by converging flows at the surface and might be driven in non-linear stage of evolution by the Lorentz force associated with the large-scale dynamo-driven magnetic field if the dynamo number is at least 2.5 times supercritical. [ABSTRACT FROM AUTHOR]

We use high resolution direct numerical simulations (DNS) to show that helical turbulence can generate significant large-scale fields even in the presence of strong small-scale dynamo action. During the kinematic stage, the unified large/small-scale dynamo grows fields with a shapeinvariant eigenfunction, with most power peaked at small scales or large k, as in Subramanian & Brandenburg. Nevertheless, the large-scale field can be clearly detected as an excess power at small k in the negatively polarized component of the energy spectrum for a forcing with positively polarized waves. Its strength B, relative to the total rms field Brms, decreases with increasing magnetic Reynolds number, ReM. However, as the Lorentz force becomes important, the field generated by the unified dynamo orders itself by saturating on successively larger scales. The magnetic integral scale for the positively polarized waves, characterizing the smallscale field, increases significantly from the kinematic stage to saturation. This implies that the small-scale field becomes as coherent as possible for a given forcing scale, which averts the ReM-dependent quenching of B/Brms. These results are obtained for 1024³ DNS with magnetic Prandtl numbers of PrM = 0.1 and 10. For PrM = 0.1, B/Brms grows from about 0.04 to about 0.4 at saturation, aided in the final stages by helicity dissipation. For PrM = 10, B/Brms grows from much less than 0.01 to values of the order the 0.2. Our results confirm that there is a unified large/small-scale dynamo in helical turbulence. [ABSTRACT FROM AUTHOR]

With the ultimate aim of using the fundamental or f mode to study helioseismic aspects of turbulence-generated magnetic flux concentrations, we use randomly forced hydromagnetic simulations of a piecewise isothermal layer in two dimensions with reflecting boundaries at top and bottom. We compute numerically diagnostic wavenumber–frequency diagrams of the vertical velocity at the interface between the denser gas below and the less dense gas above. For an Alfvén-to-sound speed ratio of about 0.1, a 5 per cent frequency increase of the f mode can be measured when kxHp = 3–4, where kx is the horizontal wavenumber and Hp is the pressure scaleheight at the surface. Since the solar radius is about 2000 times larger than Hp, the corresponding spherical harmonic degree would be 6000–8000. For weaker fields, a kx-dependent frequency decrease by the turbulent motions becomes dominant. For vertical magnetic fields, the frequency is enhanced for kxHp ≈ 4, but decreased relative to its nonmagnetic value for kxHp ≈ 9. [ABSTRACT FROM AUTHOR]

Using direct numerical simulations (DNS), we verify that in the kinematic regime, a turbulent helical dynamo grows in such a way that the magnetic energy spectrum remains to high-precision shape-invariant, i.e. at each wavenumber k the spectrum grows with the same growth rate. Signatures of large-scale dynamo action can be identified through the excess of magnetic energy at small k, of one of the two oppositely polarized constituents. Also a suitably defined planar average of the magnetic field can be chosen such that its rms value isolates the strength of the mean field. However, these different means of analysis suggest that the strength of the large-scale field diminishes with increasing magnetic Reynolds number ReM like ${\rm Re}_{\rm M}^{-1/2}$ for intermediate values and like ${\rm Re}_{\rm M}^{-3/4}$ for larger ones. Both an analysis from the Kazantsev model including helicity and the DNS show that this arises due to the magnetic energy spectrum still peaking at resistive scales, even when helicity is present. As expected, the amplitude of the large-scale field increases with increasing fractional helicity, enabling us to determine the onset of large-scale dynamo action and distinguishing it from that of the small-scale dynamo. Our DNS show that, contrary to earlier results for smaller scale separation (only 1.5 instead of now 4), the small-scale dynamo can still be excited at magnetic Prandtl numbers of 0.1 and only moderate values of the magnetic Reynolds numbers (∼160). [ABSTRACT FROM AUTHOR]

ABSTRACT We consider mean-field dynamo models with fluctuating α effect, both with and without large-scale shear. The α effect is chosen to be Gaussian white noise with zero mean and a given covariance. In the presence of shear, we show analytically that (in infinitely large domains) the mean-squared magnetic field shows exponential growth. The growth rate of the fastest growing mode is proportional to the shear rate. This result agrees with earlier numerical results of Yousef et al. and the recent analytical treatment by Heinemann, McWilliams & Schekochihin who use a method different from ours. In the absence of shear, an incoherent α2 dynamo may also be possible. We further show by explicit calculation of the growth rate of third- and fourth-order moments of the magnetic field that the probability density function of the mean magnetic field generated by this dynamo is non-Gaussian. [ABSTRACT FROM AUTHOR]

Using two- and three-dimensional hydromagnetic simulations for a range of different flows, including laminar and turbulent ones, it is shown that solutions expressing the field in terms of Euler potentials (EP) are in general incorrect if the EP are evolved with an artificial diffusion term. In three dimensions, standard methods using the magnetic vector potential are found to permit dynamo action when the EP give decaying solutions. With an imposed field, the EP method yields excessive power at small scales. This effect is more exaggerated in the dynamic case, suggesting an unrealistically reduced feedback from the Lorentz force. The EP approach agrees with standard methods only at early times when magnetic diffusivity did not have time to act. It is demonstrated that the usage of EP with even a small artificial magnetic diffusivity does not converge to a proper solution of hydromagnetic turbulence. The source of this disagreement is not connected with magnetic helicity or the three-dimensionality of the magnetic field, but is simply due to the fact that the non-linear representation of the magnetic field in terms of EP that depend on the same coordinates is incompatible with the linear diffusion operator in the induction equation. [ABSTRACT FROM AUTHOR]

The generation of mean magnetic fields is studied for a simple non-helical flow where a net cross-helicity of either sign can emerge. This flow, which is also known as the Archontis flow, is a generalization of the Arnold–Beltrami–Childress flow, but with the cosine terms omitted. The presence of cross-helicity leads to a mean-field dynamo effect that is known as the Yoshizawa effect. Direct numerical simulations of such flows demonstrate the presence of magnetic fields on scales larger than the scale of the flow. Contrary to earlier expectations, the Yoshizawa effect is found to be proportional to the mean magnetic field and can therefore lead to its exponential instead of just linear amplification for magnetic Reynolds numbers that exceed a certain critical value. Unlike α effect dynamos, it is found that the Yoshizawa effect is not notably constrained by the presence of a conservation law. It is argued that this is due to the presence of a forcing term in the momentum equation, which leads to a non-zero correlation with the magnetic field. Finally, the application to energy convergence in solar wind turbulence is discussed. [ABSTRACT FROM AUTHOR]

Using mean-field models with a dynamical quenching formalism, we show that in finite domains magnetic helicity fluxes associated with small-scale magnetic fields are able to alleviate catastrophic quenching. We consider fluxes that result from advection by a mean flow, the turbulent mixing down the gradient of mean small-scale magnetic helicity density or the explicit removal which may be associated with the effects of coronal mass ejections in the Sun. In the absence of shear, all the small-scale magnetic helicity fluxes are found to be equally strong for both large- and small-scale fields. In the presence of shear, there is also an additional magnetic helicity flux associated with the mean field, but this flux does not alleviate catastrophic quenching. Outside the dynamo-active region, there are neither sources nor sinks of magnetic helicity, so in a steady state this flux must be constant. It is shown that unphysical behaviour emerges if the small-scale magnetic helicity flux is forced to vanish within the computational domain. [ABSTRACT FROM AUTHOR]

The turbulent diffusion tensor describing the evolution of the mean concentration of a passive scalar is investigated for non-helically forced turbulence in the presence of rotation or a magnetic field. With rotation, the Coriolis force causes a sideways deflection of the flux of mean concentration. Within the magnetohydrodynamics approximation there is no analogous effect from the magnetic field because the effects on the flow do not depend on the sign of the field. Both rotation and magnetic fields tend to suppress turbulent transport, but this suppression is weaker in the direction along the magnetic field. Turbulent transport along the rotation axis is not strongly affected by rotation, except on shorter length-scales, i.e. when the scale of the variation of the mean field becomes comparable with the scale of the energy-carrying eddies. These results are discussed in the context of anisotropic convective energy transport in the Sun. [ABSTRACT FROM AUTHOR]

In the framework of mean-field electrodynamics the coefficients defining the mean electromotive force in Galloway–Proctor flows are determined. These flows show a two-dimensional pattern and are helical. The pattern wobbles in its plane. Apart from one exception a circularly polarized Galloway–Proctor flow, i.e. a circular motion of the flow pattern is assumed. This corresponds to one of the cases considered recently by Courvoisier, Hughes & Tobias. An analytic theory of the α effect and related effects in this flow is developed within the second-order correlation approximation and a corresponding fourth-order approximation. In the validity range of these approximations there is an α effect but no γ effect, or pumping effect. Numerical results obtained with the test-field method, which are independent of these approximations, confirm the results for α and show that γ is in general non-zero. Both α and γ show a complex dependency on the magnetic Reynolds number and other parameters that define the flow, that is, amplitude and frequency of the circular motion. Some results for the magnetic diffusivity and a related quantity are given, too. Finally, a result for α in the case of a randomly varying linearly polarized Galloway–Proctor flow, without the aforementioned circular motion, is presented. The flows investigated show quite interesting effects. There is, however, no straightforward way to relate these flows to turbulence and to use them for studying properties of the α effect and associated effects under realistic conditions. [ABSTRACT FROM AUTHOR]

Using direct simulations, weakly non-linear theory and non-linear mean-field theory, it is shown that the quenched velocity field of a saturated non-linear dynamo can itself act as a kinematic dynamo. The flow is driven by a forcing function that would produce a Roberts flow in the absence of a magnetic field. This result confirms an analogous finding by Cattaneo & Tobias for the more complicated case of turbulent convection, suggesting that this may be a common property of non-linear dynamos; see also the talk given online at the Kavli Institute for Theoretical Physics ( ). It is argued that this property can be used to test non-linear mean-field dynamo theories. [ABSTRACT FROM AUTHOR]

Using numerical simulations at moderate magnetic Reynolds numbers up to 220, it is shown that in the kinematic regime, isotropic helical turbulence leads to an α-effect and a turbulent diffusivity whose values are independent of the magnetic Reynolds number, Rm, provided Rm exceeds unity. These turbulent coefficients are also consistent with expectations from the first-order smoothing approximation. For small values of Rm, α and turbulent diffusivity are proportional to Rm. Over finite time-intervals, meaningful values of α and turbulent diffusivity can be obtained even when there is small-scale dynamo action that produces strong magnetic fluctuations. This suggests that the fields generated by the small-scale dynamo do not make a correlated contribution to the mean electromotive force. [ABSTRACT FROM AUTHOR]

The backreaction of the Lorentz force on the α-effect is studied in the limit of small magnetic and fluid Reynolds numbers, using the first-order smoothing approximation (FOSA) to solve both the induction and momentum equations. Both steady and time-dependent forcings are considered. In the low Reynolds number limit, the velocity and magnetic fields can be expressed explicitly in terms of the forcing function. The non-linear α-effect is then shown to be expressible in several equivalent forms in agreement with formalisms that are used in various closure schemes. On one hand, one can express α completely in terms of the helical properties of the velocity field as in traditional FOSA, or, alternatively, as the sum of two terms, a so-called kinetic α-effect and an oppositely signed term proportional to the helical part of the small-scale magnetic field. These results hold for both steady and time-dependent forcing at arbitrary strength of the mean field. In addition, the τ-approximation is considered in the limit of small fluid and magnetic Reynolds numbers. In this limit, the τ closure term is absent and the viscous and resistive terms must be fully included. The underlying equations are then identical to those used under FOSA, but they reveal interesting differences between the steady and time-dependent forcing. For steady forcing, the correlation between the forcing function and the small-scale magnetic field turns out to contribute in a crucial manner to determine the net α-effect. However for delta-correlated time-dependent forcing, this force–field correlation vanishes, enabling one to write α exactly as the sum of kinetic and magnetic α-effects, similar to what one obtains also in the large Reynolds number regime in the τ-approximation closure hypothesis. In the limit of strong imposed fields, B0, we find for delta-correlated forcing, in contrast to the well-known behaviour for the case of a steady forcing. The analysis presented here is also shown to be in agreement with numerical simulations of steady as well as random helical flows. [ABSTRACT FROM AUTHOR]

The macroscopic behaviour of cosmic rays in turbulent magnetic fields is discussed. An implementation of anisotropic diffusion of cosmic rays with respect to the magnetic field in a non-conservative, high-order, finite-difference magnetohydrodynamic code is discussed. It is shown that the standard implementation fails near singular -points of the magnetic field, which are common if the field is random. A modification to the diffusion model for cosmic rays is described and the resulting telegraph equation (implemented by solving a dynamic equation for the diffusive flux of cosmic rays) is used; it is argued that this modification may better describe the physics of cosmic ray diffusion. The present model reproduces several processes important for the propagation and local confinement of cosmic rays, including spreading perpendicular to the local large-scale magnetic field, controlled by the random-to-total magnetic field ratio, and the balance between cosmic ray pressure and magnetic tension. Cosmic ray diffusion is discussed in the context of a random magnetic field produced by turbulent dynamo action. It is argued that energy equipartition between cosmic rays and other constituents of the interstellar medium does not necessarily imply that cosmic rays play a significant role in the balance of forces. [ABSTRACT FROM AUTHOR]

In an attempt to determine the outer scale of turbulence driven by localized sources, such as supernova explosions in the interstellar medium, we consider a forcing function given by the gradient of Gaussian profiles localized at random positions. Different coherence times of the forcing function are considered. In order to isolate the effects specific to the nature of the forcing function, we consider the case of a polytropic equation of state and restrict ourselves to forcing amplitudes such that the flow remains subsonic. When the coherence time is short, the outer scale agrees with the half-width of the Gaussian. Longer coherence times can cause extra power at large scales, but this would not yield power-law behaviour at scales larger than that of the expansion waves. At scales smaller than the scale of the expansion waves the spectrum is close to power law with a spectral exponent of −2. The resulting flow is virtually free of vorticity. Viscous driving of vorticity turns out to be weak and self-amplification through the non-linear term is found to be insignificant. No evidence for small-scale dynamo action is found in cases where the magnetic induction equation is solved simultaneously with the other equations. [ABSTRACT FROM AUTHOR]

The effect of compressibility on the onset of non-helical turbulent dynamo action is investigated using both direct simulations as well as simulations with shock-capturing viscosities, keeping, however, the regular magnetic diffusivity. It is found that the critical magnetic Reynolds number increases from about 35 in the subsonic regime to about 70 in the supersonic regime. Although the shock structures are sharper in the high-resolution direct simulations compared with the low-resolution shock-capturing simulations, the magnetic field looks roughly similar in both cases and does not show any shock structures. Similarly, the onset of dynamo action is not significantly affected by the shock-capturing viscosity. [ABSTRACT FROM AUTHOR]