Structure of unsteady stably stratified turbulence with mean shear

Abstract

The statistics of unsteady turbulence with uniform stratification N (Brunt–Vaisala frequency) and shear α(=dU1/dx3) are analysed over the entire time range (0 0 and \it Ri>0.25 respectively, oscillatory momentum and positive and negative density fluxes develop. Above a critical value of \it Ri\scriptsize\it crit(∼0.3), their average values are persistently countergradient. This structural change in the turbulence is the primary mechanism whereby stable stratification reduces the fluxes and the production of variances. It is quite universal and differs from the energy and stability mechanisms of Richardson (1926) and Taylor (1931). The long-time asymptotics of the energy ratio ER(=\it PE/VKE) of the potential energy to the vertical kinetic energy generally decreases with \it Ri(≥0.25), reaching the smallest value of 3/2 when there is no shear (\it Ri→∞). For strong mean shear (\it Ri<0.25), RDT significantly overestimates ER since (as in unstratified shear flow) it underestimates the vertical kinetic energy VKE. The RDT results show that the asymptotic values of the energy ratio ER and the normalized vertical density flux are independent of the initial value of ER, in agreement with DNS. This independence of the initial condition occurs because the ratios of the contributions from the initial values PE0 and KE0 are the same for PE and VKE and can be explained by the linear processes. Stable stratification generates buoyancy oscillations in the direction of the energy propagation of the internal gravity wave and suppresses the generation of turbulence by mean shear. Because the shear distorts the wavenumber fluctuations, the low-wavenumber spectrum of the vertical kinetic energy has the general form E33(k)∝(αtk)−1, where (LXαt)−1≪k≪L−1X (LX: integral scale). The viscous decay is controlled by the shear, so that the components of larger streamwise wavenumber k1 decay faster. Then, combined with the spectrum distortion by the shear, the energy and the flux are increasingly dominated by the small-k1 components as time elapses. They oscillate at the buoyancy period π/N because even in a shear flow the components as k1→0 are weakly affected by the shear. The effects of stratification N and shear α at small scales are to reduce both VKE and PE. Even for the same \it Ri, larger N and α reduce the high-wavenumber components of VKE and PE. This supports the applicability of the linear assumption for large N and α. At large scales, the stratification and shear effects oppose each other, i.e. both VKE and PE decrease due to the stratification but they increase due to the shear. We conclude that certain of these unsteady results can be applied directly to estimate the properties of sheared turbulence in a statistically steady state, but others can only be applied qualitatively.