Gravity currents: entrainment, stratification and self-similarity

2015 ◽  
Vol 784 ◽  
pp. 130-162 ◽  
Author(s):  
Diana Sher ◽  
Andrew W. Woods

We present new experiments of the motion of a turbulent gravity current produced by the rapid release of a finite volume of dense aqueous solution from a lock of length $L$ into a channel $x>0$ filled with a finite depth, $H$, of fresh water. Using light attenuation we measure the mixing and evolving density of the flow, and, using dye studies, we follow the motion of the current and the ambient fluid. After the fluid has slumped to the base of the tank, there are two phases of the flow. When the front of the current, $x_{n}$, is within the region $2L<x_{n}<7L$, the fluid in the head of the current retains its original density and the flow travels with a constant speed. We find that approximately $0.75(\pm 0.05)$ of the ambient fluid displaced by the head mixes with the fluid in the head. The mixture rises over the head and feeds a growing stratified tail region of the flow. Dye studies show that fluid with the original density continues to reach the front of the current, at a speed which we estimate to be approximately $1.35\pm 0.05$ times that of the front, consistent with data of Berson (Q. J. R. Meteorol. Soc., vol. 84, 1958, pp. 1–16) and Kneller et al. (J. Geophys. Res. Oceans, vol. 104, 1999, pp. 5281–5291). This speed is similar to that of the ‘bore’, the trailing edge of the original lock gate fluid, as described by Rottman & Simpson (J. Fluid Mech., vol. 135, 1983, pp. 95–110). The continual mixing at the front leads to a gradual decrease of the mass of unmixed original lock gate fluid. Eventually, when the nose extends beyond $x_{n}\approx 7L$, the majority of the lock gate fluid has been diluted through the mixing. As the current continues, it adjusts to a second regime in which the position of the head increases with time as $x_{n}\approx 1.7B^{1/3}t^{2/3}$, where $B$ is the total buoyancy of the flow per unit width across the channel, while the depth-averaged reduced gravity in the head decreases through mixing with the ambient fluid according to the relation $g_{n}^{\prime }\approx 4.6H^{-1}B^{2/3}t^{-2/3}$. Measurements also show that the depth of the head $h_{n}(t)$ is approximately constant, $h_{n}\sim 0.38H$, and the reduced gravity decreases with height above the base of the current and with distance behind the front of the flow. Using the depth-averaged shallow-water equations, we derive a new class of self-similar solution which models the lateral structure of the flow by assuming the ambient fluid is entrained into the current in the head of the flow. By comparison with our data, we estimate that a fraction $0.69\pm 0.06$ of the ambient fluid displaced by the head of the current is mixed into the flow in this approximately self-similar regime, and the front of the current has a Froude number $0.9\pm 0.05$. We discuss the implications of our results for the evolution of the buoyancy in a gravity current as a function of the distance from the source.

2015 ◽  
Vol 782 ◽  
Author(s):  
Peeradon Samasiri ◽  
Andrew W. Woods

We present new experiments to measure the rate of entrainment of ambient fluid into a high Reynolds number, axisymmetric, turbulent gravity current. The current is produced by the rapid release of a finite volume of aqueous salt solution from a lock of length $r_{o}$ into a diverging channel, $r>0$, of angle $9.5^{\circ }$, filled with a finite depth, $H$, of fresh water. Using light attenuation we measure the evolving density of the flow, and using dye studies we illustrate the process of mixing between the current and ambient fluid. After an initial adjustment, a circulation develops in the head of the flow: current fluid reaches the nose of the flow, rises up and moves backwards relative to the nose. We find that, owing to the mixing, the volume of the current increases as $V\sim 0.2r_{n}^{7/4}r_{o}^{1/4}H$ while the maximum depth of the head decreases as $h_{n}\sim 0.5H(r_{o}/r_{n})^{1/4}$, where $r_{n}$ is the location of the front of the current. Combining these results, we estimate that the recirculating current fluid mixes with a fraction $E=0.33\pm 0.09$ of the ambient fluid that is directly ahead of the current and displaced upwards by it. Some of the mixed fluid supplies the tail of the flow, while the remainder recirculates into the head, which becomes progressively more dilute. In accord with Huppert & Simpson (J. Fluid Mech., vol. 99, 1980, pp. 785–799), we find that the position of the front increases with time as $r_{n}\approx (1.28\pm 0.05)B^{1/4}t^{1/2}$, where $B$ is the total buoyancy of the flow. We also find that the maximum value of the vertical integral of the buoyancy $(\overline{g^{\prime }}h)_{n}$ decreases with the position of the nose according to the relation $(\overline{g^{\prime }}h)_{n}\approx (0.89\pm 0.12)Br_{n}^{-2}$, consistent with a Froude number $0.86\pm 0.07$. We compare our measurements with a new idealised self-similar solution of the depth-averaged equations that accounts for the mixing at the nose, the vertical shear in the velocity and the lateral stratification of the buoyancy within the current.


1992 ◽  
Vol 4 (6) ◽  
pp. 1148-1155 ◽  
Author(s):  
Javier Alberto Diez ◽  
Roberto Gratton ◽  
Julio Gratton

2002 ◽  
Vol 456 ◽  
pp. 33-48 ◽  
Author(s):  
RYAN J. LOWE ◽  
P. F. LINDEN ◽  
JAMES W. ROTTMAN

