The effects of drag on turbulent gravity currents

2000 ◽  
Vol 416 ◽  
pp. 297-314 ◽  
Author(s):  
LYNNE HATCHER ◽  
ANDREW J. HOGG ◽  
ANDREW W. WOODS
Keyword(s):  
Analytical Solution ◽  
Laboratory Experiments ◽  
Gravity Current ◽  
Gravity Currents ◽  
Similarity Solutions ◽  
Flow Speed ◽  
Series Solutions ◽  
Accurate Approximation ◽  
New Class ◽  
Power Series Solutions

We model the propagation of turbulent gravity currents through an array of obstacles which exert a drag force on the flow proportional to the square of the flow speed. A new class of similarity solutions is constructed to describe the flows that develop from a source of strength q0tγ. An analytical solution exists for a finite release, γ = 0, while power series solutions are developed for sources with γ > 0. These are shown to provide an accurate approximation to the numerically calculated similarity solutions. The model is successfully tested against a series of new laboratory experiments which investigate the motion of a turbulent gravity current through a large flume containing an array of obstacles. The model is extended to account for the effects of a sloping boundary. Finally, a series of geophysical and environmental applications of the model are discussed.

Propagation of viscous currents on a porous substrate with finite capillary entry pressure

2016 ◽  
Vol 801 ◽  
pp. 65-90 ◽  
Author(s):  
Roiy Sayag ◽  
Jerome A. Neufeld
Keyword(s):  
Laboratory Experiments ◽  
Gravity Current ◽  
Fluid Pressure ◽  
Analytical Techniques ◽  
Gravity Currents ◽  
Porous Substrate ◽  
Entry Pressure ◽  
Self Similar ◽  
Capillary Entry Pressure ◽  
Theoretical Results

We study the propagation of viscous gravity currents over a thin porous substrate with finite capillary entry pressure. Near the origin, where the current is deep, propagation of the current coincides with leakage through the substrate. Near the nose of the current, where the current is thin and the fluid pressure is below the capillary entry pressure, drainage is absent. Consequently the flow can be characterised by the evolution of drainage and fluid fronts. We analyse this flow using numerical and analytical techniques combined with laboratory-scale experiments. At early times, we find that the position of both fronts evolve as $t^{1/2}$, similar to an axisymmetric gravity current on an impermeable substrate. At later times, the growing effect of drainage inhibits spreading, causing the drainage front to logarithmically approach a steady position. In contrast, the asymptotic propagation of the fluid front is quasi-self-similar, having identical structure to the solution of gravity currents on an impermeable substrate, only with slowly varying fluid flux. We benchmark these theoretical results with laboratory experiments that are consistent with our modelling assumption, but that also highlight the detailed dynamics of drainage inhibited by finite capillary pressure.

Entraining gravity currents

2013 ◽  
Vol 731 ◽  
pp. 477-508 ◽  
Author(s):  
Christopher G. Johnson ◽  
Andrew J. Hogg
Keyword(s):  
Richardson Number ◽  
Solution Structure ◽  
Large Scale ◽  
Gravity Current ◽  
Gravity Currents ◽  
Similarity Solutions ◽  
Water Model ◽  
Experimental Measurements ◽  
Entrainment Rate ◽  
Long Time

AbstractEntrainment of ambient fluid into a gravity current, while often negligible in laboratory-scale flows, may become increasingly significant in large-scale natural flows. We present a theoretical study of the effect of this entrainment by augmenting a shallow water model for gravity currents under a deep ambient with a simple empirical model for entrainment, based on experimental measurements of the fluid entrainment rate as a function of the bulk Richardson number. By analysing long-time similarity solutions of the model, we find that the decrease in entrainment coefficient at large Richardson number, due to the suppression of turbulent mixing by stable stratification, qualitatively affects the structure and growth rate of the solutions, compared to currents in which the entrainment is taken to be constant or negligible. In particular, mixing is most significant close to the front of the currents, leading to flows that are more dilute, deeper and slower than their non-entraining counterparts. The long-time solution of an inviscid entraining gravity current generated by a lock-release of dense fluid is a similarity solution of the second kind, in which the current grows as a power of time that is dependent on the form of the entrainment law. With an entrainment law that fits the experimental measurements well, the length of currents in this entraining inviscid regime grows with time approximately as ${t}^{0. 447} $. For currents instigated by a constant buoyancy flux, a different solution structure exists in which the current length grows as ${t}^{4/ 5} $. In both cases, entrainment is most significant close to the current front.

