Instantaneous Particle Image Velocimetry Research Articles

Investigation of fluid flow during flow boiling inside a horizontal rectangular channel with single-sided heating using particle image velocimetry

Subcooled flow boiling is a highly efficient cooling systems for thermal management systems. This study explores the intricate dynamics of subcooled flow boiling within a horizontal channel, investigating the impact of vapor generation on liquid-phase velocity using Particle Image Velocimetry (PIV) and advanced image processing techniques. Four mass flow rates ranging from 5–20 g/s with subcooled inlet conditions are investigated in a rectangular channel with single-sided heating. Three regions of interest along the heated channel are investigated for instantaneous PIV analysis. The PIV system captures detailed velocity profiles, illustrating the impact of varying mass flow rates and heat flux levels on flow behavior. Vapor masking techniques are introduced to enhance the precision of PIV data by mitigating interference from the vapor phase. Results demonstrate the influence of vapor bubbles on flow resistance, revealing non-uniform velocity distributions and turbulence near the liquid–vapor interface. The study emphasizes the critical role of inertia and buoyancy forces in shaping the velocity profiles. Moreover, the investigation sheds light on the effects of flow rates on the interfacial behaviors, hinting at a transition point between 10 and 15 g/s. In summary, this research contributes valuable insights into the nuanced dynamics of flow boiling, laying the foundation for future studies on turbulence, heat transfer, and phase-change phenomena in two-phase thermal management systems.

Supersonic Flow Separation in Front of a Forward-Facing Step

A detailed experimental study has been done to understand the Shock Wave-Boundary Layer Interaction (SWBLI) over a Forward-Facing Step (FFS) of height (h) in a Mach 2.5 flow using Particle Image Velocimetry (PIV). In addition to PIV, unsteady wall pressure measurements, high-speed schlieren and surface oil flow visualization have been studied. Flow separation using surface oil flow visualization was observed around 4.1 step heights upstream of the step face. Shock motions extracted from the schlieren image sequence shows a peak frequency of approximately 1000 Hz, which corresponds to the dominant (low) frequency at which the shock oscillates; this being two orders in magnitude smaller than the characteristics frequency of the incoming boundary layer. Also, spectra of the unsteady wall pressure data shows that the shock oscillates with a peak frequency of approximately 1000 Hz. The effect of shock location (xs) from instantaneous velocity PIV fields in the cross-stream plane due to an upstream influence parameter (line-averaged boundary layer streamwise velocity, ul) and a downstream parameter (separation bubble or recirculation region area, A) is studied. Instantaneous cross-stream PIV measurements show significant variations of the separation shock location, incoming boundary layer velocity fluctuations and recirculation region. Besides, conditionally averaged velocity field based on incoming boundary layer streamwise velocity fluctuations and separation bubble area in the cross-stream plane shows that the shock location (xs) is correlated to the size of the separation bubble and not to boundary layer velocity fluctuations.

Uncertainty estimation for ensemble particle image velocimetry

We present a novel approach to estimate the uncertainty in ensemble particle image velocimetry (PIV) measurements. The ensemble PIV technique is widely used when the cross-correlation signal-to-noise ratio is insufficient to perform a reliable instantaneous velocity measurement. Despite the utility of ensemble PIV, uncertainty quantification for this type of measurement has not been studied. Here, we propose a method for estimating the uncertainty directly from the probability density function of displacements found by deconvolving the ensemble cross-correlation from the ensemble autocorrelation. We then find the second moment of the probability density function and apply a scaling factor to report the uncertainty in the velocity measurement. We call this method the moment of probability of displacement (MPD). We assess MPD’s performance with synthetic and experimental images. We show that predicted uncertainties agree well with the expected root mean square (RMS) of the error in the velocity measurements over a wide range of image and flow conditions. MPD shows good sensitivity to various PIV error sources with around 86% accuracy in matching the RMS of the error in the baseline data sets. So, MPD establishes itself as a reliable uncertainty quantification algorithm for ensemble PIV. We compared the results of MPD against one of the existing instantaneous PIV uncertainty approaches, moment of correlation (MC). We adapted the MC approach for ensemble PIV, however, its primary limitations remain the assumption of the Gaussian probability density function of displacements and the Gaussian particles’ intensity profile. In addition, our analysis shows that ensemble MC consistently underestimates the uncertainty, while MPD outperforms that and removes the limiting Gaussian assumption for the particle and probability density function, thus overcoming the limitations of MC.

