Direct Numerical Simulations of Planar and Cylindrical Density Currents

The collapse of a heavy fluid column in a lighter environment is studied by direct numerical simulation of the Navier-Stokes equations using the Boussinesq approximation for small density difference. Such phenomenon occurs in many engineering and environmental problems resulting in a density current spreading over a no-slip boundary. In this work, density currents corresponding to two Grashof (Gr) numbers are investigated (⁠$105$ and $1.5×106$⁠) for two very different geometrical configurations, namely, planar and cylindrical, with the goal of identifying differences and similarities in the flow structure and dynamics. The numerical model is capable of reproducing most of the two- and three-dimensional flow structures previously observed in the laboratory and in the field. Soon after the release of the heavier fluid into the quiescent environment, a density current forms exhibiting a well-defined head with a hanging nose followed by a shallower body and tail. In the case of large Gr, the flow evolves in a three-dimensional fashion featuring a pattern of lobes and clefts in the intruding front and substantial three-dimensionality in the trailing body. For the case of the lower Gr, the flow is completely two dimensional. The dynamics of the current is visualized and explained in terms of the mean flow for different phases of spreading. The initial phase, known as slumping phase, is characterized by a nearly constant spreading velocity and strong vortex shedding from the front of the current. Our numerical results show that this spreading velocity is influenced by Gr as well as the geometrical configuration. The slumping phase is followed by a decelerating phase in which the vortices move into the body of the current, pair, stretch and decay as viscous effects become important. The simulated dynamics of the flow during this phase is in very good agreement with previously reported experiments.