Particle-driven gravity currents in non-rectangular cross section channels

Abstract

We consider a high-Reynolds-number gravity current generated by suspension of heavier particles in fluid of density ρi, propagating along a channel into an ambient fluid of the density ρa. The bottom and top of the channel are at z = 0, H, and the cross section is given by the quite general −f1(z) ≤ y ≤ f2(z) for 0 ≤ z ≤ H. The flow is modeled by the one-layer shallow-water equations obtained for the time-dependent motion which is produced by release from rest of a fixed volume of mixture from a lock. We solve the problem by the finite-difference numerical code to present typical height h(x, t), velocity u(x, t), and volume fraction of particles (concentration) ϕ(x, t) profiles. The methodology is illustrated for flow in typical geometries: power-law (f(z) = zα and f(z) = (H − z)α, where α is positive constant), trapezoidal, and circle. In general, the speed of propagation of the flows driven by suspensions decreases compared with those driven by a reduced gravity in homogeneous currents. However, the details depend on the geometry of the cross section. The runout length of suspensions in channels of power-law cross sections is analytically predicted using a simplified depth-averaged “box” model. The present approach is a significant generalization of the classical gravity current problem. The classical formulation for a rectangular channel is now just a particular case, f(z) = const., in the wide domain of cross sections covered by this new model.