Alpha effect and diffusivity in helical turbulence with shear

D. Mitra; P. J. Käpylä; R. Tavakol; A. Brandenburg

doi:10.1051/0004-6361:200810359

Free Access

Issue		A&A Volume 495, Number 1, February III 2009


Page(s)		1 - 8
Section		Astrophysical processes
DOI		https://doi.org/10.1051/0004-6361:200810359
Published online		11 December 2008

Alpha effect and diffusivity in helical turbulence with shear

D. Mitra¹ - P. J. Käpylä² - R. Tavakol¹ - A. Brandenburg³

1 - Astronomy unit, School of Mathematical Sciences, Queen Mary University of London, Mile End Road, London E1 4NS, UK
2 - Observatory, Tähtitorninmäki (PO Box 14), 00014 University of Helsinki, Finland
3 - NORDITA, Roslagstullsbacken 23, 10691 Stockholm, Sweden

Received 10 June 2008 / Accepted 28 November 2008

Abstract
Aims. We study the dependence of turbulent transport coefficients, such as the components of the $\alpha$ tensor ( $\alpha_{ij}$ ) and the turbulent magnetic diffusivity tensor ( $\eta_{ij}$ ), on shear and magnetic Reynolds number in the presence of helical forcing.
Methods. We use three-dimensional direct numerical simulations with periodic boundary conditions and measure the turbulent transport coefficients using the kinematic test field method. In all cases the magnetic Prandtl number is taken as unity.
Results. We find that with increasing shear the diagonal components of $\alpha_{ij}$ quench, whereas those of $\eta_{ij}$ increase. The antisymmetric parts of both tensors increase with increasing shear. We also propose a simple expression for the turbulent pumping velocity (or $\gamma$ effect). This pumping velocity is proportional to the kinetic helicity of the turbulence and the vorticity of the mean flow. For negative helicity, i.e. for a positive trace of $\alpha_{ij}$ , it points in the direction of the mean vorticity, i.e. perpendicular to the plane of the shear flow. Our simulations support this expression for low shear and magnetic Reynolds number. The transport coefficients depend on the wavenumber of the mean flow in a Lorentzian fashion, just as for non-shearing turbulence.

Key words: magnetohydrodynamics (MHD) - hydrodynamics - turbulence - magnetic fields

1 Introduction

Understanding the origin of non-uniform large-scale magnetic fields in stars, galaxies, and accretion discs continues to pose important challenges. Such fields are commonly thought to be the result of dynamo action that converts the kinetic energy of turbulent motions and large-scale shear into magnetic energy. The usual framework for studying such dynamo actions is mean field electrodynamics (e.g. Krause & Rädler 1980). Over the years, however, the applicability of this framework has been questioned (Piddington 1981; Vainshtein & Cattaneo 1992). In particular, an important debate in this connection revolves around the role played by magnetic helicity (Gruzinov & Diamond 1994). Magnetic helicity is conserved for ideal (non-dissipative) magnetohydrodynamics (MHD) and also in the limit of large magnetic Reynolds numbers. This conservation places severe constraints on the growth of the mean magnetic field and may regulate the quenching of the $\alpha$ effect as the magnetic Reynolds number ( $\mbox{Re}_{\rm M}$ ) increases (Brandenburg & Subramanian 2005). Under certain circumstances (no magnetic helicity fluxes), $\alpha$ quenching becomes more severe with $\alpha$ proportional to $\mbox{Re}_{\rm M}^{-1}$ . It has been suggested that this can be alleviated by expelling magnetic helicity through open boundaries (Kleeorin et al. 2000; Blackman & Field 2000), possibly mediated by shear (Subramanian & Brandenburg 2004; Vishniac & Cho 2001; Brandenburg & Sandin 2004). Furthermore, shear itself is an important ingredient in MHD dynamos in solar and stellar settings as, for example, in the solar tachocline. Hence, it is important to understand how shear affects the turbulent transport coefficients, in particular the components of the $\alpha_{ij}$ and $\eta_{ij}$ tensors.

Several studies have recently looked at various aspects of this problem (see, e.g., Rüdiger & Kitchatinov 2006; Rogachevskii & Kleeorin 2003,2004; Brandenburg 2005; Leprovost & Kim 2008; Rädler & Stepanov 2006;Brandenburg et al. 2008a). These works employ different tools and make different assumptions and are often applicable to limited regions of the parameter space. As a consequence, care must be taken in comparing these results. For example, using semi-analytical tools, which treat the nonlinear Lorentz force feedback perturbatively, Leprovost & Kim (2008) have found that shear can reduce $\alpha$ in helically forced turbulence. This is analogous to the ``rotational quenching'' of turbulent transport coefficients with increasing Coriolis or inverse Rossby numbers (Kitchatinov et al. 1994; Pipin et al. 1996). There are also other related cases in which the presence of shear enhances the growth rate of the dynamo. For example; using direct numerical simulations of the MHD equations in the presence of shear and non-helical forcing, Yousef et al. (2008a,b) have found large-scale dynamos whose growth rate increases linearly with shear. Such scaling has also been found for $\alpha$ -shear dynamos where the $\alpha$ effect is due to stratified convection with shear (Käpylä et al. 2008a). Furthermore, using the kinematic test field method (described below), Brandenburg et al. (2008a) have studied the dynamo coefficients in the presence of shear, but in the absence of helicity, and they find that Gaussian fluctuations of $\alpha$ about zero are strong enough to drive an incoherent $\alpha$ -shear dynamo (Vishniac & Brandenburg 1997; Proctor 2007). The significance of the incoherent $\alpha$ effect has also been stressed by Hughes & Proctor (2009), although their system may have also had a net $\alpha$ effect. They dismissed this on the grounds that for an imposed uniform magnetic field $\alpha$ is very small. However, this result disagrees with recent calculations of $\alpha$ using the test field method (Käpylä et al. 2008b).

In this paper we use three-dimensional direct numerical simulations of the kinematic test field equations with helical forcing in order to study the dependence of turbulent transport coefficients on shear and magnetic Reynolds number. In Sect. 2 we give a brief account of our model. Section 3 contains our results and finally we conclude in Sect. 4.

2 The model

We use the test field method (Schrinner et al. 2005,2007) to calculate the turbulent transport coefficients. This method and its modification in the presence of large-scale shear are described in Brandenburg et al. (2008a). Here we just point out the essence of the method and elaborate only on those aspects where our treatment differs from their paper.

In the presence of large-scale shear, the equations of magnetohydrodynamics are treated in the following way. Writing the velocity as the sum $\vec{ U} + \vec{ U}^{S}$ , where the large-scale shear velocity $\vec{U}^{S}=Sx\vec{\hat{y}}$ with a constant shear S, and assuming an isothermal equation of state characterised by the sound speed $c_{\rm s}$ , the momentum equation becomes

$\begin{displaymath}\frac{{\mathcal D} \vec{ U}}{{\mathcal D} t} = - \vec{ U} \... ... s}^2\vec{\nabla} \ln \rho + \vec{ f} + \vec{ F}_{{\rm visc}}. \end{displaymath}$

(1)

Here $\vec{ F}_{{\rm visc}} = \rho^{-1}\vec{\nabla}\cdot(2\rho\nu{\mbox{\boldmath${\sf S}$ } {}})$ is the viscous force, ${\sf S}_{ij}= \frac{1}{2}(U_{i,j}+U_{j,i}) - \frac{1}{3}\delta_{ij}\vec{\nabla} \cdot \vec{ U}$ is the traceless rate of strain tensor (not to be confused with the shear parameter S), $\nu$ is the kinematic viscosity, $\rho$ is the fluid density, $c_{\rm s}$ is the isothermal sound speed and

