# Triad interactions and the bidirectional turbulent cascade of magnetic helicity

###### Abstract

Using direct numerical simulations we demonstrate that magnetic helicity exhibits a bidirectional turbulent cascade at high but finite magnetic Reynolds numbers. Despite the injection of positive magnetic helicity in the flow, we observe that magnetic helicity of opposite signs is generated between large and small scales. We explain these observations by carrying out an analysis of the magnetohydrodynamic equations reduced to triad interactions using the Fourier helical decomposition. Within this framework, the direct cascade of positive magnetic helicity arises through triad interactions that are associated with small scale dynamo action, while the occurrence of negative magnetic helicity at large scales is explained through triad interactions that are related to stretch-twist-fold dynamics and small scale dynamo action, which compete with the inverse cascade of positive magnetic helicity. Our analytical and numerical results suggest that the direct cascade of magnetic helicity is a finite magnetic Reynolds number effect that will vanish in the limit .

Postprint version of the manuscript published in Phys. Rev. Fluids 2, 054605 (2017)

## I Introduction

The existence of planetary and stellar magnetic fields is currently attributed to dynamo action moffatt78 (); proctorgilbert94 (). One of the theoretical arguments to explain the generation and preservation of magnetic fields in spatial scales much larger than the outer scales of planets and stars is the inverse cascade of magnetic helicity in magnetohydrodynamic (MHD) turbulent flows Frischetal75 (). Magnetic helicity is defined as the correlation between the magnetic field and the magnetic potential , i.e., , with the angular brackets denoting a spatial average unless indicated otherwise. Magnetic helicity is considered to play a critical role in the long-term evolution of stellar and planetary magnetic fields brandenburg09 () and hence it is important to understand its dynamics across scales in order to shed light on the saturation mechanisms of the large-scale magnetic fields of planets and stars.

Early studies using turbulent closure models Pouquet76 () and mean field theory Steenbeck66 (); moffatt78 (); Krause80 () have shown (within the framework of their approximations) that magnetic helicity cascades from small scales to large scales in agreement with the prediction from equilibrium statistical mechanics for the ideal MHD equations Frischetal75 (). In these cases, the inverse cascade of magnetic helicity was associated with the effect parker55 () of large scale dynamos, where the kinetic helicity (where denotes the velocity field and the vorticity field) generates opposite signs of magnetic helicity between large and small scales. This change in sign across scales is further supported by the magnetic helicity spectra of solar wind data Brandenburgetal11 (), which show at small wavenumbers and at large wavenumbers. One way to understand this change in sign across scales is the conceptual stretch-twist-fold (STF) mechanism Childress95 (); Vainshtein72 (). This mechanism proposes that the advection of magnetic field lines by a positive helical flow leads to a positive magnetic helicity at small scales and to a negative magnetic helicity at large scales since large scale magnetic field lines are twisted in the opposite direction. However, it is presently not clear if such a mechanism is directly associated with a non-linear cascade process Mininni11 ().

The inverse cascade of magnetic helicity has been verified by direct numerical simulations (DNSs) meneguzzietal81 (); Brandenburg01 (); Alexakis06 () at moderate values of magnetic Reynolds number , where the small-scale magnetic helicity is dissipated quite fast. In this case, in contrast with the theory, a direct cascade of magnetic helicity with smaller magnitude was also observed. This bidirectional cascade of , where an inverse and a direct cascade of magnetic helicity coexist, was observed recently even at high flows Mueller12 (); ld16 (). In the limit of high concerns have been raised about the effectiveness of the effect in generating large-scale magnetic fields with strong amplitude due to the detrimental feedback that fast-growing small-scale magnetic fields have on the growth rate of the large scales Vainshtein92 (); Cattaneo96 (). These concerns are supported in some sense by the bidirectional cascade of magnetic helicity, because the magnitude of the inverse cascade is limited by the existence of the residual direct cascade, and by the nonlocality of the inverse cascade of magnetic helicity in the statistically stationary regime Alexakis06 (); ld16 ().

In this paper, we focus on the bidirectional cascade of magnetic helicity and on the mechanism of magnetic helicity to generate opposite signs of helicity across the scales. By means of DNSs we inject mean magnetic helicity in our flows using a helical electromagnetic forcing, while the velocity field is forced by a nonhelical forcing. This is in contrast to dynamo studies, which typically force only the velocity field using a helical forcing. In dynamo studies, the kinetic helicity generates opposite signs of magnetic helicity between large and small scales. This sign change across scales makes the interpretation of statistics ambiguous. With our approach we attempt to have a dominant sign of magnetic helicity across the scales in order to avoid such ambiguities as much as possible. The interpretation of our numerical results is supported by analytical work on the triadic interactions of helical modes in MHD turbulence Linkmannetal16 (); Linkmann17b (). This analysis is based on the helical decomposition of Fourier modes Craya58 (); Lesieur72 (); Herring74 (), which has been an important tool to understand the cascade dynamics of helical flows in Navier-Stokes turbulence Constantin88 (); Waleffe92 (); Chen03 (); Biferale12 (); Sahoo15 (); Alexakis17 (). Here, the analysis of triadic interactions of helical modes demonstrates that it can also be a useful tool to understand further (a) the presence of the direct cascade of magnetic helicity to small scales and (b) the mechanism that generates opposite signs of magnetic helicity across the scales.

