Helical fluid and (Hall)-MHD turbulence: a brief review.

Helicity, a measure of the breakage of reflectional symmetry representing the topology of turbulent flows, contributes in a crucial way to their dynamics and to their fundamental statistical properties. We review several of their main features, both new and old, such as the discovery of bi-directional cascades or the role of helical vortices in the enhancement of large-scale magnetic fields in the dynamo problem. The dynamical contribution in magnetohydrodynamic of the cross-correlation between velocity and induction is discussed as well. We consider next how turbulent transport is affected by helical constraints, in particular in the context of magnetic reconnection and fusion plasmas under one- and two-fluid approximations. Central issues on how to construct turbulence models for non-reflectionally symmetric helical flows are reviewed, including in the presence of shear, and we finally briefly mention the possible role of helicity in the development of strongly localized quasi-singular structures at small scale. This article is part of the theme issue 'Scaling the turbulence edifice (part 2)'.

Rotating turbulence under "precession-like" perturbation.

The effects of changing the orientation of the rotation axis on homogeneous turbulence is considered. We perform direct numerical simulations on a periodic box of 1024(3) grid points, where the orientation of the rotation axis is changed (a) at a fixed time instant (b) regularly at time intervals commensurate with the rotation time scale. The former is characterized by a dominant inverse energy cascade whereas in the latter, the inverse cascade is stymied due to the recurrent changes in the rotation axis resulting in a strong forward energy transfer and large-scale structures that resemble those of isotropic turbulence.

Small-scale behavior of Hall magnetohydrodynamic turbulence.

Decaying Hall magnetohydrodynamic (HMHD) turbulence is studied using three-dimensional (3D) direct numerical simulations with grids up to 768(3) points and two different types of initial conditions. Results are compared to analogous magnetohydrodynamic (MHD) runs and both Laplacian and Laplacian-squared dissipative operators are examined. At scales below the ion inertial length, the ratio of magnetic to kinetic energy as a function of wave number transitions to a magnetically dominated state. The transition in behavior is associated with the advection term in the momentum equation becoming subdominant to dissipation. Examination of autocorrelation functions reveals that, while current and vorticity structures are similarly sized in MHD, HMHD current structures are narrower and vorticity structures are wider. The electric field autocorrelation function is significantly narrower in HMHD than in MHD and is similar to the HMHD current autocorrelation function at small separations. HMHD current structures are found to be significantly more intense than in MHD and appear to have an enhanced association with strong alignment between the current and magnetic field, which may be important in collisionless plasmas where field-aligned currents can be unstable. When hyperdiffusivity is used, a longer region consistent with a k(-7/3) scaling is present for right-polarized fluctuations when compared to Laplacian dissipation runs.

Long-time properties of magnetohydrodynamic turbulence and the role of symmetries.

Using direct numerical simulations with grids of up to 512(3) points, we investigate long-time properties of three-dimensional magnetohydrodynamic turbulence in the absence of forcing and examine in particular the roles played by the quadratic invariants of the system and the symmetries of the initial configurations. We observe that when sufficient accuracy is used, initial conditions with a high degree of symmetries, as in the absence of helicity, do not travel through parameter space over time, whereas by perturbing these solutions either explicitly or implicitly using, for example, single precision for long times, the flows depart from their original behavior and can either become strongly helical or have a strong alignment between the velocity and the magnetic field. When the symmetries are broken, the flows evolve towards different end states, as already predicted by statistical arguments for nondissipative systems with the addition of an energy minimization principle. Increasing the Reynolds number by an order of magnitude when using grids of 64(3)-512(3) points does not alter these conclusions. Furthermore, the alignment properties of these flows, between velocity, vorticity, magnetic potential, induction, and current, correspond to the dominance of two main regimes, one helically dominated and one in quasiequipartition of kinetic and magnetic energies. We also contrast the scaling of the ratio of magnetic energy to kinetic energy as a function of wave number to the ratio of eddy turnover time to Alfvén time as a function of wave number. We find that the former ratio is constant with an approximate equipartition for scales smaller than the largest scale of the flow, whereas the ratio of time scales increases with increasing wave number.

Anisotropy and nonuniversality in scaling laws of the large-scale energy spectrum in rotating turbulence.

Rapidly rotating turbulent flow is characterized by the emergence of columnar structures that are representative of quasi-two-dimensional behavior of the flow. It is known that when energy is injected into the fluid at an intermediate scale Lf, it cascades towards smaller as well as larger scales. In this paper we analyze the flow in the inverse cascade range at a small but fixed Rossby number, Rof≈0.05. Several numerical simulations with helical and nonhelical forcing functions are considered in periodic boxes with unit aspect ratio. In order to resolve the inverse cascade range with reasonably large Reynolds number, the analysis is based on large eddy simulations which include the effect of helicity on eddy viscosity and eddy noise. Thus, we model the small scales and resolve explicitly the large scales. We show that the large-scale energy spectrum has at least two solutions: one that is consistent with Kolmogorov-Kraichnan-Batchelor-Leith phenomenology for the inverse cascade of energy in two-dimensional (2D) turbulence with a ∼k⊥-5/3 scaling, and the other that corresponds to a steeper ∼k⊥-3 spectrum in which the three-dimensional (3D) modes release a substantial fraction of their energy per unit time to the 2D modes. The spectrum that emerges depends on the anisotropy of the forcing function, the former solution prevailing for forcings in which more energy is injected into the 2D modes while the latter prevails for isotropic forcing. In the case of anisotropic forcing, whence the energy goes from the 2D to the 3D modes at low wave numbers, large-scale shear is created, resulting in a time scale τsh, associated with shear, thereby producing a ∼k-1 spectrum for the total energy with the horizontal energy of the 2D modes still following a ∼k⊥-5/3 scaling.

