Large‐scale magnetic flux concentrations from turbulent stresses

In this study we provide the first numerical demonstration of the effects of turbulence on the mean Lorentz force and the resulting formation of large‐scale magnetic structures. Using three‐dimensional direct numerical simulations (DNS) of forced turbulence we show that an imposed mean magnetic field leads to a decrease of the turbulent hydromagnetic pressure and tension. This phenomenon is quantified by determining the relevant functions that relate the sum of the turbulent Reynolds and Maxwell stresses with the Maxwell stress of the mean magnetic field. Using such a parameterization, we show by means of two‐dimensional and three‐dimensional mean‐field numerical modelling that an isentropic density stratified layer becomes unstable in the presence of a uniform imposed magnetic field. This large‐scale instability results in the formation of loop‐like magnetic structures which are concentrated at the top of the stratified layer. In three dimensions these structures resemble the appearance of bipolar magnetic regions in the Sun. The results of DNS and mean‐field numerical modelling are in good agreement with theoretical predictions. We discuss our model in the context of a distributed solar dynamo where active regions and sunspots might be rather shallow phenomena (© 2010 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

The dual role of shear in large‐scale dynamos

The role of shear in alleviating catastrophic quenching by shedding small‐scale magnetic helicity through fluxes along contours of constant shear is discussed. The level of quenching of the dynamo effect depends on the quenched value of the turbulent magnetic diffusivity. Earlier estimates that might have suffered from the force‐free degeneracy of Beltrami fields are now confirmed for shear flows where this degeneracy is lifted. For a dynamo that is saturated near equipartition field strength those estimates result in a 5‐fold decrease of the magnetic diffusivity as the magnetic Reynolds number based on the wavenumber of the energy‐carrying eddies is increased from 2 to 600. Finally, the role of shear in driving turbulence and large‐scale fields by the magneto‐rotational instability is emphasized. New simulations are presented and the 3π /4 phase shift between poloidal and toroidal fields is confirmed. It is suggested that this phase shift might be a useful diagnostic tool in identifying mean‐field dynamo action in simulations and to distinguish this from other scenarios invoking magnetic buoyancy as a means to explain migration away from the midplane. (© 2008 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Do mean‐field dynamos in nonrotating turbulent shear‐flows exist?

A plane‐shear flow in a fluid with forced turbulence is considered. If the fluid is electrically‐conducting then a mean electromotive force (EMF) results even without basic rotation and the magnetic diffusivity becomes a highly anisotropic tensor. It is checked whether in this case self‐excitation of a large‐scale magnetic field is possible (so‐called W̄ × J̄ ‐dynamo) and the answer is NO. The calculations reveal the cross‐stream components of the EMF perpendicular to the mean current having the wrong signs, at least for small magnetic Prandtl numbers. After our results numerical simulations with magnetic Prandtl number of about unity have only a restricted meaning as the Prandtl number dependence of the diffusivity tensor is rather strong. If, on the other hand, the turbulence field is strati.ed in the vertical direction then a dynamo‐active α ‐effect is produced. The critical magnetic Reynolds number for such a self‐excitation in a simple shear flow is slightly above 10 like for the other – but much more complicated – flow patterns used in existing dynamo experiments with liquid sodium or gallium. (© 2006 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)     

Toroidal field instability and eddy viscosity in Taylor‐Couette flows

Toroidal magnetic fields subject to the Tayler instability can transport angular momentum. We show that the Maxwell and Reynolds stress of the nonaxisymmetric field pattern depend linearly on the shear in the cylindrical gap geometry. Resulting angular momentum transport also scales linear with shear. It is directed outwards for astrophysical relevant flows and directed inwards for superrotating flows with dΩ/dR > 0. We define an eddy viscosity based on the linear relation between shear and angular momentum transport and show that its maximum for given Prandtl and Hartmann number depends linear on the magnetic Reynolds number Rm. For Rm ≃ 1000 the eddy viscosity is of the size of 30 in units of the microscopic value. (© 2008 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Equatorial magnetic helicity flux in simulations with different gauges

We use direct numerical simulations of forced MHD turbulence with a forcing function that produces two different signs of kinetic helicity in the upper and lower parts of the domain. We show that the mean flux of magnetic helicity from the small‐scale field between the two parts of the domain can be described by a Fickian diffusion law with a diffusion coefficient that is approximately independent of the magnetic Reynolds number and about one third of the estimated turbulent magnetic diffusivity. The data suggest that the turbulent diffusive magnetic helicity flux can only be expected to alleviate catastrophic quenching at Reynolds numbers of more than several thousands. We further calculate the magnetic helicity density and its flux in the domain for three different gauges. We consider the Weyl gauge, in which the electrostatic potential vanishes, the pseudo‐Lorenz gauge, where the speed of light is replaced by the sound speed, and the ‘resistive gauge’ in which the Laplacian of the magnetic vector potential acts as a resistive term. We find that, in the statistically steady state, the time‐averaged magnetic helicity density and the magnetic helicity flux are the same in all three gauges (© 2010 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Strong mean field dynamos require supercritical helicity fluxes

Several one and two dimensional mean field models are analyzed where the effects of current helicity fluxes and boundaries are included within the framework of the dynamical quenching model. In contrast to the case with periodic boundary conditions, the final saturation energy of the mean field decreases inversely proportional to the magnetic Reynolds number. If a nondimensional scaling factor in the current helicity flux exceeds a certain critical value, the dynamo can operate even without kinetic helicity, i.e. it is based only on shear and current helicity fluxes, as first suggested by Vishniac & Cho (2001, ApJ 550, 752). Only above this threshold is the current helicity flux also able to alleviate catastrophic quenching. The fact that certain turbulence simulations have now shown apparently non‐resistively limited mean field saturation amplitudes may be suggestive of the current helicity flux having exceeded this critical value. Even below this critical value the field still reaches appreciable strength at the end of the kinematic phase, which is in qualitative agreement with dynamos in periodic domains. However, for large magnetic Reynolds numbers the field undergoes subsequent variations on a resistive time scale when, for long periods, the field can be extremely weak. (© 2005 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Can turbulent plane‐shear flows support large‐scale dynamos?

Turbulent plane‐shear flow is found to show same basic effects of mean‐fieldMHD as rotating turbulence. In particular, the mean electromotive force (EMF) includes highly anisotropic turbulent diffusion and alpha‐effect. Only magnetic diffusion remains for spatially‐uniform turbulence. The question is addressed whether in this case a self‐excitation of a magnetic field by so‐called sher‐current dynamo is possible and the quasilinear theory provides a negative answer. The streamaligned component of the EMF has the sign opposite to that required for dynamo. If, however, the turbulence is not uniform across the flow direction then a dynamo‐active α ‐effect emerges. The critical magnetic Reynolds number for the alpha‐shear dynamo is estimated to be slightly above ten. Possibilities for cross‐checking theoretical predictions with MHD experiments are discussed. (© 2008 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Background magnetic fields during last three cycles of solar activity

This paper describes our studies of evolution of the solar magnetic field with different sign and field strength in the range from –100 G to 100 G. The structure and evolution of large‐scale magnetic fields on the Sun during the last 3 cycles of solar activity is investigated using magnetograph data from the Kitt Peak Solar Observatory. This analysis reveals two groups of the large‐scale magnetic fields evolving differently during the cycles. The first group is represented by relatively weak background fields, and is best observed in the range of 3–10 Gauss. The second group is represented by stronger fields of 75–100 Gauss. The spatial and temporal properties of these groups are described and compared with the total magnetic flux. It is shown that the anomalous behaviour of the total flux during the last cycle can be found only in the second group. (© 2008 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

The role of the Yoshizawa effect in the Archontis dynamo

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.     

Dissipation in dynamos at low and high magnetic Prandtl numbers

Using simulations of helically driven turbulence, it is shown that the ratio of kinetic to magnetic energy dissipation scales with the magnetic Prandtl number in power law fashion with an exponent of approximately 0.6. Over six orders of magnitude in the magnetic Prandtl number the magnetic field is found tobe sustained by large‐scale dynamo action of alphasquared type. This work extends a similar finding for small magnetic Prandtl numbers to the regime of large magnetic Prandtl numbers. At large magnetic Prandtl numbers, most of the energy is dissipated viscously, lowering thus the amount of magnetic energy dissipation, which means that simulations can be performed at magnetic Reynolds numbers that are large compared to the usual limits imposed by a given resolution. This is analogous to an earlier finding that at small magnetic Prandtl numbers, most of the energy is dissipated resistively, lowering the amount of kinetic energy dissipation, so simulations can then be performed at much larger fluid Reynolds numbers than otherwise. The decrease in magnetic energy dissipation at large magnetic Prandtl numbers is discussed in the context of underluminous accretion found in some quasars (© 2011 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Influence of Ohmic diffusion on the excitation and dynamics of MRI

In this paper we make an effort to understand the interaction of turbulence generated by the magnetorotational instability (MRI) with turbulence from other sources, such as supernova explosions (SNe) in galactic disks. First we perform a linear stability analysis (LSA) of non‐ideal MRI to derive the limiting value of Ohmic diffusion that is needed to inhibit the growth of the instability for different types of rotation laws. With the help of a simple analytical expression derived under first‐order smoothing approximation (FOSA), an estimate of the limiting turbulence level and hence the turbulent diffusion needed to damp the MRI is derived. Secondly, we perform numerical simulations in local cubes of isothermal nonstratified gas with external forcing of varying strength to see whether the linear result holds for more complex systems. Purely hydrodynamic calculations with forcing, rotation and shear are made for reference purposes, and as expected, non‐zero Reynolds stresses are found. In the magnetohydrodynamic calculations, therefore, the total stresses generated are a sum of the forcing and MRI contributions. To separate these contributions, we perform reference runs with MRI‐stable shear profiles (angular velocity increasing outwards), which suggest that the MRI‐generated stresses indeed become strongly suppressed as function of the forcing. The Maxwell to Reynolds stress ratio is observed to decrease by an order of magnitude as the turbulence level due to external forcing exceeds the predicted limiting value, which we interpret as a sign of MRI suppression. Finally, we apply these results to estimate the limiting radius inside of which the SN activity can suppress the MRI, arriving at a value of 14 kpc (© 2010 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Subcritical dynamos in shear flows

Identifying generic physical mechanisms responsible for the generation of magnetic fields and turbulence in differentially rotating flows is fundamental to understand the dynamics of astrophysical objects such as accretion disks and stars. In this paper, we discuss the concept of subcritical dynamo action and its hydrodynamic analogue exemplified by the process of nonlinear transition to turbulence in non‐rotating wall‐bounded shear flows. To illustrate this idea, we describe some recent results on nonlinear hydrodynamic transition to turbulence and nonlinear dynamo action in rotating shear flows pertaining to the problem of turbulent angular momentum transport in accretion disks. We argue that this concept is very generic and should be applicable to many astrophysical problems involving a shear flow and non‐axisymmetric instabilities of shearinduced axisymmetric toroidal velocity or magnetic fields, such as Kelvin‐Helmholtz, magnetorotational, Tayler or global magnetoshear instabilities. In the light of several recent numerical results, we finally suggest that, similarly to a standard linear instability, subcritical MHD dynamo processes in high‐Reynolds number shear flows could act as a large‐scale driving mechanism of turbulent flows that would in turn generate an independent small‐scale dynamo. (© 2008 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Radial mixing in the outer Milky Way disc caused by an orbiting satellite

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.     

Stability of the solar tachocline with magnetic fields

The stability of magnetic fields in the solar tachocline is investigated. We present stability limits for higher azimuthal wave numbers and results on the dependence of the stability on the location of toroidal magnetic fields in latitude. While the dependence of the wave number with the largest growth rate on the magnetic field strength and the magnetic Prandtl number is small, the dependence on the magnetic Reynolds number Rm indicates that lowest azimuthal modes are excited for very high Rm. Upon varying the latitudinal position of the magnetic field belts, we find slightly lower stability limits for high latitudes, and very large stability limits at latitudes below 10°, with little dependence on latitude in between. An increase of the maximum possible field was achieved by adding a poloidal field. The upper limit for the toroidal field which can be stored in the radiative tachocline is then 1000 G, compared to about 100 G for a purely toroidal field as was found in an earlier work. (© 2007 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Turbulent viscosity and Λ-effect from numerical turbulence models

Homogeneous anisotropic turbulence simulations are used to determine off-diagonal components of the Reynolds stress tensor and its parameterization in terms of turbulent viscosity and Λ-effect. The turbulence is forced in an anisotropic fashion by enhancing the strength of the forcing in the vertical direction. The Coriolis force is included with a rotation axis inclined relative to the vertical direction. The system studied here is significantly simpler than that of turbulent stratified convection which has often been used to study Reynolds stresses. Certain puzzling features of the results for convection, such as sign changes or highly concentrated latitude distributions, are not present in the simpler system considered here. (© 2007 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Properties of the negative effective magnetic pressure instability

As was demonstrated in earlier studies, turbulence can result in a negative contribution to the effective mean magnetic pressure, which, in turn, can cause a large‐scale instability. In this study, hydromagnetic mean‐field modelling is performed for an isothermally stratified layer in the presence of a horizontal magnetic field. The negative effective magnetic pressure instability (NEMPI) is comprehensively investigated. It is shown that, if the effect of turbulence on the mean magnetic tension force vanishes, which is consistent with results from direct numerical simulations of forced turbulence, the fastest growing eigenmodes of NEMPI are two‐dimensional. The growth rate is found to depend on a parameter β_* characterizing the turbulent contribution of the effective mean magnetic pressure for moderately strong mean magnetic fields. A fit formula is proposed that gives the growth rate as a function of turbulent kinematic viscosity, turbulent magnetic diffusivity, the density scale height, and the parameter β_*. The strength of the imposed magnetic field does not explicitly enter provided the location of the vertical boundaries are chosen such that the maximum of the eigenmode of NEMPI fits into the domain. The formation of sunspots and solar active regions is discussed as possible applications of NEMPI (© 2011 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Linear stability analysis of the Hall magnetorotational instability in a spherical domain

We investigate the stability of the Hall‐MHD system and determine its importance for neutron stars at their birth, when they still consist of differentially rotating plasma permeated by extremely strong magnetic fields. We solve the linearised HallMHD equations in a spherical shell threaded by a homogeneous magnetic field. With the fluid/flow coupling and the Hall effect included, the magnetorotational instability and the Hall effect are both acting together. Results differ for magnetic fields aligned with the rotation axis and anti‐parallel magnetic fields. For a positive alignment of the magnetic field the instability grows on a rotational time‐scale for any sufficiently large magnetic Reynolds number. Even the magnetic fields which are stable against the MRI due to the magnetic diffusion are now susceptible to the shear‐Hall instability. In contrast, the negative alignment places strong restrictions on the growth and the magnitude of the fields, hindering the effectiveness of the Hall‐MRI. While non‐axisymmetric modes of the MRI can be suppressed by strong enough rotation, there is no such restriction when the Hall effect is present. The implications for the magnitude and the topology of the magnetic field of a young neutron star may be significant (© 2011 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Probing Magnetic Fields with Square Kilometre Array and its Precursors

Origin of magnetic fields, its structure and effects on dynamical processes in stars to galaxies are not well understood. Lack of a direct probe has remained a problem for its study. The first phase of Square Kilometre Array (SKA-I), will have almost an order of magnitude higher sensitivity than the best existing radio telescope at GHz frequencies. In this contribution, we discuss specific science cases that are of interest to the Indian community concerned with astrophysical turbulence and magnetic fields. The SKA-I will allow observations of a large number of background sources with detectable polarization and measure their Faraday depths (FDs) through the Milky Way, other galaxies and their circum-galactic mediums. This will probe line-of-sight magnetic fields in these objects well and provide field configurations. Detailed comparison of observational data (e.g., pitch angles in spirals) with models which consider various processes giving rise to field amplification and maintenance (e.g., various types of dynamo models) will then be possible. Such observations will also provide the coherence scale of the fields and its random component through RM structure function. Measuring the random component is important to characterize turbulence in the medium. Observations of FDs with redshift will provide important information on magnetic field evolution as a function of redshift. The background sources could also be used to probe magnetic fields and its coherent scale in galaxy clusters and in bridges formed between interacting galaxies. Other than FDs, sensitive observations of synchrotron emission from galaxies will provide complimentary information on their magnetic field strengths in the sky plane. The core shift measurements of AGNs can provide more precise measurements of magnetic field in the sub parsec region near the black hole and its evolution. The low band of SKA-I will also be useful to study circularly polarized emission from Sun and comparing various models of field configurations with observations.     

Scaling laws and coherent structures in the solar wind

The interplanetary medium is characterized by a very high Reynolds number and is pervaded by fluctuations providing information on a wide range of scales, from fractions of second up to the solar rotation period. In the past decade or so, turbulence in the solar wind has been used as a large wind tunnel to investigate scaling laws of turbulent fluctuations and multifractal models. Moreover, new interesting insights in the theory of turbulence have been derived from the point of view which considers a turbulent flow as a complex system, a sort of benchmark for the theory of dynamical systems. Important finding like the lack of a strict self-similarity of the fluctuations with the consequent nonapplicability of strict scale invariance, the strong anisotropy of velocity and magnetic field fluctuations, the clear lack of equipartition between magnetic and kinetic fluctuations all contributed to suggest the idea that interplanetary fluctuations could possibly be due to a mixture of propagating waves and static structures convected by the wind. In this paper we further discuss this point and bring new evidence about the fact that the presence of a background magnetic field introduces not only a symmetry breaking in interplanetary space but also organizes fluctuations about its large scale orientation.