Quantitative \muPIV Measurements of Velocity Profiles

Quantitative PIV Measurements of Velocity Profiles

Bryant  P. W. Neumann  R. F. Moura  M. J. B. Steiner IBM Research – Brazil, Av. Pasteur 138 & 146, Urca, Rio de Janeiro, CEP 22290-240, Brazil  M. IBM Research – Brazil, Av. Pasteur 138 & 146, Urca, Rio de Janeiro, CEP 22290-240, Brazil Carvalho Dept. Mech. Eng., PUC–Rio, R. Marquês de São Vicente, 225, Gávea, Rio de Janeiro, Brazil  M. S. Dept. Mech. Eng., PUC–Rio, R. Marquês de São Vicente, 225, Gávea, Rio de Janeiro, Brazil Feger IBM Research – Brazil, Av. Pasteur 138 & 146, Urca, Rio de Janeiro, CEP 22290-240, Brazil  C. IBM Research – Brazil, Av. Pasteur 138 & 146, Urca, Rio de Janeiro, CEP 22290-240, Brazil
July 5, 2019

In Microscopic Particle Image Velocimetry (PIV), velocity fields in microchannels are sampled over finite volumes within which the velocity fields themselves may vary significantly. In the past, this has limited measurements often to be only qualitative in nature, blind to velocity magnitudes. In the pursuit of quantitatively useful results, one has treated the effects of the finite volume as errors that must be corrected by means of ever more complicated processing techniques. Resulting measurements have limited robustness and require convoluted efforts to understand measurement uncertainties. To increase the simplicity and utility of PIV measurements, we introduce a straightforward method, based directly on measurement, by which one can determine the size and shape of the volume over which moving fluids are sampled. By comparing measurements with simulation, we verify that this method enables quantitative measurement of velocity profiles across entire channels, as well as an understanding of experimental uncertainties. We show how the method permits measurement of an unknown flow rate through a channel of known geometry. We demonstrate the method to be robust against common sources of experimental uncertainty. We also apply the theory to model the technique of Scanning PIV, which is often used to locate the center of a channel, and we show how and why it can in fact misidentify the center. The results have general implications for research and development that requires reliable, quantitative measurement of fluid flow on the micrometer scale and below.

1 Introduction

Microscopic Particle Image Velocimetry (PIV) has become an invaluable measurement technique for probing fluid flow in systems at the micro- and the nano-scale [1, 2, 3, 4]. Well-controlled and understood PIV experiments performed on increasingly complex systems will become crucial to maintaining the rapid progress observed in the biomedical and natural resources industries [5, 6, 7]. When imaging only over a region small compared to the size of a microfluidic system, it is sometimes possible optically to restrict the volume over which the system is probed. The result is often a quantitative measurement of fluid flow in that region [8, 9]. To image simultaneously as many flow features as possible, however, one requires a field of view of size comparable to the length scales of the system. Unfortunately, this leads inevitably to having also a large depth of field along the optical axis, which is an issue [10, 11, 12] that hinders a straightforward quantification and interpretation of results. For PIV to reach its full potential as an investigative and design tool, it must provide repeatable and quantitatively useful results. It has been an open question whether or not this is possible for PIV with a depth of field of size comparable to that of the microfluidics system it measures.

There are two standard approaches to treating the effects of a large depth of field. The first is to limit PIV to qualitative comparisons of velocity profile shapes, often away from channel walls. In such cases, one treats the volumetric flow rate of the fluid as a free parameter used to match expected profiles to the measured data [10, 11]. In the second approach, one treats finite volume effects as measurement errors, particularly near the walls, and attempts to remove as many of these effects as possible via data processing. This approach has, to some extent, enabled the quantitative evaluation of PIV data [12, 13, 14, 15].

A shortcoming of the latter method, however, is the often complicated and convoluted data processing that limits or even precludes the physical understanding of the fluid flow and its analysis. It bears the risk of over-fitting results to the extent that important flow features are missed, or even misdiagnosing the causes of observed phenomena. An example is the spurious, increasing velocity near the channel walls that sometimes appears [13, 15]. In most cases noisy or uncertain experimental data are excluded from the measured flow profile altogether, without physical justification. A related problem is that excessive data processing introduces additional sources of experimental uncertainty in an obscure fashion.

In this article, we investigate the results of a series of experiments performed with a large depth of field PIV, in which we image across the entire microchannel and measure the velocity profile from wall to wall. Rather than use data processing to remove the effects of finite measurement volumes, we introduce the concept of Sampling Volume (SV) over which the velocity field can be averaged to recover the measured velocity profiles. We then systematically determine the size and shape of the SV using an error minimization procedure. For our experiment, which is described below, we find that measurements near the walls reflect in a straightforward manner the non-trivial intersection of the SV with the edges of the microchannel. Consequently, the entire measured profile contains information useful for determining the SV, and one should therefore not discard features measured near channel walls.

