Small-scale entrainment in inclined gravity currents

We investigate the effect of buoyancy on the small-scale aspects of turbulent entrainment by performing direct numerical simulation of a gravity current and a wall jet. In both flows, we detect the turbulent/nonturbulent interface separating turbulent from irrotational ambient flow regions using a range of enstrophy iso-levels spanning many orders of magnitude. Conform to expectation, the relative enstrophy isosurface velocity \(v_n\) in the viscous superlayer scales with the Kolmogorov velocity for both flow cases. We connect the integral entrainment coefficient E to the small-scale entrainment and observe excellent agreement between the two estimates throughout the viscous superlayer. The contribution of baroclinic torque to \(v_n\) is negligible, and we show that the primary reason for reduced entrainment in the gravity current as compared to the wall-jet are 1) the reduction of \(v_n\) relative to the integral velocity scale \(u_T\); and 2) the reduction in the surface area of the isosurfaces.     

The flow of high-Reynolds axisymmetric gravity currents of a stratified fluid into a stratified ambient: shallow-water and box model solutions

We consider the axisymmetric flow (in a full cylinder or a wedge) of high-Reynolds-number Boussinesq gravity currents and intrusions systems in which both the ambient and the propagating “current” are linearly stratified. The main focus is on a current of fixed volume released from a cylinder lock; the height ratio of the fluids H, and the stratification parameter of the ambient S, are quite general. We develop a one-layer shallow-water model. The internal stratification enters as a new dimensionless parameter, ${\sigma \in [0, 1]}$ . In general, the time-dependent motion is obtained by standard finite-difference solutions; a self-similar analytical solution exists for S?= 0. We show that, in general, the speed of propagation decreases when the internal stratification becomes more pronounced (σ increases). We also developed a box-model approximation, and show that the resulting radius of propagation is in good agreement with the more rigorous shallow-water prediction.     

Numerical simulation of particle-driven gravity currents

Particle-driven gravity currents frequently occur in nature, for instance as turbidity currents in reservoirs. They are produced by the buoyant forces between fluids of different density and can introduce sediments and pollutants into water bodies. In this study, the propagation dynamics of gravity currents is investigated using the FLOW-3D computational fluid dynamics code. The performance of the numerical model using two different turbulence closure schemes namely the renormalization group (RNG) ${k-\epsilon}$ scheme in a Reynold-averaged Navier-Stokes framework (RANS) and the large-eddy simulation (LES) technique using the Smagorinsky scheme, were compared with laboratory experiments. The numerical simulations focus on two different types of density flows from laboratory experiments namely: Intrusive Gravity Currents (IGC) and Particle-Driven Gravity Currents (PDGC). The simulated evolution profiles and propagation speeds are compared with laboratory experiments and analytical solutions. The numerical model shows good quantitative agreement for predicting the temporal and spatial evolution of intrusive gravity currents. In particular, the simulated propagation speeds are in excellent agreement with experimental results. The simulation results do not show any considerable discrepancies between RNG ${k-\epsilon}$ and LES closure schemes. The FLOW-3D model coupled with a particle dynamics algorithm successfully captured the decreasing propagation speeds of PDGC due to settling of sediment particles. The simulation results show that the ratio of transported to initial concentration C _o /C _i by the gravity current varies as a function of the particle diameter d _s . We classify the transport pattern by PDGC into three regimes: (1) a suspended regime (d _s is less than about 16 μm) where the effect of particle deposition rate on the propagation dynamics of gravity currents is negligible i.e. such flows behave like homogeneous fluids (IGC); (2) a mixed regime (16 μm < d _s <40 μm) where deposition rates significantly change the flow dynamics; and (3) a deposition regime (d _s ?> 40 μm) where the PDGC rapidly loses its forward momentum due to fast deposition. The present work highlights the potential of the RANS simulation technique using the RNG ${k-\epsilon}$ turbulence closure scheme for field scale investigation of particle-driven gravity currents.     

Lock-exchange gravity currents over rough bottoms

