A quasi-linear theory is presented for how randomly forced, barotropic velocity fluctuations cause an exponentially-growing, large-scale (mean) magnetic dynamo in the presence of a uniform shear flow, $\vec{U} = S x \vec{e}_y$. It is a "kinematic" theory for the growth of the mean magnetic energy from a small initial seed, neglecting the saturation effects of the Lorentz force. The quasi-linear approximation is most broadly justifiable by its correspondence with computational solutions of nonlinear magneto-hydrodynamics, and it is rigorously derived in the limit of large resistivity, $\eta \rightarrow \infty$. Dynamo action occurs even without mean helicity in the forcing or flow, but random helicity variance is then essential. In a sufficiently large domain and with small wavenumber $k_z$ in the direction perpendicular to the mean shearing plane, a positive exponential growth rate $\gamma$ can occur for arbitrary values of $\eta$, the viscosity $\nu$, and the random-forcing correlation time $t_f$ and phase angle $\theta_f$ in the shearing plane. The value of $\gamma$ is independent of the domain size. The shear dynamo is "fast", with finite $\gamma > 0$ in the limit of $\eta \rightarrow 0$. Averaged over the random forcing ensemble, the mean magnetic field grows more slowly, if at all, compared to the r.m.s. field (or magnetic energy). In the limit of small Reynolds numbers ($\eta, \ \nu \rightarrow \infty$), the dynamo behavior is related to the well-known alpha--omega {\it ansatz} when the forcing is steady ($t_f \rightarrow \infty$) and to the "incoherent" alpha--omega {\it ansatz} when the forcing is purely fluctuating.

A rigorous theory for the generation of a large-scale magnetic field by random non-helically forced motions of a conducting fluid combined with a linear shear is presented in the analytically tractable limit of low Rm and weak shear. The dynamo is kinematic and due to fluctuations in the net (volume-averaged) electromotive force. This is a minimal proof-of-concept quasilinear calculation aiming to put the shear dynamo, a new effect recently found in numerical experiments, on a firm theoretical footing. Numerically observed scalings of the wavenumber and growth rate of the fastest growing mode, previously not understood, are derived analytically. The simplicity of the model suggests that shear dynamo action may be a generic property of sheared magnetohydrodynamic turbulence.

Numerical simulations of forced turbulence in elongated shearing boxes are carried out to demonstrate that a nonhelical turbulence in conjunction with a linear shear can give rise to a mean-field dynamo. Exponential growth of magnetic field at scales larger than the outer (forcing) scale of the turbulence is found. Over a range of values of the shearing rate S spanning approximately two orders of magnitude, the growth rate of the magnetic field is proportional to the imposed shear, gamma ~ S, while the characteristic spatial scale of the field is l_b ~ S^(-1/2). The effect is quite general: earlier results for the nonrotating case by Yousef et al. 2008 (PRL 100, 184501) are extended to shearing boxes with Keplerian rotation; it is also shown that the shear dynamo mechanism operates both below and above the threshold for the fluctuation dynamo. The apparently generic nature of the shear dynamo effect makes it an attractive object of study for the purpose of understanding the generation of magnetic fields in astrophysical systems.

The feasibility of a mean-field dynamo in nonhelical turbulence with superimposed linear shear is studied numerically in elongated shearing boxes. Exponential growth of magnetic field at scales much larger than the outer scale of the turbulence is found. The charateristic scale of the field is l_B ~ S^{-1/2} and growth rate is gamma ~ S, where S is the shearing rate. This newly discovered shear dynamo effect potentially represents a very generic mechanism for generating large-scale magnetic fields in a broad class of astrophysical systems with spatially coherent mean flows.

