Weak Alfvenic turbulence in a periodic domain is considered as a mixed state of Alfven waves interacting with the two-dimensional (2D) condensate. Unlike in standard treatments, no spectral continuity between the two is assumed and indeed none is found. If the 2D modes are not directly forced, k^{-2} and k^{-1} spectra are found for the Alfven waves and the 2D modes, respectively, with the latter less energetic than the former. The wave number at which their energies become comparable marks the transition to strong turbulence. For imbalanced energy injection, the spectra are similar and the Elsasser ratio scales as the ratio of the energy fluxes in the counterpropagting Alfven waves. If the 2D modes are forced, a 2D inverse cascade dominates the dynamics at the largest scales, but at small enough scales, the same weak and then strong regimes as described above are achieved.

We investigate the anisotropy of Alfv\'enic turbulence in the inertial range of slow solar wind and in both driven and decaying reduced magnetohydrodynamic simulations. A direct comparison is made by measuring the anisotropic second-order structure functions in both data sets. In the solar wind, the perpendicular spectral index of the magnetic field is close to -5/3. In the forced simulation, it is close to -5/3 for the velocity and -3/2 for the magnetic field. In the decaying simulation, it is -5/3 for both fields. The spectral index becomes steeper at small angles to the local magnetic field direction in all cases. We also show that when using the global rather than local mean field, the anisotropic scaling of the simulations cannot always be properly measured.

Two-dimensional numerical simulations of the effect of background turbulence on 2D resistive magnetic reconnection are presented. For sufficiently small values of the resistivity ($\eta$) and moderate values of the turbulent power ($\epsilon$), the reconnection rate is found to have a much weaker dependence on $\eta$ than the Sweet-Parker scaling of $\eta^{1/2}$ and is even consistent with an $\eta-$independent value. For a given value of $\eta$, the dependence of the reconnection rate on the turbulent power exhibits a critical threshold in $\epsilon$ above which the reconnection rate is significantly enhanced.

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.

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.

This paper examines the consistency of the exact scaling laws for isotropic MHD turbulence in numerical simulations with large magnetic Prandtl numbers Pm and with Pm=1. The exact laws are used to elucidate the structure of the magnetic and velocity fields. Despite the linear scaling of certain third-order correlation functions, the situation is not analogous to the case of Kolmogorov turbulence. The magnetic field is adequately described by a model of stripy (folded) field with direction reversals at the resistive scale. At currently available resolutions, the cascade of kinetic energy is short-circuited by the direct exchange of energy between the forcing-scale motions and the stripy magnetic fields. This nonlocal interaction is the defining feature of isotropic MHD turbulence.

Decaying turbulence is studied numerically using as initial condition a random flow whose shell-integrated energy spectrum increases with wavenumber k like k^q. Alternatively, initial conditions are generated from a driven turbulence simulation by simply stopping the driving. It is known that the dependence of the decaying energy spectrum on wavenumber, time, and viscosity can be collapsed onto a unique scaling function that depends only on two parameters. This is confirmed using three-dimensional simulations and the dependence of the scaling function on its two arguments is determined.

The dimensionless kinetic energy dissipation rate C_epsilon is estimated from numerical simulations of statistically stationary isotropic box turbulence that is slightly compressible. The Taylor microscale Reynolds number Re_lambda range is 20 < Re_lambda < 220 and the statistical stationarity is achieved with a random phase forcing method. The strong Re_lambda dependence of C_epsilon abates when Re_lambda approx. 100 after which C_epsilon slowly approaches approx 0.5 a value slightly different to previously reported simulations but in good agreement with experimental results. If C_epsilon is estimated at a specific time step from the time series of the quantities involved it is necessary to account for the time lag between energy injection and energy dissipation. Also, the resulting value can differ from the ensemble averaged value by up to +-30%. This may explain the spread in results from previously published estimates of C_epsilon.

Forced turbulence simulations are used to determine the turbulent kinematic viscosity, nu_t, from the decay rate of a large scale velocity field. Likewise, the turbulent magnetic diffusivity, eta_t, is determined from the decay of a large scale magnetic field. In the kinematic regime, when the field is weak, the turbulent magnetic Prandtl number, nu_t/eta_t, is about unity. When the field is nonhelical, eta_t is quenched when magnetic and kinetic energies become comparable. For helical fields the quenching is stronger and can be described by a dynamical quenching formula.

At numerical resolutions around $512^3$ and above, three-dimensional energy spectra from turbulence simulations begin to show noticeably shallower spectra than $k^{-5/3}$ near the dissipation wavenumber (`bottleneck effect'). This effect is shown to be significantly weaker in one-dimensional spectra such as those obtained in wind tunnel turbulence. The difference can be understood in terms of the transformation between one-dimensional and three-dimensional energy spectra under the assumption that the turbulent velocity field is isotropic. Transversal and longitudinal energy spectra are similar and can both accurately be computed from the full three-dimensional spectra. Second-order structure functions are less susceptible to the bottleneck effect and may be better suited for inferring the scaling exponent from numerical simulation data.

In the past few years suggestions have emerged that the solar magnetic field might have a bi-helical component with oppositely polarized magnetic fields at large and small scales, and that the shedding of such fields may be crucial for the operation of the dynamo. It is shown that, if a bi-helical field is shed into the solar wind, positive and negative contributions of the magnetic helicity spectrum tend to mix and decay. Even in the absence of turbulence, mixing and decay can occur on a time scale faster than the resistive one provided the two signs of magnetic helicity originate from a single tube. In the presence of turbulence, positively and negatively polarized contributions mix rapidly in such a way that the ratio of magnetic helicity to magnetic energy is largest both at the largest scale and in the dissipation range. In absolute units the small scale excess of helical fields is however negligible.