气液二相流研究概述

气液二相流研究液体和气体(或蒸气)两相介质共存条件下的流动特性。二相体系可以是液体中含有气体或(蒸气)泡或者气体,(或蒸气)中含有液体微滴,也可能因其中气泡或液滴的聚并,两相间形成更复杂的分布状态。气液二相流是自然界、日常生活中常见的现象,在许      

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

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

离散型湍流多相流动的研究进展和需求

离散型多相流动,指气体-颗粒(气-固)、液体-颗粒(液-固)、液体-气泡、气体-液雾以及气泡-液体-颗粒等两相或三相流动.这种类型的多相流动广泛存在于能源, 航天和航空, 化工和冶金,交通运输, 水利, 核能等领域.本文阐述了离散型多相流动的国内外基础研究,包括颗粒/液滴/气泡在流场中受流体动力作用力的研究, 颗粒-颗粒,液滴-液滴,气泡-气泡之间以及颗粒/液滴和壁面之间碰撞和聚集规律的研究,颗粒-气体和气泡-液体湍流相互作用的研究, 和数值模拟的研究,包括多相流动的雷诺平均模拟、大涡模拟和直接数值模拟的研究进展.最后, 归纳了目前尚待研究的需求.      

分级进风燃烧室内含颗粒的湍流反应流的实验研究

应用三维激光粒子动态分析仪(PDA)，对分级进风燃烧室内含颗粒的湍流反应流的气 固两相瞬时速度场进行了实验测量，同时对单相湍流反应流的瞬时速度场也进行了实验测 量. 得到了两种情况下燃烧室内气固两相和单相气体的平均轴向与切向速度、轴向与切向脉 动速度均方根值和轴向-切向脉动速度二阶关联量的分布.     

气泡-水流两相流的激光多普勒法测量

为区分激光多普勒稀态两相流动测量中出现的流体相信号和颗粒(气泡、液滴)相信号,本文提出可见度和基座的复合判别法,并设计了专用的信号预处理器,此预处理器能将原始多普勒信号分为两路,分别代表两相,并可用计数型信号处理器测量每路的速度及其它统计量。使用本方案对矩形管道內低速气泡-水两相流的实验和测量表明,对两相的区分和测量是成功的。文中给出了两相速度、湍流度和气泡滑移速度的测量结果。     

气液两相流中气泡运动速度场的PIV分析与研究

粒子图像测速技术（PIV）作为一种无扰、瞬态、全场速度测量方法，已被广泛应用于液体或气体的单相流流速场测定。对于两相流PIV技术，目前还处于起步阶段，本文应用PIV技术的基本原理，对静止液体中的气泡运动速度进行了分析，并对有关气液两相流测量问题进行了探讨。     

节水水嘴起泡器气液两相流数值模拟

为了研究节水水嘴起泡器内部两相流的流动规律,采用Fluent软件对其内部流场进行数值模拟.根据起泡器内部流场的流动特性,采用欧拉两相流模型以及RNG(re—normalization group)κ-ε湍流模型,分析起泡器出口截面气液两相体积分数和速度的分布特点.结果表明:增大入口水流速度可以加快分散出口截面气液两相的分布,缩短流体流动的稳定时间;整流网具有分散流体,降低流速的作用;错开整流网相邻层之间的网格可以改善出口截面的液相分布;本模型中整流网采用三层网格达到较好的出水效果.     

航空发动机轴承腔润滑的气液两相均匀流研究

基于轴承腔中润滑油的两相均匀流动模型,采用湍流模式和有限差分数值方法计算轴承腔内三维定常N-S方程,对腔内润滑油的气液两相均匀流动特性进行研究,以获得气液两相均匀流条件下润滑油流场、压力场和速度场在轴承腔内的分布情况,分析转子转速、含气率和润滑油进口速度对润滑油出口压力以及润滑油与壁面之间剪切力的影响,同时对单相流和两相均匀流润滑性能差异进行比较.结果表明,转子转速、含气率和润滑油进口速度对润滑油出口压力和腔内壁面与润滑油间的平均剪切力具有不同影响,而采用2种流动模型计算出的轴承腔润滑油出口压力的差异较大,同时支持了开展航空发动机轴承腔润滑两相流动分析的必要性.     

微重力气液两相流动与池沸腾传热

