PIV measurement of flow around an arbitrarily moving body

This paper presents a PIV (particle image velocimetry) image processing method for measuring flow velocities around an arbitrarily moving body. This image processing technique uses a contour-texture analysis based on user-defined textons to determine the arbitrarily moving interface in the particle images. After the interface tracking procedure is performed, the particle images near the interface are transformed into Cartesian coordinates that are related to the distance from the interface. This transformed image always has a straight interface, so the interrogation windows can easily be arranged at certain distances from the interface. Accurate measurements near the interface can then be achieved by applying the window deformation algorithm in concert with PIV/IG (interface gradiometry). The displacement of each window is evaluated by using the window deformation algorithm and was found to result in acceptable errors except for the border windows. Quantitative evaluations of this method were performed by applying it to computer-generated images and actual PIV measurements.     

Theory of non-isotropic spatial resolution in PIV

The spatial resolution of the PIV interrogation technique is discussed from an analytical standpoint and assessed with Monte Carlo numerical simulation of particle image motion. The PIV measurement error associated with lack of spatial resolution is modelled associating the cross-correlation operator to a moving average filter. The error associated with the "low-pass filtering" effect is investigated by adopting a second-order polynomial expression for the velocity spatial distribution. According to the present error analysis, the measurement error is proportional to the second-order spatial derivative of the velocity field and increases with the square of the window linear size. The strategy for the selection of the window size and properties (aspect ratio and orientation) so as to minimize the error is discussed. The principle is based on nonisotropic interrogation windows of elliptical shape, with a constant area and elongated in the direction of the largest curvature radius. The nonisotropic parameters are defined as eccentricity and orientation, which are based on the local eigenvalues/vectors of the Hessian tensor of the displacement spatial distribution. The technique is implemented in a recursive PIV interrogation method. The performance of nonisotropic interrogation technique is assessed by means of synthetic PIV images, which simulate three situations: first, a one-dimensional sinusoidal shear displacement, which allows comparison of the cross-correlation spatial response with the transfer function of linear filters. Second, the stream-wise exponential velocity decay is simulated, which simulates the particle tracers decelerating downstream of a shock wave and gives an example of a flow with main velocity differences aligned with the velocity direction. The results show that keeping the image density fixed, the error caused by insufficient spatial resolution can be reduced by a factor two when a preferential direction is found in the flow field. Finally, a Lamb–Oseen vortex flow is presented, which shows the complex pattern formed by the interrogation windows in a two-dimensional case. In this case, the improvement in interrogation performance is limited due to the isotropic nature of the velocity spatial fluctuation.     

A novel architecture for a super-resolution PIV algorithm developed for the improvement of the resolution of large velocity gradient measurements

Three different particle image processing algorithms have been developed for the improvement of PIV velocity measurements characterized by large velocity gradients. The objectives of this study are to point out the limitations of the standard processing methods and to propose a complete algorithm to enhance the measurement accuracy. The heart of the PIV image processing is a direct cross-correlation calculation in order to obtain complete flexibility in the choice of the size and the shape of the interrogation window (IW). An iterative procedure is then applied for the reduction of the size of IW at each measurement location. This procedure allows taking into account the local particle concentration in the image. The results of this first iterative processing, applied to synthetic images, show both a significant improvement of measurement accuracy and an increase of the spatial resolution. Finally, a super-resolution algorithm is developed to further increase the spatial resolution of the measurement by determining the displacement of each particle. The computer time for a complete image processing is optimized by the introduction of original data storage in Binary Space Partitions trees. It is shown that measurement errors for large velocity gradient flows are similar to those obtained in simpler cases with uniform translation displacements. This last result validates the ability of the developed super-resolution algorithm for the aerodynamic characterization of large velocity gradient flows.     

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.     

