Determination of complete velocity gradient tensor by using cinematographic stereoscopic PIV in a turbulent jet

Cinematographic stereoscopic PIV measurements were performed in the far field of an axisymmetric co-flowing turbulent round jet (Re _T ≈ 150, where Re _T is the Reynolds number based on Taylor micro scale) to resolve small and intermediate scales of turbulence. The time-resolved three-component PIV measurements were performed in a plane normal to the axis of the jet and the data were converted to quasi-instantaneous three-dimensional (volumetric) data by using Taylor’s hypothesis. The availability of the quasi-three-dimensional data enabled the computation of all nine components of the velocity gradient tensor over a volume. The use of Taylor’s hypothesis was validated by performing a separate set of time-resolved two component “side-view” PIV measurements in a plane along the jet axis. Probability density distributions of the velocity gradients computed using Taylor’s hypothesis show good agreement with those computed directly with the spatially resolved data. The overall spatial structure of the gradients computed directly exhibits excellent similarity with that computed using Taylor’s hypothesis. The accuracy of the velocity gradients computed from the pseudo-volume was assessed by computing the divergence error in the flow field. The root mean square (rms) of the divergence error relative to the magnitude of the velocity gradient tensor was found to be 0.25, which is consistent with results based on other gradient measurement techniques. The velocity gradients, vorticity components and mean dissipation in the self-similar far field of the jet were found to satisfy the axisymmetric isotropy conditions. The divergence error present in the data is attributed to the intrinsic uncertainty associated with performing stereoscopic PIV measurements and not to the use of Taylor’s hypothesis. The divergence error in the data is found to affect areas of low gradient values and manifests as nonphysical values for quantities like the normalized eigenvalues of the strain-rate tensor. However, the high gradients are less affected by the divergence error and so it can be inferred that structural features of regions of intense vorticity and dissipation will be faithfully rendered.     

Accuracy of tomographic particle image velocimetry data on a turbulent round jet

Tomographic particle image velocimetry (Tomo-PIV) was applied on a turbulent round air jet to quantitatively assess the accuracy of velocity gradients obtained in the self-similar turbulent region. The jet Reynolds number based on the nozzle diameter (d) was Re_d = 3000. Mean velocity, turbulent intensities, and Reynolds shear stress at the center plane of the jet were measured. In addition, statistical results of Tomo-PIV along the axial direction were assessed by performing a separate set of two-dimensional two-component PIV experiments on a “side view” plane along the jet axis. Moreover, the probability distribution functions of four components of the measured velocity gradients in the axial and radial directions were validated by these “side view” planar PIV data. The root mean square of the velocity divergence values relative to the norm of the velocity gradient tensor was 0.36. Furthermore, the on- and off-diagonal components of the velocity gradients satisfied the axisymmetric isotropy conditions. The divergence error in the data affected only areas with low gradient magnitude. Therefore, turbulent structures in the regions with intense vorticity and dissipation can be closely monitored. On this basis, the joint pdfs of the invariants of the velocity gradient and strain and rotation tensor rates were produced and compared well with those in isotropic turbulence studies.     

Stereoscopic particle image velocimetry

Stereoscopic particle image velocimetry (PIV) employs two cameras to record simultaneous but distinct off-axis views of the same region of interest (illuminated plane within a flow seeded with tracer particles). Sufficient information is contained in the two views to extract the out-of-plane motion of particles, and also to eliminate perspective error which can contaminate the in-plane measurement. This review discusses the principle of stereoscopic PIV, the different stereoscopic configurations that have been used, the relative error in the out-of-plane to the in-plane measurement, and the relative merits of calibration-based methods for reconstructing the three-dimensional displacement vector in comparison to geometric reconstruction. It appears that the current trend amongst practitioners of stereoscopic PIV is to use digital cameras to record the two views in the angular displacement configuration while incorporating the Scheimpflug condition. The use of calibration methods has also gained prominence over geometric reconstruction. Received: 15 April 1999/Accepted: 1 February 2000     

3-Component acceleration field measurement by dual-time stereoscopic particle image velocimetry

In this article, a multiplane stereo-particle image velocimetry (PIV) system was implemented and validated to measure the three-component acceleration field in a plane of turbulent flows. The employed technique relies on the use of two stereoscopic particle image velocimetry (SPIV) systems to measure pairs of velocity fields superimposed in space but shifted in time. The time delay between the two velocity fields enables the implementation of a finite difference scheme to compute temporal derivatives. The use of two synchronized SPIV systems allows us to overcome the limited acquisition rate of PIV systems when dealing with highly turbulent flows. Moreover, a methodology based on the analysis of the spectral error distribution is described here to determine the optimal time delay to compute time derivatives. The present dual-time SPIV arrangement and the proposed analysis method are applied to measure three-component acceleration fields in a cross section of a subsonic plane turbulent mixing layer.     

