水气两相流中稀疏气泡流动速度场的数字图像测量初探

由空压机提供的气体通过—排微小直径的喷嘴进入静止水体，形成水气两相流流场。在单相PIV和PTV技术的基础上，研究稀疏气液两相流情况下气泡的速度场分布。PIV算法采用快速傅立叶互相关分析法，而PTV算法需要获得每幅图像中每个气泡的形心，根据连续图像中的粒子对，计算速度。用PIV和PTV两种算法处理求出气泡的速度并对两种方法进行比较，其最终研究成果可应用于流体及多相流的流量测技术，提高我们进行低密度气液两相流相关研究的测量水平。同时为水气两相流的数值分析和理论研究提供流场测试的数据。     

两相流显微PIV/PTV系统的开发

开发了一个能同时测量两相流中两相速度和细颗粒尺寸分布的显微PIV/PTV系统,其硬件系统包括大功率连续激光器、显微镜、高速摄像机;软件系统由改进的球形颗粒图像识别算法、各种图像处理算法和各种先进的PIV/PTV算法组成。其中改进的圆弧识别算法能够进行更准确地进行曲线分割而能对充满噪音并相互重叠的颗粒图像给出较好的识别结果。应用该PIV系统,可以在微秒和微米数量级上捕获细颗粒/气泡图像,并能较准确地同时得到两相速度、颗粒尺寸和浓度分布。对焚香可吸入颗粒物进行了速度和尺寸的同时测量,得到了较满意的结果。     

基于光流算法的粒子图像测速技术研究和验证

光流测量技术作为一种新的空气动力学实验技术,以其像素级分辨率的矢量场测量优势获得广泛的应用。光流测量技术使用光流约束方程,配合平滑限定条件,可以进行速度场测量,获得高分辨率的全局矢量场。本文首先通过研究积分最小化光流测速理论和算法,采用C++编写光流速度测量程序;然后通过三种典型的人工位移图像对光流计算程序进行了验证,并将结果和标准位移分布进行比对分析,以指导如何在实际应用中获得高精度光流速度场;最后进行小型风洞后向台阶实验,利用高速相机拍摄示踪粒子图像,使用光流计算程序获得速度矢量场,同采用互相关算法的粒子图像测速计算结果相比较,体现出光流计算方法像素级分辨率的矢量场测量优势。     

稠密气固两相湍流流动的实验和数值模拟

基于气固两相流动模型计算循环流化床内稠密气固两相流湍流动,颗粒动理学方法模拟颗粒相湍动能,ＳＧＳ模型模拟气相湍流,采用γ－射线密度计和非等速取样管测量局部颗粒浓度和流率,利用ＦＦＴ方法计算颗粒浓度功率谱密度。模拟计算得到上升管内气相和固相速度和浓度分布等。同时数值模拟与Ｔｓｕｊｉ等和Ｋｎｏｗｌｔｏｎ等试验结果进行了比较,结果表明数值模拟计算与实验结果相吻合。     

复杂气固两相系统的微观结构

流化床中的气固两相流动是一个高度复杂的非线性混沌系统。本文利用激光粒子动态分析仪（ＰＤＡ）得到的循环流化床中颗粒脉动速度信号,采用ＦＦＴ分析了脉动信号的宽频谱特征,在此基础上应用小波法分析了脉动信号的动态特征,得到了颗粒脉动速度的微观结构,指出颗粒脉动速度的非线性特性是流化床具有混沌特性的根源,且在不同的尺度上颗粒脉动速度表现出各向异性的特征。     

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

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

基于Harris角点检测的位移测量算法

在力学实验中,位移和应变的测量是最基本任务之一。本文提出了一种基于Harris角点检测的位移测量算法。通过检测变形前图像的角点,然后利用光流跟踪技术在变形图像中搜索其匹配点,最终计算得到位移值。算法对图像采用网格划分方式进行计算,达到了全场测量的目的。实验结果表明,该算法位移测量精度高、稳定性好,可作为基于图像处理的无损测量方法。     

混凝土搅拌车搅拌筒磨损数值分析

混凝土运输过程中搅拌筒的磨损一直是一种常见潜在危害,磨损严重时会导致叶片失效,对搅拌的质量和出料匀质性产生影响。通过实验获取搅拌筒内部的磨损费时费力,因此有必要采取一种数值分析方法对搅拌筒的磨损进行预测并提出改进。本文采用摩擦磨损实验的方法来标定颗粒与搅拌筒之间的Archard磨损常数,采用JKR接触模型表征混凝土的流动性能,采用离散元方法(DEM)对搅拌筒筒体及叶片磨损进行预测分析。通过法向接触能量与切向接触能量的对比,证明搅拌筒中的磨损主要为伴有冲击作用的磨粒磨损,搅拌车搅拌筒中搅拌叶片顶部的磨损较为严重。针对磨损比较严重的叶片顶部进行改进,采用T型耐磨结构等改进叶片结构,叶片顶部结构改进后搅拌筒使用寿命能得到明显提升。     