Not much helicity is needed to drive large-scale dynamos.

Understanding the in situ amplification of large-scale magnetic fields in turbulent astrophysical rotators has been a core subject of dynamo theory. When turbulent velocities are helical, large-scale dynamos that substantially amplify fields on scales that exceed the turbulent forcing scale arise, but the minimum sufficient fractional kinetic helicity f(h,C) has not been previously well quantified. Using direct numerical simulations for a simple helical dynamo, we show that f(h,C) decreases as the ratio of forcing to large-scale wave numbers k(F)/k(min) increases. From the condition that a large-scale helical dynamo must overcome the back reaction from any nonhelical field on the large scales, we develop a theory that can explain the simulations. For k(F)/k(min)≥8 we find f(h,C)≲3%, implying that very small helicity fractions strongly influence magnetic spectra for even moderate-scale separation.

Alfvén waves and ideal two-dimensional Galerkin truncated magnetohydrodynamics.

We investigate numerically the dynamics of two-dimensional Euler and ideal magnetohydrodynamics (MHD) flows in systems with a finite number of modes, up to 4096(2), for which several quadratic invariants are preserved by the truncation and the statistical equilibria are known. Initial conditions are the Orszag-Tang vortex with a neutral X point centered on a stagnation point of the velocity field in the large scales. In MHD, we observe that the total energy spectra at intermediate times and intermediate scales correspond to the interactions of eddies and waves, E(T)(k)~k(-3/2). Moreover, no pseudodissipative range is visible for either Euler or ideal MHD in two dimensions. In the former case, this may be linked to the existence of a vanishing turbulent viscosity whereas in MHD, the numerical resolution employed may be insufficient. When imposing a uniform magnetic field to the flow, we observe a lack of saturation of the formation of small scales together with a significant slowing down of their equilibration, with however a cutoff independent partial thermalization being reached at intermediate scales.

Lagrangian-averaged model for magnetohydrodynamic turbulence and the absence of bottlenecks.

We demonstrate that, for the case of quasiequipartition between the velocity and the magnetic field, the Lagrangian-averaged magnetohydrodynamics (LAMHD) alpha model reproduces well both the large-scale and the small-scale properties of turbulent flows; in particular, it displays no increased (superfilter) bottleneck effect with its ensuing enhanced energy spectrum at the onset of the subfilter scales. This is in contrast to the case of the neutral fluid in which the Lagrangian-averaged Navier-Stokes alpha model is somewhat limited in its applications because of the formation of spatial regions with no internal degrees of freedom and subsequent contamination of superfilter-scale spectral properties. We argue that, as the Lorentz force breaks the conservation of circulation and enables spectrally nonlocal energy transfer (associated with Alfvén waves), it is responsible for the absence of a viscous bottleneck in magnetohydrodynamics (MHD), as compared to the fluid case. As LAMHD preserves Alfvén waves and the circulation properties of MHD, there is also no (superfilter) bottleneck found in LAMHD, making this method capable of large reductions in required numerical degrees of freedom; specifically, we find a reduction factor of approximately 200 when compared to a direct numerical simulation on a large grid of 1536;{3} points at the same Reynolds number.

Highly turbulent solutions of the Lagrangian-averaged Navier-Stokes alpha model and their large-eddy-simulation potential.

We compute solutions of the Lagrangian-averaged Navier-Stokes alpha - (LANS alpha ) model for significantly higher Reynolds numbers (up to Re approximately 8300 ) than have previously been accomplished. This allows sufficient separation of scales to observe a Navier-Stokes inertial range followed by a second inertial range specific to the LANS alpha model. Both fully helical and nonhelical flows are examined, up to Reynolds numbers of approximately 1300. Analysis of the third-order structure function scaling supports the predicted l3 scaling; it corresponds to a k-1 scaling of the energy spectrum for scales smaller than alpha. The energy spectrum itself shows a different scaling, which goes as k1. This latter spectrum is consistent with the absence of stretching in the subfilter scales due to the Taylor frozen-in hypothesis employed as a closure in the derivation of the LANS alpha model. These two scalings are conjectured to coexist in different spatial portions of the flow. The l3 [E(k) approximately k-1] scaling is subdominant to k1 in the energy spectrum, but the l3 scaling is responsible for the direct energy cascade, as no cascade can result from motions with no internal degrees of freedom. We demonstrate verification of the prediction for the size of the LANS alpha attractor resulting from this scaling. From this, we give a methodology either for arriving at grid-independent solutions for the LANS alpha model, or for obtaining a formulation of the large eddy simulation optimal in the context of the alpha models. The fully converged grid-independent LANS alpha model may not be the best approximation to a direct numerical simulation of the Navier-Stokes equations, since the minimum error is a balance between truncation errors and the approximation error due to using the LANS alpha instead of the primitive equations. Furthermore, the small-scale behavior of the LANS alpha model contributes to a reduction of flux at constant energy, leading to a shallower energy spectrum for large alpha. These small-scale features, however, do not preclude the LANS alpha model from reproducing correctly the intermittency properties of the high-Reynolds-number flow.

