竖直窄矩形通道内弹状流液弹中分区特性研究

弹状流的液弹部分受气弹尾部影响,其水力特性参数沿流动方向存在分区的不一致性。本文对竖直窄矩形通道中弹状流液弹内参数的分布特性进行了研究。结果表明：液弹内气泡在近壁面附近所受径向力较为平衡,气泡频率较大;随着远离气弹尾部,管道中间气泡频率逐渐增大。根据气泡频率波动变化将液弹分为3个区域,尾流区占液弹长度的40%～45%,过渡区占10%～15%,主流区占40%～50%。尾流区和主流区内,空泡份额呈“三峰型”分布;随着气相流速的增加,尾流区内近壁面处峰值逐渐增大,管道中间峰值逐渐下降,但主流区内情况相反。气泡直径随气相流速的增大而变大,且液弹内气泡聚合和破碎现象较少。     

可变液膜厚度下垂直塞状流参数预测模型

基于流动机理的分析建立了塞状流参数预测模型;模型中考虑了液膜的厚度变化.分析了液膜厚度变化对预测结果产生的影响,并用公开发表的数据对模型进行了验证.分析表明,若忽略液膜厚度的变化,将Taylor泡简化为圆柱体,会使其长度的预测值偏小,导致压力梯度的预测出现正偏差,且偏差会随气相表观速度的增加而增大.新建模型反映了液膜的流动特性,可对不同来源的数据进行较为准确的预测.     

竖直窄矩形通道内空气-水两相流中气弹运动速度研究

以空气和水为工质,应用高速摄像仪,对竖直窄矩形通道(3.25 mm×40 mm)内气液两相弹状流进行了可视化实验研究。气、液相表观速度分别为0.1~2.51 m/s和0.16~2.62 m/s,工作压力为常压。实验中发现窄矩形通道内弹状流与圆管中存在较大差别,气弹多发生变形,高液相流速时变形更为严重。窄边液膜含气量较高,在高液相流速时窄边液膜不下落,宽边液膜中含有由气弹头部进入和气弹尾部进入的气泡。气弹速度受气弹头部形状和宽度影响较大,受气弹长度影响较小。气弹速度可由Ishii & Jones-Zuber模型计算,但在低液相折算速度时偏差较大,其主要原因为漂移速度计算值较实验值偏小。     

Modeling of bubble shape in horizontal and inclined tubes

An experimental and theoretical study on the bubble shape of intermittent flow in the horizontal and inclined pipes has been carried out. The experiment results show that the bubble shape depends on the Froude number, bubble length and pipe inclination. The bubble with staircase pattern tail is observed at low Froude numbers, which is corresponding to plug flow. A model for the prediction of the bubble shape in horizontal and inclined pipes is proposed. The model is able to predict the bubble shape, flow pattern transition between plug and slug flow regimes as well as nose-tail inversion phenomenon observed in the downwardly inclined pipe. Validation shows the model can well predict the bubble shapes in horizontal and inclined pipes. The model discloses that the transition between plug and slug flow regimes occurs within a region. The Froude number range for plug flow regime in the downwardly inclined pipe is much wider than that in the horizontal or upwardly inclined pipe. The assumption of fully developed liquid film under the long bubbles tends to under-estimate the liquid fraction in this part of the slug structure, especially, for the intermittent flow in the upwardly inclined pipe with high Froude numbers.     

Behavior of a Single Bubble in Quiescent and Flowing Liquid inside a Cylindrical Tube

The velocity of a single bubble in quiescent and flowing liquid was studied to gain information on the structure of bubble- and slug-flows of gas-liquid two-phase flow.

In the experiments, tap water was used at room temperature and the bubbles were generated by injecting air with a syringe. The bubble velocities were evaluated from photographs taken by multiflash exposure. The tube diameters adopted were 5, 1 and 0.5 cm for the flowing liquid experiments and 1, 0.5 and 0.2 cm for the quiescent liquid experiments. The average liquid velocities varied from about 0 to 2 m/sec.

The results indicated that the velocity of a single bubble in flowing liquid was the sum of the local liquid velocity in the vicinity of the bubble and of the velocity of rise of bubble in quiescent liquid. For Taylor bubbles, the local liquid velocity is given by the maximum liquid velocity near the center of the channel.