This paper is organized as follows. We begin with the description of the numerical method and the resulting database in Sec. II. Sections III and IV present the global and spectral dynamics from our numerical simulations, respectively. Our numerical results are interpreted in Sec. V in terms of the triadic interactions of helical modes, where we propose an explanation for the direct cascade of magnetic helicity and for the mechanism that generates opposite signs of magnetic helicity across the scales. We conclude in Sec. VI by summarizing our findings.

## Ii Numerical simulations

The present work focuses on the bidirectional cascade of magnetic helicity. Thus, we consider numerical simulation of three-dimensional MHD turbulent flows forced at intermediate wavenumbers in the absence of large scales condensates. Forcing at intermediate scales and aiming for a turbulent flow with high enough scale separation poses serious computational constraints, since a very large range of scales needs to be resolved. So, to circumvent the demanding scale separation requirement and to avoid the condensation of energy at large scales, we consider high-order dissipation terms acting at the small and the large scales, respectively. These dissipation terms effectively increase the extent of the inertial range simulated for a given resolution by reducing the range of scales over which dissipation is effective. Therefore, we numerically solve the equations

(1) |

where denotes the velocity field, the magnetic induction expressed in Alfvén units, the vorticity, the current density, and , with the pressure and and the external mechanical and electromagnetic forces, respectively. Energy is dissipated at the small scales by the terms proportional to and and at the large scales by and . The indices and specify the order of the Laplacian used. In order to obtain a large inertial range, we choose . In the absence of forcing and dissipation Eqs. (1) reduce to the ideal MHD equations, which conserve the total energy , the magnetic helicity and the cross-helicity . These conserved quantities are those that determine the turbulent cascades in our flows.

Equations (1) are solved numerically in a cubic periodic domain with sides of length using the standard pseudospectral method, which ensures that and . Full dealiasing is achieved by the rule and as a result the minimum and maximum wavenumbers are = 1 and = N/3, respectively, where is the number of grid points in each Cartesian coordinate. Further details of the code can be found in Refs. Yoffe12 (); BereraLinkmann14 ().

As it was mentioned above, in these simulations we choose to inject mean magnetic helicity in our flows, attempting to have a dominant sign of magnetic helicity across the scales in order to be able to interpret our results unambiguously in comparison to dynamo studies. In order to do this in a systematic way we force using a helical forcing and using a non-helical forcing. Here, we choose to force both and with the same forcing amplitude, i.e., , so that both quantities are dynamically important and has a nonlinear feedback on the flow through the Lorenz force. The forces and are constructed from a randomized superposition of eigenfunctions of the curl operator Brandenburg01 (); Mueller12 (); MalapakaMuller13 (), resulting in Gaussian distributed and correlated in time forces whose helicities and correlation can be exactly controlled. The specific random nature of the forces ensures that at steady state the total energy input rate is known a priori Novikov65 (). Therefore, can be used as a control parameter since the amplitudes of the external forces are input parameters. The external forces are chosen such that the net cross helicity in the flow is negligible by keeping the correlation . In summary, no and no are injected into the flow, while the injection of is maximal. Here, we choose to exclude the injection of cross-helicity in the flow because we want to focus on the cascade dynamics of magnetic helicity, which can be influenced by introducing correlations between the velocity and the magnetic field Frischetal75 (); Linkmannetal16 (). Finally, the initial magnetic and velocity fields are in equipartition with energy spectra peaked at the forcing wavenumber and zero helicities, i.e., .

For all simulations, we set the magnetic Prandtl number for the small scales and for the large scales to unity, viz. . The values of and are tuned such that the inverse cascade is damped before the largest scales of the system are excited while the values of and are such that is satisfied for all the simulations, where is the dissipation wavenumber BorueOrszag95 (). The magnetic Reynolds number is defined based on the control parameters of the problem as with .

The energy input and thus the dissipation rate is adjusted such that is kept constant while the scale separation increases, resulting in a set of simulations with the same Reynolds number. For simulations at with different Reynolds numbers we kept the energy input constant and we varied . All the necessary numerical parameters that make these simulations reproducible are given in Table 1 along with their total runtime of the simulations normalized by , a time scale defined based on the control parameters.

10 | 128 | 0.05 | 1 | 1 | 320 | ||

20 | 256 | 0.05 | 130 | ||||

40 | 512 | 0.05 | 2 | 2 | 100 | ||

10 | 128 | 0.05 | 1 | 1 | 90 | ||

10 | 256 | 0.05 | 1 | 1 | 110 |

## Iii Time evolution

Before discussing the interaction of the helical modes in order to understand further the dynamics of magnetic helicity across scales, we provide an overview of the statistics of the flows under study. In Fig. 0(a) we have plotted the time evolution of the normalized magnetic helicity

(2) |

which belongs to the range . For the flows have no magnetic helicity, while for the flows are fully dominated by magnetic helicity, which means that the Lorentz force is expected to be zero on average. At time , we set but instantaneously it reaches high values because of the direct injection of positive by the helical electromagnetic force. At steady state reaches a mean value of 0.7 for all the flows with scale separations , indicating that the flows are dominated to a large degree by positive mean magnetic helicity.

The strength of the direct and the inverse cascade of magnetic helicity at steady state can be quantified by the rate of dissipation in large scales and small scales as

(3) |

