Tomographic particle image velocimetry

This paper describes the principles of a novel 3D PIV system based on the illumination, recording and reconstruction of tracer particles within a 3D measurement volume. The technique makes use of several simultaneous views of the illuminated particles and their 3D reconstruction as a light intensity distribution by means of optical tomography. The technique is therefore referred to as tomographic particle image velocimetry (tomographic-PIV). The reconstruction is performed with the MART algorithm, yielding a 3D array of light intensity discretized over voxels. The reconstructed tomogram pair is then analyzed by means of 3D cross-correlation with an iterative multigrid volume deformation technique, returning the three-component velocity vector distribution over the measurement volume. The principles and details of the tomographic algorithm are discussed and a parametric study is carried out by means of a computer-simulated tomographic-PIV procedure. The study focuses on the accuracy of the light intensity field reconstruction process. The simulation also identifies the most important parameters governing the experimental method and the tomographic algorithm parameters, showing their effect on the reconstruction accuracy. A computer simulated experiment of a 3D particle motion field describing a vortex ring demonstrates the capability and potential of the proposed system with four cameras. The capability of the technique in real experimental conditions is assessed with the measurement of the turbulent flow in the near wake of a circular cylinder at Reynolds 2,700.     

Effect of particle number density in in-line digital holographic particle velocimetry

A digital in-line holographic particle tracking velocimetry (HPTV) system was developed to measure 3D (three-dimensional) velocity fields of turbulent flows. The digital HPTV (DHPTV) procedure consists of four steps: recording, numerical reconstruction, particle extraction and velocity extraction. In the recording step, a digital CCD camera was used as a recording device. Holograms contained many unwanted images or noise. To get clean holograms, digital image processing techniques were adopted. In the velocity extraction routine, we improved the HPTV algorithm to extract 3D displacement information of tracer particles. In general, the results obtained using HPTV were not fully acceptable due to technical limitations such as low spatial resolution, small volume size, and low numerical aperture (NA). The problems of spatial resolution and NA are closely related with a recording device. As one experimental parameter that can be optimized, we focused on the particle number density. Variation of the reconstruction efficiency and recovery ratio were compared quantitatively with varying particle number density to check performance of the developed in-line DHPTV system. The reconstruction efficiency represented the particle number distribution acquired through the numerical reconstruction procedure. In addition the recovery ratio showed the performance of 3D PTV algorithm employed for DHPTV measurements. The particle number density in the range of C _o = 13–17 particles/mm³ was found to be optimum for the DHPTV system tested in this study.     

用激光三角法重构混凝土断裂面三维轮廓与断裂能分析

采用基于激光三角法的测试方法对混凝土断面表面轮廓进行了三维重构,初步的断面分析结果表明,相对于传统断面分析方法,本试验方法实现了断面的无损、快速和准确的测量分析,在测量到的断面基础上,给出了本试验所得到的混凝土材料断裂能与分维的关系。     

Volume self-calibration for 3D particle image velocimetry

Planar self-calibration methods have become standard for stereo PIV to correct misalignments between laser light sheet and calibration plane. Computing cross-correlation between images from camera 1 and 2 taken at the same time, non-zero disparity vectors indicate rotational and translational misalignments relative to the coordinate system defined by a calibration plate. This approach works well for thin light sheets but fails for extended volumes recorded in 3D-PTV or tomographic PIV experiments. Here it is primarily necessary to correct calibration errors leading to triangulation errors in 3D-PTV or in degraded tomographic volume reconstruction. Tomographic PIV requires calibration accuracies of a fraction of a pixel throughout the complete volume, which is difficult to achieve experimentally. A new volumetric self-calibration technique has been developed based on the computation of the 3D position of matching particles by triangulation as in 3D-PTV. The residual triangulation error (‘disparity’) is then used to correct the mapping functions for all cameras. A statistical clustering method suitable for dense particle images has been implemented to find correct disparity map peaks from true particle matches. Disparity maps from multiple recordings are summed for better statistics. This self-calibration scheme has been validated using several tomographic PIV experiments improving the vector quality significantly. The relevance for other 3D velocimetry methods is discussed.     

Tomographic PIV measurements in a turbulent lifted jet flame

Measurements of instantaneous volumetric flow fields are required for an improved understanding of turbulent flames. In non-reacting flows, tomographic particle image velocimetry (TPIV) is an established method for three-dimensional (3D) flow measurements. In flames, the reconstruction of the particles location becomes challenging due to a locally varying index of refraction causing beam-steering. This work presents TPIV measurements within a turbulent lifted non-premixed methane jet flame. Solid seeding particles were used to provide the 3D flow field in the vicinity of the flame base, including unburned and burned regions. Four cameras were arranged in a horizontal plane around the jet flame. Following an iterative volumetric self-calibration procedure, the remaining disparity caused by the flame was less than 0.2 pixels. Comparisons with conventional two-component PIV in terms of mean and rms values provided additional confidence in the TPIV measurements.     

