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.

Using spatio-temporal spectra we show direct evidence of excitation of magnetosonic and Alfv\'en waves in three-dimensional compressible magnetohydrodynamic turbulence at small Mach numbers. For the plasma pressure dominated regime, or high $\beta$ regime (with $\beta$ the ratio between fluid and magnetic pressure), and for the magnetic pressure dominated regime, or low $\beta$ regime, we study magnetic field fluctuations parallel and perpendicular to a guide magnetic field $\textbf{B}_0$. In the low $\beta$ case we find excitation of compressible and incompressible fluctuations, with a transfer of energy towards Alfv\'enic modes and to a lesser extent towards magnetosonic modes. In particular, we find signatures of the presence of fast magnetosonic waves in a scenario compatible with that of weak turbulence. In the high $\beta$ case, fast and slow magnetosonic waves are present, with no clear trace of Alfv\'en waves, and a significant part of the energy is carried by two-dimensional turbulent eddies.

Three-dimensional direct numerical simulations are used to study the energy cascade rate in isothermal compressible magnetohydrodynamic turbulence. Our analysis is guided by a two-point exact law derived recently for this problem in which flux, source, hybrid, and mixed terms are present. The relative importance of each term is studied for different initial subsonic Mach numbers $M_S$ and different magnetic guide fields ${\bf B}_0$. The dominant contribution to the energy cascade rate comes from the compressible flux, which depends weakly on the magnetic guide field ${\bf B}_0$, unlike the other terms whose modulus increase significantly with $M_S$ and ${\bf B}_0$. In particular, for strong ${\bf B}_0$ the source and hybrid terms are dominant at small scales with almost the same amplitude but with a different sign. A statistical analysis made with an isotropic decomposition based on the SO(3) rotation group is shown to generate spurious results in presence of ${\bf B}_0$, when compared with an axisymmetric decomposition better suited to the geometry of the problem. Our numerical results are compared with previous analyses made with in-situ measurements in the solar wind and the terrestrial magnetosheath.

Recent high resolution AIA/SDO images show evidence of the development of the Kelvin-Helmholtz instability, as coronal mass ejections (CMEs) expand in the ambient corona. A large-scale magnetic field mostly tangential to the interface is inferred, both on the CME and on the background sides. However, the magnetic field component along the shear flow is not strong enough to quench the instability. There is also observational evidence that the ambient corona is in a turbulent regime, and therefore the criteria for the development of the instability are a-priori expected to differ from the laminar case. To study the evolution of the Kelvin-Helmholtz instability with a turbulent background, we perform three-dimensional simulations of the incompressible magnetohydrodynamic equations. The instability is driven by a velocity profile tangential to the CME-corona interface, which we simulate through a hyperbolic tangent profile. The turbulent background is generated by the application of a stationary stirring force. We compute the instability growth-rate for different values of the turbulence intensity, and find that the role of turbulence is to attenuate the growth. The fact that the Kelvin-Helmholtz instability is observed, sets an upper limit to the correlation length of the coronal background turbulence.

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.

We study the stability of shear flows in a fully ionized plasma. Kelvin-Helmholtz is a well known, macroscopic and ideal shear-driven instability. In sufficiently low density plasmas, also the microscopic Hall magneto-shear instability can take place. We performed three-dimensional simulations of the Hall-MHD equations where these two instabilities are present, and carried out a comparative study. We find that when the shear flow is so intense that its vorticity surpasses the ion-cyclotron frequency of the plasma, the Hall magneto-shear instability is not only non-negligible, but it actually displays growth rates larger than those of the Kelvin-Helmholtz instability.

We use direct numerical simulations to compute turbulent transport coefficients for passive scalars in turbulent rotating flows. Effective diffusion coefficients in the directions parallel and perpendicular to the rotations axis are obtained by studying the diffusion of an imposed initial profile for the passive scalar, and calculated by measuring the scalar average concentration and average spatial flux as a function of time. The Rossby and Schmidt numbers are varied to quantify their effect on the effective diffusion. It is find that rotation reduces scalar diffusivity in the perpendicular direction. The perpendicular diffusion can be estimated from mixing length arguments using the characteristic velocities and lengths perpendicular to the rotation axis. Deviations are observed for small Schmidt numbers, for which turbulent transport decreases and molecular diffusion becomes more significant.

We present a detailed study of intermittency in the velocity and magnetic field fluctuations of compressible Hall-magnetohydrodynamic turbulence with an external guide field. To solve the equations numerically, a reduced model valid when a strong guide field is present is used. Different values for the ion skin depth are considered in the simulations. The resulting data is analyzed computing field increments in several directions perpendicular to the guide field, and building structure functions and probability density functions. In the magnetohydrodynamic limit we recover the usual results with the magnetic field being more intermittent than the velocity field. In the presence of the Hall effect, field fluctuations at scales smaller than the ion skin depth show a substantial decrease in the level of intermittency, with close to monofractal scaling.

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.