基于数字粒子图像的细水雾全场速度测量

水文通过利用激光片光照明细水雾喷雾场和采用数字成像系统获取了喷雾粒子的运动轨迹图像，进而通过研制相应的图像处理和分析软件，重建了细水雾雾场粒子的速度分布，本方法是直接对流场中的粒子成像。因而测量时无需施加示踪粒子对数字成像系统的图像获取速度要求不高，且在图像数据的分析处理方面不像DPIV技术的互相关算法，对示踪粒子的数密度有一定的限制，因此，本方法可用于不宜添加示踪粒子或流场中粒子数密度较低的喷雾场或其它两相流场中粒子速度场的测量。     

一种地磁矢量测量多源误差校正方法

地磁矢量测量中的误差来源众多,而传统误差校正算法多是对其中部分误差进行校正,导致校正后的磁场数据仍存在剩余误差,在对各项误差的特性进行分析的基础上,提出了一种地磁矢量测量多源误差综合校正方法。首先,基于椭球约束法对地磁总场误差进行校正,得到了正交坐标系下的地磁矢量;然后,将椭球约束法校正后剩余的分量误差与非对准误差合并表述为一个三维旋转矩阵,利用遗传算法估算其中的旋转角度并完成误差校正。实验结果表明,所提出的方法可以对地磁矢量测量的多项误差进行有效校正,经校正后地理参考系下地磁三分量误差的改善率分别达到91.2%、95.2%和95.7%,与传统算法相比校正效果明显提升。     

Measurement and simulation of the two-phase velocity correlation in sudden-expansion gas-particle flows

In this paper the present authors measured the gas-particle two-phase velocity correlation in sudden expansion gas-particle flows with a phase Doppler particle anemometer(PDPA) and simulated the system behavior by using both a Reynolds-averaged Navier-Stokes(RANS) model and a large-eddy simulation(LES).The results of the measurements yield the axial and radial time-averaged velocities as well as the fluctuation velocities of gas and three particle-size groups(30 μ m,50 μ m,and 95 μ m) and the gas-particle velocity correlation for 30 μ m and 50 μ m particles.From the measurements,theoretical analysis,and simulation,it is found that the two-phase velocity correlation of sudden-expansion flows,like that of jet flows,is less than the gas and particle Reynolds stresses.What distinguishes the two-phase velocity correlations of sudden-expansion flow from those of jet and channel flows is the absence of a clear relationship between the two-phase velocity correlation and particle size in sudden-expansion flows.The measurements,theoretical analysis,and numerical simulation all lead to the above-stated conclusions.Quantitatively,the results of the LES are better than those of the RANS model.     

Turbulent shear stress profiles in a bubbly channel flow assessed by particle tracking velocimetry

Particle tracking velocimetry (PTV) is applied to a bubbly two-phase turbulent flow in a horizontal channel at Re = 2 × 10⁴ to investigate the turbulent shear stress profile which had been altered by the presence of bubbles. Streamwise and vertical velocity components of liquid phase are obtained using a shallow focus imaging method under backlight photography. The size of bubbles injected through a porous plate in the channel ranged from 0.3 to 1.5 mm diameter, and the bubbles show a significant backward slip velocity relative to liquid flow. After bubbles and tracer particles are identified by binarizing the image, velocity of each phase and void fraction are profiled in a downstream region. The turbulent shear stress, which consists of three components in the bubbly two-phase flow, is computed by analysis of PTV data. The result shows that the fluctuation correlation between local void fraction and vertical liquid velocity provides a negative shear stress component which promotes frictional drag reduction in the bubbly two-phase layer. The paper also deals with the source of the negative shear stress considering bubble’s relative motion to liquid.     

Simultaneous velocity field measurements in two-phase flows for turbulent mixing of sprays by means of two-phase PIV

This paper describes a novel derivative of the PIV method for measuring the velocity fields of droplets and gas phases simultaneously using fluorescence images rather than Mie scattering images. Two-phase PIV allows the simultaneous and independent velocity field measurement of the liquid phase droplets and ambient gas in the case of two-phase flow mixing. For phase discrimination, each phase is labelled by a different fluorescent dye: the gas phase is seeded with small liquid droplets, tagged by an efficient fluorescent dye while the droplets of the liquid phases are tagged by a different fluorescent dye. For each phase, the wavelength shift of fluorescence is used to separate fluorescence from Mie scattering and to distinguish between the fluorescence of each phase. With the use of two cross-correlation PIV cameras and adequate optical filters, we obtain two double frame images, one for each phase. Thus standard PIV or PTV algorithms are used to obtain the simultaneous and independent velocity fields of the two phases. Because the two-phase PIV technique relies on the ability to produce two simultaneous and independent images of the two phases, the choice of the labelling dyes and of the associated optical filter sets is relevant for the image acquisition. Thus a spectroscopic study has been carried out to choose the optimal fluorescent dyes and the associated optical filters. The method has been evaluated in a simple two-phase flow: droplets of 30–40 μm diameter, produced by an ultrasonic nozzle are injected into a gas coflow seeded with small particles. Some initial results have been obtained which demonstrate the potential of the method.     

