Search

The self-consistency of the reduced magnetohydrodynamics (RMHD) model is explored by examining whether (parallel) spectral transfer might invalidate the assumptions employed in deriving it. Using direct numerical simulations we find that transfer of energy to structures with high parallel wavenumber is in fact limited by ongoing perpendicular transfer. Thus, the dynamics associated with RMHD models remains consistent with the underlying assumptions of RMHD. In particular, in well-resolved simulations it is neither necessary nor correct to introduce additional dissipation terms that (artificially) damp spectral transfer parallel to the mean magnetic field B0.

In a number of space and astrophysical plasmas, turbulence is driven by the supply of wave energy. In the context of incompressible magnetohydrodynamics (MHD) there are basic physical reasons, associated with conservation of cross helicity, why this kind of driving may be ineffective in sustaining turbulence. Here an investigation is made into some basic requirements for sustaining steady turbulence and dissipation in the context of incompressible MHD in a weakly inhomogeneous open field line region, driven by the supply of unidirectionally propagating waves at a boundary. While such wave driving cannot alone sustain turbulence, the addition of reflection permits sustainment. Another sustainment issue is the action of the nonpropagating or quasi-two dimensional part of the spectrum; this is particularly important in setting up a steady cascade. Thus, details of the wave boundary conditions also affect the ease of sustaining a cascade. Supply of a broadband spectrum of waves can overcome the latter difficu...

Space plasma measurements, laboratory experiments, and simulations have shown that magnetohydrodynamic (MHD) turbulence exhibits a dynamical tendency towards spectral anisotropy given a sufficiently strong background magnetic field. Here the undriven decaying initial-value problem for homogeneous MHD turbulence is examined with the purpose of characterizing the variation of spectral anisotropy of the turbulent fluctuations with magnetic field strength. Numerical results for both incompressible and compressible MHD are presented. A simple model for the scaling of this spectral anisotropy as a function of the fluctuating magnetic field over total magnetic field is offered. The arguments are based on ideas from reduced MHD (RMHD) dynamics and resonant driving of certain non-RMHD modes. The results suggest physical bases for explaining variations of the anisotropy with compressibility, Reynolds numbers, and spectral width of the (isotropic) initial conditions.

We evaluate a number of simple, one‐point phenomenological models for the decay of energy‐containing eddies in magnetohydrodynamic (MHD) and hydrodynamic turbulence. The MHD models include effects of cross helicity and Alfvenic couplings associated with a constant mean magnetic field, based on physical effects well‐described in the literature. The analytic structure of three separate MHD models is discussed. The single hydrodynamic model and several MHD models are compared against results from spectral‐method simulations. The hydrodynamic model phenomenology has been previously verified against experiments in wind tunnels, and certain experimentally determined parameters in the model are satisfactorily reproduced by the present simulation. This agreement supports the suitability of our numerical calculations for examining MHD turbulence, where practical difficulties make it more difficult to study physical examples. When the triple‐decorrelation time and effects of spectral anisotropy are properly taken i...

Long‐time states of a turbulent, decaying, two‐dimensional, Navier–Stokes flow are shown numerically to relax toward maximum‐entropy configurations, as defined by the ‘‘sinh‐Poisson’’ equation. The large‐scale Reynolds number is about 14 000, the spatial resolution is (512)2, the boundary conditions are spatially periodic, and the evolution takes place over nearly 400 large‐scale eddy‐turnover times.

We study the time development of the population of X-type critical points in a two-dimensional magnetohydrodynamic model during the early stages of freely decaying turbulence. At sufficiently high magnetic Reynolds number Rem, we find that the number of neutral points increases as Rem3/2, while the rates of reconnection at the most active sites decrease. The distribution of rates remains approximately exponential. We focus in particular on delicate issues of accuracy, which arise in these numerical experiments, in that the proliferation of X-points is also a feature of under-resolved simulations. The “splitting” of neutral points at high Reynolds number appears to be a fundamental feature of the cascade that has important implications for understanding the relationship between reconnection and turbulence, an issue of considerable importance for the Magnetospheric Multiscale and Solar Probe missions as well as observation of reconnection in the solar wind.

