An analytical vorticity-based model is introduced for steady-state inviscid Boussinesq gravity currents in sheared ambients. The model enforces the conservation of mass and horizontal and vertical momentum, and it does not require any empirical closure assumptions. As a function of the given gravity current height, upstream ambient shear and upstream ambient layer thicknesses, the model predicts the current velocity as well as the downstream ambient layer thicknesses and velocities. In particular, it predicts the existence of gravity currents with a thickness greater than half the channel height, which is confirmed by direct numerical simulation (DNS) results and by an analysis of the energy loss in the flow. For high-Reynolds-number gravity currents exhibiting Kelvin–Helmholtz instabilities along the current/ambient interface, the DNS simulations suggest that for a given shear magnitude, the current height adjusts itself such as to allow for maximum energy dissipation.
Numerical simulation of gravity current propagation in a compressible atmosphere