Then, the total magnetic helicity dissipation can be computed as the sum of the dissipation at large and small scales, viz., . In Fig. 0(b), we plot the time-series of the ratio for the three scale separations that we simulated (see Table 1). The values of can range from 0 to 1. For there is no inverse cascade of magnetic helicity, while for there is no direct cascade of magnetic helicity.

In the initial transient regime the amount of magnetic helicity that is transferred to the large scales increases monotonically until large scale dissipation takes place and the ratio saturates to the value of 0.9 (see Fig. 0(b)). This clearly shows that there is a bidirectional cascade of with 90% of the magnetic helicity cascading to the large scales at steady state while 10% of the magnetic helicity cascades towards the small scales since . The bidirectional cascade of Hb had been observed qualitatively in Refs. Alexakis06 (); Mueller12 () and it was quantified in Ref. ld16 (). As scale separation increases the amount of magnetic helicity that cascades to the large and small scales remains fixed in contrast to the cascade of the total energy, where as it was shown in ld16 (), with the total energy dissipation rate being and . Thus, for we expect the ratio , implying that the total energy will cascade only toward the small scales while the direct cascade of magnetic helicity will not vanish because is expected to remain finite as our numerical simulations suggest.

However, for fixed it is plausible that the direct cascade of magnetic helicity vanishes in the high-magnetic-Reynolds-number limit. This can be inferred by considering the following upper bound derived in Ref. da15 () for the magnitude of the total dissipation rate of magnetic helicity at small scales for the MHD equations that involve dissipation terms with Laplacian operators (i.e., the exponent )

(4) |

where is the total energy dissipation rate in this case, is the magnetic resistivity, and is the kinematic viscosity. In high-Reynolds-number MHD turbulence, becomes finite and independent of and Mininni09 (); Dallas14b (); Linkmann17a (). Thus, in the limit of Eq. (4) suggests that , unless , implying infinite magnetic energy.

A similar inequality for the dissipation rate of magnetic helicity at small scales for Eqs. (1), which involve higher-order dissipation terms, can be derived as

(5) |

Note that for we recover Eq. (4), however, for the term cannot be directly related to the small-scale dissipation rate of magnetic energy . Therefore, to make any conclusions, if any, for the behavior of in the limit of , we need to compute the scaling of with as we cannot make any a priori statements like in Eq. (4). For this purpose, we performed three simulations with the same scale separation but different . These preliminary results are plotted in Fig. 2, which shows that normalized by the total dissipation of magnetic helicity has a weak power-law behavior in terms of , i.e. . This scaling suggests that the direct cascade of magnetic helicity may not persist in the limit of . However, further investigation into this issue is necessary in order to understand this scaling law and the small-scale dynamics of the magnetic helicity.

## Iv Spectral dynamics

After considering the global dynamics and the time evolution of our flows, we now elaborate on the spectral dynamics of magnetic helicity in order to have a clear picture of the dynamics across scales. This is best demonstrated by looking at the spectrum of magnetic helicity and its normalized flux for the flow with scale separation (see Fig. 3). The magnetic helicity flux is defined as

(6) |

and the magnetic helicity spectrum as , where the caret denotes the Fourier modes of the corresponding real vector fields and the asterisk the complex conjugate.

The gray curves denote instantaneous profiles at the transient regime of the flow while the green (dark gray) curves denote the time-averaged profiles at the steady state regime. Negative values of the magnetic helicity flux imply upscale transfer while positive values imply downscale transfer. Thus, Fig. 2(a) shows clearly the bidirectional cascade behavior of magnetic helicity with most of cascading towards large scales in agreement with Fig. 0(b). The gray curves indicate how the flux builds up at the transient state and of course these curves do not satisfy the balance , which is only expected to be valid on average over a statistically stationary regime.

At the transient stage of the simulation the inverse and the direct cascade of magnetic helicity is dominated by local interactions, as it has already been reported Alexakis06 (). Before the flow saturates to a steady-state solution a spectrum with a negative power-law slope can be observed at low wave numbers. The scaling of this spectrum was proposed to be Mueller12 (); Mininni09 (), which agrees with our data. However, if one integrates further in time the spectrum develops even more as the flow saturates to a statistically stationary regime. In this regime, most of the magnetic helicity is concentrated at the largest scales [see also the spectrum of normalized helicity in Fig. 5(a)] and the inverse cascade of becomes nonlocal in wave-number space, i.e., is transferred directly from the forced scale to the largest scales of the flow, while the forward cascade remains local ld16 (); Alexakis06 (); Brandenburg01 (). This transition to nonlocal interactions at steady state occurs when condensation takes place at the largest scale of the system. In this case, the large-scale vortex interacts directly with the small-scale vortices of the flow (i.e., nonlocally in wave-number space) ld16 (). This behavior is also observed in confined two-dimensional hydrodynamic turbulent flows shuklaetal16 (); xiaetal08 (). Finally, the scaling of the time-averaged spectrum of magnetic helicity at steady state is for wave numbers , suggesting that the large scales behave as if they were in absolute equilibrium ld16 (). Now, the scaling for seems to be . However, this is not conclusive from the plot due to the limited spectral range and the high-order dissipative terms that act at the largest wave numbers.

The behavior of the magnetic helicity flux at steady-state with increasing is presented in Fig. 4. As can be seen from the figure, the small-scale range where extends over a larger range of wave numbers successively with increasing as expected for a cascade process. However, the height of the plateau, i.e., the value of , decreases with increasing , as quantified in Fig. 2.

