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.

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.

The investigation of the diffusive transport of charged particles in a turbulent magnetic field remains a subject of considerable interest. Research has most frequently concentrated on determining the diffusion coefficient in the presence of a mean magnetic field. Here we consider diffusion of charged particles in fully three-dimensional isotropic turbulent magnetic fields with no mean field, which may be pertinent to many astrophysical situations. We identify different ranges of particle energy depending upon the ratio of the Larmor radius of the charged particle to the characteristic outer length scale of the turbulence. Two different theoretical models are proposed to calculate the diffusion coefficient, each applicable to a distinct range of particle energies. The theoretical results are compared with those from computer simulations, showing good agreement.

We present a full two-fluid magnetohydrodynamic (MHD) description for a completely ionized hydrogen plasma, retaining the effects of the Hall current, electron pressure and electron inertia. According to this description, each plasma species introduces a new spatial scale: the ion inertial length $\lambda_{i}$ and the electron inertial length $\lambda_{e}$, which are not present in the traditional MHD description. In the present paper, we seek for possible changes in the energy power spectrum in fully developed turbulent regimes, using numerical simulations of the two-fluid equations in two-and-a-half dimensions (2.5D). We have been able to reproduce different scaling laws in different spectral ranges, as it has been observed in the solar wind for the magnetic energy spectrum. At the smallest wavenumbers where plain MHD is valid, we obtain an inertial range following a Kolmogorov $k^{-5/3}$ law. For intermediate wavenumbers such that $\lambda_{i}^{-1} << k << \lambda_{e}^{-1}$, the spectrum is modified to a $k^{-7/3}$ power-law, as has also been obtained for Hall-MHD (HMHD) neglecting electron inertia terms. When electron inertia is retained, a new spectral region given by $k > \lambda_{e}^{-1}$ arises. The power spectrum for magnetic energy in this region is given by a $k^{-11/3}$ power law. Finally, when the terms of electron inertia are retained, we study the self-consistent electric field. Our results are discussed and compared with those obtained in solar wind observations and previous simulations.

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.

The acceleration of charged particles is relevant to the solar corona over a broad range of scales and energies. High-energy particles are usually detected in concomitance with large energy release events like solar eruptions and flares, nevertheless acceleration can occur at smaller scales, characterized by dynamical activity near current sheets. To gain insight into the complex scenario of coronal charged particle acceleration, we investigate the properties of acceleration with a test-particle approach using three-dimensional magnetohydrodynamic (MHD) models. These are obtained from direct solutions of the reduced MHD equations, well suited for a plasma embedded in a strong axial magnetic field, relevant to the inner heliosphere. A multi-box, multi-scale technique is used to solve the equations of motion for protons. This method allows us to resolve an extended range of scales present in the system, namely from the ion inertial scale of the order of a meter up to macroscopic scales of the order of $10\,$km ($1/100$th of the outer scale of the system). This new technique is useful to identify the mechanisms that, acting at different scales, are responsible for acceleration to high energies of a small fraction of the particles in the coronal plasma. We report results that describe acceleration at different stages over a broad range of time, length and energy scales.

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.

Correlation anisotropy emerges dynamically in magnetohydrodynamics (MHD), producing stronger gradients across the large-scale mean magnetic field than along it. This occurs both globally and locally, and has significant implications in space and astrophysical plasmas, including particle scattering and transport, and theories of turbulence. Properties of local correlation anisotropy are further documented here by showing through numerical experiments that the effect is intensified in more localized estimates of the mean field. The mathematical formulation of this property shows that local anisotropy mixes second-order with higher order correlations. Sensitivity of local statistical estimates to higher order correlations can be understood in connection with the stochastic coordinate system inherent in such formulations. We demonstrate this in specific cases, illustrate the connection to higher order statistics by showing the sensitivity of local anisotropy to phase randomization, and thus establish that the local structure function is not a measure of the energy spectrum. Evidently the local enhancement of correlation anisotropy is of substantial fundamental interest, and this phenomenon must be understood in terms of higher order correlations, fourth-order and above.

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.

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.

A model is presented for generation of fast solar wind in coronal holes, relying on heating that is dominated by turbulent dissipation of MHD fluctuations transported upwards in the solar atmosphere. Scale-separated transport equations include large-scale fields, transverse Alfvenic fluctuations, and a small compressive dissipation due to parallel shears near the transition region. The model accounts for proton temperature, density, wind speed, and fluctuation amplitude as observed in remote sensing and in situ satellite data.

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.

The heating of the lower solar corona is examined using numerical simulations and theoretical models of magnetohydrodynamic turbulence in open magnetic regions. A turbulent energy cascade to small length scales perpendicular to the mean magnetic field can be sustained by driving with low-frequency Alfven waves reflected from mean density and magnetic field gradients. This mechanism deposits energy efficiently in the lower corona, and we show that the spatial distribution of the heating is determined by the mean density through the Alfven speed profile. This provides a robust heating mechanism that can explain observed high coronal temperatures and accounts for the significant heating (per unit volume) distribution below two solar radius needed in models of the origin of the solar wind. The obtained heating per unit mass on the other hand is much more extended indicating that the heating on a per particle basis persists throughout all the lower coronal region considered here.

We report preliminary results from a series of numerical simulations of the reduced magnetohydrodynamic equations, used to describe the dynamics of magnetic loops in active regions of the solar corona. A stationary velocity field is applied at the photospheric boundaries to imitate the driving action of granule motions. A turbulent stationary regime is reached, characterized by a broadband power spectrum $E_k\simeq k^{-3/2}$ and heating rate levels compatible with the heating requirements of active region loops. A dimensional analysis of the equations indicates that their solutions are determined by two dimensionless parameters: the Reynolds number and the ratio between the Alfven time and the photospheric turnover time. From a series of simulations for different values of this ratio, we determine how the heating rate scales with the physical parameters of the problem, which might be useful for an observational test of this model.

We perform direct numerical simulations of an externally driven two-dimensional magnetohydrodynamic system over extended periods of time to simulate the dynamics of a transverse section of a solar coronal loop. A stationary and large-scale magnetic forcing was imposed, to model the photospheric motions at the magnetic loop footpoints. A turbulent stationary regime is reached, which corresponds to energy dissipation rates consistent with the heating requirements of coronal loops. The temporal behavior of quantities such as the energy dissipation rate show clear indications of intermittency, which are exclusively due to the strong nonlinearity of the system. We tentatively associate these impulsive events of magnetic energy dissipation to the so-called nanoflares. A statistical analysis of these events yields a power law distribution as a function of their energies with a negative slope of 1.5, which is consistent with those obtained for flare energy distributions reported from X-ray observations.