Abstract Direct numerical simulations are performed to compare the evolution of turbulent stratified shear layers with different density gradient profiles at a high Reynolds number. The density profiles include uniform stratification, two-layer hyperbolic tangent profile and a composite of these two profiles. All profiles have the same initial bulk Richardson number ( $$Ri_{b,0}$$ R i b , 0 ); however, the minimum gradient Richardson number and the distribution of density gradient across the shear layer are varied among the cases. The objective of the study is to provide a comparative analysis of the evolution of the shear layers in term of shear layer growth, turbulent kinetic energy as well as the mixing efficiency and its parameterization. The evolution of the shear layers in all cases shows the development of Kelvin–Helmholtz billows, the transition to turbulence by secondary instabilities followed by the decay of turbulence. Comparison among the cases reveals that the amount of turbulent mixing varies with the density gradient distribution inside the shear layer. The minimum gradient Richardson number and the initial bulk Richardson number do not correlate well with the integrated TKE production, dissipation and buoyancy flux. The bulk mixing efficiency for fixed $$Ri_{b,0}$$ R i b , 0 is found to be highest in the case with two-layer density profile and lowest in the case with uniform stratification. However, the parameterizations of the flux coefficient based on buoyancy Reynolds number and the ratio of Ozmidov and Ellison scales show similar scaling in all cases.
more »
« less
Turbulent shear layers in a uniformly stratified background: DNS at high Reynolds number
Direct numerical simulations are performed to investigate a stratified shear layer at high Reynolds number ( $Re$ ) in a study where the Richardson number ( $Ri$ ) is varied among cases. Unlike previous work on a two-layer configuration in which the shear layer resides between two layers with constant density, an unbounded fluid with uniform stratification is considered here. The evolution of the shear layer includes a primary Kelvin–Helmholtz shear instability followed by a wide range of secondary shear and convective instabilities, similar to the two-layer configuration. During transition to turbulence, the shear layers at low $Ri$ exhibit a period of thickness contraction (not observed at lower $Re$ ) when the momentum and buoyancy fluxes are counter-gradient. The behaviour in the turbulent regime is significantly different from the case with a two-layer density profile. The transition layers, which are zones with elevated shear and stratification that form at the shear-layer edges, are stronger and also able to support a significant internal wave flux. After the shear layer becomes turbulent, mixing in the transition layers is shown to be more efficient than that which develops in the centre of the shear layer. Overall, the cumulative mixing efficiency ( $E^C$ ) is larger than the often assumed value of 1/6. Also, $E^C$ is found to be smaller than that in the two-layer configuration at moderate Ri . It is relatively less sensitive to background stratification, exhibiting little variation for $0.08 \leqslant Ri \leqslant 0.2$ . The dependence of mixing efficiency on buoyancy Reynolds number during the turbulence phase is qualitatively similar to homogeneous sheared turbulence.
more »
« less
- Award ID(s):
- 1851390
- NSF-PAR ID:
- 10339231
- Date Published:
- Journal Name:
- Journal of Fluid Mechanics
- Volume:
- 916
- ISSN:
- 0022-1120
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
null (Ed.)Turbulence parameters in the lower troposphere (up to ~4.5 km) are estimated from measurements of high-resolution and fast-response cold-wire temperature and Pitot tube velocity from sensors onboard DataHawk Unmanned Aerial Vehicles (UAVs) operated at the Shigaraki Middle and Upper atmosphere (MU) Observatory during two ShUREX (Shigaraki UAV Radar Experiment) campaigns in 2016 and 2017. The practical processing methods used for estimating turbulence kinetic energy dissipation rate ε and temperature structure function parameter C T 2 from one-dimensional wind and temperature frequency spectra are first described in detail. Both are based on the identification of inertial (−5/3) subranges in respective spectra. Using a formulation relating ε and C T 2 valid for Kolmogorov turbulence in steady state, the flux Richardson number R f and the mixing efficiency χ m are then estimated. The statistical analysis confirms the variability of R f and χ m around ~ 0.13 − 0.14 and ~ 0.16 − 0.17 , respectively, values close to the canonical values found from some earlier experimental and theoretical studies of both the atmosphere and the oceans. The relevance of the interpretation of the inertial subranges in terms of Kolmogorov turbulence is confirmed by assessing the consistency of additional parameters, the Ozmidov length scale L O , the buoyancy Reynolds number R e b , and the gradient Richardson number Ri. Finally, a case study is presented showing altitude differences between the peaks of N 2 , C T 2 and ε , suggesting turbulent stirring at the margin of a stable temperature gradient sheet. The possible contribution of this sheet and layer structure on clear air radar backscattering mechanisms is examined.more » « less
-
more » « less
Abstract The flux Richardson number
R f , also called the mixing efficiency of stratified turbulence, is important in determining geophysical flow phenomena such as ocean circulation and air‐sea transports. MeasuringR f in the field is usually difficult, thus parameterization ofR f based on readily observed properties is essential. Here, estimates ofR f in a strongly turbulent, sediment‐stratified estuarine flow are obtained from measurements of covariance‐derived turbulent buoyancy fluxes (B ) and spectrally fitted values of the dissipation rate of turbulent kinetic energy (ε ). We test scalings forR f in terms of the buoyancy Reynolds number (Re b ), the gradient Richardson number (Ri ), and turbulent Froude number (Fr t ). Neither theRe b ‐based nor theRi ‐based scheme is able to describe the observed variations inR f , but theFr t ‐based parameterization works well. These findings support further use of theFr t ‐ based parameterization in turbulent oceanic and estuarine environments. -
null (Ed.)Ekman layers develop at the boundaries of the Earth’s fluid core in response to precession. Instabilities in these layers lead to turbulence when a local Reynolds number, Re, based on the thickness of the Ekman layer, exceeds a critical value. The transition to turbulence is often assessed using experiments for steady Ekman layers, where the interior geostrophic flow is independent of time. Precessionally driven flow varies on diurnal timescales, so the transition to turbulence may occur at a different value of Re.We use 3-D numerical calculations in a local Cartesian geometry to assess the transition to turbulence in precessional flow. Calculations retain the horizontal component of the rotation vector and account for the influence of fluid stratification. The transition to turbulence in a neutrally stratified fluid occurs near Re = 500, which is higher than the value Re = 150 usually cited for steady Ekman layers. However, it is comparable to the nominal value for precessional flow in the Earth. Complications due to fluid stratification or a magnetic field can suppress the transition to turbulence, reducing the likelihood of turbulent Ekman layers in the Earth’s core.more » « less
-
Abstract Upper-ocean turbulence is central to the exchanges of heat, momentum, and gasses across the air/sea interface, and therefore plays a large role in weather and climate. Current understanding of upper-ocean mixing is lacking, often leading models to misrepresent mixed-layer depths and sea surface temperature. In part, progress has been limited due to the difficulty of measuring turbulence from fixed moorings which can simultaneously measure surface fluxes and upper-ocean stratification over long time periods. Here we introduce a direct wavenumber method for measuring Turbulent Kinetic Energy (TKE) dissipation rates, ϵ , from long-enduring moorings using pulse-coherent ADCPs. We discuss optimal programming of the ADCPs, a robust mechanical design for use on a mooring to maximize data return, and data processing techniques including phase-ambiguity unwrapping, spectral analysis, and a correction for instrument response. The method was used in the Salinity Processes Upper-ocean Regional Study (SPURS) to collect two year-long data sets. We find the mooring-derived TKE dissipation rates compare favorably to estimates made nearby from a microstructure shear probe mounted to a glider during its two separate two-week missions for (10 −8 ) ≤ ϵ ≤ (10 −5 ) m 2 s −3 . Periods of disagreement between turbulence estimates from the two platforms coincide with differences in vertical temperature profiles, which may indicate that barrier layers can substantially modulate upper-ocean turbulence over horizontal scales of 1-10 km. We also find that dissipation estimates from two different moorings at 12.5 m, and at 7 m are in agreement with the surface buoyancy flux during periods of strong nighttime convection, consistent with classic boundary layer theory.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1017/jfm.2021.212
