Axisymmetric, constantly supplied gravity currents at high Reynolds number

2011 ◽  
Vol 675 ◽  
pp. 540-551 ◽  
Author(s):  
ANJA C. SLIM ◽  
HERBERT E. HUPPERT

We consider theoretically the long-time evolution of axisymmetric, high Reynolds number, Boussinesq gravity currents supplied by a constant, small-area source of mass and radial momentum in a deep, quiescent ambient. We describe the gravity currents using a shallow-water model with a Froude number closure condition to incorporate ambient form drag at the front and present numerical and asymptotic solutions. The predicted profile consists of an expanding, radially decaying, steady interior that connects via a shock to a deeper, self-similar frontal boundary layer. Controlled by the balance of interior momentum flux and frontal buoyancy across the shock, the front advances as (g′sQ/r1/4s)4/154/5, where g′s is the reduced gravity of the source fluid, Q is the total volume flux, rs is the source radius and is time. A radial momentum source has no effect on this solution below a non-zero threshold value. Above this value, the (virtual) radius over which the flow becomes critical can be used to collapse the solution onto the subthreshold one. We also use a simple parameterization to incorporate the effect of interfacial entrainment, and show that the profile can be substantially modified, although the buoyancy profile and radial extent are less significantly impacted. Our predicted profiles and extents are in reasonable agreement with existing experiments.

2009 ◽  
Vol 28 (5) ◽  
pp. 651-659 ◽  
Author(s):  
J.J. Monaghan ◽  
C.A. Mériaux ◽  
H.E. Huppert ◽  
J.M. Monaghan

2015 ◽  
Vol 778 ◽  
pp. 552-585 ◽  
Author(s):  
M. M. Nasr-Azadani ◽  
E. Meiburg

An analytical vorticity-based model is introduced for steady-state inviscid Boussinesq gravity currents in sheared ambients. The model enforces the conservation of mass and horizontal and vertical momentum, and it does not require any empirical closure assumptions. As a function of the given gravity current height, upstream ambient shear and upstream ambient layer thicknesses, the model predicts the current velocity as well as the downstream ambient layer thicknesses and velocities. In particular, it predicts the existence of gravity currents with a thickness greater than half the channel height, which is confirmed by direct numerical simulation (DNS) results and by an analysis of the energy loss in the flow. For high-Reynolds-number gravity currents exhibiting Kelvin–Helmholtz instabilities along the current/ambient interface, the DNS simulations suggest that for a given shear magnitude, the current height adjusts itself such as to allow for maximum energy dissipation.


2003 ◽  
Vol 494 ◽  
pp. 255-274 ◽  
Author(s):  
MARK HALLWORTH ◽  
HERBERT E. HUPPERT ◽  
MARIUS UNGARISH

2011 ◽  
Vol 673 ◽  
pp. 574-602 ◽  
Author(s):  
THOMAS BONOMETTI ◽  
MARIUS UNGARISH ◽  
S. BALACHANDAR

The time-dependent behaviour of non-Boussinesq high-Reynolds-number density currents, released from a lock of height h0 and length x0 into a deep ambient and spreading over horizontal flat boundaries, is considered. We use two-dimensional Navier–Stokes simulations to cover: (i) a wide range of current-to-ambient density ratios, (ii) a range of length-to-height aspect ratios of the initial release within the lock (termed the lock aspect ratio λ = x0/h0) and (iii) the different phases of spreading, from the initial acceleration phase to the self-similar regimes. The Navier–Stokes results are compared with predictions of a one-layer shallow-water model. In particular, we derive novel insights on the influence of the lock aspect ratio (λ) on the shape and motion of the current. It is shown that for lock aspect ratios below a critical value (λcrit), the dynamics of the current is significantly influenced by λ. We conjecture that λcrit depends on two characteristic time scales, namely the time it takes for the receding perturbation created at the lock upon release to reflect back to the front, and the time of formation of the current head. A comparison of the two with space–time diagrams obtained from the Navier–Stokes simulations supports this conjecture. The non-Boussinesq effect is observed to be significant. While the critical lock aspect ratio (λcrit) is of order 1 for Boussinesq currents, its value decreases for heavy currents and increases significantly (up to about 20) for light currents. We present a simple analytical model which captures this trend, as well as the observation that for a light current the speed of propagation is proportional to λ1/4 when λ < λcrit.


2011 ◽  
Vol 668 ◽  
pp. 351-362 ◽  
Author(s):  
MARCELLO MELDI ◽  
PIERRE SAGAUT ◽  
DIDIER LUCOR

A stochastic eddy-damped quasi-normal Markovian (EDQNM) approach is used to investigate self-similar decaying isotropic turbulence at a high Reynolds number (400 ≤ Reλ ≤ 104). The realistic energy spectrum functional form recently proposed by Meyers & Menevau (Phys. Fluids, vol. 20, 2008, p. 065109) is generalized by considering some of the model constants as random parameters, since they escape measure in most experimental set-ups. The induced uncertainty on the solution is investigated, building response surfaces for decay power-law exponents of usual physical quantities. Large-scale uncertainties are considered, the emphasis being put on Saffman and Batchelor turbulences. The sensitivity of the solution to initial spectrum uncertainties is quantified through probability density functions of the decay exponents. It is observed that the initial spectrum shape at very large scales governs the long-time evolution, even at a high Reynolds number, a parameter which is not explicitly taken into account in many theoretical works. Therefore, a universal asymptotic behaviour in which kinetic energy decays as t−1 is not detected. However, this decay law is observed at finite Reynolds numbers with low probability for some initial conditions.


Sign in / Sign up

Export Citation Format

Share Document