$\begin{displaymath}\frac{{\cal D} }{{\cal D} t} \equiv \frac{\partial}{\partial t} + Sx\frac{\partial}{\partial y}\cdot \end{displaymath}$

(2)

As our external forcing $\vec{ f}$ we employ helical, white-in-time, random forcing described in Brandenburg (2001). In this paper we consider the purely kinematic problem, so there is no feedback in Eq. (1) due to the Lorentz force. In addition to the momentum equation, we have the continuity equation

$\begin{displaymath}\frac{{\mathcal D} \rho}{{\mathcal D} t} = - \vec{ U} \cdot \vec{\nabla} \ln \rho - \vec{\nabla}\cdot \vec{ U}, \end{displaymath}$

(3)

and the uncurled induction equation in the Weyl gauge,

$\begin{displaymath}\frac{{\mathcal D}\vec{ A}}{{\mathcal D}t}=-SA_y\vec{\hat{x}}+\vec{ U}\times\vec{ B} -\eta\mu_0\vec{ J}. \end{displaymath}$

(4)

Here the magnetic field is $\vec{ B} = \vec{\nabla} \times \vec{ A}$ , the current density is $\vec{ J} = \vec{\nabla} \times \vec{ B}/\mu_0$ , and $\mu_0$ and $\eta$ are the vacuum permeability and the molecular magnetic diffusivity, respectively. In the mean field approach to MHD one usually decomposes the magnetic field (or magnetic vector potential) and velocity into mean (indicated by an overbar) and fluctuating parts respectively

$\begin{displaymath}\vec{A} = \overline{\vec{A}} + \vec{a}, ~~~~ \vec{U} = \overline{\vec{U}} + \vec{u}. \end{displaymath}$

(5)

The equations satisfied by the mean and fluctuating parts of the magnetic vector potential are given by

$\begin{displaymath}\frac{\mathcal{D} {\overline{\vec{ A}}}}{\mathcal{D} t} = -... ...{\mbox{\boldmath${\cal E}$ }}{}-\eta\mu_0{\overline{\vec{J}}}, \end{displaymath}$

(6)

and

$\begin{displaymath}\frac{{\mathcal D}\vec{ a}}{{\mathcal D}t}= -Sa_y\vec{\hat{x}... ...- \overline{\mbox{\boldmath${\cal E}$ }}{}- \eta\mu_0\vec{ j}, \end{displaymath}$

(7)

where $\vec{j} = \vec{ J} - {\overline {\vec{J}}}$ and $\overline{\mbox{\boldmath${\cal E}$ }}{}= \overline{\vec{ u} \times \vec{ b}}$ is the mean turbulent electromotive force. Specification of this second order quantity in terms of the mean field constitutes a closure problem. A common procedure is to expand $\overline{\mbox{\boldmath${\cal E}$ }}{}$ in terms of the mean field $\overline{\vec{B}}$ and its derivatives,

$\begin{displaymath}\overline{\mathcal{E}}_i = \overline{\mathcal{E}}_{0i} + \alpha_{ij} {\overline B_j} + \eta_{ijk} {\overline B_{j,k}} , \end{displaymath}$

(8)

where $\alpha_{ij}$ and $\eta_{ijk}$ are the tensorial turbulent transport coefficients, and $\overline{\mbox{\boldmath${\cal E}$ }}{}_0$ quantifies additional contributions that arise even in the absence of a mean field, owing to small-scale dynamo action, for example. Here and throughout, summation is assumed over repeated indices. In the test field method we take the mean magnetic field to be a given ``test field'' $\overline{\vec{B}}$ , and calculate these tensors by measuring $\overline{\mbox{\boldmath${\cal E}$ }}{}$ by solving Eq. (7) simultaneously with Eqs. (1) and (3), while using Eq. (5) to find $\vec{ u}$ . In order to find all the components of the $\alpha_{ij}$ and $\eta_{ij}$ tensors, and not just projections relevant to the actual fields, one has to use an orthogonal set of test fields and solve Eq. (7) for each of them.

In our simulations we employ averages over x and y directions to define our mean fields, which are therefore functions of z and t only. Thus, all other components of ${\overline B}_{j,k}$ except ${\overline B}_{x,z}$ and ${\overline B}_{y,z}$ vanish, and the $\alpha$ and $\eta$ tensors can be written as rank two tensors with $\eta_{i1} = \eta_{i23}$ and $\eta_{i2} = -\eta_{i13}$ for i,j=(1,2)(Brandenburg et al. 2008a). Now consider, as an example, the following two test fields,

$\displaystyle {\overline B}^{c1}_i = B_0(\cos kz,0,0),$		$\displaystyle \mu_0{\overline J}^{c1}_i = kB_0(0,-\sin kz,0),$	(9)
$\displaystyle {\overline B}^{s1}_i = B_0(\sin kz,0,0),$		$\displaystyle \mu_0{\overline J}^{s1}_i = kB_0(0,+\cos kz,0),$	(10)

which we use to compute the two corresponding mean electromotive forces $\overline{\mbox{\boldmath${\cal E}$ }}{}^{c1}$ and $\overline{\mbox{\boldmath${\cal E}$ }}{}^{s1}$ . The relevant components of the $\alpha$ and $\eta$ tensors are then given by

$\begin{displaymath}\left( \begin{array}{c} \alpha_{i1} \\ - \eta_{i2}k \end{arra... ...1}_i \\ {\overline{\mathcal E}}^{s1}_i \end{array} \right) . \end{displaymath}$

(11)

The rest of the 2+2 components of the $\alpha$ and $\eta$ tensors can be similarly determined by using the test fields,

$\begin{displaymath}{\overline B}^{c2}_i = B_0(0,\cos kz,0), ~~~ \/{\overline B}^{s2}_i = kB_0(0,\sin kz,0) . \end{displaymath}$

(12)

In what follows we denote a test field by ${\overline{\vec{ B}}}^{pq}$ , where p=c,s and q=1,2. A particular small-scale magnetic vector potential that develops in response to the ${\overline{\vec{ B}}}^{pq}$ is denoted by $\vec{ a}^{pq}$ and the corresponding small-scale magnetic field is given by $\vec{ b}^{pq} = \vec{\nabla} \times \vec{ a}^{pq}$ . Note that at large values of $\mbox{Re}_{\rm M}$ there will also be small-scale dynamo action that will lead to spurious time-dependencies of $\overline{\mbox{\boldmath${\cal E}$ }}{}$ . However, since neither the test fields nor $\alpha_{ij}$ or $\eta_{ij}$ depend on time, i.e.

$\begin{displaymath}\overline{\mathcal{E}}_i^{pq}(z,t) = \overline{\mathcal{E}}_{... ...ine B_j}^{pq}(z) + \eta_{ijk} {\overline B_{j,k}}^{pq}(z) ~ , \end{displaymath}$

(13)

such time dependence must be entirely due to $\overline{\mbox{\boldmath${\cal E}$ }}{}_0^{pq}(z,t)$ and can be eliminated by time averaging.

$\begin{figure} \mbox{\includegraphics[width=6cm,clip]{0359fig1a}\includegraphics[width=6cm,clip]{0359fig1b} }\end{figure}$

Figure 1:

Root-mean-square velocity $u_{\rm rms}$ normalized by the speed of sound $c_{\rm s}$ as: a) a function of $\mbox{Re}_{\rm M}$ for different values of $\mbox{Sh}$ : $\mbox{Sh}\sim 0.07 (\ast), 0.2 (\triangle), \mbox{and}\ 0.3 (\Box)$ respectively and b) a function of the shear parameter $\mbox{Sh}$ for different value of $\mbox{Re}_{\rm M}$ : $\mbox{Re}_{\rm M}\sim 1$ ( $\ast$ ), $\mbox{Re}_{\rm M}\sim 20$ ( $\triangle$ ), and $\mbox{Re}_{\rm M}\sim 72$ ( $\square$ ).