A comparative study of the MQD method and several correlation-based PIV evaluation algorithms

 The Minimum Quadratic Difference (MQD) method is compared with methods conventionally used for the evaluation of PIV recordings, i.e. correlation-based evaluation with fixed interrogation windows (auto- or cross-correlation) and correlation-based tracking. The comparison is performed by studying the evaluation accuracy achieved when applying these methods to pairs of synthetic PIV recordings for which the true displacements are known. The influence of the magnitude of the particle image displacement, evaluation window size, density of particle image distribution, and particle image size on the accuracy are investigated. In all these cases the best results in terms of a statistical error are obtained with the MQD method. The superiority of the MQD method can be explained with its potential of accounting for non-uniformities in the particle image distribution and a non-uniform illumination. It is also shown that the conventional correlation-based methods may produce principal errors that are non-existent for the MQD method. The evaluation speed achievable for the MQD method by making use of the FFT is comparable to that common for the generally used auto- or cross-correlation algorithm. Finally, a quantitative explanation is given for the often observed phenomenon that PIV velocity results tend to be smaller than the true values. Received: 15 May 1998/Accepted: 24 April 1999     

Improvements to PIV image analysis by recognizing the velocity gradients

Two iterative PIV image processing methods are introduced, which utilize displacement and deformation of the interrogation areas to maximize the correlation. The velocity gradients used for the window deformation are iteratively estimated directly from the images and no velocity values are required from neighbouring interrogation areas, as with numerical differentiation. The improved accuracy and resolution of the velocity gradient estimation compared to numerical differentiation is shown using synthetic images. The performance in a real application is shown using experimental reference images.     

A generalized processing technique in digital particle image velocimetry with direct estimation of velocity gradients

A technique is proposed for the processing of digital particle image velocimetry (PIV) images, in one single step providing direct estimates of fluid velocity, out-of-plane vorticity and in-plane shear rate tensor. The method is based on a generalization of the standard PIV cross-correlation technique and substitutes the usual discrete cross-correlation of image pairs with a correlation of interpolated two-dimensional image intensity functions, being subject to affine transformations. The correlation is implemented by using collocation points, on which image intensity values are interpolated. The resulting six-dimensional correlation function is maximized using a general purpose optimization algorithm. The use of the method is demonstrated by application to different types of synthetically generated image pairs constructed with known particle displacement functions. The resulting errors are assessed and compared with those of a representative standard PIV method as well as with those of the present technique using no differential quantities in the search of the peak location. The examples demonstrate that significant improvements in accuracy can be obtained for flow fields with regions containing strong velocity gradients.     

A vision-based hybrid particle tracking velocimetry (PTV) technique using a modified cascade correlation peak-finding method

A novel technique for particle tracking velocimetry is presented in this paper to overcome the issue of overlapping particle images encountered in the flows with high particle density or under volumetric illumination conditions. To achieve this goal, algorithms for particle identification and tracking are developed based on current methods and validated with both synthetic and experimental image sets. The results from synthetic image tests show that the particle identification algorithm is able to resolve overlapped particle images up to 50?% under noisy conditions, while keeping the root mean square peak location error under 0.07?pixels. The algorithm is also robust to the size changes up to a size ratio of 5. The tracking method developed from a classic computer vision matching algorithm is capable of capturing a velocity gradient up to 0.3 while maintaining the error under 0.2?pixels. Sensitivity tests were performed to describe the optimum conditions for the technique in terms of particle image density, particle image sizes and velocity gradients, also its sensitivity to errors of the PIV results that guide the tracking process. The comparison with other existing tracking techniques demonstrates that this technique is able to resolve more vectors out of a dense particle image field.     

Intensity Capping: a simple method to improve cross-correlation PIV results

A common source of error in particle image velocimetry (PIV) is the presence of bright spots within the images. These bright spots are characterized by grayscale intensities much greater than the mean intensity of the image and are typically generated by intense scattering from seed particles. The displacement of bright spots can dominate the cross-correlation calculation within an interrogation window, and may thereby bias the resulting velocity vector. An efficient and easy-to-implement image-enhancement procedure is described to improve PIV results when bright spots are present. The procedure, called Intensity Capping, imposes a user-specified upper limit to the grayscale intensity of the images. The displacement calculation then better represents the displacement of all particles in an interrogation window and the bias due to bright spots is reduced. Four PIV codes and a large set of experimental and simulated images were used to evaluate the performance of Intensity Capping. The results indicate that Intensity Capping can significantly increase the number of valid vectors from experimental image pairs and reduce displacement error in the analysis of simulated images. A comparison with other PIV image-enhancement techniques shows that Intensity Capping offers competitive performance, low computational cost, ease of implementation, and minimal modification to the images.     