Uncertainty Estimation for Ensemble Particle Image Velocimetry

We present a novel approach to estimate the uncertainty in ensemble particle image velocimetry (PIV) measurements. Ensemble PIV is widely used when the cross-correlation signal-to-noise ratio (SNR) is insufficient to perform a reliable instantaneous velocity measurement. Despite the utility of ensemble PIV, uncertainty quantification for this type of measurement has not been studied. The existing uncertainty quantification algorithms for PIV are developed and used only for instantaneous PIV measurement and do not account for the improved SNR in ensemble PIV. Existing instantaneous uncertainty quantification methods can be divided into direct and indirect categories. Indirect methods require calibration based on the effect of various image parameters (such as noise, particle size, density, velocity gradient, etc.) on the correlation SNR. Indirect methods have not been calibrated for error sources relevant in an Ensemble PIV measurement. Also, they have lower sensitivity to the error sources compared to direct approaches. Direct methods, such as the moment of correlation (MC) and Image Matching (IM), find the uncertainty based on the images and correlation planes without any calibration and are more reliable (Bhattacharya et al., 2018; Sciacchitano et al., 2013). Ensemble PIV is based on ensemble correlations; therefore, MC, which uses the generalized cross-correlation (GCC) plane as a measure of uncertainty, is the most suitable method to be modified to be applicable for the ensemble PIV. The GCC plane is the inverse Fourier transform of the phase correlation and represents the probability density function (PDF) of particles’ displacements (Bhattacharya et al., 2018; Eckstein and Vlachos, 2009). We replaced instantaneous GCC with ensemble GCC and modified MC’s normalization factor to account for the number of ensembles. The MC’s primary limitation is that it assumes a Gaussian shape for the PDF of displacements and estimate the standard deviation of the underlying PDF using a fitted Gaussian. However, the PDF deviates from Gaussian distribution due to velocity gradient or non-Gaussian random displacements. Therefore, MC’s reliability and applicability are reduced for flow fields with non-Gaussian PDFs. Also, our analysis shows that ensemble MC consistently underestimates the uncertainty. So, a generalized and reliable method for uncertainty quantification for ensemble PIV is needed.

Experimental investigations on the sharp leading-edge separation over a flat plate at zero incidence using particle image velocimetry

Leading-edge separated flow field over a sharp flat plate is experimentally investigated in Reynolds numbers ranging from 6.2 × 103 to 4.1 × 104, using particle image velocimetry (PIV) and its statistics. It was observed that the average reattachment length is nearly independent of Reynolds number and the small secondary bubble observed near the leading edge was found to shrink with increasing Reynolds number. The wall-normal profiles of the statistical values of kinematic quantities such as the velocity components and their fluctuations scaled well with average reattachment length lR and freestream velocity U∞. Their magnitudes compare well with previous investigations even though the current triangular shaped sharp leading edge is different from previous flat-faced or semi-circular ones. The shear layer was observed to exhibit 2 different linear growth rates over 2 distinct regions. Instantaneous PIV realizations demonstrate unsteady nature of the separation bubble, whose origins in the upstream portion of the bubble are analysed. Bimodal nature of the probability density function (PDF) of fluctuating streamwise velocity at around x/lR = 0.08–0.15 indicates successive generation and passage of vortices in the region, which subsequently interact and evolve into multiscale turbulent field exhibiting nearly Gaussian PDF. Shedding of vortices with wide range of scales are apparent in most of the instantaneous realizations. Proper Orthogonal Decomposition (POD) of the velocity fluctuation magnitude field revealed that the flow structures of the dominant modes and their relative energies are independent of Reynolds number. In each of the dominant modes (first 3 modes), the length scales corresponding to the large scale structures and their spacing are the same for all Reynolds numbers, suggesting that their Strouhal number (observed to be ~ 0.09–0.2 at Reynolds number of 6.2 × 103) of unsteadiness should also be independent of Reynolds number. A single large structure- comparable in size to lR—was apparent well before reattachment in a few instantaneous realizations, as compared to multiple small-scale structures visible in most realizations; at Reynolds number of 6.2 × 103, realizations with such large-scale structures occurred approximately after every 20–30 realizations, corresponding to non-dimensional frequency of 0.4–0.6, which is identified to be the “regular shedding”. It was possible to reconstruct the large-scale structure during the instances from just the first 3 POD modes, indicating that the Strouhal number of regular shedding too is independent of Reynolds number.Graphic abstract