We begin with a discussion of PIV for large depth of field and a motivation for our definition of the SV. We then describe our experimental setup and how our data are processed and analyzed. We also discuss the expected velocity field in the channel, and its relevant sources of uncertainty. Then we show how averaging the velocity field over a rectangular SV reproduces to within our experimental uncertainty the measured velocity profiles across the entire channel. We also demonstrate how, because of the simplicity of our method and the minimal post-processing of data, PIV can in some cases be used to infer an unknown velocity field or flow rate in a channel of known cross section. We close with a series of use cases for our approach, which include robustly measuring velocity profiles in the presence of camera misalignment and very noisy data that would typically be deemed useless, a demonstration of how data processing affects the SV, and finally a critical analysis of the procedure known as “Scanning PIV,” which is often used to place a measurement domain at the center of microchannels [12, 16]. We show why, for finite-sized SV, Scanning PIV does not necessarily locate the center of a microchannel, and should therefore be used with caution.

2 Particle Image Velocimetry

2.1 General description

Particle Image Velocimetry (PIV) is an optical method for visualizing and measuring flow fields. One seeds a fluid with tracer particles, and if the particles follow the fluid flow without disturbing it significantly, by analyzing their trajectories one can infer the velocity field of the moving fluid. To determine the particle displacements one uses consecutive snapshots of the system taken at times and , containing several particle images at various positions. The snapshots are subdivided into smaller regions called “interrogation windows” and the net particle displacement is calculated for each one of them. Finally, the velocity representative of each window is calculated simply by dividing the respective displacement by the interval, , between frames.

Analyzing and interpreting a PIV measurement is complicated for the following reasons. First, the velocity obtained from this process is an average in space and time, and its accuracy depends on the flow field itself, as well as on both spatial and temporal resolution of the experimental setup [17]. Second, tracer particles are susceptible to Brownian motion, and as a consequence their positions in images must be correlated within sub-ensembles to suppress the effect of random movements and determine net particle displacements. The application of cross-correlation techniques ultimately provides a 2D projection of the 3D velocity field of the fluid located within the imaged region, as viewed from the camera [4].

In Microscopic PIV (PIV), see Figure 1(a), fluorescent tracer particles in a fluid channel of arbitrary geometry are observed by means of an optical microscope setup. The size of the tracer particles needs to be chosen carefully such that particles emit a sufficient amount of light for detection without perturbing the flow. To perform the measurement, one focuses a laser pulse on the fluid channel to excite the tracer particles. The particles then emit light at a red-shifted fluorescence wavelength that is collected by the same objective lens, alongside with the reflected laser light. The reflected excitation laser light is then filtered out by a dichroic beam splitter while the fluoresced light passes, rendering fluorescence images of spatially resolved particles on the CCD array. As indicated in Figure 1(a), fluorescent particles (red) in the object plane are spatially resolved in the conjugate image plane with a lateral resolution in and that is determined by the specifics of the optical setup. However, a particle’s image will only be detectable if its -position falls within the finite range labeled in Figure 1(b), which depends not only on the characteristics of the optical setup but also on the data processing algorithms used in the analysis. In the following section, we will introduce a concept for the determination of the volume probed by a PIV experiment.

Figure 1: (a) Simplified schematic of the PIV setup. The flow direction of fluorescent particles in the fluidic channel (FC) is indicated by an arrow. MO: microscope objective; BS: dichroic beam splitter; L: lens; CCD: charge coupled device array detector. (b) Schematic of the fluidic channel showing the depth of the SV, . (c) Scanning electron microscope image of the cross section of the fluid channel used in the flow experiments. The length of the white scale bar represents 10 m.

2.2 Sampling Volume

In this paper, our main concern is a well-known complication inherent to PIV with a large depth of field, for which tracer particles are imaged over a finite volume with a characteristic length scale on the order of the scale of the microchannels themselves. Besides the well established depth of field and depth of focus [4, 18], several other quantities have been used to quantify the depth over which the system is probed. Among them are the laser light sheet [4, 10], the depth of measurement [19], and the depth of correlation [11, 12, 20]. Although they have various meanings and account for various effects, many of them derive from the basic definition of the depth of field.

The term “laser light sheet” [4, 10] was introduced to denote the illuminated region inside the fluid system. In most PIV applications the light sheet and the entire fluid system are of comparable sizes. When that is the case, all particles inside the channel are excited by the laser light and, in general, all fluorescent particles within the field of view of the optical setup and within the depth of field, , create in-focus images on the CCD detector. The geometrical shape of this volume can be approximated as a parallelepiped, and the focal -plane oriented perpendicular to the optical axis cuts this volume in half.