A correlation-based continuous window-shift technique to reduce the peak-locking effect in digital PIV image evaluation

In this paper the peak-locking phenomenon is investigated in the evaluation of digital PIV recordings by using a correlation-based interrogation algorithm with a discrete window shift and a correlation-based tracking algorithm. Statistical analyses indicate that nonuniformly distributed bias errors are the main cause of the peak-locking effect, and the amplitude variation of the random error is also an important source of the peak locking. Simulations and experimental examples demonstrate that very strong peak-locking effects exist for the correlation-based interrogation algorithm with discrete window shift in the cases of large particle images, small interrogation windows, and very small particle images. Very strong peak-locking effects are also observed for the correlation-based tracking algorithm when the particle images are overexposed, binarized, or very small. These strong peak-locking effects can be avoided without loss of evaluation accuracy by using a continuous window-shift technique in combination with the correlation-based interrogation algorithm. Received: 2 July 2001 / Accepted: 28 November 2001     

A digital mask technique for reducing the bias error of the correlation-based PIV interrogation algorithm

 In this paper the bias phenomenon in the evaluation of PIV recordings by using the correlation-based interrogation algorithm is discussed, and a digital mask technique, that can effectively reduce the bias error, is introduced. The correlation-based interrogation algorithm, when masked with a Gaussian window function, can achieve a higher evaluation accuracy not only for PIV recordings of flows with small velocity gradients, but also for that of flows with large gradients. Received: 14 October 1998/Accepted: 20 July 1999     

Robust evaluation of the dissimilarity between interrogation windows in image velocimetry

Image velocimetry techniques, which extract motion information by comparison of image regions, typically make use of cross-correlation to measure the degree of matching. In this work, a novel measure of the dissimilarity between interrogation windows is proposed which is based on a more robust estimator than cross-correlation. The method is validated on synthetic images and on two experimental data sets obtained from a periodically pulsed jet and a backward-facing step. The former is a basically laminar flow, whereas the latter is fully turbulent. Both of them are characterized by regions of high velocity gradients. The efficiency of the robust image velocimetry (RIV) is compared with a cross-correlation algorithm (PIV). The analysis of results shows that the RIV is less sensitive to the appearance and disappearance of particles, and to high velocity gradients and, in general, to noise, generating less spurious velocity vectors. As a consequence RIV resolves better the vorticity peaks at the center of the vortex rings generated by the pulsed jet, obtaining, for a given interrogation window size, a higher spatial resolution. Moreover, in the analysis of the flow field generated by the backward-facing step, the RIV performs better in the shear layer at the border of the recirculation region, leading to a more reliable estimation of Reynolds shear stress and horizontal velocity component.     

Simultaneous measurements of velocity and density in buoyancy-driven mixing

A variant of the particle image velocimetry (PIV) technique is described for measuring velocity and density simultaneously in a turbulent Rayleigh-Taylor mixing layer. The velocity field is computed by the usual PIV technique of cross-correlating two consecutive images, and deducing particle displacements from correlation peaks of intensity fields. Different concentrations of seed particles are used in the two streams of different temperature (density) fluids, and a local measure of the density is obtained by spatially averaging over an interrogation window. Good agreement is reported between the first- and second-order statistics for density obtained from this technique and from a thermocouple. Velocity-density correlations computed by cross-correlating individual time series are presented. The errors in the density measurements are quantified and analyzed, and the issue of spatial resolution is also discussed. Our purpose for this paper is to introduce the PIV-S method and validate its accuracy against corresponding thermocouple measurements.     

Particle relaxation and its influence on the particle image velocimetry cross-correlation function

