Direct measurement of fluid velocity gradients at a wall by PIV image processing with stereo reconstruction

Velocity gradient is typically estimated in Particle Image Velocimetry (PIV) by differentiating a measured velocity field, which amplifies noise in the measured velocities. If gradients near a boundary are sought, such noise is usually greater than in bulk fluid, because of small tracer displacement, uncertainty in the effective positions of velocity vectors, intense deformation of tracer patterns, and laser reflection. We consider here a modified form of the Particle Image Distortion (PID) method todirectly calculate velocity gradients at a fixed wall, and refer it as “PIV/IG” (“Interface Gradiometry”). Results from synthetic 2D PIV images suggest our method achieves higher SNR and accuracy than velocity differentiation. Also, we have developed a procedure to reconstruct three-dimensional velocity gradient at a fixed wall the two non-zero components from PIV/IG data obtained in stereo views; these equations simplify considerably thanks to the no-slip condition. Experimental data from the bottom wall of turbulent open channel flow appear to suffer from a form of pixel locking. As with standard PIV, this underlines the importance of adequate tracer diameter in the images, sufficient seeding density, and of dynamic range of the camera sensor.     

Simultaneous 2D flow velocity and gas temperature measurements using thermographic phosphors

In this paper a new approach for simultaneous 2D velocity and temperature measurements using phosphoric particles is presented. The phosphoric particles respond to the temperature changes in the flow while acting as tracers for velocity mapping. The temperature sensitive particles were seeded into a heated flow and were excited by a pulsed UV laser. The subsequent red shifted emission was detected and analyzed to infer temperature using calibration procedures for lifetime and emission spectra against temperature. The diameter of the temperature sensitive particles, usually in the range of 1–10 μm, makes them useful for velocity measurements using particle image velocimetry (PIV). As such, simultaneous measurement of temperature and flow velocity of a gaseous flow were performed and presented. PACS  42.62.-b; 47.80.Cb; 47.80.Fg     

Deformation of viscoelastic coating in a turbulent flow

The rate and amplitude of compliant coating deformation by turbulent pressure pulsations were calculated. Complex compliance determined by a 2D model has two components: along and across the coating. Dependence of the components of dimensionless compliance on the wavelength — coating thickness ratio was determined for 0.3 < λ/H < 30 and dependence of these components on the ratio of flow velocity to velocity of wave propagation was determined for 0.1 < V/C < 10. Deformation amplitude and rate of surface displacement for the hard compliant coatings which can be used in practice were calculated within the range of 5–55 m/s for the water and air turbulent flow. The effects of the loss tangent and Poisson’s ratio of the coating material were also studied. It is shown that the mean-square displacement of their surface does not exceed the thickness of a viscous sublayer. However, the velocity of surface motion is comparable with velocity pulsations in a boundary layer near a wall. This can be a reason for drag reduction on a compliant wall. The calculated value of ratio between energy absorbed by the wall and energy dissipated within the flow because of drag was 10⁻⁴ for water and 10⁻⁶ for air. This estimate does not confirm the hypothesis explaining drag reduction by energy takeoff from the flow.     

Structure of the channel flow behind a surface-mounted rib under conditions of laminar-turbulent transition

The PIV technique was used to measure the instantaneous vector fields of flow velocity and vorticity behind a thin cross-flow rib installed in a channel with laminar, transient, or turbulent flow. The data were treated statistically to determine the fields of mean longitudinal flow velocity and the correlations of pulsating velocity components 〈u′ν′〉 and 〈u′u′〉. Some features of the flow structure developing under conditions of laminar-turbulent transition behind the rib have been revealed.     

超声粒子图像测速技术及应用

心血管疾病的产生与动脉血流的流动状况密切相关。然而,目前普遍应用的超声多普勒成像技术不能精确测量复杂血流流场信息。本文提出了一种基于超声造影微泡的超声全流场粒子图像测速技术,能够获得多维流速速度信息,且不依赖于声束与速度向量之间的夹角。本文首先着重阐述了超声全流场粒子测速技术的基本原理以及系统组成,并对直管流和旋转流场流体动力学特性进行了实验测试研究,实验结果表明本技术能够测量全流场速度,并可作为表征复杂血流流场的有力手段。      

Analysis of jet entrainment mechanism in the transitional regime by time-resolved PIV

Abstract  

