Inertial gravity current produced by the drainage of a cylindrical reservoir from an outer or inner edge

Abstract

We consider the time-dependent flow of a fluid of density $\unicode[STIX]{x1D70C}_{1}$ in a vertical cylindrical container embedded in a fluid of density $\unicode[STIX]{x1D70C}_{2}~({<}\unicode[STIX]{x1D70C}_{1})$ whose side boundary is suddenly removed and the fluid drains freely from the edge. We show that in the inertial–buoyancy regime (large initial Reynolds number) the flow is modelled by the shallow-water equations and bears similarities to a gravity current released from a lock (the dam-break problem) driven by the reduced gravity $g^{\prime }=(1-\unicode[STIX]{x1D70C}_{2}/\unicode[STIX]{x1D70C}_{1})g$. This formulation is amenable to an efficient finite-difference solution. Moreover, we demonstrate that similarity solutions exist, and show that the flow created by the dam break approaches the predicted self-similar behaviour when the volume ratio ${\mathcal{V}}(t)/{\mathcal{V}}(0)\approx 1/2$ where $t$ is time elapsed from the dam break. We considered two cases of drainage: (i) outward from the outer boundary in a full-radius reservoir; and (ii) inward from the inner radius in an annular-shaped reservoir. For the first case the similarity solution is expressed analytically, while the second case is more complicated and requires a numerical solution. In both cases ${\mathcal{V}}(t)/{\mathcal{V}}(0)$ decays like $t^{-2}$, but the details are different. The similarity solutions admit an adjustable virtual-origin constant, which we determine by matching with the finite-difference solution. The analysis is valid for both Boussinesq and non-Boussinesq systems, and a wide range of geometric parameters (inner and outer radii, and height). The importance of the neglected viscous terms increases with time, and eventually the inertial–buoyancy model becomes invalid. An estimate for this occurrence is also provided. The predictions of the model are compared to results of direct numerical simulations; there is good agreement for the position of the interface and for the averaged radial velocity, and excellent agreement for ${\mathcal{V}}(t)/{\mathcal{V}}(0)$. A box model is used for estimating the effect of a partial (over a sector) dam break. This study is an extension of the work for a rectangular reservoir of Momen et al. (J. Fluid Mech., vol. 827, 2017, pp. 640–663). We demonstrate that there are some similarities, but also significant differences, between the rectangular and the cylindrical reservoirs concerning the velocity, shape of the interface and rate of drainage, which are of interest in applications. The overall conclusion is that this simple model captures very well the flow field under consideration.