Large-scale tomographic PIV in forced and mixed convection using a parallel SMART version

Large-scale tomographic particle image velocimetry (tomographic PIV) was used to study large-scale flow structures of turbulent convective air flow in an elongated rectangular convection cell. Three flow cases have been investigated, that is, pure forced convection and mixed convection at two different Archimedes numbers. The Reynolds number was constant at Re?=?1.04?×?10⁴ for all cases, while the Archimedes numbers were Ar?=?2.1 and 3.6 for the mixed convection cases, corresponding to Rayleigh numbers of Ra?=?1.6?×?10⁸ and 2.8?×?10⁸, respectively. In these investigations, the size of the measurement volume was as large as 840?mm?×?500?mm?×?240?mm. To allow for statistical analysis of the measured instantaneous flow fields, a large number of samples needed to be evaluated. Therefore, an efficient parallel implementation of the tomographic PIV algorithm was developed, which is based on a version of the simultaneous multiplicative reconstruction technique (SMART). Our algorithm distinguishes itself amongst other features by the fact that it does not store any weighting coefficients. The measurement of forced convection reveals an almost two-dimensional roll structure, which is orientated in the longitudinal cell direction. Its mean velocity field exhibits a core line with a wavy shape and a wavelength, which corresponds to the height and depth of the cell. In the instantaneous fields, the core line oscillates around its mean position. Under the influence of thermal buoyancy forces, the global structure of the flow field changes significantly. At lower Archimedes numbers, the resulting roll-like structure is shifted and deformed as compared to pure forced convection. Additionally, the core line oscillates much more strongly around its mean position due to the interaction of the roll structure with the rising hot air. If the Archimedes number is further increased, the roll-like structure breaks up into four counter-rotating convection rolls as a result of the increased influence of buoyancy forces. Moreover, large-scale tomographic PIV reveals that the orientation of these rolls reflects a ??W??-like shape in the horizontal X?CZ-plane of the convection cell.     

Tomographic particle image velocimetry investigation of the flow in a modeled human carotid artery bifurcation

Hemodynamic forces within the human carotid artery are well known to play a key role in the initiation and progression of vascular diseases such as atherosclerosis. The degree and extent of the disease largely depends on the prevailing three-dimensional flow structure and wall shear stress (WSS) distribution. This work presents tomographic PIV (Tomo-PIV) measurements of the flow structure and WSS in a physiologically accurate model of the human carotid artery bifurcation. The vascular geometry is reconstructed from patient-specific data and reproduced in a transparent flow phantom to demonstrate the feasibility of Tomo-PIV in a complex three-dimensional geometry. Tomographic reconstruction is performed with the multiplicative line-of-sight (MLOS) estimation and simultaneous multiplicative algebraic reconstruction (SMART) technique. The implemented methodology is validated by comparing the results with Stereo-PIV measurements in the same facility. Using a steady flow assumption, the measurement error and RMS uncertainty are directly inferred from the measured velocity field. It is shown that the measurement uncertainty increases for increasing light sheet thickness and increasing velocity gradients, which are largest near the vessel walls. For a typical volume depth of 6 mm (or 256 pixel), the analysis indicates that the velocity derived from 3D cross-correlation can be measured within ±2% of the maximum velocity (or ±0.2 pixel) near the center of the vessel and within ±5% (±0.6 pixel) near the vessel wall. The technique is then applied to acquire 3D-3C velocity field data at multiple axial locations within the carotid artery model, which are combined to yield the flow field and WSS in a volume of approximately 26 mm × 27 mm × 60 mm. Shear stress is computed from the velocity gradient tensor and a method for inferring the WSS distribution on the vessel wall is presented. The results indicate the presence of a complex and three-dimensional flow structure, with regions of flow separation and strong velocity gradients. The WSS distribution is markedly asymmetric confirming a complex swirling flow structure within the vessel. A comparison of the measured WSS with Stereo-PIV data returns an acceptable agreement with some differences in stress magnitude.     

Tomographic and time resolved PIV measurements on a finite cylinder mounted on a flat plate

