This paper presents the investigation of gravity currents by both laboratory experiments and a mathematical model. Eleven lock-exchange experiments, in which lock position, the initial current height, and density varied, were carried out to test the model validity and to compare laboratory results with previous expressions found in the literature. A two-layer shallow-water model was used to simulate all the runs. This model is new if compared with previous shallow-water models used to simulate gravity currents, because it accounts for both the entrainment and the free surface. A modified Turner’s formula is used to model the entrainment between the two fluids. The developed shallow-water models with and without entrainment are also compared, showing a better agreement when mixing is accounted for. Also, the effect of the free surface is shown by comparing the developed two-layer shallow-water model with a free surface and two different single-layer models with a rigid-lid approximation. Laboratory experim...

Entrainment and mixing in lock-exchange gravity currents are investigated by large eddy simulations. Nine cases are analysed, varying the initial excess density driving the motion and the aspect ratio r of the initial water depth to the lock length. Laboratory experiments are also performed and a fair agreement between numerical simulations and measurements is found. Mixing between the gravity current and the ambient fluid, in both the slumping and self-similar phases, is studied for a range of entrainment parameters, gravity current fractional area and using an energy budget method. The entrainment is found to increase as r decreases. The occurrence of irreversible mixing is detected during the entire development of the flow, i.e. both in the slumping and self-similar phases. A higher amount of mixing is observed as r decreases and the initial excess density increases.

Dense overflows entrain surrounding waters at specific locations, for example, sills and constrictions, but also along the descent over the continental slope. The amount of entrainment dictates the final properties of these overflows, and thus is of fundamental importance to the understanding of the formation of deep water masses. Even when resolving the overflows, coarse resolution global circulation and climate models cannot resolve the entrainment processes that are often parameterized. A new empirical parameterization is suggested, obtained using an oceanic and laboratory dataset, which includes two novel aspects. First, the parameterization depends on both the Froude number (Fr) and Reynolds number of the flow. Second, it takes into account subcritical (Fr < 1) entrainment. A weak, but nonzero, entrainment can change the final density and, consequently, the depth and location of important water masses in the open ocean. This is especially true when the dense current follows a long path over the slope in a subcritical regime, as observed in the southern Greenland Deep Western Boundary Current. A streamtube model employing this new parameterization gives results that are more consistent with previous laboratory and oceanographic observations than when a classical parameterization is used. Finally, the new parameterization predictions compare favorably to recent oceanographic measurements of entrainment and turbulent diapycnal mixing rates, using scaling arguments to relate the entrainment ratio to diapycnal diffusivities.

We discuss laboratory experiments investigating mixing in a density-driven current flowing down a sloping bottom, in a rotating homogenous fluid. A systematic study spanning a wide range of Froude, 0.8 < Fr < 10, and Reynolds, 10 < Re < 1400, numbers was conducted by varying three parameters: the bottom slope; the flow rate; and the density of the dense fluid. Different flow regimes were observed, i.e. waves (non-breaking and breaking) and turbulent regimes, while changing the above parameters. Mixing in the density-driven current has been quantified within the observed regimes, and at different locations on the slope. The dependence of mixing on the relevant non-dimensional numbers, i.e. slope, Fr and Re, is discussed. The entrainment parameter, E, was found to be dependent not only on Fr, as assumed in previous studies, but also on Re. In particular, mixing increased with increasing Fr and Re. For low Fr and Re, the magnitude of the mixing was comparable to mixing in the ocean. For large Fr and Re, mixing was comparable to that observed in previous laboratory experiments that exhibited the classic turbulent entrainment behaviour.

The dynamics of a three-dimensional gravity current is investigated by both laboratory experiments and numerical simulations. The experiments take place in a rectangular tank, which is divided into two square reservoirs with a wall containing a sliding gate of width b. The two reservoirs are filled to the same height H, one with salt water and the other with fresh water. The gravity current starts its evolution as soon as the sliding gate is manually opened. Experiments are conducted with either smooth or rough surface on the bottom of the tank. The bottom roughness is created by gluing sediment material of different diameters to the surface. Five diameter values for the surface roughness and two salinity conditions for the fluid are investigated. The mathematical model is based on shallow-water theory together with the single-layer approximation, so that the model is strictly hyperbolic and can be put into conservative form. Consequently, a finite-volume-based numerical algorithm can be applied. The Godunov formulation is used together with Roe’s approximate Riemann solver. Comparisons between the numerical and experimental results show satisfactory agreement. The behavior of the gravity current is quite unusual and cannot be interpreted using the usual model framework adopted for two-dimensional and axisymmetric gravity currents. Two main phases are apparent in the gravity current evolution; during the first phase the front velocity increases, and during the second phase the front velocity decreases and the dimensionless results, relative to the different densities, collapse onto the same curve. A systematic discrepancy is seen between the numerical and experimental results, mainly during the first phase of the gravity current evolution. This discrepancy is attributed to the limits of the mathematical formulation, in particular, the neglect of entrainment in the mathematical model. An interesting result arises from the influence of the bottom surface roughness; it both reduces the front velocity during the second phase of motion and attenuates the differences between the experimental and numerical front velocities during the first phase of motion.