In standard far-field fluorescence microscopy the total point spread function and, hence, the size of the focal volume is determined by the product of the excitation point spread function and the detection point spread function [18]. In PIV, however, where image processing is used to determine the particles’ displacements over successive frames, there are additional contributions that alter the effective volume that is being probed, such as the application of the image pair correlation function [21]. As a consequence, a “depth of correlation” (DOC) was introduced to estimate the depth of the volume that contains those particles contributing significantly to the spatial correlation function [11, 12, 20]. The calculated DOC is useful for identifying the experimental parameters relevant to the probed volume. However, it is based on an assumed intensity threshold beyond which imaged particles do not “cause a significant bias in the velocity estimate” [15]. It is thus not quantitatively well defined, and has not led to conclusive results in comparison with experiments [11, 12]. In Section 6.3 we will present a series of results that indicate the DOC should be used only as an order-of-magnitude estimate of the probed volume.

Despite significant effort to determine precisely the probed volume, a simultaneously intuitive and quantitative definition is lacking. To address this issue, we introduce a quantity with a straightforward and verifiable interpretation. We call it Sampling Volume (SV), and it is designed to account for the combined effects of illumination scheme, optical depth, and image processing. The concept, as opposed to previously defined quantities such as the DOC, is strikingly simple: the SV is defined as the volume over which one averages the velocity field of the fluid in order to recover the measured profile.

We begin with a complete description of our averaging process for one experiment performed on a simple, straight channel, for which we will ignore any variation in measured velocity down the channel ( direction). We denote the velocity field within the channel by . Because of the translational symmetry in , we take . The velocity averaged over the SV, which is compared with measurements, is then


as a function of , which is the direction across the channel, as seen below in Figure 3. In (1), the integration bounds for a rectangular SV are defined by




In (2) and (3), and denote the maximum and minimum values of the arguments. The bottom and top locations of the channel itself are and , respectively. The vertical center of the SV, as measured with respect to the channel center, is denoted by , and the vertical depth, or thickness, of the SV is represented by . The selection of the maximum or minimum values in the integration bounds is necessary to keep from continuing the average beyond the channel’s extent, where no fluid exists. In Section 6.4 is a demonstration of how failing to cut off the average properly at the edge of the channel can lead to a calculated velocity of limited applicability. Throughout Sections 4 and 6 are several demonstrations of the success of using Equation (1) to explain measurements.

3 Experimental Methods and Data Analysis

3.1 Flow Field Measurements

In this section, we briefly describe the experimental setup used for performing all the fluid flow experiments discussed in this article. Further description of the experimental setup can be found in Ref. [22]. The experiments were performed using a commercially available PIV system (TSI Inc.), composed of an optical microscope, CCD camera, two lasers, and a laser pulse synchronizer. The system is operated by a data acquisition and management software (Insight 3G, TSI Inc.). The inverted microscope (IX71S1F-3, Olympus) is equipped with a 10x/0.30 air objective (UPlanFL N, Olympus). The Peltier-cooled CCD (Sensicam 630166, PowerView) has a resolution of 1376 1024 pixels (1.4MP) and a 2x projection lens. The two pulsed Nd:YAG lasers (Gemini PIV-15, NEW WAVE) operate at 532 nm and a frequency of 15 Hz, providing a pulse energy of 50 mJ. From the specifications of the optical setup we can calculate our depth of field [18] to be approximately 20 m. The laser pulse synchronizer (610034, TSI Inc.) controls the timing sequence of laser pulses and exposure time of the CCD. For each flow experiment, 100 pairs of images were collected with a time interval, , of 500 s between successive images.

A glass microfluidic device manufactured by Dolomite Centre Ltd, UK was used to perform the flow experiments. The device has a straight channel with an elliptical cross-section, with a semi-minor axis of 50 m and semi-major axis of 55 m. To visualize the flow we use a mix of 14% aqueous solution of fluorescent particles of 1 m in diameter (ex/em 542/612 nm, Thermo Scientific Fluoro-Max Dyed Red Aqueous Fluorescent Particles) and 86% purified water (Thermo Scientific Nalgene Analytical Filter, 0.2 m). This fluid mix is injected in the device with the use of a syringe pump (11 Elite syringe pump 704501, Harvard Apparatus). The injection rates ranged from 75 l/h to 100 l/h.

3.2 Measurement Uncertainty

The extent to which we can determine the size and shape of our SV depends on our experimental uncertainties, as will be demonstrated in Sections 4 and 6. We therefore performed several calibration procedures to assess the accuracy and precision of our measurements.

