An official website of the United States government
In a stably stratified shear layer, multiple competing instabilities produce sensitivity to small changes in initial conditions, popularly called the butterfly effect (as a flapping wing may alter the weather). Three ensembles of 15 simulated mixing events, identical but for small perturbations to the initial state, are used to explore differences in the route to turbulence, the maximum turbulence level and the total amount and efficiency of mixing accomplished by each event. Comparisons show that a small change in the initial state alters the strength and timing of the primary Kelvin–Helmholtz instability, the subharmonic pairing instability and the various three-dimensional secondary instabilities that lead to turbulence. The effect is greatest in, but not limited to, the parameter regime where pairing and the three-dimensional secondary instabilities are in strong competition. Pairing may be accelerated or prevented; maximum turbulence kinetic energy may vary by up to a factor of 4.6, flux Richardson number by 12 %–15 % and net mixing by a factor of 2.
more »
« less
The effects of boundary proximity on Kelvin–Helmholtz instability and turbulence
Studies of Kelvin–Helmholtz (KH) instability have typically modelled the initial flow as an isolated shear layer. In geophysical cases, however, the instability often occurs near boundaries and may therefore be influenced by boundary proximity effects. Ensembles of direct numerical simulations are conducted to understand the effect of boundary proximity on the evolution of the instability and the resulting turbulence. Ensemble averages are used to reduce sensitivity to small variations in initial conditions. Both the transition to turbulence and the resulting turbulent mixing are modified when the shear layer is near a boundary: the time scales for the onset of instability and turbulence are longer, and the height of the KH billow is reduced. Subharmonic instability is suppressed by the boundary because phase lock is prevented due to the diverging phase speeds of the KH and subharmonic modes. In addition, the disruptive influence of three-dimensional secondary instabilities on pairing is more profound as the two events coincide more closely. When the shear layer is far from the boundary, the shear-aligned convective instability is dominant; however, secondary central-core instability takes over when the shear layer is close to the boundary, providing an alternate route for the transition to turbulence. Both the efficiency of the resulting mixing and the turbulent diffusivity are dramatically reduced by boundary proximity effects.
more »
« less
- Award ID(s):
- 1830071
- PAR ID:
- 10437708
- Date Published:
- Journal Name:
- Journal of Fluid Mechanics
- Volume:
- 966
- ISSN:
- 0022-1120
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
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 Physical processes behind flow‐topography interactions and turbulent transitions are essential for parameterization in numerical models. We examine how the Kuroshio cascades energy into turbulence upon passing over a seamount, employing a combination of shipboard measurements, tow‐yo microstructure profiling, and high‐resolution mooring. The seamount, spanning 5 km horizontally with two summits, interacts with the Kuroshio, whose flow speed ranges from 1 to 2 m s−1, modulated by tides. The forward energy cascade process is commenced by forming a train of 2–3 nonlinear lee waves behind the summit with a wavelength of 0.5–1 km and an amplitude of 50–100 m. A train of Kelvin‐Helmholtz (KH) billows develops immediately below the lee waves and extends downstream, leading to enhanced turbulence. The turbulent kinetic energy dissipation rate isO(10−7–10−4) W kg−1, varying in phase with the upstream flow speed modulated by tides. KH billows occur primarily at the lee wave's trailing edge, where the combined strong downstream shear and low‐stratification recirculation trigger the shear instability,Ri < 1/4. The recirculation also creates an overturn susceptible to gravitational instability. This scenario resembles the rotor, commonly found in atmospheric mountain waves but rarely observed in the ocean. A linear stability analysis further suggests that critical levels, where the KH instability extracts energy from the mean flow, are located predominantly at the strong shear layer of the lee wave's upwelling portion, coinciding with the upper boundary of the rotor. These novel observations may provide insights into flow‐topography interactions and improve physics‐based turbulence parameterization.more » « less
-
Abstract Kelvin‐Helmholtz Instabilities (KHI) are known to be significant drivers of atmospheric turbulence. Recent observations show KHI forming with misaligned or angled billow segments that develop connecting vortex tubes and knots (T&K); these features promote distinctive, event‐defining instability and mixing characteristics that were not accounted for in prior idealized studies. Though T&K have been shown to increase mixing in KHI events with low Richardson numbers (Ri), their influence in weakly KH‐unstable, less‐idealized environments is unknown. Here we present modeling results of KHI in the stratosphere to assess the impact of T&K dynamics in weakly KH‐unstable environments. Radiosonde wind and temperature profiles from 22 February 2006 near Lamont, Oklahoma, measured vertically offset shear and stability peaks near 16.2 km with a minimum Ri = 0.11. Direct numerical simulations (DNS) of this event reveal decreasing shear and increasing stratification, where Ri increases to 0.2 as the shear and stratification peaks move to a common altitude. The resulting KHI exhibit T&K features forming adjacent to, and in superposition with, secondary convective instabilities (CI) rather than superseding them as in prior T&K studies with Ri = 0.05. Newly discovered “crankshaft” instabilities distort the billows and give rise to secondary KHI with delayed, elevated dissipation. KHI that exhibit T&K dynamics are found to accumulate % greater mixing than axially uniform KHI with equal or lower mixing efficiency. The substantial increase in mixing suggests significant contributions of T&K dynamics to KHI events throughout the atmosphere that remain unaddressed in general circulation models' turbulence parameterizations.more » « less
-
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
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1017/jfm.2023.412