Effect of drag force on backward-facing step gas-particle turbulent flows

An improved drag force coefficient of gas-particle interaction based on the traditional Wen’s 1966 model is proposed. In this model, a two-stage continuous function is used to correct the discontinuous switch when porosity less than 0.2. Using this proposed correlation and the Wen’s 1966 model, a gas-particle kinetic energy and particle temperature model is developed to predict the hydrodynamic characteristics in backward-facing step gas-particle two-phase turbulent flows. Numerically results showed that they are in good agreement with experiment measurements and presented model are better due to a improvement of momentum transport between gas and particle phases. Particle dispersions take on the distinctively anisotropic behaviors at every directions and gas phase fluctuation velocity are about twice larger than particle phases. Particle phase has a unique transportation mechanism and completely different from the gas phase due to different density. Furthermore, the correlation values of axial–axial gas-particle are always greater than the radial–radial values at fully flow regions. The gas-particle two-phase interactions will make influence on two-phase turbulent flow behaviors.     

Theoretical and numerical studies on the flow multiplicity phenomenon for gas–solids two-phase flows in CFB risers

The dependence of the fully-developed flow profiles on the inlet flow conditions for gas–solids two-phase flows, i.e. the flow multiplicity phenomenon, in circulating fluidized bed (CFB) risers was proposed and discussed in this article. The flow multiplicity phenomenon for gas–solids two-phase flows was first proved mathematically based on the conservation equations of mass and momentum. Then the CFD model using Eulerian–Eulerian approach with k–ε turbulence model for each phase was further adopted to analyze the details of this flow multiplicity phenomenon. It is theoretically and numerically revealed that for gas–solids two-phase flows, the flow profiles in the fully-developed region are always dominated by the flow profiles at the inlet. The solids concentration profile is closely coupled with the velocity profile, and the inlet solids concentration and velocity profiles can largely influence the fully-developed concentration and velocity profiles.     

Two-phase turbulence models for simulating dense gas-particle flows

The two-fluid model is widely adopted in simulations of dense gas-particle flows in engineering facili- ties. Present two-phase turbulence models for two-fluid modeling are isotropic. However, turbulence in actual gas-particle flows is not isotropic. Moreover, in these models the two-phase velocity correlation is closed using dimensional analysis, leading to discrepancies between the numerical results, theoretical analysis and experiments. To rectify this problem, some two-phase turbulence models were proposed by the authors and are applied to simulate dense gas-particle flows in downers, risers, and horizontal channels; Experimental results validate the simulation results. Among these models the USM-O and the two-scale USM models are shown to give a better account of both anisotropic particle turbulence and particle-particle collision using the transport equation model for the two-phase velocity correlation.     

Evaluation of the effect of wall boundary conditions on numerical simulations of circulating fluidized beds

A computational fluid dynamics (CFD) modeling of the gas–solids two-phase flow in a circulating fluidized bed (CFB) riser is carried out. The Eularian–Eularian method with the kinetic theory of granular flow is used to solve the gas–solids two-phase flow in the CFB riser. The wall boundary condition of the riser is defined based on the Johnson and Jackson wall boundary theory (Johnson & Jackson, 1987) with specularity coefficient and particle–wall restitution coefficient. The numerical results show that these two coefficients in the wall boundary condition play a major role in the predicted solids lateral velocity, which affects the solid particle distribution in the CFB riser. And the effect of each of the two coefficients on the solids distribution also depends on the other one. The generality of the CFD model is further validated under different operating conditions of the CFB riser.     

Evaluation of the effect of wall boundary conditions on numerical simulations of circulating fluidized beds

A computational fluid dynamics （CFD） modeling of the gas-solids two-phase flow in a circulating fluidized bed （CFB） riser is carried out. The Eularian-Eularian method with the kinetic theory of granular flow is used to solve the gas-solids two-phase flow in the CFB riser. The wall boundary condition of the riser is defined based on the Johnson and Jackson wall boundary theory （Johnson ＆ Jackson, 1987） with specularity coefficient and particle-wall restitution coefficient.The numerical results show that these two coefficients in the wall boundary condition play a major role in the predicted solids lateral velocity, which affects the solid particle distribution in the CFB riser. And the effect of each of the two coefficients on the solids distribution also depends on the other one. The generality of the CFD model is further validated under different operatin~ conditions of the CFB riser.     

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.     

MULTIFRACTAL ANALYSIS OF PARTICLE-FLUID SYSTEM IN A CIRCULATING FLUIDIZED BED

In this paper, multifractal analysis together with wavelet transform modulus maxima (WTMM) method is used to analyze the fluctuating signals of circulating fluidized bed (CFB). Singularity spectrum shows that the gas-particle flow in CFB has multifractal character. Motion behavior of the particle-fluid system of CFB can be described by singularity spectrum. Intermittency index can be used to determine the transition of flow regime from fast fluidization to pneumatic conveying.