The major sources of uncertainty in the fluid flow itself are in the microchannel geometry and the flow rate provided by the syringe pump. With a scanning electron microscope, we imaged representative cross sections of our channels. Figure 1 contains a sample image. We digitized the boundary from two separate images and found that all digitized points lie within the region between two ellipses with semi-major and semi-minor axes of and , respectively. The syringe pump was calibrated by measuring the amount of fluid injected over time. Five measurements were taken for each of three injection rates. The measurements were also compared with the total infused volume and the time elapsed, as given on the pump’s display. We obtained a 3% uncertainty in the volumetric flow rate, which is consistent with the uncertainty of 0.5% provided by the manufacturer. For our results in Sections 4 and 6, we assume an uncertainty of 0.5% of the flow rate. We also tested the syringe pump connected to the microfluidic system to check for internal pressure effects. Again we measured the infused volume over time. No significant change in the injection rates was observed.

To address the uncertainty in lengths measured by the PIV system, we calibrated the camera’s field of view using a standard calibration target. For our experimental setup the field of view is (4405) m (3255) m, in agreement with the equipment’s specifications.

3.3 Data Processing

Image processing of the results presented in Section 4 was carried out using the Insight 3G software, starting from the raw images. We used an interrogation window of 32 32 pixels, with a 50% overlap according to the Nyquist sampling criterion. The FFT correlator and the 3-point Gaussian peak detection were employed in the determination of particle displacements. In order to increase the signal-to-noise ratio, we subtracted background noise and performed an ensemble average over correlation functions across 100 pairs of snapshots. From the camera’s pixel resolution and the field of view we know that one pixel corresponds to 0.32 m in both the and the direction. The overall resolution is therefore 5.12 m. With this resolution we obtain 85 velocity vectors along the channel length and 21 along the channel width. Because of the low Reynolds number, the resulting flow is as expected laminar and stable. Figure 2(a) shows a typical velocity profile, which will be further analyzed in Sec. 4.

Velocity profiles were also obtained independently starting from the same raw images and analyzing them with the OpenPIV [23] software 111http://www.openpiv.net. The processing options for OpenPIV, such as the geometry of the interrogation windows, the cross-correlation method, and the peak detection algorithm, were the same as those employed by Insight 3G. The difference was in the noise reduction algorithms: instead of performing an ensemble average over correlation functions, we averaged the flow field themselves over the 100 pairs of snapshots, and we did not perform background noise subtraction. The results of the two processing methods are discussed in Section 6.3.

Figure 2: Quiver plots of the measured velocity field as processed by (\subreffig:quiver_I3G) Insight 3G and (\subreffig:quiver_OPIV) OpenPIV.

4 Results

For comparison with theory we averaged the profiles located along the (streamwise) direction. The resulting average profile for an injection rate of 100 l/h is pictured in Figure 3(\subreffig:cxfit1) as dots with vertical error bars corresponding to one standard deviation. Horizontal error bars representing the uncertainty in length calibration discussed in Section 3.2 are too small to be seen. The large speeds with large error bars beyond the channel walls are numerical artifacts from data processing.

Figure 3: (\subreffig:cxfit1) Measured and theoretically expected velocity profiles for an injection rate of 100 l/h. (\subreffig:cxfit2) The magnitude of the simulated velocity field, , is overlaid by the SV (shaded rectangle) that results in the fit. The dashed line in both images corresponds to a centered SV of zero thickness.

We used Equation (1) to calculate for a given thickness, , and position, , of the SV. To compute the velocity field, , we used the Lattice Boltzmann Method [24], to which is input the injection rate of the pump, the channel geometry, and the fluid properties. To determine the position, , and the thickness, , of the SV that results in the measured profile, we searched for possible fits to the data by varying and while keeping fixed. In other words, we did not vary the pressure drop or the injection flow rate, and thus the magnitude of the theoretically expected velocity field is not a free parameter. In Figure 3(\subreffig:cxfit1), the solid blue line is calculated for an elliptical channel with semi-major and semi-minor axes of and , respectively, and a flow rate delivered by the pump of 100 l/h. The shaded area represents the combined effect of the uncertainty in the channel shape and the flow rate, as discussed in Section 3.2. The dashed line, included for reference, is the calculated velocity profile at the center of the channel, .

The fit in Figure 3 is for and . The position can be positive or negative because the channel has a reflection symmetry with respect to a horizontal plane bisecting it. The lower location is depicted in Figure 3(\subreffig:cxfit2). These two parameters that determine the disposition of the SV were found by sweeping through a range of values, for which Figure 4 shows the residual, . It is clear that no other choices for the center and depth of a rectangular SV can fit the data as well.

Figure 4: Dependence of the residual, , with respect to the fit parameters (the SV center) and (the SV thickness). The white spots around the minima depict regions in which is no more than 20% larger than its lowest value.

Based on this fit, we conclude that, to within our experimental uncertainty, for this test there is a rectangular SV over which one can average the velocity field to match the measured velocity profile. The result is a quantitatively useful PIV measurement of a velocity profile across the entire width of a channel. In Section 6 are several more examples of successful fits to experimental results. It has been claimed recently [12], however, that the simple average over a sampling region cannot reproduce velocity profiles obtained by PIV with a large depth of field. Our analysis in Section 6.4 indicates that this claim may result from the incorrect use of a technique called Scanning PIV to locate the vertical center of a channel.

