No Arabic abstract
Intense fluctuations of energy dissipation rate in turbulent flows result from the self-amplification of strain rate via a quadratic nonlinearity, with contributions from vorticity (via the vortex stretching mechanism) and the pressure Hessian tensor, which we analyze here using direct numerical simulations of isotropic turbulence in periodic domains of up to $12288^3$ grid points, and Taylor-scale Reynolds numbers in the range $140-1300$. We extract the statistics of various terms involved in amplification of strain and additionally condition them on the magnitude of strain. We find that strain is overall self-amplified by the quadratic nonlinearity, and depleted via vortex stretching; whereas pressure Hessian acts to redistribute strain fluctuations towards the mean-field and thus depleting intense strain. Analyzing the intense fluctuations of strain in terms of its eigenvalues reveals that the net amplification is solely produced by the third eigenvalue, resulting in strong compressive action. In contrast, the self-amplification terms acts to deplete the other two eigenvalues, whereas vortex stretching acts to amplify them, both effects canceling each other almost perfectly. The effect of the pressure Hessian for each eigenvalue is qualitatively similar to that of vortex stretching, but significantly weaker in magnitude. Our results conform with the familiar notion that intense strain is organized in sheet-like structures, which are in the vicinity of, but never overlap with regions of intense vorticity due to fundamental differences in their amplifying mechanisms.
An essential ingredient of turbulent flows is the vortex stretching mechanism, which emanates from the non-linear interaction of vorticity and strain-rate tensor and leads to formation of extreme events. We analyze the statistical correlations between vorticity and strain rate by using a massive database generated from very well resolved direct numerical simulations of forced isotropic turbulence in periodic domains. The grid resolution is up to $12288^3$, and the Taylor-scale Reynolds number is in the range $140-1300$. In order to understand the formation and structure of extreme vorticity fluctuations, we obtain statistics conditioned on enstrophy (vorticity-squared). The magnitude of strain, as well as its eigenvalues, is approximately constant when conditioned on weak enstrophy; whereas they grow approximately as power laws for strong enstrophy, which become steeper with increasing $R_lambda$. We find that the well-known preferential alignment between vorticity and the intermediate eigenvector of strain tensor is even stronger for large enstrophy, whereas vorticity shows a tendency to be weakly orthogonal to the most extensive eigenvector (for large enstrophy). Yet the dominant contribution to the production of large enstrophy events arises from the most extensive eigendirection, the more so as $R_lambda$ increases. Nevertheless, the stretching in intense vorticity regions is significantly depleted, consistent with the kinematic properties of weakly-curved tubes in which they are organized. Further analysis reveals that intense enstrophy is primarily depleted via viscous diffusion, though viscous dissipation is also significant. Implications for modeling are nominally addressed as appropriate.
Using exact relations between velocity structure functions (Hill, Hill and Boratav, and Yakhot) and neglecting pressure contributions in a first approximation, we obtain a closed system and derive simple order-dependent rescaling relationships between longitudinal and transverse structure functions. By means of numerical data with turbulent Reynolds numbers ranging from $Re_lambda=320$ to $Re_lambda=730$, we establish a clear correspondence between their respective scaling range, while confirming that their scaling exponents do differ. This difference does not seem to depend on Reynolds number. Making use of the Mellin transform, we further map longitudinal to (rescaled) transverse probability density functions.
Inertial particle data from three-dimensional direct numerical simulations of particle-laden homogeneous isotropic turbulence at high Reynolds number are analyzed using Voronoi tessellation of the particle positions, considering different Stokes numbers. A finite-time measure to quantify the divergence of the particle velocity by determining the volume change rate of the Voronoi cells is proposed. For inertial particles the probability distribution function (PDF) of the divergence deviates from that for fluid particles. Joint PDFs of the divergence and the Voronoi volume illustrate that the divergence is most prominent in cluster regions and less pronounced in void regions. For larger volumes the results show negative divergence values which represent cluster formation (i.e. particle convergence) and for small volumes the results show positive divergence values which represents cluster destruction/void formation (i.e. particle divergence). Moreover, when the Stokes number increases the divergence takes larger values, which gives some evidence why fine clusters are less observed for large Stokes numbers. Theoretical analyses further show that the divergence for random particles in random flow satisfies a PDF corresponding to the ratio of two independent variables following normal and gamma distributions in one dimension. Extending this model to three dimensions, the predicted PDF agrees reasonably well with Monte-Carlo simulations and DNS data of fluid particles.
Numerical simulations are made for forced turbulence at a sequence of increasing values of Reynolds number, R, keeping fixed a strongly stable, volume-mean density stratification. At smaller values of R, the turbulent velocity is mainly horizontal, and the momentum balance is approximately cyclostrophic and hydrostatic. This is a regime dominated by so-called pancake vortices, with only a weak excitation of internal gravity waves and large values of the local Richardson number, Ri, everywhere. At higher values of R there are successive transitions to (a) overturning motions with local reversals in the density stratification and small or negative values of Ri; (b) growth of a horizontally uniform vertical shear flow component; and (c) growth of a large-scale vertical flow component. Throughout these transitions, pancake vortices continue to dominate the large-scale part of the turbulence, and the gravity wave component remains weak except at small scales.
Multiscale statistical analyses of inertial particle distributions are presented to investigate the statistical signature of clustering and void regions in particle-laden incompressible isotropic turbulence. Three-dimensional direct numerical simulations of homogeneous isotropic turbulence at high Reynolds number ($Re_lambda gtrsim 200$) with up to $10^9$ inertial particles are performed for Stokes numbers ranging from $0.05$ to $5.0$. Orthogonal wavelet analysis is then applied to the computed particle number density fields. Scale-dependent skewness and flatness values of the particle number density distributions are calculated and the influence of Reynolds number $Re_lambda$ and Stokes number $St$ is assessed. For $St sim 1.0$, both the scale-dependent skewness and flatness values become larger as the scale decreases, suggesting intermittent clustering at small scales. For $St le 0.2$, the flatness at intermediate scales, i.e. for scales larger than the Kolmogorov scale and smaller than the integral scale of the flow, increases as $St$ increases, and the skewness exhibits negative values at the intermediate scales. The negative values of the skewness are attributed to void regions. These results indicate that void regions at the intermediate sales are pronounced and intermittently distributed for such small Stokes numbers. As $Re_lambda$ increases, the flatness increases slightly. For $Re_lambda ge 328$, the skewness shows negative values at large scales, suggesting that void regions are pronounced at large scales, while clusters are pronounced at small scales.