A laboratory study of the velocity structure in an intrusive gravity current

2002 ◽  
Vol 456 ◽  
pp. 33-48 ◽  
Author(s):  
RYAN J. LOWE ◽  
P. F. LINDEN ◽  
JAMES W. ROTTMAN
Keyword(s):  
Velocity Field ◽  
Velocity Structure ◽  
Fluid Velocity ◽  
Gravity Current ◽  
Maximum Speed ◽  
Wake Region ◽  
Tail Region ◽  
Front Speed ◽  
Energy Conserving ◽  
Lock Gate

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.

A Laboratory Study of the Velocity Structure in an Intrusive Gravity Current

2000 ◽  
Author(s):  
Ryan J. Lowe
Keyword(s):  
Particle Tracking ◽  
Laboratory Experiments ◽  
Velocity Structure ◽  
Gravity Current ◽  
Current Theory ◽  
Head Region ◽  
Front Speed ◽  
Energy Conserving ◽  
The Stability ◽  
Lock Gate

Abstract 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 is to determine the structure of the velocity field inside the intrusion as well as the stability characteristics of the interface. Soon after the removal of the lock-gate the speed of the front of the intrusive gravity current reached a constant speed. The observed structure of the flow inside the intrusion shows a “head region” where the flow is nearly uniform, followed by a region of intense mixing and high velocities and finally followed by another region of fairly uniform velocity with a speed slightly faster than the front speed. The results show that the maximum centerline velocity is about 50% greater than the front speed and corresponds to the position in the intrusion where the strongest Kelvin- Helmholtz billows form. Closer to the front, the relative flow within the head is weak, which explains why Benjamin’s (1968) energy-conserving gravity current theory accurately predicts the behavior of dissipative gravity currents.

scholarly journals Gravity currents and internal waves in a stratified fluid

2008 ◽  
Vol 616 ◽  
pp. 327-356 ◽  
Author(s):  
BRIAN L. WHITE ◽  
KARL R. HELFRICH
Keyword(s):  
Internal Waves ◽  
Wave Speed ◽  
Numerical Solutions ◽  
Gravity Current ◽  
Gravity Currents ◽  
Stratified Fluid ◽  
Front Speed ◽  
Long Wave ◽  
Conjugate State ◽  
Energy Conserving

A steady theory is presented for gravity currents propagating with constant speed into a stratified fluid with a general density profile. Solution curves for front speed versus height have an energy-conserving upper bound (the conjugate state) and a lower bound marked by the onset of upstream influence. The conjugate state is the largest-amplitude nonlinear internal wave supported by the ambient stratification, and in the limit of weak stratification approaches Benjamin's energy-conserving gravity current solution. When the front speed becomes critical with respect to linear long waves generated above the current, steady solutions cannot be calculated, implying upstream influence. For non-uniform stratification, the critical long-wave speed exceeds the ambient long-wave speed, and the critical-Froude-number condition appropriate for uniform stratification must be generalized. The theoretical results demonstrate a clear connection between internal waves and gravity currents. The steady theory is also compared with non-hydrostatic numerical solutions of the full lock release initial-value problem. Some solutions resemble classic gravity currents with no upstream disturbance, but others show long internal waves propagating ahead of the gravity current. Wave generation generally occurs when the stratification and current speed are such that the steady gravity current theory fails. Thus the steady theory is consistent with the occurrence of either wave-generating or steady gravity solutions to the dam-break problem. When the available potential energy of the dam is large enough, the numerical simulations approach the energy-conserving conjugate state. Existing laboratory experiments for intrusions and gravity currents produced by full-depth lock exchange flows over a range of stratification profiles show excellent agreement with the conjugate state solutions.

Gravity currents: entrainment, stratification and self-similarity

2015 ◽  
Vol 784 ◽  
pp. 130-162 ◽  
Author(s):  
Diana Sher ◽  
Andrew W. Woods
Keyword(s):  
Light Attenuation ◽  
Gravity Current ◽  
Finite Depth ◽  
Reduced Gravity ◽  
Ambient Fluid ◽  
Tail Region ◽  
Self Similar Solution ◽  
Self Similar ◽  
Two Phases ◽  
Lock Gate

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.

Numerical study of the flow of a mixture of reacting gases in an inhomogeneous catalyst layer