This article reviews recent studies of scale interactions in magnetohydrodynamic turbulence. The present day increase of computing power, which allows for the exploration of different configurations of turbulence in conducting flows, and the development of shell-to-shell transfer functions, has led to detailed studies of interactions between the velocity and the magnetic field and between scales. In particular, processes such as induction and dynamo action, the damping of velocity fluctuations by the Lorentz force, or the development of anisotropies, can be characterized at different scales. In this context we consider three different configurations often studied in the literature: mechanically forced turbulence, freely decaying turbulence, and turbulence in the presence of a uniform magnetic field. Each configuration is of interest for different geophysical and astrophysical applications. Local and non-local transfers are discussed for each case. While the transfer between scales of solely kinetic or solely magnetic energy is local, transfers between kinetic and magnetic fields are observed to be local or non-local depending on the configuration. Scale interactions in the cascade of magnetic helicity are also reviewed. Based on the results, the validity of several usual assumptions in hydrodynamic turbulence, such as isotropy of the small scales or universality, is discussed.

Much of the progress in our understanding of dynamo mechanisms has been made within the theoretical framework of magnetohydrodynamics (MHD). However, for sufficiently diffuse media, the Hall effect eventually becomes non-negligible. We present results from three dimensional simulations of the Hall-MHD equations subjected to random non-helical forcing. We study the role of the Hall effect in the dynamo efficiency for different values of the Hall parameter, using a pseudospectral code to achieve exponentially fast convergence. We also study energy transfer rates among spatial scales to determine the relative importance of the various nonlinear effects in the dynamo process and in the energy cascade. The Hall effect produces a reduction of the direct energy cascade at scales larger than the Hall scale, and therefore leads to smaller energy dissipation rates. Finally, we present results stemming from simulations at large magnetic Prandtl numbers, which is the relevant regime in hot and diffuse media such a the interstellar medium.

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.

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.

Planetary and stellar dynamos likely result from turbulent motions in magnetofluids with kinematic viscosities that are small compared to their magnetic diffusivities. Laboratory experiments are in progress to produce similar dynamos in liquid metals. This work reviews recent computations of thresholds in critical magnetic Reynolds number above which dynamo amplification can be expected for mechanically-forced turbulence (helical and non-helical, short wavelength and long wavelength) as a function of the magnetic Prandtl number $P_M$. New results for helical forcing are discussed, for which a dynamo is obtained at $P_M=5\times10^{-3}$. The fact that the kinetic turbulent spectrum is much broader in wavenumber space than the magnetic spectrum leads to numerical difficulties which are bridged by a combination of overlapping direct numerical simulations and subgrid models of magnetohydrodynamic turbulence. Typically, the critical magnetic Reynolds number increases steeply as the magnetic Prandtl number decreases, and then reaches an asymptotic plateau at values of at most a few hundred. In the turbulent regime and for magnetic Reynolds numbers large enough, both small and large scale magnetic fields are excited. The interactions between different scales in the flow are also discussed.

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.

Direct numerical simulations of turbulent Hall dynamos are presented. The evolution of an initially weak and small scale magnetic field in a system maintained in a stationary turbulent regime by a stirring force at a macroscopic scale is studied to explore the conditions for exponential growth of the magnetic energy. Scaling of the dynamo efficiency with the Reynolds numbers is studied, and the resulting total energy spectra are found to be compatible with a Kolmogorov type law. A faster growth of large scale magnetic fields is observed at intermediate intensities of the Hall effect.

Dynamo activity caused by waves in a rotating magneto-plasma is investigated. In astrophysical environments such as accretion disks and at sufficiently small spatial scales, the Hall effect is likely to play an important role. It is shown that a combination of the Coriolis force and Hall effect can produce a finite $\alpha$-effect by generating net helicity in the small scales. The shear/ion-cyclotron normal mode of the Hall plasma is the dominant contributor to the dynamo action for short scale motions.

We present 3D, gasdynamic simulations of jet/cloud collisions, with the purpose of modelling the HH 270/110 system. From the models, we obtain predictions of H$\alpha$ and H$_2$ 1-0 s(1) emission line maps, which qualitatively reproduce some of the main features of the corresponding observations of HH 110. We find that the model that better reproduces the observed structures corresponds to a jet that was deflected at the surface of the cloud $\sim 1000$ yr ago, but is now boring a tunnel directly into the cloud. This model removes the apparent contradiction between the jet/cloud collision model and the lack of detection of molecular emission in the crossing region of the HH 270 and HH 110 axes.