skip to main content


Title: Shear Instability and Turbulent Mixing in the Stratified Shear Flow Behind a Topographic Ridge at High Reynolds Number
Observations on the lee of a topographic ridge show that the turbulence kinetic energy (TKE) dissipation rate due to shear instabilities is three orders of magnitude higher than the typical value in the open ocean. Laboratory-scale studies at low Reynolds number suggest that high turbulent dissipation occurs primarily within the core region of shear instabilities. However, field-scale studies indicate that high turbulence is mainly populated along the braids of shear instabilities. In this study, a high-resolution, resolving the Ozmidov-scale, non-hydrostatic model with Large Eddy Simulation (LES) turbulent closure is applied to investigate dominant mechanisms that control the spatial and temporal scales of shear instabilities and resulting mixing in stratified shear flow at high Reynolds number. The simulated density variance dissipation rate is elevated in the cusp-like bands of shear instabilities with a specific period, consistent with the acoustic backscatter taken by shipboard echo sounder. The vertical length scale of each cusp-like band is nearly half of the vertical length scale of the internal lee wave. However, it is consistent with instabilities originating from a shear layer based on linear stability theory. The model results indicate that the length scale and/or the period of shear instabilities are the key parameters to the mixing enhancement that increases with lateral Froude number Fr L , i.e. stronger shear and/or steeper ridge.  more » « less
Award ID(s):
2128895
NSF-PAR ID:
10355953
Author(s) / Creator(s):
; ; ; ; ;
Date Published:
Journal Name:
Frontiers in Marine Science
Volume:
9
ISSN:
2296-7745
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. 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
  2. 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 = N2/S2, whereSis the vertical shear, withvscaled 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 scalev. 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 scalev, and in particular the associated turbulent velocity scale,ut. We relateutto 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.

     
    more » « less
  3. Abstract

    Thorpe analysis has been used to study turbulence in the atmosphere and ocean. It is clear that Thorpe analysis applied to individual soundings cannot be expected to give quantitatively reliable measurements of turbulence parameters because of the instantaneous nature of the measurement. A critical aspect of this analysis is the assumption of the linear relationship C = LO/LT between the Thorpe scale LT, derived from the sounding measurements, and the Ozmidov scale LO. It is the determination of LO that enables determination of the dissipation rate of turbulence kinetic energy ε. Single atmospheric and oceanic soundings cannot indicate either the source of turbulence or the stage of its evolution; different values of C are expected for different turbulence sources and stages of the turbulence evolution and thus cannot be expected to yield quantitatively reliable turbulence parameters from individual profiles. The variation of C with the stage of turbulence evolution is illustrated for direct numerical simulation (DNS) results for gravity wave breaking. Results from a DNS model of multiscale initiation and evolution of turbulence with a Reynolds number Re (which is defined using the vertical wavelength of the primary gravity wave and background buoyancy period as length and time scales, respectively) of 100 000 are sampled as in sounding of the atmosphere and ocean, and various averaging of the sounding results indicates a convergence to a well-defined value of C, indicating that applying Thorpe analysis to atmospheric or oceanic soundings and averaging over a number of profiles gives more reliable turbulence determinations. The same averaging study is also carried out when the DNS-modeled turbulence is dominated by turbulence growing from the initial instabilities, when the turbulence is fully developed, when the modeled turbulence is decaying, and when the turbulence is in a still-later decaying stage. These individual cases converge to well defined values of C, but these values of C show a large variation resulting from the different stages of turbulence evolution. This study gives guidance as to the accuracy of Thorpe analysis of turbulence as a function of the number of profiles being averaged. It also suggests that the values of C in different environments likely depend on the dominant turbulence initiation mechanisms and on the Reynolds number of the environment.

     
    more » « less
  4. 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
  5. null (Ed.)
    Abstract Mangrove swamps are extremely productive ecosystems providing many ecological services in coastal regions. The hydrodynamic interactions of mangrove roots and water flow have been proposed as a key element to mitigate erosion. Several studies reveal that precise prediction of the morphological evolution of coastal areas, in the face of global warming and the consequent sea-level rise, requires an understanding of interactions between root porosity (the fraction of the volume of void space over the total volume), water flows, and sediment transport. Water flows around the mangrove prop roots create a complex energetic process that mixes up sediments and generates a depositional region posterior to the roots. In this work, we investigated the boundary layer behind permeable arrays of cylinders (patch) that represent the mangrove roots to explore the impact of patch porosity on the onset of sediment transport. The flow measurements were performed in a vertical plane along the water depth downstream of the mangrove root models. A high-resolution Particle Image Velocimetry (PIV) was used in a flume to observe the impact of porosity on the mean flow, velocity derivatives, skin friction coefficient, and production of turbulent kinetic energy for Reynolds number of 2500 (based on patch diameter length-scale). Here, we proposed a predictive model for critical velocity for incipient motion that takes into account the mangrove roots porosity and the near-bed turbulence effect. It is found that the patch with the $$\phi =47\%$$ ϕ = 47 % porosity, has the maximum critical velocity over which the sediment transport initiates. We found the optimum porosity has the minimum sediment erosion and creates negative vorticity sources near the bed that increases the critical velocity. This signifies an optimum porosity for the onset of sediment transport consistent with the porosity of mangroves in nature. The phenomenological model is elucidated based on an analysis of the vorticity evolution equation for viscous incompressible flows. For the optimum porous patch, a sink of vorticity was formed which yielded to lower the near-bed turbulence and vorticity. The minimum velocity fluctuations were sufficient to initiate the boundary layer transition, however, the viscous dissipation dominated the turbulence production to obstruct the sediment transport. This work identified the pivotal role of mangrove root porosity in sediment transport in terms of velocity and its derivatives in wall-bounded flows. Our work also provides insight into the sediment transport and erosion processes that govern the evolution of the shapes of shorelines. 
    more » « less