Dual-plane stereoscopic particle image velocimetry: system set-up and its application on a lobed jet mixing flow

 The technical basis and system set-up of a dual-plane stereoscopic particle image velocimetry (PIV) system, which can obtain the flow velocity (all three components) fields at two spatially separated planes simultaneously, is summarized. The simultaneous measurements were achieved by using two sets of double-pulsed Nd:Yag lasers with additional optics to illuminate the objective fluid flow with two orthogonally linearly polarized laser sheets at two spatially separated planes, as proposed by Kaehler and Kompenhans in 1999. The light scattered by the tracer particles illuminated by laser sheets with orthogonal linear polarization were separated by using polarizing beam-splitter cubes, then recorded by high-resolution CCD cameras. A three-dimensional in-situ calibration procedure was used to determine the relationships between the 2-D image planes and three-dimensional object fields for both position mapping and velocity three-component reconstruction. Unlike conventional two-component PIV systems or single-plane stereoscopic PIV systems, which can only get one-component of vorticity vectors, the present dual-plane stereoscopic PIV system can provide all the three components of the vorticity vectors and various auto-correlation and cross-correlation coefficients of flow variables instantaneously and simultaneously. The present dual-plane stereoscopic PIV system was applied to measure an air jet mixing flow exhausted from a lobed nozzle. Various vortex structures in the lobed jet mixing flow were revealed quantitatively and instantaneously. In order to evaluate the measurement accuracy of the present dual-plane stereoscopic PIV system, the measurement results were compared with the simultaneous measurement results of a laser Doppler velocimetry (LDV) system. It was found that both the instantaneous data and ensemble-averaged values of the stereoscopic PIV measurement results and the LDV measurement results agree well. For the ensemble-averaged values of the out-of-plane velocity component at comparison points, the differences between the stereoscopic PIV and LDV measurement results were found to be less than 2%. Received: 18 April 2000/Accepted: 2 February 2001     

Three-dimensional temporally resolved measurements of turbulence–flame interactions using orthogonal-plane cinema-stereoscopic PIV

A new orthogonal-plane cinema-stereoscopic particle image velocimetry (OPCS-PIV) diagnostic has been used to measure the dynamics of three-dimensional turbulence–flame interactions. The diagnostic employed two orthogonal PIV planes, with one aligned perpendicular and one aligned parallel to the streamwise flow direction. In the plane normal to the flow, temporally resolved slices of the nine-component velocity gradient tensor were determined using Taylor’s hypothesis. Volumetric reconstruction of the 3D turbulence was performed using these slices. The PIV plane parallel to the streamwise flow direction was then used to measure the evolution of the turbulence; the path and strength of 3D turbulent structures as they interacted with the flame were determined from their image in this second plane. Structures of both vorticity and strain-rate magnitude were extracted from the flow. The geometry of these structures agreed well with predictions from direct numerical simulations. The interaction of turbulent structures with the flame also was observed. In three dimensions, these interactions had complex geometries that could not be reflected in either planar measurements or simple flame–vortex configurations.     

Investigation of turbulence structures in a drag-reduced turbulent channel flow with surfactant additive by stereoscopic particle image velocimetry

In the present study, we employed stereoscopic particle image velocimetry (PIV) to investigate the characteristics of turbulence structures in a drag-reduced turbulent channel flow with addition of surfactant. The tested drag-reducing fluid was a CTAC/NaSal/Water (CTAC: cetyltrimethyl ammonium chloride; NaSal: sodium salicylate) system at 25°C. The weight concentration of CTAC was 30 ppm. Stereoscopic PIV measurement was performed for a water flow (Re=1.1×10⁴) and a CTAC solution flow (Re=1.5×10⁴ with 54% drag reduction) in both the streamwise–spanwise and wall-normal-spanwise planes, respectively. The three-dimensionality of hairpin vortex structures in the near-wall region for wall-bounded turbulent flow was reproduced by conditionally averaging the stereoscopic two-dimensional-three-component velocity fields. A series of wall-normal vortex cores were found to align with the near-wall low-speed streaks with opposite vorticity signals at both sides of the streaks and with the vorticity decreased on average by about one order of magnitude in CTAC solution flow compared with water flow; the spanwise spacing between the near-wall low-speed streaks in the solution flow is increased by about 46%. The streamwise vorticity of the vortex cores appearing in the wall-normal-spanwise plane was also decreased by the use of drag-reducing surfactant additives.     