综述了近年来中国科学院微重力重点实验室(国家微重力实验室)完成的一系列微重力气液两相流动与池沸腾传热方面的地基实验、飞行实验和理论研究等方面获得的主要成果.在微重力气液两相流动方面,提出了半理论Weber数模型用于预测微重力条件下气液两相弹-环状流转换,并采用Monte Carlo方法,针对气泡初始尺寸对泡-弹状流转换的影响进行数值研究.通过俄罗斯"和平号"空间站与IL-76失重飞机实验,获得了微重力下的气液两相流型图,与此同时在地面利用小尺度毛细管模型模拟了微重力气液两相流动特征.实验测量了微重力气液两相流压降,并基于微重力流动特性建立了一个泡状流压降关联模型.在微重力池沸腾传热方面,利用我国返回式卫星完成了两次空间实验,其中,第22颗返回式卫星搭载铂丝表面R113池沸腾实验采用控制温度的稳态加热方式,而实践8号育种卫星搭载平面FC-72池沸腾实验则采用控制加热电压的准稳态加热方式.同时,还进行了地面常重力和落塔短时微重力条件下的对比实验研究.观察到丝状加热表面微重力时轻微的传热强化现象,而平板加热表面微重力核态池沸腾低热流时传热强化、高热流时传热恶化.微重力实验中观察到气泡脱落前存在横向运动现象,据此分析了气泡行为与传热之间关系,并提出了一个预测丝状加热表面气泡脱落直径的半理论模型.旨在对相关领域的进一步发展和空间两相流系统的应用提供数据及理论支持.     

绕水翼超空化流动形态与速度分布

为揭示超空化流场结构特性,利用高速全流场显示技术,观察了绕hydronautics水翼的超空化流动形态,并利用数字粒子图像测速仪(DPIV)测量了其速度分布. 在测量空穴内部流速分布时,采用空化流场中的空化泡作为示踪粒子来显示流动结构. 结果表明:随着空化数的降低,超空化流动显现出了明显的阶段特征,其中水汽混合相和汽相的分布决定了空化区域的形态与流速分布;空化区和主流区的汽液交界面处存在着较大的速度梯度;低速分布区域随着空化数的降低由水翼吸力面中后部向水翼下游移动;在空化区域内部,水汽混合区的速度相对较低,而汽相区则与主流区有着相近的速度分布.关键词超空化水翼、DPIV、高速摄像、空化形态、流速分布      

The effect of bubble-induced liquid flow on mass transfer in bubble plumes

The use of microbubbles to enhance mass transfer in a compact bubble column has become a valuable topic recently. When the liquid flow induced by the presence of microbubbles is taken into account, the behavior of the microbubbles may differ widely from simple estimations. One example is the change of the residence time, which is determined not only by slip velocity but also the velocity of the surrounding liquid. In the present study the effect of the bubble-induced liquid flow on mass transfer in microbubble plumes is analyzed numerically. A two-way coupling Eulerian–Lagrangian approach is used to simulate oxygen bubble plumes with initial bubble diameters from 100 μm to 1 mm and a maximum local void fraction of less than 2% in compact rectangular tanks. The simulations illustrate that the effect of bubble-induced liquid velocity on the residence time of microbubbles increases with the decrease of initial bubble diameters, and also increases with the reduction of initial water depth. The differences between the concentrated and uniform bubble injections are compared. The results show that the uniform injection of microbubbles provides much better mass transfer efficiency than the concentrated injection, because the bubble-induced liquid flow is suppressed when bubbles are injected uniformly over the entire bottom of the tank.     

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.     

Motion of large gas bubbles ascending in a liquid

An equation is derived for the ascent velocity of large gas bubbles in a liquid. This velocity is assumed to be governed by the propagation of a wavelike perturbation caused by the bubble in the liquid.Notation w bubble (or drop) velocity - specific gravity - dynamic viscosity - kinematic viscosity - r bubble (or drop) radius - surface tension - coefficient of friction - g gravitational acceleration - D bubble (or drop) diameter - p pressure - c propagation velocity of the wavelike perturbation - wavelength     

Physical and transport properties governing bubble column operations

In this work, physical and transport properties governing fluid dynamics in bubble columns were studied in the heterogeneous liquid circulation regime of operation based on several published experimental data. It was confirmed that the axial dispersion phenomenon of the liquid phase could be defined by the constant Peclet number predicated on the two-phase slip velocity, which comprised the time-averaged bubbly flow velocity superimposed by the velocity of large bubbles rising in the fluid. In addition, the slip velocity of mean bubbles in the bubbly flow not only indicated a constant rising velocity of a particular liquid but also apparently reduced to that for the Higbie model with respect to the bubble diameter and its terminal rising velocity, when applied to Akita and Yoshida's dimensionless correlation for mass transfer operation (Akita and Yoshida, 1974). In the course of demonstration, it was confirmed that a semi-theoretical correlation for the gas holdup proposed by Mersmann (1978) is similar to the empirical correlation proposed by Akita and Yoshida (1973) with respect to the expression and the physical perspective.     

