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.

Entrainment and Transport in Idealized Three-Dimensional Gravity Current Simulation

2006 ◽  
Vol 23 (9) ◽  
pp. 1249-1269 ◽  
Author(s):  
Yu-Heng Tseng ◽  
David E. Dietrich
Keyword(s):  
Gravity Current ◽  
Vertical Mixing ◽  
Three Dimensional ◽  
Ocean Model ◽  
Numerical Dissipation ◽  
Good Convergence ◽  
Low Viscosity ◽  
Denmark Strait ◽  
Dimensional Gravity ◽  
Fine Resolution

Abstract A purely z-coordinate Dietrich/Center for Air Sea Technology (DieCAST) ocean model is applied to the Dynamics of Overflow Mixing and Entrainment (DOME) idealized bottom density current problem that is patterned after the Denmark Strait. The numerical results show that the background viscosity plays a more important role than the chosen coordinate system in the entrainment and mixing if the background viscosity is not small enough. Both higher horizontal viscosity and coarser resolution leads to slower along-slope propagation. Reducing vertical mixing parameterization also leads to slower along-slope propagation with thicker plume size vertically. The simulation gives consistent results for the moderate- and fine-resolution runs. At a very coarse grid the dense water descends more slowly and is mainly dominated by diffusion. Time-averaged downstream transport and entrainment are not very sensitive to viscosity after the flow reaches its quasi-steady status. However, more realistic eddies and flow structures are found in low-viscosity runs. The results show good convergence of the resolved flow as expected and clarify the effects of numerical dissipation/mixing on overflow modeling. Larger numerical dissipation is not required nor recommended in z-coordinate models.

scholarly journals Axisymmetric three-dimensional gravity currents generated by lock exchange

2018 ◽  
Vol 851 ◽  
pp. 507-544 ◽  
Author(s):  
Roberto Inghilesi ◽  
Claudia Adduce ◽  
Valentina Lombardi ◽  
Federico Roman ◽  
Vincenzo Armenio
Keyword(s):  
Froude Number ◽  
Flow Rate ◽  
Rapid Development ◽  
Three Dimensional ◽  
Gravity Currents ◽  
Constant Flow ◽  
Axially Symmetric ◽  
Grashof Number ◽  
Constant Flow Rate ◽  
Dimensional Gravity

Unconfined three-dimensional gravity currents generated by lock exchange using a small dividing gate in a sufficiently large tank are investigated by means of large eddy simulations under the Boussinesq approximation, with Grashof numbers varying over five orders of magnitudes. The study shows that, after an initial transient, the flow can be separated into an axisymmetric expansion and a globally translating motion. In particular, the circular frontline spreads like a constant-flow-rate, axially symmetric gravity current about a virtual source translating along the symmetry axis. The flow is characterised by the presence of lobe and cleft instabilities and hydrodynamic shocks. Depending on the Grashof number, the shocks can either be isolated or produced continuously. In the latter case a typical ring structure is visible in the density and velocity fields. The analysis of the frontal spreading of the axisymmetric part of the current indicates the presence of three regimes, namely, a slumping phase, an inertial–buoyancy equilibrium regime and a viscous–buoyancy equilibrium regime. The viscous–buoyancy phase is in good agreement with the model of Huppert (J. Fluid Mech., vol. 121, 1982, pp. 43–58), while the inertial phase is consistent with the experiments of Britter (Atmos. Environ., vol. 13, 1979, pp. 1241–1247), conducted for purely axially symmetric, constant inflow, gravity currents. The adoption of the slumping model of Huppert & Simpson (J. Fluid Mech., vol. 99 (04), 1980, pp. 785–799), which is here extended to the case of constant-flow-rate cylindrical currents, allows reconciling of the different theories about the initial radial spreading in the context of different asymptotic regimes. As expected, the slumping phase is governed by the Froude number at the lock’s gate, whereas the transition to the viscous phase depends on both the Froude number at the gate and the Grashof number. The identification of the inertial–buoyancy regime in the presence of hydrodynamic shocks for this class of flows is important, due to the lack of analytical solutions for the similarity problem in the framework of shallow water theory. This fact has considerably slowed the research on variable-flow-rate axisymmetric gravity currents, as opposed to the rapid development of the knowledge about cylindrical constant-volume and planar gravity currents, despite their own environmental relevance.

Three-dimensional gravity-current flow within a subaqueous bend: Spatial evolution and force balance variations

Sedimentology ◽  
2013 ◽  
Vol 60 (7) ◽  
pp. 1668-1680 ◽  
Author(s):  
Taoyuan Wei ◽  
Jeff Peakall ◽  
Daniel R. Parsons ◽  
Zhongyuan Chen ◽  
Baocheng Zhao ◽  
...  
Keyword(s):  
Current Flow ◽  
Gravity Current ◽  
Three Dimensional ◽  
Force Balance ◽  
Spatial Evolution ◽  
Dimensional Gravity

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.

