Hostname: page-component-76fb5796d-25wd4 Total loading time: 0 Render date: 2024-04-26T10:46:37.320Z Has data issue: false hasContentIssue false

A non-dissipative solution of Benjamin-type gravity current for a wide range of depth ratios

Published online by Cambridge University Press:  30 June 2011

MARIUS UNGARISH*
Affiliation:
Department of Computer Science, Technion, Haifa 32000, Israel
*
Email address for correspondence: unga@cs.technion.ac.il

Abstract

We consider the steady-state propagation of a high-Reynolds-number gravity current of height h and density ρc on the bottom of a horizontal channel of height H filled with ambient fluid of density ρa(<ρc), usually known as Benjamin's current problem. The objective is to derive an analytical result for the speed of propagation, U, in the form of the dimensionless Froude number, Fr(a) = U/(gh)1/2). Here g′ = (ρca − 1)g is the reduced gravity-driving effect (g being the gravity acceleration) and a = h/H is the depth (thickness) ratio of the layer of the current to that of the ambient fluid into which the current propagates. The analysis is performed in a frame of reference attached to the current; in this frame the current is a motionless slug. The original analysis of Benjamin assumes that the speed of the ambient in the domain above the parallel-horizontal main part of the current (behind the head) is independent of the vertical coordinate z, but here we assume that a small u′(z) fluctuation about the depth-averaged speed u exists. Then, we impose the balances of volume flux, flow-force (momentum flux) and global energy conservation, for a control volume attached to the current. We show that this gives a unique analytical result for Fr as a function of a = h/H. We recall that the original counterpart solution FrB(a) of Benjamin does not satisfy the above-mentioned energy conservation condition, i.e. the system displays energy dissipation (except for the half-depth current case a = 1/2). The present dissipationless-flow Fr(a) result is valid for any a ≤ 1/2, i.e. currents of at most half-depth of the channel height. On the other hand, in agreement with Benjamin's solution, gravity currents of more than half-depth of the channel height require an energy source and are impossible in normal conditions. The new Fr(a) is slightly smaller than Benjamin's FrB(a) result for 0 < a < 1/2, and the difference vanishes at a = 1/2 and a → 0 (a current of finite height in a very deep ambient).

Type
Papers
Copyright
Copyright © Cambridge University Press 2011

Access options

Get access to the full version of this content by using one of the access options below. (Log in options will check for institutional or personal access. Content may require purchase if you do not have access.)

References

REFERENCES

Huppert, H. E. & Simpson, J. E. 1980 The slumping of gravity currents. J. Fluid Mech. 99, 785799.CrossRefGoogle Scholar
Klemp, J. B., Rotunno, R. & Skamarock, W. C. 1994 On the dynamics of gravity currents in a channel. J. Fluid Mech. 269, 169198.CrossRefGoogle Scholar
Shin, J. O., Dalziel, S. B. & Linden, P. F. 2004 Gravity currents produced by lock exchange. J. Fluid Mech. 521, 134.CrossRefGoogle Scholar
Ungarish, M. 2008 Energy balances and front speed conditions of two-layer models for gravity currents produced by lock release. Acta Mechanica 201, 6381.CrossRefGoogle Scholar
Ungarish, M. 2009 An Introduction to Gravity Currents and Intrusions. Chapman & Hall.CrossRefGoogle Scholar
Ungarish, M. 2010 Energy balances for gravity currents with a jump at the interface produced by lock release. Acta Mechanica 211, 122.CrossRefGoogle Scholar
Ungarish, M. 2011 Two-layer shallow-water dam-break solutions for non-Boussinesq gravity currents in a wide range of fractional depth. J. Fluid Mech. 675, 2759.CrossRefGoogle Scholar