A method for automatic estimation of instantaneous local uncertainty in particle image velocimetry measurements

The uncertainty of any measurement is the interval in which one believes the actual error lies. Particle image velocimetry (PIV) measurement error depends on the PIV algorithm used, a wide range of user inputs, flow characteristics, and the experimental setup. Since these factors vary in time and space, they lead to nonuniform error throughout the flow field. As such, a universal PIV uncertainty estimate is not adequate and can be misleading. This is of particular interest when PIV data are used for comparison with computational or experimental data. A method to estimate the uncertainty from sources detectable in the raw images and due to the PIV calculation of each individual velocity measurement is presented. The relationship between four error sources and their contribution to PIV error is first determined. The sources, or parameters, considered are particle image diameter, particle density, particle displacement, and velocity gradient, although this choice in parameters is arbitrary and may not be complete. This information provides a four-dimensional “uncertainty surface” specific to the PIV algorithm used. After PIV processing, our code “measures" the value of each of these parameters and estimates the velocity uncertainty due to the PIV algorithm for each vector in the flow field. The reliability of our methodology is validated using known flow fields so the actual error can be determined. Our analysis shows that, for most flows, the uncertainty distribution obtained using this method fits the confidence interval. An experiment is used to show that systematic uncertainties are accurately computed for a jet flow. The method is general and can be adapted to any PIV analysis, provided that the relevant error sources can be identified for a given experiment and the appropriate parameters can be quantified from the images obtained.     

Three-component velocity field measurement in confined liquid flows with high-speed digital image plane holography

In the last years, several techniques have been developed for the measurement of the three velocity components in a fluid plane or volume. Techniques as stereoscopic particle image velocimetry (SPIV) or tomographic PIV need a complex set-up and present serious restrictions when applied to confined liquid flows. Other like digital holographic PIV has some limitations in the particle concentration that can be measured. In this work, high-speed digital image plane holography has been applied for the measurement of the three velocity components in a complex geometry brain aneurysm model, using a two-cavity high-speed laser, one double frame camera and normal visualization, like in regular PIV. A portable and compact system has been built for adapting the high-speed laser short coherence length to the measurement of larger areas.     

Solid/free-surface juncture boundary layer and wake

 The Reynolds-averaged flow for a solid/free-surface juncture boundary layer and wake is documented. The three mean-velocity components and five of the Reynolds stresses are measured for a surface-piercing flat plate in a towing tank using a laser-Doppler velocimeter system for both boundary-layer and wake planes in regions close to the free surface. The experimental method is described, including the foil-plate model, laser-Doppler velocimeter system, conditions, and uncertainty analysis. The underlying flow data is in excellent agreement with benchmark data. Inner (near the plate and wake centerplane and below the free surface) and outer (near the free surface) regions of high streamwise vorticity of opposite sign are observed, which transport, respectively, high mean velocity and low turbulence from the outer to the inner and low mean velocity and high turbulence from the inner to the outer portions of the boundary layer and wake. For the wake, the inner region of vorticity is relatively weak. The physical mechanism for the streamwise vorticity is analyzed with regard to the Reynolds-averaged streamwise vorticity equation. The anisotropy of the crossplane normal Reynolds stresses closely correlates with the vorticity and, additionally, indicates similarity, i.e., its nature is such that it only depends on the proximity to the plate and free surface boundaries or wake centerplane symmetry plane. Free-surface effects on the Reynolds stresses are analyzed with regard to the behavior close to the free surface of the turbulent kinetic energy and the normal components of the anisotropy tensor and the anisotropy invariants. Close to the free surface, the turbulent kinetic energy is nearly constant and increases for the inner and outer portions, respectively, of the boundary layer and wake and the normal components of the anisotropy tensor and the anisotropy invariants roughly correspond to the limiting values for two-component turbulence. The similarities and differences between the present results and analysis with those from related studies are discussed. The data and analysis should have practical application with regard to the development of turbulence models for computational fluid dynamics methods for the Reynolds-averaged Navier–Stokes equations. Received: 27 May 1997/Accepted: 1 August 1997     

The influence of peak-locking errors on turbulence statistics computed from PIV ensembles

