International audienceDirect numerical simulations of stably and strongly stratified turbulent flows with Reynolds number Re " 1 and horizontal Froude number Fh Gt; 1 are presented. The results are interpreted on the basis of a scaling analysis of the governing equations. The analysis suggests that there are two different strongly stratified regimes according to the parameter R = ReFh2. When R " 1, viscous forces are nimportant and lv scales as lv ~ U/N (U is a characteristic horizontal velocity and N is the Brunt - Väis¨alä frequency) so that the dynamics of the flow is inherently three-dimensional but strongly anisotropic. When R " 1, vertical viscous shearing is important so that lv ~ lh/Re1/2 (lh is a characteristic horizontal length scale). The parameter R is further shown to be related to the buoyancy Reynolds number and proportional to (lO/?) 4/3, where lO is the Ozmidov length scale and ? the Kolmogorov length scale. This implies that there are simultaneously two distinct ranges in strongly stratified turbulence when R " 1: the scales larger than lO are strongly influenced by the stratification while those between lO and ? are weakly affected by stratification. The direct numerical simulations with forced large-scale horizontal two-dimensional motions and uniform stratification cover a wide Re and Fh range and support the main parameter controlling strongly stratified turbulence being R. The numerical results are in good agreement with the scaling laws for the vertical length scale. Thin horizontal layers are observed independently of the value of R but they tend to be smooth for R > 1, while for R > 1 small-scale three-dimensional turbulent disturbances are increasingly superimposed. The dissipation of kinetic energy is mostly due to vertical shearing for R > 1 but tends to isotropy as R increases above unity. When R > 1, the horizontal and vertical energy spectra are very steep while, when R > 1, the horizontal spectra of kinetic and potential energy exhibit an pproximate kh-5/3-power-law range and a clear forward energy cascade is observed. © 2007 Cambridge University Press

Direct numerical simulations ͑DNSs͒ and experiments of a spatially developing zeropressure-gradient turbulent boundary layer are presented up to Reynolds number Re = 2500, based on momentum thickness and free-stream velocity. For the first time direct comparisons of DNS and experiments of turbulent boundary layers at the same ͑computationally high and experimentally low͒ Re are given, showing excellent agreement in skin friction, mean velocity, and turbulent fluctuations. These results allow for a substantial reduction of the uncertainty of boundary-layer data, and cross validate the numerical setup and experimental technique. The additional insight into the flow provided by DNS clearly shows large-scale turbulent structures, which scale in outer units growing with Re , spanning the whole boundary-layer height.

No abstract

In this paper we consider the direct numerical simulation (DNS) of a spatially developing free round jet at low Reynolds numbers. Simulation of a spatially evolving flow such as the jet requires boundary conditions, which allow entrainment into the turbulent flow across the lateral boundaries of the computational domain. The boundary conditions which satisfy this requirement are so-called traction free boundary conditions. After showing that these boundary conditions lead to a correct behavior of the velocity near the lateral boundary of the jet, we will consider the DNS of the jet flow at a Reynolds number of 2.4×103 and compare the results with experimental data obtained by Hussein et al. [J. Fluid Mech. 258, 31 (1994)] and by Panchapakesan and Lumley [J. Fluid Mech. 246, 197 (1993)]. The results of our numerical simulations agree very well with the experimental data. Next we use the DNS to investigate the influence of the shape of the velocity profile at the jet orifice on the self-similarity scaling for the far-field velocity and shear stress profile. Evidence is presented in support of the suggestion by George [Advances in Turbulence (Springer, New York, 1989)] that the details of self-similarity depend on the initial conditions. This fact implies that there may exist no universally valid similarity scaling for the free jet.