The purpose of our simulations using a positively helical electromagnetic force was to impose a mean magnetic helicity in our flows with the intention to have a dominant sign of magnetic helicity across the scales. However, at the steady-state regime of our simulations there are low-wave-number modes that develop a negative sign of magnetic helicity even though the mean value of magnetic helicity is positive in the flow. This is clearly depicted in Fig. 4(a), where we plot the time-averaged spectra of the normalized magnetic helicity for flows with .

Note that all the flows generate negative at wave number , which indicates that there is a consistent mechanism generating opposite signs of helicity across the scales even in these flows that are dominated by positive net magnetic helicity. Moreover, we observe that the number of modes at with negative magnetic helicity increases as the scale separation increases.

Despite the nonhelical mechanical forcing, a whole spectrum of positive kinetic helicity is observed across the scales. This is shown in Fig. 4(b), where we plot the spectra of the relative kinetic helicity for the flows with and . The values of for wave numbers decrease as the scale separation increases, while they are particularly high close to the forcing scale where is also dominant. This correlation seems to be related to the injection of positive magnetic helicity in the flow. In order to understand further the observations from our numerical simulation on the generation of the spectrum, the direct cascade of magnetic helicity, and the mechanism that generates opposite signs of magnetic helicity between large and small scales, we study analytically the triadic interactions of helical modes in the following section.

## V Triad interactions of helical modes

The three-dimensional vector fields of the velocity and the magnetic field are solenoidal (viz., and ), hence and must lie in the plane perpendicular to the wave vector . This plane is spanned by two eigenvectors of the curl operator in Fourier space with nonzero eigenvalues. Since these eigenvectors are by definition fully helical, each Fourier mode and can be further decomposed into two modes with positive and negative helicity

(7) | ||||

(8) |

where the basis vectors
are the
orthonormal eigenvectors of the curl operator in Fourier space satisfying
with and .
This is the so-called helical decomposition Craya58 (); Lesieur72 (); Herring74 ().
To study the triad interaction of helical modes, we use the formalism that was
developed by Waleffe Waleffe92 () for the Navier-Stokes equations and extended by
Linkmann et al. Linkmannetal16 (); Linkmann17b ()
to the MHD equations.
This is essentially a linear stability analysis of a
dynamical system obtained from the MHD equations in the limits of and .
In the MHD equations that are shown here in Fourier space

(9) | ||||

(10) |

with denoting the Fourier transform as a linear operator, the convolutions that describe the respective Fourier transforms of the inertial term , the Lorentz force , and the curl of the electromotive force are reduced to single triads of wave vectors , and satisfying . For convenience and without loss of generality we impose the ordering . In order to study the dynamics in the inertial range of scales, the dissipation terms are neglected in the limits of and . Following Lessinnes et al. Lessinnes09 (), by substituting Eqs. (7) and (8) into Eqs. (9) and (10) after restricting the convolutions to single triads and subsequently taking the inner product with the helical basis vectors, one can then derive the following system of ordinary differential equations, which conserves the ideal invariants of the MHD equations Lessinnes09 ():

(11) |

where the geometric factors describe the coupling between the helical basis vectors, i.e. Waleffe92 (). Further details on the full derivation of Eqs. (V) can be found in Refs. Lessinnes09 (); Linkmannetal16 (); Linkmann17b (). A graphical representation of the system of equations (V) is shown in Fig. 6, where each convolution term corresponds to a triangle that represents a triad of wave vectors. Note that two triads are required to describe the interactions that correspond to the term due to a necessary symmetrization of the convolution in Fourier space. However, this symmetrization does not allow one to disentangle the triadic interactions of the advection term and the stretching term of the magnetic field Linkmannetal16 (). Different interactions of helical modes can now be studied via Eqs. (V) by choosing different combinations of and .

A linear stability analysis of steady solutions of Eqs. (V) has been carried out by Linkmann et al. Linkmannetal16 (), where a linear instability corresponds to a transfer of energy from the unstable helical mode (denoted in capital letters) of a single triad interaction to the perturbations (denoted in lower case letters), i.e., to the other two helical modes of the triad. Under the assumption that the statistical behavior of the flow is controlled by the stability characteristics of these isolated triads (referred to as the instability assumption) Kraichnan67 (); Waleffe92 (), it is possible to draw conclusions concerning the energy and helicity transfer in the MHD equations from the stability properties of Eqs. (V). Different physical processes such as the inverse cascade of magnetic helicity or kinematic dynamo action can then be studied on the level of triad interactions by setting up specific perturbation problems Linkmannetal16 (). Depending on the characteristic wave numbers and the sign of helicities of the interacting modes, the dominant interscale energy transfers can then be identified through a comparison of the growth rates of the perturbations Linkmann17b ().

If we consider a steady solution for the velocity field subject to magnetic perturbations, it is possible to identify all the triadic interactions of helical modes that produce small as well as large scale growth of the magnetic field (see Ref. Linkmannetal16 () for details). With this result Linkmann et al. Linkmann17b () were able to interpret how small- and large-scale dynamos operate at the level of triadic interactions. In Fig. 7 we present all the triads that lead to linear instabilities and hence to energy transfer,