Tomographic and time resolved PIV measurements were performed to examine the 3D flow topology and the flow dynamic above the upper surface of a low-aspect ratio cylinder at Re ≈ 1 ×  10⁵. This generic experiment is of fundamental interest because it represents flow features which are relevant to many applications such as laminar separation bubbles and turbulent reattachment. At Re  ≈ 1 × 10⁵, laminar separation bubbles arise on the side of the cylinder. Furthermore, on the top of the cylinder a separation with reattachment is of major interest. The tomographic PIV measurement, which allows to determine all three velocity components in a volume instantaneously, was applied to examine the flow topology and interaction between the boundary layer and wake structures on the top of the finite cylinder. In the instantaneous flow fields the tip vortices and the recirculation region becomes visible. However, it is also observed that the flow is quite unsteady due to the large separation occurring on the top of the cylinder. In order to study the temporal behaviour of the separation, time resolved PIV was applied. This technique allows capturing the dynamic processes in detail. The development of vortices in the separated shear layer is observed and in addition regions with different dominant frequencies are identified.     

Spatial resolution and dissipation rate estimation in Taylor--Couette flow for tomographic PIV

This paper assesses the spatial resolution and accuracy of tomographic particle image velocimetry (PIV). In tomographic PIV the number of velocity vectors are of the order of the number of reconstructed particle images, and sometimes even exceeds this number when a high overlap fraction between adjacent interrogations is used. This raises the question of the actual spatial resolution of tomographic PIV in relation to the various flow scales. We use a Taylor--Couette flow of a fluid between two independently rotating cylinders and consider three flow regimes: laminar flow, Taylor vortex flow and fully turbulent flow. The laminar flow has no flow structures, and the measurement results are used to assess the measurement uncertainty and to validate the accuracy of the technique for measurements through the curved wall. In the Taylor vortex flow regime, the flow contains large-scale flow structures that are much larger than the size of the interrogation volumes and are fully resolved. The turbulent flow regime contains a range of flow scales. Measurements in the turbulent flow regime are carried out for a Reynolds number Re between 3,800 and 47,000. We use the measured torque on the cylinders to obtain an independent estimate of the energy dissipation rate and estimate of the Kolmogorov length scale. The data obtained by tomographic PIV are assessed by estimating the dissipation rate and comparing the result against the dissipation rate obtained from the measured torque. The turbulent flow data are evaluated for different sizes of the interrogation volumes and for different overlap ratios between adjacent interrogation locations. The results indicate that the turbulent flow measurements for the lowest Re could be (nearly) fully resolved. At the highest Re only a small fraction of the dissipation rate is resolved, still a reasonable estimate of the total dissipation rate could be obtained by means of using a sub-grid turbulence model. The resolution of tomographic PIV in these measurements is determined by the size of the interrogation volume. We propose a range of vector spacing for fully resolving the turbulent flow scales. It is noted that the use of a high overlap ratio, that is, 75?%, yields a substantial improvement for the estimation of the dissipation rate in comparison with data for 0 and 50?% overlap. This indicates that additional information on small-scale velocity gradients can be obtained by reducing the data spacing.     

Development and validation of echo PIV

The combination of ultrasound echo images with digital particle image velocimetry (DPIV) methods has resulted in a two-dimensional, two-component velocity field measurement technique appropriate for opaque flow conditions including blood flow in clinical applications. Advanced PIV processing algorithms including an iterative scheme and window offsetting were used to increase the spatial resolution of the velocity measurement to a maximum of 1.8 mm×3.1 mm. Velocity validation tests in fully developed laminar pipe flow showed good agreement with both optical PIV measurements and the expected parabolic profile. A dynamic range of 1 to 60 cm/s has been obtained to date.     

Characterization of a system generating a homogeneous isotropic turbulence field by free synthetic jets

A facility inspired by Hwang and Eaton (2004a, b) for generating a homogeneous isotropic turbulence was built, the objective being to study evaporating droplets in the presence of turbulence. Turbulence was produced by the mixing of six synthetic jets, in ambient atmosphere. Combined PIV and LDA techniques were used to measure the statistical turbulence properties. The turbulence produced was found to be homogeneous isotropic with a small mean flow within a domain having an average size of 50 mm × 50 mm × 50 mm. The rms fluctuations were of the order of 0.9 m/s, corresponding to a Taylor Reynolds number of 240 and an integral length scale of about 40 mm. This apparatus proved to be well suited to the study of the evaporation of droplets in a controlled turbulence field.     

