On the accuracy of simulations of turbulence

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

Solar wind fluctuations and the von Kármán–Howarth equations: The role of fourth-order correlations

The von Kármán-Howarth (vKH) hierarchy of equations relate the second-order correlations of the turbulent fluctuations to the third-order ones, the third-order to the fourth-order, and so on. We recently demonstrated [1] that for MHD, self-similar solutions to the vKH equations seem to require at least two independent similarity lengthscales (one for each Elsässer energy), so that compared to hydrodynamics a richer set of behaviors seems likely to ensue. Moreover, despite the well-known anisotropy of MHD turbulence with a mean magnetic field (B₀), the equation for the second-order correlation does not contain explicit dependence on B₀. We show that there is, however, implicit dependence on B₀ via the third-order correlations, which themselves have both explicit B₀-dependence and also their own implicit dependence through fourth-order correlations. Some subtleties and consequences of this implicit-explicit balance are summarized here. In addition, we present an analysis of simulation results showing that the evolution of turbulence can depend strongly on the initial fourth-order correlations of the system. This leads to considerable variation in the energy dissipation rates. Some associated consequences for MHD turbulence are discussed

The third-order law for magnetohydrodynamic turbulence with shear: Numerical investigation

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

von Kármán self-preservation hypothesis for magnetohydrodynamic turbulence and its consequences for universality

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

Is the Kelvin Theorem Valid for High-Reynolds-Number Turbulence?

The Kelvin-Helmholtz theorem on conservation of circulations is supposed to hold for ideal inviscid fluids and is believed to be play a crucial role in turbulent phenomena, such as production of dissipation by vortex line-stretching. However, this expectation does not take into account singularities in turbulent velocity fields at infinite Reynolds number. We present evidence from numerical simulations for the breakdown of the classical Kelvin theorem in the three-dimensional turbulent energy cascade. Although violated in individual realizations, we find that circulations are still conserved in some average sense. For comparison, we show that Kelvin's theorem holds for individual realizations in the two-dimensional enstrophy cascade, in agreement with theory. The turbulent ``cascade of circulations'' is shown to be a classical analogue of phase-slip due to quantized vortices in superfluids and various applications in geophysics and astrophysics are outlined.Comment: 4 pages, 3 figure

The third-order law for magnetohydrodynamic turbulence with constant shear

The scaling laws of mixed third‐order structure functions for isotropic, homogeneous, and incompressible magnetohydrodynamic (MHD) turbulence have been recently applied in solar wind studies, even though there is recognition that isotropy is not well satisfied. Other studies have taken account of the anisotropy induced by a constant mean magnetic field. However, large‐scale shear can also cause departures from isotropy. Here we examine shear effects in the simplest case, and derive the third‐order laws for MHD turbulence with constant shear, where homogeneity is still assumed. This generalized scaling law has been checked by data from direct numerical simulations (DNS) of two‐dimensional (2D) MHD and is found to hold across the inertial range. These results suggest that third‐order structure function analysis and interpretation in the solar wind should be undertaken with some caution, since, when present, shear can change the meaning of the third‐order relations

Generation of X-points and secondary islands in 2D magnetohydrodynamic turbulence

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

Evolution of a Stratified Turbulent Cloud under Rotation

Localized turbulence is common in geophysical flows, where the roles of rotation and stratification are paramount. In this study, we investigate the evolution of a stratified turbulent cloud under rotation. Recognizing that a turbulent cloud is composed of vortices of varying scales and shapes, we start our investigation with a single eddy using analytical solutions derived from a linearized system. Compared to an eddy under pure rotation, the stratified eddy shows the physical manifestation of a known potential vorticity mode, appearing as a static stable vortex. In addition, the expected shift from inertial waves to inertial-gravity waves is observed. In our numerical simulations of the turbulent cloud, carried out at a constant Rossby number over a range of Froude numbers, stratification causes columnar structures to deviate from vertical alignment. This deviation increases with increasing stratification, slowing the expansion rate of the cloud. The observed characteristics of these columnar structures are consistent with the predictions of linear theory, particularly in their tilt angles and vertical growth rates, suggesting a significant influence of inertial-gravity waves. Using Lagrangian particle tracking, we have identified regions where wave activity dominates over turbulence. In scenarios of milder stratification, these inertial-gravity waves are responsible for a significant energy transfer away from the turbulent cloud, a phenomenon that attenuates with increasing stratification

Investigation of intermittency in magnetohydrodynamics and solar wind turbulence: scale-dependent kurtosis

The behavior of scale-dependent (or filtered) kurtosis is studied in the solar wind using magnetic field measurements from the ACE and Cluster spacecraft at 1 AU. It is also analyzed numerically with high-resolution magnetohydrodynamic spectral simulations. In each case the filtered kurtosis increases with wavenumber, implying the presence of coherent structures at the smallest scales. This phase coupling is related to intermittency in solar wind turbulence and the emergence of non-Gaussian statistics. However, it is inhibited by the presence of upstream waves and other phase-randomizing structures, which act to reduce the growth of kurtosis

Influence of molecular transport on burning rate and conditioned species concentrations in highly turbulent premixed flames

Apparent inconsistency between (i) experimental and Direct Numerical Simulation (DNS) data that show the significant influence of differential diffusion on turbulent burning rate and (ii) recent complex-chemistry DNS data that indicate mitigation of the influence of differential diffusion on conditioned profiles of various local flame characteristics at high Karlovitz numbers is explored by analyzing new DNS data obtained from lean hydrogen-air turbulent flames.\ua0Both aforementioned effects are observed by analyzing the same DNS data provided that the conditioned profiles are sampled from the entire computational domain. On the contrary, the conditioned profiles sampled at the leading edge of the mean flame brush do not indicate the mitigation, but are significantly affected by differential diffusion phenomena,e.g., because reaction zones are highly curved at the leading edge. This observation is consistent with a significant increase in the computed turbulent burning velocity with decreasing Lewis number, with all the results considered jointly being consonant with the leading point concept of premixed turbulent combustion. The concept is further supported by comparing DNS data obtained by allowing for preferential diffusion solely for a single species, either atomic or molecular hydrogen