The effects of Prandtl number on the nonlinear dynamics of Kelvin-Helmholtz instability in two dimensions

It is known that the pitchfork bifurcation of Kelvin-Helmholtz instability occurring at minimum gradient Richardson number in viscous stratified shear flows can be subcritical or supercritical depending on the value of the Prandtl number,. Here, we study stratified shear flow restricted to two dimensions at finite Reynolds number, continuously forced to have a constant background density gradient and a hyperbolic tangent shear profile, corresponding to the 'Drazin model' base flow. Bifurcation diagrams are produced for fluids with (typical for air), 3 and (typical for water). For and, steady billow-like solutions are found to exist for strongly stable stratification of beyond. Interestingly, these solutions are not a direct product of a Kelvin-Helmholtz instability, having half the wavelength of the linear instability, and arising through a superharmonic bifurcation. These short-wavelength states can be tracked down to at least and act as instigators of complex dynamics, even in strongly stratified flows. Direct numerical simulations of forced and unforced two-dimensional flows are performed, which support the results of the bifurcation analyses. Perturbations are observed to grow approximately exponentially from random initial conditions where no modal instability is predicted by a linear stability analysis.</p

Environmental Restoration Disposal Facility Lessons Learned

The purpose of lessons learned is to identify insight gained during a project – successes or failures – that can be applied on future projects. Lessons learned can contribute to the overall success of a project by building on approaches that have worked well and avoiding previous mistakes. Below are examples of lessons learned during ERDF’s ARRA-funded expansion project

AN INVESTIGATION INTO THE EFFECTS OF A SIMULATED EFFUSION IN HEALTHY SUBJECTS ON KNEE KINEMATICS AND LOWER LIMB MUSCLE ACTIVITY DURING A SINGLE LEG DROP LANDING

Arthrogenic muscle inhibition (AMI) is defined as an ongoing reflex inhibition of the musculature surrounding a joint following distension or damage to the structures of that joint [Hopkins and Ingersoll, 2000]. AMI following joint injury may affect movement and muscle recruitment which may impair rehabilitation and delay the return to activity. Knee angular displacement and velocity as well as lower limb EMG were measured in the period 250 milliseconds pre initial contact to 250 milliseconds post initial contact during a single leg drop jump in 8 healthy subjects before and after a simulated knee joint effusion of 60 millilitres. Repeated measures ANOVA and post hoc testing revealed no statistically significant differences in pre and post effusion in knee kinematic or lower limb EMG measures undertaken. A simulated knee effusion did not result in significant alterations to knee joint mechanics or lower limb muscle activation patterns during a single leg drop landing. The mechanism by which an effusion affects motor control during functional and dynamic weight bearing tasks warrants further investigation

Quantifying mixing and available potential energy in vertically periodic simulations of stratified flows

Turbulent mixing exerts a significant influence on many physical processes in the ocean. In a stably stratified Boussinesq fluid, this irreversible mixing describes the conversion of available potential energy (APE) to background potential energy (BPE). In some settings the APE framework is difficult to apply and approximate measures are used to estimate irreversible mixing. For example, numerical simulations of stratified turbulence often use triply periodic domains to increase computational efficiency. In this setup however, BPE is not uniquely defined and the method of Winters et al. (1995, J. Fluid Mech., 289) cannot be directly applied to calculate the APE. We propose a new technique to calculate APE in periodic domains with a mean stratification. By defining a control volume bounded by surfaces of constant buoyancy, we can construct an appropriate background buoyancy profile $b_\ast(z,t)$ and accurately quantify diapycnal mixing in such systems. This technique also permits the accurate calculation of a finite amplitude local APE density in periodic domains. The evolution of APE is analysed in various turbulent stratified flow simulations. We show that the mean dissipation rate of buoyancy variance $\chi$ provides a good approximation to the mean diapycnal mixing rate, even in flows with significant variations in local stratification. When quantifying measures of mixing efficiency in transient flows, we find significant variation depending on whether laminar diffusion of a mean flow is included in the kinetic energy dissipation rate. We discuss how best to interpret these results in the context of quantifying diapycnal diffusivity in real oceanographic flows.Comment: 28 pages, 10 figures, accepted to J. Fluid Mec

