Hostname: page-component-7c8c6479df-hgkh8 Total loading time: 0 Render date: 2024-03-28T09:43:00.639Z Has data issue: false hasContentIssue false

Two-layer shallow-water dam-break solutions for gravity currents in non-rectangular cross-area channels

Published online by Cambridge University Press:  12 September 2013

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

Abstract

We consider the propagation of a high-Reynolds-number gravity current in a horizontal channel along the horizontal coordinate $x$. The bottom and top of the channel are at $z= 0, H$, and the cross-section is given by the quite general $- {f}_{1} (z)\leq y\leq {f}_{2} (z)$ for $0\leq z\leq H$. We develop a two-layer shallow-water (SW) formulation, and we implement it for the solution of the dam-break problem. The dependent variables are the position of the interface, $h(x, t)$, and the velocity (averaged over the area of the current), $u(x, t)$. The non-rectangular cross-section geometry enters the formulation via $f(h)$ and integrals of $f(z)$ and $zf(z)$, where $f(z)= {f}_{1} (z)+ {f}_{2} (z)$ is the width of the channel. For a given geometry $f(z)$, the free input parameters of the lock-release problem are: (i) the height ratio $H/ {h}_{0} $ of ambient to lock; and (ii) the density ratio of ‘light’ to ‘heavy’ fluids, $R$. We show that the dam-break problem is amenable to solution by the method of characteristics, but various complications (internal jumps, critical restrictions) appear when the return flow in the ambient is significant; these features are not captured by a one-layer shallow-water model. The role of the non-standard geometry in the detection and calculation of the internal and reflected jumps, and on the dam-break flow field, is elucidated. The methodology is illustrated for Boussinesq flow in typical geometries: power-law ($f(z)= b{z}^{\alpha } $, where $b, \alpha $ are positive constants), trapezoidal and circle (full or sector); we solved for full-depth ($H= 1$) and part-depth ($H\gt 1$) lock-release configurations. The approach seems to work well: in all the tested cases, the solution was obtained with simple mathematical and numerical tools (a Runge–Kutta integrator and a secant-method solver); the insights are compatible with our previous understanding of the problem; and the standard (rectangular) case is recovered in the appropriate limit (i.e. when the side-boundaries of the trapezium cross-section are vertical or far apart). A comparison of the speed of propagation with available experiments is also presented, and the agreement is satisfactory. The present solution is a significant generalization of the classical two-layer gravity–current dam-break problem. The classical formulation for a rectangular (or laterally unbounded) channel is now just a particular case, $f(z)= $ const., in the wide domain of cross-sections covered by the new model.

Type
Papers
Copyright
©2013 Cambridge University Press 

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

Benjamin, T. B. 1968 Gravity currents and related phenomena. J. Fluid Mech. 31, 209248.Google Scholar
Bonometti, T., Ungarish, M. & Balachandar, S. 2011 A numerical investigation of constant-volume non-Boussinesq gravity currents in deep ambient. J. Fluid Mech. 673, 574602.CrossRefGoogle Scholar
Borden, Z., Meiburg, E. & Constantinescu, G. 2012 Internal bores: an improved model via a detailed analysis of the energy budget. J. Fluid Mech. 703, 279314.Google Scholar
Klemp, J. B., Rotunno, R. & Skamarock, W. C. 1994 On the dynamics of gravity currents in a channel. J. Fluid Mech. 269, 169198.Google Scholar
Klemp, J. B., Rotunno, R. & Skamarock, W. C. 1997 On the propagation of internal bores. J. Fluid Mech. 331, 81106.Google Scholar
Lowe, R. J., Rottman, J. W. & Linden, P. F. 2005 The non-Boussinesq lock exchange problem. Part 1. Theory and experiments. J. Fluid Mech. 537, 101124.Google Scholar
Marino, B. M. & Thomas, L. P. 2009 Front condition for gravity currents in channels of nonrectangular symmetric cross-section shapes. Trans. ASME J. Fluids Engng 131 (5), 051201.Google Scholar
Marino, B. M. & Thomas, L. P. 2011 Dam-break release of a gravity current in a power-law channel section. J. Phys.: Conference Series 296, 012008 doi:10.1088/1742-6596/296/1/012008.Google Scholar
Monaghan, J., Mériaux, C., Huppert, H. & Monaghan, J. 2009 High Reynolds number gravity currents along V-shaped valleys. Eur. J. Mech. – B/Fluids 28 (5), 651659.CrossRefGoogle Scholar
Rottman, J. & Simpson, J. 1983 Gravity currents produced by instantaneous release of a heavy fluid in a rectangular channel. J. Fluid Mech. 135, 95110.CrossRefGoogle Scholar
Rotunno, R., Klemp, J. B., Bryan, G. H. & Muraki, D. J. 2011 Models of non-Boussinesq lock-exchange flow. J. Fluid Mech. 675, 126.Google Scholar
Simpson, J. E. 1997 Gravity Currents in the Environment and the Laboratory. Cambridge University Press.Google Scholar
Ungarish, M. 2009 An Introduction to Gravity Currents and Intrusions. Chapman & Hall/CRC Press.Google Scholar
Ungarish, M. 2010 Energy balances for gravity currents with a jump at the interface produced by lock release. Acta Mech. 211, 122.Google 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.Google Scholar
Ungarish, M. 2012 A general solution of Benjamin-type gravity current in a channel of non-rectangular cross-section. Environ. Fluid Mech. 12 (3), 251263.Google Scholar
Zemach, T. & Ungarish, M. 2013 Gravity currents in non-rectangular cross-section channels: analytical and numerical solutions of the one-layer shallow-water model for high-Reynolds-number propagation. Phys. Fluids 25, 026601–1–24.Google Scholar