Open with DEXTER

In the following we shall, to begin with, use k=k₁, the wavenumber corresponding to the box size, to study the dependence of $\alpha_{ij}$ and $\eta_{ij}$ on shear and magnetic Reynolds number. We shall discuss the dependence of the $\alpha_{ij}$ and $\eta_{ij}$ tensors on k in Sect. 3.3. Note that the usual approach of using uniform applied fields for calculating $\alpha$ (e.g., Courvoisier et al. 2006) corresponds to a special case of the test field method for k=0. However, dynamos generate large-scale fields with non-zero k, so it is important to relax this restriction. It is then also important to calculate $\eta_{ij}$ . The test field method allows the simultaneous calculation of all the components of the $\alpha_{ij}$ and $\eta_{ij}$ tensors for arbitrary values of k.

The test field method has recently been criticised by Cattaneo & Hughes (2008) on the grounds that the test fields are arbitrary predetermined mean fields. They argue that the resulting turbulent transport coefficients will only be approximations to the true values unless the test fields are close to the actual mean fields - a criticism equally applicable to other methods using arbitrary uniform fields. This concern has been addressed by Tilgner & Brandenburg (2008), who argue that Eq. (7) can instead be applied to any mean field. This statement has been numerically verified in three cases that we describe below.

Firstly, the test field method correctly reproduces a vanishing growth rate in saturated nonlinear cases (Brandenburg et al. 2008b). Secondly, in the time-dependent case, the test field method correctly reproduces also a non-vanishing growth rate, but in that case it is no longer permissible to express $\overline{\mbox{\boldmath${\cal E}$ }}{}$ in terms of a multiplication of turbulent transport coefficients with the mean field and its spatial derivatives. One must therefore write Eq. (8) as a convolution in time (Hubbard & Brandenburg 2008). Finally, the success of the test field method becomes particularly clear when it is applied to a passive vector field that obeys a separate induction equation with a velocity field from a saturated dynamo (Tilgner & Brandenburg 2008). This question was originally posed by Cattaneo & Tobias (2008). In particular for the Roberts flow with a mean field of Beltrami type, e.g. one that is proportional to $(\cos k_1z,\sin k_1z, 0)$ , the $\alpha_{ij}$ tensor is anisotropic and has an additional component proportional to $\meanB_i\meanB_j$ that tends to quench the components of the isotropic part of $\alpha_{ij}$ . The fastest growing passive vector field is then proportional to $(\sin k_1z,-\cos k_1z, 0)$ . This result has been confirmed both numerically and using weakly nonlinear theory (Tilgner & Brandenburg 2008).

In the following we ignore the complications involving time-dependent mean fields and restrict ourselves to transport coefficients that apply strictly speaking only to the time-independent or marginally excited case. For our numerical simulations we use the PENCIL CODE, where the test field algorithm has already been implemented. All our numerical simulations are performed in a periodic cubic box. The forcing scale is chosen to have the wavenumber $k_{\rm f}/k_1=5$ . This gives enough scale separation for a large-scale field to develop (Haugen et al. 2004), and is still not too big to reduce the resulting Reynolds numbers too much. We use units such that $c_{\rm s}=k_1=\rho_0=\mu_0=1$ and arrange the forcing amplitude so that the Mach number is around 0.1. For the magnetic Prandtl number we choose $\mbox{Pr}_{\rm M}=\nu/\eta=1$ , where $\nu$ lies in the range $3\times10^{-3}$ to $2\times10^{-4}$ (in units of $c_{\rm s}/k_1$ ). We choose the shear S such that the parameter

$\begin{displaymath}\mbox{Sh}\equiv S/(u_{\rm rms} k_{\rm f}) \end{displaymath}$

(14)

takes values in the range -0.02 to -0.9, where $u_{\rm rms}$ is the root-mean-square velocity. All the runs are started with uniform density, $\rho=\rho_0$ , $\vec{ U}=0$ , and $\vec{ a}^{pq}=0$ . Depending upon the parameters of a particular run we use up to 256³ grid points. For each test field ${\overline{\vec{ B}}}^{pq}$ we calculate the time averages of the $\alpha_{ij}$ and $\eta_{ij}$ tensors over time intervals over which $u_{\rm rms}$ is statistically stationary. In runs with higher magnetic Reynolds numbers,

$\begin{displaymath}\mbox{Re}_{\rm M}\equiv u_{\rm rms}/(\eta k_{\rm f}), \end{displaymath}$

(15)

we obtain an exponential growth of small-scale magnetic field (see below). We interpret this as being associated with the $\overline{\mbox{\boldmath${\cal E}$ }}{}_0^{pq}(z,t)$ term. This often gives rise to large fluctuations in all components of the $\alpha_{ij}$ and $\eta_{ij}$ tensors at late times. In such cases we confine our calculations of the time averages to time intervals over which $\vec{b}^{pq}_{\rm rms}$ does not exceed ${\overline{\vec{ B}}}^{pq}$ by more than a factor of about 20. Up until this point, the components of the $\alpha_{ij}$ and $\eta_{ij}$ tensors show only their intrinsic fluctuations, but at later times these will be swamped by additional contributions from the small-scale dynamo that grows exponentially in time. If necessary, we repeat our calculations over several independent realizations by resetting $\vec{ a}^{pq}=0$ at regular time intervals.

3 Results

Our principal aim in this paper is to study the effects of varying shear ( $\mbox{Sh}$ ) and magnetic Reynolds number ( $\mbox{Re}_{\rm M}$ ) on the components of the $\alpha_{ij}$ and $\eta_{ij}$ tensors. In the subsection below we summarize our results concerning the different components of these tensors.

3.1 Diagonal components of the transport tensors

The isotropic parts of the $\alpha_{ij}$ and $\eta_{ij}$ tensors are respectively characterised as