Visual hull method for tomographic PIV measurement of flow around moving objects

Tomographic particle image velocimetry (PIV) is a recently developed method to measure three components of velocity within a volumetric space. We present a visual hull technique that automates identification and masking of discrete objects within the measurement volume, and we apply existing tomographic PIV reconstruction software to measure the velocity surrounding the objects. The technique is demonstrated by considering flow around falling bodies of different shape with Reynolds number?~1,000. Acquired image sets are processed using separate routines to reconstruct both the volumetric mask around the object and the surrounding tracer particles. After particle reconstruction, the reconstructed object mask is used to remove any ghost particles that otherwise appear within the object volume. Velocity vectors corresponding with fluid motion can then be determined up to the boundary of the visual hull without being contaminated or affected by the neighboring object velocity. Although the visual hull method is not meant for precise tracking of objects, the reconstructed object volumes nevertheless can be used to estimate the object location and orientation at each time step.     

Flow field investigation in rotating rib-roughened channel by means of particle image velocimetry

The turbulent velocity field over the rib-roughened wall of an orthogonally rotating channel is investigated by means of two-dimensional particle image velocimetry (PIV). The flow direction is outward, with a bulk Reynolds number of 1.5 × 10⁴ and a rotation number ranging from 0.3 to 0.38. The measurements are obtained along the wall-normal/streamwise plane at mid-span. The PIV system rotates with the channel, allowing to measure directly the relative flow velocity with high spatial resolution. Coriolis forces affect the stability of the boundary layer and free shear layer. Due to the different levels of shear layer entrainment, the reattachment point is moved downstream (upstream) under stabilizing (destabilizing) rotation, with respect to the stationary case. Further increase in rotation number pushes further the reattachment point in stabilizing rotation, but does not change the recirculation length in destabilizing rotation. Turbulent activity is inhibited along the leading wall, both in the boundary layer and in the separated shear layer; the opposite is true along the trailing wall. Coriolis forces affect indirectly the production of turbulent kinetic energy via the Reynolds shear stresses and the mean shear. Two-point correlation is used to characterize the coherent motion of the separated shear layer. Destabilizing rotation is found to promote large-scale coherent motions and accordingly leads to larger integral length scales; on the other hand, the spanwise vortices created in the separating shear layer downstream of the rib are less organized and tend to be disrupted by the three-dimensional turbulence promoted by the rotation. The latter observation is consistent with the distributions of span-wise vortices detected in instantaneous flow realizations.     

Pulsed, high-power LED illumination for tomographic particle image velocimetry

This paper investigates the use of high-power light-emitting diode (LED) illumination for tomographic particle image velocimetry (PIV) as an alternative to traditional laser-based illumination. Modern solid-state LED devices can provide averaged radiant power in excess of 10 W and by operating the LED with short high current pulses theoretical pulse energies up to several tens of mJ can be achieved. In the present work, a custom-built drive circuit is used to drive a Luminus PT-120 high-power LED at pulsed currents of up to 150 A and 1 μs duration. Volumetric illumination is achieved by directly projecting the LED into the flow to produce a measurement volume of ≈3–4 times the size of the LED die. The feasibility of the volumetric LED illumination is assessed by performing tomographic PIV of homogenous, grid-generated turbulence. Two types of LEDs are investigated, and the results are compared with measurements of the same flow using pulsed Nd:YAG laser illumination and DNS data of homogeneous isotropic turbulence. The quality of the results is similar for both investigated LEDs with no significant difference between the LED and Nd:YAG illumination. Compared with the DNS, some differences are observed in the power spectra and the probability distributions of the fluctuating velocity and velocity gradients. These differences are attributed to the limited spatial resolution of the experiments and noise introduced during the tomographic reconstruction (i.e. ghost particles). The uncertainty in the velocity measurements associated with the LED illumination is estimated to approximately 0.2–0.3 pixel for both LEDs, which compares favourably with similar tomographic PIV measurements of turbulent flows. In conclusion, the proposed high-power, pulsed LED volume illumination provides accurate and reliable tomographic PIV measurements in water and presents a promising technique for flow diagnostics and velocimetry.     

