Decaying Magnetohydrodynamic Turbulence Research Articles

Helios 2 observations of solar wind turbulence decay in the inner heliosphere

Aims. A linear scaling of the mixed third-order moment of the magnetohydrodynamic (MHD) fluctuations is used to estimate the energy transfer rate of the turbulent cascade in the expanding solar wind. Methods. In 1976, the Helios 2 spacecraft measured three samples of fast solar wind originating from the same coronal hole, at different distances from the Sun. Along with the adjacent slow solar wind streams, these intervals represent a unique database for studying the radial evolution of turbulence in samples of undisturbed solar wind. A set of direct numerical simulations of the MHD equations performed with the Lattice-Boltzmann code FLAME was also used for interpretation. Results. We show that the turbulence energy transfer rate decays approximately as a power law of the distance and that both the amplitude and decay law correspond to the observed radial temperature profile in the fast wind case. Results from MHD numerical simulations of decaying MHD turbulence show a similar trend for the total dissipation, suggesting an interpretation of the observed dynamics in terms of decaying turbulence and that multi-spacecraft studies of the solar wind radial evolution may help clarify the nature of the evolution of the turbulent fluctuations in the ecliptic solar wind.

MHD turbulence: a biased review

This review of scaling theories of magnetohydrodynamic (MHD) turbulence aims to put the developments of the last few years in the context of the canonical time line (from Kolmogorov to Iroshnikov–Kraichnan to Goldreich–Sridhar to Boldyrev). It is argued that Beresnyak's (valid) objection that Boldyrev's alignment theory, at least in its original form, violates the Reduced-MHD rescaling symmetry can be reconciled with alignment if the latter is understood as an intermittency effect. Boldyrev's scalings, a version of which is recovered in this interpretation, and the concept of dynamic alignment (equivalently, local 3D anisotropy) are thus an example of a physical theory of intermittency in a turbulent system. The emergence of aligned structures naturally brings into play reconnection physics and thus the theory of MHD turbulence becomes intertwined with the physics of tearing, current-sheet disruption and plasmoid formation. Recent work on these subjects by Loureiro, Malletet al.is reviewed and it is argued that we may, as a result, finally have a reasonably complete picture of the MHD turbulent cascade (forced, balanced, and in the presence of a strong mean field) all the way to the dissipation scale. This picture appears to reconcile Beresnyak's advocacy of the Kolmogorov scaling of the dissipation cutoff (as$\mathrm {Re}^{3/4}$) with Boldyrev's aligned cascade. It turns out also that these ideas open the door to some progress in understanding MHD turbulence without a mean field – MHD dynamo – whose saturated state is argued to be controlled by reconnection and to contain, at small scales, a tearing-mediated cascade similar to its strong-mean-field counterpart (this is a new result). On the margins of this core narrative, standard weak-MHD-turbulence theory is argued to require some adjustment – and a new scheme for such an adjustment is proposed – to take account of the determining part that a spontaneously emergent 2D condensate plays in mediating the Alfvén-wave cascade from a weakly interacting state to a strongly turbulent (critically balanced) one. This completes the picture of the MHD cascade at large scales. A number of outstanding issues are surveyed: imbalanced turbulence (for which a new, tentative theory is proposed), residual energy, MHD turbulence at subviscous scales, and decaying MHD turbulence (where there has been dramatic progress recently, and reconnection again turned out to feature prominently). Finally, it is argued that the natural direction of research is now away from the fluid MHD theory and into kinetic territory – and then, possibly, back again. The review lays no claim to objectivity or completeness, focusing on topics and views that the author finds most appealing at the present moment.

Reynolds-number dependence of the dimensionless dissipation rate in homogeneous magnetohydrodynamic turbulence.

This paper examines the behavior of the dimensionless dissipation rate C_{ɛ} for stationary and nonstationary magnetohydrodynamic (MHD) turbulence in the 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_{ɛ} 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_{ɛ,∞} lie between 0.193≤C_{ɛ,∞}≤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_{ɛ,∞}=0.223.