The entrainment mechanism in the near field of daisy-shaped and circular orifice jets have been investigated in the transitional regime using time-resolved 2D PIV measurements. The objective is to improve the knowledge from one previous investigation at initial Reynolds number of 800, based on the construction of a pseudo-time resolved PIV fields using the combination of non time-resolved PIV measurements and time-resolved visualizations (Nastase and Meslem J Vis 11(4):309–318, 2008). As expected in the previous work, the entrainment in the circular jet is correlated to the periodic Kelvin–Helmholtz (K–H) ring passing and the entrainment is produced in the braid region where the streamwise structures develop. In the daisy jet, we found that the entrainment rate is not correlated to the periodic K–H vortex passing. The observed small variation of the entrainment rate amplitude in the daisy jet could be related to the K–H dynamics. However, at the studied low Reynolds number the contribution of the K–H vortex on the daisy jet entrainment seems negligible comparing with the streamwise structures role. Furthermore, the real-time resolved measurements allow an indepth analysis of the role played by the K–H ring in the entrainment of circular jet. It is shown that the entrainment is not only produced in the braid region but is also present in the upstream part of the K–H ring. In the downstream part of the ring, the entrainment is dramatically reduced. This new observation opens a question which still has to be answered with time-resolved 3D PIV measurements. The question is “Whether the depression formed due to the ring passing or the streamwise structures rolled-up on the ring is responsible for entrainment at the upstream part of the ring?”     

Electrohydrodynamic flow in a wire-plate non-thermal plasma reactor measured by 3D PIV method

This work was aimed at measurements of the electrohydrodynamic (EHD) secondary flow in a non-thermal plasma reactor using three-dimensional particle image velocimetry (3D PIV) method. The wide-type non-thermal plasma reactor used in this work was an acrylic box with a wire discharge electrode and two plate collecting electrodes. The positive DC voltage was applied to the wire electrode through a 10 MΩ resistor. The collecting electrodes were grounded. The voltage applied to the wire electrode was 28 kV. Air flow seeded with a cigarette smoke was blown along the reactor duct with an average velocity of 0.6 m/s. The 3D PIV velocity fields measurements were carried out in four parallel planes stretched along the reactor duct, perpendicularly to the wire electrode and plate electrodes. The measured flow velocity fields illustrate complex nature of the EHD induced secondary flow in the non-thermal plasma reactor.     

对称槽道涡波流场流动特征的POD分析

为深入分析层流状态下对称槽道内涡波流场的流动特性及其变化规律,对流场进行了二维粒子图像测速(2DPIV)测量获取瞬态速度矢量数据,利用本征正交分解(POD)技术进行模态分解以及涡波流场的重构,然后根据重构的流场对对称槽道内涡波流场进行了平均速度剖面、流场脉动强度以及特征点的速度和频谱分布等方面的分析。结果表明:POD的前15阶模态能够表征涡波流场的主导结构,第1,3阶模态主要表现为一对旋向相反的涡对特征,第2阶模态具有涡旋和波状主流的特征;提取了5个涡旋涡核的位置作为流场流动特性的特征点;根据POD重构流场分析发现流向平均速度呈抛物线形状分布,法向平均速度呈对称分布特征;流向脉动强度受壁面的影响较大,法向脉动强度呈现抛物线形状分布;距离中心主流较近的1#,4#,5#特征点的速度脉动程度受主流的脉动强度影响较大,速度的脉动主频0.15 Hz与次频、流场的自然频率0.35 Hz共同影响特征点的速度分布;2#,3#特征点的流向速度呈衰减趋势,法向速度在初期幅度变化较大。     

Rapid imaging of fluid flow patterns in a narrow packed bed using MRI

A gradient echo rapid velocity and acceleration imaging sequence (GERVAIS) has been developed and implemented to image liquid flow within a narrow packed bed. Two-dimensional velocity images have been acquired with an in-plane pixel size of 781 microm x 781 microm, with a data acquisition time of 20 ms for a single velocity component. Images of the x, y and z velocity vectors are reported. Data are reported for Reynolds numbers (based on particle diameter) of 200 and 300. In each case, GERVAIS images are compared with the results of a standard spin-echo phase-encoding velocity measurement. At Re = 200, steady-state flow is expected and the velocity images acquired using both techniques are consistent. At Re = 300, the GERVAIS sequence is able to image the unsteady-state flow field within this system. In contrast, the standard phase-encoding velocity measurement contains significant artefacts.     

A comparison of 2D and 3D PIV measurements in an oblique jet

A stereo PIV (SPIV) acquisition and analysis system was developed to measure three velocity components in planar flow fields. The analysis software is based on a third order mapping function method. The system was calibrated by imaging a square grid in three measurement planes with two Kodak Megaplus cameras oriented at 30 to the bisector between them. The camera images were dewarped into real coordinates by employing a set of transform matrices computed for each calibration plane. Bias and rms errors were determined by comparing displacements measured directly with displacements estimated from the dewarping and recombination algorithm. The bias errors in the directions parallel with the measurement plane were negligible while the bias in thez direction was about 0.6 pixel. The rms errors, 0.2–0.3 pixels, were largest in thez direction. These errors were thought to result from limitations in the calibration method. The SPIV system was tested in a two-dimensional oblique jet with Reynolds number of 1800. The three dimensional results were taken in a vertical (x, y) plane parallel with the jet span. The SPIV results were compared with LDV data and two-dimensional PIV data obtained in a vertical (y, z) plane of the same jet. The SPIV measurements yielded accurate values for the in-plane mean and rms velocity components. The measured out-of-plane mean component was underestimated due to the bias error mentioned above. The rms component was accurate in part of the field but overestimated in another part due to local variations in rms error. It is expected that in the future, the out-of-plane errors can be minimized by improving the calibration and transformation procedures.     