Direct numerical simulations of incompressible nonhelical randomly forced MHD turbulence are used to demonstrate for the first time that the fluctuation dynamo exists in the limit of large magnetic Reynolds number Rm>>1 and small magnetic Prandtl number Pm<<1. The dependence of the critical Rm_c for dynamo on the hydrodynamic Reynolds number Re is obtained for 1<Re<6700. In the limit Pm<<1, Rm_c is about three times larger than for the previously well established dynamo at large and moderate Prandtl numbers: Rm_c<200 for Re>6000 compared to Rm_c~60 for Pm>=1. Is is not as yet possible to determine numerically whether the growth rate of the magnetic energy is ~Rm^{1/2} in the limit Rm->infinity, as it should if the dynamo is driven by the inertial-range motions at the resistive scale.

This paper is a detailed report on a programme of simulations used to settle a long-standing issue in the dynamo theory and demonstrate that the fluctuation dynamo exists in the limit of large magnetic Reynolds number Rm>>1 and small magnetic Prandtl number Pm<<1. The dependence of the critical Rm_c vs. the hydrodynamic Reynolds number Re is obtained for 1<Re<6700. In the limit Pm<<1, Rm_c is ~3 times larger than for Pm>1. The stability curve Rm_c(Re) (and, it is argued, the nature of the dynamo) is substantially different from the case of the simulations and liquid-metal experiments with a mean flow. It is not as yet possible to determine numerically whether the growth rate is ~Rm^{1/2} in the limit Re>>Rm>>1, as should be the case if the dynamo is driven by the inertial-range motions. The magnetic-energy spectrum in the low-Pm regime is qualitatively different from the Pm>1 case and appears to develop a negative spectral slope, although current resolutions are insufficient to determine its asymptotic form. At 1<Rm<Rm_c, the magnetic fluctuations induced via the tangling by turbulence of a weak mean field are investigated and the possibility of a k^{-1} spectrum above the resistive scale is examined. At low Rm<1, the induced fluctuations are well described by the quasistatic approximation; the k^{-11/3} spectrum is confirmed for the first time in direct numerical simulations.

We study numerically the dependence of the critical magnetic Reynolds number Rmc for the turbulent small-scale dynamo on the hydrodynamic Reynolds number Re. The turbulence is statistically homogeneous, isotropic, and mirror--symmetric. We are interested in the regime of low magnetic Prandtl number Pm=Rm/Re<1, which is relevant for stellar convective zones, protostellar disks, and laboratory liquid-metal experiments. The two asymptotic possibilities are Rmc->const as Re->infinity (a small-scale dynamo exists at low Pm) or Rmc/Re=Pmc->const as Re->infinity (no small-scale dynamo exists at low Pm). Results obtained in two independent sets of simulations of MHD turbulence using grid and spectral codes are brought together and found to be in quantitative agreement. We find that at currently accessible resolutions, Rmc grows with Re with no sign of approaching a constant limit. We reach the maximum values of Rmc~500 for Re~3000. By comparing simulations with Laplacian viscosity, fourth-, sixth-, and eighth-order hyperviscosity and Smagorinsky large-eddy viscosity, we find that Rmc is not sensitive to the particular form of the viscous cutoff. This work represents a significant extension of the studies previously published by Schekochihin et al. 2004, PRL 92, 054502 and Haugen et al. 2004, PRE, 70, 016308 and the first detailed scan of the numerically accessible part of the stability curve Rmc(Re).

We report an extensive numerical study of the small-scale turbulent dynamo at large magnetic Prandtl numbers Pm. A Pm scan is given for the model case of low-Reynolds-number turbulence. We concentrate on three topics: magnetic-energy spectra and saturation levels, the structure of the field lines, and the field-strength distribution. The main results are (1) the folded structure (direction reversals at the resistive scale, field lines curved at the scale of the flow) persists from the kinematic to the nonlinear regime; (2) the field distribution is self-similar and appears to be lognormal during the kinematic regime and exponential in the saturated state; and (3) the bulk of the magnetic energy is at the resistive scale in the kinematic regime and remains there after saturation, although the spectrum becomes much shallower. We propose an analytical model of saturation based on the idea of partial two-dimensionalization of the velocity gradients with respect to the local direction of the magnetic folds. The model-predicted spectra are in excellent agreement with numerical results. Comparisons with large-Re, moderate-Pm runs are carried out to confirm the relevance of these results. New features at large Re are elongation of the folds in the nonlinear regime from the viscous scale to the box scale and the presence of an intermediate nonlinear stage of slower-than-exponential magnetic-energy growth accompanied by an increase of the resistive scale and partial suppression of the kinetic-energy spectrum in the inertial range. Numerical results for the saturated state do not support scale-by-scale equipartition between magnetic and kinetic energies, with a definite excess of magnetic energy at small scales. A physical picture of the saturated state is proposed.