Cancellation exponent and multifractal structure in two-dimensional magnetohydrodynamics: direct numerical simulations and Lagrangian averaged modeling.

We present direct numerical simulations and Lagrangian averaged (also known as alpha model) simulations of forced and free decaying magnetohydrodynamic turbulence in two dimensions. The statistics of sign cancellations of the current at small scales is studied using both the cancellation exponent and the fractal dimension of the structures. The alpha model is found to have the same scaling behavior between positive and negative contributions as the direct numerical simulations. The alpha model is also able to reproduce the time evolution of these quantities in free decaying turbulence. At large Reynolds numbers, an independence of the cancellation exponent with the Reynolds numbers is observed.

Shell-to-shell energy transfer in magnetohydrodynamics. II. Kinematic dynamo.

We study the transfer of energy between different scales for forced three-dimensional magnetohydrodynamics turbulent flows in the kinematic dynamo regime. Two different forces are examined: a nonhelical Taylor-Green flow with magnetic Prandtl number P(M) = 0.4 and a helical ABC flow with P(M) = 1. This analysis allows us to examine which scales of the velocity flow are responsible for dynamo action and identify which scales of the magnetic field receive energy directly from the velocity field and which scales receive magnetic energy through the cascade of the magnetic field from large to small scales. Our results show that the turbulent velocity fluctuations in the inertial range are responsible for the magnetic field amplification at small scales (small-scale dynamo) while the large-scale field is amplified mostly due to the large-scale flow. A direct cascade of the magnetic field energy from large to small scales is also presented and is a complementary mechanism for the increase of the magnetic field at small scales. The input of energy from the inertial range velocity field into the small magnetic scales dominates over the energy cascade up to the wave number where the magnetic energy spectrum peaks. At even smaller scales, most of the magnetic energy input is from the cascading process.

Shell-to-shell energy transfer in magnetohydrodynamics. I. Steady state turbulence.

We investigate the transfer of energy from large scales to small scales in fully developed forced three-dimensional magnetohydrodynamics (MHD) turbulence by analyzing the results of direct numerical simulations in the absence of an externally imposed uniform magnetic field. Our results show that the transfer of kinetic energy from large scales to kinetic energy at smaller scales and the transfer of magnetic energy from large scales to magnetic energy at smaller scales are local, as is also found in the case of neutral fluids and in a way that is compatible with the Kolmogorov theory of turbulence. However, the transfer of energy from the velocity field to the magnetic field is a highly nonlocal process in Fourier space. Energy from the velocity field at large scales can be transferred directly into small-scale magnetic fields without the participation of intermediate scales. Some implications of our results to MHD turbulence modeling are also discussed.

Numerical solutions of the three-dimensional magnetohydrodynamic alpha model.

We present direct numerical simulations and alpha -model simulations of four familiar three-dimensional magnetohydrodynamic (MHD) turbulence effects: selective decay, dynamic alignment, inverse cascade of magnetic helicity, and the helical dynamo effect. The MHD alpha model is shown to capture the long-wavelength spectra in all these problems, allowing for a significant reduction of computer time and memory at the same kinetic and magnetic Reynolds numbers. In the helical dynamo, not only does the alpha model correctly reproduce the growth rate of magnetic energy during the kinematic regime, it also captures the nonlinear saturation level and the late generation of a large scale magnetic field by the helical turbulence.

Continuum analysis of an avalanche model for solar flares.

We investigate the continuum limit of a class of self-organized critical lattice models for solar flares. Such models differ from the classical numerical sandpile model in their formulation of stability criteria in terms of the curvature of the nodal field, and are known to belong to a different universality class. A fourth-order nonlinear hyperdiffusion equation is reverse engineered from the discrete model's redistribution rule. A dynamical renormalization-group analysis of the equation yields scaling exponents that compare favorably with those measured in the discrete lattice model within the relevant spectral range dictated by the sizes of the domain and the lattice grid. We argue that the fourth-order nonlinear diffusion equation that models the behavior of the discrete model in the continuum limit is, in fact, compatible with magnetohydrodynamics (MHD) of the flaring phenomenon in the regime of strong magnetic field and the effective magnetic diffusivity characteristic of strong MHD turbulence.