Small-scale statistics in high-resolution direct numerical simulation of turbulence: Reynolds number dependence of one-point velocity gradient statistics

Abstract

One-point statistics of velocity gradients and Eulerian and Lagrangian accelerations are studied by analysing the data from high-resolution direct numerical simulations (DNS) of turbulence in a periodic box, with up to 40963 grid points. The DNS consist of two series of runs; one is with kmaxη ~ 1 (Series 1) and the other is with kmaxη ~ 2 (Series 2), where kmax is the maximum wavenumber and η the Kolmogorov length scale. The maximum Taylor-microscale Reynolds number Rλ in Series 1 is about 1130, and it is about 675 in Series 2. Particular attention is paid to the possible Reynolds number (Re) dependence of the statistics. The visualization of the intense vorticity regions shows that the turbulence field at high Re consists of clusters of small intense vorticity regions, and their structure is to be distinguished from those of small eddies. The possible dependence on Re of the probability distribution functions of velocity gradients is analysed through the dependence on Rλ of the skewness and flatness factors (S and F). The DNS data suggest that the Rλ dependence of S and F of the longitudinal velocity gradients fit well with a simple power law: S ~ −0.32Rλ0.11 and F ~ 1.14Rλ0.34, in fairly good agreement with previous experimental data. They also suggest that all the fourth-order moments of velocity gradients scale with Rλ similarly to each other at Rλ > 100, in contrast to Rλ < 100. Regarding the statistics of time derivatives, the second-order time derivatives of turbulent velocities are more intermittent than the first-order ones for both the Eulerian and Lagrangian velocities, and the Lagrangian time derivatives of turbulent velocities are more intermittent than the Eulerian time derivatives, as would be expected. The flatness factor of the Lagrangian acceleration is as large as 90 at Rλ ≈ 430. The flatness factors of the Eulerian and Lagrangian accelerations increase with Rλ approximately proportional to RλαE and RλαL, respectively, where αE ≈ 0.5 and αL ≈ 1.0, while those of the second-order time derivatives of the Eulerian and Lagrangian velocities increases approximately proportional to RλβE and RλβL, respectively, where βE ≈ 1.5 and βL ≈ 3.0.