where the red arrows indicate large-scale dynamo action (i.e., energy is transferred to small-wavenumber helical modes of the magnetic field), the blue arrows indicate small-scale dynamo action (i.e., energy is transferred to large-wavenumber helical modes of the magnetic field), and the thickness of the arrows indicates the magnitude of the transfer (see Refs. Linkmannetal16 (); Linkmann17b () for details). More precisely, the thickness of the arrows is qualitative, reflecting that for any nontrivial triad geometry the growth rates of the perturbations have a consistent ordering [see, e.g., Eq. (20)].

Similarly, if we consider the linear stability analysis of a steady solution for the magnetic field subject to velocity and magnetic perturbations, one can obtain the triadic interactions involving unstable helical modes shown in Fig. 8,

where the green arrows indicate magnetic energy transfer, the magenta arrows indicate the conversion of magnetic to kinetic energy due to the action of the Lorentz force, and the thickness of the arrows indicates the magnitude of the transfer. As can be seen in Fig. 8, triadic interactions of helical modes leading to energy transfer occur only if (a) the steady solution of the magnetic field and the magnetic perturbation have the same sign of magnetic helicity, and (b) the characteristic wave number of the magnetic perturbation is smaller than the characteristic wavenumber of the steady solution. In other words, the energy transfer between the magnetic modes occurs exclusively from large to small wavenumbers and only between magnetic modes with the same sign of helicity, which implies an inverse cascade of magnetic helicity Linkmannetal16 (). As indicated by the thickness of the green arrows, the inverse cascade of magnetic helicity is stronger if the magnetic and kinetic helicities are of the same sign Linkmann17b (). Following the analysis that was carried out in Ref. Linkmann17b () for the triadic interactions corresponding to kinematic dynamo action and the inverse transfer of magnetic helicity, we carry out an analysis of the relative magnitude of the growth rates corresponding to triad interactions of the Lorentz force in the Appendix. Our analysis shows that the conversion of magnetic to kinetic energy through the Lorentz force occurs mainly between helical modes of the velocity and the magnetic field that have the same sign of helicity, as indicated by the thick magenta arrows in Fig. 8.

### v.1 Direct cascade of magnetic helicity

In order to understand the presence of the direct cascade of positive magnetic helicity, observed in our simulations, at the level of triadic interactions, we aim to identify interactions among helical modes that transfer energy to helical modes of positive magnetic helicity at large wavenumbers. The only possible way to get such interactions is via the triads in Figs. 6(a), 6(b), and 6(d). From these three triad interactions, those that dominate the small scales are shown in Figs. 6(a) and 6(b), where the unstable helical modes and transfer energy to via transfers that represent a small scale dynamo. In all these cases, energy is transferred from a small-wave-number mode to a large-wave-number mode and this transfer originates from a positively helical unstable mode of the velocity field. As we have already observed from our simulations, kinetic helicity has a well-defined spectrum [see Fig. 4(b)] and thus positively helical modes of the velocity field are present and can enable such triadic interactions.

The generation of kinetic helicity in these flows can be understood via the triadic interactions of Fig. 7(a)-7(c). These three triads are the dominant interactions in this linear stability analysis and in all these cases modes of positive magnetic helicity generate modes of positive kinetic helicity. In other words, the analysis of triadic interactions suggests that the magnetic helicity generates kinetic helicity through triad interactions associated with the Lorentz force. As discussed in Sec. III, our flows are dominated by positively helical modes of the magnetic field because we maintain positive net magnetic helicity [see Fig. 0(a)] via the positively helical electromagnetic force .

In summary, since there is no direct transfer of energy between modes that have opposite sign of magnetic helicity, the only way for a direct cascade of to occur is via the Lorentz force, which generates modes with positive kinetic helicity. These in turn excite modes with both positive and negative magnetic helicity with the dominant interactions creating positively helical magnetic-field modes at small scales (see Figs. 6(a), 6(b), and 6(d)).

### v.2 Negative magnetic helicity at large scales

Here we discuss the triad interactions that lead to the generation of modes with negative magnetic helicity at small wave numbers. As stated above, there is no direct transfer of energy between modes that have opposite sign of magnetic helicity. Therefore, the triad interactions that generate opposite signs of magnetic helicity between large and small scales have to involve a transfer of energy from positively helical modes of the velocity field. According to the results of the stability analysis of the triadic interactions shown in Figs. 7 and 8, we can conclude that there are two types of instabilities leading to the occurrence of large-scale magnetic helicity, namely, triadic instabilities that represent large-scale dynamo, including STF dynamics [see Figs. 6(b) and 6(c)], and those that can be directly associated with the inverse cascade of positive magnetic helicity [see Figs. 7(b) and 7(c)]. Here we expect a competition between triad interactions in Figs. 6(b) and 6(c) and those in Figs. 7(b) and 7(c). Based on our numerical simulations, we observe the persistent generation of negative magnetic helicity at wave number [see Fig. 4(a)] despite the positive mean value of magnetic helicity that is maintained by the positively helical electromagnetic forcing. This observation indicates that the triad interactions of Figs. 6(b) and 6(c) dominate over Figs. 7(b) and 7(c) at low wavenumbers as the scale separation increases. Thus, the generation of opposite signs of magnetic helicity between large and small scales can only be partially explained by STF dynamics, which is essentially what the triad in Fig. 6(b) represents. This is because the triadic interactions of the helical modes in Fig. 6(c), which represent a small scale dynamo, can play an important role. Moreover, the persistence of the triads in Figs. 6(b) and 6(c) can also explain why the scales larger that the forcing scale do not become fully helical and hence the relative helicity at wavenumbers .

