3D Flow reconstruction using ultrasound PIV

Ultrasound particle image velocimetry (PIV) can be used to obtain velocity fields in non-transparent geometries and/or fluids. In the current study, we use this technique to document the flow in a curved tube, using ultrasound contrast bubbles as flow tracer particles. The performance of the technique is first tested in a straight tube, with both steady laminar and pulsatile flows. Both experiments confirm that the technique is capable of reliable measurements. A number of adaptations are introduced that improve the accuracy and applicability of ultrasound PIV. Firstly, due to the method of ultrasound image acquisition, a correction is required for the estimation of velocities from tracer displacements. This correction accounts for the fact that columns in the image are recorded at slightly different instances. The second improvement uses a slice-by-slice scanning approach to obtain three-dimensional velocity data. This approach is here demonstrated in a strongly curved tube. The resulting flow profiles and wall shear stress distribution shows a distinct asymmetry. To meaningfully interpret these three-dimensional results, knowledge of the measurement thickness is required. Our third contribution is a method to determine this quantity, using the correlation peak heights. The latter method can also provide the third (out-of-plane) component if the measurement thickness is known, so that all three velocity components are available using a single probe.     

3D holographic PIV with a forward-scattering laser sheet and stereoscopic analysis

 A new 3D PIV system combining holography and stereoscopic PIV is presented. The double pulsed holographic recording relies on the forward scattering of particles in the laser sheet. The holographic images of the particles are used for a stereoscopic PIV analysis. An imaging system with a rightangle prism is used to acquire a stereoscopic pair of images. The application of the system to the vortex flow from an inclined delta wing shows the prospects and limitations of the technique. Received: 23 December 1996/Accepted: 1 June 1997     

Error quantification of 3D homogeneous and isotropic turbulence measurements using 2D PIV

Particle image velocimetry (PIV) has become a popular non-intrusive tool for measuring various types of flows. However, when measuring three dimensional flows with 2D PIV, there is inherent measurement error due to out-of-plane motion. Errors in the measured velocity field propagate to turbulence statistics. Since this can distort the overall flow characteristics, it is important to understand the effect of this out-of-plane error. In this study, the effect of out-of-plane motion on turbulence statistics is quantified. Using forced isotropic turbulence direct numerical simulation (DNS) flow field data provided by the Johns Hopkins turbulence database (JHTDB), synthetic image tests are performed. Turbulence statistics such as turbulence kinetic energy, dissipation rate, Taylor microscale, Kolmogorov scale, and velocity correlations are calculated. Various test cases were simulated while controlling three main parameters which affect the out-of-plane motion: PIV interrogation window size, camera inter-frame time, and laser sheet thickness. The amount of out-of-plane motion was first quantified, and then the error variation according to these parameters was examined. This information can be useful when examining fully three dimensional flows such as homogeneous and isotropic turbulence via 2D PIV.     

Time resolved PIV and flow visualization of 3D sheet cavitation

Time-resolved PIV was applied to study fully developed sheet cavitation on a hydrofoil with a spanwise varying angle of attack. The hydrofoil was designed to have a three-dimensional cavitation pattern closely related to propeller cavitation, studied for its adverse effects as vibration, noise, and erosion production. For the PIV measurements, fluorescent tracer particles were applied in combination with an optical filter, in order to remove the reflections of the laser lightsheet by the cavitation. An adaptive mask was developed to find the interface between the vapor and liquid phase. The velocity at the interface of the cavity was found to be very close to the velocity predicted by a simple streamline model. For a visualization of the global flow dynamics, the laser beam was expanded and used to illuminate the entire hydrofoil and cavitation structure. The time-resolved recordings reveal the growth of the attached cavity and the cloud shedding. Our investigation proves the viability of accurate PIV measurements around developed sheet cavitation. The presented results will further be made available as a benchmark for the validation of numerical simulations of this complicated flow.     

Digital direct analysis of a multiexposed photograph in PIV

In Particle Image Velocimetry the analysis of the multi-exposed photographic plate can be carried out either by means of the Young's fringe analysis or by a direct technique. With this second approach the image of each interrogation area, enlarged by a microscope, is recorded by means of a digital image acquisition system. After some preliminary filtering operations the digitized image can be analyzed with numerical tools. In this paper a review of the digital techniques used in the analysis of these images is carried out. Both twodimensional and onedimensional techniques are analyzed and the problems involved in their use are discussed.     

Volumetric correlation PIV: a new technique for 3D velocity vector field measurement