Volumetric correlation PIV: a new technique for 3D velocity vector field measurement

A method is proposed that allows three-dimensional (3D) two-component measurements to be made by means of particle image velocimetry (PIV) in any volume illuminated over a finite thickness. The method is based on decomposing the cross-correlation function into various contributions at different depths. Because the technique is based on 3D decomposition of the correlation function and not reconstruction of particle images, there is no limit to particle seeding density as experienced by 3D particle tracking algorithms such as defocusing PIV and tomographic PIV. Correlations from different depths are differentiated by the variation in point spread function of the lens used to image the measurement volume over that range of depths. A number of examples are demonstrated by use of synthetic images which simulate micro-PIV (μPIV) experiments. These examples vary from the trivial case of Couette flow (linear variation of one velocity component over depth) to a general case where both velocity components vary by different complex functions over the depth. A final validation—the measurement of a parabolic velocity profile over the depth of a microchannel flow—is presented. The same method could also be applied using a thick light sheet in macro-scale PIV and in a stereo configuration for 3D three-component PIV.     

Turbulent free convection in a vertical converging channel

This paper presents the results of experimental and numerical investigations of the problem of turbulent natural convection in a converging-plate vertical channel. The channel has two isothermally heated inclined walls and two adiabatic vertical side walls. The parameters involved in this study are the channel geometry represented by the channel width at exit, the inclination of the heated walls and the temperature difference between the heated walls and the ambient. The investigation covered modified Rayleigh numbers up to 10⁸ in the computational study and up to 9.3 × 10⁶ in the experimental work. The experimental measurements focused on the velocity field and were carried out using a PIV system and included measurements of the mean velocity profiles as well as the root-mean-square velocity and shear stress profiles. The experiments were conducted for an inclination angle of 30°, a gap width of 10 mm and two temperature differences (∆T=25.4°C and 49.8°C). The velocity profiles in the lower part of the channel indicated the presence of two distinct layers. The first layer is adjacent to the heated plate and driven by buoyancy forces while the second layer extends from the point of maximum velocity to the channel center plane and driven mainly by shear forces. The velocity profile at the upper portion of the channel has shown the merging of the two boundary layers growing over the two heated walls. The measured values of the Reynolds shear stress and root mean square of the horizontal and vertical velocity fluctuation components have reached their maximum near the wall while having smaller values in the core region. The computational results have shown that the average Nusselt number increases approximately linearly with the increase of the modified Rayleigh number when plotted on log–log scale. The variation of the local Nusselt number indicated infinite values at the channel inlet (leading edge effect) and high values at the channel exit (trailing edge effect). For a fixed value of the top channel opening, the increase of the inclination angle tended to reduce flow velocity at the inlet section while changing the flow structure near the heated plates in such a way to create boundary-layer type flow. The maximum value of the average Nusselt number occurs when θ = 0 and decreases with the increase of the inclination angle. On the other hand, the increase of the channel width at exit for the same inclination angle caused a monotonic increase in the flow velocity at the channel inlet.     

Investigation of thermal convection in water columns using particle image velocimetry

Particle image velocity measurements were applied on thermally driven convection at low Rayleigh numbers. In a model experiment using a water column heated from bottom and cooled from above, the velocity field was studied at different vertical temperature gradients. In the testing facility with high aspect ratio (about 19) representing a 1-m-long column with 5?cm diameter, occurrence of free convection was verified for destabilizing temperature gradients of 0.1–2?K/m. The PIV results revealed that significant flow exists already at low vertical temperature gradients. The velocity of the stable large-scale circulations increased linearly with temperature gradient (<1?K/m) from 8?×?10^?5 to 1?×?10^?3?m/s. At higher temperature gradients (1–2?K/m), a transition from quasi-stationary into time-dependent flow was observed, where convection cells changed position, number, and form temporarily. The motivation of this research was to gain more insight into density-driven convection in boreholes and groundwater monitoring wells.     

