We study the spatio-temporal behavior of the Els\"asser variables describing magnetic and velocity field fluctuations, using direct numerical simulations of three-dimensional magnetohydrodynamic turbulence. We consider cases with relatively small, intermediate, and large values of a mean background magnetic field, and with null, small, and high cross-helicity (correlations between the velocity and the magnetic field). Wavenumber-dependent time correlation functions are computed for the different simulations. From these correlation functions, the decorrelation time is computed and compared with different theoretical characteristic times: the local non-linear time, the random-sweeping time, and the Alfv\'enic time. It is found that decorrelation times are dominated by sweeping effects for low values of the mean magnetic field and for low values of the cross-helicity, while for large values of the background field or of the cross-helicity and for wave vectors sufficiently aligned with the guide field, decorrelation times are controlled by Alfv\'enic effects. Finally, we observe counter-propagation of Alfv\'enic fluctuations due to reflections produced by inhomogeneities in the total magnetic field. This effect becomes more prominent in flows with large cross-helicity, strongly modifying the propagation of waves in turbulent magnetohydrodynamic flows.

This paper briefly reviews a number of fundamental measurements that need to be made in order to characterize turbulence in space plasmas such as the solar wind. It has long been known that many of these quantities require simultaneous multipoint measurements to attain a proper characterization that would reveal the fundamental physics of plasma turbulence. The solar wind is an ideal plasma for such an investigation, and it now appears to be technologically feasible to carry out such an investigation, following the pioneering Cluster and MMS missions. Quantities that need to be measured using multipoint measurements include the two-point, two-time second correlation function of velocity, magnetic field and density, and higher order statistical objects such as third and fourth order structure functions. Some details of these requirements are given here, with a eye towards achieving closure on fundamental questions regarding the cascade rate, spectral anisotropy, characteristic coherent structures, intermittency, and dissipation mechanisms that describe plasma turbuelence, as well as its variability with plasma parameters in the solar wind. The motivation for this discussion is the current planning for a proposed Helioswarm mission that would be designed to make these measurements,leading to breakthrough understanding of the physics of space and astrophysical turbulence.

We live in an age in which high-performance computing is transforming the way we do science. Previously intractable problems are now becoming accessible by means of increasingly realistic numerical simulations. One of the most enduring and most challenging of these problems is turbulence. Yet, despite these advances, the extreme parameter regimes encountered in space physics and astrophysics (as in atmospheric and oceanic physics) still preclude direct numerical simulation. Numerical models must take a Large Eddy Simulation (LES) approach, explicitly computing only a fraction of the active dynamical scales. The success of such an approach hinges on how well the model can represent the subgrid-scales (SGS) that are not explicitly resolved. In addition to the parameter regime, heliophysical and astrophysical applications must also face an equally daunting challenge: magnetism. The presence of magnetic fields in a turbulent, electrically conducting fluid flow can dramatically alter the coupling between large and small scales, with potentially profound implications for LES/SGS modeling. In this review article, we summarize the state of the art in LES modeling of turbulent magnetohydrodynamic (MHD) flows. After discussing the nature of MHD turbulence and the small-scale processes that give rise to energy dissipation, plasma heating, and magnetic reconnection, we consider how these processes may best be captured within an LES/SGS framework. We then consider several specific applications in heliophysics and astrophysics, assessing triumphs, challenges, and future directions.

Employing a simple ideal magnetohydrodynamic model in spherical geometry,we show that the presence of either rotation or finite magnetic helicity is sufficient to induce dynamical reversals of the magnetic dipole moment. The statistical character of the model is similar to that of terrestrial magnetic field reversals, with the similarity being stronger when rotation is present.The connection between long time correlations, $1/f$ noise, and statistics of reversals is supported, consistent with earlier suggestions.