The influence of peak-locking errors on turbulence statistics computed from ensembles of PIV data is considered. PIV measurements are made in the streamwise–wall-normal plane of turbulent channel flow. The PIV images are interrogated in three distinct ways, generating ensembles of velocity fields with absolute, moderate, and minimal peak locking. Turbulence statistics computed for all three ensembles of data indicate a general sensitivity to peak locking in the single-point statistics, except for the mean velocity profile. Peak-locking errors propagate into the fluctuations of velocity, rendering single-point statistics inaccurate when severe peak locking is present. Multi-point correlations of both streamwise and wall-normal velocity are also found to be influenced by severe levels of peak locking. The displacement range of the measurement, defined by the PIV time delay, appears to affect the influence of peak-locking errors on turbulence statistics. Smaller displacement ranges, particularly those that produce displacement fluctuations that are less than one pixel in magnitude, yield inaccurate turbulence statistics in the presence of peak locking.     

Stereoscopic PIV: validation and application to an isotropic turbulent flow

 A new stereoscopic PIV system to measure the three velocity components is developed and applied to grid turbulence flows. This system uses two CCD cameras coupled with an accurate cross-correlation calculation method. An experimental test (based upon three-dimensional displacements) has been carried out to demonstrate the capability of this process to locate the maximum of correlation, and to detect accurately the 3D displacements. Experiments in a well-established turbulent flow have validated the method for quantitative measurements and a comparison with LDV results showed a good agreement in terms of mean and fluctuating velocities. Combined PIV and stereoscopic PIV measurements on a turbulent flow revealed the need to the stereoscopic systems to measure accurate 2D velocity fields. It has been shown that an error of up to 10% in the velocity fluctuation measured by conventional PIV could be attained due to 3D effects in highly turbulent cases. Finally, the digital cross-correlation technique adapted to the determination of small displacements seems to be the most suitable technique for stereoscopic PIV. Received: 22 July 1997/Accepted: 27 January 1998     

A kilohertz frame rate cinemagraphic PIV system for laboratory-scale turbulent and unsteady flows

A kilohertz frame rate cinemagraphic particle image velocimetry (PIV) system has been developed for acquiring time-resolved image sequences of laboratory-scale gas and liquid-phase turbulent flows. Up to 8000 instantaneous PIV images per second are obtained, with sequence lengths exceeding 4000 images. The two-frame cross-correlation method employed precludes directional ambiguity and has a higher signal-to-noise ratio than single-frame autocorrelation or cross-correlation methods, facilitating acquisition of long uninterrupted sequences of valid PIV images. Low and high velocities can be measured simultaneously with similar accuracy by adaptively cross-correlating images with the appropriate time delay. Seed particle illumination is provided by two frequency-doubled Nd:YAG lasers producing Q-switched pulses at the camera frame rate. PIV images are acquired using a 16 mm high-speed rotating prism camera. Frame-to-frame registration is accomplished by imaging two pairs of crossed lines onto each frame and aligning the digitized image sequence to these markers using image processing algorithms. No flow disturbance is created by the markers because only their image is projected to the PIV imaging plane, with the physical projection device residing outside the flow field. The frame-to-frame alignment uncertainty contributes 2% to the overall velocity measurement uncertainty, which is otherwise comparable to similar film-based PIV methods. Received: 11 July 2000 / Accepted: 21 June 2001 Published online: 29 November 2001     

Three-dimensional magnetic resonance velocimetry measurements of turbulence quantities in complex flow

A magnetic resonance velocimetry (MRV) experimental technique based on magnetic resonance imaging and capable of measuring the turbulent Reynolds stresses in a 3D flow domain is described. Results are presented in backward facing step flow in a square channel with a Reynolds number of 48,000 based on step height and freestream velocity at the step. MRV results are compared to particle image velocimetry (PIV) measurements in the centerplane containing the streamwise and cross-stream axes. MRV and PIV mean velocity measurements show excellent agreement. MRV measurements for Reynolds normal stresses compare to within ±20% of the PIV results while results for the turbulent shear are less accurate.     

An experimental study of a turbulent vortex ring: a three-dimensional representation

This paper presents a reconstruction of the three-dimensional velocity field of a turbulent vortex ring by means of Taylor’s hypothesis. Stereoscopic PIV is used to acquire three velocity component information on a plane. The accuracy of the Taylor’s hypothesis for this particular flow pattern is first discussed, and the three-dimensional velocity and vorticity information are then presented. This study also introduces an azimuthally averaging method in order to give a mean structure in cylindrical coordinates from a single realization and from which turbulent stresses and production can be estimated. The azimuthally averaged quantities are then compared with the ensemble-averaged results from the previous planar (two-dimensional and stereoscopic) PIV experiments.     