A method is proposed that allows three-dimensional (3D) two-component measurements to be made by means of particle image velocimetry (PIV) in any volume illuminated over a finite thickness. The method is based on decomposing the cross-correlation function into various contributions at different depths. Because the technique is based on 3D decomposition of the correlation function and not reconstruction of particle images, there is no limit to particle seeding density as experienced by 3D particle tracking algorithms such as defocusing PIV and tomographic PIV. Correlations from different depths are differentiated by the variation in point spread function of the lens used to image the measurement volume over that range of depths. A number of examples are demonstrated by use of synthetic images which simulate micro-PIV (μPIV) experiments. These examples vary from the trivial case of Couette flow (linear variation of one velocity component over depth) to a general case where both velocity components vary by different complex functions over the depth. A final validation—the measurement of a parabolic velocity profile over the depth of a microchannel flow—is presented. The same method could also be applied using a thick light sheet in macro-scale PIV and in a stereo configuration for 3D three-component PIV.     

3D pressure imaging of an aircraft propeller blade-tip flow by phase-locked stereoscopic PIV

The flow field at the tip region of a scaled DHC Beaver aircraft propeller, running at transonic speed, has been investigated by means of a multi-plane stereoscopic particle image velocimetry setup. Velocity fields, phase-locked with the blade rotational motion, are acquired across several planes perpendicular to the blade axis and merged to form a 3D measurement volume. Transonic conditions have been reached at the tip region, with a revolution frequency of 19,800 rpm and a relative free-stream Mach number of 0.73 at the tip. The pressure field and the surface pressure distribution are inferred from the 3D velocity data through integration of the momentum Navier-Stokes equation in differential form, allowing for the simultaneous flow visualization and the aerodynamic loads computation, with respect to a reference frame moving with the blade. The momentum and pressure data are further integrated by means of a contour-approach to yield the aerodynamic sectional force components as well as the blade torsional moment. A steady Reynolds averaged Navier-Stokes numerical simulation of the entire propeller model has been used for comparison to the measurement data.     

Fast and accurate PIV computation using highly parallel iterative correlation maximization

Our contribution deals with fast computation of dense two-component (2C) PIV vector fields using Graphics Processing Units (GPUs). We show that iterative gradient-based cross-correlation optimization is an accurate and efficient alternative to multi-pass processing with FFT-based cross-correlation. Density is meant here from the sampling point of view (we obtain one vector per pixel), since the presented algorithm, folki, naturally performs fast correlation optimization over interrogation windows with maximal overlap. The processing of 5 image pairs (1,376 × 1,040 each) is achieved in less than a second on a NVIDIA Tesla C1060 GPU. Various tests on synthetic and experimental images, including a dataset of the 2nd PIV challenge, show that the accuracy of folki is found comparable to that of state-of-the-art FFT-based commercial softwares, while being 50 times faster.     

Spray analysis of a gasoline direct injector by means of two-phase PIV

The hollow-cone spray of a high-pressure swirl injector for a direct-injection spark-ignition (DISI) engine was investigated inside a pressure vessel by means of particle image velocimetry (PIV). As the interaction between the spray droplets and the ambient air is of particular interest for the mixture preparation process, two-phase PIV techniques were applied. To allow phase discrimination, fluorescent seeding particles were used to trace the gas phase. Because of the periodicity of piston engine injection, a statistical evaluation of ensemble-averaged fields to reduce cycle-to-cycle variations and to provide more general information about the two-phase flow was performed. Besides the general spray/air interaction process the investigation of the spray collapse at elevated ambient pressures was the main focus of the study. Future investigations of transient interaction processes require simultaneous techniques in combination with a high-speed camera to resolve the transient interaction phenomena. Therefore, optical filters that attenuate Mie-scattered light and transmit fluorescent light were used to collect both phases on the same image. Consequently, phase separation techniques were employed for data analysis. A masking and a peak separation technique are described and a comparison between the results of an instantaneous two-phase flow field in the spray cone of a DISI injector is presented in the paper.     

Assessment of tomographic PIV in wall-bounded turbulence using direct numerical simulation data

Simulations of tomographic particle image velocimetry (Tomo-PIV) are performed using direct numerical simulation data of a channel flow at Reynolds number of Re _τ = 934, to investigate the influence of experimental parameters such as camera position, seeding density, interrogation volume size and spatial resolution. The simulations employ camera modelling, a Mie scattering illumination model, lens distortion effects and calibration to realistically model a tomographic experiment. Results are presented for camera position and orientation in three-dimensional space, to obtain an optimal reconstruction quality. Furthermore, a quantitative analysis is performed on the accuracy of first and second order flow statistics, at various voxel sizes normalised using the viscous inner length scale. This enables the result to be used as a general reference for wall-bounded turbulent experiments. In addition, a ratio relating seeding density and the interrogation volume size is proposed to obtain an optimal reference value that remains constant. This can be used to determine the required seeding density concentration for a certain interrogation volume size.     

