Energy decay laws in strongly anisotropic magnetohydrodynamic turbulence

We investigate the influence of a uniform magnetic field B(0)=B(0)e( parallel) on energy decay laws in incompressible magnetohydrodynamic (MHD) turbulence. The nonlinear transfer reduction along B(0) is included in a model that distinguishes parallel and perpendicular directions, following a phenomenology of Kraichnan. We predict a slowing down of the energy decay due to anisotropy in the limit of strong B(0), with distinct power laws for energy decay of shear- and pseudo-Alfvén waves. Numerical results from the kinetic equations of Alfvén wave turbulence recover these predictions, and MHD numerical results clearly tend to follow them in the lowest perpendicular planes.     

Effect of current sheets on the solar wind magnetic field power spectrum from the Ulysses observation: from Kraichnan to Kolmogorov scaling

The MHD turbulence theory developed by Iroshnikov and Kraichnan predicts a k(-1.5) power spectrum. Solar wind observations, however, often show a k(-5/3) Kolmogorov scaling. Based on 3 years worth of Ulysses magnetic field data where over 28,000 current sheets are identified, we propose that the current sheet is the cause of the Kolmogorov scaling. We show that for 5 longest current-sheet-free periods the magnetic field power spectra are all described by the Iroshnikov-Kraichnan scaling. In comparison, for 5 periods that have the most number of current sheets, the power spectra all exhibit Kolmogorov scaling. The implication of our results is discussed.     

Observation of inertial energy cascade in interplanetary space plasma

Direct evidence for the presence of an inertial energy cascade, the most characteristic signature of hydromagnetic turbulence (MHD), is observed in the solar wind by the Ulysses spacecraft. After a brief rederivation of the equivalent of Yaglom's law for MHD turbulence, a linear relation is indeed observed for the scaling of mixed third-order structure functions involving Els?sser variables. This experimental result firmly establishes the turbulent character of low-frequency velocity and magnetic field fluctuations in the solar wind plasma.     

Scaling properties of three-dimensional magnetohydrodynamic turbulence

The scaling properties of three-dimensional magnetohydrodynamic turbulence with finite magnetic helicity are obtained from direct numerical simulations using 512(3) modes. The results indicate that the turbulence does not follow the Iroshnikov-Kraichnan phenomenology. The scaling exponents of the structure functions can be described by a modified She-Leveque model zeta(p) = p/9+1-(1/3)(p/3), corresponding to basic Kolmogorov scaling and sheetlike dissipative structures. In particular, we find zeta(2) approximately 0.7, consistent with the energy spectrum E(k) approximately k(-5/3) as observed in the solar wind, and zeta(3) approximately 1, confirming a recent analytical result.     

Quasi-Gaussian statistics of hydrodynamic turbulence in 4 / 3 +epsilon dimensions

The statistics of two-dimensional turbulence exhibit a riddle: the scaling exponents in the regime of inverse energy cascade agree with the K41 theory of turbulence far from equilibrium, but the probability distribution functions are close to Gaussian-like in equilibrium. The skewness S identical with S3(R)/S(3/2)(2)(R) was measured as S (exp) approximately 0.03. This contradiction is lifted by understanding that two-dimensional turbulence is not far from a situation with equipartition of enstrophy, which exists as true thermodynamic equilibrium with K41 exponents in space dimension of d= 4 / 3. We evaluate the skewness S( d) for 4 / 3 < or =d< or =2, showing that S(d)=0 at d= 4 / 3, and that it remains as small as S (exp) in two dimensions.     

Spectral energy dynamics in magnetohydrodynamic turbulence