In nature, density driven currents often flow over or within a bottom roughness: a sea breeze encountering tall buildings, a shallow flow encountering aquatic vegetation, or a dense oceanic current flowing over a rough bottom. Laboratory experiments investigating the mechanisms by which bottom roughness enhances or inhibits entrainment and dilution in a lock-exchange dense gravity current have been conducted. The bottom roughness has been idealized by an array of vertical, rigid cylinders. Both spacing (sparse vs. dense configuration) and height of the roughness elements compared with the height of the current have been varied. Two-dimensional density fields have been obtained. Experimental results suggest that enhancement of the entrainment/dilution of the current can occur due to two different mechanisms. For a sparse configuration, the dense current propagates between the cylinders and the entrainment is enhanced by the vortices generated in the wake of the cylindrical obstacles. For a dense configuration, the dense current rides on top of the cylinders and the dilution is enhanced by the onset of convective instability between the dense current above the cylinders and the ambient lighter water between the cylinders. For low values of the ratio of the cylinder to lock height \(\lambda \) the dense current behavior approaches that of a current over a smooth bottom, while the largest deviations from the smooth bottom case are observed for large values of \(\lambda \).     

Dynamics of the head of gravity currents

The present work experimentally investigates the dynamics of unsteady gravity currents produced by lock-release of a saline mixture into a fresh water tank. Seven different experimental runs were performed by varying the density of the saline mixture in the lock and the bed roughness. Experiments were conducted in a Perspex flume, of horizontal bed and rectangular cross section, and recorded with a CCD camera. An image analysis technique was applied to visualize and characterize the current allowing thus the understanding of its general dynamics and, more specifically, of the current head dynamics. The temporal evolution of both head length and mass shows repeated stretching and breaking cycles: during the stretching phase, the head length and mass grow until reaching a limit, then the head becomes unstable and breaks. In the instants of break, the head aspect ratio shows a limit of 0.2 and the mass of the head is of the order of the initial mass in the lock. The average period of the herein called breaking events is seen to increase with bed roughness and the spatial periodicity of these events is seen to be approximately constant between runs. The rate of growth of the mass at the head is taken as a measure to assess entrainment and it is observed to occur at all stages of the current development. Entrainment rate at the head decreases in time suggesting this as a phenomenon ruled by local buoyancy and the similarity between runs shows independence from the initial reduced gravity and bed roughness.     

Lock-release gravity currents propagating over roughness elements

Environmental Fluid Mechanics - The influence of bottom roughness on the development of lock-release gravity currents is investigated through laboratory experiments using Particle Image...     

Experiments on gravity currents propagating on unbounded uniform slopes

Gravity currents propagating on \(12^\circ \), \(9^\circ \), \(6^\circ \), \(3^\circ \) unbounded uniform slopes and on an unbounded horizontal boundary are reported. Results show that there are two stages of the deceleration phase. In the early stage of the deceleration phase, the front location history follows \({(x_f+x_0)}^2 = {(K_I B)}^{1/2} (t+t_{I})\), where \((x_f+x_0)\) is the front location measured from the virtual origin, \(K_I\) an experimental constant, B the total buoyancy, t time and \(t_I\) the t-intercept. In the late stage of the deceleration phase for the gravity currents on \(12^\circ \), \(9^\circ \), \(6^\circ \) unbounded uniform slopes, the front location history follows \({(x_f+x_0)}^{8/3} = K_{VS} {{B}^{2/3} V^{2/9}_0 }{\nu }^{-1/3} ({t+t_{VS}})\), where \(K_{VS}\) is an experimental constant, \(V_0\) the initial volume of heavy fluid, \(\nu \) the kinematic viscosity and \(t_{VS}\) the t-intercept. In the late stage of the deceleration phase for the gravity currents on a \(3^\circ \) unbounded uniform slope and on an unbounded horizontal boundary, the front location history follows \({(x_f+x_0)}^{4} = K_{VM} {{B}^{2/3} V^{2/3}_0 }{\nu }^{-1/3} ({t+t_{VM}})\), where \(K_{VM}\) is an experimental constant and \(t_{VM}\) the t-intercept. Two qualitatively different flow morphologies are identified in the late stage of the deceleration phase. For the gravity currents on \(12^\circ \), \(9^\circ \), \(6^\circ \) unbounded uniform slopes, an ‘active’ head separates from the body of the current. For the gravity currents on a \(3^\circ \) unbounded uniform slope and on an unbounded horizontal boundary, the gravity currents maintain an integrated shape throughout the motion. Results indicate two possible routes to the final stage of the gravity currents on unbounded uniform slopes.

     

Mass-based depth and velocity scales for gravity currents and related flows

