Particle Image Velocimetry
Laboratory for Applied Fluid
Dynamics
Northwestern University, Evanston, Illinois
Particle Image
Velocimetry (PIV) is a non-intrusive, two-dimensional, instantaneous velocity
measurement technique. The measurements are performed by seeding the flow of
a fluid with small particles and illuminating a region of interest with a thin
sheet of powerful laser light. Then, sequential images of illuminated particles
are captured either in pairs or in a double-exposed form by a CCD camera, depending
on the method that the images are interrogated. In the case where the images
are captured in pairs one obtains the velocity field via cross-correlation,
whereas in the other case via auto-correlation. Below, the cross-correlation
operation is described briefly by using numerically simulated PIV images that
were created with Matlab. Image 1 and Image 2 show the simulated particle images
at time instants t = t_1 and t = t_2, respectively. Particles are shifted by
4 pixels in positive x and y directions as they move from Image 1 to Image 2.
|
|
The average displacement, and consequently
the average velocity, of the particles in the images are obtained by performing
cross-correlation, which is a statistical method to determine the relation between
the data contained in the images. As mentioned above, the prescribed relation
between data contained in Image 1 and Image 2 is that particles move 4 pixels
in positive x and y directions, so the resulting correlation function has a
maximum value at this particular location. Below, contour and three-dimensional
plots of the correlation function are shown. The peak is shifted by 4 pixels
in each direction with respect to the center of the plotted area.
|
|
In the case of an actual PIV image however, the number of particles and the size of the image are much larger than the numerically simulated images shown above. Hence the image is divided into many small regions from which one velocity vector representing the average velocity of the particles contained only in that small region is obtained. This can be best described by an actual image and measured velocity field.
The image below is a real image of
wavy vortex flow in a cylindrical Couette device investigated in our laboratory.
The bright arcs at the top and the bottom (the bottom one is not as visible)
are the outer and inner cylinder boundaries that become visible as the laser
light sheet passes through the flow region of interest. The inner cylinder rotates
in the CW direction and the outer cylinder is stationary. The white spots are
individual tracer particles in the flow. The size of the image is 1K x 1K pixels.
The average velocity of the tracer particles in a small 64 x 128 pixel interrogation
region (subset of the image) is found using the cross-correlation for this region.
Repeating this process over a grid of locations in the flow field provides an
array of green velocity vectors. The length of the green vectors is proportional
to the magnitude of the velocity. The orientation of the green vectors shows
the direction of the flow. The velocity field shown below is obtain from a single
pair of images. Ensemble averaging several realizations of velocity data produces
even better results.
|