On the interaction of Taylor bubbles rising in two-phase co-current slug flow in vertical columns: turbulent wakes

An experimental study on the interaction between Taylor bubbles rising through a co-current flowing liquid in a vertical tube with 32 mm of internal diameter is reported. The flow pattern in the bubble's wake was turbulent and the flow regime in the liquid slug was either turbulent or laminar. When the flow regime in the liquid slug is turbulent (i) the minimum distance between bubbles above which there is no interaction is 5D-6D; (ii) the bubble's rising velocity is in excellent agreement with the Nicklin relation; (iii) the experimental values of the bubble length compare well with theoretical predictions (Barnea 1990); (iv) the distance between consecutive bubbles varied from 13D to 16D and is insensitive to the liquid Reynolds number. When the flow regime in the liquid slug is laminar (i) the wake length is about 5D-6D; (ii) the minimum distance between bubbles above which there is no interaction is higher than 25D; (iii) the bubble's rising velocity is significantly smaller than theoretical predictions. These results were explained in the light of the findings of Pinto et al. (1998) on coalescence of two Taylor bubbles rising through a co-current liquid. Received: 2 February 2000 / Accepted: 15 March 2001     

The minimum in-line coalescence height of bubbles in non-Newtonian fluid

The minimum in-line coalescence height of bubbles generated from a submerged nozzle was investigated experimentally in shear thinning non-Newtonian fluid at lower Reynolds number (2∼60). Carboxymethyl cellulose sodium (CMC) aqueous solution and carbon dioxide were used as the liquid phase and the gas phase, respectively. The process of the formation, movement and in-line coalescence of bubbles was visualized and recorded by a high-speed digital camera. The influences of bubble size, bubble generation frequency and liquid property on the minimum in-line coalescence height of bubbles were investigated by changing nozzle diameter, gas flow rate and the mass concentration of CMC aqueous solutions. For a given liquid, the generating frequency and size of bubbles increased but the minimum coalescence height of in-line bubbles decreased when the nozzle diameter and gas flow rate were increased. When the nozzle diameter and gas flow rate were fixed, the shear-thinning effect of CMC aqueous solution became stronger with increasing CMC mass concentration, which led to the increase in both the terminal rise velocity and average acceleration of the trailing bubble, consequently, the minimum in-line coalescence height of bubbles decreased. An empirical correlation for estimating the minimum in-line bubble coalescence height was proposed, the calculating values accords well with experimental data with a mean relative deviation only 7.6%.     

Interaction between Taylor bubbles rising in stagnant non-Newtonian fluids

The interaction between Taylor bubbles rising in stagnant non-Newtonian solutions was studied. Aqueous solutions of carboxymethylcellulose (CMC) and polyacrylamide (PAA) polymers were used to study the effect of different rheological properties: shear viscosity and viscoelasticity. The solutions studied covered a range of Reynolds numbers between 10 and 714, and Deborah numbers up to 14. The study was performed with pairs of Taylor bubbles rising in a vertical column (0.032 m internal diameter) filled with stagnant liquid. The velocities of the leading and trailing bubbles were measured by sets of laser diodes/photocells placed along the column. The velocity of the trailing bubble was analysed together with the liquid velocity profile in the wake of a single rising bubble (Particle Image Velocimetry data obtained from the literature). For the less concentrated CMC solutions, with moderate shear viscosity and low viscoelasticity, the interaction between Taylor bubbles was similar to that found in Newtonian fluids. For the most concentrated CMC solution, which has high shear viscosity and moderate viscoelasticity, a negative wake forms behind the Taylor bubbles, inhibiting coalescence since the bubbles maintain a minimum distance of about 1D between them. For the PAA solutions, with moderate shear viscosity but higher viscoelasticity than the CMC solutions, longer wake lengths are seen, which are responsible for trailing bubble acceleration at greater distances from the leading bubble. Also in the PAA solutions, the long time needed for the fluid to recover its initial shear viscosity after the passage of the first bubble makes the fluid less resistant to the trailing bubble flow. Hence, the trailing bubble can travel at a higher velocity than the leading bubble, even at distances above 90D.     

Experimental investigation of horizontal air–water bubbly-to-plug and bubbly-to-slug transition flows in a 3.81 cm ID pipe

