On the propagation of gravity currents over and through a submerged array of circular cylinders

2017 ◽  
Vol 831 ◽  
pp. 394-417 ◽  
Author(s):  
Jian Zhou ◽  
Claudia Cenedese ◽  
Tim Williams ◽  
Megan Ball ◽  
Subhas K. Venayagamoorthy ◽  
...  
Keyword(s):  
Volume Fraction ◽  
Gravity Current ◽  
Gravity Currents ◽  
Front Velocity ◽  
Circular Cylinders ◽  
Dimensional Parameter ◽  
Flow Interruption ◽  
Line Array ◽  
Through Flow ◽  
Array Density

The propagation of full-depth lock-exchange bottom gravity currents past a submerged array of circular cylinders is investigated using laboratory experiments and large eddy simulations. Firstly, to investigate the front velocity of gravity currents across the whole range of array density $\unicode[STIX]{x1D719}$ (i.e. the volume fraction of solids), the array is densified from a flat bed ($\unicode[STIX]{x1D719}=0$) towards a solid slab ($\unicode[STIX]{x1D719}=1$) under a particular submergence ratio $H/h$, where $H$ is the flow depth and $h$ is the array height. The time-averaged front velocity in the slumping phase of the gravity current is found to first decrease and then increase with increasing $\unicode[STIX]{x1D719}$. Next, a new geometrical framework consisting of a streamwise array density $\unicode[STIX]{x1D707}_{x}=d/s_{x}$ and a spanwise array density $\unicode[STIX]{x1D707}_{y}=d/s_{y}$ is proposed to account for organized but non-equidistant arrays ($\unicode[STIX]{x1D707}_{x}\neq \unicode[STIX]{x1D707}_{y}$), where $s_{x}$ and $s_{y}$ are the streamwise and spanwise cylinder spacings, respectively, and $d$ is the cylinder diameter. It is argued that this two-dimensional parameter space can provide a more quantitative and unambiguous description of the current–array interaction compared with the array density given by $\unicode[STIX]{x1D719}=(\unicode[STIX]{x03C0}/4)\unicode[STIX]{x1D707}_{x}\unicode[STIX]{x1D707}_{y}$. Both in-line and staggered arrays are investigated. Four dynamically different flow regimes are identified: (i) through-flow propagating in the array interior subject to individual cylinder wakes ($\unicode[STIX]{x1D707}_{x}$: small for in-line array and arbitrary for staggered array; $\unicode[STIX]{x1D707}_{y}$: small); (ii) over-flow propagating on the top of the array subject to vertical convective instability ($\unicode[STIX]{x1D707}_{x}$: large; $\unicode[STIX]{x1D707}_{y}$: large); (iii) plunging-flow climbing sparse close-to-impermeable rows of cylinders with minor streamwise intrusion ($\unicode[STIX]{x1D707}_{x}$: small; $\unicode[STIX]{x1D707}_{y}$: large); and (iv) skimming-flow channelized by an in-line array into several subcurrents with strong wake sheltering ($\unicode[STIX]{x1D707}_{x}$: large; $\unicode[STIX]{x1D707}_{y}$: small). The most remarkable difference between in-line and staggered arrays is the non-existence of skimming-flow in the latter due to the flow interruption by the offset rows. Our analysis reveals that as $\unicode[STIX]{x1D719}$ increases, the change of flow regime from through-flow towards over- or skimming-flow is responsible for increasing the gravity current front velocity.

Front dynamics and entrainment of finite circular gravity currents on an unbounded uniform slope

2016 ◽  
Vol 801 ◽  
pp. 322-352 ◽  
Author(s):  
N. Zgheib ◽  
A. Ooi ◽  
S. Balachandar
Keyword(s):  
Froude Number ◽  
Temporal Evolution ◽  
Gravity Current ◽  
Gravity Currents ◽  
Front Velocity ◽  
Acceleration Phase ◽  
Simple Method ◽  
Entrainment Coefficient ◽  
Converging Flow ◽  
Front Dynamics

We report on the dynamics of circular finite-release Boussinesq gravity currents on a uniform slope. The study comprises a series of highly resolved direct numerical simulations for a range of slope angles between $5^{\circ }$ and $20^{\circ }$. The simulations were fixed at Reynolds number $Re=5000$ for all slopes considered. The temporal evolution of the front is compared to available experimental data. One of the interesting aspects of this study is the detection of a converging flow towards the centre of the gravity current. This converging flow is a result of the finite volume of the release coupled with the presence of a sloping boundary, which results in a second acceleration phase in the front velocity of the current. The details of the dynamics of this second acceleration and the redistribution of material in the current leading to its development will be discussed. These finite-release currents are invariably dominated by the head where most of the mixing and ambient entrainment occurs. We propose a simple method for defining the head of the current from which we extract various properties including the front Froude number and entrainment coefficient. The Froude number is seen to increase with steeper slopes, whereas the entrainment coefficient is observed to be weakly dependent on the bottom slope.

