Three-component micro-PIV using the continuity equation and a comparison of the performance with that of stereoscopic measurements

Recently, a number of techniques have been presented for the determination of the third “out-of-plane” velocity component in micro particle image velocimetry (micro-PIV) data. In particular, the conventional macroscopic stereo-PIV technique has been converted to the micro scale by the use of stereo-microscopy. In this work a different technique is investigated, which uses conventional, two-component micro-PIV to generate velocity data on a number of planes. The in-plane velocity gradients are then calculated, which can be used in the continuity equation to produce the out-of-plane velocity gradients. These, together with the no-penetration boundary condition, can then be used to calculate the out-of-plane velocities. An algorithm is presented that is capable of handling up to one invalid vector per column of data by using a combination of first order and second order projections of the velocity. The advantage of the continuity based technique is that it uses the existing two-component micro-PIV technology, which at present is in a more advanced stage of development then stereo-microscopy based micro-PIV. The technique is investigated using a flow similar to one used previously to assess stereoscopic micro-PIV (Meas Sci Technol 17:2175–2185, 2006). This allows a comparison of the performance of the two techniques. The results show that the continuity based data agrees well with an independent computational fluid dynamics solution and has a smaller experimental uncertainty than the stereoscopic technique at a better spatial resolution. There are, however, potential limitations to the continuity based technique. These include the two-dimensionality of the data, which is limited to the planes on which the original images were taken, and the dependence of the technique on the data close to surfaces, where the experimental errors are often greatest. Stereoscopic micro-PIV does not have these limitations so, whilst at present it appears that continuity based techniques may be more accurate, there is sufficient potential for stereoscopic techniques to justify their further development.     

Turbulent thermal diffusion in a multi-fan turbulence generator with imposed mean temperature gradient

We studied experimentally the effect of turbulent thermal diffusion in a multi-fan turbulence generator which produces a nearly homogeneous and isotropic flow with a small mean velocity. Using particle image velocimetry and image processing techniques, we showed that in a turbulent flow with an imposed mean vertical temperature gradient (stably stratified flow) particles accumulate in the regions with the mean temperature minimum. These experiments detected the effect of turbulent thermal diffusion in a multi-fan turbulence generator for relatively high Reynolds numbers. The experimental results are in compliance with the results of the previous experimental studies of turbulent thermal diffusion in oscillating grid turbulence (Buchholz et al. 2004; Eidelman et al. 2004). We demonstrated that the turbulent thermal diffusion is an universal phenomenon. It occurs independently of the method of turbulence generation, and the qualitative behavior of particle spatial distribution in these very different turbulent flows is similar. Competition between turbulent fluxes caused by turbulent thermal diffusion and turbulent diffusion determines the formation of particle inhomogeneities.     

Velocity fields and streamline patterns of miscible displacements in cylindrical tubes

Displacements of a viscous fluid by a miscible fluid of a lesser viscosity and density in cylindrical tubes were investigated experimentally. Details of velocity and Stokes streamline fields in vertical tubes were measured using a DPIV (digital particle image velocimetry) technique. In a reference frame moving with the fingertip, the streamline patterns around the fingertip obtained from the present measurements confirm the hypothesis of Taylor (1961) for the external patterns, and that of Petitjeans and Maxworthy (1996) for the internal patterns. As discussed in these papers, the dependent variable, m, a measure of the volume of viscous fluid left on the tube wall after the passage of the displacing finger, is a parameter that determines the flow pattern. When m>0.5 there is one stagnation point at the tip of the finger; when m<0.5 there are two stagnation points on the centerline, one at the tip and the other inside the fingertip, and a stagnation ring on the finger surface with a toroidal recirculation in the fingertip between the two stagnation points. The finger profile is obtained from the zero streamline of the streamline pattern.An erratum to this article can be found at     

Experimental Investigation on Turbulent Structure of Drag Reducing Channel Flow with Blowing Polymer Solution from the Wall