2003 ◽  
Vol 3 ◽  
pp. 246-254
Author(s):  
C.I. Mikhaylenko ◽  
S.F. Urmancheev
Keyword(s):  
Velocity Field ◽  
Chemical Reactions ◽  
Porous Layer ◽  
Computational Experiment ◽  
Numerical Study ◽  
Fluid Velocity ◽  
Catalyst Layer ◽  
Granular Catalyst ◽  
Fluid Velocity Field

The behavior of a liquid flowing through a fixed bulk porous layer of a granular catalyst is considered. The effects of the nonuniformity of the fluid velocity field, which arise when the surface of the layer is curved, and the effect of the resulting inhomogeneity on the speed and nature of the course of chemical reactions are investigated by the methods of a computational experiment.

Wave Motion of a Compressible Viscous Fluid Contained in a Cylindrical Shell

1993 ◽  
Vol 115 (3) ◽  
pp. 302-312 ◽  
Author(s):  
J. H. Terhune ◽  
K. Karim-Panahi
Keyword(s):  
Velocity Field ◽  
Viscous Fluid ◽  
Imaginary Part ◽  
Fluid Velocity ◽  
Wave Number ◽  
Mode Shapes ◽  
Infinite Length ◽  
Complex Frequency ◽  
Compressible Viscous Fluid ◽  
Fluid Velocity Field

The free vibration of cylindrical shells filled with a compressible viscous fluid has been studied by numerous workers using the linearized Navier-Stokes equations, the fluid continuity equation, and Flu¨gge ’s equations of motion for thin shells. It happens that solutions can be obtained for which the interface conditions at the shell surface are satisfied. Formally, a characteristic equation for the system eigenvalues can be written down, and solutions are usually obtained numerically providing some insight into the physical mechanisms. In this paper, we modify the usual approach to this problem, use a more rigorous mathematical solution and limit the discussion to a single thin shell of infinite length and finite radius, totally filled with a viscous, compressible fluid. It is shown that separable solutions are obtained only in a particular gage, defined by the divergence of the fluid velocity vector potential, and the solutions are unique to that gage. The complex frequency dependence for the transverse component of the fluid velocity field is shown to be a result of surface interaction between the compressional and vortex motions in the fluid and that this motion is confined to the boundary layer near the surface. Numerical results are obtained for the first few wave modes of a large shell, which illustrate the general approach to the solution. The axial wave number is complex for wave propagation, the imaginary part being the spatial attenuation coefficient. The frequency is also complex, the imaginary part of which is the temporal damping coefficient. The wave phase velocity is related to the real part of the axial wave number and turns out to be independent of frequency, with numerical value lying between the sonic velocities in the fluid and the shell. The frequency dependencies of these parameters and fluid velocity field mode shapes are computed for a typical case and displayed in non-dimensional graphs.

Gravity currents and related phenomena

1968 ◽  
Vol 31 (2) ◽  
pp. 209-248 ◽  
Author(s):  
T. Brooke Benjamin
Keyword(s):  
Salt Water ◽  
Gravity Current ◽  
Circular Cross Section ◽  
Gravity Currents ◽  
Head Wave ◽  
Physical Problems ◽  
Energy Conserving ◽  
Fluid Theory ◽  
Hydrodynamical Problem ◽  
Upper Boundary

This paper presents a broad investigation into the properties of steady gravity currents, in so far as they can be represented by perfect-fluid theory and simple extensions of it (like the classical theory of hydraulic jumps) that give a rudimentary account of dissipation. As usually understood, a gravity current consists of a wedge of heavy fluid (e.g. salt water, cold air) intruding into an expanse of lighter fluid (fresh water, warm air); but it is pointed out in § 1 that, if the effects of viscosity and mixing of the fluids at the interface are ignored, the hydrodynamical problem is formally the same as that for an empty cavity advancing along the upper boundary of a liquid. Being simplest in detail, the latter problem is treated as a prototype for the class of physical problems under study: most of the analysis is related to it specifically, but the results thus obtained are immediately applicable to gravity currents by scaling the gravitational constant according to a simple rule.In § 2 the possible states of steady flow in the present category between fixed horizontal boundaries are examined on the assumption that the interface becomes horizontal far downstream. A certain range of flows appears to be possible when energy is dissipated; but in the absence of dissipation only one flow is possible, in which the asymptotic level of the interface is midway between the plane boundaries. The corresponding flow in a tube of circular cross-section is found in § 3, and the theory is shown to be in excellent agreement with the results of recent experiments by Zukoski. A discussion of the effects of surface tension is included in § 3. The two-dimensional energy-conserving flow is investigated further in § 4, and finally a close approximation to the shape of the interface is obtained. In § 5 the discussion turns to the question whether flows characterized by periodic wavetrains are realizable, and it appears that none is possible without a large loss of energy occurring. In § 6 the case of infinite total depth is considered, relating to deeply submerged gravity currents. It is shown that the flow must always feature a breaking ‘head wave’, and various properties of the resulting wake are demonstrated. Reasonable agreement is established with experimental results obtained by Keulegan and others.

