Stratified shear instability in a field of pre-existing turbulence

2019 ◽  
Vol 862 ◽  
pp. 639-658 ◽  
Author(s):  
A. K. Kaminski ◽  
W. D. Smyth

Turbulent mixing of heat and momentum in the stably-stratified ocean interior occurs in discrete events driven by vertical variations of the horizontal velocity. Typically, these events have been modelled assuming an initially laminar stratified shear flow which develops wavelike instabilities, becomes fully turbulent, and then relaminarizes into a stable state. However, in the real ocean there is always some level of turbulence left over from previous events. Using direct numerical simulations, we show that the evolution of a stably-stratified shear layer may be significantly modified by pre-existing turbulence. The classical billow structure associated with Kelvin–Helmholtz instability is suppressed and eventually eliminated as the strength of the initial turbulence is increased. A corresponding energetics analysis shows that potential energy changes and dissipation of kinetic energy depend non-monotonically on initial turbulence strength, with the largest effects when initial turbulence is present but insufficient to prevent billow formation. The mixing efficiency decreases with increasing initial turbulence amplitude as the development of the Kelvin–Helmholtz billow, with its large pre-turbulent mixing efficiency, is arrested.

2017 ◽  
Vol 826 ◽  
pp. 522-552 ◽  
Author(s):  
A. Mashayek ◽  
C. P. Caulfield ◽  
W. R. Peltier

Turbulent mixing plays a major role in enabling the large-scale ocean circulation. The accuracy of mixing rates estimated from observations depends on our understanding of basic fluid mechanical processes underlying the nature of turbulence in a stratified fluid. Several of the key assumptions made in conventional mixing parameterizations have been increasingly scrutinized in recent years, primarily on the basis of adequately high resolution numerical simulations. We add to this evidence by compiling results from a suite of numerical simulations of the turbulence generated through stratified shear instability processes. We study the inherently intermittent and time-dependent nature of wave-induced turbulent life cycles and more specifically the tight coupling between inherently anisotropic scales upon which small-scale isotropic turbulence grows. The anisotropic scales stir and stretch fluid filaments enhancing irreversible diffusive mixing at smaller scales. We show that the characteristics of turbulent mixing depend on the relative time evolution of the Ozmidov length scale $L_{O}$ compared to the so-called Thorpe overturning scale $L_{T}$ which represents the scale containing available potential energy upon which turbulence feeds and grows. We find that when $L_{T}\sim L_{O}$, the mixing is most active and efficient since stirring by the largest overturns becomes ‘optimal’ in the sense that it is not suppressed by ambient stratification. We argue that the high mixing efficiency associated with this phase, along with observations of $L_{O}/L_{T}\sim 1$ in oceanic turbulent patches, together point to the potential for systematically underestimating mixing in the ocean if the role of overturns is neglected. This neglect, arising through the assumption of a clear separation of scales between the background mean flow and small-scale quasi-isotropic turbulence, leads to the exclusion of an highly efficient mixing phase from conventional parameterizations of the vertical transport of density. Such an exclusion may well be significant if the mechanism of shear-induced turbulence is assumed to be representative of at least some turbulent events in the ocean. While our results are based upon simulations of shear instability, we show that they are potentially more generic by making direct comparisons with $L_{T}-L_{O}$ data from ocean and lake observations which represent a much wider range of turbulence-inducing physical processes.


2017 ◽  
Vol 837 ◽  
pp. 341-380 ◽  
Author(s):  
Peter P. Sullivan ◽  
James C. McWilliams

