Search

In an earlier paper (Wan et al. 2012), the authors showed that a similarity solution for anisotropic incompressible 3D magnetohydrodynamic (MHD) turbulence, in the presence of a uniform mean magnetic field $\vB_0$, exists if the ratio of parallel to perpendicular (with respect to $\vB_0$) similarity length scales remains constant in time. This conjecture appears to be a rather stringent constraint on the dynamics of decay of the energy-containing eddies in MHD turbulence. However, we show here, using direct numerical simulations, that this hypothesis is indeed satisfied in incompressible MHD turbulence. After an initial transient period, the ratio of parallel to perpendicular length scales fluctuates around a steady value during the decay of the eddies. We show further that a Taylor--K\'arm\'an-like similarity decay holds for MHD turbulence in the presence of a mean magnetic field. The effect of different parameters, including Reynolds number, DC field strength, and cross-helicity, on the nature of similarity decay is discussed., Comment: Accepted for publication in Journal of Fluid Mechanics

When magnetohydrodynamic turbulence evolves in the presence of a large-scale mean magnetic field, an anisotropy develops relative to that preferred direction. The well-known tendency is to develop stronger gradients perpendicular to the magnetic field, relative to the direction along the field. This anisotropy of the spectrum is deeply connected with anisotropy of estimated timescales for dynamical processes, and requires reconsideration of basic issues such as scale locality and spectral transfer. Here analysis of high-resolution three-dimensional simulations of unforced magnetohydrodynamic turbulence permits quantitative assessment of the behavior of theoretically relevant timescales in Fourier wavevector space. We discuss the distribution of nonlinear times, Alfv\'en times, and estimated spectral transfer rates. Attention is called to the potential significance of special regions of the spectrum, such as the two-dimensional limit and the "critical balance" region. A formulation of estimated spectral transfer in terms of a suppression factor supports a conclusion that the quasi two-dimensional fluctuations (characterized by strong nonlinearities) are not a singular limit, but may be in general expected to make important contributions., Comment: In Press at Physics of Plasmas

In presence of an externally supported, mean magnetic field a turbulent, conducting medium, such as plasma, becomes anisotropic. This mean magnetic field, which is separate from the fluctuating, turbulent part of the magnetic field, has considerable effects on the dynamics of the system. In this paper, we examine the dissipation rates for decaying incompressible magnetohydrodynamic (MHD) turbulence with increasing Reynolds number, and in the presence of a mean magnetic field of varying strength. Proceeding numerically, we find that as the Reynolds number increases, the dissipation rate asymptotes to a finite value for each magnetic field strength, confirming the K\'arm\'an-Howarth hypothesis as applied to MHD. The asymptotic value of the dimensionless dissipation rate is initially suppressed from the zero-mean-field value by the mean magnetic field but then approaches a constant value for higher values of the mean field strength. Additionally, for comparison, we perform a set of two-dimensional (2DMHD) and a set of reduced MHD (RMHD) simulations. We find that the RMHD results lie very close to the values corresponding to the high mean-field limit of the three-dimensional runs while the 2DMHD results admit distinct values far from both the zero mean field cases and the high mean field limit of the three-dimensional cases. These findings provide firm underpinnings for numerous applications in space and astrophysics wherein von K\'arm\'an decay of turbulence is assumed., Comment: Accepted for publication in Physical Review X

We argue that the hypothesis of preservation of shape of dimensionless second- and third-order correlations during decay of incompressible homogeneous magnetohydrodynamic (MHD) turbulence requires, in general, at least two independent similarity length scales. These are associated with the two Elsässer energies. The existence of similarity solutions for the decay of turbulence with varying cross-helicity implies that these length scales cannot remain in proportion, opening the possibility for a wide variety of decay behaviour, in contrast to the simpler classic hydrodynamics case. Although the evolution equations for the second-order correlations lack explicit dependence on either the mean magnetic field or the magnetic helicity, there is inherent implicit dependence on these (and other) quantities through the third-order correlations. The self-similar inertial range, a subclass of the general similarity case, inherits this complexity so that a single universal energy spectral law cannot be anticipated, even though the same pair of third-order laws holds for arbitrary cross-helicity and magnetic helicity. The straightforward notion of universality associated with Kolmogorov theory in hydrodynamics therefore requires careful generalization and reformulation in MHD.

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...