Gravity driven flows on inclines can be caused by cold, saline or turbid inflows into water bodies. Another example are cold downslope winds, which are caused by cooling of the atmosphere at the lower boundary. In a well-known contribution, Ellison and Turner (ET) investigated such flows by making use of earlier work on free shear flows by Morton, Taylor and Turner (MTT). Their entrainment relation is compared here with a spread relation based on a diffusion model for jets by Prandtl. This diffusion approach is suitable for forced plumes on an incline, but only when the channel topography is uniform, and the flow remains supercritical. A second aspect considered here is that the structure of ET’s entrainment relation, and their shallow water equations, agrees with the one for open channel flows, but their depth and velocity scales are those for free shear flows, and derived from the velocity field. Conversely, the depth of an open channel flow is the vertical extent of the excess mass of the liquid phase, and the average velocity is the (known) discharge divided by the depth. As an alternative to ET’s parameterization, two sets of flow scales similar to those of open channel flows are outlined for gravity currents in unstratified environments. The common feature of the two sets is that the velocity scale is derived by dividing the buoyancy flux by the excess pressure at the bottom. The difference between them is the way the volume flux is accounted for, which—unlike in open channel flows—generally increases in the streamwise direction. The relations between the three sets of scales are established here for gravity currents by allowing for a constant co-flow in the upper layer. The actual ratios of the three width, velocity, and buoyancy scales are evaluated from available experimental data on gravity currents, and from field data on katabatic winds. A corresponding study for free shear flows is referred to. Finally, a comparison of mass-based scales with a number of other flow scales is carried out for available data on a two-layer flow over an obstacle. Mass-based flow scales can also be used for other types of flows, such as self-aerated flows on spillways, water jets in air, or bubble plumes.     

Front propagation of gravity currents on inclined bottoms in linearly stratified fluids

Environmental Fluid Mechanics - Gravity currents propagating on an inclined bottom into stratified environment can be frequently encountered in nature or engineering fields. However, theoretical...     

On gravity currents confined by porous boundaries in containers of general cross-sections

Environmental Fluid Mechanics - We consider an inertial (large Reynolds number) gravity current (GC) released from a lock of length $$x_0$$ and height $$h_0$$ into an ambient fluid of height...     

Propagation,deposition, and suspension characteristics of constant-volume particle-driven gravity currents

Environmental Fluid Mechanics - In this laboratory study, propagation behaviour, particle deposition patterns, and suspension characteristics of non-cohesive particle-driven gravity currents formed...     

On the front conditions for gravity currents in channels of general cross-section

