Search for: All records

We conduct three-dimensional direct numerical simulations to investigate the mixing, entrainment and energy budgets of gravity currents emerging from two-layer stratified locks. Depending on the density and layer thickness ratios, we find that either the upper layer or lower layer fluid can propagate faster, and that the density structure of the overall gravity current can range from strongly stratified to near-complete mixing. We furthermore observe that intermediate values of the density ratio can maximise mixing between the gravity current layers. Based on the vorticity budget, we propose a theoretical model for predicting the overall gravity current height, along with the front velocity of the two layers, for situations in which the lower layer moves faster than the upper layer. The model identifies the role of the height and thickness ratios in determining the velocity structure of the current, and it clarifies the dynamics of the ambient counter-current. A detailed analysis of the energy budget quantifies the conversion of potential into kinetic energy as a function of the governing parameters, along with the energy transfer between the different layers of the gravity current and the ambient fluid. Depending on the values of the density and layer thickness ratios, we find that the lower lock layer can gain or lose energy, whereas the upper layer always loses energy. 

We investigate the submerged collapse of weakly polydisperse, loosely packed cohesive granular columns, as a function of aspect ratio and cohesive force strength, via grain-resolving direct numerical simulations. The cohesive forces act to prevent the detachment of individual particles from the main body of the collapsing column, reduce its front velocity, and yield a shorter and thicker final deposit. All of these effects can be captured accurately across a broad range of parameters by piecewise power-law relationships. The cohesive forces reduce significantly the amount of available potential energy released by the particles. For shallow columns, the particle and fluid kinetic energy decreases for stronger cohesion. For tall columns, on the other hand, moderate cohesive forces increase the maximum particle kinetic energy, since they accelerate the initial free-fall of the upper column section. Only for larger cohesive forces does the peak kinetic energy of the particles decrease. Computational particle tracking indicates that the cohesive forces reduce the mixing of particles within the collapsing column, and it identifies the regions of origin of those particles that travel the farthest. The simulations demonstrate that cohesion promotes aggregation and the formation of aggregates. Furthermore, they provide complete information on the temporally and spatially evolving network of cohesive and direct contact force bonds. While the normal contact forces are aligned primarily in the vertical direction, the cohesive bonds adjust their preferred spatial orientation throughout the collapse. They result in a net macroscopic stress that counteracts deformation and slows the spreading of the advancing particle front.