Experimental investigations of the asymmetric vortex formation over a slender body of revolution at high angles of attack

This paper describes an experimental investigation of the initial growth of flow asymmetries over a slender body of revolution at high angles of attack with natural and disturbed noses. Time-resolved particle image velocimetry (PIV) is used to investigate the flow field around the body. Flow visualization clearly shows the formation of the asymmetric vortices. Instantaneous PIV shows that the amplified asymmetric disturbances lead to Kelvin–Helmholtz instability appearing first on one side, which increases the momentum exchange crossing the layer. As a result, the separation region shrinks which creates the initial vortex asymmetry.

Large-scale motions in a plane wall jet

The plane wall jet (PWJ) is a wall-bounded flow in which a wall shear layer develops in the presence of extremely energetic flow structures of the outer free-shear layer. The structure of a PWJ, developing in still air, was studied with the focus on the large scales in the flow. Wall-normal hot-wire anemometry (HWA) measurements along with double-frame particle image velocimetry (PIV) measurements (wall-normal–streamwise plane) were carried out at streamwise distances up to $162b$, where $b$ is the slot width of the PWJ exit. The nominal PWJ Reynolds number based on exit parameters was $Re_{j}\approx 5940$. Comparisons with a zero-pressure-gradient boundary layer (ZPGBL) at nominally matched friction Reynolds number $Re_{\unicode[STIX]{x1D70F}}$ were also carried out as appropriate, to highlight key features of the PWJ structure. Consistent with previous work, the PWJ showed a dependence of the peak turbulent stresses on the jet exit Reynolds number. The turbulent production showed a peak corresponding to the near-wall cycle similar to the peak seen in the ZPGBL. However, another turbulent production peak was observed in the outer free-shear layer that was an order of magnitude larger than the inner one. Along with the change in sign of the viscous and Reynolds shear stresses, the PWJ was shown to have a region of very low turbulent production between these two peaks. The dissipation rate increased over the PWJ layer with a peak also in the outer region. Visualizations of the flow and two-point correlations reveal that the most energetic large-scale structures within a PWJ are vortical motions in the wall-normal–streamwise plane similar to those structures seen in free-shear layers. These structures are referred to as J (for jet) type structures. In addition two-point correlations reveal the existence of large-scale structures in the wall region which have a signature similar to those structures seen in canonical boundary layers. These structures are referred to as W (for wall) type structures. Instantaneous PIV realizations and flow visualizations reveal that these W type large-scale features are consistent with the paradigm of hairpin vortex packets in the wall region. The J type structures were seen to intrude well into the wall region while the W type structures were also seen to extend into the outer shear layer. Further, these large-scale structures were shown to modulate the amplitude of the finer scales of the flow.

Insights into the dynamics of conical breakdown modes in coaxial swirling flow field