Nonuniversality and Finite Dissipation in Decaying Magnetohydrodynamic Turbulence.

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ϵ=Cϵ,∞+C/R-+O(1/R-(2)), where R- is a generalized Reynolds number. The constant Cϵ,∞ 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ϵ,∞ 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.

A method for spectral DNS of low Rm channel flows based on the least dissipative modes

We put forward a new type of spectral method for the direct numerical simulation of flows where anisotropy or very fine boundary layers are present. The main idea is to take advantage of the fact that such structures are dissipative and that their presence should reduce the number of degrees of freedom of the flow, when paradoxically, their fine resolution incurs extra computational cost in most current methods. The principle of this method is to use a functional basis with elements that already include these fine structures so as to avoid these extra costs. This leads us to develop an algorithm to implement a spectral method for arbitrary functional bases, and in particular, non-orthogonal ones. We construct a basic implementation of this algorithm to simulate magnetohydrodynamic (MHD) channel flows with an externally imposed, transverse magnetic field, where very thin boundary layers are known to develop along the channel walls. In this case, the sought functional basis can be built out of the eigenfunctions of the dissipation operator, which incorporate these boundary layers, and it turns out to be non-orthogonal. We validate this new scheme against numerical simulations of freely decaying MHD turbulence based on a finite volume code and it is found to provide accurate results. Its ability to fully resolve wall-bounded turbulence with a number of modes close to that required by the dynamics is demonstrated on a simple example. This opens the way to full-blown simulations of MHD turbulence under very high magnetic fields. Until now such simulations were too computationally expensive. In contrast to traditional methods the computational cost of the proposed method, does not depend on the intensity of the magnetic field.

Nonhelical inverse transfer of a decaying turbulent magnetic field.

In the presence of magnetic helicity, inverse transfer from small to large scales is well known in magnetohydrodynamic (MHD) turbulence and has applications in astrophysics, cosmology, and fusion plasmas. Using high resolution direct numerical simulations of magnetically dominated self-similarly decaying MHD turbulence, we report a similar inverse transfer even in the absence of magnetic helicity. We compute for the first time spectral energy transfer rates to show that this inverse transfer is about half as strong as with helicity, but in both cases the magnetic gain at large scales results from velocity at similar scales interacting with smaller-scale magnetic fields. This suggests that both inverse transfers are a consequence of universal mechanisms for magnetically dominated turbulence. Possible explanations include inverse cascading of the mean squared vector potential associated with local near two dimensionality and the shallower k^{2} subinertial range spectrum of kinetic energy forcing the magnetic field with a k^{4} subinertial range to attain larger-scale coherence. The inertial range shows a clear k^{-2} spectrum and is the first example of fully isotropic magnetically dominated MHD turbulence exhibiting weak turbulence scaling.

Quantifying energetics and dissipation in magnetohydrodynamic turbulence

We perform a suite of two- and three-dimensional magnetohydrodynamic (MHD) simulations with the Athena code of the non-driven Kelvin-Helmholtz instability in the subsonic, weak magnetic field limit. Focusing the analysis on the non-linear turbulent regime, we quantify energy transfer on a scale-by-scale basis and identify the physical mechanisms responsible for energy exchange by developing the diagnostic known as spectral energy transfer function analysis. At late times when the fluid is in a state of MHD turbulence, magnetic tension mediates the dominant mode of energy injection into the magnetic reservoir, whereby turbulent fluid motions twist and stretch the magnetic field lines. This generated magnetic energy turbulently cascades to smaller scales, while being exchanged backwards and forwards with the kinetic energy reservoir, until finally being dissipated. Incorporating explicit dissipation pushes the dissipation scale to larger scales than if the dissipation were entirely numerical. For scales larger than the dissipation scale, we show that the physics of energy transfer in decaying MHD turbulence is robust to numerical effects.