Visualization of low Reynolds boundary-driven cavity flows in thin liquid shells

Abstract  

Classic examples of low-Reynolds recirculating cavity flows are typically generated from lid-driven boundary motion at a solid–fluid interface, or alternatively may result from shear flow over cavity openings. Here, we are interested in an original family of boundary-driven cavity flows occurring, in contrast to classic setups, at fluid–fluid interfaces. Particle image velocimetry (PIV) is used to investigate the structure of internal convective flows observed in thin liquid shells. Under the specific configuration investigated, the soap bubble’s liquid shell is in fact in motion and exhibits sporadic local “bursts”. These bursts induce transient flow motion within the cavity of order Re ∼ O(1). The combination of PIV and proper orthogonal decomposition (POD) is used to extract dominant flow structures present within bubble cavities. Next, we show that thermally induced Marangoni flows in the liquid shell can lead to forced, (quasi) steady-state, internal recirculating flows. The present findings illustrate a novel example of low-Reynolds boundary-driven cavity flows.     

3D imaging of evaporating fuel droplets by stereoscopic PIV

A gun-type burner is a widely used oil burner for industrial and domestic applications. The oil is pressure-atomized and mixed with air generating a recirculating, swirling flow. Because of the surrounding flame, fuel droplets evaporate, being difficult to obtain information on droplets’ dynamics. Several laser techniques have been applied to this burner for spray diagnosis. PDA provides information about droplet size and velocity but can say little about the instantaneous spatial structures in the flow. Planar laser techniques as PIV can describe the 2D instantaneous spatial structures, but cannot provide information about the 3D structures in the flow. Then Stereoscopic PIV was applied. This technique allows us to measure the full 3D velocity vector map in a whole fluid plane. This paper has a double purpose. Firstly, to visualize the 3D structures which are present in the burner; secondly, to show that Stereoscopic PIV is an applicable technique for the diagnosis of an evaporating spray.     

Imaging of subcutaneous blood vessels and flow velocity profiles by optical coherence tomography

We have applied a compact low power rapid scanning Doppler Optical Coherence Tomography system to monitor multi-dimensional velocity profiles within the complex vessels and simultaneous real-time non-invasive imaging of skin tissues morphology in vivo, in the wavelength range of 1.3–1.5 nm. Optical clearing of skin tissues has been utilized to achieve depth of OCT images up to 1.7 mm. Current approach enables applying low-power (0.4–0.5 mW) and low-noise broadband near-infrared light sources and obtaining OCT images with down to 12 μm spatial resolution. Two-dimensional time-domain OCT images of complex flow velocity profiles in blood vessel phantom and in vivo subcutaneous human skin tissues are presented. The effect of optical clearing on in vivo images is demonstrated and discussed.     

A comparative study of the PIV and LDV measurements on a self-induced sloshing flow

Particle Imaging Velocimetry (PIV) and Laser Doppler Velocimetry (LDV) measurements on a self-induced sloshing flow in a rectangular tank had been conducted in the present study. The PIV measurement result was compared with LDV measurement result quantitatively in order to evaluate the accuracy level of the PIV measurement. The comparison results show that the PIV and LDV measurement results agree with each other well in general for both mean velocity and fluctuations of the velocity components. The average disagreement level of the mean velocity between PIV and LDV measurement results was found to be within 3% of the target velocity for the PIV system parameter selection. Bigger disagreements between the PIV and LDV measurement results were found to concentrate at high shear regions. The spatial resolution and temporal resolution differences of the PIV and LDV measurements and the limited frames of the PIV instantaneous results were suggested to be the main reasons for the disagreement.     

Schlieren “PIV” for turbulent flows

The possibility of using commercial PIV equipment combined with schlieren optics to measure the velocity fields of turbulent flows is explored. Given a sufficiently high Reynolds number and adequate refractive flow differences, turbulent eddies can serve as the PIV “particles” in a schlieren image or shadowgram. The PIV software analyzes motion between consecutive schlieren or shadowgraph frames to obtain velocity fields. Velocimetry examples of an axisymmetric sonic helium jet in air and a 2D turbulent boundary layer at Mach 3 are shown. Due to optical path integration, axisymmetric flows require the inverse Abel transform to extract center-plane velocity data. Conditions for optimum schlieren sensitivity are examined. In its present embodiment, “schlieren PIV” is not useful for laminar flows nor for fully 3D flows. Otherwise it functions much like standard PIV under conditions where individual particles are not resolved and velocimetry is instead based on correlation of the motion of turbulent structures. “Schlieren PIV” shows significant promise for general refractive turbulent flow velocimetry if its integrative nature can be overcome through sharp-focusing optics.     

