The behaviour of a gravity current propagating into a two-layer stratified ambient fluid is described in detail. A comprehensive description is given of the different flow regimes, with particular emphasis on the front condition linking the thickness of the gravity current to its speed of propagation and the transfer of energy to upstream disturbances in the form of internal bores and nonlinear solitary waves. Hydraulic theory analogous to that of two-layer flow over topography (Baines, J. Fluid Mech., vol. 146, 1984, pp. 127–167) is extended to the gravity current problem to classify frontal behaviour into the following regimes: Type I, subcritical currents; Type II, currents that generate upstream undular bores; Type III, currents that generate an upstream monotonic bore connected by a rarefaction; Type IV, supercritical fronts with a large-amplitude trapped solitary-wave-like disturbance; and Type V, supercritical gravity currents. Over 200 two-dimensional Boussinesq–Euler simulations spanning a range of gravity current properties demonstrate good agreement, for both the behavioural regime and the front condition , with hydraulic theory that extends original work by Rottman & Simpson (Q. J. R. Meteorol. Soc., vol. 115, 1989, pp. 941–963) to arbitrary ambient layer thickness, and uses an improved closure for the upstream bore that correctly predicts the behaviour in the limit of large bore amplitude. In addition, the energy balance is analysed, and it is shown that the energy transfer from the gravity current to upstream disturbances is significant, and consistent with the hydraulic theory. The results demonstrate a clear connection to the problem of upstream resonance in two-layer flow over topography, and have significant implications for interpreting field observations of nonlinear internal waves generated by atmospheric density currents and coastal river plumes.

/pdf/a-general-description-of-a-gravity-current-front-propagating-sx2dd55cwi.pdf

A general description of a gravity current front propagating in a two-layer stratified fluid

Three-dimensional 3-D Large eddy simulation (LES) has become a powerful tool to investigate evolution and structure of gravity currents, especially for cases (e.g., high Reynolds number flows, flows with massive separation) where 3-D Direct numerical simulation using non-dissipative viscous solvers is computationally too expensive. In this paper we briefly review some important results obtained based on high-resolution 3-D LES of bottom-propagating compositional Boussinesq currents in lock-exchange configurations. LES was used to provide a detailed description of the structure of the current, to discuss the role of the large-scale coherent structures, and to predict the evolution of the front velocity over the different stages of the current propagation. Three main types of lock-exchange flows are considered: (1) currents with a high volume of release (HVR) and a low volume of release (LVR) propagating in a channel with a smooth horizontal bed; (2) HVR and LVR currents propagating in a horizontal channel containing a porous layer; and (3) currents propagating in a horizontal channel containing an array of bottom obstacles (2-D dunes and ribs). The simulations are performed using non-dissipative numerical algorithms and sub-grid scale models that predict a zero eddy viscosity in regions where the turbulence is negligible. Experimental data is used to validate LES predictions. LES results show that in most cases the evolution of the front velocity is consistent with that predicted based on shallow-flow theory. LES flow fields are then used to estimate important quantities (e.g., bed friction velocity, sediment entrainment capacity) that are very difficult to obtain from experiments and to understand how the structure and evolution of the current change because of the additional drag induced by obstacles present within the channel or at the channel bed. The paper also discusses how the evolution and structure of the current change as the Reynolds number is increased to values that are relevant for gravity currents encountered in geosciences and environmental engineering applications.

LES of lock-exchange compositional gravity currents: a brief review of some recent results

Internal solitary waves are hump-shaped, large-amplitude waves that are physically analogous to surface waves except that they propagate within the fluid, along density steps that typically characterize the layered vertical structure of lakes, oceans and the atmosphere. As do surface waves, internal solitary waves may overturn and break, and the process is thought to provide a globally significant source of turbulent mixing and energy dissipation. Although commonly observed in geophysical fluids, the origins of internal solitary waves remain unclear. Here we report a rarely observed natural case of the birth of internal solitary waves from a frontally forced interfacial gravity current intruding into a two-layer and vertically sheared background environment. The results of the analysis carried out suggest that fronts may represent additional and unexpected sources of internal solitary waves in regions of lakes, oceans and atmospheres that are dynamically similar to the situation examined here in the Saguenay Fjord, Canada.

Generation of internal solitary waves by frontally forced intrusions in geophysical flows

We consider the propagation of a gravity current of density ρc at the bottom of a two-layer stratified ambient in a horizontal channel of height H, in the high-Reynolds number Boussinesq domain. The study emphasizes theoretical-analytical modeling, however, experimental and Navier-Stokes simulation data are also presented and their comparison with theory is discussed. The stratification parameters are S = (ρ1 − ρ2)/(ρc − ρ2) where ρ is the fluid density, and φ = h1R/H where h1R is the (unperturbed) ambient interface height. Here, 1 and 2 denote, respectively, the lower and upper layer and c denotes the gravity current. The reduced gravity is defined as g′ = (ρc/ρ2 − 1)g. Rigorous results are obtained for the steady-state analogue of the classical problem of Benjamin [J. Fluid Mech. 31, 209 (1968)]10.1017/S0022112068000133, in which the half-infinite gravity current has thickness h and speed U. We thereby demonstrate that the Froude number F=U/(g′h)1/2 is a function of a = h/H, S, and φ. In general, two so...

/pdf/gravity-currents-in-a-two-layer-stratified-ambient-the-bj2ar8i8kz.pdf

Gravity currents in a two-layer stratified ambient: The theory for the steady-state (front condition) and lock-released flows, and experimental confirmations

We develop a vorticity-based approach for modelling quasi-steady, supercritical gravity currents propagating into a finite-height channel with arbitrary density and velocity stratification. The model enforces the conservation of mass, horizontal and vertical momentum. In contrast to previous approaches, it does not rely on empirical, energybased closure assumptions. Instead, the effective energy loss of the flow can be calculated a posteriori. The present model results in the formulation of a second-order, nonlinear ordinary differential equation (ODE) that can be solved in a straightforward fashion to determine the gravity-current velocity, along with the downstream ambient velocity and density profiles. Comparisons between model predictions and direct numerical simulations (DNS) show excellent agreement. Furthermore, they indicate that, for high Reynolds numbers, the gravity-current height adjusts itself so as to maximize the loss of energy.

Gravity currents propagating into ambients with arbitrary shear and density stratification: vorticity-based modelling

http://absimage.aps.org/image/DFD10/MWS_DFD10-2010-020099.pdf

A.W. Tan

Papers

Gravity currents in two-layer stratified media