Results are presented of direct numerical simulations of incompressible, homogeneous magnetohydrodynamic turbulence without a mean magnetic field, subject to different mechanical forcing functions commonly used in the literature. Specifically, the forces are negative damping (which uses the large-scale velocity field as a forcing function), a nonhelical random force, and a nonhelical static sinusoidal force (analogous to helical ABC forcing). The time evolution of the three ideal invariants (energy, magnetic helicity and cross helicity), the time-averaged energy spectra, the energy ratios and the dissipation ratios are examined. All three forcing functions produce qualitatively similar steady states with regards to the time evolution of the energy and magnetic helicity. However, differences in the cross helicity evolution are observed, particularly in the case of the static sinusoidal method of energy injection. Indeed, an ensemble of sinusoidally-forced simulations with identical parameters shows significant variations in the cross helicity over long time periods, casting some doubt on the validity of the principle of ergodicity in systems in which the injection of helicity cannot be controlled. Cross helicity can unexpectedly enter the system through the forcing function and must be carefully monitored.

In this paper we investigate the dynamics of magnetic helicity in magnetohydrodynamic (MHD) turbulent flows focusing at scales larger than the forcing scale. Our results show a nonlocal inverse cascade of magnetic helicity, which occurs directly from the forcing scale into the largest scales of the magnetic field. We also observe that no magnetic helicity and no energy is transferred to an intermediate range of scales sufficiently smaller than the container size and larger than the forcing scale. Thus, the statistical properties of this range of scales, which increases with scale separation, is shown to be described to a large extent by the zero flux solutions of the absolute statistical equilibrium theory exhibited by the truncated ideal MHD equations.

We present a numerical and analytical study of incompressible homogeneous conducting fluids using a helical Fourier representation. We analytically study both small- and large-scale dynamo properties, as well as the inverse cascade of magnetic helicity, in the most general minimal subset of interacting velocity and magnetic fields on a closed Fourier triad. We mainly focus on the dependency of magnetic field growth as a function of the distribution of kinetic and magnetic helicities among the three interacting wavenumbers. By combining direct numerical simulations of the full magnetohydrodynamics equations with the helical Fourier decomposition we numerically confirm that in the kinematic dynamo regime the system develops a large-scale magnetic helicity with opposite sign compared to the small-scale kinetic helicity, a sort of triad-by-triad $\alpha$-effect in Fourier space. Concerning the small-scale perturbations, we predict theoretically and confirm numerically that the largest instability is achieved for the magnetic component with the same helicity of the flow, in agreement with the Stretch-Twist-Fold mechanism. Vice versa, in presence of a Lorentz feedback on the velocity, we find that the inverse cascade of magnetic helicity is mostly local if magnetic and kinetic helicities have opposite sign, while it is more nonlocal and more intense if they have the same sign, as predicted by the analytical approach. Our analytical and numerical results further demonstrate the potential of the helical Fourier decomposition to elucidate the entangled dynamics of magnetic and kinetic helicities both in fully developed turbulence and in laminar flows.

This paper examines the behavior of the dimensionless dissipation rate $C_{\varepsilon}$ for stationary and nonstationary magnetohydrodynamic (MHD) turbulence in presence of external forces. By combining with previous studies for freely decaying MHD turbulence, we obtain here both the most general model equation for $C_{\varepsilon}$ applicable to homogeneous MHD turbulence and a comprehensive numerical study of the Reynolds number dependence of the dimensionless total energy dissipation rate at unity magnetic Prandtl number. We carry out a series of medium to high resolution direct numerical simulations of mechanically forced stationary MHD turbulence in order to verify the predictions of the model equation for the stationary case. Furthermore, questions of nonuniversality are discussed in terms of the effect of external forces as well as the level of cross- and magnetic helicity. The measured values of the asymptote $C_{\varepsilon,\infty}$ lie between $0.193 \leqslant C_{\varepsilon,\infty} \leqslant 0.268$ for free decay, where the value depends on the initial level of cross- and magnetic helicities. In the stationary case we measure $C_{\varepsilon,\infty} = 0.223$.

Spectral transfer processes in magnetohydrodynamic (MHD) turbulence are investigated analytically by decomposition of the velocity and magnetic fields in Fourier space into helical modes. Steady solutions of the dynamical system which governs the evolution of the helical modes are determined, and a stability analysis of these solutions is carried out. The interpretation of the analysis is that unstable solutions lead to energy transfer between the interacting modes while stable solutions do not. From this, a dependence of possible interscale energy and helicity transfers on the helicities of the interacting modes is derived. As expected from the inverse cascade of magnetic helicity in 3D MHD turbulence, mode interactions with like helicities lead to transfer of energy and magnetic helicity to smaller wavenumbers. However, some interactions of modes with unlike helicities also contribute to an inverse energy transfer. As such, an inverse energy cascade for nonhelical magnetic fields is shown to be possible. Furthermore, it is found that high values of the cross-helicity may have an asymmetric effect on forward and reverse transfer of energy, where forward transfer is more quenched in regions of high cross-helicity than reverse transfer. This conforms with recent observations of solar wind turbulence. For specific helical interactions the relation to dynamo action is established.

A model equation for the Reynolds number dependence of the dimensionless dissipation rate in freely decaying homogeneous magnetohydrodynamic turbulence in the absence of a mean magnetic field is derived from the real-space energy balance equation, leading to $C_{\varepsilon}=C_{\varepsilon, \infty}+C/R_- +O(1/R_-^2))$, where $R_-$ is a generalized Reynolds number. The constant $C_{\varepsilon, \infty}$ describes the total energy transfer flux. This flux depends on magnetic and cross helicities, because these affect the nonlinear transfer of energy, suggesting that the value of $C_{\varepsilon,\infty}$ is not universal. Direct numerical simulations were conducted on up to $2048^3$ grid points, showing good agreement between data and the model. The model suggests that the magnitude of cosmological-scale magnetic fields is controlled by the values of the vector field correlations. The ideas introduced here can be used to derive similar model equations for other turbulent systems.

Ensemble averaged high resolution direct numerical simulations of reverse spectral transfer are presented, extending on the many single realization numerical studies done up to now. This identifies this type of spectral transfer as a statistical property of magnetohydrodynamic turbulence and thus permits reliable numerical exploration of its dynamics. The magnetic energy decay exponent from these ensemble runs has been determined to be $n_E = (0.47 \pm 0.03) + (13.9 \pm 0.8)/R_{\lambda}$ for initially helical magnetic fields. We show for the first time that even after removing the Lorentz force term in the momentum equation, thus decoupling it from the induction equation, reverse spectral transfer still persists. The induction equation is now linear with an externally imposed velocity field, thus amenable to numerous analysis techniques. A new door has opened for analyzing reverse spectral transfer, with various ideas discussed.