Universal statistics of non-linear energy transfer in turbulent models

A class of shell models for turbulent energy transfer at varying the inter-shell separation, $\lambda$, is investigated. Intermittent corrections in the continuous limit of infinitely close shells ($\lambda \rightarrow 1$) have been measured. Although the model becomes, in this limit, non-intermittent, we found universal aspects of the velocity statistics which can be interpreted in the framework of log-poisson distributions, as proposed by She and Waymire (1995, Phys. Rev. Lett. 74, 262). We suggest that non-universal aspects of intermittency can be adsorbed in the parameters describing statistics and properties of the most singular structure. On the other hand, universal aspects can be found by looking at corrections to the monofractal scaling of the most singular structure. Connections with similar results reported in other shell models investigations and in real turbulent flows are discussed.Comment: 4 pages, 2 figures available upon request to [email protected]

Bifractality of the Devil's staircase appearing in the Burgers equation with Brownian initial velocity

It is shown that the inverse Lagrangian map for the solution of the Burgers equation (in the inviscid limit) with Brownian initial velocity presents a bifractality (phase transition) similar to that of the Devil's staircase for the standard triadic Cantor set. Both heuristic and rigorous derivations are given. It is explained why artifacts can easily mask this phenomenon in numerical simulations.Comment: 12 pages, LaTe

A new scaling property of turbulent flows

We discuss a possible theoretical interpretation of the self scaling property of turbulent flows (Extended Self Similarity). Our interpretation predicts that, even in cases when ESS is not observed, a generalized self scaling, must be observed. This prediction is checked on a number of laboratory experiments and direct numerical simulations.Comment: Plain Latex, 1 figure available upon request to [email protected]

Developed turbulence: From full simulations to full mode reductions

Developed Navier-Stokes turbulence is simulated with varying wavevector mode reductions. The flatness and the skewness of the velocity derivative depend on the degree of mode reduction. They show a crossover towards the value of the full numerical simulation when the viscous subrange starts to be resolved. The intermittency corrections of the scaling exponents of the pth order velocity structure functions seem to depend mainly on the proper resolution of the inertial subrange. Universal scaling properties (i.e., independent of the degree of mode reduction) are found for the relative scaling exponents rho which were recently defined by Benzi et al.Comment: 4 pages, 5 eps-figures, replaces version from August 5th, 199

Turbulence and Multiscaling in the Randomly Forced Navier Stokes Equation

