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.

A wall-layer model is proposed for large eddy simulation of high Reynolds number turbulent flows in conjunction with immersed boundaries. The model is based on two main steps: the reconstruction of the velocity field at the first grid point off the immersed body and the modelization of the actual wall shear stress at the immersed boundary through imposition of a Reynolds averaged Navier–Stokes-like eddy viscosity obtained by means of analytical considerations. The model is tested in a turbulent plane channel flow with walls reproduced by immersed boundaries considering both Cartesian and curvilinear grids. Even with coarse and distorted grids the proposed methodology is able to reproduce accurately both first- and second-order turbulent statistics.

Lock-exchange gravity currents propagating up a slope are investigated by large eddy simulations, focusing on the entrainment and mixing processes occurring between the dense current and the ambient fluid. Relevant parameters, such as the aspect ratio of the initial volume of dense fluid in the lock R, the angle between the bottom boundary and the horizontal direction θ and the depth aspect ratio ϕ, are varied. The numerical results are compared with laboratory experiments and a good agreement is found. Entrainment and mixing in a lock-release gravity current are studied using different entrainment parameters and an energy budget method. The entrainment is found to depend on both Froude, Fr, and Reynolds, Re, numbers. In addition, the dependence of both entrainment and mixing on the parameters varied is discussed. The entrainment decreases with increasing steepness of the bottom and R. Irreversible mixing is not affected by the varied parameters during the slumping phase, while during the successive phases of motion, it is found to decrease with the increase of θ and R. Low entrainment and mixing occur for ϕ < 1.

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.