An official website of the United States government
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
Shear Instabilities and Stratified Turbulence in an Estuarine Fluid Mud
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.
more »
« less
- Award ID(s):
- 1830071
- PAR ID:
- 10381491
- Date Published:
- Journal Name:
- Journal of Physical Oceanography
- Volume:
- 52
- Issue:
- 10
- ISSN:
- 0022-3670
- Page Range / eLocation ID:
- 2257 to 2271
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
We perform a direct numerical simulation (DNS) of interacting Kelvin–Helmholtz instabilities (KHI) that arise at a stratified shear layer where KH billow cores are misaligned or exhibit varying phases along their axes. Significant evidence of these dynamics in early laboratory shear-flow studies by Thorpe ( Geophys. Astrophys. Fluid Dyn. , vol. 34, 1985, pp. 175–199) and Thorpe ( J. Geophys. Res. , vol. 92, 1987, pp. 5231–5248), in observations of KH billow misalignments in tropospheric clouds (Thorpe, Q. J. R. Meteorol. Soc. , vol. 128, 2002, pp. 1529–1542) and in recent direct observations of such events in airglow and polar mesospheric cloud imaging in the upper mesosphere reveals that these dynamics are common. More importantly, the laboratory and mesospheric observations suggest that these dynamics lead to more rapid and more intense instabilities and turbulence than secondary convective instabilities in billow cores and secondary KHI in stratified braids between and around adjacent billows. To date, however, no simulations exploring the dynamics and energetics of interacting KH billows (apart from pairing) have been performed. Our DNS performed for Richardson number $Ri=0.10$ and Reynolds number $Re=5000$ demonstrates that KHI tubes and knots (i) comprise strong and complex vortex interactions accompanying misaligned KH billows, (ii) accelerate the transition to turbulence relative to secondary instabilities of individual KH billows, (iii) yield significantly stronger turbulence than secondary KHI in billow braids and secondary convective instabilities in KHI billow cores and (iv) expand the suite of secondary instabilities previously recognized to contribute to KHI dynamics and breakdown to turbulence in realistic geophysical environments.more » « less
-
We perform a direct numerical simulation (DNS) of interacting Kelvin–Helmholtz instabilities (KHI) that arise at a stratified shear layer where KH billow cores are misaligned or exhibit varying phases along their axes. Significant evidence of these dynamics in early laboratory shear-flow studies by Thorpe (Geophys. Astrophys. Fluid Dyn., vol. 34, 1985, pp. 175–199) and Thorpe (J. Geophys. Res., vol. 92, 1987, pp. 5231–5248), in observations of KH billow misalignments in tropospheric clouds (Thorpe, Q. J. R. Meteorol. Soc., vol. 128, 2002, pp. 1529–1542) and in recent direct observations of such events in airglow and polar mesospheric cloud imaging in the upper mesosphere reveals that these dynamics are common. More importantly, the laboratory and mesospheric observations suggest that these dynamics lead to more rapid and more intense instabilities and turbulence than secondary convective instabilities in billow cores and secondary KHI in stratified braids between and around adjacent billows. To date, however, no simulations exploring the dynamics and energetics of interacting KH billows (apart from pairing) have been performed. Our DNS performed for Richardson number Ri = 0.10 and Reynolds number Re = 5000 demonstrates that KHI tubes and knots (i) comprise strong and complex vortex interactions accompanying misaligned KH billows, (ii) accelerate the transition to turbulence relative to secondary instabilities of individual KH billows, (iii) yield significantly stronger turbulence than secondary KHI in billow braids and secondary convective instabilities in KHI billow cores and (iv) expand the suite of secondary instabilities previously recognized to contribute to KHI dynamics and breakdown to turbulence in realistic geophysical environments.more » « less
-
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
-
Abstract Flow over a surface can be stratified by imposing a fixed mean vertical temperature (density) gradient profile throughout or via cooling at the surface. These distinct mechanisms can act simultaneously to establish a stable stratification in a flow. Here, we perform a series of direct numerical simulations of open-channel flows to study adaptation of a neutrally stratified turbulent flow under the combined or independent action of the aforementioned mechanisms. We force the fully developed flow with a constant mass flow rate. This flow forcing technique enables us to keep the bulk Reynolds number constant throughout our investigation and avoid complications arising from the acceleration of the bulk flow if a constant pressure gradient approach were to be adopted to force the flow instead. When both stratification mechanisms are active, the dimensionless stratification perturbation number emerges as an external flow control parameter, in addition to the Reynolds, Froude, and Prandtl numbers. We demonstrate that significant deviations from the Monin–Obukhov similarity formulation are possible when both types of stratification mechanisms are active within an otherwise weakly stable flow, even when the flux Richardson number is well below 0.2. An extended version of the similarity theory due to Zilitinkevich and Calanca shows promise in predicting the dimensionless shear for cases where both types of stratification mechanisms are active, but the extended theory is less accurate for gradients of scalar. The degree of deviation from neutral dimensionless shear as a function of the vertical coordinate emerges as a qualitative measure of the strength of stable stratification for all the cases investigated in this study.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1175/JPO-D-21-0230.1