Order and disorder in particle–liquid flows in a pipe

The spherical expanded polystyrene particle–oil two-phase flow in a vertical pipe was used to simulate the dispersed phase distribution in laminar bubbly flows. A three-dimensional particle image tracking technique was used to track the particles in the flow to study the ordered structure of dispersed phase distribution and its transition to disorder. The ordered structures behaved as particle strings aligned in the flow direction as induced by the flow shear. The structures were quite durable in high liquid velocity flows and dispersed gradually as the liquid velocity decreased. In lower velocity flows, the particles tended to form clusters in the horizontal direction, as predicted by potential theory for spherical bubbles rising in a quiescent inviscid liquid and as observed in experiments on non-shear bubbly water flows.     

Three-dimensional dynamic imaging of sand particles under wheel via gamma-ray camera system

In recent years, attempts have been made to deploy robots for use in various activities such as planetary exploration, post-tsunami seashore reconnaissance, and volcano investigations. These robots may have to move on soft terrain. The movement of sand or soil particles under the wheels or tracks greatly affects the robot’s ability to maneuver. There is a simple but difficult problem with measuring particle movement: the sand and soil particles beneath the surface are not visible. Only 2D visualization techniques that take a surface picture of the ground or use transparent boards are available. A nuclear 3D imaging technique called positron emission particle tracking (PEPT) was developed at the University of Birmingham for this purpose. PEPT detects pairs of gamma rays emitted by a positron-emitting radionuclide of a tracer particle, which produces an image of the tracer. Thus, the overarching goal of this study was to explore the 3D terramechanics between terrain particles and a wheel or track using PEPT. As an initial step, this paper introduces an imaging technique for standard sand under a rotating wheel using PEPT and presents some images of sand particles under various conditions. Absolute displacements along the longitudinal, vertical, and lateral axes are presented.     

Rod-like particles matching algorithm based on SOM neural network in dispersed two-phase flow measurements

A matching algorithm based on self-organizing map (SOM) neural network is proposed for tracking rod-like particles in 2D optical measurements of dispersed two-phase flows. It is verified by both synthetic images of elongated particles mimicking 2D suspension flows and direct numerical simulations-based results of prolate particles dispersed in a turbulent channel flow. Furthermore, the potential benefit of this algorithm is evaluated by applying it to the experimental data of rod-like fibers tracking in wall turbulence. The study of the behavior of elongated particles suspended in turbulent flows has a practical importance and covers a wide range of applications in engineering and science. In experimental approach, particle tracking velocimetry of the dispersed phase has a key role together with particle image velocimetry of the carrier phase to obtain the velocities of both phases. The essential parts of particle tracking are to identify and match corresponding particles correctly in consecutive images. The present study is focused on the development of an algorithm for pairing non-spherical particles that have one major symmetry axis. The novel idea in the algorithm is to take the orientation of the particles into account for matching in addition to their positions. The method used is based on the SOM neural network that finds the most likely matching link in images on the basis of feature extraction and clustering. The fundamental concept is finding corresponding particles in the images with the nearest characteristics: position and orientation. The most effective aspect of this two-frame matching algorithm is that it does not require any preliminary knowledge of neither the flow field nor the particle behavior. Furthermore, using one additional characteristic of the non-spherical particles, namely their orientation, in addition to its coordinate vector, the pairing is improved both for more reliable matching at higher concentrations of dispersed particles and for higher robustness against loss of particle pairs between image frames.     

Holographic particle image velocimetry system for measurements of hairpin vortices in air channel flow

A holographic particle image velocimetry system for investigating hairpin vortices, artificially generated in a subcritical plane Poiseuille air flow, is presented. The optical setup is a modified version of the hybrid scheme, previously employed in turbulent water flows. Accordingly, separate reconstruction of holograms, successively recorded on the same photoplate, is provided by using two reference beams. The positioning of the photoplate within the image of the sample volume accompanied by special alignment procedures, minimizes the apparent displacement caused by the misalignment of the reconstruction waves. A novel method is employed for detecting in-focus particles. Testing the system with a fixed 5 μm diameter wire, results in a corresponding 3D wire image having a diameter of ≈25 μm. Finally, the instantaneous topology and 3D distribution of the two velocity components associated with the hairpin vortex are presented.     

A multiresolution continuum simulation of the ductile fracture process