Three-component micro-PIV using the continuity equation and a comparison of the performance with that of stereoscopic measurements

Recently, a number of techniques have been presented for the determination of the third “out-of-plane” velocity component in micro particle image velocimetry (micro-PIV) data. In particular, the conventional macroscopic stereo-PIV technique has been converted to the micro scale by the use of stereo-microscopy. In this work a different technique is investigated, which uses conventional, two-component micro-PIV to generate velocity data on a number of planes. The in-plane velocity gradients are then calculated, which can be used in the continuity equation to produce the out-of-plane velocity gradients. These, together with the no-penetration boundary condition, can then be used to calculate the out-of-plane velocities. An algorithm is presented that is capable of handling up to one invalid vector per column of data by using a combination of first order and second order projections of the velocity. The advantage of the continuity based technique is that it uses the existing two-component micro-PIV technology, which at present is in a more advanced stage of development then stereo-microscopy based micro-PIV. The technique is investigated using a flow similar to one used previously to assess stereoscopic micro-PIV (Meas Sci Technol 17:2175–2185, 2006). This allows a comparison of the performance of the two techniques. The results show that the continuity based data agrees well with an independent computational fluid dynamics solution and has a smaller experimental uncertainty than the stereoscopic technique at a better spatial resolution. There are, however, potential limitations to the continuity based technique. These include the two-dimensionality of the data, which is limited to the planes on which the original images were taken, and the dependence of the technique on the data close to surfaces, where the experimental errors are often greatest. Stereoscopic micro-PIV does not have these limitations so, whilst at present it appears that continuity based techniques may be more accurate, there is sufficient potential for stereoscopic techniques to justify their further development.     

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.     

PIV measurements of turbulent flow in planar mixing layer

A turbulent mixing layer consists of two different flow types, i.e. shear layer (shear-flow turbulence) and free stream regions (nearly homogeneous turbulence). The inherent non-uniform seeding tracer distributions observed around the interfaces between the shear layer and two free stream regions usually lead to a difficulty in particle image velocimetry (PIV) measurements. A parametric study on the application of PIV to the measurement of velocity field in a planar mixing layer is made by means of six factors, including interrogation window size, aspect ratio of interrogation window, interrogation window offset, threshold of data validation, sharpening spatial filters (Prewitt and Sobel masks), and smoothing spatial filter (median mask). The objective of this study is to obtain accurate turbulent measurements in both mean and fluctuating velocities using PIV under an appropriate parametric setting. The optimal levels, which are trade-off in between the accuracy and fine spatial resolution of velocity field measurements, are determined with the aid of the Taguchi method. It is shown that the PIV measurements made with this optimal set of parameters are in good agreement with the measurements made by a two-component hot-wire anemometer. Case independency of the proposed optimal set of parameters on the flow condition of the mixing layer is validated through the applications to two additional tests under the different experimental conditions in changing solely either velocity ratio of high-speed to low-speed free stream velocities or Reynolds number.     

Investigation of the flow in a circular cavity using stereo and tomographic particle image velocimetry

The turbulent flow over a circular cavity with an aspect ratio of D/H = 2 is investigated by multi-planar stereoscopic particle image velocimetry and with tomographic particle image velocimetry (PIV). The main aim of the study is the flow topology and the turbulent structure of the asymmetrical flow pattern that forms inside the cavity at these specific conditions. The flow field is measured in the vertical symmetry plane to describe the overall recirculation pattern in the cavity and the turbulent shear layer developing from the separation point. In this specific regime the shear layer fluctuations are recognized as those caused by instabilities together with the effect of the incoming boundary layer turbulence. Additional observations performed at several wall-parallel planes at different height inside the cavity allow to further evaluate the secondary flow circulation generated by this asymmetric regime. The observed flow pattern consists of a steady vortex, occupying the entire cavity volume and placed diagonally inside the cavity such to entrain the external flow from one side, capture it into a circulatory motion and eject it from the opposite side of the cavity. The spatial distribution of the turbulent fluctuations also reveals the same structure. The tomographic PIV measurement returns a visual inspection to the instantaneous three-dimensional structure of the turbulent fluctuations, which at the investigated height exhibit a low level of coherence with slightly elongated vortices in the recirculating flow inside the cavity.