We consider the propagation of a high-Reynolds-number gravity current in a horizontal channel with general cross-section whose width is \(f(z), 0 \le z\le H\), and the gravity acceleration g acts in \(-z\) direction. (The classical rectangular cross-section is covered by the particular case \(f(z) =\) const.) We assume a two-layer system of homogeneous fluids of constant densities \(\rho _{c}\) (current, of height \(h < H \)) and smaller \(\rho _{a}\) (ambient, filling the remaining part of the channel). We focus attention on the calculation and assessment of the nose Froude-number condition \(Fr = U/(g' h)^{1/2}\); here U is the speed of propagation of the current and \(g' = (\rho _{c}/\rho _{a}-1) g\) is the reduced gravity. We first revisit the steady-state current, and derive compact insightful expressions of Fr and energy dissipation as a function of \(\varphi \) (\(=\) area fraction occupied by the current in the cross-section). We show that the head loss \(\delta _0\) on the stagnation line is formally a degree of freedom in the determination of \(Fr(\varphi )\), and we clarify the strong connections with the head loss \(\delta \) in the ambient fluid, and with the overall rate of dissipation \(\dot{{\mathcal{D}}}\). We demonstrate that the closure \(\delta _0 = 0\) [suggested by Benjamin (J Fluid Mech 31, 209–248, 1968) for the rectangular cross-section] produces in general the smallest Fr for a given \(\varphi \); the results are valid for a significant range \([0, \varphi _{\max }]\), in which the current is dissipative, except for the point \(\varphi _{\max }\) where \(\delta = \dot{{\mathcal{D}}} = 0\). We show that imposing the closure \(\delta = \dot{{\mathcal{D}}} = 0\), which corresponds to an energy-conserving or non-dissipative current, produces in general unacceptable restrictions of the range of validity, and large values of Fr; in particular, deep currents (\(\varphi < 0.3\) say) must be excluded because they are inherently dissipative. On the other hand, the compromise closure \(\delta (\varphi ) =\delta _0(\varphi )\) produces the simple \(Fr(\varphi ) = \sqrt{2}(1 - \varphi )\) formula whose values and dissipation properties are very close, and the range of validity is identical, to these obtained with Benjamin’s closure (moreover, we show that this corresponds to circulation-conservation solutions). The results are illustrated for practical cross-section geometries (rectangle, \(\Delta \) and \(\nabla \) triangle, circle, and the general power-law \(f(z) = b z ^\alpha \) (\(b>0, \alpha \ge 0, 0< z \le H\)). Next, we investigate the connection of the steady-state results with the time-dependent current, and show that in a lock-released current the rate of dissipation of the system is equal to, or larger than, that obtained for Fr corresponding to the conditions at the nose of the current. The results and insights of this study cover a wide range of cross-section geometry and apply to both Boussinesq and non-Boussinesq systems; they reveal a remarkable robustness of Fr as a function of \(\varphi \).     

Vortical structures,entrainment and mixing process in the lateral discharge of the gravity current

Environmental Fluid Mechanics - Lateral gravity currents can play a critical role in the exchange of materials between terrestrial and marine ecosystems. The three-dimensional flow structure and...     

A steady-state model for asymmetric intrusive gravity currents in a linearly stratified ambient

The behavior of the steady intrusive gravity current of thickness h and density ρ _c which propagates with speed U at the neutral buoyancy level of a long horizontal channel of height H into a stratified ambient fluid whose density increases linearly from ρ _o to ρ _b is investigated. The intrusive and the ambient fluids are assumed to be asymmetric with respect to the neutral-buoyancy level. The Boussinesq, high-Reynolds number two-dimensional configuration is considered. Long’s model combined with the flow-force balance over the width of the channel and the pressure balances over a density current are used to obtain the desired results. It is shown that the intrusion velocity decreases with decreasing the asymmetry of the system and approaches its minimum for the symmetric configuration (however, the difference of speed between asymmetric and symmetric configurations shows no significant differences).     

Thin-layer models for gravity currents in channels of general cross-section area,a review

We present a brief review of the recent investigations on gravity currents in horizontal channels with non-rectangular cross-section area (such as triangle, \(\bigvee \)-valley, circle/semi-circle, trapezoid) which occur in nature (e.g., rivers) and constructed environment (tunnels, reservoirs, canals). To be specific, we discuss the propagation of a gravity current (GC) in a horizontal channel along the horizontal coordinate x, with gravity g acting in the \(-z\) direction, and y the horizontal–lateral coordinate. The bottom and top of the channel are at \(z=0,H\). The “standard” problem is concerned with 2D flow in a channel with rectangular (or laterally unbounded) cross-section area (CSA). Recent investigations have successfully extended the standard knowledge to the channels of CSA given by the quite general \(-f_1(z)\le y \le f_2(z)\) for \(0 \le z \le H\). This includes the practical \(\bigvee \)-valley, triangle, circle/semi-circle and trapezoid; these geometries may be in “up” or “down” setting with respect to gravity, e.g., \(\bigtriangleup \) and \(\bigtriangledown \). The major objective of the extended theory is to predict the height of the interface \(z=h(x,t)\) and the velocity (averaged over the CSA) u(x, t), where t is time; the prediction includes the speed and position of the nose \(u_N(t), x_N(t)\). We show that the motion is governed by a set of simplified equations, called “model,” that provides versatile and insightful solutions and trends. The emphasis in on a high-Reynolds-number current whose motion is dominated by buoyancy–inertia balance; in particular a GC released from a lock, which also contains general effects such as front and internal jumps (shocks), and reflected bore. We discuss two-layer, one-layer, and box models; Boussinesq and non-Boussinesq systems; compositional and particle-driven cases; and the effect of stratification of the ambient fluid. The models are self-contained, and admit realistic initial and boundary conditions. The governing equations are amenable to analytical solutions in some special circumstances. Some salient features of the buoyancy-viscous regime, and the estimate for the length at which transition to this regime takes place, are also presented. Some experimental support to the theory, and open questions for further investigations, are also mentioned. The major conclusions are (1) The CSA geometry has significant influence on the motion of the GC; and (2) The new theory is a useful, very significant, extension of the standard two-dimensional GC problem. The standard current is just a particular case, \(f_{1,2} =\) constants, among many other covered by the new theory.