On the interaction of an internal wavepacket with its induced mean flow and the role of streaming

2018 ◽  
Vol 838 ◽  
Author(s):  
Boyu Fan ◽  
T. Kataoka ◽  
T. R. Akylas
Keyword(s):  
Nonlinear Interaction ◽  
Laboratory Experiments ◽  
Modulation Instability ◽  
Three Dimensional ◽  
Mean Flow ◽  
Interaction Dynamics ◽  
Dimensional Gravity ◽  
Combined Action ◽  
Beam Bending

The coupled nonlinear interaction of three-dimensional gravity–inertia internal wavepackets, in the form of beams with nearly monochromatic profile, with their induced mean flow is discussed. Unlike general three-dimensional wavepackets, such modulated nearly monochromatic beams are not susceptible to modulation instability from their inviscid, purely modulation-induced mean flow. However, streaming – the induced mean flow associated with the production of mean potential vorticity via the combined action of dissipation and nonlinearity – can cause cross-beam bending, transverse broadening and increased along-beam decay of the beam profile, in qualitative agreement with earlier laboratory experiments. For wavepackets with general three-dimensional modulations, by contrast, streaming does arise, but plays a less prominent role in the interaction dynamics.

On the formation and propagation of nonlinear internal boluses across a shelf break

2007 ◽  
Vol 577 ◽  
pp. 137-159 ◽  
Author(s):  
SUBHAS K. VENAYAGAMOORTHY ◽  
OLIVER B. FRINGER
Keyword(s):  
High Resolution ◽  
Gravity Waves ◽  
Three Dimensional ◽  
Internal Gravity Waves ◽  
Gravity Currents ◽  
Interaction Process ◽  
Shelf Break ◽  
Ambient Fluid ◽  
Dense Fluid ◽  
Incident Waves

High-resolution two- and three-dimensional numerical simulations are performed of first-mode internal gravity waves interacting with a shelf break in a linearly stratified fluid. The interaction of nonlinear incident waves with the shelf break results in the formation of upslope-surging vortex cores of dense fluid (referred to here as internal boluses) that propagate onto the shelf. This paper primarily focuses on understanding the dynamics of the interaction process with particular emphasis on the formation, structure and propagation of internal boluses onshelf. A possible mechanism is identified for the excitation of vortex cores that are lifted over the shelf break, from where (from the simplest viewpoint) they essentially propagate as gravity currents into a linearly stratified ambient fluid.

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.

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.

Gravity currents impinging on bottom-mounted square cylinders: flow fields and associated forces

2009 ◽  
Vol 631 ◽  
pp. 65-102 ◽  
Author(s):  
E. GONZALEZ-JUEZ ◽  
E. MEIBURG ◽  
G. CONSTANTINESCU
Keyword(s):  
Drag Coefficient ◽  
Flow Structure ◽  
Gravity Current ◽  
Three Dimensional ◽  
Gravity Currents ◽  
Constant Density ◽  
Two Dimensional ◽  
Drag And Lift ◽  
The Impact ◽  
Three Dimensional Simulations

The unsteady drag and lift generated by the interaction of a gravity current with a bottom-mounted square cylinder are investigated by means of high-resolution Navier–Stokes simulations. Two-dimensional simulations for Reynolds numbers (Re) O(1000) and three-dimensional simulations for Re = O(10000) demonstrate that the drag coefficient increases exponentially towards a maximum as the current meets the cylinder, then undergoes strong fluctuations and eventually approaches a quasi-steady value. The simulation results show that the maximum drag coefficient can reach a value of 3, with the quasi-steady value being O(1), which should aid in selecting a design drag coefficient for submarine structures under the potential impact of gravity currents. The transient drag and lift fluctuations after impact are associated with the Kelvin–Helmholtz vortices in the mixing layer between the gravity current and the ambient fluid. As these vortices pass over the cylinder, they cause the convection of separated flow regions along the bottom wall towards the cylinder. In two-dimensional simulations at Re = O(10000), these flow structures are seen to be unrealistically coherent and to persist throughout the interaction, thus resulting in a noticeable overprediction of the drag and lift fluctuations. On the other hand, the impact of the current on the cylinder is seen to be very well captured by two-dimensional simulations at all Re values. Three-dimensional simulations lead to excellent agreement with available experimental data throughout the flow/structure interaction. They show that the spanwise variation of the drag is determined by the gravity current's lobe-and-cleft structure at impact and by an unsteady cellular flow structure similar to that found in constant-density flows at later times. A comparison between gravity-current flows and corresponding constant-density flows shows the hydrostatic drag component to be important for gravity currents.

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.

Sign in / Sign up

Export Citation Format

Share Document