The main idea of this paper is to understand the fundamental vortex breakdown mechanisms in the coaxial swirling flow field. In particular, the interaction dynamics of the flow field is meticulously addressed with the help of high fidelity laser diagnostic tools. Time-resolved particle image velocimetry (PIV) (${\sim}1500~\text{frames}~\text{s}^{-1}$) is employed in $y{-}r$ and multiple $r{-}\unicode[STIX]{x1D703}$ planes to precisely delineate the flow dynamics. Experiments are carried out for three sets of co-annular flow Reynolds number $Re_{a}=4896$, 10 545, 17 546. Furthermore, for each $Re_{a}$ condition, the swirl number ‘$S_{G}$’ is varied independently from $0\leqslant S_{G}\leqslant 3$. The global evolution of flow field across various swirl numbers is presented using the time-averaged PIV data. Three distinct forms of vortex breakdown namely, pre-vortex breakdown (PVB), central toroidal recirculation zone (CTRZ; axisymmetric toroidal bubble type breakdown) and sudden conical breakdown are witnessed. Among these, the conical form of vortex breakdown is less explored in the literature. In this paper, much attention is therefore focused on exploring the governing mechanism of conical breakdown. It is should be interesting to note that, unlike other vortex breakdown modes, conical breakdown persists only for a very short band of $S_{G}$. For any small increase/decrease in $S_{G}$ beyond a certain threshold, the flow spontaneously reverts back to the CTRZ state. Energy ranked and frequency-resolved/ranked robust structure identification methods – proper orthogonal decomposition (POD) and dynamic mode decomposition (DMD) respectively – are implemented over instantaneous time-resolved PIV data sets to extract the dynamics of the coherent structures associated with each vortex breakdown mode. The dominant structures obtained from POD analysis suggest the dominance of the Kelvin–Helmholtz (KH) instability (axial $+$ azimuthal; accounts for ${\sim}80\,\%$ of total turbulent kinetic energy, TKE) for both PVB and CTRZ while the remaining energy is contributed by shedding modes. On the other hand, shedding modes contribute the majority of the TKE in conical breakdown. The frequency signatures quantified from POD temporal modes and DMD analysis reveal the occurrence of multiple dominant frequencies in the range of ${\sim}10{-}400~\text{Hz}$ with conical breakdown. This phenomenon may be a manifestation of high energy contribution by shedding eddies in the shear layer. Contrarily, with PVB and CTRZ, the dominant frequencies are observed in the range of ${\sim}20{-}40~\text{Hz}$ only. We have provided a detailed exposition of the mechanism through which conical breakdown occurs. In addition, the current work explores the hysteresis (path dependence) phenomena of conical breakdown as functions of the Reynolds and Rossby numbers. It has been observed that the conical mode is not reversible and highly dependent on the initial conditions.

Real-time processing methods to characterize streamwise vortices

One of the bottlenecks of an active control of turbulent boundary layers based on particle image velocimetry (PIV) measurements relates to the ability to process PIV images and extract useful for control information in real-time, on a time scale relevant for the selected control strategy. We propose two methodologies to extract useful information from PIV measurements in the y−z spanwise - wall normal cross-section, based on characterisation of properties such as size and strength of streamwise vortices. The vortices are created by active vortex generators embedded in a fully turbulent boundary layer. The proposed methodologies combine vortex identification and characterisation methods that use instantaneous PIV realisations to extract centers and strengths of streamwise vortices. The main purpose is to compare a standard vortex identification method that provides a robust and accurate estimate with another, ad-hoc method that is less robust or accurate but has a large computational speed gain potential. For demonstration purposes we use PIV measurements obtained at the wind tunnel facility downstream active vortex generators (Carlier and Stanislas, 2005; Foucaut et al., 2014). The first algorithm uses the Q-criteria and the integration of vorticity of each extracted vortex. This robust and accurate method requires the full spatially resolved PIV field followed by the computationally expensive spatial derivatives calculations and an integration. The second algorithm utilises the prior knowledge about the presence and shape of streamwise vortices expected in the measurement, therefore uses only a single line velocity profile estimate and can be parallelized to obtain few horizontal velocity profiles at different wall-normal distances. We compare the two methods and discuss their potential in term of computation times and robustness. Results show that in a turbulent boundary layer, a moving window average of a number of instantaneous fields is needed for both methods.

Symmetry breaking and instabilities of conical vortex pairs over slender delta wings

An investigation of symmetry breaking and naturally occurring instabilities over thin slender delta wings with sharp leading edges was carried out in a water tunnel using particle image velocimetry (PIV) measurements. Time-averaged location, strength and core radius of conical vortices vary almost linearly with chordwise distance for three delta wings with $75^{\circ }$, $80^{\circ }$ and $85^{\circ }$ sweep angles over a wide range of angles of attack. Properties of the time-averaged vortex pairs depend only on the similarity parameter, which is a function of the angle of attack and the sweep angle. It is shown that time-averaged vortex pairs develop asymmetry gradually with increasing values of the similarity parameter. Vortex asymmetry can develop in the absence of vortex breakdown on the wing. Instantaneous PIV snapshots were analysed using proper orthogonal decomposition and dynamic mode decomposition, revealing the shear layer and vortex instabilities. The shear layer mode is the most periodic and more dominant for lower values of the similarity parameter. The Strouhal number based on the free stream velocity component in the cross-flow plane is a function of only the similarity parameter. The dominant frequency of the shear layer mode decreases with the increasing similarity parameter. The vortex modes reveal the fluctuations of the vorticity magnitude and helical displacement of the cores, but with little periodicity. There is little correlation between the fluctuations in the cores of the vortices.

