No Arabic abstract
Gravity currents modify their flow characteristics by entraining ambient fluid, which depends on a variety of governing parameters such as the initial density, $Delta rho$, the total initial height of the fluid, $H$, and the slope of the terrain, $alpha$, from where it is released. Depending on these parameters, the gravity current may be designated as sub-critical, critical, or super-critical. It is imperative to study the entrainment dynamics of a gravity current in order to have a clear understanding of mixing transitions that govern the flow physics, the shear layer thickness, $delta_{u}$, and the mixing layer thickness, $delta_{rho}$. Experiments were conducted in a lock-exchange facility in which the dense fluid was separated from the ambient lighter fluid using a gate. As the gate is released instantaneously, an energy conserving gravity current is formed, for which the only governing parameter is the Reynolds number defined as $Re=frac{Uh}{ u}$, where $U$ is the front velocity of the gravity current, and $h$ is the height of the current. In our study, the bulk Richardson number, $Ri_{b}$=$frac{g^{}H}{U_{b}^{2}}$=1, takes a constant value for all the experiments, with $U_{b}$ being the bulk velocity of the layer defined as $U_{b}$=$sqrt{g^{}H}$. Simultaneous Particle Image Velocimetry (PIV) and Planar Laser Induced Fluorescence (PLIF) measurement techniques are employed to get the velocity and density statistics. A flux-based method is used to calculate the entrainment coefficient, E$_{F}$, for a Reynolds number range of $Reapprox$400-13000 used in our experiments. The result shows a mixing transition at $Reapprox$2700 that is attributed to the flow transitioning from weak Holmboe waves to Kelvin-Helmholtz type instabilities.
A series of laboratory experiments on energy conserving gravity currents in a lock-exchange facility are conducted for a range of Reynolds numbers, $Re= frac{U_Fh}{ u} =$ 485-12270. The velocity and density fields are captured simultaneously using a PIV-PLIF system. A moving average method is employed to compute the mean field and a host of turbulence statistics, namely, turbulent kinetic energy ($K$), shear production ($P$), buoyancy flux ($B$), and energy dissipation ($epsilon$) during the slumping phase of the current. The subsequent findings are used to ascertain the quantitative values of mixing efficiency, $Ri_{f}$, Ozmidov length-scale ($L_O$), Kolmogorov length-scale ($L_kappa$), and eddy diffusivities of momentum ($kappa_m$) and scalar ($kappa_rho$). Two different forms of $Ri_{f}$ are characterized in this study, denoted by $Ri_{f}^I=frac{B}{P}$ and $Ri_{f}^{II}=frac{B}{B+epsilon}$. The results cover the entire diffusive regime (3 $<Re_b<$ 10) and a portion of the intermediate regime (10 $<Re_b<$ 50), where $Re_b=frac{epsilon}{ u N^2}$ is the buoyancy Reynolds number that measures the level of turbulence in a shear-stratified flow. The values of $P$, $B$, and $epsilon$ show a marked increase at the interface of the ambient fluid and the current, owing to the development of a shear-driven mixed layer. Based on the changes in the turbulence statistics and the length scales, it is inferred that the turbulence decays along the length of the current. The mixing efficiency monotonically increases in the diffusive regime ($Re_{b}<$10), and is found to have an upper bound of $Ri_{f}^{I}approx$ 0.15 and $Ri_{f}^{II}approx$ 0.2 in the intermediate regime. Using the values of $Ri_{f}$, the normalized eddy diffusivity of momentum is parameterized as $frac{kappa_m}{ u.Ri_{g}}$=1.2$Re_{b}$ and normalized eddy diffusivity of scalar as $frac{kappa_{rho}}{ u}$=0.2$Re_{b}$
We present experimental measurements of a wall-bounded gravity current, motivated by characterizing natural gravity currents such as oceanic overflows. We use particle image velocimetry and planar laser-induced fluorescence to simultaneously measure the velocity and density fields as they evolve downstream of the initial injection from a turbulent channel flow onto a plane inclined at 10$^circ$ with respect to horizontal. The turbulence level of the input flow is controlled by injecting velocity fluctuations upstream of the output nozzle. The initial Reynolds number based on Taylor microscale of the flow, R$_lambda$, is varied between 40 and 120, and the effects of the initial turbulence level are assessed. The bulk Richardson number $Ri$ for the flow is about 0.3 whereas the gradient Richardson number $Ri_g$ varies between 0.04 and 0.25, indicating that shear dominates the stabilizing effect of stratification. Kelvin-Helmholtz instability results in vigorous vertical transport of mass and momentum. We present baseline characterization of standard turbulence quantities and calculate, in several different ways, the fluid entrainment coefficient $E$, a quantity of considerable interest in mixing parameterization for ocean circulation models. We also determine properties of mixing as represented by the flux Richardson number $Ri_f$ as a function of $Ri_g$ and diapycnal mixing parameter $K_rho$ versus buoyancy Reynolds number $Re_b$. We find reasonable agreement with results from natural flows.
We study entrainment in dry thermals in neutrally and unstably stratified ambients, and moist thermals in dry-neutrally stratified ambients using direct numerical simulations (DNS). We find, in agreement with results of Lecoanet and Jeevanjee [1] that turbulence plays a minor role in entrainment in dry thermals in a neutral ambient for Reynolds numbers $Re < 10^4$ . We then show that the net entrainment rate increases when the buoyancy of the thermals increases, either by condensation heating or because of an unstably stratified ambient. This is in contrast with the findings of Morrison et al. [2]. We also show that the role of turbulence is greater in these cases than in dry thermals and, significantly, that the combined action of condensation heating and turbulence creates intense small scale vorticity, destroying the vortex ring that is seen in dry and moist laminar thermals. These findings suggest that fully resolved simulations at Reynolds numbers significantly larger than the mixing transition Reynolds number $Re = 10^4$ are necessary to understand the role of turbulence in the entrainment in growing cumulus clouds, which consist of a series of thermals rising and decaying in succession.
In this paper, we employ Lagrangian coherent structures (LCSs) theory for the three dimensional vortex eduction and investigate the effect of large-scale vortical structures on the turbulent/non-turbulent interface (TNTI) and entrainment of a gravity current. The gravity current is realized experimentally and different levels of stratification are examined. For flow measurements, we use a multivolume three-dimensional particle tracking velocimetry technique. To identify vortical LCSs (VLCSs), a fully automated 3D extraction algorithm for multiple flow structures based on the so-called Lagrangian-Averaged Vorticity Deviation method is implemented. The size, the orientation and the shape of the VLCSs are analyzed and the results show that these characteristics depend only weakly on the strength of the stratification. Through conditional analysis, we provide evidence that VLCSs modulate the average TNTI height, affecting consequently the entrainment process. Furthermore, VLCSs influence the local entrainment velocity and organize the flow field on both the turbulent and non-turbulent sides of the gravity current boundary.
The presence of stratified layer in atmosphere and ocean leads to buoyant vertical motions, commonly referred to as plumes. It is important to study the mixing dynamics of a plume at a local scale in order to model their evolution and growth. Such a characterization requires measuring the velocity and density of the mixing fluids simultaneously. Here, we present the results of a buoyant plume propagating in a linearly stratified medium with a density difference of 0.5%, thus yielding a buoyancy frequency of N=0.15 s^{-1}. To understand the plume behaviour, statistics such as centerline and axial velocities along varying downstream locations, turbulent kinetic energy, Reynolds stress, and buoyancy flux were measured. The centerline velocity was found to decrease with increase in height. The Reynolds stress and buoyancy flux profiles showed the presence of a unstable layer and the mixing associated within that layer.