2 results
Energy balances for the collision of gravity currents of equal strengths
- Albert Dai, Yu-Lin Huang, Ching-Sen Wu
-
- Journal:
- Journal of Fluid Mechanics / Volume 959 / 25 March 2023
- Published online by Cambridge University Press:
- 17 March 2023, A20
-
- Article
-
- You have access Access
- Open access
- HTML
- Export citation
-
Collision of two counterflowing gravity currents of equal densities and heights was investigated by means of three-dimensional high-resolution simulations with the goal of understanding the flow structures and energetics in the collision region in more detail. The lifetime of collision is approximately $3 \tilde {H}/\tilde {u}_f$, where $\tilde {H}$ is the depth of heavy and ambient fluids, and $\tilde {u}_f$ is the front velocity of the approaching gravity currents, and the lifetime of collision can be divided into three phases. During Phase I, $-0.2 \leqslant (\tilde {t}-\tilde {t}_c) \tilde {u}_f/\tilde {H} \leqslant 0.5$, where $\tilde {t}$ is the time, and $\tilde {t}_c$ is the time instance at which the two colliding gravity currents have fully osculated, geometric distortions of the gravity current fronts result in stretching of pre-existing vorticity in the wall-normal direction inside the fronts, and an array of vertical vortices extending throughout the updraught fluid column develop along the interface separating the two colliding gravity currents. The array of vertical vortices is responsible for the mixing between the heavy fluids of the two colliding gravity currents and for the production of turbulent kinetic energy in the collision region. The presence of the top boundary deflects the updraughts into the horizontal direction, and a number of horizontal streamwise vortices are generated close to the top boundary. During Phase II, $0.5 \leqslant (\tilde {t}-\tilde {t}_c) \tilde {u}_f/\tilde {H} \leqslant 1.2$, the horizontal streamwise vortices close to the top boundary induce turbulent buoyancy flux and break up into smaller structures. While the production of turbulent kinetic energy weakens, the rate of transfer of energy to turbulent flow due to turbulent buoyancy flux reaches its maximum and becomes the primary supply in the turbulent kinetic energy in Phase II. During Phase III, $1.2 \leqslant (\tilde {t}-\tilde {t}_c) \tilde {u}_f/\tilde {H} \leqslant 2.8$, the collided fluid slumps away from the collision region, while the production of turbulent kinetic energy, turbulent buoyancy flux and dissipation of energy attenuate. From the point of view of energetics, the production of turbulent kinetic energy and turbulent buoyancy flux transfers energy away from the mean flow to the turbulent flow during the collision. Our study complements previous experimental investigations on the collision of gravity currents in that the flow structures, spatial distribution and temporal evolution of the mean flow and turbulent flow characteristics in the collision region are presented clearly. It is our understanding that such complete information on the energy budgets in the collision region can be difficult to attain in laboratory experiments.
High-resolution simulations of cylindrical gravity currents in a rotating system
- Albert Dai, Ching-Sen Wu
-
- Journal:
- Journal of Fluid Mechanics / Volume 806 / 10 November 2016
- Published online by Cambridge University Press:
- 29 September 2016, pp. 71-101
-
- Article
- Export citation
-
Cylindrical gravity currents, produced by a full-depth lock release, in a rotating system are investigated by means of three-dimensional high-resolution simulations of the incompressible variable-density Navier–Stokes equations with the Coriolis term and using the Boussinesq approximation for a small density difference. Here, the depth of the fluid is chosen to be the same as the radius of the cylindrical lock and the ambient fluid is non-stratified. Our attention is focused on the situation when the ratio of Coriolis to inertia forces is not large, namely $0.1\leqslant {\mathcal{C}}\leqslant 0.3$, and the non-rotating case, namely ${\mathcal{C}}=0$, is also briefly considered. The simulations reproduce the major features observed in the laboratory and provide more detailed flow information. After the heavy fluid contained in a cylindrical lock is released in a rotating system, the influence of the Coriolis effects is not significant during the initial one-tenth of a revolution of the system. During the initial one-tenth of a revolution of the system, Kelvin–Helmholtz vortices form and the rotating cylindrical gravity currents maintain nearly perfect axisymmetry. Afterwards, three-dimensionality of the flow quickly develops and the outer rim of the spreading heavy fluid breaks away from the body of the current, which gives rise to the maximum dissipation rate in the system during the entire adjustment process. The detached outer rim of heavy fluid then continues to propagate outward until a maximum radius of propagation is attained. The body of the current exhibits a complex contraction–relaxation motion and new outwardly propagating pulses form regularly in a period slightly less than half-revolution of the system. Depending on the ratio of Coriolis to inertia forces, such a contraction–relaxation motion may be initiated after or before the attainment of a maximum radius of propagation. In the contraction–relaxation motion of the heavy fluid, energy is transformed between potential energy and kinetic energy, while it is mainly the kinetic energy that is consumed by the dissipation. As a new pulse initially propagates outward, the potential energy in the system increases at the expense of decreasing kinetic energy, until a local maximum of potential energy is reached. During the latter part of the new pulse propagation, the kinetic energy in the system increases at the expense of decreasing potential energy, until a local minimum of potential energy is reached and another new pulse takes form. With the use of three-dimensional high-resolution simulations, the lobe-and-cleft structure at the advancing front can be clearly observed. The number of lobes is maintained only for a limited period of time before merger between existing lobes occurs when a maximum radius of propagation is approached. The high-resolution simulations complement the existing shallow-water formulation, which accurately predicts many important features and provides insights for rotating cylindrical gravity currents with good physical assumptions and simple mathematical models.