A review of some of the issues that have arisen over the years concerning the energy distribution among scales for magnetohydrodynamics (MHD) turbulence is given here. A variety of tools are employed to that effect, and a central role is played by taking into consideration the ideal (non-dissipative) invariants, namely the total energy, the magnetic helicity and the cross-correlations between the velocity and the magnetic field (concentrating on the three-dimensional case). These concepts, based mostly on theory, models and direct numerical simulations, are briefly put in the context of observations, in particular the solar wind, and some of the remaining open questions are delineated as well. New results on ideal MHD dynamics in three dimensions on equivalent grids of up to $6144^3$ points using the Taylor-Green flow generalized to MHD are also mentioned.

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 backreaction from any non-helical 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.

This study shows the connection between three previously observed but seemingly unrelated phenomena in hydrodynamic (HD) and magnetohydrodynamic (MHD) turbulent flows, involving the emergence of fluctuations occurring on very long time scales: the low-frequency 1/f noise in the power frequency spectrum, the delayed ergodicity of complex valued amplitude fluctuations in wavenumber space, and the spontaneous flippings or reversals of large scale fields. Direct numerical simulations of ideal MHD and HD are employed in three space dimensions, at low resolution, for long periods of time, and with high accuracy to study several cases: Different geometries, presence of rotation and/or a uniform magnetic field, and different values of the associated conserved global quantities. It is conjectured that the origin of all these long-time phenomena is rooted in the interaction of the longest wavelength fluctuations available to the system with fluctuations at much smaller scales. The strength of this non-local interaction is controlled either by the existence of conserved global quantities with a back-transfer in Fourier space, or by the presence of a slow manifold in the dynamics.

We demonstrate that, for the case of quasi-equipartition between the velocity and the magnetic field, the Lagrangian-averaged magnetohydrodynamics alpha-model (LAMHD) reproduces well both the large-scale and small-scale properties of turbulent flows; in particular, it displays no increased (super-filter) bottleneck effect with its ensuing enhanced energy spectrum at the onset of the sub-filter-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 super-filter-scale spectral properties. No such regions are found in LAMHD, making this method capable of large reductions in required numerical degrees of freedom; specifically, we find a reduction factor of 200 when compared to a direct numerical simulation on a large grid of 1536^3 points at the same Reynolds number.

We present an analysis of data stemming from numerical simulations of decaying magnetohydrodynamic (MHD) turbulence up to grid resolution of 1536^3 points and up to Taylor Reynolds number of 1200. The initial conditions are such that the initial velocity and magnetic fields are helical and in equipartition, while their correlation is negligible. Analyzing the data at the peak of dissipation, we show that the dissipation in MHD seems to asymptote to a constant as the Reynolds number increases, thereby strengthening the possibility of fast reconnection events in the solar environment for very large Reynolds numbers. Furthermore, intermittency of MHD flows, as determined by the spectrum of anomalous exponents of structure functions of the velocity and the magnetic field, is stronger than for fluids, confirming earlier results; however, we also find that there is a measurable difference between the exponents of the velocity and those of the magnetic field, as observed recently in the solar wind. Finally, we discuss the spectral scaling laws that arise in this flow.

A recently developed spectral-element adaptive refinement incompressible magnetohydrodynamic (MHD) code [Rosenberg, Fournier, Fischer, Pouquet, J. Comp. Phys. 215, 59-80 (2006)] is applied to simulate the problem of MHD island coalescence instability (MICI) in two dimensions. MICI is a fundamental MHD process that can produce sharp current layers and subsequent reconnection and heating in a high-Lundquist number plasma such as the solar corona [Ng and Bhattacharjee, Phys. Plasmas, 5, 4028 (1998)]. Due to the formation of thin current layers, it is highly desirable to use adaptively or statically refined grids to resolve them, and to maintain accuracy at the same time. The output of the spectral-element static adaptive refinement simulations are compared with simulations using a finite difference method on the same refinement grids, and both methods are compared to pseudo-spectral simulations with uniform grids as baselines. It is shown that with the statically refined grids roughly scaling linearly with effective resolution, spectral element runs can maintain accuracy significantly higher than that of the finite difference runs, in some cases achieving close to full spectral accuracy.