## Vi Conclusions

In this paper we investigated the bidirectional cascade of magnetic helicity in homogeneous MHD turbulence by means of direct numerical simulations and analysis of triad interactions of helical modes. In order to be able to simulate statistically stationary high- flows with large enough scale separation we considered high-order dissipation terms acting at small and large scales of the velocity and the magnetic field. To maintain positive mean magnetic helicity in our flows we used a helical electromagnetic forcing, while the velocity field was forced by a nonhelical forcing.

Our numerical simulations demonstrated the following main results. (a) Magnetic helicity exhibits a clear bidirectional cascade even in these high- flows. This bidirectional cascade consists of a dominant inverse cascade and a residual direct cascade of smaller magnitude, which remains invariant as the scale separation increases. (b) Despite the mean positive magnetic helicity measured in all our simulations, we consistently observe the generation of opposite signs of magnetic helicity between large and small scales.

Our numerical results were interpreted in terms of the triad interactions of the helical Fourier modes. The residual forward cascade of magnetic helicity may occur due to a nonlinear small-scale dynamo, while the occurrence of large-scale negative helicity may result from a competition between triadic interactions associated with dynamo action and those corresponding to the inverse cascade of magnetic helicity.

The analysis of triadic interactions of helical modes predicts that the kinetic helicity in our flows is generated at all scales by the unstable modes of positive magnetic helicity. This is deduced from the stability properties of triadic interactions associated with the Lorentz force, and the predictions are in agreement with our numerical simulations, where kinetic helicity becomes significant by developing a whole spectrum of predominantly positive helicity across the scales. When the modes of positive kinetic helicity are unstable they preferentially generate modes of positive magnetic helicity at small scales Linkmannetal16 (); Linkmann17b (). This explains the existence of the residual direct cascade of magnetic helicity, which is associated with a triad interaction representing the Lorentz force and a triad that represents STF dynamics (i.e. the generation of opposite signs of magnetic helicity between between large and small scales due to a helical velocity field). Moreover, since triad interactions between modes that have opposite sign of magnetic helicity are prohibited, the only way for the magnetic helicity to have opposite signs between large and small scales is via the combined effect of triad interactions associated with STF dynamics and small-scale dynamo action, which involve unstable modes of positive kinetic helicity.

A large network of triadic interactions given by the MHD equations can behave differently from a collection of isolated systems of triads Moffatt14 (). For example, we cannot comment just from the analysis of triad interactions whether the direct cascade of magnetic helicity will vanish at the limit . Therefore, numerical investigations of the effect of different couplings of helical modes in DNSs like in Refs. Biferale12 (); Alexakis17 () would complement our analysis on the interscale transfers of magnetic helicity.

Preliminary results from numerical simulations with fixed scale separation and increasing suggest that the direct cascade of magnetic helicity may be a finite-Reynolds-number effect. This is because the scaling of the normalized dissipation rate of at small scales is . This weak power law implies that the direct cascade of magnetic helicity will persist as increases, but it will eventually vanish in the limit of . This claim is supported by deriving upper bounds on the absolute value of the dissipation rate of at small scales. However, further investigation is required in order to understand better the small scale dynamics of magnetic helicity. In principle, this could be achieved through a series of high-resolution DNSs of helical MHD turbulence covering a large range of magnetic Reynolds numbers from which one could extrapolate towards the high-Reynolds-number asymptotic limit. However, such a series of simulations will be computationally challenging as sufficient scale separation is required between the forcing scale and the box size, while the use of the Laplacian operator requires a large number of collocation points in order to achieve adequate small scale resolution.

Magnetic helicity is considered to play a fundamental role in the dynamo action of planetary, stellar, and galactic magnetic fields. The time scale and the level of the saturation of the magnetic field in models of the solar activity cycle depend on assumptions about the diffusion of magnetic helicity by small-scale velocity fluctuations kleeorinetal03 (). Our results improve our understanding of the magnetic helicity dynamics on the formation and evolution of large-scale magnetic fields in astrophysical objects, which is one of the challenging problems in modern astrophysical fluid dynamics.

Note added. Recently, it has come to our attention that similar conclusions concerning the behavior of the direct cascade of magnetic helicity in the limit of infinite have been reached Aluie17 ().

## Acknowledgements

We thank Luca Biferale and Ganapati Sahoo for helpful conversations. This work has made use of the resources provided by ARCHER, made available through the Edinburgh Compute and Data Facility (ECDF). V. D. acknowledges support from the Royal Society and the British Academy of Sciences (Newton International Fellowship, Grant No. NF140631). The research leading to these results has received funding from the European Union’s Seventh Framework Programme under grant agreement No. 339032.

## Appendix A Triadic analysis of the Lorentz force

Following the work of Refs. Linkmannetal16 (); Linkmann17b (), we present here the linear stability analysis of the triadic interactions that correspond
to the Lorentz force as an example. This analysis allows us to derive the additional results concerning
the magnitude of the magnetic-to-kinetic energy transfers depicted in
Fig. 8. We consider a positively helical magnetic steady solution
of Eqs. (V)
at wave vector
subject to magnetic and velocity field perturbations , , , and ,

