An official website of the United States government
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
A comparative study of turbulent stratified shear layers: effect of density gradient distribution
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
- Award ID(s):
- 1851390
- PAR ID:
- 10356816
- Date Published:
- Journal Name:
- Environmental Fluid Mechanics
- ISSN:
- 1567-7419
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
High Reynolds number wall-bounded turbulent flows subject to buoyancy forces are fraught with complex dynamics originating from the interplay between shear generation of turbulence ( $$S$$ ) and its production or destruction by density gradients ( $$B$$ ). For horizontal walls, $$S$$ augments the energy budget of the streamwise fluctuations, while $$B$$ influences the energy contained in the vertical fluctuations. Yet, return to isotropy remains a tendency of such flows where pressure–strain interaction redistributes turbulent energy among all three velocity components and thus limits, but cannot fully eliminate, the anisotropy of the velocity fluctuations. A reduced model of this energy redistribution in the inertial (logarithmic) sublayer, with no tuneable constants, is introduced and tested against large eddy and direct numerical simulations under both stable ( $B<0$ ) and unstable ( $B>0$ ) conditions. The model links key transitions in turbulence statistics with flux Richardson number (at $$Ri_{f}=-B/S\approx$$ $-2$ , $-1$ and $-0.5$ ) to shifts in the direction of energy redistribution. Furthermore, when coupled to a linear Rotta-type closure, an extended version of the model can predict individual variance components, as well as the degree of turbulence anisotropy. The extended model indicates a regime transition under stable conditions when $$Ri_{f}$$ approaches $$Ri_{f,max}\approx +0.21$$ . Buoyant destruction $$B$$ increases with increasing stabilizing density gradients when $$Ri_{f}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
-
null (Ed.)Abstract Large-eddy simulations are used to investigate the influence of a horizontal frontal zone, represented by a stationary uniform background horizontal temperature gradient, on the wind- and wave-driven ocean surface boundary layers. In a frontal zone, the temperature structure, the ageostrophic mean horizontal current, and the turbulence in the ocean surface boundary layer all change with the relative angle among the wind and the front. The net heating and cooling of the boundary layer could be explained by the depth-integrated horizontal advective buoyancy flux, called the Ekman Buoyancy Flux (or the Ekman-Stokes Buoyancy Flux if wave effects are included). However, the detailed temperature profiles are also modulated by the depth-dependent advective buoyancy flux and submesoscale eddies. The surface current is deflected less (more) to the right of the wind and wave when the depth-integrated advective buoyancy flux cools (warms) the ocean surface boundary layer. Horizontal mixing is greatly enhanced by submesoscale eddies. The eddy-induced horizontal mixing is anisotropic and is stronger to the right of the wind direction. Vertical turbulent mixing depends on the superposition of the geostrophic and ageostrophic current, the depth-dependent advective buoyancy flux, and submesoscale eddies.more » « less
-
Abstract The flux Richardson numberRf, also called the mixing efficiency of stratified turbulence, is important in determining geophysical flow phenomena such as ocean circulation and air‐sea transports. MeasuringRfin the field is usually difficult, thus parameterization ofRfbased on readily observed properties is essential. Here, estimates ofRfin 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 forRfin terms of the buoyancy Reynolds number (Reb), the gradient Richardson number (Ri), and turbulent Froude number (Frt). Neither theReb‐based nor theRi‐based scheme is able to describe the observed variations inRf, but theFrt‐based parameterization works well. These findings support further use of theFrt‐ based parameterization in turbulent oceanic and estuarine environments.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1007/s10652-022-09873-2