Scanning liquid-crystal thermometry and stereo velocimetry for simultaneous three-dimensional measurement of temperature and velocity field in a turbulent Rayleigh-Bérnard convection

In this paper, we report on an experimental technique for the simultaneous measurement of temperature and three components of velocity in a three-dimensional thermal flow using scanning liquid-crystal thermometry and stereo velocimetry. The temperature is measured by the color image analysis of the liquid-crystal particles suspended in a fluid, while the three velocity components are measured by stereo particle image velocimetry (stereo PIV) with the aid of tracer particles. The measurement is carried out by scanning the light-sheet plane while capturing the sequential color images of the liquid crystals and tracer particles. This measurement allows the reconstruction of the three-dimensional distribution of temperature and full velocity field simultaneously. The present experimental technique is applied to the horizontal fluid layer of a turbulent Rayleigh-Bérnard convection and the three-dimensional structures of thermal plumes are evaluated. The experimental results indicate that the structures of plumes are often correlated with the vertical velocity of the fluid, but they behave randomly in space, influenced by the large-scale turbulence evident in the middle of the fluid layer.     

Measurements of a supersonic turbulent boundary layer by focusing schlieren deflectometry

 Some novel, non-intrusive, high-frequency, localized optical measurements of turbulence in compressible flows are described. The technique is based upon focusing schlieren optics coupled with high-speed quantitative measurement of light intensity fluctuations in the schlieren image. Measurements of density gradient fluctuations confined to a thin slice of the flowfield are thus obtained. The new instrument was used to investigate the structure of a two-dimensional, adiabatic, wind tunnel wall boundary layer at a Mach number of 3. The measurements were compared to data obtained using hot-wire anemometry and good agreement was found between the two. Distributions of broadband convection velocity of large-scale structures through the boundary later were also measured. In marked contrast to earlier results, it is shown here that the convection velocity is essentially identical to the local mean velocity. Further, results obtained using the VITA conditional sampling technique shed new light on the turbulent boundary layer structure. Overall, the data presented herein serve to validate the new measurement technique. Received: 12 February 1997/Accepted: 31 January 1998     

Planar velocity measurements of a two-dimensional compressible wake

The present study describes the application of particle image velocimetry (PIV) to investigate the compressible flow in the wake of a two-dimensional blunt base at a freestream Mach number M_X=2. The first part of the study addresses specific issues related to the application of PIV to supersonic wind tunnel flows, such as the seeding particle flow-tracing fidelity and the measurement spatial resolution. The seeding particle response is assessed through a planar oblique shock wave experiment. The measurement spatial resolution is enhanced by means of an advanced image-interrogation algorithm. In the second part, the experimental results are presented. The PIV measurements yield the spatial distribution of mean velocity and turbulence. The mean velocity distribution clearly reveals the main flow features such as expansion fans, separated shear layers, flow recirculation, reattachment, recompression and wake development. The turbulence distribution shows the growth of turbulent fluctuations in the separated shear layers up to the reattachment location. Increased velocity fluctuations are also present downstream of reattachment outside of the wake due to unsteady flow reattachment and recompression. The instantaneous velocity field is analyzed seeking coherent flow structures in the redeveloping wake. The instantaneous planar velocity and vorticity measurements return evidence of large-scale turbulent structures detected as spatially coherent vorticity fluctuations. The velocity pattern consistently shows large masses of fluid in vortical motion. The overall instantaneous wake flow is organized as a double row of counter-rotating structures. The single structures show vorticity contours of roughly elliptical shape in agreement with previous studies based on spatial correlation of planar light scattering. Peak vorticity is found to be five times higher than the mean vorticity value, suggesting that wake turbulence is dominated by the activity of large-scale structures. The unsteady behavior of the reattachment phenomenon is studied. Based on the instantaneous flow topology, the reattachment is observed to fluctuate mostly in the streamwise direction suggesting that the unsteady separation is dominated by a pumping-like motion.     

Particle image velocimetry investigation of the flow-field of a 3D turbulent annular swirling decaying flow induced by means of a tangential inlet