Magnetohydrodynamic (MHD) systems can be strongly nonlinear (turbulent) when their kinetic and magnetic Reynolds numbers are high, as is the case in many astrophysical and space plasma flows. Unfortunately these high Reynolds numbers are typically much greater than those currently attainable in numerical simulations of MHD turbulence. A natural question to ask is how can researchers be sure that their simulations have reproduced all of the most influential physics of the flows and magnetic fields? In this paper, a metric is defined to indicate whether the necessary physics of interest has been captured. It is found that current computing resources will typically not be sufficient to achieve this minimum state metric.

The widely recognized issue of adequate spatial resolution in numerical simulations of turbulence is studied in the context of two-dimensional magnetohydrodynamics. The familiar criterion that the dissipation scale should be resolved enables accurate computation of the spectrum, but fails for precise determination of higher-order statistical quantities. Examination of two straightforward diagnostics, the maximum of the kurtosis and the scale-dependent kurtosis, enables the development of simple tests for assessing adequacy of spatial resolution. The efficacy of the tests is confirmed by examining a sample problem, the distribution of magnetic reconnection rates in turbulence. Oversampling the Kolmogorov dissipation scale by a factor of 3 allows accurate computation of the kurtosis, the scale-dependent kurtosis, and the reconnection rates. These tests may provide useful guidance for resolution requirements in many plasma computations involving turbulence and reconnection.

The scaling laws of third-order structure functions for isotropic, homogeneous, and incompressible magnetohydrodynamic (MHD) turbulence relate the observable structure function with the energy dissipation rate. Recently [Wan et al. Phys. Plasmas 16, 090703 (2009)], the theory was extended to the case in which a constant velocity shear is present, motivated by the application of the third-order law to the solar wind. We use direct numerical simulations of two-dimensional MHD with shear to confirm this new generalization of the theory. The presence of the shear effect broadens the circumstances in which the law can be applied. Important implications for laboratory and space plasmas are discussed.

We extend the theory for third-order structure functions in homogeneous incompressible magnetohydrodynamic (MHD) turbulence to the case in which a constant velocity shear is present. A generalization is found of the usual relation [Politano and Pouquet, Phys. Rev. E, 57 21 (1998)] between third-order structure functions and the dissipation rate in steady inertial range turbulence, in which the shear plays a crucial role. In particular, the presence of shear leads to a third-order law which is not simply proportional to the relative separation. Possible implications for laboratory and space plasmas are discussed.

Spectral method simulations of ideal magnetohydrodynamics are used to investigate production of coherent small scale structures, a feature of fluid models that is usually associated with inertial range signatures of nonuniform dissipation, and the associated emergence of non-Gaussian statistics. The near-identical growth of non-Gaussianity in ideal and nonideal cases suggests that generation of coherent structures and breaking of self-similarity are essentially ideal processes. This has important implications for understanding the origin of intermittency in turbulence.

Initial states in incompressible two-dimensional magnetohydrodynamics that are known to lead to strong current sheets and (laminar) magnetic reconnection are modified by the addition of small-scale turbulent perturbations of various energies. The evolution of these states is computed with the aim of ascertaining the influence of the turbulence on the underlying laminar solution. Two main questions are addressed here: (1) What effect does small-scale turbulence have on the energy dissipation rate of the underlying solution? (2) What is the threshold turbulent perturbation level above which the original laminar reconnective dynamics is no longer recognizable. The simulations show that while the laminar dynamics persist the dissipation rates are largely unaffected by the turbulence, other than modest increases attributable to the additional small length scales present in the new initial condition. The solutions themselves are also remarkably insensitive to small-scale turbulent perturbations unless the perturbations are large enough to undermine the integrity of the underlying cellular flow pattern. Indeed, even initial states that lead to the evolution of small-scale microscopic sheets can survive the addition of modest turbulence. The role of a large-scale organizing background magnetic field is also addressed.

A one-point closure model for energy decay in three-dimensional magnetohydrodynamic (MHD) turbulence is developed. The model allows for influence of a large-scale magnetic field that may be of strength sufficient to induce Alfven wave propagation effects, and takes into account components of turbulence in which either the wave-like character is negligible or is dominant. This two-component model evolves energy and characteristic length scales, and may be useful as a simple description of homogeneous MHD turbulent decay. In concert with spatial transport models, it can form the basis for approximate treatment of low-frequency plasma turbulence in a variety of solar, space, and astrophysical contexts.