Laboratory experiments were performed in which an intrusive gravity current was observed using shadowgraph and particle tracking methods. The intrusion was generated in a two-layer fluid with a sharp interface by mixing the fluid behind a vertical lock gate and then suddenly withdrawing the gate from the tank. The purpose of the experiments was to determine the structure of the velocity field inside the intrusion and the stability characteristics of the interface. Soon after the removal of the lock gate, the front of the intrusive gravity current travelled at a constant speed close to the value predicted by theory for an energy-conserving gravity current. The observed structure of the flow inside the intrusion can be divided into three regions. At the front of the intrusion there is an energy-conserving head region in which the fluid velocity is nearly uniform with speed equal to the front speed. This is followed by a dissipative wake region in which large billows are present with their associated mixing and in which the fluid velocity is observed to be non-uniform and have a maximum speed approximately 50% greater than the front speed. Behind the wake region is a tail region in which there is very little mixing and the velocity field is nearly uniform with a speed slightly faster than the front speed.


2015 ◽  
Vol 778 ◽  
pp. 669-690 ◽  
Author(s):  
Zhong Zheng ◽  
Sangwoo Shin ◽  
Howard A. Stone

We study the propagation of viscous gravity currents along a thin permeable substrate where slow vertical drainage is allowed from the boundary. In particular, we report the effect of this vertical fluid drainage on the second-kind self-similar solutions for the shape of the fluid–fluid interface in three contexts: (i) viscous axisymmetric gravity currents converging towards the centre of a cylindrical container; (ii) viscous gravity currents moving towards the origin in a horizontal Hele-Shaw channel with a power-law varying gap thickness in the horizontal direction; and (iii) viscous gravity currents propagating towards the origin of a porous medium with horizontal permeability and porosity gradients in power-law forms. For each of these cases with vertical leakage, we identify a regime diagram that characterizes whether the front reaches the origin or not; in particular, when the front does not reach the origin, we calculate the final location of the front. We have also conducted laboratory experiments with a cylindrical lock gate to generate a converging viscous gravity current where vertical fluid drainage is allowed from various perforated horizontal substrates. The time-dependent position of the propagating front is captured from the experiments, and the front position is found to agree well with the theoretical and numerical predictions when surface tension effects can be neglected.


2001 ◽  
Vol 444 ◽  
pp. 23-47 ◽  
Author(s):  
DAVID PRITCHARD ◽  
ANDREW W. WOODS ◽  
ANDREW J. HOGG

We examine the gravitational dispersal of dense fluid through a horizontal permeable layer, which is separated from a second underlying layer by a narrow band of much lower permeability. We derive a series of analytical solutions which describe the propagation of the fluid through the upper layer and the draining of the fluid into the underlying region. The model predicts that the current initially spreads according to a self-similar solution. However, as the drainage becomes established, the spreading slows, and in fact the fluid only spreads a finite distance before it has fully drained into the underlying layer. We examine the sensitivity of the results to the initial conditions through numerical solution of the governing equations. We find that for sources of sufficiently large initial aspect ratio (defined as the ratio of height to length), the solution converges rapidly to the initially self-similar regime. For longer and shallower initial source conditions, this convergence does not occur, but we derive estimates for the run-out length of the current, which compare favourably with our numerical data. We also present some preliminary laboratory experiments, which support the model.


2007 ◽  
Vol 577 ◽  
pp. 363-383 ◽  
Author(s):  
M. A. HESSE ◽  
H. A. TCHELEPI ◽  
B. J. CANTWEL ◽  
F. M. ORR

We investigate the evolution of a finite release of fluid into an infinite, two-dimensional, horizontal, porous slab saturated with a fluid of different density and viscosity. The vertical boundaries of the slab are impermeable and the released fluid spreads as a gravity current along a horizontal boundary. At early times the released fluid fills the entire height of the layer, and the governing equation admits a self-similar solution that is a function of the viscosity ratio between the two fluids. This early similarity solution describes a tilting interface with tips propagating as x ∝ t1/2. At late times the released fluid has spread along the boundary and the height of the current is much smaller than the thickness of the layer. The governing equation simplifies and admits a different similarity solution that is independent of the viscosity ratio. This late similarity solution describes a point release of fluid in a semi-infinite porous half-space, where the tip of the interface propagates as x ∝ t1/3. The same simplification of the governing equation occurs if the viscosity of the released fluid is much higher than the viscosity of the ambient fluid. We have obtained an expression for the time when the solution transitions from the early to the late self-similar regime. The transition time increases monotonically with increasing viscosity ratio. The transition period during which the solution is not self-similar also increases monotonically with increasing viscosity ratio, for mobility ratios larger than unity. Numerical computations describing the full evolution of the governing equation show good agreement with the theoretical results. Estimates of the spreading of injected fluids over long times are important for geological storage of CO2, and for the migration of pollutants in aquifers. In all cases it is important to be able to anticipate when the spreading regime transitions from x ∝ t1/2 to x ∝ t1/3.


1992 ◽  
Vol 3 (4) ◽  
pp. 319-341 ◽  
Author(s):  
S. P. Hastings ◽  
L. A. Peletier

We discuss the self-similar solutions of the second kind associated with the propagation of turbulent bursts in a fluid at rest. Such solutions involve an eigenvalue parameter μ, which cannot be determined from dimensional analysis. Existence and uniqueness are established and the dependence of μ on a physical parameter λ in the problem is studied: estimates are obtained and the asymptotic behaviour as λ → ∞ is established.


Sign in / Sign up

Export Citation Format

Share Document