5 Inferring Flow Rate, Channel Geometry, or Velocity Field

In Section 4, we demonstrated how to determine the SV for a system in which the flow rate, channel geometry, and fluid properties are known. An important benefit of minimizing the post-processing of data and instead working with the SV is that one may be able to reverse the analysis and infer either the flow rate, the velocity field, or the channel geometry in systems where one of them is not known. For example, a common approach when one uses PIV to measure the flow rate of, for example, blood in an artery or oil in porous media [7], is to measure the velocity profile and then match the maximum velocity to the peak of the theoretically expected profile, such as the Hagen-Poiseuille profile for cylindrical channels. As is clear from our analysis of the SV, however, this approach will underestimate the experimental flow rate.

One can already obtain a better evaluation from the same data by acknowledging the presence of the SV. If one first calibrates the PIV apparatus for a well-controlled experiment in which the thickness, , of the SV is determined for a given set of optical and processing parameters, as described in Section 4, one can use the knowledge of to understand the degree to which a measured velocity profile will underestimate the actual flow rate.

Perhaps more interestingly, by retaining all measured features, such as the kinks near the walls in the profiles seen in Figure 3(\subreffig:cxfit1), for some channel shapes one can determine with a single experiment the size and position of the SV, without fitting to the measured profile. As shown in Figure 5 for our experiment, the values where the measured profile goes to zero uniquely determine where the top of the SV intersects the channel’s walls, where the velocity field itself is zero. Similarly, the profile’s kinks determine where the bottom of the SV intersects the walls. With and uniquely determined, and given the measured velocity magnitudes, one can infer the flow rate or even the velocity field of the flowing fluid. For example, if we had known the channel geometry but not the flow rate for the measurement described above, we could have recovered, to within experimental uncertainty, the entire velocity field by performing several simulations with various flow rates.

Figure 5: Correspondence between the intersections of the channel walls with the SV and the position of the near-wall features.

In this scenario, one must know the fluid properties and have enough knowledge of possible flow rates to choose correctly the PIV parameters such as the between images. For some channel shapes and SV sizes, a unique determination of and may not be possible from a single measured profile. In that case, the theory can still be used to identify a set of possible solutions and to judge the relative likelihood of the possibilities. Furthermore, by parallel reasoning, one can reverse the process: If one knows the flow rate but not the channel geometry, one can in principle iterate over possible channel shapes to infer as much as possible about the unknown geometry of the microchannel itself.

Another intriguing possibility is to verify flow models for complex fluids that present non-Newtonian behavior, for which theoretical models or numerical solutions may not be accurate. In this case, one must either determine the SV as in Figure 5, or else by means of a calibration experiment using known optical parameters and a well-understood fluid, such as we have described in Section 4. Then one can measure the velocity profile of the complex fluid, and the result will be a PIV measurement that can provide quantitative insight into the various theoretical treatments of the complex fluid.

6 Use Cases

Here we describe several use cases for the simple theory of the SV presented in Section 4. These use cases demonstrate its robustness against common sources of experimental error and against noisy data, its utility for quantifying the effect of processing variables, and how the theory has led to a deeper understanding of the Scanning PIV technique.

6.1 Robustness Against Camera Misalignment

In some experiments, we measured a consistent and smooth drop in the velocity along the direction, which is the direction of flow. One example can be seen in Figure 6, for the speed at the center (in ) of the channel as it changes in the direction. Once an experiment was prepared and fluid was flowing, the change in velocity down the channel was consistent, even for other locations along the length of the channel. In other experiments, such as the one presented in Figure 3, there was no appreciable change in velocity down the channel. The simplest explanation for this occasional smooth change in velocity is an occasional misalignment between the camera and the microchannel.

Figure 6: The points represent the speed measured at the center of the cross section for several positions along , the direction of flow. The solid line shows the theoretically expected speed, calculated at the center of the profile using Equations (4), (5), and (6).

With our method, one can easily treat a misaligned camera. The translational symmetry of the channel is unaffected, so we continue to use for the velocity field inside the channel. The position and shape of the SV do depend on the optical system, however, so we model a misalignment of the camera by varying in the position of the center of the SV: . The result is an expected measured velocity that also varies in , , which is calculated as


Here the integration bounds must be defined as




For the fit in Figure 6 the thickness of the SV, , is a constant . The center of the SV, , varies linearly from to . The channel depth is , so the fit in Figure 6 suggests a camera misalignment of approximately , which is reasonable for our experimental setup. Even though the camera was misaligned, we were able to verify the flow in the microchannel.

6.2 Robustness Against Noisy Data