where each set of equations represents a particular helical combination of velocity- and magnetic-field modes from the system of equations (V) for a particular choice of helicities. The procedure is conceptually similar to a linear stability analysis of rigid body rotation Waleffe92 (), that is, we differentiate Eqs. (12)-(A) further in time and replace any occurrence of a first-order time derivative by Eqs. (V). In order to isolate the dynamical effects of the Lorentz force on the flow we focus on the evolution of the velocity perturbations

(16) | ||||

(17) | ||||

(18) | ||||

(19) |

and we observe that Eqs. (18) and (19) cannot have exponentially growing solutions as the term is always negative. Concerning Eqs. (16) and (17), the wavenumber ordering results in , hence these two equations admit exponentially growing solutions. In other words, linear instabilities occur only in Eqs. (16) and (17).

From Eqs. (12) and (13) we observe that the coefficient corresponds to and coupling to . Therefore, the term in Eq. (16) describes the coupling of magnetic-field modes with the same signs of helicity, i.e., and to . Similarly, the term in Eq. (18) describes the coupling of magnetic-field modes with opposite signs of helicity to , because couples to and . Equations. (17) and (19) are analyzed analogously. Hence linear instabilities only occur in triadic interactions where the unstable magnetic mode and the magnetic perturbation have the same sign of helicity Linkmannetal16 (). In terms of the instability assumption this result implies that the Lorentz force can only convert magnetic to kinetic energy if the interaction proceeds by triads involving magnetic-field modes with the same sign of helicity Linkmannetal16 ().

The energy transfer due the linear instability governed by Eq. (16) corresponds to the magenta arrow in Fig. 7(b), while the linear instability described by Eq. (18) corresponds to the magenta arrow in Fig. 7(f). The remaining energy transfers depicted in Fig. 8 can be derived analogously by changing the characteristic wavenumber of the steady magnetic-field mode. The magnitude of the transfers indicated by the thickness of the arrows in Fig. 8 is determined by a comparison of the growth rate of the mechanical perturbations

(20) |

since the magnitude of the geometric factors can be written as Waleffe92 (); Linkmann17b (). Hence the energy transfer from a positively helical magnetic field into a positively helical flow is stronger than that occurring from a positively helical magnetic field into a negatively helical flow. This result is indicated by a thicker magenta arrow in Fig. 7(c) than in Fig. 7(f). The relative magnitudes of the remaining energy transfers connected with dynamo action as shown in Figs. 7 and the inverse cascade of magnetic helicity shown by the green arrows in Fig. 8 have been derived by the same method in Ref. Linkmann17b ().

## References

