Search for: All records

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. 

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. 

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. 

Abstract Acrucial region of the ocean surface boundary layer (OSBL) is the strongly-sheared and -stratified transition layer (TL) separating the mixed layer from the upper pycnocline, where a diverse range of waves and instabilities are possible. Previous work suggests that these different waves and instabilities will lead to different OSBL behaviours. Therefore, understanding which physical processes occur is key for modelling the TL. Here we present observations of the TL from a Lagrangian float deployed for 73 days near Ocean Weather Station Papa (50°N, 145°W) during Fall 2018. The float followed the vertical motion of the TL, continuously measuring profiles across it using an ADCP, temperature chain and salinity sensors. The temperature chain made depth/time images of TL structures with a resolution of 6cm and 3 seconds. These showed the frequent occurrence of very sharp interfaces, dominated by temperature jumps of O(1)°C over 6cm or less. Temperature inversions were typically small (≲ 10cm), frequent, and strongly-stratified; very few large overturns were observed. The corresponding velocity profiles varied over larger length scales than the temperature profiles. These structures are consistent with scouring behaviour rather than Kelvin-Helmholtz-type overturning. Their net effect, estimated via a Thorpe-scale analysis, suggests that these frequent small temperature inversions can account for the observed mixed layer deepening and entrainment flux. Corresponding estimates of dissipation, diffusivity, and heat fluxes also agree with previous TL studies, suggesting that the TL dynamics is dominated by these nearly continuous 10cm-scale mixing structures, rather than by less frequent larger overturns.