Shear-induced breaking of internal gravity waves

Motivated by observations of turbulence in the strongly stratified ocean thermocline, we use direct numerical simulations to investigate the interaction of a sinusoidal shear flow and a large-amplitude internal gravity wave. Despite strong nonlinearities in the flow and a lack of scale separation, we find that linear ray tracing theory is qualitatively useful in describing the early development of the flow as the wave is refracted by the shear. Consistent with the linear theory, the energy of the wave accumulates in regions of negative mean shear where we observe evidence of convective and shear instabilities. Streamwise-aligned convective rolls emerge the fastest, but their contribution to irreversible mixing is dwarfed by shear-driven billow structures that develop later. Although the wave strongly distorts the buoyancy field on which these billows develop, the mixing efficiency of the subsequent turbulence is similar to that arising from Kelvin-Helmholtz instability in a stratified shear layer. We run simulations at Reynolds numbers of 5000 and 8000, and vary the initial amplitude of the internal gravity wave. For high values of initial wave amplitude, the results are qualitatively independent of $Re$. Smaller initial wave amplitudes delay the onset of the instabilities, and allow for significant laminar diffusion of the internal wave, leading to reduced turbulent activity. We discuss the complex interaction between the mean flow, internal gravity wave and turbulence, and its implications for internal wave-driven mixing in the ocean.Comment: 27 pages, 12 figures, accepted to J. Fluid. Mec

Optimal perturbation growth on a breaking internal gravity wave

The breaking of internal gravity waves in the abyssal ocean is thought to be responsible for much of the mixing necessary to close oceanic buoyancy budgets. The exact mechanism by which these waves break down into turbulence remains an active area of research and can have significant implications on the mixing efficiency. Recent evidence has suggested that both shear instabilities and convective instabilities play a significant role in the breaking of an internal gravity wave in a high Richardson number mean shear flow. We perform a systematic analysis of the stability of a configuration of an internal gravity wave superimposed on a background shear flow first considered by Howland et al. (J. Fluid Mech., vol. 921, 2021, A24), using direct–adjoint looping to find the perturbation giving maximal energy growth on this evolving flow. We find that three-dimensional, convective mechanisms produce greater energy growth than their two-dimensional counterparts. In particular, we find close agreement with the direct numerical simulations of Howland et al. (J. Fluid Mech., 2021, in press), which demonstrated a clear three-dimensional mechanism causing breakdown to turbulence. The results are shown to hold at realistic Prandtl numbers. At low mean Richardson numbers, two-dimensional, shear-driven mechanisms produce greater energy growth

Optimal perturbation growth on a breaking internal gravity wave

The breaking of internal gravity waves in the abyssal ocean is thought to be responsible for much of the mixing necessary to close oceanic buoyancy budgets. The exact mechanism by which these waves break down into turbulence remains an active area of research and can have significant implications on the mixing efficiency. Recent evidence has suggested that both shear instabilities and convective instabilities play a significant role in the breaking of an internal gravity wave in a high Richardson number mean shear flow. We perform a systematic analysis of the stability of a configuration of an internal gravity wave superimposed on a background shear flow first considered by Howland et al. (J. Fluid Mech., vol. 921, 2021, A24), using direct–adjoint looping to find the perturbation giving maximal energy growth on this evolving flow. We find that three-dimensional, convective mechanisms produce greater energy growth than their two-dimensional counterparts. In particular, we find close agreement with the direct numerical simulations of Howland et al. (J. Fluid Mech., 2021, in press), which demonstrated a clear three-dimensional mechanism causing breakdown to turbulence. The results are shown to hold at realistic Prandtl numbers. At low mean Richardson numbers, two-dimensional, shear-driven mechanisms produce greater energy growth