Turbulent drag reducing flow with blowing polymer solution from the channel wall was investigated experimentally using particle image velocimetry (PIV). Experiments were carried out with varying conditions of blowing polymer solution (e.g. weight concentration of polymer solution). Reynolds number based on the channel height and mean velocity was set to 20000 and 40000. When the polymer solution was blown from the channel wall, streamwise velocity fluctuation little increased, but wall-normal velocity fluctuation, Reynolds shear stress and correlation coefficient decreased significantly only near the blower wall. This behavior corresponds to the decrease of the ejection and sweep in the near-wall region observed by the investigation of instantaneous velocity map. On the contrary, this characteristic behavior was not observed at a position away from the blower wall (y/(H/2) > 0.4) and the scatter plot was almost the same as that of the water flow in this region. These results suggest that there are two regions in the drag reducing flow with blowing polymer solution from the wall; one is a non-Newtonian region which exists near the blower wall, and the other is a Newtonian region at a distance from the wall. The non-Newtonian region plays a key role in the drag reduction by the blowing polymer solution.     

Reliable treatment of a new analytical method for solving MHD boundary-layer equations

The purpose of this study is to implement a new analytical method which is a combination of the homotopy analysis method (HAM) and the Padé approximant for solving magnetohydrodynamic boundary-layer flow. The solution is compared with the numerical solution. Comparisons between the HAM–Padé and the numerical solution reveal that the new technique is a promising tool for solving MHD boundary-layer equations. The effects of the various parameters on the velocity and temperature profiles are presented graphically form. Favorable comparisons with previously published works (Crane, J. Appl. Math. Phys. 21:645–647, 1970, and Vajravelu and Hadjinicolaou, Int. J. Eng. Sci. 35:1237–1244, 1997) are obtained. It is predicted that HAM–Padé can have wide application in engineering problems (especially for boundary-layer and natural convection problems).     

Use of particle image velocimetry to study heterogeneous drag reduction

Large polymer filaments can form when drag reducing polymers are injected through wall slots. The presence of these structures enhances the performance of the drag reducing function by mechanisms which are not understood. This paper shows how particle image velocimetry (PIV) techniques can be used to study changes in the configuration of the injected polymer and in the structure of the velocity field with increasing drag reduction. The filaments are found to behave as solid bodies which break up in high shear regions close to a boundary. The breakup process provides an explanation of why the filaments are not observed close to a wall and offers the possibility of providing a heterogeneous distribution of small aggregates of polymers which could be more effective than uniformly distributed molecules as suggested by Hoyer and Gyr (J Non-Newton Fluid Mech 65:221–240, 1996; J Fluids Eng 120:818–823, 1998), Dunlop and Cox (Phys Fluids 20:203–213, 1977) and Vlachogiannis et al. (Phys Fluid 15:3786–3794, 2004). PIV measurements show dramatic qualitative changes in the velocity patterns at maximum drag reduction.     

Slow motions of a circular cylinder experiencing slip near a plane wall

A theoretical study is presented for the two-dimensional creeping flow caused by a long circular cylindrical particle translating and rotating in a viscous fluid near a large plane wall parallel to its axis. The fluid is allowed to slip at the surface of the particle. The Stokes equations for the fluid velocity field are solved in the quasi-steady limit using cylindrical bipolar coordinates. Semi-analytical solutions for the drag force and torque acting on the particle by the fluid are obtained for various values of the slip coefficient associated with the particle surface and of the relative separation distance between the particle and the wall. The results indicate that the translation and rotation of the confined cylinder are not coupled with each other. For the motion of a no-slip cylinder near a plane wall, our hydrodynamic drag force and torque results reduce to the closed-form solutions available in the literature. The boundary-corrected drag force and torque acting on the particle decrease with an increase in the slip coefficient for an otherwise specified condition. The plane wall exerts the greatest drag on the particle when its migration occurs normal to it, and the least in the case of motion parallel to it. The enhancement in the hydrodynamic drag force and torque on a translating and rotating particle caused by a nearby plane wall is much more significant for a cylinder than for a sphere.     

Multilayer nano-particle image velocimetry