We present an extensive pseudospectral study of the randomly forced Navier-Stokes equation (RFNSE) stirred by a stochastic force with zero mean and a variance $\sim k^{4-d-y}$, where $k$ is the wavevector and the dimension $d = 3$. We present the first evidence for multiscaling of velocity structure functions in this model for $y \geq 4$. We extract the multiscaling exponent ratios $\zeta_p/\zeta_2$ by using extended self similarity (ESS), examine their dependence on $y$, and show that, if $y = 4$, they are in agreement with those obtained for the deterministically forced Navier-Stokes equation ($3d$NSE). We also show that well-defined vortex filaments, which appear clearly in studies of the $3d$NSE, are absent in the RFNSE.Comment: 4 pages (revtex), 6 figures (postscript

Universality in fully developed turbulence

We extend the numerical simulations of She et al. [Phys.\ Rev.\ Lett.\ 70, 3251 (1993)] of highly turbulent flow with $15 \le$ Taylor-Reynolds number $Re_\lambda\le 200$ up to $Re_\lambda \approx 45000$, employing a reduced wave vector set method (introduced earlier) to approximately solve the Navier-Stokes equation. First, also for these extremely high Reynolds numbers $Re_\lambda$, the energy spectra as well as the higher moments -- when scaled by the spectral intensity at the wave number $k_p$ of peak dissipation -- can be described by {\it one universal} function of $k/k_p$ for all $Re_\lambda$. Second, the ISR scaling exponents $\zeta_m$ of this universal function are in agreement with the 1941 Kolmogorov theory (the better, the large $Re_\lambda$ is), as is the $Re_\lambda$ dependence of $k_p$. Only around $k_p$ viscous damping leads to slight energy pileup in the spectra, as in the experimental data (bottleneck phenomenon).Comment: 14 pages, Latex, 5 figures (on request), 3 tables, submitted to Phys. Rev.

Supersonic turbulence and structure of interstellar molecular clouds

The interstellar medium (ISM) provides a unique laboratory for highly supersonic, driven hydrodynamics turbulence. We present a theory of such turbulence, confirm it by numerical simulations, and use the results to explain observational properties of interstellar molecular clouds, the regions where stars are born.Comment: 5 pages, 3 figures include

Finite size corrections to scaling in high Reynolds number turbulence

We study analytically and numerically the corrections to scaling in turbulence which arise due to the finite ratio of the outer scale $L$ of turbulence to the viscous scale $\eta$, i.e., they are due to finite size effects as anisotropic forcing or boundary conditions at large scales. We find that the deviations \dzm from the classical Kolmogorov scaling $\zeta_m = m/3$ of the velocity moments \langle |\u(\k)|^m\rangle \propto k^{-\zeta_m} decrease like $\delta\zeta_m (Re) =c_m Re^{-3/10}$. Our numerics employ a reduced wave vector set approximation for which the small scale structures are not fully resolved. Within this approximation we do not find $Re$ independent anomalous scaling within the inertial subrange. If anomalous scaling in the inertial subrange can be verified in the large $Re$ limit, this supports the suggestion that small scale structures should be responsible, originating from viscosity either in the bulk (vortex tubes or sheets) or from the boundary layers (plumes or swirls)

On the stationarity of linearly forced turbulence in finite domains

A simple scheme of forcing turbulence away from decay was introduced by Lundgren some time ago, the `linear forcing', which amounts to a force term linear in the velocity field with a constant coefficient. The evolution of linearly forced turbulence towards a stationary final state, as indicated by direct numerical simulations (DNS), is examined from a theoretical point of view based on symmetry arguments. In order to follow closely the DNS the flow is assumed to live in a cubic domain with periodic boundary conditions. The simplicity of the linear forcing scheme allows one to re-write the problem as one of decaying turbulence with a decreasing viscosity. Scaling symmetry considerations suggest that the system evolves to a stationary state, evolution that may be understood as the gradual breaking of a larger approximate symmetry to a smaller exact symmetry. The same arguments show that the finiteness of the domain is intimately related to the evolution of the system to a stationary state at late times, as well as the consistency of this state with a high degree of isotropy imposed by the symmetries of the domain itself. The fluctuations observed in the DNS for all quantities in the stationary state can be associated with deviations from isotropy. Indeed, self-preserving isotropic turbulence models are used to study evolution from a direct dynamical point of view, emphasizing the naturalness of the Taylor microscale as a self-similarity scale in this system. In this context the stationary state emerges as a stable fixed point. Self-preservation seems to be the reason behind a noted similarity of the third order structure function between the linearly forced and freely decaying turbulence, where again the finiteness of the domain plays an significant role.Comment: 15 pages, 7 figures, changes in the discussion at the end of section VI, formula (60) correcte

Circulation Statistics in Three-Dimensional Turbulent Flows

We study the large $\lambda$ limit of the loop-dependent characteristic functional $Z(\lambda)=$, related to the probability density function (PDF) of the circulation around a closed contour $c$. The analysis is carried out in the framework of the Martin-Siggia-Rose field theory formulation of the turbulence problem, by means of the saddle-point technique. Axisymmetric instantons, labelled by the component $\sigma_{zz}$ of the strain field -- a partially annealed variable in our formalism -- are obtained for a circular loop in the $xy$ plane, with radius defined in the inertial range. Fluctuations of the velocity field around the saddle-point solutions are relevant, leading to the lorentzian asymptotic behavior $Z(\lambda) \sim 1/{\lambda^2}$. The ${\cal O}(1 / {\lambda^4})$ subleading correction and the asymmetry between right and left PDF tails due to parity breaking mechanisms are also investigated.Comment: Computations are discussed in a more detailed way; accepted for publication in Physical Review