With the advancement in computational science that is stepping into the Exascale era and experimental techniques that enable rapid reconstruction of the 3D microstructure, quantitative microstructure simulations at an unprecedented fidelity level are giving rise to new possibilities for linking microstructure to property. This paper presents recent advances in 3D computational modeling of ductile fracture in high toughness steels. Ductile fracture involves several concurrent and mutually interactive mechanisms at multiple length scales of microstructure. With serial sectioning tomographic techniques, a digital dataset of microstructure features associated with the fracture process has been experimentally reconstructed. In this study, primary particles are accurately and explicitly modeled while the secondary particles are modeled by a two scale multiresolution continuum model. The present numerical simulation captures detailed characteristics of the fracture process, such as zigzag crack morphology, critical void growth ratios, local stress triaxiality variation, and intervoid ligament structure. For the first time, fracture toughness is linked to multiscale microstructures in a realistic large 3D model.     

Volumetric-correlation PIV to measure particle concentration and velocity of microflows

Volumetric-correlation particle image velocimetry (VPIV) is a new technique that provides a 3-dimensional 2-component velocity field from a single image plane. This single camera technique is simpler and cheaper to implement than multi-camera systems and has the capacity to measure time-varying flows. Additionally, this technique has significant advantages over other 3D PIV velocity measurement techniques, most notably in the capacity to measure highly seeded flows. Highly seeded flows, often unavoidable in industrial and biological flows, offer considerable advantages due to higher information density and better overall signal-to-noise ratio allowing for optimal spatial and temporal resolution. Here, we further develop VPIV adding the capability to measure concentration and increasing the robustness and accuracy of the technique. Particle concentrations are calculated using volumetric auto-correlations, and subsequently the velocities are calculated using volumetric cross-correlation corrected for variations in particle concentration. Along with the ability to calculate the particle concentration profile, our enhanced VPIV produces significant improvement in the accuracy of velocity measurements. Furthermore, this technique has been demonstrated to be insensitive to out-of-plane flows. The velocity measurement accuracy of the enhanced VPIV exceeds that of standard micro-PIV measurements, especially in near-wall regions. The 3D velocity and particle-concentration measurement capability of VPIV are demonstrated using both synthetic and experimental results.     

Numerical investigation of particle lateral migration in straight channel flows using a direct-forcing immersed boundary method

Inertia-induced cross-stream migration has been recently exploited for precise position of particles in confined channel flows. In this work, a three-dimensional finite volume based immersed boundary method has been developed to study the lateral migration and hydrodynamic self-assembly of neutrally-buoyant particles in pressure-driven flows. Simulation results show that, in 2D channel flows, the equilibrium position for a circular particle is closer to the centreline for larger particle Reynolds number due to the increasing flow rate, while in 3D square duct flow, the equilibrium position for a spherical particle is near a face centre and is closer to the wall for larger particle Reynolds number. Self-assembly of a pair of particles is observed in 3D square duct flows but not in 2D channel flows. Mechanisms for the self-assembly are discussed.     

两相流PIV粒子图像处理方法的研究

本文在单相PIV技术的基础上研究了两相流动PIV图像处理方法，采用摸板匹配法和灰度加权标定法对两相粒子进行了识别、区分和标定，采用灰度互相关法对区分后的单相粒子图像进行了处理，应用基于以上方法编制的Windows应用软件，首先对由美国Minnesota大学复杂流动实验室提供的两相流动粒子图片进行了处理，通过对比分析可见，应用本文所采用的方法能对两相粒子进行有效的识别和区分，然后以搅拌槽内液固两相流场为例对此方法进行了应用。     

Target-free Stereo PIV: a novel technique with inherent error estimation and improved accuracy

A novel, accurate and simple stereo particle image velocimetry (SPIV) technique utilising three cameras is presented. The key feature of the new technique is that there is no need of a separate calibration phase. The calibration data are measured concurrently with the PIV data by a third paraxial camera. This has the benefit of improving ease of use and reducing the time taken to obtain data. This third camera also provides useful velocity information, considerably improving the accuracy of the resolved 3D vectors. The additional redundancy provided by this third perspective in the stereo reconstruction equations suggests a least-squares approach to their solution. The least-squares process further improves the utility of the technique by means of the reconstruction residual. Detailed error analysis shows that this residual is an accurate predictor of resolved vector errors. The new technique is rigorously validated using both pure translation and rotation test cases. However, while this kind of validation is standard, it is shown that such validation is substantially flawed. The case of the well-known confined vortex breakdown flow is offered as an alternative validation. This flow is readily evaluated using CFD methods, allowing a detailed comparison of the data and evaluation of PIV errors in their entirety for this technique.     