Nano-particle image velocimetry (nPIV), based on evanescent-wave illumination of fluorescent colloidal tracers, measures the two velocity components parallel to the wall averaged over the first few hundred nanometers next to the wall. The intensity of the evanescent wave decays exponentially with z, or the distance normal to the wall. Illuminated tracers closer to the wall therefore have images that are brighter than those farther from the wall. This nonuniform illumination presents the possibility to extend the technique to “multilayer nPIV,” where the two velocity components parallel to the wall can be estimated at different z-locations within the illuminated region. In this paper, the variation of tracer image intensity with distance from the wall was predicted using diffraction optics-based approaches. The predictions, which were validated by calibration experiments, show that particle image intensity decays exponentially with distance normal to the wall. The feasibility of multilayer nPIV was evaluated using artificial images of plane Couette flow that incorporate evanescent-wave illumination, hindered Brownian diffusion and image noise. Each image was divided into three sub-images based on tracer image intensity, and standard techniques were then used to extract temporally and spatially averaged velocities at three different z-locations. In these simulations, velocity data were obtained within 80 nm of the wall, a threefold improvement over previous measurements. The results demonstrate that multilayer nPIV is feasible if appropriate classification techniques are developed and used to separate tracer images into different layers.

---

Application of multivariate outlier detection to fluid velocity measurements

A statistical-based approach to detect outliers in fluid-based velocity measurements is proposed. Outliers are effectively detected from experimental unimodal distributions with the application of an existing multivariate outlier detection algorithm for asymmetric distributions (Hubert and Van der Veeken, J Chemom 22:235–246, 2008). This approach is an extension of previous methods that only apply to symmetric distributions. For fluid velocity measurements, rejection of statistical outliers, meaning erroneous as well as low probability data, via multivariate outlier rejection is compared to a traditional method based on univariate statistics. For particle image velocimetry data, both tests are conducted after application of the current de facto standard spatial filter, the universal outlier detection test (Westerweel and Scarano, Exp Fluids 39:1096–1100, 2005). By doing so, the utility of statistical outlier detection in addition to spatial filters is demonstrated, and further, the differences between multivariate and univariate outlier detection are discussed. Since the proposed technique for outlier detection is an independent process, statistical outlier detection is complementary to spatial outlier detection and can be used as an additional validation tool.     

Full 3D correlation tensor computed from double field stereoscopic PIV in a high Reynolds number turbulent boundary layer

The turbulence structure near a wall is a very active subject of research and a key to the understanding and modeling of this flow. Many researchers have worked on this subject since the fifties Hama et al. (J Appl Phys 28:388–394, 1957). One way to study this organization consists of computing the spatial two-point correlations. Stanislas et al. (C R Acad Sci Paris 327(2b):55–61, 1999) and Kahler (Exp Fluids 36:114–130, 2004) showed that double spatial correlations can be computed from stereoscopic particle image velocimetry (SPIV) fields and can lead to a better understanding of the turbulent flow organization. The limitation is that the correlation is only computed in the PIV plane. The idea of the present paper is to propose a new method based on a specific stereoscopic PIV experiment that allows the computation of the full 3D spatial correlation tensor. The results obtained are validated by comparison with 2D computation from SPIV. They are in very good agreement with the results of Ganapthisubramani et al. (J Fluid Mech 524:57–80, 2005a).     

The two-dimensional velocity shift caused by the use of a rotating mirror in PIV flow field measurements

A theoretical analysis of the two-dimensional velocity shift by the use of a rotating mirror in PIV flow field measurements has been carried out by the application of ray optics. The velocity shift has been calculated directly in the flow field co-ordinates. In dimensionless form the results are also available in the image plane. It is shown that the distribution of the velocity shift over the entire observation field may not usually be assumed to be uniform (and may vary by up to 20%). The effect of the layout and arrangement of the rotating mirror on the velocity shift and its distribution has been analysed. In addition, the sensitivity of the velocity shift to the design tolerances and imperfections of the rotating mirror is given.     

The measurement of square channel velocity profiles using a microcomputer-based image analysis system