The velocities of spherical and ellipsoidal bubbles are determined by the balance of forces acting on them, but that of the Taylor bubble appears to be determined by a different mechanism. A good explanation of such a mechanism is Taylor instability at the gas-liquid interface well removed from the channel wall.     

Modeling of bubble coalescence and disintegration in confined upward two-phase flow

This paper presents the modeling of bubble interaction mechanisms in the two-group interfacial area transport equation (IATE) for confined gas–liquid two-phase flow. The transport equation is applicable to bubbly, cap-turbulent, and churn-turbulent flow regimes. In the two-group IATE, bubbles are categorized into two groups: spherical/distorted bubbles as Group 1 and cap/slug/churn-turbulent bubbles as Group 2. Thus, two sets of equations are used to describe the generation and destruction rates of bubble number density, void fraction, and interfacial area concentration for the two groups of bubbles due to bubble expansion and compression, coalescence and disintegration, and phase change. Five major bubble interaction mechanisms are identified for the gas–liquid two-phase flow of interest, and are analytically modeled as the source/sink terms for the transport equation in the confined flow. These models include both intra-group and inter-group bubble interactions.     

Two-group interfacial area transport in vertical air-water flow: I. Mechanistic model

The prediction of the dynamical evolution of interfacial area concentration is one of the most challenging tasks in two-fluid model application. This paper is focused on developing theoretical models for interfacial area source and sink terms for a two-group interfacial area transport equation. Mechanistic models of major fluid particle interaction phenomena involving two bubble groups are proposed, including the shearing-off of small bubbles from slug/cap bubbles, the wake entrainment of spherical/distorted bubble group into slug/cap bubble group, the wake acceleration and coalescence between slug/cap bubbles, and the breakup of slug/cap bubbles due to turbulent eddy impacts. The existing one-group interaction terms are extended in considering the generation of cap bubbles, as well as different parametric dependences when these terms are applied to the slug flow regime. The complete set of modeling equations is closed and continuously covers the bubbly flow, slug flow, and churn-turbulent flow regimes. Prediction of the interfacial area concentration evolution using a one-dimensional two-group transport equation and evaluation with experimental results are described in a companion paper.     

基于竖直圆管空气-水两相流实验的相间曳力模型研究

研究两相流相间阻力特性对系统程序关键本构模型封闭具有重要意义。本文基于竖直圆管开展了空气-水两相流实验,采用四探头电导探针对空泡份额、气泡弦长和界面面积浓度等气泡参数的径向分布进行了测量。结果表明空泡份额和气泡弦长呈现“核峰型”分布,而界面面积浓度并没有表现出随流速的单调关系。进一步开发了泡状流和弹状流的相间曳力模型,考虑了液相表观流速与管径对气泡尺寸分布的影响,建立了临界韦伯数与不同液相流速的关系。计算得到的空泡份额和界面面积浓度与实验数据整体符合较好,验证了模型的可靠性,为两相流相间阻力特性研究提供参考意义。     

Analysis of Gas-liquid Two-phase Flow Pattern Transition in Inclined Rising Pipe under Fluctuant Nonlinear Vibration Condition

The flow pattern and transition law of gas-liquid two-phase flow in inclined rising pipe under different fluctuant nonlinear vibration conditions were studied experimentally. The flow patterns of gas-liquid two-phase flow under fluctuant nonlinear vibration condition were identified by using high-speed camera. The results show that there are four kinds of patterns in inclined rising pipe with diffuse bubble flow, fluctuant slug flow, proto slug flow and liquid-ring annular flow. The mechanisms of the transition from diffuse bubble flow to fluctuant slug flow, proto slug flow to liquid-ring annular flow were derived, the vibration parameter was added on the basis of the steady state transition mechanism and the flow pattern transition formulae considering the vibration acceleration were established. The results show that the flow pattern transition formulae established in this paper agree well with the experimental results.     

起伏非线性振动下倾斜上升管内气液两相流流型转变分析

本文对不同起伏非线性振动条件下倾斜上升管内气液两相流流型及转变规律进行实验研究,借助高速摄影仪对起伏非线性振动状态下气液两相流的流型进行分类。结果表明,倾斜上升管气液两相流有弥散泡状流、起伏弹状流、准弹状流和液环式环状流4种。对弥散泡状流向起伏弹状流和准弹状流向液环式环状流的转变机理进行分析,在稳定状态转变机理的基础上引入振动参数,建立了考虑振动加速度的关系式。本文建立的流型转变关系式与实验结果吻合较好。     