LONG-TERM EVOLUTION OF DECAYING MAGNETOHYDRODYNAMIC TURBULENCE IN THE MULTIPHASE INTERSTELLAR MEDIUM

Supersonic turbulence in the interstellar medium (ISM) is believed to decay rapidly within a flow crossing time irrespective of the degree of magnetization. However, this general consensus of decaying magnetohydrodynamic (MHD) turbulence relies on local isothermal simulations, which are unable to take into account the roles of global structures of magnetic fields and the ISM. Utilizing three-dimensional MHD simulations including interstellar cooling and heating, we investigate decaying MHD turbulence within cold neutral medium sheets embedded in warm neutral medium. Early evolution of turbulent kinetic energy is consistent with previous results for decaying compressible MHD turbulence characterized by rapid energy decay with a power-law form of $E\propto t^{-1}$ and by short decay time compared to a flow crossing time. If initial magnetic fields are strong and perpendicular to the sheet, however, long-term evolution of the kinetic energy shows that a significant amount of turbulent energy ($\sim 0.2E_0$) still remains even after ten flow crossing times for models with periodic boundary conditions. The decay rate is also greatly reduced as the field strength increases for such initial and boundary conditions, but not if the boundary conditions are that for a completely isolated sheet. We analyze velocity power spectra of the remaining turbulence to show that in-plane, incompressible motions parallel to the sheet dominate at later times.

Broken ergodicity in magnetohydrodynamic turbulence

Turbulent magnetofluids appear in various geophysical and astrophysical contexts, in phenomena associated with planets, stars, galaxies and the universe itself. In many cases, large-scale magnetic fields are observed, though a better knowledge of magnetofluid turbulence is needed to more fully understand the dynamo processes that produce them. One approach is to develop the statistical mechanics of ideal (i.e. non-dissipative), incompressible, homogeneous magnetohydrodynamic (MHD) turbulence, known as “absolute equilibrium ensemble” theory, as far as possible by studying model systems with the goal of finding those aspects that survive the introduction of viscosity and resistivity. Here, we review the progress that has been made in this direction. We examine both three-dimensional (3-D) and two-dimensional (2-D) model systems based on discrete Fourier representations. The basic equations are those of incompressible MHD and may include the effects of rotation and/or a mean magnetic field B o. Statistical predictions are that Fourier coefficients of the velocity and magnetic field are zero-mean random variables. However, this is not the case, in general, for we observe non-ergodic behavior in very long time computer simulations of ideal turbulence: low wavenumber Fourier modes that have relatively large means and small standard deviations, i.e. coherent structure. In particular, ergodicity appears strongly broken when B o = 0 and weakly broken when B o ≠ 0. Broken ergodicity in MHD turbulence is explained by an eigenanalysis of modal covariance matrices. This produces a set of modal eigenvalues inversely proportional to the expected energy of their associated eigenvariables. A large disparity in eigenvalues within the same mode (identified by wavevector k ) can occur at low values of wavenumber k = | k |, especially when B o = 0. This disparity breaks the ergodicity of eigenvariables with smallest eigenvalues (largest energies). This leads to coherent structure in models of ideal homogeneous MHD turbulence, which can occur at lowest values of wavenumber k for 3-D cases, and at either lowest or highest k for ideal 2-D magnetofluids. These ideal results appear relevant for unforced, decaying MHD turbulence, so that broken ergodicity effects in MHD turbulence survive dissipation. In comparison, we will also examine ideal hydrodynamic (HD) turbulence, which, in the 3-D case, will be seen to differ fundamentally from ideal MHD turbulence in that coherent structure due to broken ergodicity can only occur at maximum k in numerical simulations. However, a nonzero viscosity eliminates this ideal 3-D HD structure, so that unforced, decaying 3-D HD turbulence is expected to be ergodic. In summary, broken ergodicity in MHD turbulence leads to energetic, large-scale, quasistationary magnetic fields (coherent structures) in numerical models of bounded, turbulent magnetofluids. Thus, broken ergodicity provides a large-scale dynamo mechanism within computer models of homogeneous MHD turbulence. These results may help us to better understand the origin of global magnetic fields in astrophysical and geophysical objects.