In the past, when normalized velocities have been measured in a more qualitative fashion [10, 11], interesting features in velocity profiles, such as the kinks highlighted in Figure 5, have typically been ignored. In the more quantitative evaluations [12, 14, 15], one has tried to eliminate any kinks or other features near walls via post-processing meant to reduce the effects of a finite SV.

With Figure 7 we demonstrate how maintaining such features can in fact lead to more robust PIV measurements, particularly when one has noisy data.

Figure 7: Flow profile for an experiment in which some of the velocity vectors could not be determined correctly. The pump rate for this experiment was 75 l/h. (\subreffig:qhole) Measured flow field in which at least four of the profiles have incorrect velocities near the center of the channel. (\subreffig:fhole) Fit to the data, demonstrating the robustness of our method. The large error bars near the center of the profile result from the incorrect vectors.

In this experiment, raw particle images included a large, stationary bright spot that resulted, after processing, in the area with zero velocity shown in Figure 7(\subreffig:qhole). The average profile, shown in Figure 7(\subreffig:fhole), has large uncertainty in the center of the channel. Rather than discard the results of the experiment, by fitting the entire measured speed profile, from wall to wall, we were able to quantify the flow in the microchannel. Because we fit the entire profile, and thereby incorporate more data that have been minimally processed, we believe that quantitative evaluation using the SV is significantly more robust and more straightforward than with other methods.

6.3 Varying Experimental and Processing Parameters

In this section we show fits from a variety of other tests, from which can be seen how variables such as flow rate, focus, and processing options can affect the size, , and position, , of the SV. We also discuss how the measured SV differs from the calculated DOC [20], and how the two have complimentary roles. As in Section 4, error bars represent one standard deviation in the set of all profiles measured along the direction. The blue shaded area represents the combined uncertainty in the actual fluid flow, which results from our uncertainty in the channel cross section and the flow rate of the pump. Note that, to within our experimental uncertainties, a simple rectangular SV suffices to explain all of the results.

Figure 8 shows fits for two measurements on the same channel and for the same flow rate. The velocity field of the fluid within the channel is therefore the same for both measurements. The velocity profiles measured by PIV differ, however, because, between the measurement in Figure 8(\subreffig:751) and the measurement in Figure 8(\subreffig:752), the pump was briefly disconnected from and then reconnected to the microchannel, and the microscope was refocused. All processing options were consistent between the two experiments. Though quality of focus is largely subjective, averaging over the SV can explain both measurements. In this case, refocusing changes the position of the SV but not its thickness.

Figure 8: Fits to the experimental profiles obtained with fluid injection rate of 75 l/h. (\subreffig:751) Before disconnecting from the pump and (\subreffig:752) after reconnecting to the pump and refocusing.

In Figure 9 and Table 1 one can see the effects on the SV of differing processing options and flow rates.

Figure 9: Flow profile obtained after processing the same set of CCD images with (\subreffig:100norm) Insight 3G and (\subreffig:100Open) OpenPIV software. The processing options are detailed in Section 3.3. Figure 9(\subreffig:100norm) is a reproduction of the fit shown in Figure 3. We have included it again here for easier comparison.

Figure 9 shows two measured results for the same experiment, in which the fluid injection rate was 100 l/h. Both measured velocity profiles were obtained from the same set of CCD images. They differ only in the image processing options, which are discussed in Section 3.3. The positions of the SV center, , and the SV thickness, , resulting from the various processing options are listed in Table 1 for two fluid injection rates: 100 l/h and 75 l/h. For both flow rates, processing with OpenPIV results in a substantially thinner SV that is located slightly farther from the center of the channel. The preferred processing options may vary from case to case, and it is important to note that their effects on the SV can be quantified.

Insight 3G OpenPIV
100 l/h m m
m m
75 l/h m m
m m
Table 1: Comparison of the SV for two processing options and flow rates.

As discussed in Section 2.2, the DOC derived in [20] is often used as an estimate for the thickness of the SV, . We calculate the DOC for our experiments to be approximately m. This value is based on an assumed intensity threshold, , required for a particle image to affect significantly the processing. We have used the typical [15] value of , which means we assume that contributions are made by out-of-focus particles with a minimum of 1% of the intensity of an in-focus particle. Note that this calculated value of the DOC is constant for all of our experiments presented in Sections 4 and 6, for which we have measured that the actual varies between m and m. The calculated DOC can be reconciled with the measured if one assumes that in fact varies between and , but it is not clear how one could determine these thresholds without first measuring the SV for the various focus levels and processing options. While the DOC gives important intuition into the variables affecting , it does not replace a quantitatively determined SV.

6.4 Analysis of Scanning PIV