We show that local directional alignment of the velocity and magnetic field fluctuations occurs rapidly in magnetohydrodynamics for a variety of parameters. This is observed both in direct numerical simulations and in solar wind data. The phenomenon is due to an alignment between the magnetic field and either pressure gradients or shear-associated kinetic energy gradients. A similar alignment, of velocity and vorticity, occurs in the Navier Stokes fluid case. This may be the most rapid and robust relaxation process in turbulent flows, and leads to a local weakening of the nonlinear terms in the small scale vorticity and current structures where alignment takes place.

We present a numerical analysis of an incompressible decaying magnetohydrodynamic turbulence run on a grid of 1536^3 points. The Taylor Reynolds number at the maximum of dissipation is ~1100, and the initial condition is a superposition of large scale ABC flows and random noise at small scales, with no uniform magnetic field. The initial kinetic and magnetic energies are equal, with negligible correlation. The resulting energy spectrum is a combination of two components, each moderately resolved. Isotropy obtains in the large scales, with a spectral law compatible with the Iroshnikov-Kraichnan theory stemming from the weakening of nonlinear interactions due to Alfven waves; scaling of structure functions confirms the non-Kolmogorovian nature of the flow in this range. At small scales, weak turbulence emerges with a k_{\perp}^{-2} spectrum, the perpendicular direction referring to the local quasi-uniform magnetic field.

We present an extension of the Karman-Howarth theorem to the Lagrangian averaged magnetohydrodynamic (LAMHD-alpha) equations. The scaling laws resulting as a corollary of this theorem are studied in numerical simulations, as well as the scaling of the longitudinal structure function exponents indicative of intermittency. Numerical simulations for a magnetic Prandtl number equal to unity are presented both for freely decaying and for forced two dimensional MHD turbulence, solving directly the MHD equations, and employing the LAMHD-alpha equations at 1/2 and 1/4 resolution. Linear scaling of the third-order structure function with length is observed. The LAMHD-alpha equations also capture the anomalous scaling of the longitudinal structure function exponents up to order 8.

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.

A three-dimensional numerical computation of magnetohydrodynamic dynamo behavior is described. The dynamo is mechanically forced with a driving term of the Taylor-Green type. The magnetic field development is followed from negligibly small levels to saturated values that occur at magnetic energies comparable to the kinetic energies. Though there is locally a helicity density, there is no overall integrated helicity in the system. Persistent oscillations are observed in the saturated state for not-too-large mechanical Reynolds numbers, oscillations in which the kinetic and magnetic energies vary out of phase but with no reversal of the magnetic field. The flow pattern exhibits considerable geometrical structure in this regime. As the Reynolds number is raised, the oscillations disappear and the energies become more nearly stationary, but retain some unsystematically fluctuating turbulent time dependence. The regular geometrical structure of the fields gives way to a more spatially disordered distribution. The injection and dissipation scales are identified and the different components of energy transfer in Fourier space are analyzed, in particular in the context of clarifying the role played by different flow scales in the amplification of the magnetic field.

The observations of magnetic field variations as a signature of flaring activity is one of the main goal in solar physics. Some efforts in the past give apparently no unambiguous observations of changes. We observed that the scaling laws of the current helicity inside a given flaring active region change clearly and abruptly in correspondence with the eruption of big flares at the top of that active region. Comparison with numerical simulations of MHD equations, indicates that the change of scaling behavior in the current helicity, seems to be associated to a topological reorganization of the footpoint of the magnetic field loop, namely to dissipation of small scales structures in turbulence. It is evident that the possibility of forecasting in real time high energy flares, even if partially, has a wide practical interest to prevent the effects of big flares on Earth and its environment.

We derive a weak turbulence formalism for incompressible MHD. Three-wave interactions lead to a system of kinetic equations for the spectral densities of energy and helicity. We find energy spectra solution of the kinetic equations. The constants of the spectra are computed exactly and found to depend on the amount of correlation between the velocity and the magnetic field. Comparison with several numerical simulations and models is also made.