A phenomenological theory describes radial evolution of plasma turbulence in the solar wind from 1 to 50 astronomical units. The theory includes a simple closure for local anisotropic magnetohydrodynamic turbulence, spatial transport, and driving by large-scale shear and pickup ions. Results compare well to plasma and magnetic field data from the Voyager 2 spacecraft, providing a basis for a concise, tractable description of turbulent energy transport in a variety of astrophysical plasmas.

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.

In an evolving turbulent medium, a natural timescale can be defined in terms of the energy decay time. The time evolution may be complicated by other effects such as energy supply due to driving, and spatial inhomogeneity. In the solar wind the turbulence appears not to be simply engaging in free decay, but rather the energy level observed at a particular position in the hellosphere is affected by expansion, "mixing," and driving by stream shear. Here we discuss a new approach for estimating the "age" of solar wind turbulence as a function of heliocentric distance, using the local turbulent decay rate as the natural dock, but taking into account expansion and driving effects. The simplified formalism presented here is appropriate to low cross hellcity (non-Alfvnic) turbulence in the outer hellosphere especially at low hello-latitudes. We employ Voyager data to illustrate our method, which improves upon the familiar estimates in terms of local eddy turnover times.

The properties and structure of second-order ~Cartesian! correlation tensors are derived for the general case of two solenoidal random vector fields. The theory is intended to describe homogeneous magnetohydrodynamic turbulence, with no assumed rotational or reflectional symmetries. Each correlation tensor can be written in terms of four scalar generating functions and the relationship of these functions to the potentials that generate the poloidal and toroidal components of the underlying vector fields is derived. The physical nature of the scalar functions is investigated and their true or pseudoscalar character is ascertained. In our general discussion we clarify several misleading statements dating back to Robertson’s original paper in the field @Proc. Camb. Philos. Soc. 36, 209 ~1940!#. It is also shown that using the one-dimensional correlation function, it is possible to obtain spectral information on the induced electric field in directions perpendicular to the measurement direction. @S1063-651X~97!09208-8#

Solar wind fluctuations admit well-documented anisotropies of the variance matrix, or polarization, related to the mean magnetic field direction. Typically, one finds a ratio of perpendicular variance to parallel variance of the order of 9:1 for the magnetic field. Here we study the question of whether a kinetic plasma spontaneously generates and sustains parallel variances when initiated with only perpendicular variance. We find that parallel variance grows and saturates at about 5% of the perpendicular variance in a few nonlinear times irrespective of the Reynolds number. For sufficiently large systems (Reynolds numbers) the variance approaches values consistent with the solar wind observations.

Building on results from two-dimensional magnetohydrodynamic (MHD) turbulence (Shebalin, Matthaeus & Montgomery 1983), the development of anisotropic states from initially isotropic ones is investigated numerically for fully three-dimensional incompressible MHD turbulence. It is found that when an external d.c. magnetic field (B0) is imposed on viscous and resistive MHD systems, excitations are preferentially transferred to modes with wavevectors perpendicular to B0). The anisotropy increases with increasing mechanical and magnetic Reynolds numbers, and also with increasing wavenumber. The tendency of B0 to inhibit development of turbulence is also examined.

We observe the anisotropy of the power spectral tensor of magnetic field fluctuations in the fast solar wind for the first time. In heliocentric RTN coordinates the power in each element of the tensor has a unique dependence on the angle between the magnetic field and velocity of the solar wind (\theta) and the angle of the vector in the plane perpendicular to the velocity (\phi). We derive the geometrical effect of the high speed flow of the solar wind past the spacecraft on the power spectrum in the frame of the plasma P(k) to arrive at the observed power spectrum P(f,\theta,\phi) based on a scalar field description of turbulence theory. This allows us to predict the variation in the \phi direction and compare it to the data. We then transform the observations from RTN coordinates to magnetic-field-aligned coordinates. The observed reduced power spectral tensor matches the theoretical predictions we derive in both RTN and field-aligned coordinates which means that the local magnetic field we calculate with wavelet envelope functions is an accurate representation of the physical axis of symmetry for the turbulence and implies that on average the turbulence is axi-symmetric. We also show that we can separate the dominant toroidal component of the turbulence from the smaller but significant poloidal component and that these have different power anisotropy. We also conclude that the magnetic helicity is anisotropic and mostly two-dimensional, arising from wavevectors largely confined to the plane perpendicular to B., Comment: 21 pages, 8 figures, 1 table

Numerical solution of two-dimensional incompressible hydrodynamics shows that states of a near-minimal ratio of enstrophy to energy can be attained in times short compared with the flow decay time, confirming the simplest turbulent selective decay conjecture, and suggesting that coherent vortex structures do not terminate nonlinear processes. After all possible vortex mergers occur, the vorticity attains a particlelike character, suggested by the late-time similarity of the streamlines to Ewald potential contours.