Vortex Rope Formation in a High Head Model Francis Turbine

Francis turbine working at off-design operating condition experiences high swirling flow at the runner outlet. In the present study, a high head model Francis turbine was experimentally investigated during load rejection, i.e., best efficiency point (BEP) to part load (PL), to detect the physical mechanism that lies in the formation of vortex rope. For that, a complete measurement system of dynamic pressure, head, flow, guide vanes (GVs) angular position, and runner shaft torque was setup with corresponding sensors at selected locations of the turbine. The measurements were synchronized with the two-dimensional (2D) particle image velocimetry (PIV) measurements of the draft tube. The study comprised an efficiency measurement and maximum hydraulic efficiency of 92.4 ± 0.15% was observed at BEP condition of turbine. The severe pressure fluctuations corresponding to rotor–stator interaction (RSI), standing waves, and rotating vortex rope (RVR) have been observed in the draft tube and vaneless space of the turbine. Moreover, RVR in the draft tube has been decomposed into two different modes; rotating and plunging modes. The time of occurrence of both modes was investigated in pressure and velocity data and results showed that the plunging mode appears 0.8 s before the rotating mode. In the vaneless space, the plunging mode was captured before it appears in the draft tube. The physical mechanism behind the vortex rope formation was analyzed from the instantaneous PIV velocity vector field. The development of stagnation region at the draft tube center and high axial velocity gradients along the draft tube centerline could possibly cause the formation of vortex rope.

Vortex dynamics for flow over a circular cylinder in proximity to a wall

The dynamics of vortical structures in flow over a circular cylinder in the vicinity of a flat plate is investigated using particle image velocimetry (PIV). The cylinder is placed above the flat plate with its axis parallel to the wall and normal to the flow direction. The Reynolds number $Re_{D}$ based on the cylinder diameter $D$ is 1072 and the gap $G$ between the cylinder and the flat plate is varied from gap-to-diameter ratio $G/D=0$ to $G/D=3.0$. The flow statistics and vortex dynamics are strongly dependent on the gap ratio $G/D$. Statistics show that as the cylinder comes close to the wall ($G/D\leqslant 2.0$), the cylinder wake becomes more and more asymmetric and a boundary layer separation is induced on the flat plate downstream of the cylinder. The wake vortex shedding frequency increases with decreasing $G/D$ until a critical gap ratio (about $G/D=0.25$) below which the vortex shedding is irregular. The deflection of the gap flow away from the wall and its following interaction with the upper shear layer may be the cause of the higher shedding frequency. The vortex dynamics is investigated based on the phase-averaged flow field and virtual dye visualization in the instantaneous PIV velocity field. It is revealed that when the cylinder is close to the wall ($G/D=2.0$), the cylinder wake vortices can periodically induce secondary spanwise vortices near the wall. As the cylinder approaches the wall ($G/D=1.0$) the secondary vortex can directly interact with the lower wake vortex, and a further approaching of the cylinder ($G/D=0.5$) can result in more complex interactions among the secondary vortex, the lower wake vortex and the upper wake vortex. The breakdown of vortices into filamentary debris during vortex interactions is clearly revealed by the coloured virtual dye visualizations. For $G/D<0.25$, the lower shear layer is strongly inhibited and only the upper shear layer can shed vortices. Investigation of the vortex formation, evolution and interaction in the flow promotes the understanding of the flow physics for different gap ratios.

Insights into the dynamics of spray–swirl interactions