A technique was developed to perform automated velocity measurements from a sequence of particle images. A very thin sheet of laser light allows determination of essentially two-dimensional velocity profiles in very small conduits. A four image sequence was captured by a microcomputer-based frame grabber. After thresholding to eliminate particles not perfectly centered in the sheet of light, the sequence of pseudo-colored images of a given particle is used to determine its velocity. By measuring several two-dimensional velocity profiles across a square channel, the complete three-dimensional velocity profile was assembled. The experimentally measured velocity profile agrees closely with the known theoretical velocity profile for flow in a square channel.     

Three-dimensional particle tracking method using FPGA-based real-time image processing and four-view image splitter

We present a cost-effective solution of the three-dimensional particle tracking velocimetry (3D-PTV) system based on the real-time image processing method (Kreizer et al. Exp Fluids 48:105–110, 2010) and a four-view image splitter. The image processing algorithm, based on the intensity threshold and intensity gradients estimated using the fixed-size Sobel kernel, is implemented on the field-programmable gate array integrated into the camera electronics. It enables extracting positions of tracked objects, such as tracers or large particles, in real time. The second major component of this system is a four-view split-screen device that provides four views of the observation volume from different angles. An open-source ray-tracing software package allows for a customized optical setup for the given experimental settings of working distances and camera parameters. The specific design enables tracking in larger observation volumes when compared to the designs published up to date. The present cost-effective solution is complemented with open-source particle tracking software that receives raw data acquired by the real-time image processing system and returns trajectories of the identified particles. The combination of these components simplifies the 3D-PTV technique by reducing the size and increasing recording speed and storage capabilities. The system is capable to track a multitude of particles at high speed and stream the data over the computer network. The system can provide a solution for the remotely controlled tracking experiments, such as in microgravity, underwater or in applications with harsh experimental conditions.     

Propeller tip and hub vortex dynamics in the interaction with a rudder

In the present paper, the interaction mechanisms of the vortices shed by a single-screw propeller with a rudder installed in its wake are addressed; in particular, following the works by Felli et al. (Exp Fluids 6(1):1–11, 2006a, Exp Fluids 46(1):147–1641, 2009a, Proceedings of the 8th international symposium on particle image velocimetry: Piv09, Melbourne, 2009b), the attention is focused on the analysis of the evolution, instability, breakdown and recovering mechanisms of the propeller tip and hub vortices during the interaction with the rudder. To investigate these mechanisms in detail, a wide experimental activity consisting in time-resolved visualizations, velocity measurements by particle image velocimetry (PIV) and laser Doppler velocimetry (LDV) along horizontal chordwise, vertical chordwise and transversal sections of the wake have been performed in the Cavitation Tunnel of the Italian Navy. Collected data allows to investigate the major flow features that distinguish the flow field around a rudder operating in the wake of a propeller, as, for example, the spiral breakdown of the vortex filaments, the rejoining mechanism of the tip vortices behind the rudder and the mechanisms governing the different spanwise misalignment of the vortex filaments in the pressure and suction sides of the appendage.     

Simulation of particles released near the wall in a turbulent boundary layer

A direct numerical simulation was used along with a Lagrangian particle tracking technique to study particle motion in a horizontal, spatially developing turbulent boundary layer along an upper-wall (with terminal velocity directed away from the wall). The objective of the research was to study particle diffusion, dispersion, reflection, and mean velocity in the context of two parametric studies: one investigated the effect of the drift parameter (the ratio of particle terminal velocity to fluid friction velocity) for a fixed and finite particle inertia, and the second varied the drift parameter and particle inertia by the same amount (i.e. for a constant Froude number). A range of drift parameters from 10⁻⁴ to 10⁰ were considered for both cases. The particles were injected into the simulation at a height of four wall units for several evenly distributed points across the span and a perfectly elastic wall collision was specified at one wall unit.Statistics collected along the particle trajectories demonstrated a transition in particle movement from one that is dominated by diffusion to one that is dominated by gravity. For small and intermediate sized particles (i.e. ones with outer Stokes numbers and drift parameters much less than unity) transverse diffusion away from the wall dominated particle motion. However, preferential concentration is seen near the wall for intermediate-sized particles due to inhomogeneous turbulence effects (turbophoresis), consistent with previous channel flow studies. Particle–wall collision statistics indicated that impact velocities tended to increase with increasing terminal velocity for small and moderate inertias, after which initial conditions become important. Finally, high relative velocity fluctuations (compared to terminal velocity) were found as particle inertia increased, and were well described with a quasi-one-dimensional fluctuation model.     

