A hybrid digital particle tracking velocimetry technique
|
|
- Julian Gallagher
- 5 years ago
- Views:
Transcription
1 Experiments in Fluids (1997) Springer-Verlag 1997 A hybrid digital particle tracking velocimetry technique E. A. Cowen, S. G. Monismith 199 Abstract A novel approach to digital particle tracking velocimetry (DPTV) based on cross-correlation digital particle image velocimetry (DPIV) is presented that eliminates the need to interpolate the randomly located velocity vectors (typical of tracking techniques) and results in significantly improved resolution and accuracy. In particular, this approach allows for the direct measurement of mean squared fluctuating gradients, and thus several important components of the turbulent dissipation. The effect of various parameters (seeding density, particle diameter, dynamic range, out-of-plane motion, and gradient strength) on accuracy for both DPTV and DPIV are investigated using a Monte Carlo simulation and optimal values are reported. Validation results are presented from the comparison of measurements by the DPTV technique in a turbulent flat plate boundary layer to laser Doppler anemometer (LDA) measurements in the same flow as well as direct numerical simulation (DNS) data. The DPIV analysis of the images used for the DPTV validation is included for comparison. 1 Introduction In the study of turbulent flows, particle image velocimetry (PIV) is a valuable tool capable of measuring the significant spatial structure of turbulence. Liu et al. (1991) have successfully demonstrated the use of PIV to probe the structure of a turbulent channel flow. Previous experimental studies of Received: 9 August 199/Accepted: 31 May 199 E. A. Cowen, S. G. Monismith Environmental Fluid Mechanics Laboratory Stanford University, Stanford, CA , USA Correspondence to: E. A. Cowen The authors gratefully acknowledge the support of this work by the Fluid Dynamics Program, Office of Naval Research (Scientific Director: Dr. E. Rood, grant N ) and the High Performance Computing and Communication Program, National Science Foundation (Scientific Director: Dr. A. Thaler, grant ASC 93181). The authors wish to thank Jonathan Harris, whose PIV work laid the foundation from which our technique grew and who never gave up on PTV, Bob Street, Jeff Koseff and Chris Rehmann for their steady supply of ideas and assistance, and Bob Brown for his mechanical wizardry. We are particularly indebted to John Crimaldi for making the detailed LDA measurements reported in Sect. 5. turbulent flows have generally relied on point measurement techniques laser Doppler anemometry (LDA) and hot wire velocimetry. While providing high temporal resolution, these methods have extremely limited spatial resolution. Until recently, high resolution accurate PIV techniques, such as that of Liu et al., who report spatial resolution of 0.3 mm 1.0 mm, were based on interrogating photographic images. This process is time consuming, as the photographs must first be developed and then interrogated to determine if the experiment has been successful. Willert and Gharib (1991) demonstrated that it is possible to use a CCD based camera to acquire digital images eliminating the time consuming development and interrogation of photographs. Their technique (DPIV), however, must be considered low resolution as they report a spatial resolution of only 5. mm 5. mm (we note that this is more a function of their chosen optical configuration and not their technique). Westerweel et al. (199) demonstrated that DPIV can operate at resolutions on the order of conventional photographic PIV (and for that matter DNS) with comparable accuracy. They interrogated digital images at a spatial resolution of 1.3 mm 1.3 mm which is comparable to that achieved by Liu et al. However, correlation based analysis does not extract the maximum amount of velocity information from an image. Recently Keane et al. (1995)1 demonstrated that tracking a significant number of particles in a photographic image is feasible. Essentially, they followed Guezennec and Kiritsis (1990) algorithm and used an auto-correlation analysis estimate of the velocity field to guide the particle tracking algorithm. Their results suggest an improvement in resolution by a factor of.5 (to a dimensional grid spacing of about 0.1 mm) over their conventional PIV approach. While they recognized the power of such an algorithm to improve spatial resolution, it appears that the hybrid approach can offer substantial gains in accuracy over PIV as well. As we will demonstrate, PTV algorithms are inherently more accurate than correlation-based PIV algorithms since they are relatively unaffected by the presence of displacement gradients. We believe this behavior is typically overlooked because the common practice of interpolating the randomly located tracked velocity vectors onto a regular grid introduces noise that obviates the noise reduction obtained with PTV relative to PIV (Agüí and Jiménez 1987; Hesselink 1988). 1 The authors wish to thank R. J. Adrian for making us aware of his group s efforts to use correlation analysis to guide a particle tracking algorithm and for a pre-print of their work.
2 00 In this paper we present a DPTV technique capable of both high spatial resolution and high accuracy. Our technique is able to capture the significant spatial structure of turbulence, including both mean velocity statistics and instantaneous spatial gradients. Its essential feature is the use of the randomly located tracked velocity vectors in place to determine all turbulent statistics. This results in resolution and accuracy comparable to that found by Keane et al. (1995) under conventional photographic conditions but through fully digital means. Similar to Keane et al. and Guezennec and Kiritsis (1990), it relies on a correlation analysis of images to estimate the displacement of particles as a function of their position. It uses this estimate to track individual particle displacements from which the instantaneous velocity vectors are determined. A major advance in our approach relative to other DPTV approaches is the elimination of the need to interpolate the random vector fields. Since the approach is digital, we are able to collect greater than a thousand velocity field realizations with relative ease and determine the turbulent statistics by binning the random velocity vectors into small measurement volumes, thus avoiding the need to interpolate the data. We have also developed a technique for directly determining the instantaneous gradients from the randomly located data. In order to optimize our DPTV algorithm we used a Monte Carlo simulation to investigate the effects of various parameters on its accuracy. We present results and conclusions drawn from these simulations. The technique was validated against both LDA and DNS measurements of a flat plate boundary layer and the measurements of mean quantities, turbulence intensities, stresses, as well as dissipation are presented. We include both the DPIV and the DPTV results for this validation experiment allowing for a direct comparison of the two algorithms. PIV/PTV hybrid PTV has been employed in many ways but essentially there are two distinct types: multiply-exposed single images and singly-exposed multiple images. Multiply-exposed single image techniques (e.g. Agüí and Jiménez 1987) rely on relatively sparse seeding of the flow to avoid overlapping images. Because we were interested in a high-resolution technique, we chose to base our technique on singly-exposed multiple images (e.g. Hassan et al. 199). The challenge in developing a singlyexposed multiple image PTV technique is tracking a given particle through sequential images which contain a relatively high density of particles. Previous PTV implementations have relied on large search windows in the second image of an image pair in order to track particles in the pair (Hassan et al. 199). This necessitates low seeding densities to avoid pairing ambiguity. However, the size of the search window can be significantly reduced if an accurate estimate can be made of the particle s location in the second sub-window (e.g. Guezennec and Kiritsis 1990; Keane et al. 1995). As pointed out by a reviewer the method of binning is really zero-order interpolation. However, due to the small length scales associated with the binning process, typically one sixth the length scales associated with a 3 3 pixel sub-window, the interpolative noise incurred by this zero-order interpolation is minimal. Cross-correlation based PIV is probably the most accurate form of correlation based PIV (Keane and Adrian 199). By employing cross-correlation based PIV to two images, an accurate function to estimate a particle s position in the second image can be determined. Using this function the search window for the particle in the second image can be minimized, allowing PTV to perform at very high seeding densities (we have demonstrated the algorithm successfully at seeding densities greater than 0 particles per 3 3 sub-window in simulation). There are two significant advantages in using PIV as a guide to track all of the particles in an image pair: increased accuracy and improved resolution. Traditionally the main arguments against PTV have been that the data is randomly located and low seeding densities are required to avoid ambiguity in determining particle pairs (Hesselink 1988; Adrian 1991; Westerweel 1993a). Using the results of PIV as an estimate to a particle s displacement allows a PTV algorithm to search a small region for the particle s pair, markedly increasing the seeding densities at which PTV can successfully operate. Researchers have previously interpolated the randomly located data onto a grid in order to calculate turbulent statistics (e.g. Agüí and Jiménez 1987; Keane et al. 1995) but this results in interpolative errors (Agüı and Jiménez 1987; Hesselink 1988). We developed a DPTV technique that was originally conceived as the merger of Willert and Gharib s (1991) DPIV technique with the PTV method of Hassan et al. (199). Our method can be summarized as follows: 1. Precondition the images (i.e. remove any mean background noise and non-uniformity).. Make a coarse DPIV analysis pass of two images (nonoverlapping 3 3 pixel sub-windows). 3. Remove stray vectors from this vector set.. Develop a displacement estimating function from which the displacement at any pixel location in the first image can be estimated. 5. Threshold the images and search the result for patterns that meet a user definition of a particle.. Search for particle pairs. 7. Remove stray vectors from this vectors set. 8. Calculate the turbulent statistics. The specifics of each step are presented in the sections that follow..1 Preconditioning the images An advantage of digital images is the ability to easily manipulate them to remove the effects of non-ideal aspects of the imaging system. Preconditioning of the images for both PIV and PTV algorithms is generally beneficial (Westerweel 1993a; Hesselink 1988). Our pre-processing consists of analyzing the entire set of images to determine the minimum value at each pixel. These minimum values are assembled into a minimum image which is then subtracted from all of the original images. This pre-processing removes any constant noise source as well as the effects of non-uniform illumination (for example from the Gaussian intensity profile of a light sheet
3 expanded with a cylindrical lens) while enhancing the signalto-noise ratio of the images. Non-uniform background levels will cause problems in the thresholding of the image to find the particles; it is particularly important to correct for them.. The DPIV pass Our DPIV algorithm is a modification of the algorithm outlined in Willert and Gharib (1991). We use non-overlapping 3 3 pixel sub-windows, however, the location of the sub-window in the second image is determined dynamically. Through an iterative process the DPIV algorithm determines the amount to shift the second image sub-window. At each iteration the mean displacement between the sub-windows is determined. When the nearest integer pixel displacement is non-zero, the second image sub-window is shifted by this displacement in the appropriate direction. The process is repeated until the nearest integer pixel displacement is zero. If the displacement does not converge in three iterations, the sub-window is flagged as invalid. In this way our DPIV algorithm accomplishes two tasks: areas with insufficient seeding density are flagged as invalid, and the sub-window in the second image is displaced to ensure that a majority of particles in the first sub-window are imaged in the second sub-window. This latter step reduces a principal source of mean bias associated with in-plane loss of correlation (discussed further in Sect. 3.1) by correlating mean displacements that are limited to 0.5 pixels. We note that our iterative technique to locate the second image sub-window does not eliminate the possibility of particles with large displacements, relative to the mean (e.g. in strong shear), not residing in this sub-window; it merely reduces it. Keane and Adrian (199) use a different approach to ensure all particles in the first sub-window are in the second sub-window. They allow the second sub-window to be larger than the first. This approach requires the second sub-window to be pixels to accommodate the FFT. We prefer our method because it has slightly lower computational costs and the important ability to detect sub-windows that are not able to produce a valid velocity vector. A small gradient bias is incurred (a maximum of about 0.03 pixels) which will be removed when the individual particles are tracked..3 Removal of stray vectors The PIV and the PTV algorithm will, on occasion, produce erroneous results. We developed a routine to remove stray vectors based on the statistics of the displacements being analyzed. The probability density function of turbulent fluctuations can be approximated as Gaussian (Lesieur 1987). While this assumption is tentative near boundaries of the flow, where significant skewness can occur, it turns out to be adequate if used only to set maximum possible deviations from mean values. Our DPTV system is designed to measure turbulent statistics and hence we generally work with hundreds if not thousands of velocity field measurements. We invoke the Gaussian assumption to predict the maximum expected deviation from the mean quantities given the sample size. For the DPIV analysis of the data this entails calculating the ensemble statistics at each 3 3 sub-window location. Based Count a Count b Displacement (pixels) Displacement (pixels) Fig. 1a, b. Velocity density functions. a DPIV pass; b DPTV pass, pre-filter, post-filter on a sample s size and variance, the maximum expected fluctuation is determined and fluctuations beyond this are removed. The statistics are then recalculated and the thresholding of the fluctuations is repeated. This whole process is repeated until it has converged (i.e. no velocity fluctuations lie outside the determined threshold values). The traditional estimator of the mean, n i /N is not robust to skewness in the outliers (Snedecor and Cochran 1989), therefore we use the median value to estimate the mean in the first iteration. For all subsequent iterations the traditional estimator of the mean is used. Typical results for a DPIV data set are shown in Fig. 1a. The filter rejected about % of the raw DPIV data which has relatively high uncertainty because the seeding density was below optimal, about particles per 3 3 pixel sub-window. This vector filtering process is first done on the entire DPIV data set which is then used to construct the displacement estimator functions for each image pair. The DPTV algorithm is applied and this data set is filtered in the same manner. We divide the images into small two-dimensional measurement areas and bin the randomly located vectors into these areas. The areas are sized to allow the assumption of homogeneous statistics within the measurement area. This assumption is easily checked by looking at the variation of statistics in the neighboring measurement volumes. Enough images are collected to allow each binned measurement area to contain thousands of velocity vectors. Typical results for a DPTV data set are shown in Fig. 1b. The filtered DPIV data is clearly providing an excellent estimate of the particle displacements 01
4 0 as only about 0.5% of the data is rejected by the filter. Note Fig. 1a, b were calculated from the same data set (the validation data discussed in Sect. 5).. The displacement estimating function Using the DPIV data we develop a displacement estimating function based on the Hardy multiquadratic scheme (Hassan et al. 199; Hardy 1971), a global interpolation scheme given by V i (x, y) (u i, v i ) (1) u (x, y ) a 1 d () i i i j ij j 1 v (x, y ) b 1 d (3) i i i j ij j 1 where d ij (x i x j ) (y i y j ) () V i the velocity vector at coordinates x i, y i with components u i, v i N the number of vectors to be used in the interpolation x i,j, y i,j the location of the vectors a j, b j constants to be determined Using the determined u i and v i at x i and y i the a j and b j are solved for. Given a j and b j and the location of the original data the velocity at any point (u i, v i ) can be determined. We choose this scheme over other schemes because it exactly reproduces the original data and because using non-overlapping subwindows leaves N relatively small, so that we can afford to use the order N3 calculation implied by the inversion of Eqs. (1) (). If computational expense becomes a concern it may prove necessary to switch to an order N process such as adaptive Gaussian window (see Spedding and Rignot 1993)..5 The PTV algorithm Our PTV method can be summarized as follows: 1. Threshold the two images at a user-defined value above the local mean intensity of the image. This local mean intensity is determined by calculating the mode of each image along a column (or row) perpendicular to the direction of illumination. The modes (a function of the direction of the optical light path) are smoothed using a uniformly weighted six nearest neighbor kernel.. Store the thresholded image as a binary image. 3. Identify and label the particles in each image. This is achieved by scanning the binary image line by line and clustering the binary l s into particles based on their proximity to other binary 1 s. A pixel is considered part of a particle if it is adjacent horizontally or vertically to another pixel with value 1. Each pixel is given a number corresponding to the particle of which it is a part.. The location of each particle is determined with a Gaussian sub-pixel fit estimator (this is expanded upon in the Sect. 3.). 5. Particles in the first image are paired with particles in the second image. The position of a particle in the second image is estimated using the above global interpolation scheme and a9 9 sub-window is centered on each position. A local correlation is performed and the estimate of the second image particle location is refined based on this determined displacement. Finally a 3 3 pixel sub-window (found to be adequately large for estimates based on the refined DPIV results) is centered at the estimated position in the second image. If, and only if, exactly one image two particle is located in this 3 3 pixel sub-window it is considered the pairing particle.. Determining statistics, including instantaneous gradients The stray vector filter described in Sect..3 is applied to the entire data set of tracked particles. The statistics on the valid vectors are calculated at each measurement bin. A differentiating capability of full-field methods is the ability to measure instantaneous velocity gradients, and thus mean squared gradients, yet these are never reported. The reason for this is likely the increased noise level that results from differentiating velocity data with uncertainty in it. If the uncertainty (or noise) in the data could be filtered out, or at least reduced to low levels, it would be trivial to calculate accurate mean squared gradients for DPIV or DPTV type data. We have developed a dynamic process for calculating the spatial derivatives of randomly located PTV data. Because we are interested in turbulent phenomenon the smallest length scale that needs to be resolved is the Kolmogorov scale, η (ν3/ε)1/. In fact, this scale generally does not need to be resolved in order to fully characterize the flow. Using the proposed universal spectrum of Pao (195) to integrate the dissipation spectrum, (kη)φ(kη), one can show that 99% of the dissipation takes place for k 5.5 η 1 (as shown in Fig. ). Our algorithm determines derivatives by identifying each image one particle, say at (x 1, y 1 ), and constructing a small measurement area in the following image centered at (x 1 cη, y 1 ) and (x 1, y 1 cη) where c is a constant set to 5.5 (the sampling frequency required to measure 99% of the dissipation). A pictorial representation of this process is shown in Fig. 3. If a single particle is found in either of these two small measurement areas, with typical side length of about 0.05η, the appropriate derivative is calculated. The process is dynamic because η is estimated directly from the data based on the estimation for the dissipation ε c ν[ ( u / x) ( v / y) ( u / y) ( v / x)] (5) where c is an O(1) constant that takes into account the missing mean squared gradients. The first iteration requires that the user start the algorithm with an estimate of η which is taken as a constant throughout the domain for this iteration. Fortunately, since η depends weakly on ε (η ε 1/), this process converges relatively quickly to final estimates of the gradients (and η). As with the PTV vectors, the gradient quantities themselves are arrayed irregularly in space, and at lower densities than the velocity vectors. The in-plane mean squared gradients are determined in the same manner as the PTV statistics, namely the randomly located gradients are binned into small measurement areas and the mean squares are determined. The vorticity and divergence of the velocity field can be determined
5 (k η) Φ(k η) kη Fig.. Fraction of total dissipation as a function of resolved wavenumber dissipation, fraction of total dissipation Fraction of total dissipation Gradient (pixels/pixel) Fig.. Monte Carlo simulations of four sub-pixel estimators rms error: DPIV three point center-of-mass, DPTV three point center-of-mass, DPIV five point center-of-mass, DPTV five point center-of-mass, DPIV Gaussian, DPTV Gaussian, DPIV parabolic, DPTV parabolic 03 Fig. 3. Schematic drawing of the dynamic derivative algorithm in a similar manner by calculating both the / x and / y gradients when vectors randomly occur in the x and y components of the stencil shown in Fig. 3 simultaneously. 3 The Monte Carlo simulations To investigate the important parameters of the PIV and PTV processes, we developed a numerical model to simulate them. This approach has been taken by many researchers, including Westerweel (1993a), Keane et al. (1995), Keane and Adrian (199) and Guezennec and Kiritsis (1990). Our model generates images with randomly located, uniformly distributed particles that have a Gaussian intensity profile (which we acknowledge to be somewhat ideal). Gaussian noise can be superimposed on the entire image. The parameters in the model include background noise intensity, particle size, dynamic range, particle density, and out-of-plane motions. The CCD camera with which we typically work (see Sect. ) has a background noise standard deviation of about counts, where a count is defined as one least significant bit, so we held this value for all of our tests. We fixed the correlation subwindow size at a relatively small (Prasad et al. 199) 3 3 pixels used by Willert and Gharib (1991). This sub-window size was chosen because it allows significant spatial resolution on a CCD camera without significant loss of accuracy (Westerweel 1993b). Westerweel argues that it is the bandwidth of the spectral density that we are interested in since we are determining only the position of the particles (given by their low wave-number components) and not their exact shape (which would require the sampling of the high wave number range as well). His analysis shows that for particles in the μm range sampled in a 1 mm sub-window that sub-window sample rates of 3 3 to pixels are sufficient to obtain high levels of accuracy. Since we do not require that our DPIV analysis of the data have particularly high accuracy ( 0.5 pixels as an upper bound is sufficient), 3 3 pixels is more than adequate. The model enables us to apply an X and a Y displacement as well as displacement gradients X i / x j, i, j 1,. Using this model, we were able to systematically identify the effects of each of the parameters on the accuracy of our DPTV and DPIV processes. In order to insure that the statistics generated from the simulations converged, sub-windows were simulated for each set of parameter combinations (i.e. each data point shown in Figs. and 5a f). 3.1 Error type and structure Before reviewing the results of our simulations it is important to consider the type of error expected in PIV and PTV measurements. PIV relies on using the displacement correlation peak as a measure of displacement. However, when a gradient exists in the particle displacement field the second image sub-window can not contain all of the particles in the first image sub-window and the result is an in-plane loss-ofcorrelation (Adrian 1988; Keane and Adrian 199; Westerweel 1993a, b). This causes the correlation peak to be a biased estimate of the displacement. It is biased toward lower displacements since smaller displaced particle pairs have a higher probability of remaining in the second image subwindow. PTV does not suffer this type of mean bias since individual particles are tracked and a particle s image is not affected by displacement gradients. After Keane and Adrian, (199), we will refer to this as gradient biasing.
6 a Particle diameter (pixels) b Dynamic range (counts) Particle Diameter 18 = 1.8 pixels c Number of particles (per 3 3 sub-window) Vectors per 3 3 sub-window Vectors per 3 3 sub-window Vectors per 3 3 sub-window e Gradient (pixels/pixel) 0. f 0.0 Interpolation rms Dynamic range (counts) Vectors per 3 3 sub-window d Out-of-plane fraction (pixels/pixel) Vectors per 3 3 sub-window Fig. 5a f. Monte Carlo simulations a Particle size; b dynamic range; c seeding density d out-of-plane fraction: e gradient magnitude; f interpolation. DPIV mean error, DPTV mean error, DPIV rms error, DPTV rms error, DPIV valid vectors, DPTV valid vectors Both DPIV and DPTV rely on sub-pixel fit estimators to locate the centers of correlation peaks and particle images, respectively, to sub-pixel accuracy. Various subpixel fit estimators have been proposed (e.g. Willert and Gharib 1991; Adrian 1988), however, to date none of them provide an unbiased estimate (Westerweel 1993a). Following Westerweel, we will call these tracking biases. Both PIV and PTV suffer from tracking bias as both techniques rely on sub-pixel fit estimators. The goal becomes the minimization of this bias through optimal particle size and dynamic range. Finally, the imaging process is subject to random noise from various sources (e.g. light quantization, CCD dark current, particle blocking, etc.). This noise results in a random uncertainty in locating both the correlation peak and particle images. This uncertainty can be thought of as the rms fluctuation in the position returned by an algorithm for particles displaced exactly the same amount repeatedly. We will refer to this as rms error. The total error is the sum of all of these errors. In the results discussed below we plot the rms error and mean error, the sum of the gradient bias and tracking bias described above.
7 3. The optimal sub-pixel fit algorithm To improve the resolution and accuracy of our CCD based implementation we required an algorithm capable of tracking displacements to sub-pixel accuracy. Sub-pixel determination of displacements are calculated using intensity variation information over the correlation peak in PIV and over the particle image in PTV. Willert and Gharib (1991) report that a three-point exponential curve fit outperforms a centroid based fit while R.D. Keane concludes from his simulations that in the presence of gradients the center-of-mass method performs best (private communication to Prasad et al. 199). Westerweel (1993a) performed a detailed analytical and experimental study of three sub-pixel fit estimators: center-ofmass, parabolic, and Gaussian. He found that, in general, center-of-mass sub-pixel fit schemes are strongly biased towards integer displacements. His work showed that the Gaussian sub-pixel fit scheme performed best. We investigated the same sub-pixel fit schemes sensitivity to displacement gradients and our results are shown in Fig.. While Westerweel confines his investigation of the center-of-mass estimator to a three point model we also investigate the five point model, in particular to see if this may help to improve the results under highly strained conditions for PIV. The parameters used for our simulations were: particle diameter. pixels, dynamic range 00 counts, 1 particles per 3 3 sub-window. We find for both correlation peak-fitting (PIV) and particle image fitting (PTV) that the Gaussian sub-pixel fit performs the best. We note that we are primarily interested in relatively small particle image diameters (1.5 3 pixels), our Monte Carlo simulations have indicated that the optimal sub-pixel fit algorithm may be a function of particle image diameter. Not surprisingly, the three point center-of-mass estimator outperforms the five point estimator for particle image fitting (the particle diameter was. pixels) while the opposite holds for correlation peak fitting. The apparent reduced rms error at a gradient of 3% for the center-of-mass sub-pixel estimator results from the near-integer displacement when interrogating a 3%displacement gradient with a 3 pixel interrogation spacing. It was pointed out above that these schemes are strongly biased toward integer displacements and hence perform better when displacements are restricted to nearinteger. We recognize that the choice of Gaussian particle intensity profiles skews the analysis in favor of the Gaussian sub-pixel estimator, however, given Westerweel s (1993a) analytic work and the reality that particles imaging near the diffraction limit of an optical system have intensity profiles that are well approximated by a Gaussian we feel the results are useful. Given our work and that of Westerweel we choose to use the Gaussian sub-pixel estimator for both the correlation peak and particle image fits in our DPTV algorithm. For reference the Gaussian sub-pixel estimator is x i (ξ ξ 1 )ln(γ /γ 3 ) (ξ 3 ξ )ln(γ 1 /γ ) [(ξ ξ 1 )ln(γ /γ 3 ) (ξ 3 ξ )ln(γ 1 /γ )] where x i the sub-pixel estimate of position () ξ 1, ξ, ξ 3 the pixel locations of the left of peak, peak, and right of peak intensity γ 1, γ, γ 3 the intensity value at ξ 1, ξ, ξ 3 We note that it may be possible to improve the performance of correlation peak fitting under highly strained conditions (displacement gradients greater than 3%) through the use of asymetric, two dimensional peak fits (Westerweel et al. 1995). This type of fit may also give modest improvements to particle center peak fits, particularly for non-spherical particles. 3.3 Particle size To investigate the effect of particle size on error, we generated images with an average seeding density of twelve particles per 3 3 pixel sub-window, no out-of-plane motions, and a dynamic range of 00 counts. A linear displacement gradient of 3% (0.03 pixels/pixel) was applied to the data. The results of this simulation are shown in Fig. 5a. Note that the right axis is discontinuous for the Figs. 5a e. The mean error for DPTV is essentially zero. For DPIV it is also essentially constant but at a bias of about 0.03 pixels. This is principally due to the gradient bias effects on the PIV algorithm. The rms error for DPTV drops off initially with increasing particle diameter, reaching a minimum error at a particle diameter of about. pixels before increasing again for larger diameters. The DPIV rms error s functional dependence on particle diameter is similar to DPTV s, however, the error is roughly twice as large and reaches its minimum closer to a diameter of 3.0 pixels. The significant increase in resolution is clearly reflected in the number of valid vectors found per 3 3 pixel sub-window for each technique. As expected, DPIV produces very close to one vector per sub-window. However, DPTV resolves nearly an order of magnitude greater valid vector density with a particle diameter of close to pixels. It is not surprising that the valid vector density decays strongly with increasing diameter as particle image overlapping increases with increasing diameter. DPIV does not suffer this loss of valid vector density with increasing particle diameter nearly as severely since overlapping particle images still contribute to the correlation while they are discarded by the DPTV algorithm. Reconciling the desire for a high valid vector density with the need for low error, the optimal particle size for DPTV is taken to be between 1.8. pixels. Our results for DPIV indicate an optimal particle size in the broad range of.0 to nearly.0 pixels. This is larger than the slightly less than.0 pixels suggested by Prasad et al. (199) for the centerof-mass estimator and 1.0 pixels determined by Westerweel (1993a) for Gaussian and parabolic estimators. However, both Prasad et al. and Westerweel looked at linear displacements in the absence of shear. In contrast, our results suggest that the presence of a gradient increases the optimal particle size for the correlation analysis. This is not surprising as the correlation algorithm will be more robust for larger particles since there will be a higher probability of sub-window particle correlations in the presence of gradients for larger particles. 3. Dynamic range Fixing the particle diameter at. pixels we varied dynamic range, defined as the difference in mean maximum particle 05
8 0 intensity level and the mean background intensity level, to investigate its effect on error. From experience gathered working with several traditional 8-bit CCD cameras we find that typically we are able to produce dynamic ranges of about counts. As Fig. 5b shows, the accuracy of DPTV is strongly tied to the dynamic range while DPIV accuracy is relatively independent of dynamic range, in agreement with Westerweel (1993a). As expected, the DPTV mean error is essentially zero while the DPIV mean error again shows the negative bias due to the mean gradient. The DPTV rms error decays monotonically with increasing dynamic range, while for dynamic ranges above 50 counts the DPIV rms error is essentially constant at a level greater than twice the asymptotic limit of the DPTV rms error. The valid vector density for both DPIV and DPTV is essentially independent of dynamic range. This simulation suggests that the optimal dynamic range for DPTV is at minimum 00 counts with clear gains out to a dynamic range of 500. These results lead us to the conclusion that while 7-bit cameras are sufficient to minimize error for DPIV, true 10-bit or better cameras (potential dynamic range of 10) offer significant advantages over lower bit depth devices for DPTV. We note that we have been working with rather pedestrian particles, Pliolite VT-AC-L, which sell for around two dollars per pound. It is possible to more effectively use the dynamic range of a camera by imaging fluorescent particles (Willert 199) and bandpass filtering the scattered light. This may allow dynamic ranges in excess 00 to be obtained with an 8-bit camera, however, these particles are expensive (about $500 per pound) and our experimental facilities are large, O(10 m3), making a 10-bit camera the less expensive alternative. 3.5 Seeding density Maintaining the particle diameter of. pixels and dynamic range of 00, we varied the particle seeding density, defined as the mean number of particles per 3 3 sub-window. Figure 5c shows that, as expected, the DPTV mean error is essentially zero while the DPIV mean error shows the negative bias discussed above. The DPTV rms error increases nearly linearly with increasing seeding density but remains below the rms error for DPIV for densities up to 5 particles per sub-window. DPTV rms error is half that of DPIV for seeding densities below 1 particles per sub-window. The DPIV rms error decays monotonically with increasing seeding density and can be driven down to the level of optimal DPTV for very large seeding densities (greater than 0 particles per subwindow). We have assumed, however, moderate velocity gradients and zero out-of-plane motion which is tenuous in turbulent flows. It is shown in the next section that for even moderate out-of-plane motions DPIV rms error increases significantly relative to DPTV rms error. The DPTV valid vector density initially increases almost linearly with increasing seeding density but this growth rate decays, a result of particle image overlapping. The dashed lines are the results under the same conditions but with a particle diameter of 1.8 pixels. Clearly much larger valid vector densities can be achieved at relatively small rms error cost if necessary. Our simulation suggests that the optimal seeding density for DPTV is 5 0 particles per sub-window (balancing the desire for greater valid vector density with low error). Our DPIV simulations are in agreement with Keane and Adrian (199), suggesting that seeding density should be greater than 10 particles per 3 3 sub-window. For seeding densities greater than 10 particles per sub-window the DPIV valid vector density is essentially constant at Out-of-plane motions Single-camera particle imaging techniques are inherently two-dimensional measurement techniques while turbulence is fully three-dimensional and thus it is important to consider the robustness of DPTV and DPIV to out-of-plane motions. Maintaining the previously determined optimal parameter settings: particle diameter of. pixels, dynamic range of 00, and 1 particles per sub-window, we varied the out-of-plane fraction, where 0.0 corresponds to all particles in the first image remaining in the second image and 1.0 corresponds to no original particles remaining in the second image. Figure 5d demonstrates that DPTV error is relatively insensitive to out-of-plane fraction, while the valid vector density drops linearly as expected. In constrast, DPIV is extremely sensitive to out-of-plane fraction, with rms error doubling for fractions of about This agrees with Keane et al. (1991) who find that the valid detection probability is a strong function of out-of-plane fraction for a constant in-plane displacement. They recommend keeping the out-of-plane fraction below We find that the presence of an in-plane gradient requires a more severe restriction for DPIV (3 3 sub-window); out-of-plane fraction should be kept to below This is significant since all of the above simulations have been carried out at zero out-of-plane fraction, almost never a characteristic of turbulent flows of interest. Typically, turbulence intensities can exceed 10% near walls or in moderate shear conditions. Therefore, the simulations of rms error discussed above must be considered low estimates for DPIV yet reasonable estimates for DPTV. 3.7 Gradients Another significant advantage of DPTV over DPIV is its increased robustness in the presence of gradients. Figure 5e shows that while DPIV error (both rms and mean) increases sharply with gradient strength, DPTV error is nearly independent of gradient strength. Since the DPTV algorithm only requires that the DPIV estimate of position be good to 0.5 pixels, the DPTV algorithm is capable of tracking nearly constant valid vector density up to a displacement gradient of about 7%. This compares with a nearly linear increase in DPIV rms error of pixels/%. Thus it is not recommended that DPIV be used for displacement gradients in excess of about 3% which agrees with the parameters set by Keane and Adrian (199) while 7% is the limit for DPTV. 3.8 Interpolation error As discussed above, previous researchers have relied on interpolating randomly located PTV vectors onto a uniform grid in order to calculate turbulent statistics. We used the AGW interpolator with the optimal parameters suggested by Agüı
9 and Jime nez (1987) to interpolate the simulation data (this was the interpolator used by Keane et al. 1995). The results are shown for the dynamic range sensitivity test in Fig. 5f. Only the rms error is shown since the mean results are unaffected (indicating that AGW is an unbiased interpolator). Clearly the accuracy gains made by using the DPTV process are more than offset by the interpolation noise incurred using the AGW interpolator. We have tried tweaking the optimal parameters suggested by Agüí and Jime nez but with little improvement. Spedding and Rignot (1993) demonstrate that a spline thin shell (STS) interpolator may perform in a superior manner to AGW. We have only worked briefly with this interpolator and were unable to find significant improvement relative to AGW but STS deserves further exploration. 3.9 Summary of PTV findings and comparison to PIV The results of our Monte Carlo simulation are summarized in Table 1. Looking at Figs. 5a f and Table 1 it is apparent that the DPTV technique has performance characteristics that are superior to those of DPIV. In general DPTV eliminates mean error and reduces rms errors to one half to one third the associated rms DPIV error. Significantly, as we demonstrate, DPTV determines an order of magnitude more valid vectors from an image than does DPIV. It should also be noted that our PTV technique is much more robust in the face of out-of-plane motions and strong gradients than correlation-based PIV. It is comparably insensitive to out-of-plane fraction relative to PIV and can operate accurately at better than twice the gradient strength relative to PIV. Specifications of our PTV hardware We designed our PTV system using the above analysis. The centerpiece of the imaging system is the camera. In order to measure velocity fields resolved to order 100 μm using 3 3 pixel sub-windows for the correlation analysis in laboratory scale water flows, it is necessary to capture two images between 3 and 0 ms apart. Since our technique needs to be digital we must use a CCD type camera. These cameras typically run at framing rates between 5 and 0 Hz (interimage times of 17 0 ms). There are some specialty cameras that are capable of framing rates greater than 500 Hz but these are quite small format, not greater than 5 5 pixels. Faster cameras appear to be on the horizon. Rood (199) states that Princeton Scientific Instruments is working on a one-million frame per second CCD camera. Table 1. Optimal values for DPTV and DPIV parameters Parameter DPTV DPIV Sub-pixel fit estimator Gaussian Gaussian Particle diameter (pixels) Dynamic range (counts) Seeding density (per 3 3 sub-window) Out-of-plane fraction Displacement gradient 7% 3% When only two images are required, as in our case, there is a solution to the CCD camera s limited framing rate (Willert 199). A few CCD chips are built in a format known as fullframe transfer a chip which captures an image and transfers it to an on-chip storage area, typically a masked off portion of the CCD array. Since the transfer is accomplished on chip, it occurs quite rapidly. Transfer times on the order of 1 ms are possible. A camera based on a full-frame transfer chip can be used to collect two images, spaced on the order of 1 ms apart, by strobing the image area, shifting the image into on-chip storage, and strobing the image area again. To take full advantage of the velocity field measurement capabilities of PIV requires that the image area be as large as possible. Since full-frame transfer chips are available in formats as large as pixels, we chose to use a full-frame transfer CCD camera. The camera we are presently working with is a cooled 1- bit digital slow-scan camera manufactured by Patterson Electronics of Tustin, CA. The camera is a full-frame transfer camera with a transfer time of about 3.5 ms. It has the advantage that it can house essentially any full-frame transfer CCD chip. Because of economics we began with a Texas Instruments TC-17 chip which has an active area of pixels. This chip has a fill factor of 100% but is not truly a scientific-grade chip. Its electron well is not quite deep enough to support a full 1-bit dynamic range, but it is capable of dynamic ranges in excess of 11 bits. It has a pixel size of 7.8 μm 13. μm and thus has an aspect ratio of 1.7. This causes a problem in trying to maintain an ideal particle diameter since the particle will appear roughly twice the diameter in one direction as the other. Despite these limitations the chip performs quite well as will be demonstrated in Sect. 5. Cooling the chip to 0 C reduces most of the thermal noise associated with the chip (the rms noise is reduced from 3. counts at room temperature to 1.8 counts) and any non-linearities and spatial gradients in sensitivity. We illuminate the flow with two Continuum Minilite Nd : YAG pulsed lasers. The lasers output wavelengths are frequency doubled to 53 nm, have typical power of 10 mj/ pulse, and have typical pulse length of 5 ns. They are combined into a single optical path via a polarizing beam splitter and a half-wave retarder before passing through a cylindrical lens to form a 1 mm thick light sheet. The laser and camera triggers are controlled with National Instruments LabVIEW software running on a Macintosh IIfx computer. The digital images are collected and processed on an Intel i80 based processing card, manufactured by Alacron, and hosted by a Pentium based PC equipped with 3.5 gigabytes of hard drive space. The images are written to the host s hard drives and stored on DAT tapes for later analysis. 5 Validation experiment the flat-plate boundary layer To validate our DPTV technique we made measurements in a partially developed free-surface channel flow with both the DPTV technique and a two-component LDA. We also include the results from a DPIV analysis of the data set for comparison. This flow is a good approximation of a zero pressure gradient flat-plate boundary layer. The Reynolds number of the flow 07
10 08 based on momentum thickness (Re ) was 1300 allowing the data to be compared directly to the zero pressure gradient flat-plate boudnary layer DNS at Re 110 of Spalart (198). 5.1 Experimental facility and LDA The flow facility used was a constant head type recirculating flume fitted with a smooth bottom, described in detail in O Riordan et al. (1993). The facility is 10 m in length and the flow is driven through an inlet section by a m constant head. The inlet section is equipped with a baffle to minimize secondary flow structures developing at the inlet jet. The baffle is followed by three screens of decreasing mesh size to further homogenize the flow. The flow enters the main channel section through a.5 : 1 two-dimensional contraction. To trip the boundary layer, there is a 3 mm diameter cylindrical rod located just downstream of the contraction. The main channel is m long and is equipped with glass side walls over the center 3 m. The flume was operated at a depth of 5.5 cm which produced a free-stream velocity of 11. cm/s. All data was taken with the image area (or LDA measurement volume) centered.0 m downstream of the boundary layer trip. We defined a right-hand coordinate system such that x is taken positive in the streamwise direction and z is the vertical distance above the smooth bottom. The two component LDA system included two Dantec trackers and a Spectraphysics argon-ion laser operating at 1 W with output beam wavelength of 88 nm. The beams were electronically shifted to differentiate the two components. The optical configuration resulted in a measurement volume with in-plane dimensions of mm and an out-of-plane (transverse) length of 1 mm. Doppler frequencies determined by the trackers were sampled at 80 Hz to yield data records of between and points, depending on turbulence intensity. The LDA was used to measure the entire boundary layer profile from which the momentum thickness was determined to be 1.1 cm, giving a Reynolds number based on momentum thickness, Re, of This Re was intentionally chosen close to Spalart s (198) DNS at Re 110 to allow direct comparison. To validate the DPTV technique LDA measurements were made at the same location as the DPTV measurements. 5. DPTV calibration, data collection and processing The DPTV system was calibrated by imaging a very fine grid etched in Plexiglas that is placed in the light sheet. The grid is 10 cm 10 cm with a grid spacing of 1.00 cm. A single image was manually interrogated to determine to the nearest pixel the location of the grid elements. In this way the distance per pixel in the horizontal and vertical directions are determined to better than 0.5% (1 part in 00 pixels). The calibration returned a horizontal value of 5. μm per pixel and vertical value of 31.1 μm per pixel which gives a magnification of 0.51 and yields an image area of mm. The flow was seeded with Pliolite VT-AC-L, manufactured by Goodyear Tire & Rubber Company. Pliolite has a specific gravity of 1.03 and good reflectivity. We used a mechanical shaker and sieves to acquire particles with a diameter of between 5 μm and 75 μm. We expected that particles in this size range would follow the fluid motion which we confirm using Eq. (7), derived by Adrian (1991) for the difference between particle and fluid motion, in Sect v u ρ d ρ d r vr 3ν where v the measured particle velocity u the exact velocity per the velocity field ρ d /ρ the specific gravity of the particle d r the particle diameter vr the particle acceleration We will estimate the peak acceleration, vr, asu /Lwhere L is an appropriate length scale. In order to maintain a uniform distribution of particles across the boundary layer we seeded our entire reservoir of water, about 8 m3, with approximately 13 g of particles. Elgoboshi (199) reports that particle laden flows show evidence of two-way interaction at volume fractions as low as 10. Our volume fraction was about However, Elgoboshi s findings are for monodisperse particles while our seed is polydisperse, with a factor of 8 difference in particle volume over the diameter distribution. The number density function for particle diameter is probably skewed heavily toward the low diameter range (due to the sieving process) which would give an effective void fraction closer to Elgoboshi s cutoff. (Westerweel et al. (1995a) use a similar argument to justify a void fraction in the same range). To measure the boundary layer flow, 109 image pairs were collected. The camera was mounted in a portrait orientation (horizontal axis of CCD area in the vertical, z, direction). This was done to take advantage of the asymmetry of the pixel dimensions (Liu et al. (1991) also use an asymmetric interrogation window in their boundary layer measurements). To obtain optimal signal-to-noise ratio the maximum mean displacement in each direction should be equal to about 35 0% of the sub-window length, about 1 pixels for a 3 3 sub-window. In a boundary layer flow the dominant velocity is the mean free-stream velocity and thus to obtain maximal signal-to-noise ratio we aligned this velocity component with the long side of the pixels. The time between images, Δt, was set at ms. The area of each image above the flat plate was pixels and thus the area imaged was mm. The DPIV pass on the data was performed on a 99 8 sub-region with non-overlapping 3 3 pixel sub-windows yielding 3 measurement sites. Due to a relatively low seeding density the DPIV algorithm converged on average in only 80 sub-windows. Again due to the low seeding density the stray vector filter determined that only 75.9% of these vectors were valid (13 sub-windows). The images were thresholded such that a particle contained at least two pixels with intensity greater than 15 counts above the local mean intensity level. This process yielded a mean particle density of 3. and. particles per sub-window for the image one and image two data sets respectively (the limiting factor was the volume of water that needed to be seeded). The image two set contains more particles due to a slightly greater output intensity of the second laser. The median dynamic range for the image set was (7)
The Elimination of Correlation Errors in PIV Processing
9 th International Symposium on Applications of Laser Techniques to Fluid Mechanics, Lisbon, Portugal, July, 1998 The Elimination of Correlation Errors in PIV Processing Douglas P. Hart Massachusetts Institute
More informationTime-resolved PIV measurements with CAVILUX HF diode laser
Time-resolved PIV measurements with CAVILUX HF diode laser Author: Hannu Eloranta, Pixact Ltd 1 Introduction Particle Image Velocimetry (PIV) is a non-intrusive optical technique to measure instantaneous
More information9.9 Coherent Structure Detection in a Backward-Facing Step Flow
9.9 Coherent Structure Detection in a Backward-Facing Step Flow Contributed by: C. Schram, P. Rambaud, M. L. Riethmuller 9.9.1 Introduction An algorithm has been developed to automatically detect and characterize
More informationModule 3: Velocity Measurement Lecture 14: Analysis of PIV data. The Lecture Contains: Flow Visualization. Test Cell Flow Quality
The Lecture Contains: Flow Visualization Test Cell Flow Quality Influence of End-Plates Introduction To Data Analysis Principle of Operation of PIV Various Aspects of PIV Measurements Recording of the
More informationParticle Image Velocimetry for Fluid Dynamics Measurements
Particle Image Velocimetry for Fluid Dynamics Measurements Lyes KADEM, Ph.D; Eng kadem@encs.concordia.ca Laboratory for Cardiovascular Fluid Dynamics MIE Concordia University Presentation - A bit of history
More informationThe effect of a discrete window offset on the accuracy of cross-correlation analysis of digital PIV recordings
Experiments in Fluids 23 (1997) 20 28 Springer-Verlag 1997 The effect of a discrete window offset on the accuracy of cross-correlation analysis of digital PIV recordings J. Westerweel, D. Dabiri, M. Gharib
More informationFLOW VISUALISATION OF POLYMER MELT CONTRACTION FLOWS FOR VALIDATION OF NUMERICAL SIMULATIONS
FLOW VISUALISATION OF POLYMER MELT CONTRACTION FLOWS FOR VALIDATION OF NUMERICAL SIMULATIONS R Spares, T Gough, M T Martyn, P Olley and P D Coates IRC in Polymer Science & Technology, Mechanical & Medical
More informationDriven Cavity Example
BMAppendixI.qxd 11/14/12 6:55 PM Page I-1 I CFD Driven Cavity Example I.1 Problem One of the classic benchmarks in CFD is the driven cavity problem. Consider steady, incompressible, viscous flow in a square
More informationParticle Velocimetry Data from COMSOL Model of Micro-channels
Particle Velocimetry Data from COMSOL Model of Micro-channels P.Mahanti *,1, M.Keebaugh 1, N.Weiss 1, P.Jones 1, M.Hayes 1, T.Taylor 1 Arizona State University, Tempe, Arizona *Corresponding author: GWC
More informationLecture # 11: Particle image velocimetry
AerE 344 Lecture Notes Lecture # 11: Particle image velocimetry Dr. Hui Hu Dr. Rye M Waldman Department of Aerospace Engineering Iowa State University Ames, Iowa 50011, U.S.A Sources/ Further reading:
More informationDoppler Global Velocimetry: A Potential Velocity Measurement Method for General Aviation Applications
Doppler Global Velocimetry: A Potential Velocity Measurement Method for General Aviation Applications L. Scott Miller The Wichita State University Wichita, Kansas and James F. Meyers, and Jimmy W. Usry
More informationSYNTHETIC SCHLIEREN. Stuart B Dalziel, Graham O Hughes & Bruce R Sutherland. Keywords: schlieren, internal waves, image processing
8TH INTERNATIONAL SYMPOSIUM ON FLOW VISUALIZATION (998) SYNTHETIC SCHLIEREN Keywords: schlieren, internal waves, image processing Abstract This paper outlines novel techniques for producing qualitative
More informationHydrodynamic Instability and Particle Image Velocimetry
Hydrodynamic Instability and Particle Image Velocimetry Instabilities in lid-driven cavities First important investigations of hydrodynamic instabilities were published by v. Helmholtz (1868), Lord Rayleigh
More informationIMAGE PROCESSING IN PIV
ISTP-,, PRAGUE TH INTERNATIONAL SYMPOSIUM ON TRANSPORT PHENOMENA IMAGE PROCESSING IN PIV Milan Pěta, Jan Novotný Czech Technical University in Prague, Fakulty of Mechanical Engineering, Division of Fluid
More informationTurbulencja w mikrokanale i jej wpływ na proces emulsyfikacji
Polish Academy of Sciences Institute of Fundamental Technological Research Turbulencja w mikrokanale i jej wpływ na proces emulsyfikacji S. Błoński, P.Korczyk, T.A. Kowalewski PRESENTATION OUTLINE 0 Introduction
More informationRASNIK Image Processing with a Steepest Ascent Algorithm
ATLAS Internal Note MUON-No-092 RASNIK Image Processing with a Steepest Ascent Algorithm Kevan S. Hashemi and James R. Bensinger Brandeis University August 1995 Introduction The RASNIK alignment instrument
More informationMeasurement Techniques. Digital Particle Image Velocimetry
Measurement Techniques Digital Particle Image Velocimetry Heat and Mass Transfer Laboratory (LTCM) Sepideh Khodaparast Marco Milan Navid Borhani 1 Content m Introduction m Particle Image Velocimetry features
More informationspecular diffuse reflection.
Lesson 8 Light and Optics The Nature of Light Properties of Light: Reflection Refraction Interference Diffraction Polarization Dispersion and Prisms Total Internal Reflection Huygens s Principle The Nature
More informationParticle Image Velocimetry Part - 1
AerE 545X class notes #23 Particle Image Velocimetry Part - 1 Hui Hu Department of Aerospace Engineering, Iowa State University Ames, Iowa 50011, U.S.A Announcement Room 1058, Sweeney Hall for Lab#4 (LDV
More informationMeasurements of Three-Dimensional Velocity Fields Under Breaking Waves
The Journal of Undergraduate Research Volume 11 Journal of Undergraduate Research, Volume 11: 2013 Article 3 2013 Measurements of Three-Dimensional Velocity Fields Under Breaking Waves Matthew Auch South
More informationLaser speckle based background oriented schlieren measurements in a fire backlayering front
Laser speckle based background oriented schlieren measurements in a fire backlayering front Philipp Bühlmann 1*, Alexander H. Meier 1, Martin Ehrensperger 1, Thomas Rösgen 1 1: ETH Zürich, Institute of
More informationPulsating flow around a stationary cylinder: An experimental study
Proceedings of the 3rd IASME/WSEAS Int. Conf. on FLUID DYNAMICS & AERODYNAMICS, Corfu, Greece, August 2-22, 2 (pp24-244) Pulsating flow around a stationary cylinder: An experimental study A. DOUNI & D.
More informationRange Imaging Through Triangulation. Range Imaging Through Triangulation. Range Imaging Through Triangulation. Range Imaging Through Triangulation
Obviously, this is a very slow process and not suitable for dynamic scenes. To speed things up, we can use a laser that projects a vertical line of light onto the scene. This laser rotates around its vertical
More informationDispersion of rod-like particles in a turbulent free jet
Test case Dispersion of rod-like particles in a turbulent free jet 1. MOTIVATION: Turbulent particle dispersion is a fundamental issue in a number of industrial and environmental applications. An important
More informationCharacterization of Single Pixel Evaluation PIV for Precise Flow Measurement
Characterization of Single Pixel Evaluation PIV for Precise Flow Measurement Han-Sheng Chuang 1, Steve T. Wereley 1, Lichuan Gui 2 1: School of Mechanical Engineering, Purdue University, West Lafayette,
More informationFlow Visualization around Generic Bridge Shapes using Particle Image Velocimetry
Flow Visualization around Generic Bridge Shapes using Particle Image Velocimetry by Harold Bosch 1 and Kornel Kerenyi 2 ABSTRACT This paper examines the flow field around generic bridge shape models using
More informationAn Experimental and Computational Investigation of a 3D, l/h=5 Transonic Cavity Flow. Prof Kevin Knowles Dr Simon Ritchie Dr Nick Lawson
An Experimental and Computational Investigation of a 3D, l/h=5 Transonic Cavity Flow Prof Kevin Knowles Dr Simon Ritchie Dr Nick Lawson Overview Background Experimental Studies Computational Studies Results
More informationDEVELOPMENT OF ROBUST CORRELATION ALGORITHMS FOR IMAGE VELOCIMETRY USING ADVANCED FILTERING
DEVELOPMENT OF ROBUST CORRELATION ALGORITHMS FOR IMAGE VELOCIMETRY USING ADVANCED FILTERING by Adric Eckstein Thesis submitted to the faculty of the Virginia Polytechnic Institute and State University
More informationApplication of fluorescent particles for particle tracking velocimetry in wind tunnels
Application of fluorescent particles for particle tracking velocimetry in wind tunnels Tamara Guimarães 1,*, K. Todd Lowe 1 1: Dept. of Mechanical Engineering, Virginia Polytechnic Institute and State
More informationFLOW VISUALISATION AROUND A SOLID SPHERE ON A ROUGH BED UNDER REGULAR WAVES
FLOW VISUALISATION AROUND A SOLID SPHERE ON A ROUGH BED UNDER REGULAR WAVES H.P.V.Vithana 1, Richard Simons 2 and Martin Hyde 3 Flow visualization using Volumetric Three-component Velocimetry (V3V) was
More informationAir Assisted Atomization in Spiral Type Nozzles
ILASS Americas, 25 th Annual Conference on Liquid Atomization and Spray Systems, Pittsburgh, PA, May 2013 Air Assisted Atomization in Spiral Type Nozzles W. Kalata *, K. J. Brown, and R. J. Schick Spray
More informationWake Vortex Tangential Velocity Adaptive Spectral (TVAS) Algorithm for Pulsed Lidar Systems
Wake Vortex Tangential Velocity Adaptive Spectral (TVAS) Algorithm for Pulsed Lidar Systems Hadi Wassaf David Burnham Frank Wang Communication, Navigation, Surveillance (CNS) and Traffic Management Systems
More informationA three-dimensional velocimetry approach using a combination of tomographic reconstruction and triangulation for double-frame particle tracking
A three-dimensional velocimetry approach using a combination of tomographic reconstruction and triangulation for double-frame particle tracking Thomas Fuchs *, Rainer Hain, Christian J. Kähler Institute
More informationFlow structure and air entrainment mechanism in a turbulent stationary bore
Flow structure and air entrainment mechanism in a turbulent stationary bore Javier Rodríguez-Rodríguez, Alberto Aliseda and Juan C. Lasheras Department of Mechanical and Aerospace Engineering University
More informationPARTICLE IMAGE VELOCIMETRY (PIV) AND VOLUMETRIC VELOCIMETRY (V3V) SYSTEMS
PARTICLE IMAGE VELOCIMETRY (PIV) AND VOLUMETRIC VELOCIMETRY (V3V) SYSTEMS VERSATILE, UPGRADEABLE FLUID MECHANICS MEASUREMENT SOLUTIONS UNDERSTANDING, ACCELERATED FULL SPECTRUM OF GLOBAL VELOCITY SYSTEMS
More informationPIV Data Validation Software Package
NASA Contractor Report 2070 PIV Data Validation Software Package James L. Blackshire ViGYAN, Inc., Hampton, Virginia CONTRACT NAS-9505 June 997 National Aeronautics and Space Administration Langley Research
More informationMOMENTUM AND HEAT TRANSPORT INSIDE AND AROUND
MOMENTUM AND HEAT TRANSPORT INSIDE AND AROUND A CYLINDRICAL CAVITY IN CROSS FLOW G. LYDON 1 & H. STAPOUNTZIS 2 1 Informatics Research Unit for Sustainable Engrg., Dept. of Civil Engrg., Univ. College Cork,
More informationPIV and LDV measurements behind a backward facing step
PIV and LDV measurements behind a backward facing step M.T. Pilloni, C. Schram, M.L. Riethmulle/^ ^ Mechanical Engineering Department, 09123 Cagliari, Italy ^ von Karman Institute for Fluid Dynamics, 1640
More informationExperiments with Edge Detection using One-dimensional Surface Fitting
Experiments with Edge Detection using One-dimensional Surface Fitting Gabor Terei, Jorge Luis Nunes e Silva Brito The Ohio State University, Department of Geodetic Science and Surveying 1958 Neil Avenue,
More informationWorkhorse ADCP Multi- Directional Wave Gauge Primer
Acoustic Doppler Solutions Workhorse ADCP Multi- Directional Wave Gauge Primer Brandon Strong October, 2000 Principles of ADCP Wave Measurement The basic principle behind wave the measurement, is that
More informationDigital Volume Correlation for Materials Characterization
19 th World Conference on Non-Destructive Testing 2016 Digital Volume Correlation for Materials Characterization Enrico QUINTANA, Phillip REU, Edward JIMENEZ, Kyle THOMPSON, Sharlotte KRAMER Sandia National
More informationMulti-Scale Stereoscopic PIV measurements of a Jet in a Cross-Flow
Multi-Scale Stereoscopic PIV measurements of a Jet in a Cross-Flow Nicolas Lanitis 1,*, James R. Dawson 1 1: Department of Engineering, University of Cambridge, Cambridge, UK * correspondent author: nl258@cam.ac.uk
More informationELEC Dr Reji Mathew Electrical Engineering UNSW
ELEC 4622 Dr Reji Mathew Electrical Engineering UNSW Review of Motion Modelling and Estimation Introduction to Motion Modelling & Estimation Forward Motion Backward Motion Block Motion Estimation Motion
More informationMotion Estimation for Video Coding Standards
Motion Estimation for Video Coding Standards Prof. Ja-Ling Wu Department of Computer Science and Information Engineering National Taiwan University Introduction of Motion Estimation The goal of video compression
More informationExperimental Data Confirms CFD Models of Mixer Performance
The Problem Over the years, manufacturers of mixing systems have presented computational fluid dynamics (CFD) calculations to claim that their technology can achieve adequate mixing in water storage tanks
More informationPredictive Interpolation for Registration
Predictive Interpolation for Registration D.G. Bailey Institute of Information Sciences and Technology, Massey University, Private bag 11222, Palmerston North D.G.Bailey@massey.ac.nz Abstract Predictive
More informationPIV Sleuth: Integrated Particle Image Velocimetry (PIV) Interrogation/Validation Software
PIV Sleuth: Integrated Particle Image Velocimetry (PIV) Interrogation/Validation Software Kenneth T. Christensen Steven M. Soloff Ronald J. Adrian Laboratory for Turbulence and Complex Flow Department
More informationDirect Numerical Simulation of Turbulent Boundary Layers at High Reynolds Numbers.
Direct Numerical Simulation of Turbulent Boundary Layers at High Reynolds Numbers. G. Borrell, J.A. Sillero and J. Jiménez, Corresponding author: guillem@torroja.dmt.upm.es School of Aeronautics, Universidad
More informationParticle Image Velocimetry Part - 3
AerE 545X class notes #5 Particle Image Velocimetry Part - 3 Hui Hu Department of Aerospace Engineering, Iowa State University Ames, Iowa 50011, U.S.A PIV System Setup Particle tracers: Illumination system:
More informationProf. Fanny Ficuciello Robotics for Bioengineering Visual Servoing
Visual servoing vision allows a robotic system to obtain geometrical and qualitative information on the surrounding environment high level control motion planning (look-and-move visual grasping) low level
More informationChapter 6 : Results and Discussion
Refinement and Verification of the Virginia Tech Doppler Global Velocimeter (DGV) 86 Chapter 6 : Results and Discussion 6.1 Background The tests performed as part of this research were the second attempt
More informationA Direct Simulation-Based Study of Radiance in a Dynamic Ocean
A Direct Simulation-Based Study of Radiance in a Dynamic Ocean Lian Shen Department of Civil Engineering Johns Hopkins University Baltimore, MD 21218 phone: (410) 516-5033 fax: (410) 516-7473 email: LianShen@jhu.edu
More informationCOMPUTATIONAL FLUID DYNAMICS ANALYSIS OF ORIFICE PLATE METERING SITUATIONS UNDER ABNORMAL CONFIGURATIONS
COMPUTATIONAL FLUID DYNAMICS ANALYSIS OF ORIFICE PLATE METERING SITUATIONS UNDER ABNORMAL CONFIGURATIONS Dr W. Malalasekera Version 3.0 August 2013 1 COMPUTATIONAL FLUID DYNAMICS ANALYSIS OF ORIFICE PLATE
More informationDevelopment of automated ultraviolet laser beam profiling system using fluorometric technique
Development of automated ultraviolet laser beam profiling system using fluorometric technique BB Shrivastava*, NS Benerji, P Bhatnagar, HS Vora a and U Nundy Chemical and Excimer Laser Section a Laser
More informationGood Practice guide to measure roundness on roller machines and to estimate their uncertainty
Good Practice guide to measure roundness on roller machines and to estimate their uncertainty Björn Hemming, VTT Technical Research Centre of Finland Ltd, Finland Thomas Widmaier, Aalto University School
More informationTutorial: Instantaneous Measurement of M 2 Beam Propagation Ratio in Real-Time
Tutorial: Instantaneous Measurement of M 2 Beam Propagation Ratio in Real-Time By Allen M. Cary, Jeffrey L. Guttman, Razvan Chirita, Derrick W. Peterman, Photon Inc A new instrument design allows the M
More informationCFD modelling of thickened tailings Final project report
26.11.2018 RESEM Remote sensing supporting surveillance and operation of mines CFD modelling of thickened tailings Final project report Lic.Sc.(Tech.) Reeta Tolonen and Docent Esa Muurinen University of
More informationData Partitioning. Figure 1-31: Communication Topologies. Regular Partitions
Data In single-program multiple-data (SPMD) parallel programs, global data is partitioned, with a portion of the data assigned to each processing node. Issues relevant to choosing a partitioning strategy
More information1.12 Optimal Filters (Wiener Filters)
Random Data 75 1.12 Optimal Filters (Wiener Filters) In calculating turbulent spectra we frequently encounter a noise tail just as the spectrum roles down the inertial subrange (-5/3 slope region) toward
More informationHigh spatial resolution measurement of volume holographic gratings
High spatial resolution measurement of volume holographic gratings Gregory J. Steckman, Frank Havermeyer Ondax, Inc., 8 E. Duarte Rd., Monrovia, CA, USA 9116 ABSTRACT The conventional approach for measuring
More informationMeasurement of droplets temperature by using a global rainbow technique with a pulse laser
, 23rd Annual Conference on Liquid Atomization and Spray Systems, Brno, Czech Republic, September 2010 Measurement of droplets temperature by using a global rainbow technique with a pulse laser S. Saengkaew,
More informationThree Dimensional Numerical Simulation of Turbulent Flow Over Spillways
Three Dimensional Numerical Simulation of Turbulent Flow Over Spillways Latif Bouhadji ASL-AQFlow Inc., Sidney, British Columbia, Canada Email: lbouhadji@aslenv.com ABSTRACT Turbulent flows over a spillway
More informationHuman beings are extremely interested in the observation of nature, as this was and still is of utmost importance for their survival.
Historical Background Human beings are extremely interested in the observation of nature, as this was and still is of utmost importance for their survival. (www.copyright-free-images.com) 1 Historical
More informationQuantifying the Dynamic Ocean Surface Using Underwater Radiometric Measurement
DISTRIBUTION STATEMENT A. Approved for public release; distribution is unlimited. Quantifying the Dynamic Ocean Surface Using Underwater Radiometric Measurement Lian Shen Department of Mechanical Engineering
More informationBackward facing step Homework. Department of Fluid Mechanics. For Personal Use. Budapest University of Technology and Economics. Budapest, 2010 autumn
Backward facing step Homework Department of Fluid Mechanics Budapest University of Technology and Economics Budapest, 2010 autumn Updated: October 26, 2010 CONTENTS i Contents 1 Introduction 1 2 The problem
More informationComputational Simulation of the Wind-force on Metal Meshes
16 th Australasian Fluid Mechanics Conference Crown Plaza, Gold Coast, Australia 2-7 December 2007 Computational Simulation of the Wind-force on Metal Meshes Ahmad Sharifian & David R. Buttsworth Faculty
More informationTop Layer Subframe and Node Analysis
Top Layer Subframe and Node Analysis By Paul Rasmussen 2 August, 2012 Introduction The top layer of the CCAT backing structure forms a critical interface between the truss and the primary subframes. Ideally
More information[HALL PROBE GRADIOMETRY ]
2008 [HALL PROBE GRADIOMETRY ] A novel Scanning Hall probe gradiometer has been developed and a new method to image x, y & z components of the magnetic field on the sample surface has been demonstrated
More informationCalculate a solution using the pressure-based coupled solver.
Tutorial 19. Modeling Cavitation Introduction This tutorial examines the pressure-driven cavitating flow of water through a sharpedged orifice. This is a typical configuration in fuel injectors, and brings
More informationACTIVE SEPARATION CONTROL WITH LONGITUDINAL VORTICES GENERATED BY THREE TYPES OF JET ORIFICE SHAPE
24 TH INTERNATIONAL CONGRESS OF THE AERONAUTICAL SCIENCES ACTIVE SEPARATION CONTROL WITH LONGITUDINAL VORTICES GENERATED BY THREE TYPES OF JET ORIFICE SHAPE Hiroaki Hasegawa*, Makoto Fukagawa**, Kazuo
More informationQuantifying Three-Dimensional Deformations of Migrating Fibroblasts
45 Chapter 4 Quantifying Three-Dimensional Deformations of Migrating Fibroblasts This chapter presents the full-field displacements and tractions of 3T3 fibroblast cells during migration on polyacrylamide
More informationCFD MODELING FOR PNEUMATIC CONVEYING
CFD MODELING FOR PNEUMATIC CONVEYING Arvind Kumar 1, D.R. Kaushal 2, Navneet Kumar 3 1 Associate Professor YMCAUST, Faridabad 2 Associate Professor, IIT, Delhi 3 Research Scholar IIT, Delhi e-mail: arvindeem@yahoo.co.in
More informationTracePro Stray Light Simulation
TracePro Stray Light Simulation What Is Stray Light? A more descriptive term for stray light is unwanted light. In an optical imaging system, stray light is caused by light from a bright source shining
More informationFlow Field of Truncated Spherical Turrets
Flow Field of Truncated Spherical Turrets Kevin M. Albarado 1 and Amelia Williams 2 Aerospace Engineering, Auburn University, Auburn, AL, 36849 Truncated spherical turrets are used to house cameras and
More informationAP* Optics Free Response Questions
AP* Optics Free Response Questions 1978 Q5 MIRRORS An object 6 centimeters high is placed 30 centimeters from a concave mirror of focal length 10 centimeters as shown above. (a) On the diagram above, locate
More informationDAMAGE INSPECTION AND EVALUATION IN THE WHOLE VIEW FIELD USING LASER
DAMAGE INSPECTION AND EVALUATION IN THE WHOLE VIEW FIELD USING LASER A. Kato and T. A. Moe Department of Mechanical Engineering Chubu University Kasugai, Aichi 487-8501, Japan ABSTRACT In this study, we
More informationVISUALIZATION OF MULTI-SCALE TURBULENT STRUCTURE IN LOBED MIXING JET USING WAVELETS
9TH. INTERNATIONAL SYMPOSIUM ON FLOW VISUALIZATION, 000 VISUALIZATION OF MULTI-SCALE TURBULENT STRUCTURE IN LOBED MIXING JET USING WAVELETS Hui LI, Hui HU, Toshio KOBAYASHI, Tetsuo SAGA, Nobuyuki TANIGUCHI
More informationDirect Measurements of Reynolds Stresses and Turbulence in the Bottom Boundary Layer
LONG-TERM GOALS Direct Measurements of Reynolds Stresses and Turbulence in the Bottom Boundary Layer Joseph Katz Department of Mechanical Engineering The Johns Hopkins University Baltimore, MD 21218 katz@titan.me.jhu.edu
More informationMinimizing Noise and Bias in 3D DIC. Correlated Solutions, Inc.
Minimizing Noise and Bias in 3D DIC Correlated Solutions, Inc. Overview Overview of Noise and Bias Digital Image Correlation Background/Tracking Function Minimizing Noise Focus Contrast/Lighting Glare
More informationHui Hu Department of Aerospace Engineering, Iowa State University Ames, Iowa 50011, U.S.A
AerE 311L & AerE343L Lecture Notes Lecture # 14: Advanced Particle Image Velocimetry Technique Hui Hu Department of Aerospace Engineering, Iowa State University Ames, Iowa 511, U.S.A Particle-based techniques:
More informationFLOW STRUCTURE AROUND A HORIZONTAL CYLINDER AT DIFFERENT ELEVATIONS IN SHALLOW WATER
FLOW STRUCTURE AROUND A HORIZONTAL CYLINDER AT DIFFERENT ELEVATIONS IN SHALLOW WATER 1 N.Filiz OZDIL, * 2 Huseyin AKILLI 1 Department of Mechanical Engineering, Adana Science and Technology University,
More informationContinued Investigation of Small-Scale Air-Sea Coupled Dynamics Using CBLAST Data
Continued Investigation of Small-Scale Air-Sea Coupled Dynamics Using CBLAST Data Dick K.P. Yue Center for Ocean Engineering Department of Mechanical Engineering Massachusetts Institute of Technology Cambridge,
More informationADVANCED MEASUREMENT TECHNIQUES IN HYDRODYNAMICS. Chittiappa Muthanna
ADVANCED MEASUREMENT TECHNIQUES IN HYDRODYNAMICS Chittiappa Muthanna Outline Why use these techniques? Constant temperature anemometry Laser Doppler Velocimetry PIV Measuring Shapes and Deformations 2
More informationChapter 3 Image Registration. Chapter 3 Image Registration
Chapter 3 Image Registration Distributed Algorithms for Introduction (1) Definition: Image Registration Input: 2 images of the same scene but taken from different perspectives Goal: Identify transformation
More informationContrast Optimization A new way to optimize performance Kenneth Moore, Technical Fellow
Contrast Optimization A new way to optimize performance Kenneth Moore, Technical Fellow What is Contrast Optimization? Contrast Optimization (CO) is a new technique for improving performance of imaging
More informationMeasurement Principle What is VidPIV? An Overview of the VidPIV Interface Data Management
Measurement Principle What is VidPIV? An Overview of the Data Management Contents General Overview PIV Project 1 2 3 4 Measuring Principle 3 What is VidPIV? 4 An Overview of the 5 Data Management 6 4.1
More informationPhysical Optics. You can observe a lot just by watching. Yogi Berra ( )
Physical Optics You can observe a lot just by watching. Yogi Berra (1925-2015) OBJECTIVES To observe some interference and diffraction phenomena with visible light. THEORY In a previous experiment you
More informationFLUOSTAR Rhodamine B- encapsulating microspheres are seeding particles optimized for Particle Image Velocimetry. Technical handbook ver.
www.ebm.vc FLUOSTAR Rhodamine B- encapsulating microspheres are seeding particles optimized for Particle Image Velocimetry Technical handbook ver.1, Feb, 2010 CONTENTS 1) Introduction of fluorescent PIV
More informationSimulation of Flow Development in a Pipe
Tutorial 4. Simulation of Flow Development in a Pipe Introduction The purpose of this tutorial is to illustrate the setup and solution of a 3D turbulent fluid flow in a pipe. The pipe networks are common
More informationFlat-Plate As stated earlier, the purpose of the flat-plate study was to create an opportunity for side-by-side comparison of ÒfastÓ RNG and
Chapter Six: Comparison of Turbulence Models Performance Comparisons There are a number of criteria by which to judge turbulence models. One criterion sometimes important to mathematicallyminded model
More informationLarge Eddy Simulation of Flow over a Backward Facing Step using Fire Dynamics Simulator (FDS)
The 14 th Asian Congress of Fluid Mechanics - 14ACFM October 15-19, 2013; Hanoi and Halong, Vietnam Large Eddy Simulation of Flow over a Backward Facing Step using Fire Dynamics Simulator (FDS) Md. Mahfuz
More informationindex of refraction-light speed
AP Physics Study Guide Chapters 22, 23, 24 Reflection, Refraction and Interference Name Write each of the equations specified below, include units for all quantities. Law of Reflection Lens-Mirror Equation
More informationAccurately measuring 2D position using a composed moiré grid pattern and DTFT
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 Accurately measuring 2D position using a composed moiré grid pattern and DTFT S. Van
More informationOn the flow and noise of a two-dimensional step element in a turbulent boundary layer
On the flow and noise of a two-dimensional step element in a turbulent boundary layer Danielle J. Moreau 1, Jesse L. Coombs 1 and Con J. Doolan 1 Abstract This paper presents results of a study on the
More informationAnalysis of Cornell Electron-Positron Storage Ring Test Accelerator's Double Slit Visual Beam Size Monitor
Analysis of Cornell Electron-Positron Storage Ring Test Accelerator's Double Slit Visual Beam Size Monitor Senior Project Department of Physics California Polytechnic State University San Luis Obispo By:
More informationANALYSIS AND MEASUREMENT OF SCARF-LAP AND STEP-LAP JOINT REPAIR IN COMPOSITE LAMINATES
16 TH INTERNATIONAL CONFERENCE ON COMPOSITE MATERIALS ANALYSIS AND MEASUREMENT OF SCARF-LAP AND STEP-LAP JOINT REPAIR IN COMPOSITE LAMINATES David H. Mollenhauer*, Brian Fredrickson*, Greg Schoeppner*,
More informationMode-Field Diameter and Spot Size Measurements of Lensed and Tapered Specialty Fibers
Mode-Field Diameter and Spot Size Measurements of Lensed and Tapered Specialty Fibers By Jeffrey L. Guttman, Ph.D., Director of Engineering, Ophir-Spiricon Abstract: The Mode-Field Diameter (MFD) and spot
More information1. Estimation equations for strip transect sampling, using notation consistent with that used to
Web-based Supplementary Materials for Line Transect Methods for Plant Surveys by S.T. Buckland, D.L. Borchers, A. Johnston, P.A. Henrys and T.A. Marques Web Appendix A. Introduction In this on-line appendix,
More informationFinite Element Method. Chapter 7. Practical considerations in FEM modeling
Finite Element Method Chapter 7 Practical considerations in FEM modeling Finite Element Modeling General Consideration The following are some of the difficult tasks (or decisions) that face the engineer
More informationTwo-dimensional Gaussian regression for sub-pixel displacement estimation in particle image velocimetry or particle position
Experiments in Fluids (25) 38: 5 55 DOI.7/s348-5-942-3 ORIGINALS H. Nobach Æ M. Honkanen Two-dimensional Gaussian regression for sub-pixel displacement estimation in particle image velocimetry or particle
More information