We investigate the properties of highly compressible turbulence and its ability to produce self-gravitating structures. The compressibility is parameterized by an effective polytropic exponent gama-eff. In the limit of small gama-eff, the density jump at shocks is shown to be of the order of e^{M^2}, and the production of vorticity by the nonlinear terms appears to be negligible. In the presence of self-gravity, we suggest that turbulence can produce bound structures for gama-eff < 2(1-1/n), where 'n' is the typical dimensionality of the turbulent compressions. We show, by means of numerical simulations, that, for sufficiently small gama-eff, small-scale turbulent density fluctuations eventually collapse even though the medium is globally stable. This result is preserved in the presence of a magnetic field for supercritical mass-to-flux ratios.

We investigate the properties of highly compressible turbulence, the compressibility arising from a small effective polytropic exponent $\gamma_e$ due to cooling. In the limit of small $\gamma_e$, the density jump at shocks is shown to be of the order of $e^{M^2}$. Without self-gravity, the density structures arising in the moderately compressible case consist mostly of patches separated by shocks and behaving like waves, while in the highly compressible case clearly defined long-lived object-like clouds emerge. When the forcing in the momentum equation is purely compressible, the rotational energy decays monotonically in time, indicating that the vortex-stretching term is not efficient in transferring energy to rotational modes. This property may be at the origin of the low amount of rotation found in interstellar clouds. Vorticity production is found to rely heavily on the presence of additional terms in the equations. In the presence of self-gravity, we suggest that turbulence can produce bound structures for $\gamma_e < 2(1-1/n)$, where $n$ is the typical dimensionality of the turbulent compressions. We support this result by means of numerical simulations in which, for sufficiently small $\gamma_e$, small-scale turbulent density fluctuations eventually collapse even though the medium is globally stable. This result is preserved in the presence of a magnetic field for supercritical mass-to-flux ratios. At larger polytropic exponents, turbulence alone is not capable of producing bound structures, and collapse can only occur when the medium is globally unstable. This mechanism is a plausible candidate for the differentiation between primordial and present-day stellar-cluster formation and for the low efficiency of star formation.

We discuss the role of turbulence in cloud and star formation, as observed in numerical simulations of the interstellar medium. Turbulent compression at the interfaces of colliding gas streams is responsible for the formation of intermediate ($\simlt 100$ pc) and small clouds (a few tens of pc), although the smallest clouds can also form from fragmentation of expanding shells around stellar heating centers. The largest cloud complexes (several hundred pc) seem to form by slow, gravitational instability-driven merging of individual clouds, which can actually be described as a large-scale tendency towards homogenization of the flow due to gravity rather than cloud collisions. These mechanisms operate as well in the presence of a magnetic field and rotation, although slight variations on the compressibility and cloud morphology are present which depend on the strength and topology of the field. In summary, the role of turbulence in the life-cycle of clouds appears to be twofold: small-scale modes contribute to cloud support, while large-scale modes can both form or disrupt clouds.

We present results from two-dimensional numerical simulations of a supersonic turbulent flow in the plane of the galactic disk, incorporating shear, thresholded and discrete star formation (SF), self-gravity, rotation and magnetic fields. A test of the model in the linear regime supports the results of the linear theory of Elmegreen (1991). In the fully nonlinear turbulent regime, while some results of the linear theory persist, new effects also emerge. Some exclusively nonlinear effects are: a) Even though there is no dynamo in 2D, the simulations are able to maintain or increase their net magnetic energy in the presence of a seed uniform azimuthal component. b) A well-defined power-law magnetic spectrum and an inverse magnetic cascade are observed in the simulations, indicating full MHD turbulence. Thus, magnetic field energy is generated in regions of SF and cascades up to the largest scales. c) The field has a slight but noticeable tendency to be aligned with density features. d) The magnetic field prevents HII regions from expanding freely, as in the recent results of Slavin \& Cox (1993). e) A tendency to exhibit {\it less} filamentary structures at stronger values of the uniform component of the magnetic field is present in several magnetic runs. f) For fiducial values of the parameters, the flow in general appears to be in rough equipartition between magnetic and kinetic energy. There is no clear domination of either the magnetic or the inertial forces. g) A median value of the magnetic field strength within clouds is $\sim 12\mu$G, while for the intercloud medium a value of $\sim 3\mu$G is found. Maximum contrasts of up to a factor of $\sim 10$ are observed.