The boundary layer instability beneath internal solitary waves and its sensitivity to vortex wakes

We investigated the stability of the bottom boundary layer (BBL) beneath periodic internal solitary waves (ISWs) of depression over a flat bottom through two-dimensional direct numerical simulations. We explored the convective versus absolute/global nature of the BBL instability in response to changes in Reynolds number, and the sensitivity of the instability to seeding noise in the front of the ISW – spanning laboratory to geophysical scales. The BBL was laminar at $Re_{ISW}=90$ and convectively unstable at $Re_{ISW}=300$ . At laboratory-scale $Re_{ISW}=300$ , the convective wave packet was periodically amplified by each successive ISW, until vortex shedding occurred. The associated noise-amplification behaviour potentially explains the discrepancies of the critical $Re_{ISW}$ between the lock–release laboratory experiments and our Dubreil–Jacotin–Long-initialized numerical simulations as the result of the difference in background noise. Instability energy decreased under the front shoulder of the ISW, analogous to flow relaminarization under a favourable pressure gradient. At geophysical-scale $Re_{ISW}=900$ , the BBL was initially convectively unstable, and then the instability tracked with the ISW, appearing phenomenologically similar to a global instability. The simulated initial convective instability at both $Re_{ISW}=300$ and $Re_{ISW}=900$ is in agreement with local linear stability analysis which predicts that the instability group speed is always lower than the ISW celerity. Increased free stream perturbations in front of the ISW and larger $Re_{ISW}$ shift the location of vortex shedding (and enhanced bed shear stress) beneath the wave, closer to the ISW trough, thereby potentially changing the location of maximum sediment resuspension, in agreement with field observations at higher $Re_{ISW}$ .