The present study seeks to investigate horizontal bubbly-to-plug and bubbly-to-slug transition flows. The two-phase flow structures and transition mechanisms in these transition flows are studied based on experimental database established using the local four-sensor conductivity probe in a 3.81 cm inner diameter pipe. While slug flow needs to be distinguished from plug flow due to the presence of large number of small bubbles (and thus, large interfacial area concentration), both differences and similarities are observed in the evolution of interfacial structures in bubbly-to-plug and bubbly-to-slug transitions. The bubbly-to-plug transition is studied by decreasing the liquid flow rate at a fixed gas flow rate. It is found that as the liquid flow rate is lowered, bubbles pack near the top wall of the pipe due to the diminished role of turbulent mixing. As the flow rate is lowered further, bubbles begin to coalesce and form the large bubbles characteristic of plug flow. Bubble size increases while bubble velocity decreases as liquid flow rate decreases, and the profile of the bubble velocity changes its shape due to the changing interfacial structure. The bubbly-to-slug transition is investigated by increasing the gas flow rate at a fixed liquid flow rate. In this transition, gas phase becomes more uniformly distributed throughout the cross-section due to the formation of large bubbles and the increasing bubble-induced turbulence. The size of small bubbles decreases while bubble velocity increases as gas flow rate increases. The distributions of bubble size and bubble velocity become more symmetric in this transition. While differences are observed in these two transitions, similarities are also noticed. As bubbly-to-plug or bubbly-to-slug transition occurs, the formation of large elongated bubbles is observed not in the uppermost region of bubble layer, but in a lower region. At the beginning of transitions, relative differences in phase velocities near the top of the pipe cross-section to those near the pipe center become larger for both gas and liquid phases, because more densely packed bubbles introduce more resistance to both phases.     

Interfacial area concentration in gas–liquid bubbly to churn flow regimes in large diameter pipes

This study performed a survey on existing correlations for interfacial area concentration (IAC) prediction and collected an IAC experimental database of two-phase flows taken under various flow conditions in large diameter pipes. Although some of these existing correlations were developed by partly using the IAC databases taken in the low-void-fraction two-phase flows in large diameter pipes, no correlation can satisfactorily predict the IAC in the two-phase flows changing from bubbly, cap bubbly to churn flow in the collected database of large diameter pipes. So this study presented a systematic way to predict the IAC for the bubbly-to-churn flows in large diameter pipes by categorizing bubbles into two groups (group 1: spherical or distorted bubble, group 2: cap bubble). A correlation was developed to predict the group 1 void fraction by using the void fraction for all bubble. The group 1 bubble IAC and bubble diameter were modeled by using the key parameters such as group 1 void fraction and bubble Reynolds number based on the analysis of Hibiki and Ishii (2001, 2002) using one-dimensional bubble number density and interfacial area transport equations. The correlations of IAC and bubble diameter for group 2 cap bubbles were developed by taking into account the characteristics of the representative bubbles among the group 2 bubbles and the comparison between a newly-derived drift velocity correlation for large diameter pipes and the existing drift velocity correlation of Kataoka and Ishii (1987) for large diameter pipes. The predictions from the newly-developed two-group IAC correlation were compared with the collected experimental data in gas–liquid bubbly to churn flow regimes in large diameter pipes and their mean absolute relative deviations were obtained to be 28.1%, 54.4% and 29.6% for group 1, group 2 and all bubbles respectively.     

Comprehensive experimental investigation of the hydrodynamics of large-scale, 3D,oscillating bubble plumes

An extensive study of the most important hydrodynamic characteristics of fairly large-scale bubble plumes was conducted using several measurement techniques and a variety of tools to analyze the data. Particle image velocimetry (PIV), double-tip optical probes (OP) and photographic techniques were extensively applied to measure bubble and liquid velocities, void-fraction and bubble sizes. PIV measurements in a vertical plane crossing the centre of the injector provided the instantaneous velocity fields for both phases, as well as hydrodynamic parameters, such as the movement of the axis of the plume and its instantaneous shape. Statistical studies were performed using image processing to determine the distribution of the apparent instantaneous plume diameter and centreline position. An important finding was that there is little instantaneous spreading of the bubble plume core; the spreading of the time-averaged plume width (as measured from the time-averaged void-fraction and time-averaged liquid velocity fields) is largely due to plume meandering and oscillations. The liquid-phase stress tensor distributions obtained from the instantaneous velocity data indicate that, for the continuous phase, these stresses scale linearly with the local void-fraction in the range of 0.5% < α < 2.5%. The bubbles were found to be ellipsoidal, with shape factor e ≈ 0.5.