The near-field breakup and interaction of a hollow-cone liquid sheet with coannular swirling air flow have been examined using high-speed diagnostics. Time-resolved PIV (particle image velocimetry;$3500~\text{frames}~\text{s}^{-1}$) is employed to capture the spatio-temporal behaviour of the swirling air flow field. The combined liquid–gas phase interaction is visualized with the help of high-speed ($20\,000~\text{frames}~\text{s}^{-1}$) shadowgraphy. In this study, the transition from weak to strong spray–swirl interaction is explained based on the momentum ratio. Proper orthogonal decomposition (POD) is implemented on instantaneous PIV and shadowgraphy images to extract the energetic spatial eigenmodes and characteristic modal frequencies. The POD results suggest the dominance of the KH (Kelvin–Helmholtz) instability mechanism (pure axial shear, axial plus azimuthal shear) in swirl–spray interaction. In addition, linear stability analysis also shows the destabilization of the liquid–air interface caused by KH waves ($\unicode[STIX]{x1D706}_{p}$), which arises from the formation of a vorticity layer of thickness$\unicode[STIX]{x1D6FF}_{g}$near the liquid–air interface. The frequency values obtained from the primary KH wavelength ($\unicode[STIX]{x1D706}_{p}$) exhibit good agreement with the POD modal frequencies. Scaling laws are proposed to elucidate the relationships between the global length scales (breakup length, spray spread) and the primary wavelength. The breakup length scale and liquid sheet oscillations are meticulously analysed in the time domain to reveal the breakup dynamics of the liquid sheet. Furthermore, the large-scale coherent structures of the swirl flow exhibit different sheet breakup phenomena in the spatial domain. For instance, flapping breakup is induced by the central toroidal recalculation zone in the swirling flow field. Finally, the ligament formation mechanism and its diameter, i.e. the size of first-generation droplets, are measured with phase Doppler interferometry. The measured sizes scale reasonably with KH waves.

Low-frequency behavior of the turbulent axisymmetric near-wake

The turbulent wake past an axisymmetric body is investigated with time-resolved stereoscopic particle image velocimetry (PIV) at a Reynolds number ReD = 6.7 × 104 based on the object diameter. The azimuthal organization of the near-wake is studied at different locations downstream of the trailing edge. The time-averaged velocity field features a circular shear layer bounding a region of recirculating flow. Inspection of instantaneous PIV snapshots reveals azimuthal meandering of the reverse flow region with a significant radial offset with respect to the time-averaged position. The backflow meandering appears as the major contribution to the near-wake dynamics in proximity of the base, whereas closer to the rear-stagnation point, the shear layer fluctuations become important. For x/D ≤ 0.75, the time-history and probability distributions of the backflow centroid position allow to identify this motion with an irregular precession about the model symmetry axis occurring at time scales in the order of 103 D/U∞ or higher. The first two modes obtained by snapshot proper orthogonal decomposition of the velocity fluctuations can be related to an anti-symmetric mode of azimuthal wave-number m = 1 reflecting a radial displacement of the separated flow region, while the third and fourth proper orthogonal decomposition modes are identified with a second mode pair m = 2 and are representing wake ovalization. Close to the base, a third axisymmetric mode m = 0 is identified, corresponding to a streamwise pulsation of the reverse flow region. Based on the analysis of the spatial eigen-functions and frequency spectra of the time-coefficients, it is concluded that the anti-symmetric mode m = 1 is associated with the backflow instability in the very-low frequency range StD = 10−4/10−3 close to separation, whereas more downstream it reflects the fluctuations related to the shear layer development.

Coherent Structures Formation During Wake-Boundary Layer Interaction on a LP Turbine Blade

Particle Image Velocimetry (PIV) measurements have been analyzed in order to characterize the dynamics of coherent structures (eddies and streaks) within the suction side boundary layer of a low pressure turbine cascade perturbed by impinging wakes. To this end, the instantaneous flow fields at low Reynolds number and elevated free-stream turbulence intensity level (simulating the real condition of the blade row within the engine) were investigated in two orthogonal planes (a blade-to-blade and a wall-parallel plane). Proper Orthogonal Decomposition (POD) has been employed to filter the instantaneous flow maps allowing a better visualization of the structures involved in the transition process of the boundary layer. For the unsteady case properly selected POD modes have been also used to sort the instantaneous PIV images in the wake passage period. This procedure allows computing phase-averaged data and visualizing structures size and intensity in the different parts of the boundary layer during the different wake passage phases. The contributions to the whole shear stress due to the largest spanwise oriented scales at the leading and trailing boundaries of the wake-jet structures and those associated with streaky structures observed in the bulk of the wake are discussed. Instantaneous images in the wall-parallel plane are filtered with POD and they allow us to further highlight the occurrence of low and high speed traveling streaks (Klebanoff mode). The periodic advection along the suction side of the high turbulent content regions carried by the wakes anticipates both formation and sinuous instability of the streaks inside the boundary layer as compared with the steady case. The dynamics driving the breakdown of the streaks and the consequent formation of nuclei with high wall-normal vorticity have been found to be almost the same in the steady and the unsteady cases. Auto-correlation of the instantaneous images are also presented in order to highlight analogies and differences in the size and spacing of streaks in the two cases. These results are also compared with the available literature concerning simplified geometries (i.e flat plate) operating under steady inflow.