Kloosterman et al. [12] have recently claimed that a simple average, such as Equation (1), cannot explain the velocity measured by PIV with a large depth of field. Rather than treat the entire profile from wall to wall, they have considered the difference between measured and predicted maximum speeds at the center of the microfluidic channel. In our Figure 3, for example, this corresponds to explaining why the maximum measured speed is approximately 4.4 mm/s (in symbols) rather than the theoretical maximum for a centered SV with zero thickness, 6.4 mm/s (dashed line).

Their claim is based on a common experimental procedure called “Scanning PIV”, in which one systematically varies the location of the focal plane along the optical () axis and measures the velocity profile for each location. Their analysis is for Hagen-Poiseuille flow in a cylindrical channel, and because the maximum velocity in the channel occurs at its vertical center, the authors claim that from all the measured profiles from the scan, the one with the highest peak velocity is the centered profile. Having centered the region over which particle images are resolved, then, they find that for large depths of field that cover the entire channel, the measured velocity is larger than what is predicted from averaging over the entire channel.

Using the average over the SV in Equation (1), however, we can show that this use of Scanning PIV will not find the center of the channel for large SV. To analyze the results published in Ref. [12], we use the Hagen-Poiseuille profile for flow in a cylindrical channel,


where is the maximum speed in the channel, is the radial coordinate, is the channel radius, and at the center of the channel. Let us change to Cartesian coordinates with at the vertical center of the channel, and substitute for the channel’s diameter. In this section we are interested in the maximum speed, which occurs at the fixed value of the center of the channel, , while the SV changes. We will therefore fix and emphasize the dependence of the speed on the position, , and size, , of the SV. Let denote our prediction of the maximum speed. Using Equation (1) with the Poiseuille profile (7), we find


where the dependence on and is in the integration bounds and , as before.

In terms of the notation in Equation (8), when one uses scanning PIV to find the center of a channel for a given , one varies until is maximized. The principal assumption, then, is that this maximum always occurs at . To check this assumption, we solve for and find multiple solutions. One solution is for a centered SV, with , for which the maximum speed is


Another is for , for which the SV is not centered in the channel, and the maximum speed is


Figure 10 is a plot of for these two maximizing values of . The curves cross at . Because the experimenter selects the value of that maximizes the measured speed, for in a cylindrical channel, this scanning PIV method will position the SV away from the channel’s center, at . Thus, one must use Scanning PIV with caution.

Figure 10: Maximum speed measured by the scanning PIV procedure as a function of . The dashed line is for a DOC centered at , while the solid line represents the maximum speed when the DOC is centered at .

Finally, we note that in Ref. [12] the authors base their claim on a method for calculating the simple average that is different from Equation (1). Because of the assumed ability of scanning PIV to locate the center of the cylindrical channel, they fix . Furthermore, they do not stop the integration at the bottom and top of the channel, and instead use and , where is the Depth of Correlation we have discussed in Section 6.3. The result is


which has positive values only for . Though the does not equal , the two can be roughly compared when both are as large as or larger than the channel diameter, , and images cover the entire depth of the channel. In this case, the authors predict from Equation (11) that , but they measure and find that it remains nearly constant for . We see from Figure 10 and Equation (10), however, that the simple average, combined with this use of Scanning PIV, does predict for large .

7 Conclusion

We show that a simple average over the Sampling Volume (SV) suffices to explain quantitatively the velocity profiles measured in PIV experiments for which the depth of field has a size comparable to that of the microfluidics system itself. This allows rigorous verification of the velocity fields of fluids flowing in microchannels. With this simple approach, one can easily understand and physically interpret the measured features of velocity profiles. In some cases, it allows the complete determination of an unknown flow rate or velocity field based only on the PIV measurement. Furthermore, in a straightforward manner one can determine measurement uncertainty, which is essential for any quantitatively useful experiment. Because of the presence of a finite sampling region, in the past one has either ignored the magnitude of the measured velocity or else treated the effects of the finite volume as errors to be corrected by complicated processing.

By minimally processing data and retaining all measured features, such as kinks in the velocity profile, we have developed a method that is especially robust against noise and other common sources of measurement error. By comparing the measured velocity profiles to the theoretical predictions, we are able to quantify the effect on the SV of processing options and focus. We also provide a critical evaluation of Scanning PIV, which is often used to locate the center of a microchannel. We show how and why it in fact fails to locate the center for large SV. In general, the concept of the SV has implications for quantitative measurement of flow within channels on the micrometer scale and below.