scholarly journals Stratified gravity currents in porous media

2016 ◽  
Vol 791 ◽  
pp. 329-357 ◽  
Author(s):  
Samuel S. Pegler ◽  
Herbert E. Huppert ◽  
Jerome A. Neufeld
Keyword(s):  
Porous Media ◽  
Density Distribution ◽  
Vertical Structure ◽  
Gravity Current ◽  
Gravity Currents ◽  
Similarity Solutions ◽  
Density Field ◽  
Vertical Position ◽  
Uniform Density ◽  
Convex Shape

We consider theoretically and experimentally the propagation in porous media of variable-density gravity currents containing a stably stratified density field, with most previous studies of gravity currents having focused on cases of uniform density. New thin-layer equations are developed to describe stably stratified fluid flows in which the density field is materially advected with the flow. Similarity solutions describing both the fixed-volume release of a distributed density stratification and the continuous input of fluid containing a distribution of densities are obtained. The results indicate that the density distribution of the stratification significantly influences the vertical structure of the gravity current. When more mass is distributed into lighter densities, it is found that the shape of the current changes from the convex shape familiar from studies of the uniform-density case to a concave shape in which lighter fluid accumulates primarily vertically above the origin of the current. For a constant-volume release, the density contours stratify horizontally, a simplification which is used to develop analytical solutions. For currents introduced continuously, the horizontal velocity varies with vertical position, a feature which does not apply to uniform-density gravity currents in porous media. Despite significant effects on vertical structure, the density distribution has almost no effect on overall horizontal propagation, for a given total mass. Good agreement with data from a laboratory study confirms the predictions of the model.

Gravity currents from a line source in an ambient flow

2008 ◽  
Vol 606 ◽  
pp. 1-26 ◽  
Author(s):  
ANJA C. SLIM ◽  
HERBERT E. HUPPERT
Keyword(s):  
Line Source ◽  
Source Region ◽  
Gravity Currents ◽  
Similarity Solutions ◽  
Flow Speed ◽  
Moving Frame ◽  
Ambient Flow ◽  
Constant Rate ◽  
Long Time ◽  
High Reynolds

We present a mainly theoretical study of high-Reynolds-number planar gravity currents in a uniformly flowing deep ambient. The gravity currents are generated by a constant line source of fluid, and may also be supplied with a source of horizontal momentum and a source of particles. We model the motion using a shallow-water approximation and represent the effects of the ambient flow by imposing a Froude-number condition in a moving frame. We present analytic and numerical expressions for the threshold ambient flow speed above which no upstream propagation can occur at long times. For homogeneous gravity currents in an ambient flow below threshold, we find similarity solutions in which the up- and downstream fronts spread at a constant rate and the current propagates indefinitely in both directions. For gravity currents consisting of both interstitial fluid of a different density to the ambient and a sedimenting particle load, we find long-time asymptotic solutions for ambient flow strengths below threshold. These consist of a steady particle-rich near-source region, in which settling and advection of particles balance, and an effectively particle-free frontal region. The homogeneous behaviour of the fronts ensures that they also spread at a constant rate and therefore can propagate upstream indefinitely. For gravity currents driven solely by a sedimenting particle load, we find numerically that a single regime exists for ambient flow strengths below threshold. In these solutions, settling balances advection near the source leading to a steady region, which joins on to a complex frontal boundary layer. The upstream front progressively decelerates. Our solutions for homogeneous and particle-driven gravity currents compare well with published experimental results.

The effect of confining impermeable boundaries on gravity currents in a porous medium

2010 ◽  
Vol 649 ◽  
pp. 1-17 ◽  
Author(s):  
MADELEINE J. GOLDING ◽  
HERBERT E. HUPPERT
Keyword(s):  
Free Surface ◽  
Gravity Current ◽  
Mass Conservation ◽  
Propagation Rate ◽  
Gravity Currents ◽  
Similarity Solutions ◽  
Constant Flux ◽  
Channel Shape ◽  
Dimensional Gravity ◽  
Uniform Cross Section