Stability analysis in solid–liquid fluidized beds: Experimental and computational

Abstract In this study the transition from homogeneous to heterogeneous flow in a solid–liquid fluidized bed (SLFB) is examined both experimentally and numerically. The experimental apparatus comprised a refractive index-matched SLFB, comprising 5 mm diameter borosilicate glass and sodium iodine solution, which allowed for both instantaneous particle image velocimetry of the liquid flow field and solids hold-up measurements to be undertaken for superficial liquid velocities in the range of 0.06–0.22 m/s. The motion of individual, spherical steel balls (with diameters 6, 7, 8, 9 mm) was then tracked as it settled through the fluidized bed for differing superficial liquid velocities. It was observed that, for all the steel balls covered in this work, there was a change in slope in their respective classification velocity curves at a superficial liquid velocity of 0.08 m/s. This value was very close to the critical velocity of 0.085 m/s predicted from 1-D linear stability analysis; and therefore deemed to be the critical condition that marked the transition from homogeneous to non-homogenous flow. It is proposed that the change in slope of the classification velocity curve is due to the encounter of the settling foreign particles with liquid bubbles whose presence marks the onset of heterogeneous flow. Additional computational analysis, involving both Eulerian–Eulerian (E–E) and Eulerian–Lagrangian (E–L) approaches, is used to confirm the presence of liquid bubbles at a critical liquid hold-up of 0.54, which corresponds to that predicted from 1-D linear stability analysis. In summary, the study has highlighted that experimentally the transition condition for a SLFB can be obtained simply by observing the behavior of the classification velocity of a single foreign particle at different superficial liquid velocities. This transition condition was found to agree with the 1D linear stability criterion, Eulerian–Eulerian CFD (3D) and Eulerian–Lagrangian DEM (3D) approaches.

Multispectral processing for color particle image velocimetry

Since the adoption of digital video cameras and cross-correlation methods for particle image velocimetry (PIV), the use of color images has largely been abandoned. Recently, however, with the re-emergence of color-based stereo and volumetric techniques, and the extensive use of color microscopy, color imaging for PIV has again become relevant. In this work, we explore the potential advantages of color PIV processing by developing and proposing new methods for handling multi-color images. The first method uses cross-correlation of every color channel independently to build a color vector cross-correlation plane. The vector cross-correlation can then be searched for one or more peaks corresponding to either the average displacement of several flow components using a color ensemble operation, or for the individual motion of colored particles, each with a different behavior. In the latter case, linear unmixing is used on the correlation plane to separate each known particle type as captured by the different color channels. The second method introduces the use of quaternions to encode the color data, and the cross-correlation is carried out simultaneously on all colors. The resulting correlation plane can be searched either for a single peak, corresponding to the mean flow or for multiple peaks, with velocity phase separation to determine which velocity corresponds to which particle type. Each of these methods was tested using synthetic images simulating the color recording of noisy particle fields both with and without the use of a Bayer filter and demosaicing operation. It was determined that for single-phase flow, both color methods decreased random errors by approximately a factor of two due to the noise signal being uncorrelated between color channels, while maintaining similar bias errors as compared to traditional monochrome PIV processing. In multi-component flows, the color vector correlation technique was able to successfully resolve displacements of two distinct yet coupled flow components with errors similar to traditional grayscale PIV processing of a single phase. It should be noted that traditional PIV processing is bound to fail entirely under such processing conditions. In contrast, the quaternion methods frequently failed to properly identify the correct velocity and phase and showed significant cross talk in the measurements between particle types. Finally, the color vector method was applied to experimental color images of a microchannel designed for contactless dielectrophoresis particle separation, and good results were obtained for both instantaneous and ensemble PIV processing. However, in both the synthetic color images that were generated using a Bayer filter and the experimental data, a significant peak-locking effect with a period of two pixels was observed. This effect is attributed to the inherent architecture of the Bayer filter. In order to mitigate this detrimental artifact, it is suggested that improved image interpolation or demosaicing algorithms tuned for use in PIV be developed and applied on the color images before processing, or that cameras that do not use a Bayer filter and therefore do not require a demosaicing algorithm be used for color PIV.