Spectral direct numerical simulations of incompressible MHD turbulence at a resolution of up to 1024(3) collocation points are presented for a statistically isotropic system as well as for a setup with an imposed strong mean magnetic field. The spectra of residual energy, E(R)k=|E(M)k - E(K)k|, and total energy, Ek=E(K)k+E(M)k, are observed to scale self-similarly in the inertial range as E(R)k approximately k(-7/3), E(k)approximately k(-5/3) (isotropic case) and E(R)(k(perpendicular) approximately k(-2)(perpendicular), E(k(perpendicular))approximately k(-3/2)(perpendicular) (anisotropic case, perpendicular to the mean field direction). A model of dynamic equilibrium between kinetic and magnetic energy, based on the corresponding evolution equations of the eddy-damped quasinormal Markovian closure approximation, explains the findings. The assumed interplay of turbulent dynamo and Alfvén effect yields E(R)k approximately kE2(k), which is confirmed by the simulations.     

Universal spectrum of two-dimensional turbulence on a rotating sphere and some basic features of atmospheric circulation on giant planets

The Kolmogorov-Batchelor-Kraichnan (KBK) theory of two-dimensional turbulence is generalized for turbulence on the surface of a rotating sphere. The energy spectrum develops considerable anisotropy; a steep -5 slope emerges in the zonal direction, while in all others the classical KBK scaling prevails. This flow regime in robust steady state is reproduced in simulations with linear drag. The conditions favorable for this regime may be common for giant planets' atmospheric circulations; the same steep spectra are found in their observed zonal velocity profiles and utilized to explain their basic characteristics.     

Magnetic discontinuities in magnetohydrodynamic turbulence and in the solar wind

Recent measurements of solar wind turbulence report the presence of intermittent, exponentially distributed angular discontinuities in the magnetic field. In this Letter, we study whether such discontinuities can be produced by magnetohydrodynamic (MHD) turbulence. We detect the discontinuities by measuring the fluctuations of the magnetic field direction, Δθ, across fixed spatial increments Δx in direct numerical simulations of MHD turbulence with an imposed uniform guide field B(0). A large region of the probability density function (pdf) for Δθ is found to follow an exponential decay, proportional to exp(-Δθ/θ(*)), with characteristic angle θ(*)≈(14°)(b(rms)/B(0))(0.65) for a broad range of guide-field strengths. We find that discontinuities observed in the solar wind can be reproduced by MHD turbulence with reasonable ratios of b(rms)/B(0). We also observe an excess of small angular discontinuities when Δx becomes small, possibly indicating an increasing statistical significance of dissipation-scale structures. The structure of the pdf in this case closely resembles the two-population pdf seen in the solar wind. We thus propose that strong discontinuities are associated with inertial-range MHD turbulence, while weak discontinuities emerge from dissipation-range turbulence. In addition, we find that the structure functions of the magnetic field direction exhibit anomalous scaling exponents, which indicates the existence of intermittent structures.     

Restricted scaling range models for turbulent velocity and scalar energy transfers in decaying turbulence

The effect of finite Reynolds numbers and/or internal intermittency on the total kinetic energy and scalar energy transfers is examined in detail. For this purpose, two distinct models for velocity and scalar energy transfer are proposed in the specific context of freely decaying isotropic turbulence. The first one extends the already existing dynamical models (hereafter DYM, i.e. based on transport equations originated in Navier–Stokes and advection-diffusion transport equations). The second one relies on the characteristic time of the strain at a specific scale (hereafter SBM). Both models account for the Reynolds number dependence of the scaling exponent of the second-order structure functions, over a range of scales where such exponents may be defined, i.e. a restricted scaling range (RSR). Therefore, the models developed aim at reproducing the energy transfer over the RSR. The predicted energy transfer is very sensible to variations of the scaling exponent, especially at low Reynolds numbers. The approach towards the asymptotic 4/3 law is closely reproduced by the two models. The dynamical model reproduces the experimental results accurately especially in the vicinity of the Taylor microscale, while the SBM agrees almost perfectly with measurements at nearly all scales.     

Collisionless magnetohydrodynamic turbulence in two dimensions

In this paper we give a formulation of two-dimensional (2D) collisionless magnetohydrodynamic (MHD) turbulence that includes the effects of both electron inertia and electron pressure (or parallel electron compressibility) and is applicable to strongly magnetized collisionless plasmas. We place particular emphasis on the departures from the 2D classical MHD turbulence results produced by the collisionless MHD effects. We investigate the fractal/multi-fractal aspects of spatial intermittency. The fractal model for intermittent collisionless MHD turbulence appears to be able to describe the observed k⁻¹ spectrum in the solar wind. Multi-fractal scaling behaviors in the inertial range are first deduced, and are then extrapolated down to the dissipative microscales. We then consider a parabolic-profile model for the singularity spectrum f (α), as an explicit example of a multi-fractal scenario. These considerations provide considerable insights into the basic mechanisms underlying spatial intermittency in 2D fully developed collisionless MHD turbulence.     

Exact relation for correlation functions in compressible isothermal turbulence

Compressible isothermal turbulence is analyzed under the assumption of homogeneity and in the asymptotic limit of a high Reynolds number. An exact relation is derived for some two-point correlation functions which reveals a fundamental difference with the incompressible case. The main difference resides in the presence of a new type of term which acts on the inertial range similarly as a source or a sink for the mean energy transfer rate. When isotropy is assumed, compressible turbulence may be described by the relation -2/3ε(eff)r = F(r)(r), where F(r) is the radial component of the two-point correlation functions and ε(eff) is an effective mean total energy injection rate. By dimensional arguments, we predict that a spectrum in k(-5/3) may still be preserved at small scales if the density-weighted fluid velocity ρ(1/3)u is used.     

Exchange-Driven Growth with Birth Rate Less Than Death

We further study the kinetic behavior of the exchange-driven growth with birth and death for the case of birth rate kernel being less than that of death based on the mean-field theory. The symmetric exchange rate kernel is K(k,j)=K'(k,j)=Ikj^υ, and the birth and death rates are proportional to the aggregate's size. The long time asymptotic behavior of the aggregate size distribution a_k(t) is found to obey a much unusual scaling law with an exponentially growing scaling function Φ(x)=exp(x).     

非均衡湍流研究进展

Richardson-Kolmogorov能量级串理论是湍流研究中最重要的基础理论,其中一个推论是能量的传输和耗散应当是均衡的,对应于耗散系数${C_\mathit{\boldsymbol{\epsilon}}}$为常数.然而近些年的实验及数值模拟都发现了不符合Richardson-Kolmogorov能量级串理论的非均衡耗散律,即使在此区域内Reynolds数足够高,能谱满足Kolmogorov的-5/3标度律,${C_\mathit{\boldsymbol{\epsilon}}}$也不为常数,而满足${C_\mathit{\boldsymbol{\epsilon}}}$~Re_Im/Re_Ln,其中m≈1≈n,Re_I为入口Reynolds数,Re_L为以积分尺度为特征长度的当地Reynolds数.近三年来又发现流向速度梯度扭率S_k和Lagrange速度梯度自相关的时间演化Φ'_ijij也可以用来度量非均衡湍流现象,为非均衡湍流的研究开辟了新路.      

Unifying scaling theory for vortex dynamics in two-dimensional turbulence

We present a scaling theory for unforced inviscid two-dimensional turbulence. Our model unifies existing spatial and temporal scaling theories. The theory is based on a self-similar distribution of vortices of different sizes A. Our model uniquely determines the spatial and temporal scaling of the associated vortex number density which allows the determination of the energy spectra and the vortex distributions. We find that the vortex number density scales as n(A,t)-t(-2/3)/A, which implies an energy spectrum E-k(-5), significantly steeper than the classical Batchelor-Kraichnan scaling. High-resolution numerical simulations corroborate the model.     

Pressure spectrum in homogeneous turbulence

The pressure spectrum in homogeneous steady turbulence is studied using direct numerical simulation with resolution up to 1024(3) and the Reynolds number R(lambda) between 38 and 478. The energy spectrum is found to have a finite inertial range with the Kolmogorov constant K = 1.65+/-0.05 followed by a bump at large wave numbers. The pressure spectrum in the inertial range is found to be approximately P(k) = B(p)epsilon;(4/3)k(-7/3) with B(p) = 8.0+/-0.5, and followed by a bump of nearly k(-5/3) at higher wave numbers. Universality and a new scaling of the pressure spectrum are discussed.     

Self-similar signature of the active solar corona within the inertial range of solar-wind turbulence

We quantify the scaling of magnetic energy density in the inertial range of solar-wind turbulence seen in situ at 1 AU with respect to solar activity. At solar maximum, when the coronal magnetic field is dynamic and topologically complex, we find self-similar scaling in the solar wind, whereas at solar minimum, when the coronal fields are more ordered, we find multifractality. This quantifies the solar-wind signature that is of direct coronal origin and distinguishes it from that of local MHD turbulence, with quantitative implications for coronal heating of the solar wind.     

On Yaglom’s Law for the Interplanetary Proton Density and Temperature Fluctuations in Solar Wind Turbulence

In the past decades, there has been an increasing literature on the presence of an inertial energy cascade in interplanetary space plasma, being interpreted as the signature of Magnetohydrodynamic turbulence (MHD) for both fields and passive scalars. Here, we investigate the passive scalar nature of the solar wind proton density and temperature by looking for scaling features in the mixed-scalar third-order structure functions using measurements on-board the Ulysses spacecraft during two different periods, i.e., an equatorial slow solar wind and a high-latitude fast solar wind, respectively. We find a linear scaling of the mixed third-order structure function as predicted by Yaglom’s law for passive scalars in the case of slow solar wind, while the results for fast solar wind suggest that the mixed fourth-order structure function displays a linear scaling. A simple empirical explanation of the observed difference is proposed and discussed.     

Anisotropy of imbalanced Alfvénic turbulence in fast solar wind

We present the first measurement of the scale-dependent power anisotropy of Elsasser variables in imbalanced fast solar wind turbulence. The dominant Elsasser mode is isotropic at lower spacecraft frequencies but becomes increasingly anisotropic at higher frequencies. The subdominant mode is anisotropic throughout. There are two distinct subranges exhibiting different scalings within what is normally considered the inertial range. The low Alfvén ratio and the different scaling of the Elsasser modes suggests an interpretation of the observed discrepancy between the velocity and magnetic field scalings, the total energy is dominated by the latter. These results do not appear to be fully explained by any of the current theories of incompressible imbalanced MHD turbulence.     

Measurement of the electric fluctuation spectrum of magnetohydrodynamic turbulence

Magnetohydrodynamic (MHD) turbulence in the solar wind is observed to show the spectral behavior of classical Kolmogorov fluid turbulence over an inertial subrange and departures from this at short wavelengths, where energy should be dissipated. Here we present the first measurements of the electric field fluctuation spectrum over the inertial and dissipative wave number ranges in a Beta > or approximately = 1 plasma. The k(-5/3) inertial subrange is observed and agrees strikingly with the magnetic fluctuation spectrum; the wave phase speed in this regime is shown to be consistent with the Alfvén speed. At smaller wavelengths krho(i) > or = 1 the electric spectrum is enhanced and is consistent with the expected dispersion relation of short-wavelength kinetic Alfvén waves. Kinetic Alfvén waves damp on the solar wind ions and electrons and may act to isotropize them. This effect may explain the fluidlike nature of the solar wind.