No Arabic abstract
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.
Viscoelastic flow through an abrupt planar contraction geometry above a certain Weissenberg number (Wi) is well known to become unstable upstream of the contraction plane via a central jet separating from the walls and forming vortices in the salient corners. Here, for the first time we consider three-dimensional (3D) viscoelastic contraction flows in a microfabricated glass square-square contraction geometry. We employ state-of-the-art microtomographic particle image velocimetry to produce time-resolved and volumetric quantification of the 3D viscoelastic instabilities arising in a dilute polymer solution driven through the geometry over a wide range of Wi but at negligible Reynolds number. Based on our observations, we describe new insights into the growth, propagation, and transient dynamics of an elastic vortex formed upstream of the 3D micro-contraction due to flow jetting towards the contraction. At low Wi we observe vortex growth for increasing Wi, followed by a previously unreported vortex growth plateau region. In the plateau region, the vortex circulates around the jet with a period that decreases with Wi but an amplitude that is independent of Wi. In addition, we report new out-of-plane asymmetric jetting behaviour with a phase-wise dependence on Wi. Finally, we resolve the rate-of-strain tensor D and ascribe local gradients in D as the underlying driver of circulation via strain-hardening of the fluid in the wake of the jet.
It is argued herein that when PIV is used to measure turbulence, it can be treated as a time-dependent signal. The `output velocity consists of three primary contributions: the time-dependent velocity, a noise arising from the quantization (or pixelization), and a noise contribution from the fact that the velocity is not uniform inside the interrogation volume. For both of the latter their variances depend inversely on the average number of particles or images) in this interrogation volume. All three of these are spatially filtered by the finite extent of the interrogation window. Since the above noises are associated directly with the individual particles (or particle images), the noise between different realizations and different interrogation volumes is statistically independent.
In the past decade, advances in electronics technology have made larger imaging sensors available to the experimental fluid mechanics community. These advancements have enabled the measurement of 2-component 2-dimensional (2C-2D) velocity fields using particle image velocimetry (PIV) with much higher spatial resolution than previously possible. However, due to the large size of the sensor, the lens distortion needs to be taken into account as it will now have a more significant effect on the measurement quality that must be corrected to ensure accurate high-fidelity 2C-2D velocity field measurements. In this paper, two dewarping models, a second-order rational function (R2) and a bicubic polynomial (P3) are investigated with regards to 2C-2D PIV measurements of a turbulent boundary layer (TBL) using a large imaging sensor. Two approaches are considered and compared: (i) dewarping the images prior to the PIV cross-correlation analysis and (ii) undertaking the PIV cross-correlation analysis using the original recorded distorted images then followed by using the mapping functions derived for image dewarping to provide the correct spatial location of the velocity measurement point. The results demonstrate that the use of P3 dewarping model to correct lens distortion yields better results than the R2 dewarping model. Furthermore, both approaches for the P3 dewarping model yield results which are statistically indistinguishable.
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.
The physical characteristics and evolution of a large-scale helium plume are examined through a series of numerical simulations with increasing physical resolution using adaptive mesh refinement (AMR). The five simulations each model a 1~m diameter circular helium plume exiting into a (4~m)$^3$ domain, and differ solely with respect to the smallest scales resolved using the AMR, spanning resolutions from 15.6~mm down to 0.976~mm. As the physical resolution becomes finer, the helium-air shear layer and subsequent Kelvin-Helmholtz instability are better resolved, leading to a shift in the observed plume structure and dynamics. In particular, a critical resolution is found between 3.91~mm and 1.95~mm, below which the mean statistics and frequency content of the plume are altered by the development of a Rayleigh-Taylor instability near the centerline in close proximity to the base of the plume. This shift corresponds to a plume puffing frequency that is slightly higher than would be predicted using empirical relationships developed for buoyant jets. Ultimately, the high-fidelity simulations performed here are intended as a new validation dataset for the development of subgrid-scale models used in large eddy simulations of real-world buoyancy-driven flows.