scholarly journals Buoyant gravity currents along a sloping bottom in a rotating fluid

2002 ◽  
Vol 464 ◽  
pp. 251-278 ◽  
Author(s):  
STEVEN J. LENTZ ◽  
KARL R. HELFRICH
Keyword(s):  
Vertical Wall ◽  
Gravity Current ◽  
Classical Problem ◽  
Propagation Speed ◽  
Gravity Currents ◽  
Scaling Theory ◽  
Bottom Slope ◽  
Dimensional Parameter ◽  
Density Anomaly ◽  
Sloping Bottom

The dynamics of buoyant gravity currents in a rotating reference frame is a classical problem relevant to geophysical applications such as river water entering the ocean. However, existing scaling theories are limited to currents propagating along a vertical wall, a situation almost never realized in the ocean. A scaling theory is proposed for the structure (width and depth), nose speed and flow field characteristics of buoyant gravity currents over a sloping bottom as functions of the gravity current transport Q, density anomaly g′, Coriolis frequency f, and bottom slope α. The nose propagation speed is cp ∼ cw/ (1 + cw/cα) and the width of the buoyant gravity current is Wp ∼ cw/ f(1 + cw/cα), where cw = (2Qg′ f)1/4 is the nose propagation speed in the vertical wall limit (steep bottom slope) and cα = αg/f is the nose propagation speed in the slope-controlled limit (small bottom slope). The key non-dimensional parameter is cw/cα, which indicates whether the bottom slope is steep enough to be considered a vertical wall (cw/cα → 0) or approaches the slope-controlled limit (cw/cα → ∞). The scaling theory compares well against a new set of laboratory experiments which span steep to gentle bottom slopes (cw/cα = 0.11–13.1). Additionally, previous laboratory and numerical model results are reanalysed and shown to support the proposed scaling theory.

Lock-exchange gravity currents propagating in a channel containing an array of obstacles

2015 ◽  
Vol 765 ◽  
pp. 544-575 ◽  
Author(s):  
Ayse Yuksel Ozan ◽  
George Constantinescu ◽  
Andrew J. Hogg
Keyword(s):  
Reynolds Number ◽  
Numerical Simulations ◽  
Drag Coefficient ◽  
Gravity Current ◽  
Gravity Currents ◽  
Front Velocity ◽  
Release Time ◽  
Exchange Flow ◽  
Layer Model ◽  
Shallow Layer

AbstractLarge eddy simulation (LES) is used to investigate the evolution of Boussinesq gravity currents propagating through a channel of height $H$ containing a staggered array of identical cylinders of square cross-section and edge length $D$. The cylinders are positioned with their axes horizontal and perpendicular to the (streamwise) direction along which the lock-exchange flow develops. The effects of the volume fraction of solids, ${\it\phi}$, the Reynolds number and geometrical parameters describing the array of obstacles on the structure of the lock-exchange flow, total drag force acting on the gravity current, front velocity and global energy budget are analysed. Simulation results show that the currents rapidly transition to a state in which the extra resistance provided by the cylinders strongly retards the motion and dominates the dissipative processes. A shallow layer model is also formulated and similarity solutions for the motion are found in the regime where the driving buoyancy forces are balanced by the drag arising from the interaction with the cylinders. The numerical simulations and this shallow layer model show that low-Reynolds-number currents transition to a drag-dominated regime in which the resistance is linearly proportional to the flow speed and, consequently, the front velocity, $U_{f}$, is proportional to $t^{-1/2}$, where $t$ is the time measured starting at the gate release time. By contrast, high-Reynolds-number currents, for which the cylinder Reynolds number is sufficiently high that the drag coefficient for most of the cylinders can be considered constant, transition first to a quadratic drag-dominated regime in which the front speed determined from the simulations is given by $U_{f}\sim t^{-0.25}$, before undergoing a subsequent transition to the aforementioned linear drag regime in which $U_{f}\sim t^{-1/2}$. Meanwhile, away from the front, the depth-averaged gravity current velocity is proportional to $t^{-1/3}$, a result that is in agreement with the shallow water model. It is suggested that the difference between these two is due to mixing processes, which are shown to be significant in the numerical simulations, especially close to the front of the motion. Direct estimation of the drag coefficient $C_{D}$ from the numerical simulations shows that the combined drag parameter for the porous medium, ${\it\Gamma}_{D}=C_{D}{\it\phi}(H/D)/(1-{\it\phi})$, is the key dimensionless grouping of variables that determines the speed of propagation of the current within arrays with different $C_{D},{\it\phi}$ and $D/H$.