$\begin{displaymath}\alpha = {\textstyle\frac{1}{2}}\langle \alpha_{11}+\alpha_{2... ...= {\textstyle\frac{1}{2}} \langle \eta_{11}+\eta_{22} \rangle, \end{displaymath}$

(16)

where $\langle \cdot \rangle$ denotes an average taken over z and t. We normalize these quantities by $\alpha_{0}=-{\textstyle\frac{1}{3}}u_{\rm rms}$ and $\eta_{{\rm t}0}={\textstyle\frac{1}{3}}u_{\rm rms}/k_{\rm f}$ which are their respective expressions obtained using the First Order Smoothing Approximation (FOSA) for $\mbox{Re}_{\rm M}\ll1$ , which has previously been confirmed with the test field method in simulations of helical turbulence without shear (Sur et al. 2008).

Since $u_{\rm rms}$ enters the normalization of $\alpha_{ij}$ and $\eta_{ij}$ , it is useful to first look at how it changes as a function of $\mbox{Re}_{\rm M}$ (Fig. 1a) and $\mbox{Sh}$ (Fig. 1b). As can be seen, $u_{\rm rms}$ increases as a function of $\mbox{Re}_{\rm M}$ for small $\mbox{Re}_{\rm M}\sim 1$ and then reaches a plateau for high $\mbox{Re}_{\rm M}$ . On the other hand, $u_{\rm rms}$ is almost a constant as a function of shear except for high $\mbox{Sh}$ , where we observe excitation of the vorticity dynamo discussed further in Sect. 3.4.

In Fig. 2 we show $\alpha$ as a function of shear for three different values of the magnetic Reynolds number $\mbox{Re}_{\rm M}\approx 1, 20$ and 72. Figure 2b shows the corresponding results for $\eta _{\rm t}$ . We note that for small shear the turbulent transport coefficients are close to their values for zero shear. As $\vert\mbox{Sh}\vert$ increases, $\alpha$ decreases and $\eta _{\rm t}$ increases up to four times $\eta_{\rm t0}$ . As can be seen there is a clear reduction (quenching) of $\alpha$ with increasing shear in all these cases. In order to examine the possible convergence of the results with $\mbox{Re}_{\rm M}$ we plot $\alpha$ and $\eta _{\rm t}$ as functions of $\mbox{Re}_{\rm M}$ for three different values of the shear parameter $\mbox{Sh}$ , see Fig. 3. For small $\mbox{Re}_{\rm M}\sim 1$ we observe an increase in both $\alpha$ and $\eta _{\rm t}$ with $\mbox{Re}_{\rm M}$ . A similar initial increase of $\alpha$ and $\eta _{\rm t}$ was also seen in earlier simulations of helical turbulence without shear (Sur et al. 2008) and in non-helical shear flow turbulence (Brandenburg et al. 2008a). For higher values of $\mbox{Re}_{\rm M}$ and $\mbox{Sh}$ , both $\alpha$ and $\eta _{\rm t}$ show large variations. In these kinematic simulations, however, we expect them to tend to constant values asymptotically at high $\mbox{Re}_{\rm M}$ (Sur et al. 2008).

$\begin{figure} \par\includegraphics[width=12cm]{0359fig2} \end{figure}$	Figure 2: Turbulent transport coefficients a) $\alpha /\alpha _0$ and b) $\eta _{\rm t}/\eta _{{\rm t}0}$ as functions of shear parameter $\mbox{Sh}$ for different values of $\mbox{Re}_{\rm M}$ : $\mbox{Re}_{\rm M}\sim 1$ ( $\ast$ ), $\mbox{Re}_{\rm M}\sim 20$ ( $\triangle$ ), and $\mbox{Re}_{\rm M}\sim 72$ ( $\square$ ).
Open with DEXTER

$\begin{figure} \par\includegraphics[width=12cm]{0359fig3} \end{figure}$	Figure 3: a) $\alpha /\alpha _0$ and b) $\eta _{\rm t}/\eta _{{\rm t}0}$ as functions of $\mbox{Re}_{\rm M}$ , for different values of $\mbox{Sh}$ : $\mbox{Sh}\sim 0.07~ (\ast), 0.2~ (\triangle)$ , and $0.3~ (\Box)$ respectively. Horizontal dashed lines at $\alpha /\alpha _0 = 1$ and $\eta _{\rm t}/\eta _{{\rm t}0}=1$ are added to facilitate comparison.
Open with DEXTER

In an earlier study of shear flow turbulence with non-helical forcing by Brandenburg et al. (2008a), the diagonal components of $\eta_{ij}$ were found to be the same. Obviously, in the absence of helicity all the components of $\alpha_{ij}$ are zero. It turns out that in the presence of helicity the two diagonal components of $\eta_{ij}$ are still the same, but those of $\alpha_{ij}$ are now non-zero and not equal to each other. This is best shown by considering the quantities

$\begin{displaymath}\epsilon_{\alpha} = {\textstyle\frac{1}{2}}\langle \alpha_{11... ... {\textstyle\frac{1}{2}} \langle \eta_{11}-\eta_{22} \rangle . \end{displaymath}$

(17)

The results are shown in Fig. 4. Note that especially for large values of $\mbox{Re}_{\rm M}$ the values of $\epsilon _{\alpha }$ are predominantly negative. Since both $\alpha_{11}$ and $\alpha_{22}$ are negative, this means that $\vert\alpha_{11}\vert$ is larger than $\vert\alpha_{22}\vert$ , although the relative difference is only of the order of at most 5 per cent.

$\begin{figure} \par\includegraphics[width=12cm]{0359fig4} \end{figure}$	Figure 4: Normalized a) $\epsilon _{\alpha }$ and b) $\epsilon _{\eta })$ , as defined in Eq. (17), as functions of $\mbox{Re}_{\rm M}$ for constant value of the shear parameter $\mbox{Sh}$ : $\mbox{Sh}\sim 0.07~ (\ast),~ 0.2~ (\triangle)$ , and $0.3~ (\Box)$ respectively. Horizontal dashed lines are added to facilitate comparison.
Open with DEXTER

$\begin{figure} \par\includegraphics[width=12cm]{0359fig5} \end{figure}$	Figure 5: Normalized off-diagonal components of the $\alpha$ and $\eta$ tensors as functions of the shear parameter $\mbox{Sh}$ for different values of $\mbox{Re}_{\rm M}$ : $\mbox{Re}_{\rm M}\sim 1$ ( $\ast$ ), $\mbox{Re}_{\rm M}\sim 20$ ( $\triangle$ ), and $\mbox{Re}_{\rm M}\sim 72$ ( $\square$ ). Horizontal dashed lines are added to facilitate comparison.
Open with DEXTER

$\begin{figure} \par\includegraphics[width=12cm]{0359fig6} \end{figure}$	Figure 6: Normalized off-diagonal components of the $\alpha$ and $\eta$ tensors as functions of $\mbox{Re}_{\rm M}$ for constant values of $\mbox{Sh}$ : $\mbox{Sh}\sim 0.07~(\ast),~ 0.2~(\triangle),~ \mbox{and} \ 0.3 ~(\Box)$ . Horizontal dashed lines are added to facilitate comparison.
Open with DEXTER

3.2 Off-diagonal components

We now consider the off-diagonal components of $\alpha_{ij}$ and $\eta_{ij}$ . The results are depicted in Figs. 5 and 6. Of particular interest among these is the component $\eta_{21}$ , whose dependence on $\mbox{Sh}$ and $\mbox{Re}_{\rm M}$ is shown in Figs. 5d and 6d, respectively. This component can indicate the possible presence of a shear-current dynamo that may operate when $\eta_{21}S/(\eta_{\rm T} k_1)^2>1$ (Rogachevskii & Kleeorin 2003,2004). Here $\eta_{\rm T}=\eta_{\rm t}+\eta$ is the total (sum of turbulent and microscopic) magnetic diffusivity. In our case we have S < 0, which implies that the necessary condition for the shear-current dynamo to operate is $\eta_{21} < 0$ . As can be seen from Fig. 6d, for the range of $\mbox{Re}_{\rm M}$ values considered here $\eta_{21}$ is positive for small shear but becomes negative for strong shear and certain values of $\mbox{Re}_{\rm M}$ . Earlier work of Brandenburg et al. (2008a) without helicity did indicate a similar sign change, although only for larger $\mbox{Re}_{\rm M}$ . However, the error bars were so large that this result could not be regarded as significant. For the run with the strongest shear ( $-\mbox{Sh}\approx0.3$ ) and for $\mbox{Re}_{\rm M}=40$ we now find $\eta_{21}$ to be more clearly negative, but for smaller and larger values of $\mbox{Re}_{\rm M}$ the results are again, at least within error bars, compatible with zero. Also, of course, the present results apply to the case with helicity and are therefore not really comparable with those of Brandenburg et al. (2008a), where the helicity is zero. The antisymmetric contributions to the $\alpha_{ij}$ and $\eta_{ij}$ tensors are characterised by the vectors

