3D Magnetohydrodynamic (MHD) turbulent flows with initially magnetic and kinetic energies at equipartition spontaneously develop a magnetic excess (or residual energy), as well in numerical simulations and in the solar wind. Closure equations obtained in 1983 describe the residual spectrum as being produced by a dynamo source proportional to the total energy spectrum, balanced by a linear Alfv\'en damping term. A good agreement was found in 2005 with incompressible simulations; however, recent solar wind measurements disagree with these results. The previous dynamo-Alfv\'en theory is generalized to a family of models, leading to simple relations between residual and total energy spectra. We want to assess these models in detail against MHD simulations and solar wind data. The family of models is tested against compressible decaying MHD simulations with low Mach number, low cross-helicity, zero mean magnetic field, without or with expansion terms (EBM or expanding box model). A single dynamo-Alfv\'en model is found to describe correctly both solar wind scalings and compressible simulations without or with expansion. It is equivalent to the 1983-2005 closure equation but with critical balance of nonlinear turnover and linear Alfv\'en times, while the dynamo source term remains unchanged. The discrepancy with previous incompressible simulations is elucidated. The model predicts a linear relation between the spectral slopes of total and residual energies $m_R = -1/2 + 3/2 m_T$. Examining the solar wind data as in \cite{2013ApJ...770..125C}, our relation is found to be valid whatever the cross-helicity, even better so at high cross-helicity, with the total energy slope varying from $1.7$ to $1.55$.