The effect of confining boundaries on gravity currents in porous media is investigated theoretically and experimentally. Similarity solutions are derived for currents when the volume increases as tα in horizontal channels of uniform cross-section with boundary height b satisfying b ~ a|y/a|n, where y is the cross-channel coordinate and a is a length scale of the channel width. Experiments were carried out in V-shaped and semicircular channels for the case of gravity currents with constant volume (α=0) and constant flux (α=1). These showed generally good agreement with the theory.Typically, we find that the propagation of the current is well described by L ~ tc for some scalar c. We study the dependence of c on the time exponent of the volume of fluid in the current, α, and the geometry of the channel, parameterized by n. For all channel shapes, there exists a critical value of α, αc = 1/2, above which increasing n causes an increase in c and below which increasing n causes a decrease in c, where increasing n corresponds to opening up the channel boundary to the horizontal. The current height increases or decreases with respect to time depending on whether α is greater or less than αc. It is this fact, along with global mass conservation, which explains why varying the channel shape n affects the propagation rate c in different ways depending on α.We also consider channels inclined at an angle θ to the horizontal. When the slope of the channel is much greater than the slope of the free surface of the current, the component of gravity parallel to the slope dominates, causing the current to move with a constant velocity, Vf say, regardless of channel shape n and flux parameter α, in agreement with results for a two-dimensional gravity current obtained by Huppert & Woods (1995) and some initially axisymmetric gravity currents presented by Vella & Huppert (2006). If the effect of the component of gravity perpendicular to the channel may not be neglected, i.e. if the slopes of the channel and free surface of the current are comparable, we find that, in a frame moving with speed Vf, the form of the governing equation for the height of a current in an equivalent horizontal channel is recovered. We calculate that the height of a constant flux gravity current down an inclined channel will tend to a fixed depth, which is determined by the channel shape, n, and the physical properties of the fluid and rock. Experimental and numerical results for inclined V-shaped channels agree very well with this theory.

Particle-driven gravity currents down planar slopes

1999 ◽  
Vol 390 ◽  
pp. 75-91 ◽  
Author(s):  
ROGER T. BONNECAZE ◽  
JOHN R. LISTER
Keyword(s):  
Gravity Current ◽  
Gravity Currents ◽  
Similarity Solutions ◽  
Bottom Friction ◽  
Turbidity Currents ◽  
Scaling Analysis ◽  
Ambient Fluid ◽  
Buoyancy Driven Flows ◽  
Depositional Patterns ◽  
Planar Slope

Particle-driven gravity currents, as exemplified by either turbidity currents in the ocean or ignimbrite flows in the atmosphere, are buoyancy-driven flows due to a suspension of dense particles in an ambient fluid. We present a theoretical study on the dynamics of and deposition from a turbulent current flowing down a uniform planar slope from a constant-flux point source of particle-laden fluid. The flow is modelled using the shallow-water equations, including the effects of bottom friction and entrainment of ambient fluid, coupled to an equation for the transport and settling of the particles. Two flow regimes are identified. Near the source and for mild slopes, the flow is dominated by a balance between buoyancy and bottom friction. Further downstream and for steeper slopes, entrainment also affects the behaviour of the current. Similarity solutions are also developed for the simple cases of homogeneous gravity currents with no settling of particles in the friction-dominated and entrainment-dominated regimes. Estimates of the width and length of the deposit from a monodisperse particle-driven gravity current with settling are derived from scaling analysis for each regime, and the contours of the depositional patterns are determined from numerical solution of the governing equations.

scholarly journals Gravity current propagation up a valley

2014 ◽  
Vol 762 ◽  
pp. 417-434 ◽  
Author(s):  
Catherine S. Jones ◽  
Claudia Cenedese ◽  
Eric P. Chassignet ◽  
P. F. Linden ◽  
Bruce R. Sutherland
Keyword(s):  
Numerical Simulations ◽  
Laboratory Experiments ◽  
Rectangular Channel ◽  
Gravity Current ◽  
Gravity Currents ◽  
Momentum Transport ◽  
Front Speed

AbstractThe advance of the front of a dense gravity current propagating in a rectangular channel and V-shaped valley both horizontally and up a shallow slope is examined through theory, full-depth lock–release laboratory experiments and hydrostatic numerical simulations. Consistent with theory, experiments and simulations show that the front speed is relatively faster in the valley than in the channel. The front speed measured shortly after release from the lock is 5–22 % smaller than theory, with greater discrepancy found in upsloping V-shaped valleys. By contrast, the simulated speed is approximately 6 % larger than theory, showing no dependence on slope for rise angles up to ${\it\theta}=8^{\circ }$. Unlike gravity currents in a channel, the current head is observed in experiments to be more turbulent when propagating in a V-shaped valley. The turbulence is presumably enhanced due to the lateral flows down the sloping sides of the valley. As a consequence, lateral momentum transport contributes to the observed lower initial speeds. A Wentzel–Kramers–Brillouin like theory predicting the deceleration of the current as it runs upslope agrees remarkably well with simulations and with most experiments, within errors.

