On the dynamics of gravity currents in a channel

Abstract

We attempt to clarify the factors that regulate the propagation and structure of gravity currents through evaluation of idealized theoretical models along with two-dimensional numerical model simulations. In particular, we seek to reconcile research based on hydraulic theory for gravity currents evolving from a known initial state with analyses of gravity currents that are assumed to be at steady state, and to compare these approaches with both numerical simulations and laboratory experiments. The time-dependent shallow-water solution for a gravity current propagating in a channel of finite depth reveals that the flow must remain subcritical behind the leading edge of the current (in a framework relative to the head). This constraint requires that hf/d ≤ 0.347, where hf is the height of the front and d is the channel depth. Thus, in the lock-exchange problem, inviscid solutions corresponding to hf/d = 0.5 are unphysical, and the actual currents have depth ratios of less than one half near their leading edge and require dissipation or are not steady. We evaluate the relevance of Benjamin's (1968) well-known formula for the propagation of steady gravity currents and clarify discrepancies with other theoretical and observed results. From two-dimensional simulations with a frictionless lower surface, we find that Benjamin's idealized flow-force balance provides a good description of the gravity-current propagation. Including surface friction reduces the propagation speed because it produces dissipation within the cold pool. Although shallow-water theory over-estimates the propagation speed of the leading edge of cold fluid in the ‘dam-break’ problem, this discrepancy appears to arise from the lack of mixing across the current interface rather than from deficiencies in Benjamin's front condition. If an opposing flow restricts the propagation of a gravity current away from its source, we show that the propagation of the current relative to the free stream may be faster than predicted by Benjamin's formula. However, in these situations the front propagation remains dependent upon the specific source conditions and cannot be generalized.