We consider the problem of incompressible, forced, nonhelical, homogeneous, isotropic MHD turbulence with no mean magnetic field. This problem is essentially different from the case with externally imposed uniform mean field. There is no scale-by-scale equipartition between magnetic and kinetic energies as would be the case for the Alfven-wave turbulence. The isotropic MHD turbulence is the end state of the turbulent dynamo which generates folded fields with small-scale direction reversals. We propose that the statistics seen in numerical simulations of isotropic MHD turbulence could be explained as a superposition of these folded fields and Alfven-like waves that propagate along the folds.

We consider the problem of incompressible, forced, nonhelical, homogeneous and isotropic MHD turbulence with no mean magnetic field and large magnetic Prandtl number. This type of MHD turbulence is the end state of the turbulent dynamo, which generates folded fields with small-scale direction reversals. We propose a model in which saturation is achieved as a result of the velocity statistics becoming anisotropic with respect to the local direction of the magnetic folds. The model combines the effects of weakened stretching and quasi-two-dimensional mixing and produces magnetic-energy spectra in remarkable agreement with numerical results at least in the case of a one-scale flow. We conjecture that the statistics seen in numerical simulations could be explained as a superposition of these folded fields and Alfv\'en-like waves that propagate along the folds.

The amplification of magnetic fields in a highly conducting fluid is studied numerically. During growth, the magnetic field is spatially intermittent: it does not uniformly fill the volume, but is concentrated in long thin folded structures. Contrary to a commonly held view, intermittency of the folded field does not increase indefinitely throughout the growth stage if diffusion is present. Instead, as we show, the probability-density function (PDF) of the field strength becomes self-similar. The normalized moments increase with magnetic Prandtl number in a powerlike fashion. We argue that the self-similarity is to be expected with a finite flow scale and system size. In the nonlinear saturated state, intermittency is reduced and the PDF is exponential. Parallels are noted with self-similar behavior recently observed for passive-scalar mixing and for map dynamos.

We report a series of numerical simulations showing that the critical magnetic Reynolds number Rm_c for the nonhelical small-scale dynamo depends on the Reynolds number Re. Namely, the dynamo is shut down if the magnetic Prandtl number Pr=Rm/Re is less than some critical value Pr_c<1 even for Rm for which dynamo exists at Pr>=1. We argue that, in the limit of Re->infinity, a finite Pr_c may exist. The second possibility is that Pr_c->0 as Re->infinity, while Rm_c tends to a very large constant value inaccessible at current resolutions. If there is a finite Pr_c, the dynamo is sustainable only if magnetic fields can exist at scales smaller than the flow scale, i.e., it is always effectively a large-Pr dynamo. If there is a finite Rm_c, our results provide a lower bound: Rm_c<220 for Pr<=1/8. This is larger than Rm in many planets and in all liquid-metal experiments.