A 3D direct coupling method for steady ship hydroelastic analysis

Asides from the influence of incoming waves, ships can experience steady motions, such as rigid-body sinkage and trim motions, and flexible-body vertical bending motions, due to a constant forward speed even under calm water conditions. In this paper, a novel approach to analyze steady-ship hydroelasticity, particularly for the steady-ship motions and surrounding steady-wave disturbances, is proposed using a three-dimensional (3D) direct coupling method, based on a higher-order boundary element method (HOBEM) and a higher-order shell finite element method (FEM). Within the linearized framework, a solution method is proposed based on a two-step procedure, using two types of Neumann–Kelvin (NK) linear flow models for the fluid part and a virtual work equilibrium equation for the structural part. The first step is to compute a mean position wave-resistance problem using the modified NK equation, the second step is to solve a perturbed position wave-resistance problem, by employing a classical NK model and a virtual work equation based on the first step’s solution. Detailed mathematical formulation and numerical procedures are described, and a few numerical results are illustrated. These include both rigid and flexible steady-ship motions, Von-Mises stress distributions, and wave-resistance coefficients for Froude numbers ranging from 0.15 to 0.5. Furthermore, the numerical results obtained using the present direct coupling method and a modal-based one are compared.     

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).     

Very high Reynolds number boundary layers over 3D sparse roughness and obstacles: the mean flow

Hot-wire and oil-film interferometry measurements are taken for 3D rough wall boundary layers at very high Reynolds numbers (61,000 < Re θ < 120,000) with low blockage ratios, 10 < δ/H < 135, and high roughness, 100 < H ⁺ < 4,900. The results cover flows over both rough walls and over obstacles and are compared with and provide extension to recent lower Reynolds number results. The validity of the Townsend ‘wall similarity hypothesis’ in relation to consistently increasing 3D roughness is interrogated. In agreement with recent work, Schultz and Flack (J Fluid Mech 580:381–405, 2007) and Castro (J Fluid Mech 585:469–485, 2007) found that, for relatively low roughness, Townsend’s hypothesis holds for the mean velocity field. With increasing roughness, the equilibrium layer diminishes and gradually vanishes. The viscous component of the wall shear stress decreases, while the turbulent component increases as the roughness effects extend across the boundary layer.     

3-D PIV via spatial correlation in a color-coded light-sheet

Coding of the light-sheet in depth with different colors and recording with color-sensitive films or CCD cameras enables three-dimensional correlation analysis to obtain the out-of-plane velocity component in 3-D PIV. In the system used, two overlapping light-sheets of different color are recorded simultaneously and the particles' images of the separate sheets are discernible by color splitting. For only two successive exposures as usual in cross-correlation PIV, the resulting images allow for cross-wise plane-to-plane correlations between the separated sheets. This yields altogether three independent correlations. In addition to the usual procedure to obtain the in-plane component, one can determine the out-of-plane velocity component from the three correlation peak values by an appropriate fit of the correlation profile in depth to determine the maximum location with higher accuracy compared to previous methods. In addition, there is no need of a third exposure at a third moment which increases the accuracy for time-varying flows.     

PIV through moving shocks with refracting curvature

Particle image velocimetry (PIV) is applied to moving millimeter shock waves whose density jump and small radii of curvature make refraction significant. The motion of the shock front is also much larger than the motion of the corresponding mass at the front. A Lagrangian model of particle displacement in response to a moving shock is developed to investigate the relationship between particle displacements and the actual mass velocity behind the shock. Errors in PIV measurements due to light refraction across a curved, moving shock are investigated in terms of both position and velocity errors using a refraction model developed from geometrical optics. The model is experimentally validated and applied to 1-D slices of data extracted from PIV vector fields, and the resulting measurement errors are quantified.     

Endoscopic 2D particle image velocimetry (PIV) flow field measurements in IC engines

The paper describes camera and laser endoscopes designed for particle image velocimetry (PIV) applications like measurements in IC engines or turbomachinery. Endoscopic PIV measurements through 8-mm optical access on an IC engine are presented and compared with the measurements using standard optical access through a window.