A study of some aspects of tracer particle responses to step changes in fluid velocity is presented. The effect of size distribution within a seed material on measured relaxation time is examined, with polydisperse particles of the same median diameter shown to possess a significantly higher relaxation time than their monodisperse counterparts when measured via a particle image velocimetry algorithm. The influence of a shock wave–induced velocity gradient within a PIV interrogation window on the correlation function is also examined using the noiseless cross-correlation function of Soria (Turbulence and coherent structures in fluids, plasmas and nonlinear media. World Scientific, Singapore, 2006). The presence of a shock is shown to introduce an artificial fluctuation into the measurement of velocity. This fluctuation is a function of the shock position, shock strength, spatial ratio and particle distribution. When the shock is located at the middle of the window, the magnitude of the fluctuation increases monotonically with increasing spatial ratio, increases asymptotically with shock strength, and decreases for increasing particle polydispersity. When the shock is located at the left-hand edge of the window, the magnitude of the artificial fluctuation is highest for intermediate spatial ratios, going to zero at infinitely high and low values. In this instance, particle polydispersity acts to increase the magnitude of fluctuations in measured velocity. In both cases, particle polydispersity serves to broaden the PDF of measured velocity. For the cases presented herein, with a shock located within the interrogation window, the root mean square of the artificial velocity fluctuations reaches values in excess of 30% of the freestream velocity.     

Limits and accuracy of the Stereo-LFC PIV technique and its application to flows of industrial interest

Particle image velocimetry with local field correction (LFC PIV) has been tested in the past to obtain two components of velocity in a two dimensional domain (2D2C). When compared to conventional correlation based algorithms, this advanced technique has shown improvements in three important aspects: robustness, resolution and ability to cope with large displacements gradients. A further step in the development of PIV algorithms consists in the combination of LFC with the stereo technique, which is able to obtain three components of velocity in a plane (2D3C PIV). In this work this combination is implemented and its performance is evaluated carrying out the following two different tasks:

A new particle tracking algorithm based on deterministic annealing and alternative distance measures

We describe a new particle tracking algorithm for the interrogation of double frame single exposure data, which is obtained with particle image velocimetry. The new procedure is based on an algorithm which has recently been proposed by Gold et al. (Gold et al., 1998) for solving point matching problems in statistical pattern recognition. For a given interrogation window, the algorithm simultaneously extracts: (i) the correct correspondences between particles in both frames and (ii) an estimate of the local flow-field parameters. Contrary to previous methods, the algorithm determines not only the local velocity, but other local components of the flow field, for example rotation and shear. This makes the new interrogation method superior to standard methods in particular in regions with high velocity gradients (e.g. vortices or shear flows). We perform benchmarks with three standard particle image velocimetry (PIV) and particle tracking velocimetry (PTV) methods: cross-correlation, nearest neighbour search, and image relaxation. We show that the new algorithm requires less particles per interrogation window than cross-correlation and allows for much higher particle densities than the other PTV methods. Consequently, one may obtain the velocity field at high spatial resolution even in regions of very fast flows. Finally, we find that the new algorithm is more robust against out-of-plane noise than previously proposed methods. Received: 1 March 1999 / Accepted: 29 July 1999     

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.     

三曝光三相关PIV技术研究

对不等间隔三次曝光单幅记录的PIV粒子图像,本文提出一和中三相关位移诊断方法,三相关函数有两个次大峰,且不对称地人布在最大峰两侧,本方法可同时获得位移值及判别位移方向,可有效解决位移方向二义性问题,对线性涡流模拟粒子图像,应用本方法进行粒子位移诊断,结果证实了本方法对含涡复杂流动速度方向判别的有效性。     

Reynolds stress estimation up to single-pixel resolution using PIV-measurements

This article derives a theory for estimating Reynolds normal and shear stresses from PIV images with single-pixel resolution. The main idea is the analysis of the correlation function to identify the probability density function from which the Reynolds stresses can be derived in a 2-D regime. The work establishes a theoretical framework including the influence of the particle image diameter and the velocity gradients on the shape of the correlation function. Synthetic data sets are used for the validation of the proposed method. The application of the evaluation method on two experimental data sets shows that high resolution and accuracy are also obtained with experimental data. The approach is very general and can also be applied to correlation peaks that are obtained from sum-of-correlation PIV evaluations.     

Limitation and improvement of PIV

The deformation of particle image patterns by strong velocity gradients and out-of-pattern motions is a major source of error for the PIV (Particle Image Velocimetry) technique. This deformation is investigated and its effect on conventional PIV techniques is quantified for 2D flows. Simulations and comparisons with independent experiments verify the results.