$\begin{displaymath}\gamma_k= -{\textstyle{1\over2}}\epsilon_{ijk}\alpha_{ij},\quad \delta_k= -{\textstyle{1\over2}}\epsilon_{ijk}\eta_{ij}. \end{displaymath}$

(18)

Since our averages depend only on z, the z components of these tensors are irrelevant and therefore only the z components of the $\vec{\gamma}$ and $\vec{\delta}$ vectors are of interest. We denote those simply by $\gamma$ and $\delta$ , with

$\begin{displaymath}\gamma = {\textstyle\frac{1}{2}}\langle \alpha_{21}-\alpha_{1... ...= {\textstyle\frac{1}{2}} \langle \eta_{21}-\eta_{12} \rangle. \end{displaymath}$

(19)

A time series of both quantities is shown in Fig. 7 for positive and negative signs of the kinetic helicity. The results show that, in our case with S<0, $\gamma$ is positive (negative) for positive (negative) kinetic helicity, whilst $\delta$ has always the same sign.

$\begin{figure} \par\includegraphics[width=11cm]{0359fig7.ps} \end{figure}$

Figure 7:

Normalized time series of turbulent transport coefficients, a) $\alpha /\vert\alpha _0\vert$ , b) $\eta _t/\eta _{t0}$ , c) $\gamma /\vert\alpha _0\vert$ and d) $\delta /\eta _{t0}$ , from two runs with exactly the same parameters ( $\mbox{Re}_{\rm M}\sim 1.4$ , $\mbox{Sh}=-0.5$ ), but different signs of helicity for the external force. Continuous and broken lines denote results from runs with positive and negative helicities, respectively.

Open with DEXTER

Some idea about the functional forms of $\gamma$ and $\delta$ can be obtained from symmetry considerations. The vectors $\vec{\gamma}$ and $\vec{\delta}$ enter the electromotive force thus

$\begin{displaymath}\overline{\mbox{\boldmath${\cal E}$ }}{}=...+\vec{\gamma}\times\overline{\vec{B}}-\mu_0\vec{\delta}\times\overline{\vec{J}}, \end{displaymath}$

(20)

so we see that $\vec{\gamma}$ must be a polar vector and $\vec{\delta}$ must be an axial vector. Using the shear flow $\overline{\vec{ U}}^{S}$ , the only axial vector that can be constructed is the mean vorticity, $\overline{\vec{W}}=\vec{\nabla}\times\overline{\vec{ U}}^{S}$ , so we expect $\vec{\delta}$ to have a component that is proportional to $\overline{\vec{W}}$ . Likewise, since $\vec{\gamma}$ is a polar vector which points in the direction of $\overline{\vec{W}}$ , the two must be related via a pseudoscalar. In the present case the only pseudoscalar available is the kinetic helicity, $\overline{\vec{\omega}\cdot\vec{u}}$ . Based on these symmetry arguments we write

$\begin{displaymath}\vec{\gamma}=C_\gamma\tau^2\overline{\vec{\omega}\cdot\vec{ u... ...delta}=C_\delta\tau^2\overline{\vec{ u}^2}~\overline{\vec{W}}, \end{displaymath}$

(21)

where we have introduced two non-dimensional quantities, $C_\gamma$ and $C_\delta$ which could be either positive or negative, and $\tau$ is a correlation time that we approximate here by

$\begin{displaymath}\tau=(u_{\rm rms}k_{\rm f})^{-1}. \end{displaymath}$

(22)

In the present case, $\overline{\vec{\omega}\cdot\vec{u}}$ itself is a negative multiple of $\alpha$ . We therefore expect $\vec{\gamma}$ to have a component proportional to $\alpha\overline{\vec{W}}$ , multiplied by the correlation time $\tau$ . By similar arguments we expect $\vec{\delta}$ to have a component proportional to $\eta_{\rm t}\overline{\vec{W}}$ , multiplied by $\tau$ . Based on these arguments we can also write

$\begin{displaymath}\vec{\gamma}=\tilde{C}_\gamma~\alpha\overline{\vec{W}}\tau,\q... ...c{\delta}=\tilde{C}_\delta~\eta_{\rm t}\overline{\vec{W}}\tau, \end{displaymath}$

(23)

with new non-dimensional quantities, $\tilde{C}_\gamma\approx C_\gamma/3$ and $\tilde{C}_\delta\approx C_\delta/3$ , that we expect to be of order unity. Our simulations confirm this reasoning and show that both coefficients are of order unity with positive $\tilde{C}_\gamma$ and negative $\tilde{C}_\delta$ (e.g., $\tilde{C}_\gamma= 0.5$ and $\tilde{C}_\delta= 0.25$ for the run shown in Fig. 8) for runs with small $\mbox{Re}_{\rm M}$ and $\mbox{Sh}$ ; see Figs. 8 and 9, respectively. The sign in Eq. (23) has been chosen a posteriori from our simulations. For high shear and Reynolds number this simple reasoning is no longer accurate and, at least in one case (Fig. 9a), $C_\gamma$ even becomes negative. For lower values of $\mbox{Sh}$ we find an almost constant negative $\gamma$ of the order of $0.05\alpha_0$ (or less). This result is fairly independent of $\mbox{Re}_{\rm M}$ . The moduli of $\gamma$ and $\delta$ increase with increasing shear parameter, and for larger values of $\mbox{Re}_{\rm M}$ both quantities may approach an asymptotic value. For small $\mbox{Re}_{\rm M}$ we observe a roughly linear increase of $\gamma$ with shear (see inset of Fig 8a), further verifying Eq. (21).

The coefficient $\gamma$ can be interpreted as turbulent pumping, i.e. advection of the magnetic field by means other than the mean velocity field. In strongly stratified convection, turbulent pumping has been seen to be directed from higher to lower turbulence intensity (Käpylä et al. 2006; Tobias et al. 1998; Ossendrijver et al. 2002; Tobias et al. 2001), which is usually in the downward direction. Thus, turbulent pumping is likely to play an important role in convection zones of stars where it can overcome the buoyancy of the magnetic field. In the present case where the turbulence is homogeneous, however, stratification does not play a role and the pumping is just due to the combined action of shear and helical turbulence.

$\begin{figure} \par\includegraphics[width=12cm]{0359fig8.ps} \end{figure}$	Figure 8: Normalized a) $\gamma$ and b) $\delta$ as functions of $\mbox{Sh}$ for $\mbox{Re}_{\rm M}\approx 1$ ( $\ast$ ), $\mbox{Re}_{\rm M}\approx 20$ ( $\triangle$ ), and $\mbox{Re}_{\rm M}\approx 72$ ( $\square$ ). Horizontal dashed lines are added to facilitate comparison. The inset shows $\gamma /\vert\alpha _0\vert$ versus $-\mbox{Sh}$ for $\mbox{Re}_{\rm M}\sim 1$ .
Open with DEXTER

$\begin{figure} \par\includegraphics[width=12cm]{0359fig9.ps} \end{figure}$	Figure 9: Normalized a) $\gamma$ and b) $\delta$ as functions of $\mbox{Re}_{\rm M}$ for constant values of the shear parameter $\mbox{Sh}$ : $\mbox{Sh}\sim 0.07~(\ast),~ 0.2~ (\triangle), \mbox{and} \ 0.3~ (\Box)$ . Horizontal dashed lines are added to facilitate comparison.
Open with DEXTER

3.3 Scale-dependence