The measure of the third-order structure function, Y, is employed in the solar wind to compute the cascade rate of turbulence. In the absence of a mean field B0=0, Y is expected to be isotropic (radial) and independent of the direction of increments, so its measure yields directly the cascade rate. For turbulence with mean field, as in the solar wind, Y is expected to become more two dimensional (2D), that is, to have larger perpendicular components, loosing the above simple symmetry. To get the cascade rate one should compute the flux of Y, which is not feasible with single-spacecraft data, thus measurements rely upon assumptions about the unknown symmetry. We use direct numerical simulations (DNS) of magneto-hydrodynamic (MHD) turbulence to characterize the anisotropy of Y. We find that for strong guide field B0=5 the degree of two-dimensionalization depends on the relative importance of shear and pseudo polarizations (the two components of an Alfv\'en mode in incompressible MHD). The anisotropy also shows up in the inertial range. The more Y is 2D, the more the inertial range extent differs along parallel and perpendicular directions. We finally test the two methods employed in observations and find that the so-obtained cascade rate may depend on the angle between B0 and the direction of increments. Both methods yield a vanishing cascade rate along the parallel direction, contrary to observations, suggesting a weaker anisotropy of solar wind turbulence compared to our DNS. This could be due to a weaker mean field and/or to solar wind expansion.

Forced, weak MHD turbulence with guide field is shown to adopt different regimes, depending on the magnetic excess of the large forced scales. When the magnetic excess is large enough, the classical perpendicular cascade with $5/3$ scaling is obtained, while when equipartition is imposed, an isotropic $3/2$ scaling appears in all directions with respect to the mean field (\cite{2010PhRvE..82b6406G} or GM10). We show here that the $3/2$ scaling of the GM10 regime is not ruled by a small-scale cross-helicity cascade, and propose that it is a 3D extension of a perpendicular weak Iroshnikov-Kraichnan (IK) cascade. We analyze in detail the structure functions in real space and show that they closely follow the critical balance relation both in the local frame and the global frame: we show that there is no contradiction between this and the isotropic $3/2$ scaling of the spectra. We propose a scenario explaining the spectral structure of the GM10 regime, that starts with a perpendicular weak IK cascade and extends to 3D by using quasi-resonant couplings. The quasi-resonance condition happens to reduce the energy flux in the same way as is done in the weak perpendicular cascade, so leading to a $3/2$ scaling in all directions. We discuss the possible applications of these findings to solar wind turbulence.

The nonlinear dynamics of magnetic helicity, $H^M$, which is responsible for large-scale magnetic structure formation in electrically conducting turbulent media is investigated in forced and decaying three-dimensional magnetohydrodynamic turbulence. This is done with the help of high resolution direct numerical simulations and statistical closure theory. The numerically observed spectral scaling of $H^M$ is at variance with earlier work using a statistical closure model [Pouquet et al., J. Fluid Mech. \textbf{77} 321 (1976)]. By revisiting this theory a universal dynamical balance relation is found that includes effects of kinetic helicity, as well as kinetic and magnetic energy on the inverse cascade of $H^M$ and explains the above-mentioned discrepancy. Considering the result in the context of mean-field dynamo theory suggests a nonlinear modification of the $\alpha$-dynamo effect important in the context of magnetic field excitation in turbulent plasmas.

Investigations of the inverse cascade of magnetic helicity are conducted with pseudospectral, three-dimensional direct numerical simulations of forced and decaying incompressible magnetohydrodynamic turbulence. The high-resolution simulations which allow for the necessary scale-separation show that the observed self-similar scaling behavior of magnetic helicity and related quantities can only be understood by taking the full nonlinear interplay of velocity and magnetic fluctuations into account. With the help of the eddy-damped quasi-normal Markovian approximation a probably universal relation between kinetic and magnetic helicities is derived that closely resembles the extended definition of the prominent dynamo pseudoscalar $\alpha$. This unexpected similarity suggests an additional nonlinear quenching mechanism of the current-helicity contribution to $\alpha$.

The inverse cascade of magnetic helicity in 3D-MHD turbulence is believed to be one of the processes responsible for large scale magnetic structure formation in astrophysical systems. In this work we present an exhaustive set of high resolution direct numerical simulations (DNS) of both forced and decaying 3D-MHD turbulence, to understand this structure formation process. It is first shown that an inverse cascade of magnetic helicity in small-scale driven turbulence does not necessarily generate coherent large-scale magnetic structures. The observed large-scale magnetic field, in this case, is severely perturbed by magnetic fluctuations generated by the small-scale forcing. In the decaying case, coherent large-scale structure form similar to those observed astronomically. Based on the numerical results the formation of large-scale magnetic structures in some astrophysical systems, is suggested to be the consequence of an initial forcing which imparts the necessary turbulent energy into the system, which, after the forcing shuts off, decays to form the large-scale structures. This idea is supported by representative examples e.g. cluster of galaxies.

We propose a phenomenological theory of incompressible magnetohydrodynamic turbulence in the presence of a strong large-scale magnetic field, which establishes a link between the known anisotropic models of strong and weak MHD turbulence We argue that the Iroshnikov-Kraichnan isotropic cascade develops naturally within the plane perpendicular to the mean field, while oblique-parallel cascades with weaker amplitudes can develop, triggered by the perpendicular cascade, with a reduced flux resulting from a quasi-resonance condition. The resulting energy spectrum $E(k_\parallel,k_\bot)$ has the same slope in all directions. The ratio between the extents of the inertial range in the parallel and perpendicular directions is equal to $b_{rms}/B_0$. These properties match those found in recent 3D MHD simulations with isotropic forcing reported in [R. Grappin and W.-C. M\"uller, Phys. Rev. E \textbf{82}, 26406 (2010)].

We quantify possible differences between turbulent dynamo action in the Sun and the dynamo action studied in idealized simulations. For this purpose we compare Fourier-space shell-to-shell energy transfer rates of three incrementally more complex dynamo simulations: an incompressible, periodic simulation driven by random flow, a simulation of Boussinesq convection, and a simulation of fully compressible convection that includes physics relevant to the near-surface layers of the Sun. For each of the simulations studied, we find that the dynamo mechanism is universal in the kinematic regime because energy is transferred from the turbulent flow to the magnetic field from wavenumbers in the inertial range of the energy spectrum. The addition of physical effects relevant to the solar near-surface layers, including stratification, compressibility, partial ionization, and radiative energy transport, does not appear to affect the nature of the dynamo mechanism. The role of inertial-range shear stresses in magnetic field amplification is independent from outer-scale circumstances, including forcing and stratification. Although the shell-to-shell energy transfer functions have similar properties to those seen in mean-flow driven dynamos in each simulation studied, the saturated states of these simulations are not universal because the flow at the driving wavenumbers is a significant source of energy for the magnetic field.

The dynamics of spectral transport in two-dimensional turbulent convection of electrically conducting fluids is studied by means of direct numerical simulations (DNS) in the frame of the magnetohydrodynamic (MHD) Boussinesq approximation. The system performs quasi-oscillations between two different regimes of small-scale turbulence: one dominated by nonlinear MHD interactions, the other governed by buoyancy forces. The self-excited change of turbulent states is reported here for the first time. The process is controlled by the ideal invariant cross-helicity, $H^\mathrm{C}=\int_S \mathrm{d}S \mathbf{v}\cdot\mathbf{b}$. The observations are explained by the interplay of convective driving with the nonlinear spectral transfer of total MHD energy and cross-helicity.