The evolution of upper ocean currents involves a set of complex, poorly understood interactions between submesoscale turbulence (e.g. density fronts and filaments and coherent vortices) and smaller-scale boundary-layer turbulence. Here we simulate the lifecycle of a cold (dense) filament undergoing frontogenesis in the presence of turbulence generated by surface stress and/or buoyancy loss. This phenomenon is examined in large-eddy simulations with resolved turbulent motions in large horizontal domains using${\sim}10^{10}$grid points. Steady winds are oriented in directions perpendicular or parallel to the filament axis. Due to turbulent vertical momentum mixing, cold filaments generate a potent two-celled secondary circulation in the cross-filament plane that is frontogenetic, sharpens the cross-filament buoyancy and horizontal velocity gradients and blocks Ekman buoyancy flux across the cold filament core towards the warm filament edge. Within less than a day, the frontogenesis is arrested at a small width,${\approx}100~\text{m}$, primarily by an enhancement of the turbulence through a small submesoscale, horizontal shear instability of the sharpened filament, followed by a subsequent slow decay of the filament by further turbulent mixing. The boundary-layer turbulence is inhomogeneous and non-stationary in relation to the evolving submesoscale currents and density stratification. The occurrence of frontogenesis and arrest are qualitatively similar with varying stress direction or with convective cooling, but the detailed evolution and flow structure differ among the cases. Thus submesoscale filament frontogenesis caused by boundary-layer turbulence, frontal arrest by frontal instability and frontal decay by forward energy cascade, and turbulent mixing are generic processes in the upper ocean.


2001 ◽  
Vol 428 ◽  
pp. 349-386 ◽  
Author(s):  
E. J. STRANG ◽  
H. J. S. FERNANDO

The results of a laboratory experiment designed to study turbulent entrainment at sheared density interfaces are described. A stratified shear layer, across which a velocity difference ΔU and buoyancy difference Δb is imposed, separates a lighter upper turbulent layer of depth D from a quiescent, deep lower layer which is either homogeneous (two-layer case) or linearly stratified with a buoyancy frequency N (linearly stratified case). In the parameter ranges investigated the flow is mainly determined by two parameters: the bulk Richardson number RiB = ΔbD/ΔU2 and the frequency ratio fN = ND=ΔU.When RiB > 1.5, there is a growing significance of buoyancy effects upon the entrainment process; it is observed that interfacial instabilities locally mix heavy and light fluid layers, and thus facilitate the less energetic mixed-layer turbulent eddies in scouring the interface and lifting partially mixed fluid. The nature of the instability is dependent on RiB, or a related parameter, the local gradient Richardson number Rig = N2L/ (∂u/∂z)2, where NL is the local buoyancy frequency, u is the local streamwise velocity and z is the vertical coordinate. The transition from the Kelvin–Helmholtz (K-H) instability dominated regime to a second shear instability, namely growing Hölmböe waves, occurs through a transitional regime 3.2 < RiB < 5.8. The K-H activity completely subsided beyond RiB ∼ 5 or Rig ∼ 1. The transition period 3.2 < RiB < 5 was characterized by the presence of both K-H billows and wave-like features, interacting with each other while breaking and causing intense mixing. The flux Richardson number Rif or the mixing efficiency peaked during this transition period, with a maximum of Rif ∼ 0.4 at RiB ∼ 5 or Rig ∼ 1. The interface at 5 < RiB < 5.8 was dominated by ‘asymmetric’ interfacial waves, which gradually transitioned to (symmetric) Hölmböe waves at RiB > 5:8.Laser-induced fluorescence measurements of both the interfacial buoyancy flux and the entrainment rate showed a large disparity (as large as 50%) between the two-layer and the linearly stratified cases in the range 1.5 < RiB < 5. In particular, the buoyancy flux (and the entrainment rate) was higher when internal waves were not permitted to propagate into the deep layer, in which case more energy was available for interfacial mixing. When the lower layer was linearly stratified, the internal waves appeared to be excited by an ‘interfacial swelling’ phenomenon, characterized by the recurrence of groups or packets of K-H billows, their degeneration into turbulence and subsequent mixing, interfacial thickening and scouring of the thickened interface by turbulent eddies.Estimation of the turbulent kinetic energy (TKE) budget in the interfacial zone for the two-layer case based on the parameter α, where α = (−B + ε)/P, indicated an approximate balance (α ∼ 1) between the shear production P, buoyancy flux B and the dissipation rate ε, except in the range RiB < 5 where K-H driven mixing was active.