Systematics of the magnetic-Prandtl-number dependence of homogeneous, isotropic magnetohydrodynamic turbulence

We present the results of our detailed pseudospectral direct numerical simulation (DNS) studies, with up to 10243 collocation points, of incompressible, magnetohydrodynamic (MHD) turbulence in three dimensions, without a mean magnetic field. Our study concentrates on the dependence of various statistical properties of both decaying and statistically steady MHD turbulence on the magnetic Prandtl number PrM over a large range, namely 0.01⩽PrM⩽10. We obtain data for a wide variety of statistical measures, such as probability distribution functions (PDFs) of the moduli of the vorticity and current density, the energy dissipation rates, and velocity and magnetic-field increments, energy and other spectra, velocity and magnetic-field structure functions, which we use to characterize intermittency, isosurfaces of quantities, such as the moduli of the vorticity and current density, and joint PDFs, such as those of fluid and magnetic dissipation rates. Our systematic study uncovers interesting results that have not been noted hitherto. In particular, we find a crossover from a larger intermittency in the magnetic field than in the velocity field, at large PrM, to a smaller intermittency in the magnetic field than in the velocity field, at low PrM. Furthermore, a comparison of our results for decaying MHD turbulence and its forced, statistically steady analogue suggests that we have strong universality in the sense that, for a fixed value of PrM, multiscaling exponent ratios agree, at least within our error bars, for both decaying and statistically steady homogeneous, isotropic MHD turbulence.

Broken symmetries and magnetic dynamos

Phase space symmetries inherent in the statistical theory of ideal magnetohydrodynamic (MHD) turbulence are known to be broken dynamically to produce large-scale coherent magnetic structure. Here, results of a numerical study of decaying MHD turbulence are presented that show large-scale coherent structure also arises and persists in the presence of dissipation. Dynamically broken symmetries in MHD turbulence may thus play a fundamental role in the dynamo process.

Large-eddy simulation of very large kinetic and magnetic Reynolds number isotropic magnetohydrodynamic turbulence using a spectral subgrid model

A spectral eddy viscosity and magnetic resistivity subgrid-scale model based on the eddy-damped quasi-normal Markovian (EDQNM) kinetic and magnetic energy transfers is used in large-eddy simulation (LES) of asymptotically large kinetic and magnetic Reynolds number magnetohydrodynamic (MHD) turbulence. The model is assessed a posteriori on three-dimensional, incompressible, isotropic, nonhelical, freely decaying MHD turbulence. Using LES initialized with spectra such that the Alfvén ratio of kinetic to magnetic energy equals unity, it is shown that the kinetic energy and magnetic energy spectra exhibit k−5∕3 Kolmogorov inertial subrange scalings consistent with the EDQNM model.

The residual energy in freely decaying magnetohydrodynamic turbulence

Within the framework of the eddy-damped quasi-normal Markovian approximation for incompressible isotropic magnetohydrodynamic (MHD) turbulence a prediction for the inertial range scaling of the residual energy spectrum, , is obtained. This scaling, while in contradiction to earlier theoretical results, is shown to be in agreement with high-resolution direct numerical simulations of nonhelical decaying MHD turbulence. The underlying phenomenology states a dynamic quasi-equilibrium of the small-scale turbulent dynamo and the Alfvén effect.

Large-eddy simulation of magnetohydrodynamic turbulence

Pseudospectral large-eddy simulations (LES) of incompressible magnetohydrodynamic (MHD) turbulence are performed with a maximum of 643 collocation points. Different gradient-diffusion type subgrid models are evaluated. This is done by comparing the obtained LES results to reference data stemming from high-resolution direct numerical simulations of decaying MHD turbulence with up to 5123 spectral modes. A newly developed MHD subgrid model, based on the ideal cross-helicity invariant of MHD turbulence, is tested and observed to perform excellently well.