An experimental investigation of jet flows through a row of forward expanded holes into a mainstream over a concave surface

This study performed detailed measurements of jet flows through a row of forward expanded holes into a mainstream over a concave surface using digital particle image velocimetry. Each of ejected holes had a streamwise inclined angle of 35° bounded on a concave surface with constant radius of 382 mm. The spacing of adjacent holes is 1.5D. The density and the momentum flux ratio of the mainstream to the jet flow were 1.0. Results show detailed 2D mean velocity maps on several horizontal and vertical planes and a 3D streamline pattern of jet mean velocity. The streamlines of 3D mean velocity clearly display different flow characteristics of the ejected jet flow along the transverse direction. In addition, the particle trajectory of a ring enclosing an ejected jet above the injection hole was also presented to show movement of jet.     

Quantitative visualization of compressible turbulent shear flows using condensate-enhanced Rayleigh scattering

This paper describes several flow visualization experiments carried out in Mach 3 and Mach 8 turbulent shear flows. The experimental technique was based on laser scattering from particles of H₂O or CO₂ condensate that form in the wind tunnel nozzle expansion process. The condensate particles vaporize extremely rapidly on entering the relatively hot fluid within a turbulent structure, so that a sharp vaporization interface marks the outer edge of the rotational shear layer fluid. Calculations indicate that the observed thin interface corresponds to a particle size of 10 nm or less, which is consistent with optical measurements, and that particles of this size track the fluid motions well. Further, calculations and experiments show that the freestream concentration of condensate required for flow visualization has only a small effect on the wind tunnel pressure distribution. Statistics based on the image data were compared to corresponding results from probe measurements and agreement was obtained in statistical measures of speed, scale, and orientation of the large-scale structures in the shear layer turbulence. The condensate-enhanced Rayleigh scattering technique is judged to be a useful tool for quantitative studies of shear layer structure, particularly for identifying the instantaneous boundary layer edge and for extracting comparative information on the large-scale structures represented there.     

Speckle tomography of turbulent flows with density fluctuations

The speckle tomography technique is used for reconstructing both large-scale structures in turbulent flows and the microstructure of turbulence. The technique is based on multi-projectional line-of-sight speckle photography measurements with a subsequent computer-assisted tomographic reconstruction of the interior structure of the flowfield. The large-scale structure is reconstructed using the Radon integral equation, and the microstructure is analysed using a statistical approach and a novel Erbeck-Merzkirch integral transform. Digital speckle photography and speckle tomography methods are described. Numerical simulation of the optical technique is performed using digital ray tracing through a turbulent flowfield. The methods are illustrated by the 3D "averaged" temperature fields in turbulent convective flows obtained earlier and by the recent reconstruction of 3D correlation functions of density variations in turbulent flows. Local values of turbulence (Kolmogorov) microscale are evaluated using these correlation functions and the Erbeck-Merzkirch integral transform The precision of the reconstruction and the spatial resolution achieved are analysed.     

基于未标定序列图的形貌测量方法

提出一种形貌测量方法,从未标定序列图中同时提取出特征点的三维坐标和摄像机的运动矩阵,实现对目标物体的三维测量。首先用齐次坐标表达了空间点的投影变换矩阵,并根据该矩阵推导出投影比例式,运用投影深度和投影矩阵之间的循环相关性进行迭代运算,得到比普通双目视觉更加精确的深度数据。其次对分解得到的摄像机运动矩阵和空间点坐标施加度量约束。运用基于因式分解的归一化算法,恢复出欧氏空间里的摄像机运动矩阵和空间点坐标。基于这种方法,实现了曲面标记点的坐标测量,得到了三维形貌数据,误差小于0.16mm。     

3D scanning particle tracking velocimetry

In this article, we present an experimental setup and data processing schemes for 3D scanning particle tracking velocimetry (SPTV), which expands on the classical 3D particle tracking velocimetry (PTV) through changes in the illumination, image acquisition and analysis. 3D PTV is a flexible flow measurement technique based on the processing of stereoscopic images of flow tracer particles. The technique allows obtaining Lagrangian flow information directly from measured 3D trajectories of individual particles. While for a classical PTV the entire region of interest is simultaneously illuminated and recorded, in SPTV the flow field is recorded by sequential tomographic high-speed imaging of the region of interest. The advantage of the presented method is a considerable increase in maximum feasible seeding density. Results are shown for an experiment in homogenous turbulence and compared with PTV. SPTV yielded an average 3,500 tracked particles per time step, which implies a significant enhancement of the spatial resolution for Lagrangian flow measurements.