The flux Richardson number
- Award ID(s):
- 1830071
- PAR ID:
- 10370191
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- Geophysical Research Letters
- Volume:
- 49
- Issue:
- 13
- ISSN:
- 0094-8276
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
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
-
Abstract Kelvin‐Helmholtz (KH) instability plays an important role in turbulent mixing in deep oceans, coastal seas, and estuaries. Though widely observed and studied in thermohaline‐stratified waters, KH instability has only rarely been observed in sediment‐stratified environments. For the first time, we present direct observations of KH billows on an estuarine lutocline by combining echosounder images with velocity and density measurements. The interaction between velocity shear and the density stratification induced by suspended sediments initiated shear instabilities near the bed, indicated by gradient Richardson number (Ri) < 0.25 in the early stages of the observed billows. Once formed, the instabilities enhanced the vertical mixing of momentum, reducing vertical shear and elevating Ri. Linear instability analysis using measured velocity and density profiles well predicts the vertical location and spatial characteristics of the observed billows. These instabilities are believed to contribute to the vertical mixing, entrainment, and transport of estuarine and coastal sediments.
-
Abstract Anisotropic turbulence is ubiquitous in atmospheric and oceanic boundary layers due to differences in energy injection mechanisms. Unlike mechanical production that injects energy in the streamwise velocity component, buoyancy affects only the vertical velocity component. This anisotropy in energy sources, quantified by the flux Richardson number
, is compensated by a “return to isotropy” (RTI) tendency of turbulent flows. Describing RTI in Reynolds‐averaged models and across scales continues to be a challenge in stratified turbulent flows. Using phenomenological models for spectral energy transfers, the necessary conditions for which the widely‐used Rotta model captures RTI across variousR i f and eddy sizes are discussed for the first time. This work unravels adjustments to the Rotta constant, withR i f and scale, necessary to obtain consistency between RTI models and the measured properties of the atmospheric surface layer for planar‐homogeneous and stationary flows in the absence of subsidence. A range ofR i f and eddy sizes where the usage of a conventional Rotta model is prohibited is also found. Those adjustments lay the groundwork for new closure schemes.R i f -
Abstract We formulate an expression for the turbulent kinetic energy dissipation rate,
ϵ , associated with shear‐generated turbulence in terms of quantities in the ocean or atmosphere that, depending on the situation, may be measurable or resolved in models. The expression depends on the turbulent vertical length scale,ℓ v , the inverse time scaleN , and the Richardson numberRi =N 2/S 2, whereS is the vertical shear, withℓ v scaled in a way consistent with theories and observations of stratified turbulence. Unlike previous studies, the focus is not so much on the functional form ofRi , but the vertical variation of the length scaleℓ v . Using data from two ∼7‐day time series in the western equatorial Pacific, the scaling is compared with the observedϵ . The scaling works well with the estimatedϵ capturing the differences in amplitude and vertical distribution of the observedϵ between the two times series. Much of those differences are attributable to changes in the vertical distribution of the length scaleℓ v , and in particular the associated turbulent velocity scale,u t . We relateu t to a measure of the fine‐scale variations in velocity,. Our study highlights the need to consider the length scale and its estimation in environmental flows. The implications for the vertical variation of the associated turbulent diffusivity are discussed. -
Abstract We performed laboratory experiments to investigate the influence of sand content on the dynamics of wave‐supported gravity flows in mud‐dominant environments. The experiments were carried out in an oscillatory water tunnel with a sediment bed of either 1% or 13% sand. Low and high energy regimes are differentiated based on a Stokes Reynolds number
Re Δ ≈ 500. In the low energy regime, the sand fraction influences flow dynamics primarily through ripple formation; no ripples form in the 1% sand experiments, whereas ripples form in the 13% experiments that increase turbulence and the wave boundary layer thickness,δ m . In the high energy regime, small ripples form in both the 1% and 13% sand experiments and we observe high near‐bed suspended sediment concentrations. The influence of stratification on the boundary layer flow is characterized in terms of the gradient Richardson numberRi g . The flow is weakly stratified inside the boundary layer for all runs and critically stratified at or above the top of the boundary layer. In the lower energy regime, the sand content reduces the relative influence of stratification in the boundary layer, shifting the elevation of critical stratification,L B , from approximately 1.3δ m to 2.5δ m in the 1% and 13% experiments, respectively. In both sets of experimentsL B ≈δ m at the strongest wave energy, indicating a transition to strongly stratified dynamics.