Rare negative streamwise velocities and extreme wall-normal velocity fluctuations near the wall are investigated for turbulent channel flow at a series of Reynolds numbers based on friction velocity up to Re τ = 1000. Probability density functions of the wall-shear stress and velocity components are presented as well as joint probability density functions of the velocity components and the pressure. Backflow occurs more often (0.06% at the wall at Re τ = 1000) and further away (up to y + = 8.5) from the wall for increasing Reynolds number. The regions of backflow are circular with an average diameter, based on ensemble averages, of approximately 20 viscous units independent of Reynolds number. A strong oblique vortex outside the viscous sublayer is found to cause this backflow. Extreme wall-normal velocity events occur also more often for increasing Reynolds number. These extreme fluctuations cause high flatness values near the wall (F(v) = 43 at Re τ = 1000). Positive and negative velocity spikes appear in pairs, located on the two edges of a strong streamwise vortex as documented by Xu et al. [Phys. Fluids 8, 1938 (1996)] for Re τ = 180. The spikes are elliptical and orientated in streamwise direction with a typical length of 25 and a typical width of 7.5 viscous units at y + ≈ 1. The negative spike occurs in a high-speed streak indicating a sweeping motion, while the positive spike is located in between a high and low-speed streak. The joint probability density functions of negative streamwise and extreme wall-normal velocity events show that these events are largely uncorrelated. The majority of both type of events can be found lying underneath a large-scale structure in the outer region with positive sign, which can be understood by considering the more intense velocity fluctuations due to amplitude modulation of the inner layer by the outer layer. Simulations performed at different resolutions give only minor differences. Results from experiments and recent turbulent boundary layer simulations show similar results indicating that these rare events are universal for wall-bounded flows. In order to detect these rare events in experiments, measurement techniques have to be specifically tuned.

We consider mixing of the density field in stratified turbulence and argue that, at sufficiently high Reynolds numbers, stationary turbulence will have a mixing efficiency and closely related mixing coefficient described solely by the turbulent Froude number$Fr={\it\epsilon}_{k}/(Nu^{2})$, where${\it\epsilon}_{k}$is the kinetic energy dissipation,$u$is a turbulent horizontal velocity scale and$N$is the Brunt–Väisälä frequency. For$Fr\gg 1$, in the limit of weakly stratified turbulence, we show through a simple scaling analysis that the mixing coefficient scales as${\it\Gamma}\propto Fr^{-2}$, where${\it\Gamma}={\it\epsilon}_{p}/{\it\epsilon}_{k}$and${\it\epsilon}_{p}$is the potential energy dissipation. In the opposite limit of strongly stratified turbulence with$Fr\ll 1$, we argue that${\it\Gamma}$should reach a constant value of order unity. We carry out direct numerical simulations of forced stratified turbulence across a range of$Fr$and confirm that at high$Fr$,${\it\Gamma}\propto Fr^{-2}$, while at low$Fr$it approaches a constant value close to${\it\Gamma}=0.33$. The parametrization of${\it\Gamma}$based on$Re_{b}$due to Shihet al.(J. Fluid Mech., vol. 525, 2005, pp. 193–214) can be reinterpreted in this light because the observed variation of${\it\Gamma}$in their study as well as in datasets from recent oceanic and atmospheric measurements occurs at a Froude number of order unity, close to the transition value$Fr=0.3$found in our simulations.

Direct numerical simulations of subcritical rotating, stratified and magneto-hydrodynamic wall-bounded flows are performed in large computational domains, focusing on parameters where laminar and turbulent flow can stably coexist. In most cases, a regime of large-scale oblique laminar-turbulent patterns is identified at the onset of transition, as in the case of pure shear flows. The current study indicates that this oblique regime can be shifted up to large values of the Reynolds number $\mathit{Re}$ by increasing the damping by the Coriolis, buoyancy or Lorentz force. We show evidence for this phenomenon in three distinct flow cases: plane Couette flow with spanwise cyclonic rotation, plane magnetohydrodynamic channel flow with a spanwise or wall-normal magnetic field, and open channel flow under stable stratification. Near-wall turbulence structures inside the turbulent patterns are invariably found to scale in terms of viscous wall units as in the fully turbulent case, while the patterns themselves remain large-scale with a trend towards shorter wavelength for increasing $\mathit{Re}$. Two distinct regimes are identified: at low Reynolds numbers the patterns extend from one wall to the other, while at large Reynolds number they are confined to the near-wall regions and the patterns on both channel sides are uncorrelated, the core of the flow being highly turbulent without any dominant large-scale structure.

No abstract