Uncertainty Quantification in Particle Image Velocimetry

Sayantan Bhattacharya, Purdue University

Abstract

Particle Image Velocimetry (PIV) is a non-invasive measurement technique which resolves the flow velocity by taking instantaneous snapshots of tracer particle motion in the flow and uses digital image cross-correlation to estimate the particle shift up to subpixel accuracy. The measurement chain incorporates numerous sets of parameters, such as the particle displacements, the particle image size, the flow shear rate, the out-of-plane motion for planar PIV and image noise to name a few, and these parameters are interrelated and influence the final velocity estimate in a complicated way. In the last few decades, PIV has become widely popular by virtue of developments in both the hardware capabilities and correlation algorithms, especially with the scope of 3-component (3C) and 3-dimensional (3D) velocity measurements using stereo-PIV and tomographic-PIV techniques, respectively. The velocity field measurement not only leads to other quantities of interest such as Pressure, Reynold stresses, vorticity or even diffusion coefficient, but also provides a reference field for validating numerical simulations of complex flows. However, such a comparison with CFD or applicability of the measurement to industrial design requires one to quantify the uncertainty in the PIV estimated velocity field. Even though the PIV community had a strong impetus in minimizing the measurement error over the years, the problem of uncertainty estimation in local instantaneous PIV velocity vectors have been rather unnoticed. A typical norm had been to assign an uncertainty of 0.1 pixels for the whole field irrespective of local flow features and any variation in measurement noise. The first article on this subject was published in 2012 and since then there has been a concentrated effort to address this gap. The current dissertation is motivated by such a requirement and aims to compare the existing 2D PIV uncertainty methods, propose a new method to directly estimate the planar PIV uncertainty from the correlation plane and subsequently propose the first comprehensive methods to quantify the measurement uncertainty in stereo-PIV and 3D Particle Tracking Velocimetry (PTV) measurements. The uncertainty quantification in a PIV measurement is, however, non-trivial due to the presence of multitude of error sources and their non-linear coupling through the measurement chain transfer function. In addition, the advanced algorithms apply iterative correction process to minimize the residual which increases the complexity of the process and hence, a simple datareduction equation for uncertainty propagation does not exist. Furthermore, the calibration or a reconstruction process in a stereo or volumetric measurement makes the uncertainty estimation more challenging. Thus, current uncertainty quantification methods develop a-posterior models utilizing the evaluated displacement information and combine it with either image information, correlation plane information or even calibration “disparity map” information to find the desired uncertainties in the velocity estimates. The first chapter evaluates the existing 2D PIV uncertainty models prior to the ones developed in this dissertation work. The performance of four different methods (Peak to Peak Ratio (PPR), Mutual Information (MI), Image Matching (IM) and Correlation Statistics (CS)) were evaluated for two experimental cases of a cylinder wake and a free jet flow and the prescribed uncertainty was compared against the expected uncertainty or RMS error. The ground truth was established by a high magnification PIV measurement, which in turn was validated against a high frequency LDV measurement. The results indicated that the PPR and MI methods were less sensitive to the RMS error while the IM and CS methods, although had a better response to the variation in RMS error, underpredicted the uncertainty level. This work highlighted the need for further development in the field of planar PIV uncertainty estimation. The second chapter proposes a novel method to estimate the 2D PIV uncertainty directly from the cross-correlation plane. PIV being a statistical estimate, the evaluated displacement from the cross-correlation peak represents the most probable displacement for the interrogation window. However, individual particles contributing to the correlation peak do not shift by the exact same amount, due to measurement noise and thus the PDF of displacement can be obtained by extracting the phase information from the correlation plane. This idea is used to estimate the standard deviation of the PDF or the standard uncertainty. This method, namely Moment of Correlation (MC), showed a strong sensitivity to the primary error sources in PIV, as well as an improved uncertainty coverage, especially for window resolution of 64 pixels and above, for a wide range of flow fields. However, for synthetic images with very low image noise and for smaller interrogation windows the predicted uncertainty incurred a bias, of the order of 0.02 pixels, due to the resolution limit in estimating the PDF standard deviation.

Degree

Ph.D.

Advisors

Slabaugh, Purdue University.

Subject Area

Computer science|Fluid mechanics|Mechanics|Optics

COinS