Numerical simulations of lock-exchange compositional gravity current

2009 ◽  
Vol 635 ◽  
pp. 361-388 ◽  
Author(s):  
SENG KEAT OOI ◽  
GEORGE CONSTANTINESCU ◽  
LARRY WEBER
Keyword(s):  
High Resolution ◽  
Dissipation Rate ◽  
Friction Velocity ◽  
Gravity Current ◽  
Gravity Currents ◽  
Inviscid Limit ◽  
Grashof Number ◽  
Front Speed ◽  
Initial Volume ◽  
Bed Friction

Compositional gravity current flows produced by the instantaneous release of a finite-volume, heavier lock fluid in a rectangular horizontal plane channel are investigated using large eddy simulation. The first part of the paper focuses on the evolution of Boussinesq lock-exchange gravity currents with a large initial volume of the release during the slumping phase in which the front of the gravity current propagates with constant speed. High-resolution simulations are conducted for Grashof numbers $\sqrt {Gr}$ = 3150 (LGR simulation) and $\sqrt {Gr}$ = 126000 (HGR simulation). The Grashof number is defined with the channel depth h and the buoyancy velocity ub = $\sqrt {g'h}$ (g′ is the reduced gravity). In the HGR simulation the flow is turbulent in the regions behind the two fronts. Compared to the LGR simulation, the interfacial billows lose their coherence much more rapidly (over less than 2.5h behind the front), which results in a much faster decay of the large-scale content and turbulence intensity in the trailing regions of the flow. A slightly tilted, stably stratified interface layer develops away from the two fronts. The concentration profiles across this layer can be approximated by a hyperbolic tangent function. In the HGR simulation the energy budget shows that for t > 18h/ub the flow reaches a regime in which the total dissipation rate and the rates of change of the total potential and kinetic energies are constant in time. The second part of the paper focuses on the study of the transition of Boussinesq gravity currents with a small initial volume of the release to the buoyancy–inertia self-similar phase. When the existence of the back wall is communicated to the front, the front speed starts to decrease, and the current transitions to the buoyancy–inertia phase. Three high-resolution simulations are performed at Grashof numbers between $\sqrt {Gr}$ = 3 × 104 and $\sqrt {Gr}$ = 9 × 104. Additionally, a calculation at a much higher Grashof number ($\sqrt {Gr}$ = 106) is performed to understand the behaviour of a bottom-propagating current closer to the inviscid limit. The three-dimensional simulations correctly predict a front speed decrease proportional to t−α (the time t is measured from the release time) over the buoyancy–inertia phase, with the constant α approaching the theoretical value of 1/3 as the current approaches the inviscid limit. At Grashof numbers for which $\sqrt {Gr}$ > 3 × 104, the intensity of the turbulence in the near-wall region behind the front is large enough to induce the formation of a region containing streaks of low and high streamwise velocities. The streaks are present well into the buoyancy–inertia phase before the speed of the front decays below values at which the streaks can be sustained. The formation of the velocity streaks induces a streaky distribution of the bed friction velocity in the region immediately behind the front. This distribution becomes finer as the Grashof number increases. For simulations in which the only difference was the value of the Grashof number ($\sqrt {Gr}$ = 4.7 × 104 versus $\sqrt {Gr}$ = 106), analysis of the non-dimensional bed friction velocity distributions shows that the capacity of the gravity current to entrain sediment from the bed increases with the Grashof number. Past the later stages of the transition to the buoyancy–inertia phase, the temporal variations of the potential energy, the kinetic energy and the integral of the total dissipation rate are logarithmic.

Sign in / Sign up

Export Citation Format

Share Document