- [1] H. K. Moffatt. Magnetic Field Generation in Electrically Conducting Fluids. (Cambridge University Press, Cambridge, 1978).
- [2] M. R. E. Proctor and A. D. Gilbert. Lectures on Solar and Planetary Dynamos. (Cambridge University Press, Cambridge, 1994).
- [3] U. Frisch, A. Pouquet, J. Léorat, and A. Mazure. Possibility of an inverse cascade of magnetic helicity in magnetohydrodynamic turbulence. J. Fluid Mech., 68:769–778, 1975.
- [4] A. Brandenburg. The critical role of magnetic helicity in astrophysical large-scale dynamos. Plasma Physics and Contr. Fusion, 51(12):124043, 2009.
- [5] A. Pouquet, U. Frisch, and J. Léorat. Strong MHD helical turbulence and the nonlinear dynamo effect. J. Fluid Mech., 77:321–354, 1976.
- [6] M. Steenbeck, F. Krause, and K.-H. Rädler. Berechnung der mittleren Lorentz-Feldstärke v x B für ein elektrisch leitendes Medium in turbulenter, durch Coriolis-Kräfte beeinflußter Bewegung. Z. Naturforsch. A, 21:369, 1966.
- [7] F. Krause and K. Rädler. Mean-Field Magnetohydrodynamics and Dynamo Theory. (Pergamon, Oxford, 1980).
- [8] E. N Parker. Hydromagnetic dynamo models. The Astrophysical Journal, 122:293, 1955.
- [9] A. Brandenburg, K. Subramanian, A. Balogh, and M. L. Goldstein. Scale dependence of magnetic helicity in the solar wind. The Astrophysical Journal, 734(1):9, 2011.
- [10] S. Childress and A. D. Gilbert. Stretch, Twist, Fold: The Fast Dynamo. (Springer, Berlin, 1995).
- [11] S. I. Vainshtein and Y. B. Zeldovich. Origin of magnetic fields in astrophysics. Sov. Phys. Usp., 15:159–172, 1972.
- [12] P. D. Mininni. Scale interactions in magnetohydrodynamic turbulence. Annu. Rev. Fluid Mech., 43:377–397, 2011.
- [13] M. Meneguzzi, U. Frisch, and A. Pouquet. Helical and nonhelical turbulent dynamos. Phys. Rev. Lett., 47:1060–1064, 1981.
- [14] A. Brandenburg. The inverse cascade and nonlinear alpha-effect in simulations of isotropic helical magnetohydrodynamic turbulence. Astrophys. Journal, 550:824–840, 2001.
- [15] A. Alexakis, P. D. Mininni, and A. Pouquet. On the inverse cascade of magnetic helicity. Astrophys. J., 640:335–343, 2006.
- [16] W. C. Müller, S. K. Malapaka, and A. Busse. Inverse cascade of magnetic helicity in magnetohydrodynamic turbulence. Phys. Rev. E, 85:015302, 2012.
- [17] M. Linkmann and V. Dallas. Large-scale dynamics of magnetic helicity. Phys. Rev. E, 94:053209, 2016.
- [18] S. I. Vainshtein and F. Cattaneo. Nonlinear restrictions on dynamo action. Astrophys. J., 393:165–171, 1992.
- [19] F. Cattaneo and D. W. Hughes. Nonlinear saturation of the turbulent -effect. Phys. Rev. E, 54:R4532, 1996.
- [20] M. F. Linkmann, A. Berera, M. E. McKay and J. Jäger. Helical mode interactions and spectral energy transfer in magnetohydrodynamic turbulence. J. Fluid Mech., 791:61–96, 2016.
- [21] M. Linkmann, G. Sahoo, M. McKay, A. Berera, and L. Biferale. Effects of magnetic and kinetic helicities on the growth of magnetic fields in laminar and turbulent flows by helical-Fourier decomposition. Astrophys. Journal, 836:26, 2017.
- [22] A. Craya. Contribution à l’analyse de la turbulence associée à des vitesses moyennes. Ministere de l’Air report, 1958 (unpublished).
- [23] M Lesieur. Décomposition d’un champ de vitesse non divergent en ondes d’hélicité. Observatoire de Nice report, 1972 (unpublished).
- [24] J. R. Herring. Approach of axisymmetric turbulence to isotropy. Physics of Fluids, 17(5):859–872, 1974.
- [25] P. Constantin and A. Majda. The Beltrami spectrum for incompressible flows. Commun. Math. Phys., 115:435–456, 1988.
- [26] F. Waleffe. The nature of triad interactions in homogeneous turbulence. Phys. Fluids A, 4:350–363, 1992.
- [27] Q. Chen, S. Chen, and G. L. Eyink. The joint cascade of energy and helicity in three-dimensional turbulence. Phys. Fluids, 15:361, 2003.
- [28] L. Biferale, S. Musacchio, and F. Toschi. Inverse energy cascade in three-dimensional isotropic turbulence. Phys. Rev. Lett., 108:164501, 2012.
- [29] G. Sahoo, F. Bonaccorso, and L. Biferale. Role of helicity for large- and small-scale turbulent fluctuations. Phys. Rev. E, 92:051002, 2015.
- [30] Alexandros Alexakis. Helically decomposed turbulence. Journal of Fluid Mechanics, 812:752–770, 2017.
- [31] S. R. Yoffe. Investigation of the transfer and dissipation of energy in isotropic turbulence. PhD thesis, University of Edinburgh, 2012.
- [32] A. Berera and M. F. Linkmann. Magnetic helicity and the evolution of decaying magnetohydrodynamic turbulence. Phys. Rev. E, 90:041003(R), 2014.
- [33] S. K. Malapaka and W.-C. Müller. Large-scale magnetic structure formation in three-dimensional magnetohydrodynamic turbulence. Astrophys. J., 778:21, 2013.
- [34] E. A. Novikov. Functionals and the random-force method in turbulence theory. Soviet Physics JETP, 20:1290–1294, 1965.
- [35] V. Borue and S. A. Orszag. Self-similar decay of three-dimensional homogeneous turbulence with hyperviscosity. Phys. Rev. E, 51:2859(R), 1995.
- [36] V. Dallas and A. Alexakis. Self-organisation and non-linear dynamics in driven magnetohydrodynamic turbulent flows. Physics of Fluids, 27(4):045105, 2015.
- [37] P. D. Mininni and A. G. Pouquet. Finite dissipation and intermittency in magnetohydrodynamics. Phys. Rev. E., 80:025401, 2009.
- [38] V. Dallas and A. Alexakis. The signature of initial conditions on magnetohydrodynamic turbulence. The Astrophysical Journal, 788:L36, 2014.
- [39] M. Linkmann, A. Berera, and E. E. Goldstraw. Reynolds-number dependence of the dimensionless dissipation rate in homogeneous magnetohydrodynamic turbulence. Phys. Rev. E, 95:013102, 2017.
- [40] Vishwanath Shukla, Stephan Fauve, and Marc Brachet. Statistical theory of reversals in two-dimensional confined turbulent flows. Phys. Rev. E, 94:061101, 2016.
- [41] H. Xia, H. Punzmann, G. Falkovich, and M. G. Shats. Turbulence-condensate interaction in two dimensions. Phys. Rev. Lett., 101:194504, 2008.
- [42] T. Lessinnes, F. Plunian, and D. Carati. Helical shell models for MHD. Theor. Comput. Fluid Dyn., 23:439–450, 2009.
- [43] R. H. Kraichnan. Inertial ranges in two-dimensional turbulence. Phys. Fluids, 10(10):1417, 1967.
- [44] H. K. Moffatt. Note on the triad interactions of homogeneous turbulence. J. Fluid Mech., 741:R3, 2014.
- [45] N Kleeorin, K Kuzanyan, D Moss, I Rogachevskii, D Sokoloff, and H Zhang. Magnetic helicity evolution during the solar activity cycle: observations and dynamo theory. Astronomy & Astrophysics, 409(3):1097–1105, 2003.
- [46] H. Aluie. Coarse-Grained incompressible magnetohydrodynamics: Analyzing the turbulent cascades. New J. Phys., 19:025008, 2017.