2006 ◽  
Vol 128 (4) ◽  
pp. 874-879 ◽  
Author(s):  
Roberto C. Aguirre ◽  
Jennifer C. Nathman ◽  
Haris C. Catrakis

Flow geometry effects are examined on the turbulent mixing efficiency quantified as the mixture fraction. Two different flow geometries are compared at similar Reynolds numbers, Schmidt numbers, and growth rates, with fully developed turbulence conditions. The two geometries are the round jet and the single-stream planar shear layer. At the flow conditions examined, the jet exhibits an ensemble-averaged mixing efficiency which is approximately double the value for the shear layer. This substantial difference is explained fluid mechanically in terms of the distinct large-scale entrainment and mixing-initiation environments and is therefore directly due to flow geometry effects.


2011 ◽  
Vol 41 (7) ◽  
pp. 1364-1384 ◽  
Author(s):  
W. D. Smyth ◽  
S. Kimura

Abstract Mixing due to sheared salt fingers is studied by means of direct numerical simulations (DNS) of a double-diffusively unstable shear layer. The focus is on the “moderate shear” case, where shear is strong enough to produce Kelvin–Helmholtz (KH) instability but not strong enough to produce the subharmonic pairing instability. This flow supports both KH and salt-sheet instabilities, and the objective is to see how the two mechanisms work together to flux heat, salt, and momentum across the layer. For observed values of the bulk Richardson number Ri and the density ratio Rρ, the linear growth rates of KH and salt-sheet instabilities are similar. These mechanisms, as well as their associated secondary instabilities, lead the flow to a fully turbulent state. Depending on the values of Ri and Rρ, the resulting turbulence may be driven mainly by shear or mainly by salt fingering. Turbulent mixing causes the profiles of temperature, salinity, and velocity to spread; however, in salt-sheet-dominated cases, the net density (or buoyancy) layer thins over time. This could be a factor in the maintenance of the staircase and is also an argument in favor of an eventual role for Holmboe instability. Fluxes are scaled using both laboratory scalings for a thin layer and an effective diffusivity. Fluxes are generally stronger in salt-sheet-dominated cases. Shear instability disrupts salt-sheet fluxes while adding little flux of its own. Shear therefore reduces mixing, despite providing an additional energy source. The dissipation ratio Γ is near 0.2 for shear-dominated cases but is much larger when salt sheets are dominant, supporting the use of Γ in the diagnosis of observed mixing phenomena. The profiler approximation Γz, however, appears to significantly overestimate the true dissipation ratio.


2011 ◽  
Vol 354-355 ◽  
pp. 559-563
Author(s):  
Lei Shi ◽  
Shen Jie Zhou ◽  
Feng Ling Yang ◽  
Fan Jin Hu

Mixing efficiency is an important parameter in the design of many industrial processes in stirred tanks. In this study, CFD technology was used to simulate the mixing process inside the stirred tank with dislocated blades and standard turbine. Calculations were performed to study the effects of agitator speed and the configuration of impellers on mixing efficiency. The results showed that the flow field in the stirred tank with the dislocated blades is better than the standard turbine, and the flow number of the dislocated blades had been improved while the power number had been reduced. According to calculation results of Wr, we found the mixing efficiency of the dislocated blades had been improved about 4 times than that of standard turbine.


2017 ◽  
Vol 815 ◽  
pp. 243-256
Author(s):  
Philippe Odier ◽  
Robert E. Ecke