We would like to thank Michael Engel from IBM Research – Watson as well as the LABNANO–CBPF for the SEM images, Diney Ether and the LPO–UFRJ for help with the calibration, Angelo Gobbi at the LMF–LNNano for profilometer measurements, and José Florián at PUC–Rio for help with the PIV equipment.


  • [1] Stuart J Williams, Choongbae Park, and Steven T Wereley. Advances and applications on microfluidic velocimetry techniques. Microfluid Nanofluid, 8(6):709–726, 2010.
  • [2] Ralph Lindken, Massimiliano Rossi, Sebastian Große, and Jerry Westerweel. Micro-particle image velocimetry (PIV): Recent developments, applications, and guidelines. Lab Chip, 9(17):2551–2567, 2009.
  • [3] Steven T Wereley and Carl D Meinhart. Recent advances in micro-particle image velocimetry. Annu. Rev. Fluid Mech., 42:557–576, 2010.
  • [4] Ronald J Adrian and Jerry Westerweel. Particle image velocimetry, volume 30. Cambridge University Press, 2010.
  • [5] Eric K Sackmann, Anna L Fulton, and David J Beebe. The present and future role of microfluidics in biomedical research. Nature, 507(7491):181–189, 2014.
  • [6] Daniel Mark, Stefan Haeberle, Günter Roth, Felix von Stetten, and Roland Zengerle. Microfluidic lab-on-a-chip platforms: requirements, characteristics and applications. Chem. Soc. Rev., 39(3):1153–1182, 2010.
  • [7] Naga Siva Kumar Gunda, Bijoyendra Bera, Nikolaos K Karadimitriou, Sushanta K Mitra, and S Majid Hassanizadeh. Reservoir-on-a-Chip (ROC): A new paradigm in reservoir engineering. Lab Chip, 11(22):3785–3792, 2011.
  • [8] Pierre Joseph and Patrick Tabeling. Direct measurement of the apparent slip length. Phys. Rev. E, 71(3):035303, 2005.
  • [9] Peichun Tsai, Alisia M Peters, Christophe Pirat, Matthias Wessling, Rob GH Lammertink, and Detlef Lohse. Quantifying effective slip length over micropatterned hydrophobic surfaces. Phys. Fluids, 21(11):112002, 2009.
  • [10] Carl D Meinhart, Steve T Wereley, and Juan G Santiago. PIV measurements of a microchannel flow. Exp. Fluids, 27(5):414–419, 1999.
  • [11] Kristin Sott, Tobias Gebäck, Maria Pihl, Niklas Lorén, Anne-Marie Hermansson, Alexei Heintz, and Anders Rasmuson. PIV methodology using model systems for flow studies in heterogeneous biopolymer gel microstructures. J. Colloid Interf. Sci., 398:262–269, 2013.
  • [12] A Kloosterman, C Poelma, and J Westerweel. Flow rate estimation in large depth-of-field micro-PIV. Exp. Fluids, 50(6):1587–1599, 2011.
  • [13] Massimiliano Rossi, Rodrigo Segura, Christian Cierpka, and Christian J Kähler. On the effect of particle image intensity and image preprocessing on the depth of correlation in micro-PIV. Exp. Fluids, 52(4):1063–1075, 2012.
  • [14] J Westerweel, PF Geelhoed, and R Lindken. Single-pixel resolution ensemble correlation for micro-PIV applications. Exp. Fluids, 37(3):375–384, 2004.
  • [15] C Cierpka and CJ Kähler. Particle imaging techniques for volumetric three-component (3d3c) velocity measurements in microfluidics. J. Vis., 15(1):1–31, 2012.
  • [16] Boris Chayer, Katie L Pitts, Guy Cloutier, and Marianne Fenech. Velocity measurement accuracy in optical microhemodynamics: experiment and simulation. Physiol. Meas., 33(10):1585, 2012.
  • [17] J Westerweel. Fundamentals of digital particle image velocimetry. Meas. Sci. Technol., 8(12):1379, 1997.
  • [18] Lukas Novotny and Bert Hecht. Principles of nano-optics. Cambridge university press, 2012.
  • [19] CD Meinhart, ST Wereley, and MHB Gray. Volume illumination for two-dimensional particle image velocimetry. Meas. Sci. Technol., 11(6):809, 2000.
  • [20] MG Olsen and RJ Adrian. Out-of-focus effects on particle image visibility and correlation in microscopic particle image velocimetry. Exp. Fluids, 29(1):S166–S174, 2000.
  • [21] CD Meinhart, AK Prasad, and RJ Adrian. A parallel digital processor system for particle image velocimetry. Meas. Sci. Technol., 4(5):619, 1993.
  • [22] José Angel Florián Gutiérrez and Marcio S Carvalho. Flow of oil drops through micro-capillaries. In 22nd International Congress of Mechanical Engineering (COBEM 2013), pages 6458–6469, 2013.
  • [23] Zachary J Taylor, Roi Gurka, Gregory A Kopp, and Alex Liberzon. Long-duration time-resolved PIV to study unsteady aerodynamics. Instrumentation and Measurement, IEEE Transactions on, 59(12):3262–3269, 2010.
  • [24] Shiyi Chen and Gary D Doolen. Lattice boltzmann method for fluid flows. Annu. Rev. Fluid Mech., 30(1):329–364, 1998.
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description