Previous analyses of the flow of low-Reynolds-number, viscous gravity currents down inclined planes are investigated further and extended. Particular emphasis is on the motion of the fluid front and tail, which previous analyses treated somewhat cavalierly. We obtain reliable, approximate, analytic solutions in these regions, the accuracies of which are satisfactorily tested against our numerical evaluations. The solutions show that the flow has several time scales determined by the inclination angle, $\alpha$ . At short times, the influence of initial and boundary conditions is important and the flow is governed by both the pressure gradient and the direct action of gravity due to inclination. Later on, the areas where the boundary conditions are important shrink. This fact explains why previous solutions, being inaccurate near the front and the tail, described experimental data with high accuracy. At larger times, of the order of $\alpha ^{-5/2}$ , the influence of the pressure gradient may be neglected and the fluid profile converges to the square-root shape predicted in previous works. Important extensions of our approach are also outlined.
Many flows, including those containing suspended particles, are kept turbulent by the action of the bottom stress, and this turbulence is also responsible for maintaining sedimenting particles in suspension and in some cases entraining more particles from the bed. A convenient one-dimensional analogue of these processes is provided by laboratory experiments conducted in a mixing box, where a characterizable turbulence is generated by the vertical oscillation of a horizontal grid. In the present paper we report the results of a series of experiments with a grid located close to the bottom boundary to simulate the action of stresses acting at a rough boundary, and compare the results with those obtained using the more extensively studied geometry in which a similar grid is located in the interior of a stirred fluid layer. Experiments have been conducted both with dense, particle-free fluid layers and with layers containing sufficiently high concentrations of dense particles to have a significant effect on the bulk density. In the fluid case, the interface at the top of the stirred dense layer continues to rise as lighter fluid is entrained across the interface. Sediment layers are distinctly different, because the particles responsible for the density difference between the layers can fall out of the suspension as it changes in thickness. The work done in keeping particles in suspension and the effect of this on the turbulence above the grid must be taken into account. The mechanism of resuspension of particles depends on the level of turbulence near the bottom boundary, below the grid. As the stirring rate, and thus the intensity of turbulence, are increased three possible equilibrium states can be attained sequentially: the particles eventually all precipitate; or some particles precipitate while the remainder are held indefinitely in suspension; or all the particles are suspended. In the last two cases a stable, self-limited suspension layer is produced, separated from the overlying fluid by a sharp density interface at a fixed height. Theoretical arguments are presented which provide a satisfactory scaling of the experimental data. These are compared with previous theories and numerical experiments aimed at modelling both the one-dimensional problem and the corresponding processes in turbulent gravity currents. Comparisons are also made with sediment-laden channel flows and convecting layers containing sedimenting particles. Similar results will hold for light, positively buoyant particles or non-coalescing bubbles.
We consider theoretically the long-time evolution of axisymmetric, high Reynolds number, Boussinesq gravity currents supplied by a constant, small-area source of mass and radial momentum in a deep, quiescent ambient. We describe the gravity currents using a shallow-water model with a Froude number closure condition to incorporate ambient form drag at the front and present numerical and asymptotic solutions. The predicted profile consists of an expanding, radially decaying, steady interior that connects via a shock to a deeper, self-similar frontal boundary layer. Controlled by the balance of interior momentum flux and frontal buoyancy across the shock, the front advances as (g′sQ/r1/4s)4/154/5, where g′s is the reduced gravity of the source fluid, Q is the total volume flux, rs is the source radius and is time. A radial momentum source has no effect on this solution below a non-zero threshold value. Above this value, the (virtual) radius over which the flow becomes critical can be used to collapse the solution onto the subthreshold one. We also use a simple parameterization to incorporate the effect of interfacial entrainment, and show that the profile can be substantially modified, although the buoyancy profile and radial extent are less significantly impacted. Our predicted profiles and extents are in reasonable agreement with existing experiments.
The propagation at high Reynolds number of a heavy, two-dimensional gravity current of given initial volume at the base of a uniform flow is considered. An experimental setup is described for which a known volume of fluid is rapidly introduced halfway down a 9 m channel in which there is a uniform flow of water. The density excess of the released fluid is produced by either dissolving salt or suspending particles in water. The upstream and downstream propagation of the current was measured for different initial salt concentrations, particle sizes and concentrations. A simple box model for the motion of and deposit from the gravity current is constructed. The analytical results obtained compare well with our numerical solutions of one-layer and two-layer models incorporating the appropriate shallow-water equations. Both sets of results are in very good agreement with the experimental data.
Gravity currents created by the release of a fixed volume of a suspension into a lighter ambient fluid are studied theoretically and experimentally. The greater density of the current and the buoyancy force driving its motion arise primarily from dense particles suspended in the interstitial fluid of the current. The dynamics of the current are assumed to be dominated by a balance between inertial and buoyancy forces; viscous forces are assumed negligible. The currents considered are two-dimensional and flow over a rigid horizontal surface. The flow is modelled by either the single- or the two-layer shallow-water equations, the two-layer equations being necessary to include the effects of the overlying fluid, which are important when the depth of the current is comparable to the depth of the overlying fluid. Because the local density of the gravity current depends on the concentration of particles, the buoyancy contribution to the momentum balance depends on the variation of the particle concentration. A transport equation for the particle concentration is derived by assuming that the particles are vertically well-mixed by the turbulence in the current, are advected by the mean flow and settle out through the viscous sublayer at the bottom of the current. The boundary condition at the moving front of the current relates the velocity and the pressure head at that point. The resulting equations are solved numerically, which reveals that two types of shock can occur in the current. In the late stages of all particle-driven gravity currents, an internal bore develops that separates a particle-free jet-like flow in the rear from a dense gravity-current flow near the front. The second type of bore occurs if the initial height of the current is comparable to the depth of the ambient fluid. This bore develops during the early lock-exchange flow between the two fluids and strongly changes the structure of the current and its transport of particles from those of a current in very deep surroundings. To test the theory, several experiments were performed to measure the length of particle-driven gravity currents as a function of time and their deposition patterns for a variety of particle sizes and initial masses of sediment. The comparison between the theoretical predictions, which have no adjustable parameters, and the experimental results are very good.
Recent studies have modelled the flow of particle–driven gravity currents over horizontal boundaries using either shallow–water equations or simple 'box' (integral) models. The shallow–water equations are typically integrated numerically, whereas box models admit analytical solutions. However, the theoretical validity of the latter models has not been fully established. In this paper a novel mathematical technique is developed which permits the derivation of analytical solutions to the shallow–water model for gravity current motion. These solutions, confirmed by comparison with the results of numerical integration, are in good agreement with experimental observations. They also indicate why the simplified box models have been so successful. Moreover, they reveal how the internal dynamics of particle–driven flows are different from gravity currents arising solely due to compositional density differences. While compositionally driven gravity currents, which have a fixed density difference between the intruding and ambient fluids, may be modelled using similarity solutions to the governing equations, particle–driven gravity currents do not possess such solutions because their density is progressively reduced by particle sedimentation. Instead the new analysis determines how their behaviour progressively diverges from the similarity solution. By a change of independent variables, it is possible to develop convergent series expansions for each of the dependent variables which characterize the motion. It is suggested that this approach may find application to a number of other problems in which the dynamics are initially governed by a simple dynamical balance which is progressively lost as extra physical effects begin to influence the motion.
'%3e%3cpath fill='%23012169' d='M0 0v30h60V0z'/%3e%3cpath stroke='%23fff' stroke-width='6' d='m0 0 60 30m0-30L0 30'/%3e%3cpath stroke='%23C8102E' stroke-width='4' d='m0 0 60 30m0-30L0 30' clip-path='url(%23b)'/%3e%3cpath stroke='%23fff' stroke-width='10' d='M30 0v30M0 15h60'/%3e%3cpath stroke='%23C8102E' stroke-width='6' d='M30 0v30M0 15h60'/%3e%3c/g%3e%3c/svg%3e)
'/%3e%3c/defs%3e%3cpath fill='%23012169' d='M0 0h10080v5040H0z'/%3e%3cpath stroke='%23fff' stroke-width='.6' d='m0 0 6 3m0-3L0 3' clip-path='url(%23b)' transform='scale(840)'/%3e%3cpath stroke='%23e4002b' stroke-width='.4' d='m0 0 6 3m0-3L0 3' clip-path='url(%23c)' transform='scale(840)'/%3e%3cpath stroke='%23fff' stroke-width='840' d='M2520 0v2520M0 1260h5040'/%3e%3cpath stroke='%23e4002b' stroke-width='504' d='M2520 0v2520M0 1260h5040'/%3e%3cg fill='%23fff'%3e%3cuse xlink:href='%23d' x='2520' y='3780'/%3e%3cuse xlink:href='%23a' x='7560' y='4200'/%3e%3cuse xlink:href='%23a' x='6300' y='2205'/%3e%3cuse xlink:href='%23a' x='7560' y='840'/%3e%3cuse xlink:href='%23a' x='8680' y='1869'/%3e%3cuse xlink:href='%23e' x='8064' y='2730'/%3e%3c/g%3e%3c/svg%3e)