The growth and saturation of magnetic field in conducting turbulent media with large magnetic Prandtl numbers are investigated. This regime is very common in low-density hot astrophysical plasmas. During the early (kinematic) stage, weak magnetic fluctuations grow exponentially and concentrate at the resistive scale, which lies far below the hydrodynamic viscous scale. The evolution becomes nonlinear when the magnetic energy is comparable to the kinetic energy of the viscous-scale eddies. A physical picture of the ensuing nonlinear evolution of the MHD dynamo is proposed. Phenomenological considerations are supplemented with a simple Fokker--Planck model of the nonlinear evolution of the magnetic-energy spectrum. It is found that, while the shift of the bulk of the magnetic energy from the subviscous scales to the velocity scales may be possible, it occurs very slowly -- at the resistive, rather than dynamical, time scale (for galaxies, this means that generation of large-scale magnetic fields cannot be explained by this mechanism). The role of Alfvenic motions and the implications for the fully developed isotropic MHD turbulence are discussed.

We study the intermittency and field-line structure of the MHD turbulence in plasmas with very large magnetic Prandtl numbers. In this regime, which is realized in the interstellar medium, some accretion disks, protogalaxies, galaxy-cluster gas, early Universe, etc., magnetic fluctuations can be excited at scales below the viscous cutoff. The salient feature of the resulting small-scale magnetic turbulence is the folded structure of the fields. It is characterized by very rapid transverse spatial oscillation of the field direction, while the field lines remain largely unbent up to the scale of the flow. Quantitatively, the fluctuation level and the field-line geometry can be studied in terms of the statistics of the field strength and of the field-line curvature. In the kinematic limit, the distribution of the field strength is an expanding lognormal, while that of the field-line curvature K is stationary and has a power tail K^{-13/7}. The field strength and curvature are anticorrelated, i.e. the growing fields are mostly flat, while the sharply curved fields remain relatively weak. The field, therefore, settles into a reduced-tension state. Numerical simulations demonstrate three essential features of the nonlinear regime. First, the total magnetic energy is equal to the total kinetic energy. Second, the intermittency is partially suppressed compared to the kinematic case, as the fields become more volume-filling and their distribution develops an exponential tail. Third, the folding structure of the field is unchanged from the kinematic case: the anticorrelation between the field strength and the curvature persists and the distribution of the latter retains the same power tail. We propose a model of back reaction based on the folding picture that reproduces all of the above numerical results.

This is a brief review of the main results of our recent studies of the nonlinear evolution of the small-scale MHD dynamo in the high-Prandtl-number regime and of the structure of the resulting saturated state of the isotropic homogeneous MHD turbulence. It is emphasized that the MHD regime without a uniform mean field (as is the case in our studies) is fundamentally different from the one in which such a field is externally imposed. The ability of the turbulence to bend and fold the magnetic-field lines leads to the emergence of a distinctive small-scale structure. The fields are organized in folds of characteristic length comparable to the size of the largest turbulent eddies with spatial-direction reversals at the resistive scale. These folds are very hard to destroy. In the nonlinear regime, the folding structure coexists with Alfven waves propagating along the folds. The turbulent energy injected by the forcing is dissipated in part resistively via the small-scale magnetic fields, and in part viscously via the Alfven waves.

The effect of magnetic shear and shear flow on local gravitationally induced instabilities is investigated. A simple model is constructed allowing for an arbitrary entropy gradient and a shear plasma flow in the Boussinesq approximation. A transformation to shearing magnetic coordinates achieves a model with plasma flow along the magnetic field lines where the coordinate lines are coincident with the field lines. The solution for the normal modes of the system depends on two parameters: the Alfven Mach number of the plasma flow and the entropy gradient. The behavior of the unstable normal modes of this system is summarized by a stability diagram. Important characteristics of this stability diagram are the following: magnetic shear is stabilizing and the entropy gradient must exceed a threshold value for unstable mode growth to occur; flow acts to suppress mode growth in a substantially unstable regime as expected, yet near marginal stability it can lessen the stabilizing effect of magnetic shear and enhance the growth rates of the instability; and, as the Alfven Mach number approaches one, the instability is completely stabilized. Analytical work is presented supporting the characteristics of the stability diagram and illuminating the physical mechanisms controlling the behavior of the model. The implications of this work for astrophysical and fusion applications and the potential for future research extending the results to include compressibility are discussed.