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

A cascade hypothesis for a strongly stratified fluid is developed on the basis of the Boussinesq equations. According to this hypothesis, kinetic and potential energy are transferred from large to small scales in a highly anisotropic turbulent cascade. A relation for the ratio, $ l_{v}/l_{h} $, between the vertical and horizontal length scale is derived, showing how this ratio decreases with increased stratification. Similarity expressions are formulated for the horizontal and vertical spectra of kinetic and potential energy. A series of box simulations of the Boussinesq equations are carried out and a good agreement between the proposed hypothesis and the simulations is seen. The simulations with strongest stratification give horizontal kinetic and potential energy spectra of the form $ E_{K_{h}} \,{=}\, C_{1} \epsilon_{K}^{2/3} k_{h}^{-5/3} $ and $ E_{P_{h}} \,{=}\, C_{2} \epsilon_{P} k_{h}^{-5/3}/\epsilon_{K}^{1/3} $, where $ k_{h} $ is the horizontal wavenumber, $ \epsilon_{K} $ and $ \epsilon_{P} $ are the dissipation of kinetic and potential energy, respectively, and $ C_{1} $ and $ C_{2} $ are two constants. Within the given numerical accuracy, it is found that these two constants have the same value: $ C_{1} \approx C_{2} \,{=}\, 0.51 \pm 0.02 $.

The statistical features of turbulence can be studied either through spectral quantities, such as the kinetic energy spectrum, or through structure functions, which are statistical moments of the difference between velocities at two points separated by a variable distance. In this paper structure function relations for two-dimensional turbulence are derived and compared with calculations based on wind data from 5754 airplane flights, reported in the MOZAIC data set. For the third-order structure function two relations are derived, showing that this function is generally positive in the two-dimensional case, contrary to the three-dimensional case. In the energy inertial range the third-order structure function grows linearly with separation distance and in the enstrophy inertial range it grows cubically with separation distance. A Fourier analysis shows that the linear growth is a reflection of a constant negative spectral energy flux, and the cubic growth is a reflection of a constant positive spectral enstrophy flux. Various relations between second-order structure functions and spectral quantities are also derived. The measured second-order structure functions can be divided into two different types of terms, one of the form r2/3, giving a k−5/3-range and another, including a logarithmic dependence, giving a k−3-range in the energy spectrum. The structure functions agree better with the two-dimensional isotropic relation for larger separations than for smaller separations. The flatness factor is found to grow very fast for separations of the order of some kilometres. The third-order structure function is accurately measured in the interval [30, 300] km and is found to be positive. The average enstrophy flux is measured as Πω≈1.8×10−13 s−3 and the constant in the k−3-law is measured as [Kscr ]≈0.19. It is argued that the k−3-range can be explained by two-dimensional turbulence and can be interpreted as an enstrophy inertial range, while the k−5/3-range can probably not be explained by two-dimensional turbulence and should not be interpreted as a two-dimensional energy inertial range.

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.

Several existing sets of smaller-scale ocean and atmospheric data appear to display Kolmogorov–Obukov–Corrsin inertial ranges in horizontal spectra for length scales up to at least a few hundred meters. It is argued here that these data are inconsistent with the assumptions for these inertial range theories. Instead, it is hypothesized that the dynamics of stratified turbulence explain these data. If valid, these dynamics may also explain the behavior of strongly stratified flows in similar dynamic ranges of other geophysical flows.