In stably stratified turbulent flows, the mixing efficiency associated with eddy diffusivity for heat, or equivalently the turbulent Prandtl number (Pr(t)), is fraught with complex dynamics originating from the scalewise interplay between shear generation of turbulence and its dissipation by density gradients. A large corpus of data and numerical simulations agree on a near-universal relation between Pr(t) and the Richardson number (R(i)), which encodes the relative importance of buoyancy dissipation to mechanical production of turbulent kinetic energy. The Pr(t)-R(i) relation is shown to be derivable solely from the cospectral budgets for momentum and heat fluxes if a Rotta-like return to isotropy closure for the pressure-strain effects and Kolmogorov's theory for turbulent cascade are invoked. The ratio of the Kolmogorov to the Kolmogorov-Obukhov-Corrsin phenomenological constants, and a constant associated with isotropization of the production whose value (= 3/5) has been predicted from Rapid Distortion Theory, explain all the macroscopic nonlinearities.

Direct numerical simulations of stably stratified Ekman layers are conducted to study the effect of increasing static stability on turbulence dynamics and modelling in wall-bounded flows at three moderate Reynolds numbers. The flow field is analysed by examining the mean profiles of wind speed, potential temperature and momentum flux, as well as streamwise velocity and temperature spectra. The maximum stabilizing buoyancy flux that a flow can sustain while remaining fully turbulent is found to depend on the Reynolds number. The flows with the highest Reynolds number display a relatively well-developed inertial range and logarithmic layer, and are found to bear similarities to much higher-Reynolds-number flows like the ones encountered in the atmospheric boundary layer. In particular, the near-wall mean profiles follow the Monin-Obukhov similarity theory. However, several flow features, such as the critical Richardson number and the stress-strain alignment, are found to maintain significant dependence on the Reynolds number. The budgets of turbulence kinetic energy (TKE), vertical velocity variance, momentum and buoyancy fluxes, and temperature variance are analysed. The results indicate that the effect of stability on turbulence is first directly manifested in the vertical velocity variance budget, and results in damping of vertical motions. This then leads to a reduction in the downward transport of horizontal momentum components towards the surface, and consequently to a decrease in the shear production term in the TKE budget: changes in the vertical profile of TKE shear production with increasing Richardson number are significant and have a larger impact on TKE than direct buoyancy destruction. The reduction in vertical velocity variance also results in significant drops in the production terms in the momentum flux, buoyancy flux and temperature variance budgets. Various assumptions and parameters related to low-order turbulence closures are investigated. The results suggest that the vertical velocity variance is a more appropriate parameter than the full TKE on which to base eddy-diffusivity and viscosity models.

No abstract

No abstract

No abstract