So far we have confined the calculations of the $\alpha_{ij}$ and $\eta_{ij}$ tensors to test fields whose characteristic length scale is the largest scale in the domain, i.e., test fields of the form $\sin kz$ or $\cos kz$ with k=k₁, where $k_1=2\pi/L$ and L is the box size of our simulations. We have also done similar calculations for other values of k, which correspond to test fields with different characteristic length scales. In Fig. 10 we show the dependence of $\alpha$ and $\eta_t$ on k for $\mbox{Re}_{\rm M}\approx3$ . The decrease can be modelled by a Lorentzian peaked at k=0 (Fig. 10). This result is not surprising, because Lorentzian fits have been obtained earlier also in the absence of shear (Brandenburg et al. 2008c).

$\begin{figure} \par\includegraphics[width=12cm]{0359fig10.ps} \end{figure}$	Figure 10: a) $\alpha (k)/\alpha _0$ and b) $\eta _t(k)/\eta _{\rm t0}$ as a function of the wavenumber k of the large-scale magnetic field for $\mbox{Re}_{\rm M}\approx3$ . The solid lines show the Lorentzian fits of width 1.4 $k_{\rm f}$ .
Open with DEXTER

$\begin{figure} \par\mbox{\includegraphics[width=6cm,clip]{0359f11a}\includegraphics[width=6cm,clip]{0359f11b} }\end{figure}$

Figure 11:

Panel a) shows the root-mean-square values of the velocity for one of our runs (with ${\rm Re}\approx 45$ and ${\rm Sh}\approx -0.3$ ) which exhibits a vorticity dynamo. The inset shows the time-series for $\eta _{\rm t}$ . Panel b) shows a pronounced pattern in U_y on the periphery of the box at the end of the same run. Similar structures are also found in the x component of the velocity.

Open with DEXTER

3.4 Effects from the vorticity dynamo

In the runs with relatively high values of shear and Reynolds number, the root-mean-square velocity, $U_{\rm rms}$ , which initially reaches a steady state, shows an exponential growth at late times, as can be seen in Fig. 11a. Similar behaviour is seen for turbulent transport coefficients, which also show large fluctuations for later times. Large fluctuations of the turbulent transport coefficients at later times are also seen by Sur et al. (2008), but those fluctuations are more irregular and have a different origin. They are interpreted as being due to the development of a small-scale dynamo at late times and are not due to the more systematic increase in the rms velocity, which in turn is associated with the shear. However, we emphasize that even before $\vec{ b}^{pq}$ becomes dominated by this type of dynamo action, the temporal and spatial fluctuations of $\alpha_{ij}$ and $\eta_{ij}$ are of the order of $\alpha_0$ and $\eta_{\rm t0}$ , respectively. This is true even for large values of $\mbox{Re}_{\rm M}$ ; see Brandenburg et al. (2008a), who have argued that these fluctuations also contribute to large-scale dynamo action via the incoherent $\alpha$ effect (Vishniac & Brandenburg 1997; Proctor 2007). For example, the time-series of $\eta _{\rm t}$ shown in the inset of Fig. 11 has a plateau even beyond the range over which $U_{\rm rms}$ is steady.

The vertical spikes in the inset of Fig. 11 come from resetting $\vec{ a}^{pq}$ to zero in regular time intervals; see the end of Sect. 2. Note, however, that the mean values and also the upper envelope trace the evolution of $U_{\rm rms}$ reasonably well, including the increased rise after t=300. The late time behaviour is accompanied by the formation of large-scale vortical structures, as seen in Fig. 11b, which is a signature of the vorticity dynamo proposed by Elperin et al. (2003), see also Yousef et al. (2008a,b). A detailed numerical study of the vorticity dynamo has been performed by Käpylä et al. (2008).

The presence of the vorticity dynamo and the resulting systematic variation in $U_{\rm rms}$ as well as the turbulent transport coefficients for late times limit the lengths of time over which average values of the components of the $\alpha_{ij}$ and $\eta_{ij}$ tensors can be calculated. In the above case this interval lies between the two arrows in Fig. 11a. This limits the range of $\mbox{Sh}$ and $\mbox{Re}_{\rm M}$ that we have probed and explains the larger error bars shown for example in Fig. 6, and why they cannot be reduced by simply running our simulations for longer times. This problem would be avoided in the presence of magnetic fields, because then the resulting Lorentz force would quench the vorticity dynamo (Käpylä & Brandenburg 2008). This is however beyond the scope of this paper.

4 Conclusions

We have studied the effects of varying shear and magnetic Reynolds number on the turbulent transport coefficients in the presence of helicity in the kinematic limit. We have shown that for fixed $\mbox{Re}_{\rm M}$ , $\alpha$ is reduced (quenched) with increasing shear. Despite the differences in the assumptions made, this quenching is qualitatively similar to the recent results obtained by Leprovost & Kim (2008). To the best of our knowledge this is the first numerical study to show quenching of $\alpha$ as a function of shear in helical turbulence. We find that $\eta _{\rm t}$ increases with increasing shear in the range of $\mbox{Re}_{\rm M}$ values that we have considered here. A similar behaviour for $\eta _{\rm t}$ was also seen in Brandenburg et al. (2008a) where the forcing was non-helical.

We also compute the off-diagonal components of the $\alpha_{ij}$ and $\eta_{ij}$ tensors. The antisymmetric part of $\alpha_{ij}$ corresponds to a turbulent pumping velocity $\vec{\gamma}$ in the direction perpendicular to the plane of the shear flow. It shows a roughly linear increase with shear for small magnetic Reynolds number. We propose simple expressions for $\gamma$ and $\delta$ in Eq. (21) which show reasonable agreement with our numerical results for small $\mbox{Sh}$ and $\mbox{Re}_{\rm M}$ . Our expression shows that for negative helicity, $\vec{\gamma}$ points in the direction of the vorticity of the mean flow. Regarding the component $\eta_{21}$ we find indications that, at least in one or two cases, this component changes sign and becomes negative. This could be of significance in connection with the shear-current effect.

We also find that all the turbulent transport coefficients depend on the wavenumber of the mean flow in a Lorentzian fashion, just as in the case of non-shearing turbulence (Brandenburg et al. 2008c). This means that the kinematic values of $\alpha$ and $\eta _{\rm t}$ for k=k₁ are close to the values obtained for $k\to0$ . This is not the case for certain non-turbulent flows such as the Galloway-Proctor flow (Rädler & Brandenburg 2009; Courvoisier 2008). In an earlier paper, Courvoisier et al. (2006) considered only the limiting case k=0 for this flow.

Several aspects of the present investigations could be of astrophysical relevance. Turbulence in celestial bodies is helical and exhibits an $\alpha$ effect. In addition, shear (S) can be an important ingredient in that the efficiency of large-scale dynamo action is determined by the product of $\alpha$ and S. However, as S increases, $\alpha$ itself becomes quenched when Sbecomes comparable with the inverse turnover time, i.e. $S\tau=O(1)$ . Furthermore, the turbulent diffusivity becomes enhanced, suppressing the dynamo even further. Finally, it is found that the combined action of helicity and shear gives rise to a pumping velocity of mean magnetic field perpendicular to the plane of the shear flow. The existence of such a pumping velocity has not been emphasized before. On the other hand, it is well known that $\alpha\Omega$ (or rather $\alpha S$ ) dynamos can have travelling wave solutions (Brandenburg et al. 2001). When the product of $\alpha$ and S is positive, these waves travel in the positive z direction, which agrees with the direction of pumping. In the near-surface shear layer of the Sun this pumping would therefore support the equatorward migration in that layer.