Particle image velocimetry (PIV) has been used in order to measure the three mean components and turbulence intensities of the velocity vector in a swirling decaying flow induced by a tangential inlet in an annulus. This kind of flow motion is found to be very complex, exhibiting three-dimensional and non-axisymmetric characteristics coupled with a free decay of the swirling intensity along the flow path, thereby making it difficult to study. A method allowing the measurement of the three components of the velocity flow-field with a standard PIV system with two-dimensional acquisitions, is presented. The evolution of each velocity component between the inlet and the outlet of the annulus is obtained. Furthermore, the PIV technique is extended to the measurement of turbulent characteristics such as turbulent intensities and dimensionless turbulent energy. The main characteristics of the swirling flow are discussed and the swirl number is estimated as a function of the axial distance from the tangential inlet. Received: 6 July 1998/Accepted: 20 March 1999     

Volume self-calibration for 3D particle image velocimetry

Planar self-calibration methods have become standard for stereo PIV to correct misalignments between laser light sheet and calibration plane. Computing cross-correlation between images from camera 1 and 2 taken at the same time, non-zero disparity vectors indicate rotational and translational misalignments relative to the coordinate system defined by a calibration plate. This approach works well for thin light sheets but fails for extended volumes recorded in 3D-PTV or tomographic PIV experiments. Here it is primarily necessary to correct calibration errors leading to triangulation errors in 3D-PTV or in degraded tomographic volume reconstruction. Tomographic PIV requires calibration accuracies of a fraction of a pixel throughout the complete volume, which is difficult to achieve experimentally. A new volumetric self-calibration technique has been developed based on the computation of the 3D position of matching particles by triangulation as in 3D-PTV. The residual triangulation error (‘disparity’) is then used to correct the mapping functions for all cameras. A statistical clustering method suitable for dense particle images has been implemented to find correct disparity map peaks from true particle matches. Disparity maps from multiple recordings are summed for better statistics. This self-calibration scheme has been validated using several tomographic PIV experiments improving the vector quality significantly. The relevance for other 3D velocimetry methods is discussed.     

Stereoscopic and tomographic PIV of a pitching plate

This paper applies particle image velocimetry (PIV) to a simplified, canonical, pitch-hold-return problem of a pitching plate in order to gain some understanding of how three dimensionality develops in such flows. Data from a progression of PIV studies, from stereoscopic PIV yielding three-component, two-dimensional (3C-2D) data to tomographic PIV yielding three-component, three-dimensional (3C-3D) data are presented thus providing progressively more detailed information. A comparison of results is made between the two techniques. The PIV study is performed in a water tunnel facility with cross-sectional area 500 × 500 mm, and involves a full-span (nominally two-dimensional) plate, suspended between a wall end boundary condition and a free surface, pitching at a dimensionless pitch rate of K _c  = 0.93 in flow at Re = 7,500. Results demonstrate the existence of spanwise flows in both the leading edge and trailing edge vortices, but with strong directionality in the leading edge vortex towards the wall end boundary condition. Observations of instantaneous flow patterns suggest also the existence of three-dimensional coherent vortex filament structures in the outer regions of the leading edge vortex.     

A scanning PIV method for fine-scale turbulence measurements

A hybrid technique is presented that combines scanning PIV with tomographic reconstruction to make spatially and temporally resolved measurements of the fine-scale motions in turbulent flows. The technique uses one or two high-speed cameras to record particle images as a laser sheet is rapidly traversed across a measurement volume. This is combined with a fast method for tomographic reconstruction of the particle field for use in conjunction with PIV cross-correlation. The method was tested numerically using DNS data and with experiments in a large mixing tank that produces axisymmetric homogeneous turbulence at \(R_\lambda \simeq 219\) . A parametric investigation identifies the important parameters for a scanning PIV set-up and provides guidance to the interested experimentalist in achieving the best accuracy. Optimal sheet spacings and thicknesses are reported, and it was found that accurate results could be obtained at quite low scanning speeds. The two-camera method is the most robust to noise, permitting accurate measurements of the velocity gradients and direct determination of the dissipation rate.