scholarly journals Dynamics of viscous grounding lines

2010 ◽  
Vol 648 ◽  
pp. 363-380 ◽  
Author(s):  
ROSALYN A. V. ROBISON ◽  
HERBERT E. HUPPERT ◽  
M. GRAE WORSTER
Keyword(s):  
Laboratory Experiments ◽  
Gravity Current ◽  
Ice Sheets ◽  
Gravity Currents ◽  
Solid Base ◽  
Leading Order ◽  
Ice Shelves ◽  
Grounding Line ◽  
Theoretical Predictions ◽  
Good Agreement

We have used viscous fluids in simple laboratory experiments to explore the dynamics of grounding lines between marine ice sheets and the freely floating ice shelves into which they develop. We model the ice sheets as shear-dominated gravity currents, and the ice shelves as extensional gravity currents having zero shear to leading order. We consider the flow of viscous fluid down an inclined plane into a dense inviscid ‘ocean’. A fixed flux of fluid is supplied at the top of the plane, which is at ‘sea level’. The fluid forms a gravity current flowing down and attached to the plane for some distance before detaching to form a freely floating extensional current. We have derived a mathematical model of the flow that incorporates a new dynamic boundary condition for the position of the grounding line, where the gravity current loses contact with the solid base. The grounding line initially advances and eventually reaches a steady position. Good agreement between our theoretical predictions and experimental measurements and observations gives confidence in the fundamental assumptions of our model.

The effect of confining boundaries on viscous gravity currents

2007 ◽  
Vol 577 ◽  
pp. 495-505 ◽  
Author(s):  
DAISUKE TAKAGI ◽  
HERBERT E. HUPPERT
Keyword(s):  
Laboratory Experiments ◽  
Gravity Currents ◽  
Similarity Solutions ◽  
Mathematical Equations ◽  
Inclined Channels

Newtonian viscous gravity currents propagating along horizontal and inclined channels with semicircular and V-shaped boundaries are examined. Similarity solutions are obtained from the governing mathematical equations and compared with closely matching data from laboratory experiments in which the propagation of glycerine along different channels was recorded. Geological applications of the results are discussed briefly.

A study of three-dimensional gravity currents on a uniform slope

2002 ◽  
Vol 453 ◽  
pp. 239-261 ◽  
Author(s):  
ANDREW N. ROSS ◽  
P. F. LINDEN ◽  
STUART B. DALZIEL
Keyword(s):  
Laboratory Experiments ◽  
Gravity Current ◽  
Three Dimensional ◽  
Gravity Currents ◽  
Integral Model ◽  
Conductivity Measurements ◽  
Dense Fluid ◽  
Wedge Shape ◽  
Dimensional Gravity ◽  
Shallow Water Models

In many geophysical, environmental and industrial situations, a finite volume of fluid with a density different to the ambient is released on a sloping boundary. This leads to the formation of a gravity current travelling up, down and across the slope. We present novel laboratory experiments in which the dense fluid spreads both down-slope (and initially up-slope) and laterally across the slope. The position, shape and dilution of the current are determined through video and conductivity measurements for moderate slopes (5° to 20°). The entrainment coefficient for different slopes is calculated from the experimental results and is found to depend very little on the slope. The value agrees well with previously published values for entrainment into gravity currents on a horizontal surface. The experimental measurements are compared with previous shallow-water models and with a new wedge integral model developed and presented here. It is concluded that these simplified models do not capture all the significant features of the flow. In the models, the current takes the form of a wedge which travels down the slope, but the experiments show the formation of a more complicated current. It is found that the wedge integral model over-predicts the length and width of the gravity current but gives fair agreement with the measured densities in the head. The initial stages of the flow, during which time the wedge shape develops, are studied. It is found that although the influence of the slope is seen relatively quickly for moderate slopes, the time taken for the wedge to develop is much longer. The implications of these findings for safety analysis are briefly discussed.

Sign in / Sign up

Export Citation Format

Share Document