skip to main content

This content will become publicly available on May 18, 2023

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 more » the mixing enhancement that increases with lateral Froude number Fr L , i.e. stronger shear and/or steeper ridge. « less
; ; ; ; ;
Award ID(s):
Publication Date:
Journal Name:
Frontiers in Marine Science
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 ,more »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.« less
  2. 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,more »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.

    « less
  3. 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$ )more »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.« less
  4. 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).more »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.« less
  5. Abstract This study presents field observations of fluid mud and the flow instabilities that result from the interaction between mud-induced density stratification and current shear. Data collected by shipborne and bottom-mounted instruments in a hyperturbid estuarine tidal channel reveal the details of turbulent sheared layers in the fluid mud that persist throughout the tidal cycle. Shear instabilities form during periods of intense shear and strong mud-induced stratification, particularly with gradient Richardson number smaller than or fluctuating around the critical value of 0.25. Turbulent mixing plays a significant role in the vertical entrainment of fine sediment over the tidal cycle. The vertical extent of the billows identified seen in the acoustic images is the basis for two useful parameterizations. First, the aspect ratio (billow height/wavelength) is indicative of the initial Richardson number that characterizes the shear flow from which the billows grew. Second, we describe a scaling for the turbulent dissipation rate ε that holds for both observed and simulated Kelvin–Helmholtz billows. Estimates for the present observations imply, however, that billows growing on a lutocline obey an altered scaling whose origin remains to be explained.