Stereoscopic PIV measurements of flow around a marine propeller

A stereoscopic PIV (Particle Image Velocimetry) technique has been employed to measure the 3 dimensional flow structure of turbulent wake behind a marine propeller with 5 blades. The out-of-plane velocity component was measured using particle images captured simultaneously by two CCD cameras installed in the angular displacement configuration. 400 instantaneous velocity fields were acquired for each of four different blade phases of 0°, 18°, 36° and 54°. They were ensemble averaged to investigate the spatial evolution of propeller wake in the near wake region up to one propeller diameter (D) downstream. The phase-averaged velocity fields show clearly the viscous wake formed by the boundary layers developed along both surfaces of the blade. Tip vortices were generated periodically and the slipstream contraction occurs in the near-wake region. The out-of-plane velocity component has large values at the locations of tip and trailing vortices. With going downstream, the axial turbulence intensity and the strength of tip vortices were decreased due to the viscous dissipation, turbulence diffusion and blade-to-blade interaction. The difference in the mean velocity fields measured by SPIV and 2-D PIV methods was about 5% ≈ 10%. However, the 2-D PIV results also give sufficient information on propeller wake beyond the region of X/D=0.2.     

A generalised optimisation method for double pulsed particle image velocimetry

A new method is presented which allows for the optimisation of a double pulsed particle image velocimetry (PIV) system given a priori information on the dynamic range and velocity gradients at a chosen local region in the flow. The method is developed from considerations of the worst case signal-to-noise characteristics corresponding to the autocorrelation function obtained from this region. Optimisation permits the user to define a window of operation for a PIV system which, providing a minimum level of seeding is present, will ensure that from a series of experiments on average a minimum of 50% valid vectors will be obtained from the chosen region. Validation of the method was carried out by using a Monte Carlo simulation based on the theoretical model. Optimisation and Monte Carlo prediction were shown to agree to within 1.5%.     

Specific features of waves in magnetic crystals

On the basis of the asymmetric velocity tensor introduced we develop an orthogonal representation of wave propagation in anisotropic magnetic dielectrics. Simple expressions are obtained for the direct dependence of polarization vectors, refractive indices, and phase and group velocities of isonormal waves in a magnetic dielectric on the vectors of principal velocities. We show that the presence of magnetic properties is manifested in the optical range in the fact that the velocities of wave propagation in uniaxial crystals are determined by three main refractive indices, and therefore, both waves except for specific directions become extraordinary, whereas in biaxial crystals (also because of a large number of parameters) the cones of internal and external refraction become elliptic. Translated from Zhurnal Prikladnoi Spektroskopii, Vol. 66, No. 5, pp. 597–600, September–October, 1999     

Quantitative visualization of micro-tube flow using micro-PIV

This paper describes PIV measurements ofthe flow field in a micro round tube with an internal diameter of 100 μm in order to examine micro-scale effects. Since the refractive index of the micro tube almost corresponds to that of water, the inner flow in the tube can be observed clearly. The micro PIV system has been developed using a microscope, a high sensitive CCD camera, a double pulsed Nd:YAG laser and optics. Applying the micro PIV technique to the flow, the velocity distributions with spatial resolution of 1.8 × 1.8 μm were measured even near the wall in the center plane of the round tube. It was found that the velocities near the tube wall were smaller than the theoretical values calculated by using Poiseuille’s law. It is believed that this disparity is due to micro-scale effects such as interference between particles and the wall, friction at the wall, surface tension and so on.     

Application of PIV to over-expanded supersonic flows: Possibilities and limits

Two-dimensional velocity distributions outside a Mach 2.0 supersonic nozzle have been investigated using a digital particle im age velocimetry (PIV). Mean velocities , vor ticity field and volume dilatation field were obtained from PIV images using 0 .33 μm titanium dioxide (TiO₂) particle. The seeding particle of larger size , 1.4 μrn Ti0₂, was also used for the experimental comparison of velocity lag downstream of shock waves. The results have been compared and analyzed with schlieren photographs for the locations of shock waves and over-expanded shock structure to inspect possibilities and limits of a PIV technique to over-expanded supersonic flows. It is found that although the quantitative velocity measurement using PIV on over-expanded supersonic flows with large velocity and pressure gradients is limited, the locations of normal shock and oblique shock waves can be resolved by the axial/radial velocity fields, and over-expanded shock structure can be predicted by vorticity field and volume dilatation field which are acquired from the spatial differential of the velocity field.