Stratified shear flows occur in many geophysical contexts, from oceanic overflows and river estuaries to wind-driven thermocline layers. We explore a turbulent wall-bounded shear flow of lighter miscible fluid into a quiescent fluid of higher density with a range of Richardson numbers$0.05\lesssim Ri\lesssim 1$. In order to find a stability parameter that allows close comparison with linear theory and with idealized experiments and numerics, we investigate different definitions of$Ri$. We find that a gradient Richardson number defined on fluid interface sections where there is no overturning at or adjacent to the maximum density gradient position provides an excellent stability parameter, which captures the Miles–Howard linear stability criterion. For small$Ri$the flow exhibits robust Kelvin–Helmholtz instability, whereas for larger$Ri$interfacial overturning is more intermittent with less frequent Kelvin–Helmholtz events and emerging Holmboe wave instability consistent with a thicker velocity layer compared with the density layer. We compute the perturbed fraction of interface as a quantitative measure of the flow intermittency, which is approximately 1 for the smallest$Ri$but decreases rapidly as$Ri$increases, consistent with linear theory. For the perturbed regions, we use the Thorpe scale to characterize the overturning properties of these flows. The probability distribution of the non-zero Thorpe length yields a universal exponential form, suggesting that much of the overturning results from increasingly intermittent Kelvin–Helmholtz instability events. The distribution of turbulent kinetic energy, conditioned on the intermittency fraction, has a similar form, suggesting an explanation for the universal scaling collapse of the Thorpe length distribution.


2007 ◽  
Vol 577 ◽  
pp. 127-136 ◽  
Author(s):  
E. GUYEZ ◽  
J.-B. FLOR ◽  
E. J. HOPFINGER

Experiments conducted on mixing across a stable density interface in a turbulent Taylor–Couette flow show, for the first time, experimental evidence of an increase in mixing efficiency at large Richardson numbers. With increasing buoyancy gradient the buoyancy flux first passes a maximum, then decreases and at large values of the buoyancy gradient the flux increases again. Thus, the curve of buoyancy flux versus buoyancy gradient tends to be N-shaped (rather than simply bell shaped), a behaviour suggested by the model of Balmforth et al. (J. Fluid Mech. vol. 428, 1998, p. 349). The increase in mixing efficiency at large Richardson numbers is attributed to a scale separation of the eddies active in mixing at the interface; when the buoyancy gradient is large mean kinetic energy is injected at scales much smaller than the eddy size fixed by the gap width, thus decreasing the eddy turnover time. Observations show that there is no noticeable change in interface thickness when the mixing efficiency increases; it is the mixing mechanism that changes. The curves of buoyancy flux versus buoyancy gradient also show a large variability for identical experimental conditions. These variations occur at time scales one to two orders of magnitude larger than the eddy turnover time scale.


2017 ◽  
Vol 837 ◽  
pp. 129-146 ◽  
Author(s):  
Jason Olsthoorn ◽  
Stuart B. Dalziel

The study of vortex-ring-induced mixing has been significant for understanding stratified turbulent mixing in the absence of a mean flow. Renewed interest in this topic has prompted the development of a one-dimensional model for the evolution of a stratified system in the context of isolated mixing events. This model is compared to numerical simulations and physical experiments of vortex rings interacting with a stratification. Qualitative agreement between the evolution of the density profiles is observed, along with close quantitative agreement of the mixing efficiency. This model highlights the key dynamical features of such isolated mixing events.


2020 ◽  
Author(s):  
Jiahao Wang ◽  
Xi Chen ◽  
Kefeng Mao ◽  
Kelan Zhu

Abstract. Satellite measurements during April to June in 2019 and direct observations from 28th to 30th May in 2019 about the Kuroshio Extension Front are conducted. The former shows the front experience a process of stable-unstable-stable state caused by the movement of the Kuroshio Extension’s second meander and a pinched-off eddy. The latter indicates the steep upward slopes of the isopycnals tilt northward in the strong frontal zone as well as several over 100 m thick blobs of cold and fresh water in the salinity minimum zone of North Pacific Intermediate Water. Using isopycnal anomaly method and diapycnal spiciness curvature method, characteristic interleaving layers are shown primarily in σθ = 26.3–26.9 kg/m3, which corresponds to large variations of potential spiciness in intermediate layers. Further analysis indicates the development of thermohaline intrusions may be driven by the double diffusive instability and the velocity anomalies. Besides, we find the turbulence mixing attributed to symmetric instability and shear instability is very strong in intermediate layer.


Sign in / Sign up

Export Citation Format

Share Document