We present a quantitative estimate of the anisotropic power and scaling of magnetic field fluctuations in inertial range magnetohydrodynamic turbulence, using a novel wavelet technique applied to spacecraft measurements in the solar wind. We show for the first time that, when the local magnetic field direction is parallel to the flow, the spacecraft-frame spectrum has a spectral index near 2. This can be interpreted as the signature of a population of fluctuations in field-parallel wavenumbers with a k −2 k spectrum but is also consistent with the presence of a “critical balance” style turbulent cascade. We also find, in common with previous studies, that most of the power is contained in wavevectors at large angles to the local magnetic field and that this component of the turbulence has a spectral index of 5/3. Magnetised plasmas fill most of the Universe and in many regions, turbulence plays important roles in the transport of energy and momentum and the acceleration and scattering of charged particles. Many aspects of plasma turbulence remain poorly understood, however. Here we present results on one of these, the anisotropy of the energy spectrum of magnetohydrodynamic (MHD) turbulence with respect to the magnetic field. In classical hydrodynamics, velocity fluctuations δuk with a wavenumber k decay and transfer energy to smaller scales on the shear timescale, τS ≈ 1/(kδuk). Within the steady inertial range, far from the energy input (“outer”) and dissipation scales, this leads to the dimensional result (δuk) 3 ∝ ǫ/k, where ǫ is the energy dissipation rate per unit mass. This gives the familiar Kolmogorov energy spectrum P(k) ∝ k −5/3 , widely observed in hydrodynamic turbulence. In a plasma, fluc⊥ ! .

[1] We examine numerically the influence of Hall effect corrections to Ohm's law upon the decay of homogeneous compressible magnetohydrodynamic turbulence and conclude that there are no significant differences in global decay rate associated with the Hall effect. This affirms expectations that energy decay is controlled by the large-scale eddies.

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., 19 pages, 5 figures. Accepted for publication in ApJ

We study the statistics of coherent current sheets, the population of X-type critical points, and reconnection rates in a coronal loop geometry, via numerical simulations of reduced magnetohydrodynamic turbulence. Current sheets and sites of reconnection (magnetic X-points) are identified in two-dimensional planes of the three-dimensional simulation domain. The geometry of the identified current sheets—including area, length, and width—and the magnetic dissipation occurring in the current sheets are statistically characterized. We also examine the role of magnetic reconnection, by computing the reconnection rates at the identified X-points and investigating their association with current sheets.

The random walk of magnetic field lines is examined numerically and analytically in the context of reduced magnetohydrodynamic (RMHD) turbulence, which provides a useful description of plasmas dominated by a strong mean field, such as in the solar corona. A recently developed non-perturbative theory of magnetic field line diffusion is compared with the diffusion coefficients obtained by accurate numerical tracing of magnetic field lines for both synthetic models and direct numerical simulations of RMHD. Statistical analysis of an ensemble of trajectories confirms the applicability of the theory, which very closely matches the numerical field line diffusion coefficient as a function of distance z along the mean magnetic field for a wide range of the Kubo number R. This theory employs Corrsin's independence hypothesis, sometimes thought to be valid only at low R. However, the results demonstrate that it works well up to R = 10, both for a synthetic RMHD model and an RMHD simulation. The numerical results from the RMHD simulation are compared with and without phase randomization, demonstrating a clear effect of coherent structures on the field line random walk for a very low Kubo number.

Evolution of magnetohydrodynamic turbulence is often discussed in terms of second-order statistics like the energy spectra, but consideration of the structure of the von K´ arm´ an–Howarth hierarchy for MHD indicates that higherorder statistical correlations occupy an influential role. Here we show that both spectral anisotropy and energy decay are strongly associated with higher-order statistics. Dynamical emergence of spectral anisotropy must occur at a higher order in the statistical hierarchy, while numerical evidence suggests that strong variations in energy decay are connected with variations in higher-order statistics.

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.