Measurement of void fraction in bubbly-slug flow with a constant electric current method

In several void fraction measurement methods, a constant electric current method which is one of conductance methods is focused in the present study. By using this method, void fraction can be measured with higher temporal resolution. However, it has been mainly applied to annular flow in previous studies. In the present study, Maxwell's estimation, Bruggemann's estimation, low void fraction approximation and new estimations which consider the bubble shape are applied in order to measure more accurately void fraction of dispersed bubbly flow and slug flow. To understand the effect of bubble shapes and flow patterns, void fraction was measured by the constant electric current method for a rising single spherical bubble and a rising single slug bubble without a forced convection. In addition, void fraction was also measured in bubbly flow and bubbly-slug flow with a forced convection. Then, effects of flow patterns on the proposed estimations of void fraction and the accuracy of their estimations were discussed with the measurement results. From the result, the new estimations which consider a bubble shape are more accurate than the previous estimation in a slug bubble and bubbly-slug flow.     

摇摆状态下气液两相流流型转变的实验研究

通过可视化观察和数码照片对摇摆状态下光滑有机玻璃管内气液两相流流型进行分类和定义,并分析了不同管径、摇摆角度以及摇摆周期对流型之间转变的影响.结果表明:在液相折算流速一样的情况下,管径增加、摇摆周期缩短或摇摆角度减小会使得环状流形成需要更高的气体折算流速;弹状流向搅混流转变所需气相流量则随着管径的减小、摇摆周期的增加或摇摆角度的减小而增加.而在气相折算流速一样的条件下,管径增加、摇摆周期缩短或摇摆角度增大会使泡状流产生需要更高的液相流量.     

竖直窄矩形通道内弹状流中液膜特性研究

气液两相弹状流广泛存在于工程领域,弹状流中液膜特性对弹状流模型的建立具有重要意义。为此利用高速摄像系统,对竖直窄矩形通道(3.25 mm×40 mm)内弹状流中液膜进行了可视化研究。实验中发现窄矩形通道中气弹左右两侧窄边液膜厚度不等且存在波动,但其对两侧液膜速度影响较小,两侧液膜速度相等。液膜脱离厚度主要受两相流速及气弹长度影响。液膜脱离速度随液相折算速度增加而增大;在低液相流速时,随气相折算速度增加而减小;当液相流速≥1.204 m/s时,液膜不下落,液膜脱离速度随气相速度变化较小,主要受液相流速影响。     

竖直管内弹状流向搅混流转变界限的判定

本文对两种不同管径的竖直管内弹状流向搅混流转变界限的判定进行了实验研究.实验发现,在弹状流向搅混流转变的过程中,管内压力会随气体流量的增加出现先降低后上升的非单调性变化.通过对管内压力和摩阻变化的理论分析,并在可视化观察的基础上,提出压力变化曲线的转折点即为搅混流产生的起始点.然后根据实验数据,绘出两种不同管径下的搅混流转变界限.     

Two-group interfacial area transport in vertical air-water flow -II. Model evaluation

In a companion paper, mechanistic models of major fluid particle interaction phenomena involving two bubble groups have been proposed. The prediction of interfacial area concentration evolution using the one-dimensional two-group transport equation and evaluation with experimental results are performed in the paper. These evaluations are based on solid databases for a 2-inch air-water loop with sufficient information on the axial development and the radial distribution of the local parameters. Model evaluation strategies are systematically analyzed. The predictions for the interfacial area concentration evolution demonstrate satisfactory accuracy. The proposed model predicts a smooth transition across the bubbly-to-slug flow regime boundary and demonstrates mechanisms for the generation and development of the cap/slug bubble group. The two-group interfacial area transport equation covers a wide range from bubbly, slug, to churn turbulent flow regimes for adiabatic air-water upward flow in moderate diameter pipes. The generality of the interfacial transport model is also discussed.     

Isotopic separation of Lithium using Inorganic Ion Exchangers

The authors have been investigating geysering, which may appear during start-up of a natural circulation boiling reactor. From the results, it became clear that a large slug bubble covering the entire channel and its condensation in the upper plenum were indispensable for geysering to occur. In this work, an image processing unit has been developed in order to measure the coalescence process of multiple slug bubbles. The proposed system was used to measure the coalescence process of multiple slug bubbles flow in vertical circular channel on air-water system. It was shown that a following slug bubble was influenced by the preceding slug bubble as the distance between the slug bubbles became 5 times inner channel diameter. The rise velocity of the following slug bubble increased due to the pressure drop induced by the wake behind the preceding slug bubble. The bubbles coalesced together in a fraction of a second.     