Large Eddy Simulation of Gravity Current Flow Past a Circular Cylinder Using Immersed Boundary Method

2015 ◽  
Author(s):  
Jae Hwan Wang ◽  
Hyun Sik Yoon
Keyword(s):  
Circular Cylinder ◽  
Large Eddy Simulation ◽  
Gravity Current ◽  
Numerical Approach ◽  
Gravity Currents ◽  
Hydrodynamic Forces ◽  
Circular Cylinders ◽  
Immersed Boundary ◽  
Eddy Simulation ◽  
Large Eddy

We investigated the flow of gravity currents past circular cylinders above the smooth bed. In order to simulate the gravity current interacting with cylinder, large eddy simulation (LES) coupled with a direct forcing/fictitious domain (DF/FD) method was employed. The hydrodynamic forces induced by the gravity current on a circular cylinder and the flow features were investigated according to a gap distance from location below the cylinder to the bed and Reynolds number. When the gravity current encounters the circular cylinder, the maximum drag force occurs regardless of gap distance. In addition, Von Karman vortex shedding emerges behind the circular cylinder for larger gap distance. In this regard, the results can be divided into impact, transient and quasi-steady stages based on the characteristics of the hydrodynamic forces varying with time, which is consistent with previous observations. Consequently, our numerical approach well simulated the flow of gravity currents past circular cylinders in good agreement with those of previous studies.

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.

scholarly journals Front Velocity and Front Location of Lock-Exchange Gravity Currents Descending a Slope in a Linearly Stratified Environment

2018 ◽  
Vol 144 (11) ◽  
pp. 04018068 ◽  
Author(s):  
Liang Zhao ◽  
Zhiguo He ◽  
Yafei Lv ◽  
Ying-Tien Lin ◽  
Peng Hu ◽  
...  
Keyword(s):  
Gravity Currents ◽  
Front Velocity

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.

scholarly journals Entrainment in Shallow Rotating Gravity Currents: A Modeling Study

2010 ◽  
Vol 40 (8) ◽  
pp. 1819-1834 ◽  
Author(s):  
Lars Umlauf ◽  
Lars Arneborg ◽  
Richard Hofmeister ◽  
Hans Burchard
Keyword(s):  
Froude Number ◽  
Gravity Current ◽  
Vertical Mixing ◽  
Spatial Separation ◽  
Gravity Currents ◽  
Transverse Jet ◽  
Entrainment Process ◽  
Set Up ◽  
Rotation Set ◽  
Bulk Parameters

Abstract The physics of shallow gravity currents passing through a rotating channel at subcritical Froude number is investigated here with a series of idealized numerical experiments. It is found that the combined effects of friction and rotation set up a complex transverse circulation that has some crucial implications for the entrainment process. A key component of this secondary circulation is a geostrophically balanced transverse jet in the interface that laterally drains fluid from the interface. This effect is shown to result in a strong cross-channel asymmetry and a spatial separation of the entrainment process: drained interfacial fluid is partly replaced by entrained ambient fluid on the deep side of the gravity current, whereas the downward mixing of buoyant fluid occurs on the shallow side. These results, closely corresponding to recent measurements in a shallow, channelized gravity current in the western Baltic Sea, illustrate that the description of entrainment as a strictly vertical mixing process with the help of local bulk parameters like the Froude number is not generally applicable in rotating gravity currents.

scholarly journals Axisymmetric viscous gravity currents flowing over a porous medium

2009 ◽  
Vol 622 ◽  
pp. 135-144 ◽  
Author(s):  
MELISSA J. SPANNUTH ◽  
JEROME A. NEUFELD ◽  
J. S. WETTLAUFER ◽  
M. GRAE WORSTER
Keyword(s):  
Porous Medium ◽  
Steady State ◽  
Gravity Current ◽  
Gravity Currents ◽  
Constant Input ◽  
Radial Extent ◽  
Vertically Aligned ◽  
Analytical Expressions ◽  
Permeability Structure ◽  
Good Agreement

We study the axisymmetric propagation of a viscous gravity current over a deep porous medium into which it also drains. A model for the propagation and drainage of the current is developed and solved numerically in the case of constant input from a point source. In this case, a steady state is possible in which drainage balances the input, and we present analytical expressions for the resulting steady profile and radial extent. We demonstrate good agreement between our experiments, which use a bed of vertically aligned tubes as the porous medium, and the theoretically predicted evolution and steady state. However, analogous experiments using glass beads as the porous medium exhibit a variety of unexpected behaviours, including overshoot of the steady-state radius and subsequent retreat, thus highlighting the importance of the porous medium geometry and permeability structure in these systems.

Sign in / Sign up

Export Citation Format

Share Document