It is important to understand the quenching of turbulent transport coefficients in the presence of shear beyond the kinematic approximation. In that case one needs to include the induction equation (in addition to the test field equations) and incorporate the resulting Lorentz force in the momentum equation (Brandenburg et al. 2008b). This would also help in alleviating problems of strong late-time fluctuations arising from the vorticity dynamo, because the vorticity dynamo tends to be suppressed by magnetic fields of equipartition strengths (Käpylä & Brandenburg 2008). Similarly, given their potential importance in allowing the escape of magnetic helicity, the effects of open boundary conditions also needs to be considered. These questions are under study and will be reported elsewhere.

Acknowledgements

The authors acknowledge the hospitality of Nordita during the programme ``Turbulence and Dynamos''. A.B. and P.J.K. thank Astronomy Unit, Queen Mary University of London, for hospitality. Computational resources were granted by CSC (Espoo, Finland), UKMHD, and QMUL HPC facilities purchased under the SRIF initiative. This work was supported by the the Leverhulme Trust (DM, RT), the Academy of Finland grant No. 121431 (PJK), and the Swedish Research Council (AB).

References

Blackman, E., & Field, G. B. 2000, ApJ, 534, 984 [NASA ADS] [CrossRef]
Blackman, E., & Brandenburg, A. 2002, ApJ, 579, 359 [NASA ADS] [CrossRef]
Brandenburg, A. 2001, ApJ, 550, 824 [NASA ADS] [CrossRef] (In the text)
Brandenburg, A. 2005, ApJ, 625, 539 [NASA ADS] [CrossRef]
Brandenburg, A. 2008, ApJ, submitted [arXiv:0808.0961]
Brandenburg, A., & Sandin, C. 2004, A&A, 427, 13 [NASA ADS] [CrossRef] [EDP Sciences]
Brandenburg, A., & Subramanian, K. 2005, Phys. Rep., 417, 1 [NASA ADS] [CrossRef] (In the text)
Brandenburg, A., Bigazzi, A., & Subramanian, K. 2001, MNRAS, 325, 685 [NASA ADS] [CrossRef] (In the text)
Brandenburg, A., Rädler, K.-H., Rheinhardt, M., & Käpylä, P. J. 2008a, ApJ, 676, 740 [NASA ADS] [CrossRef] (In the text)
Brandenburg, A., Rädler, K.-H., Rheinhardt, M., & Subramanian, K. 2008b, ApJ, 687, L49 [NASA ADS] [CrossRef] (In the text)
Brandenburg, A., Rädler, K.-H., & Schrinner, M. 2008c, A&A, 482, 739 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
Cattaneo, F., & Hughes, D. W. 2008 [arXiv:0805.2138] (In the text)
Cattaneo, F., & Tobias, S. M. 2008, J. Fluid Mech. [arXiv:0809.1801] (In the text)
Courvoisier, A. 2008, GAFD, 102, 217 [CrossRef]
Courvoisier, A., Hughes, D. W., & Tobias, S. M. 2006, Phys. Rev. Lett., 96, 034503 [NASA ADS] [CrossRef] (In the text)
Elperin, T., Kleeorin, N., & Rogachevskii, I. 2003, Phys. Rev. E, 68, 016311 [NASA ADS] [CrossRef] (In the text)
Field, G., & Blackman, E. 2002, ApJ, 572, 685 [NASA ADS] [CrossRef]
Gruzinov, A. V., & Diamond, P. H. 1994, Phys. Rev. Lett., 72, 1651 [NASA ADS] [CrossRef] (In the text)
Haugen, N. E. L., Brandenburg, A., & Dobler, W. 2004, Phys. Rev. E 70, 016308 [NASA ADS]
Hubbard, A., & Brandenburg, A. 2008 [arXiv:0811.2561] (In the text)
Hughes, D. W., & Proctor, M. R. E. 2009, Phys. Rev. Lett. 102, 044501 (In the text)
Käpylä, P. J., & Brandenburg, A. 2008, ApJ, submitted [arXiv:0810.2298] (In the text)
Käpylä, P. J., Korpi, M. J., Ossendrijver, M., & Stix, M. 2006, A&A, 455, 401 [NASA ADS] [CrossRef] [EDP Sciences]
Käpylä, P. J., Korpi, M. J., & Brandenburg, A. 2008a, A&A, 491, 353 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
Käpylä, P. J., Korpi, M. J., & Brandenburg, A. 2008b, A&A, submitted [arXiv:0812.1792] (In the text)
Käpylä, P. J., Mitra, D., & Brandenburg, A. 2008, Phys. Rev. E, 79, 016302 (In the text)
Kitchatinov, L. L., Rüdiger, G., & Pipin, V. V. 1994, AN, 315, 157 [NASA ADS]
Kleeorin, N., Moss, D., Rogachevskii, I., & Sokoloff, D. 2000, A&A, 361, L5 [NASA ADS]
Krause, F., & Rädler, K.-H. 1980, Mean-field magnetohydrodynamics and dynamo theory (Oxford: Pergamon Press) (In the text)
Leprovost, N., & Kim, E.-J. 2008, Phys. Rev. Lett., 100, 144502 [NASA ADS] [CrossRef]
Ossendrijver, M., Stix, M., Brandenburg, A., & Rüdiger, G. 2002, A&A, 394, 735 [NASA ADS] [CrossRef] [EDP Sciences]
Piddington, J. H. 1981 ApJ, 247, 291
Pipin, V. V., Rüdiger G., & Kitchatinov, L. L. 1996, GAFD, 83, 119 [CrossRef]
Proctor, M. R. E. 2007, MNRAS, 382, L38 [NASA ADS]
Rädler, K.-H., & Stepanov, R. 2006, Phys. Rev. E, 73, 056311 [NASA ADS] [CrossRef]
Rädler, K.-H., & Brandenburg, A. 2009, MNRAS, 393, 113 [CrossRef]
Rogachevskii, I., & Kleeorin, N. 2003, Phys. Rev. E, 68, 036301 [NASA ADS] [CrossRef]
Rogachevskii, I., & Kleeorin, N. 2004, Phys. Rev. E, 70, 046310 [NASA ADS] [CrossRef]
Rüdiger, G., & Kitchatinov, L. 2006, Astron. Nachr., 327, 298 [NASA ADS] [CrossRef]
Schrinner, M., Rädler, K.-H., Schmitt, D., Rheinhardt, M., & Christensen, U. 2005, Astron. Nachr., 326, 245 [NASA ADS] [CrossRef]
Schrinner, M., Rädler, K.-H., Schmitt, D., Rheinhardt, M., & Christensen, U. 2007, GAFD, 101, 81 [CrossRef]
Subramanian, K., & Brandenburg, A. 2004, Phys. Rev. Lett., 93, 205001 [NASA ADS] [CrossRef]
Sur, S., Brandenburg, A., & Subramanian, K. 2008, MNRAS, 385, L15 [NASA ADS] (In the text)
Tilgner, A., & Brandenburg, A. 2008, MNRAS, 391, 1477 [NASA ADS] [CrossRef] (In the text)
Tobias, S. M., Brummell, N. H., Clune, T. L., & Toomre, J. 1998, ApJ, 502, L177 [NASA ADS] [CrossRef]
Tobias, S. M., Brummell, N. H., Clune, T. L., & Toomre, J. 2001, ApJ, 549, 1183 [NASA ADS] [CrossRef]
Vishniac, E. T., & Brandenburg, A. 1997, ApJ, 584, L99
Vishniac, E. T., & Cho, J. 2001, ApJ, 550, 752 [NASA ADS] [CrossRef]
Vainshtein, E. T., & Cattaneo, J. 1992, ApJ, 393, 165 [NASA ADS] [CrossRef]
Yousef, T. A., Heinemann, T., Schekochihin, A. A., et al. 2008a, Phys. Rev. Lett., 100, 184501 [NASA ADS] [CrossRef]
Yousef, T. A., Heinemann, T., Rincon, F., et al. 2008b, AN, 329, 737 [NASA ADS]