Modeling of microburst outflows using impinging jet and cooling source approaches and their comparison

Microbursts have been simulated and studied using different physical and numerical modeling methods. In the present study, the steady impinging jet model was comprehensively studied by using a 2-foot-diameter (0.61m) microburst simulator available in the Department of Aerospace Engineering at Iowa State University. Point measurements and Particle Image Velocimetry (PIV) results revealed a detailed picture of the overall flow and distribution of velocity and turbulence in the outflow of the steady impinging jet. Comparisons suggested that the average wind velocity profile of the steady impinging jet matched well with those derived from field data and previous research. FFT of the velocity time-history and instantaneous PIV results implied that the outflow consisted of low-frequency periodic shedding of vortices and the steady impinging jet model could be seen as an ensemble average of a series of simulated microburst events. Due to lack of time-dependent evolutionary information of the steady impinging jet model, a transient impinging jet model was studied to capture the transient features which were then compared with those of the cooling-source model by performing numerical simulations. Transient features of the transient impinging jet model and cooling source model showed several differences mainly related to the different formation and transportation process of the primary vortex. Ground surface pressure distributions were found to be different due to different forcing parameter of the two models. Comparison with the field data suggested that both models resembled the dynamic features of a real microburst outflow. However, results showed that the cooling source model could produce a reasonable instantaneous radial velocity profile at maximum wind condition, while the transient impinging jet model resulted in a large deviation. Finally, merits and demerits of each modeling methods were discussed.

Turbulent near wake of an Ahmed vehicle model

The lasting high fuel cost has recently inspired resurgence in drag reduction research for vehicles, which calls for a thorough understanding of the vehicle wake. The simplified Ahmed vehicle model is characterized by controllable flow separation, thus especially suitable for this purpose. In spite of a considerable number of previous investigations, our knowledge of flow around this model remains incomplete. This work aims to revisit turbulent flow structure behind this model. Two rear slant angles, i.e., α = 25o and 35o, of the model were examined, representing two distinct flow regimes. The Reynolds number was 5.26 × 104 based on the model height (H) and incident flow velocity. Using particle image velocimetry (PIV), flow was measured with and without a gap (g/H = 0.174) between the vehicle underside and ground in three orthogonal planes, viz. the x–z, x–y and y–z planes, where x, y, and z are the coordinates along longitudinal, transverse, and spanwise directions, respectively. The flow at g/H = 0 serves as an important reference for the understanding of the highly complicated vehicle wake (g/H ≠ 0). While reconfirming the well-documented major characteristics of the mean flow structure, both instantaneous and time-averaged PIV data unveil a number of important features of the flow structure, which have not been previously reported. As such, considerably modified flow structure models are proposed for both regimes. The time-averaged velocities, second moments of fluctuating velocities, and vorticity components are presented and discussed, along with their dependence on g/H in the two distinct flow regimes.

Transient and time-averaged characteristics of a compressible ground vortex flow

Particle image velocimetry data are presented from a scaled jet-lift aircraft ground vortex compressible flow. The scaled ground vortex is generated by a vertical compressible jet in cross-flow impinging on a moving ground plane. Particle image velocimetry is used to generate both transient and time-averaged flow statistics from the ground vortex region over a range of nozzle pressure ratios from 2.3 to 3.7, nozzle height-to-diameter ratios ( h/ dn) from 3 to 10 (where dn = 12.7 mm) and cross-flow velocities ( V∞) from 10 to 20 m/s. These conditions correspond to effective (jet-to-cross flow) velocity ratios of 15 < [Formula: see text] < 60. For each condition, mean and root mean square ground vortex core position was analysed from sets of 72 instantaneous particle image velocimetry vector maps. Over the range of effective velocity ratios, [Formula: see text], the particle image velocimetry results showed that the ground vortex mean streamwise position varied from 5 dn to 16 dn and the root mean square fluctuation in this position varied from 0.7 dn to 1.5 dn. Further analysis of the ground vortex temporal characteristics did not reveal any dominant non-dimensional frequencies.