The magnetic and velocity field fluctuations in magnetohydrodynamic turbulence can be characterized by their directional alignment and induced electric field. These manifest as coherent spatial correlations which are measures of Alfvenicity and turbulence cascade strength, respectively. Solar wind observations and direct numerical simulations find that these distinctive correlations, caused by rapid relaxation processes that act to suppress nonlinearity, occur in localized spatial patches. This cellularization of magnetofluid turbulence is inconsistent with a superposition of Gaussian fields and could be related to spatial intermittency or other non-Gaussian statistics.

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 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 present various numerical and analytical solutions for the transport of solar wind turbulence. The model used takes into account the effects of convection, expansion, and wave propagation, as well as the recently illuminated effects of (non-WKB) 'mixing' terms. The radial evolution of the fluctuating kinetic energy, magnetic energy, and normalized cross helicity is computed, and, it is demonstrated that in appropriate limits the solutions converge to the WKB forms. In the general case, solutions which differ substantially from those predicted by WKB theory are obtained. The degree of turbulent 'mixing' shows considerable dependence on the nature of the turbulence, giving rise to varying levels, at 1 AU, of the ratio of 'inward' and 'outward' fluctuation energies and the ratio of kinetic to magnetic energies in the fluctuations. The transport properties described here may provide at least a partial explanation for the observed mixing of cross helicities with increasing heliocentric distance in the solar wind.

An evaluation is conducted of recent theories of small-scale MHD turbulence transport in an inhomogeneous background that are pertinent to the evolution of solar wind turbulence. Attention is given to the WKB formalism that has been used in many solar wind-related physics applications, with a view to its shortcomings. Also discussed are the structure of two-scale transport theories, and their relationship to WKB theory in light of multiple-scales analysis.

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.

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.

Previous theoretical treatments of fluid-scale turbulence in the solar wind have concentrated on describing the state and dynamical evolution of fluctuations in the inertial range, which are characterized by power law energy spectra. In the present paper a model for the evolution of somewhat larger, more energetic magnetohydrodynamic (MHD) fluctuations is developed by analogy with classical hydrodynamic turbulence in the quasi-equilibrium range. The model is constructed by assembling and extending existing phenomenologies of homogeneous MHD turbulence, as well as simple two-length-scale models for transport of MHD turbulence in a weekly inhomogeneous medium. A set of equations is presented for the evolution of the turbulence, including the transport and nonlinear evolution of magnetic and kinetic energy, cross helicity, and their correlation scales. Two versions of the model are derived, depending on whether the fluctuations are distributed isotropically in three dimensions or restricted to the two-dimensional plane perpendicular to the mean magnetic field. This model includes a number of potentially important physical effects that have been neglected in previous discussions of transport of solar wind turbulence.

An alternative approach, based on a multiple scale analysis, is presented in order to reconcile the traditional Wentzel-Kramer-Brillouin (WKB) approach to the modeling of interplanetary fluctuations in a mildly inhomogeneous large-scale flow with a more recently developed transport theory. This enables us to compare directly, at a formal level, the inherent structure of the two models. In the case of noninteracting, incompressible (Alven) waves, the principle difference between the two models is the presence of leading-order couplings (called 'mixing effects') in the non-WKB turbulence model which are absent in a WKB development. Within the context of linearized MHD, two cases have been identified for which the leading order non-WJB 'mixing term' does not vanish at zero wavelength. For these cases the WKB expansion is divergent, whereas the multiple-scale theory is well behaved. We have thus established that the WKB results are contained within the multiple-scale theory, but leading order mixing effects, which are likely to have important observational consequences, can never be recovered in the WKB style expansion. Properties of the higher-order terms in each expansion are also discussed, leading to the conclusion that the non-WKB hierarchy may be applicable even when the scale separation parameter is not small.

A candidate mechanism for the heating of the solar corona in open field line regions is described. The interaction of Alfven waves, generated in the photosphere or chromosphere, with their reflections and the subsequent driving of quasi-two-dimensional MHD turbulence is considered. A nonlinear cascade drives fluctuations toward short wavelengths which are transverse to the mean field, thereby heating at rates insensitive to restrictive Alfven timescales. A phenomenology is presented, providing estimates of achievable heating efficiency that are most favorable.

A computational investigation of magnetic helicity of the fluctuating magnetic field Hm in ideal and freely decaying three‐dimensional (3‐D) magnetohydrodynamics (MHD) in the presence of a uniform mean magnetic field is performed. It is shown that for ideal 3‐D MHD Hm, which is a rugged invariant in the absence of a mean magnetic field [Frisch et al., J. Fluid Mech. 77, 796 (1975)], decays from its initial value and proceeds to oscillate about zero. The decay of Hm is shown to result from the presence of a new ‘‘generalized’’ helicity invariant, which includes contributions from the uniform magnetic field. The loss of invariance of Hm will diminish the effects of inverse transfer of Hm on freely decaying turbulence. This is demonstrated in a discussion of the selective decay relaxation process.