Two-group interfacial area transport in vertical air–water flow: I. Mechanistic model

In a companion paper, mechanistic models of major fluid particle interaction phenomena involving two bubble groups have been proposed. The prediction of interfacial area concentration evolution using the one-dimensional two-group transport equation and evaluation with experimental results are performed in the paper. These evaluations are based on solid databases for a 2-inch air–water loop with sufficient information on the axial development and the radial distribution of the local parameters. Model evaluation strategies are systematically analyzed. The predictions for the interfacial area concentration evolution demonstrate satisfactory accuracy. The proposed model predicts a smooth transition across the bubbly-to-slug flow regime boundary and demonstrates mechanisms for the generation and development of the cap/slug bubble group. The two-group interfacial area transport equation covers a wide range from bubbly, slug, to churn turbulent flow regimes for adiabatic air–water upward flow in moderate diameter pipes. The generality of the interfacial transport model is also discussed.     

Liquid entrainment by an expanding core disruptive accident bubble—a Kelvin/Helmholtz phenomenon

The final stage of a postulated energetic core disruptive accident (CDA) in a liquid metal fast breeder reactor is believed to involve the expansion of a high-pressure core-material bubble against the overlying pool of sodium. Some of the sodium will be entrained by the CDA bubble which may influence the mechanical energy available for damage to the reactor vessel. The following considerations of liquid surface instability indicate that the Kelvin–Helmholtz (K–H) mechanism is primarily responsible for liquid entrainment by the expanding CDA bubble. First, an instability analysis is presented which shows that the K–H mechanism is faster than the Taylor acceleration mechanism of entrainment at the high fluid velocities expected within the interior of the expanding CDA bubble. Secondly, a new model of liquid entrainment by the CDA bubble is introduced which is based on spherical-core-vortex motion and entrainment via the K–H instability along the bubble surface. The model is in agreement with new experimental results presented here on the reduction of nitrogen-gas-simulant CDA bubble work potential. Finally, a one-dimensional air-over-water parallel flow experiment was undertaken which demonstrates that the K–H instability results in sufficiently rapid and fine liquid atomization to account for observed CDA gas-bubble work reductions. An important byproduct of the theoretical and experimental work is that the liquid entrainment rate is well described by the Ricou–Spalding entrainment law.     

三面加热窄矩形通道内气液两相流流型研究

以去离子水为工质,对截面为3 mm×43 mm的三面加热窄矩形通道内流型转化过程进行可视化实验研究。借助高速摄影仪记录可视化数据,观察到泡状流、弹状流、搅拌流和气膜塞状流等4种主要流型,并详细描述了各种流型发生时通道内气泡转化的过程。记录不同流型转化时的临界点,绘制出三面加热窄矩形通道的流型图,分析流型图中流型转化边界曲线的趋势及形成机理。将本实验流型图与现有相似通道尺寸流型图进行对比,结果表明：三面加热条件下的流型转化过程与绝热条件下的空气 水流型转化过程差异很大,某些流型转化曲线存在趋势上的不同;由于窄边加热部分的影响,与单面加热通道的流型转化过程也存在明显差异。气膜塞状流在绝热条件和单面加热条件下均未出现。     

基于WMS的无旋和螺旋气液两相流流型及空泡份额实验研究

为探究气液两相流流型从无旋状态转变为螺旋状态前后的流型特征及空泡份额时空分布特性,基于高速摄影仪和自主开发的丝网传感器（WMS）测量技术,对内径为30 mm的水平管内起旋装置作用下空气-水两相流的相态时空演变特性进行了可视化实验研究。结果表明,在起旋器诱导的离心力作用下,流场内存在明显的气泡聚并行为和液滴沉积现象,其中,泡状流将转变为螺旋气柱流,塞状流转变为螺旋间歇流,弹状流转变为螺旋环状流,环状流转变为螺旋丝带流;相比于弹状流和环状流,泡状流和塞状流的截面平均空泡份额在起旋器出口波动幅值明显减弱,但离心力场并未明显改变各流型从无旋状态转变为螺旋状态前后的截面平均空泡份额。