Footnotes

...ODE: http://www.nordita.org/software/pencil-code

All Figures

$\begin{figure} \mbox{\includegraphics[width=6cm,clip]{0359fig1a}\includegraphics[width=6cm,clip]{0359fig1b} }\end{figure}$	Figure 1: Root-mean-square velocity $u_{\rm rms}$ normalized by the speed of sound $c_{\rm s}$ as: a) a function of $\mbox{Re}_{\rm M}$ for different values of $\mbox{Sh}$ : $\mbox{Sh}\sim 0.07 (\ast), 0.2 (\triangle), \mbox{and}\ 0.3 (\Box)$ respectively and b) a function of the shear parameter $\mbox{Sh}$ for different value of $\mbox{Re}_{\rm M}$ : $\mbox{Re}_{\rm M}\sim 1$ ( $\ast$ ), $\mbox{Re}_{\rm M}\sim 20$ ( $\triangle$ ), and $\mbox{Re}_{\rm M}\sim 72$ ( $\square$ ).
Open with DEXTER
In the text

$\begin{figure} \par\includegraphics[width=12cm]{0359fig2} \end{figure}$	Figure 2: Turbulent transport coefficients a) $\alpha /\alpha _0$ and b) $\eta _{\rm t}/\eta _{{\rm t}0}$ as functions of shear parameter $\mbox{Sh}$ for different values of $\mbox{Re}_{\rm M}$ : $\mbox{Re}_{\rm M}\sim 1$ ( $\ast$ ), $\mbox{Re}_{\rm M}\sim 20$ ( $\triangle$ ), and $\mbox{Re}_{\rm M}\sim 72$ ( $\square$ ).
Open with DEXTER
In the text

$\begin{figure} \par\includegraphics[width=12cm]{0359fig3} \end{figure}$	Figure 3: a) $\alpha /\alpha _0$ and b) $\eta _{\rm t}/\eta _{{\rm t}0}$ as functions of $\mbox{Re}_{\rm M}$ , for different values of $\mbox{Sh}$ : $\mbox{Sh}\sim 0.07~ (\ast), 0.2~ (\triangle)$ , and $0.3~ (\Box)$ respectively. Horizontal dashed lines at $\alpha /\alpha _0 = 1$ and $\eta _{\rm t}/\eta _{{\rm t}0}=1$ are added to facilitate comparison.
Open with DEXTER
In the text

$\begin{figure} \par\includegraphics[width=12cm]{0359fig4} \end{figure}$	Figure 4: Normalized a) $\epsilon _{\alpha }$ and b) $\epsilon _{\eta })$ , as defined in Eq. (17), as functions of $\mbox{Re}_{\rm M}$ for constant value of the shear parameter $\mbox{Sh}$ : $\mbox{Sh}\sim 0.07~ (\ast),~ 0.2~ (\triangle)$ , and $0.3~ (\Box)$ respectively. Horizontal dashed lines are added to facilitate comparison.
Open with DEXTER
In the text

$\begin{figure} \par\includegraphics[width=12cm]{0359fig5} \end{figure}$	Figure 5: Normalized off-diagonal components of the $\alpha$ and $\eta$ tensors as functions of the shear parameter $\mbox{Sh}$ for different values of $\mbox{Re}_{\rm M}$ : $\mbox{Re}_{\rm M}\sim 1$ ( $\ast$ ), $\mbox{Re}_{\rm M}\sim 20$ ( $\triangle$ ), and $\mbox{Re}_{\rm M}\sim 72$ ( $\square$ ). Horizontal dashed lines are added to facilitate comparison.
Open with DEXTER
In the text

$\begin{figure} \par\includegraphics[width=12cm]{0359fig6} \end{figure}$	Figure 6: Normalized off-diagonal components of the $\alpha$ and $\eta$ tensors as functions of $\mbox{Re}_{\rm M}$ for constant values of $\mbox{Sh}$ : $\mbox{Sh}\sim 0.07~(\ast),~ 0.2~(\triangle),~ \mbox{and} \ 0.3 ~(\Box)$ . Horizontal dashed lines are added to facilitate comparison.
Open with DEXTER
In the text

$\begin{figure} \par\includegraphics[width=11cm]{0359fig7.ps} \end{figure}$	Figure 7: Normalized time series of turbulent transport coefficients, a) $\alpha /\vert\alpha _0\vert$ , b) $\eta _t/\eta _{t0}$ , c) $\gamma /\vert\alpha _0\vert$ and d) $\delta /\eta _{t0}$ , from two runs with exactly the same parameters ( $\mbox{Re}_{\rm M}\sim 1.4$ , $\mbox{Sh}=-0.5$ ), but different signs of helicity for the external force. Continuous and broken lines denote results from runs with positive and negative helicities, respectively.
Open with DEXTER
In the text

$\begin{figure} \par\includegraphics[width=12cm]{0359fig8.ps} \end{figure}$	Figure 8: Normalized a) $\gamma$ and b) $\delta$ as functions of $\mbox{Sh}$ for $\mbox{Re}_{\rm M}\approx 1$ ( $\ast$ ), $\mbox{Re}_{\rm M}\approx 20$ ( $\triangle$ ), and $\mbox{Re}_{\rm M}\approx 72$ ( $\square$ ). Horizontal dashed lines are added to facilitate comparison. The inset shows $\gamma /\vert\alpha _0\vert$ versus $-\mbox{Sh}$ for $\mbox{Re}_{\rm M}\sim 1$ .
Open with DEXTER
In the text

$\begin{figure} \par\includegraphics[width=12cm]{0359fig9.ps} \end{figure}$	Figure 9: Normalized a) $\gamma$ and b) $\delta$ as functions of $\mbox{Re}_{\rm M}$ for constant values of the shear parameter $\mbox{Sh}$ : $\mbox{Sh}\sim 0.07~(\ast),~ 0.2~ (\triangle), \mbox{and} \ 0.3~ (\Box)$ . Horizontal dashed lines are added to facilitate comparison.
Open with DEXTER
In the text

$\begin{figure} \par\includegraphics[width=12cm]{0359fig10.ps} \end{figure}$	Figure 10: a) $\alpha (k)/\alpha _0$ and b) $\eta _t(k)/\eta _{\rm t0}$ as a function of the wavenumber k of the large-scale magnetic field for $\mbox{Re}_{\rm M}\approx3$ . The solid lines show the Lorentzian fits of width 1.4 $k_{\rm f}$ .
Open with DEXTER
In the text

$\begin{figure} \par\mbox{\includegraphics[width=6cm,clip]{0359f11a}\includegraphics[width=6cm,clip]{0359f11b} }\end{figure}$	Figure 11: Panel a) shows the root-mean-square values of the velocity for one of our runs (with ${\rm Re}\approx 45$ and ${\rm Sh}\approx -0.3$ ) which exhibits a vorticity dynamo. The inset shows the time-series for $\eta _{\rm t}$ . Panel b) shows a pronounced pattern in U_y on the periphery of the box at the end of the same run. Similar structures are also found in the x component of the velocity.
Open with DEXTER
In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.