Axisymmetric creeping motion of a slip spherical particle in a nonconcentric spherical cavity

A semianalytical study of the creeping flow caused by a spherical fluid or solid particle with a slip surface translating in a viscous fluid within a spherical cavity along the line connecting their centers is presented in the quasisteady limit of small Reynolds number. In order to solve the Stokes equations for the flow field, a general solution is constructed from the superposition of the fundamental solutions in the two spherical coordinate systems based on both the particle and cavity. The boundary conditions on the particle surface and cavity wall are satisfied by a collocation technique. Numerical results for the hydrodynamic drag force exerted on the particle are obtained with good convergence for various values of the ratio of particle-to-cavity radii, the relative distance between the centers of the particle and cavity, the relative viscosity or slip coefficient of the particle, and the slip coefficient of the cavity wall. In the limits of the motions of a spherical particle in a concentric cavity and near a cavity wall with a small curvature, our drag results are in good agreement with the available solutions in the literature. As expected, the boundary-corrected drag force exerted on the particle for all cases is a monotonic increasing function of the ratio of particle-to-cavity radii, and becomes infinite in the touching limit. For a specified ratio of particle-to-cavity radii, the drag force is minimal when the particle is situated at the cavity center and increases monotonically with its relative distance from the cavity center to infinity in the limit as it is located extremely away from the cavity center. The drag force acting on the particle, in general, increases with an increase in its relative viscosity or with a decrease in its slip coefficient for a given configuration, but surprisingly, there are exceptions when the ratio of particle-to-cavity radii is large.     

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

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

A color-coded backlighted defocusing digital particle image velocimetry system

Defocusing digital particle image velocimetry (DDPIV), as a true three-dimensional (3D) measurement technique, allows for the measurement of 3D velocities within a volume. Initially designed using a single CCD and 3-pinhole mask (Willert and Gharib in Exp Fluids 12:353–358, 1992), it has evolved into a multi-camera system in order to overcome the limitations of image saturation due to multiple exposures of each particle. In order to still use a single camera and overcome this limitation, we have modified the original single CCD implementation by placing different color filters over each pinhole, thus color-coding each pinhole exposure, and using a 3-CCD color camera for image acquisition. Due to the pinhole mask, there exists the problem of a significant lack of illumination in a conventional lighting setup, which we have solved by backlighting the field-of-view and seeding the flow with black particles. This produces images with a white background superimposed with colored triple exposures of each particle. A color space linear transformation is used to allow for accurate identification of each pinhole exposure when the color filters’ spectrum does not match those of the 3-CCD color camera. Because the imaging is performed with a multi-element lens instead of a single-element lens, an effective pinhole separation, d _e, is defined when using a pinhole mask within a multi-element lens. Calibration results of the system with and without fluid are performed and compared, and a correction of the effective pinhole separation, d _e, due to refraction through multiple surfaces is proposed. Uncertainty analyses are also performed, and the technique is successfully applied to a buoyancy-driven flow, where a 3D velocity field is extracted.     

Suggestive correctional methods for PIV image biasing caused by a rotating mirror system

A rotating mirror is widely used to generate the velocity shift that can resolve the directional ambiguities of PIV (particle image velocimetry) measurements. The rotating mirror system inevitably creates the normal displacement of the resulting PIV images and causes systematic image errors. Two corrective methods are proposed to eliminate or reduce the image biasing in PIV system. The use of two linearly traversing mirrors, instead of a single rotating mirror, shows that image biasing can be eliminated and the velocity shift well generated. As a second option, two co-rotating mirrors, instead of one, can reduce the image biasing with a maximum velocity shift available. Detailed imaging kinematics of the two suggestive methods are presented to lead to designing of practical devices that improve the PIV capabilities by reducing the systematic image errors. The authors like to acknowledge the partial financial support provided from the Turbo and Power Machinery Research Center (TPMRC) of Seoul National University, Seoul, Korea.