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

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

水平管内气液两相螺旋流实验研究

为研究水平管内气液两相螺旋流的流动特性,开展了以空气和水为实验介质,含气率为10％～90％,气相折算速度为0.01～3.4m/s,液相折算速度为0.05～2.7m/s的气液两相螺旋流实验.利用高速摄影机记录并参考借鉴相关研究结果分析和划分了不同工况下的流型;给出了水平管内气液两相螺旋流的流型图;研究了不同流速、不同起旋参数对流动特性(压降、流型衰减、螺距、螺旋直径以及流型转换边界等)的影响.实验结论如下:将水平管内气液两相螺旋流的流型划分为螺旋波状分层流、螺旋泡状流、螺旋团状流、螺旋线状流、螺旋轴状流、螺旋弥散流6种;将绘制的流型图与经典Mandhane流型图进行对比,出现了线状流、弥散流和轴状流3种新的流型;泡状流的分布基本不变,层状流的分布发生变化,当气相流速在2m/s以内时是线状流和轴状流,而不是层状流;随着液相流速的提高,管内两相流动的损失逐渐变大,流型的衰减程度变弱,螺旋扭矩逐渐变大,螺旋直径逐渐变小.另外,随着叶轮角度的增大或者随着叶片面积的减小,流型转换边界均向进气量增大的方向推移.而当进气量一定时,随着叶轮角度的增大或者随着叶片面积的减小,同样流型转换边界趋于进水量增大的方向.最后,随着起旋角度的增大或者随着叶片面积的减小,压降均有逐渐变大的趋势.     

水平管气液两相间歇流含气率研究

对大直径水平管内气液二相流动进行了实验，实验结果表明大直径水平管与小直径水平管具有不同的气液二相流动特征。分析和提出了适用于大直径水平管内间歇流液弹含气率模型，其计算值与实验值非常接近。     

下倾管-立管水气严重段塞流数值模拟

针对海洋油气传输中常见的下倾管-立管系统, 采用Brackbill模型模拟气液相界面间表面张力, VOF方法追踪气液两相运动界面, 提出了管内气液两相流数值模拟方法. 在低气液相进口折算速度下, 数值模拟了该种管型下的严重段塞流动现象, 分析了相关物理参数的变化特性. 结果表明, 在严重段塞流下, 管内流型流态、压力、液塞运动速度、立管出口气液相平均速度、下倾管及立管内含气率等均具有明显周期性特征, 而且一个周期内严重段塞流可分为4个阶段, 进而给出了各阶段中相关参数的变化特性. 数值模拟结果与相关文献中的实验结果吻合良好,表明了该数值模拟方法的有效性.      

不同重力下90°弯管内气液两相流流型及流动特性研究

研究了不同重力条件下90°弯管内气液两相流流型分布形态及流动特性. 通过建立90°弯管内气液两相流流动的三维数学物理模型,采用VOF 方法,对10-6g₀, 10-4g₀, 10-2g₀, 1g₀ (g₀= 9.8m/s2) 重力下的90°弯管内气液两相流流型分布特征、截面空隙率、滑速比及气相尾部最大斜向角进行了比较分析. 研究结果表明:所建立的模型能够正确模拟不同重力条件下90°弯管内气液两相流流型和截面空隙率,并得到气液两相弯管二次流与单相二次流的不同特性. 随着重力水平的提高,90°弯管对气相流型的影响作用减弱,气相整体向弯管内侧积聚靠拢,弯管对尾部的斜向作用减弱.      

垂直向下管内两相流泡状-弹状流型转换研究

建立实验系统,在维持管道出口压力为0.2MPa的条件下,对内径分别为15mm、25mm、40mm、65mm的垂直向下管内空气-水气液两相流动进行了实验研究,获得了两相流泡状-弹状流型分布。实验研究发现:管径对于泡状流与弹状流流型特征有较大影响,并且进一步影响流型转换边界,随着管径增加,泡状流-弹状流的流型转换边界向折算气速减小的方向移动。基于理论推导及实验数据,建立了垂直向下管内气液两相流泡状流-弹状流流型转换预测模型,该模型对本文实验工况条件下的垂直向下管内空气-水气液两相流流型转换具有良好的预测效果,预测模型的计算结果与实验数据之间的误差小于10%。     

气液二相流研究概述

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

微重力条件下气/液两相流流型的研究进展

气／液两相流流型是两相流研究领域最基本的课题之一,至今已有数十年的研究历史．但是,由于气／液两相流动现象极为复杂,目前还没有得到一致的结论．近十多年来,利用微重力环境减弱甚至完全消除重力的影响,简化流动中各种不同因素间的相互作用及流型特征,大大促进了对气／液两相流动特征及流型产生与转换机理的研究．同时,微重力条件下的气／液两相流动是空间技术领域必须解决的关键技术问题之一,具有重要的学术意义和重大的应用价值．本文简要总结了微重力条件下气／液两相流流型研究的基本方法以及实验结果和理论进展,指出今后研究中应该注意的一些方向．      

一种骨小管中液体流动产生的流量及切应力模型

骨组织受力变形后其内部液体就会流动,同时在其微观结构——骨单元壁中扩散,并进一步产生一系列与骨液流动相关的物理效应,如流体剪切应力、流动电位等,这些物理效应被细胞感知并做出破骨或成骨等反应,来使骨适应外部载荷环境.鉴于骨组织产生的内部液体流动很难实验测定,理论模拟是目前的主要研究手段.基于骨单元的多孔弹性性质建立了骨小管内部液体的流动模型,该模型将骨单元所受的外部载荷与骨小管内部液体的压力、流速、流量和切应力联系起来,并进一步可以研究其力传导与力电传导机制.骨小管模型的建立分别基于中空和考虑哈弗液体的骨单元模型,并考虑了骨单元外壁的弹性约束和刚性位移约束两种边界条件.最终得到骨单元在外部轴向载荷作用下,骨小管内部液体的流量及流体切应力的解析解.结果表明:骨小管中的液体流量与流体切应力都正比于应变载荷幅值和频率,并由载荷的应变率决定.因此应变率可以作为控制流量和流体切应力的一种生理载荷因素.流量随着骨小管半径的增大而非线性增大,而流体切应力则随着骨小管半径的增大而线性增大.此外,在相同的载荷下,含哈弗液体的骨单元的模型中,骨小管中液体的流量和切应力均大于中空骨单元模型.     

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

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

Experimental study of R134A condensation heat transfer inside the horizontal micro-fin tubes

The Investigation of the two-phase flow patterns and their transitions during the condensation has gained increasing interest and importance from the well-known phenomenon that the heat transfer characteristics are strongly dependent on the flow patterns. Therefore, it is very important to study on which heat transfer enhancement approach is suitable for an individual flow pattern inside a condenser, so that an accurate heat transfer mechanism can be understood that is consistent with the flow patterns. The condensation heat transfer for R134a in the two kinds of in-tube three-dimensional (3-D) micro-fin tubes with different geometries is experimentally investigated. Based on the flow pattern observations, the flow patterns in the Soliman flow regime map are divided into two-flow regimes; one with the vapor-shear-dominant annular regime and the other with the gravitational-force-dominant stratified-wavy regime. The flow regime transition criterion between the annular regime and the stratified-wavy regime is at Fr equal to 2. In the annular regime, the heat transfer coefficients h of the two kinds of in-tube 3-D micro-fin tubes decreases as the vapor quality x decreases. The regressed condensation heat transfer correlation from the experimental data of the annular flow region is obtained. The dispersibility of the experimental data is inside the limits of ±25%. In the stratified-wavy regime, the average heat transfer coefficient h of the two kinds of in-tube 3-D micro-fin tubes increases as the mass flux increases and the number of micro fins in the 3-D micro-fin tube is not the controlling factor for the performance of a condensation heat transfer. The regressed condensation heat transfer correlation of the stratified-wavy flow regime is experimentally obtained. The dispersibility of the experimental data is inside the limits of ±22%. Combined with the criteria of flow pattern transitions, the correlations can be used for the design of a condenser with 3-D micro-fin tubes.     

An analysis of film flow and its application to condensation in a horizontal tube

An analysis is presented for turbulent flow in a liquid film which is being dragged along a horizontal tube by an axial shear at the film surface. While being dragged along, the film drains down the wall due to gravity. The analysis can be applied to a number of horizontal, annular, gas-liquid flow problems and an example is given here of its use in analysing condensation in a horizontal tube. For this problem, the predictions show limited agreement with the little experimental data presently available.     

Effect of void fraction and friction factor models on the prediction of pressure drop of R134a during downward condensation in a vertical tube

The theoretical flow models of homogeneous and separated flow are applied to in-tube condensation to predict the pressure drop characteristics of R134a. The homogeneous flow model is modified by ten different dynamic viscosity correlations and various alternative correlations of total, frictional and momentum pressure drops to take account of the partial condensation inside the tube. Numerical analyses were performed to determine the average and local homogeneous wall shear stresses and friction factors by means of a CFD program. The equivalent Reynolds number model is modified by six different two-phase friction factors to determine the total condensation pressure drop in the separated flow model. The refrigerant side total pressure drops, frictional pressure drops, friction factors and wall shear stresses are determined within a ±30% error band. The importance of using the alternative total, momentum and frictional pressure drop correlations for the homogeneous flow model is also shown.     

Two-phase flow in inclined tubes with specific reference to condensation: A review

Tilting influences the flow patterns and thus the heat transfer and pressure drop during condensation in smooth tubes. However, few studies are available on diabatic two-phase flows in inclined tubes. The purpose of the present paper is to review two-phase flow in inclined tubes, with specific reference to condensation. Firstly, the paper reviews convective condensation in horizontal tubes. Secondly, an overview is given of two-phase flow in inclined tubes. Thirdly, a review is conducted on condensation in inclined tubes. It is shown for convective condensation in inclined tubes that the inclination angle influences the heat transfer coefficient. The heat transfer coefficient can be increased or decreased depending on the experimental conditions, and especially the flow pattern. Under certain conditions, an inclination angle may exist, which leads to an optimum heat transfer coefficient. Furthermore, this paper highlights the lack of experimental studies for the prediction of the inclination angle effect on the flow pattern, the heat transfer coefficient and the pressure drop in two-phase flows during phase change.     

Convective instability of Marangoni-Poiseuille flow under a longitudinal temperature gradient

An exact solution is obtained for the problem of steady flow in a system of two horizontal layers of immiscible fluids with a common interface. The stability of the flow is studied by a linearization method. It is shown that the occurrence of instabilities is due to the different governing parameters of the fluids (thickness, heating conditions, viscous and thermal conductivity of the fluids). It is found that under constant gravity conditions, the perturbations are monotonic, and in zero gravity, oscillatory thermocapillary instability occurs.     

Developed cavitation behind a disk in a vertical tube

The theoretical and experimental investigation of developed cavities in a vertical flow shows [1–3] that the size and shape of the cavities depend in this case to a greater extent on the Froude number than for horizontal flow. Experiments in a longitudinal gravity field are usually made under conditions of a confined flow, which requires analysis of the influence of this circumstance on their results. Such an analysis has been made for horizontal flow confined by walls in both an approximate [4–6] and an exact [7–10] formulation; the problems have been solved for both two- and three-dimensional flows. In the present paper we solve in a nonlinear formulation the axisymmetric cavitation problem of determining the influence of the restriction of flow in a vertical circular tube on a developed cavity. The obtained results describe well the cavities realized under experimental conditions even when there is an appreciable deviation of the shape of the tube section from a circle.     

Dynamic instabilities of boiling two-phase flow in a single horizontal channel

Dynamic instabilities of two-phase flow associated with refrigerant R-11 in a uniformly heated horizontal in-tube boiling system were experimentally investigated. An experimental setup was designed and built to work in a wide range of mass fluxes G [75–1050 kg/(m² s)], heat fluxes q (0–100 kW/m²), and fluid inlet temperatures T_inlet (2–24°C). Dynamic instability data were obtained under various working conditions. The dependence of oscillation amplitude and period on system parameters is discussed, and the boundaries of various oscillations are located on the steady-state characteristic curves.     

Stratification of boiling two-phase flow in a single horizontal channel

The objective of this study is to investigate experimentally the stratification phenomena of boiling two-phase flow in a uniformly heated horizontal channel. Two-phase flow stratification due to gravity effects, and consequently its thermal and hydrodynamic behavior, under steady state conditions, have been determined by measuring 16 top and 16 bottom wall temperatures. Six distinct wall temperature profiles are found, and the corresponding flow patterns are discussed. A dimensionless number has been formulated for the prediction of the occurrence of different flow patterns.     

Viscous Spreading of Non-Newtonian Gravity Currents on a Plane

A gravity current originated by a power-law viscous fluid propagating on a horizontal rigid plane below a fluid of lower density is examined. The intruding fluid is considered to have a pure Ostwald power-law constitutive equation. The set of equations governing the flow is presented, under the assumption of buoyancy-viscous balance and negligible inertial forces. The conditions under which the above assumptions are valid are examined and a self-similar solution in terms of a nonlinear ordinary differential equation is derived. For the release of a time-variable volume of fluid, the shape of the gravity current is determined numerically using an approximate analytical solution derived close to the current front as a starting condition. A closed-form analytical expression is derived for the special case of the release of a fixed volume of fluid. The space-time development of the gravity current is discussed for different flow behavior indexes.     

Instability of Hadley–Prats Flow with Viscous Heating in a Horizontal Porous Layer

The onset of convective rolls instability in a horizontal porous layer subject to a basic temperature gradient inclined with respect to gravity is investigated. The basic velocity has a linear profile with a non-vanishing mass flow rate, i.e., it is the superposition of a Hadley-type flow and a uniform flow. The influence of the viscous heating contribution on the critical conditions for the onset of the instability is assessed. There are four governing parameters: a horizontal Rayleigh number and a vertical Rayleigh number defining the intensity of the inclined temperature gradient, a Péclet number associated with the basic horizontal flow rate, and a Gebhart number associated with the viscous dissipation effect. The critical wave number and the critical vertical Rayleigh number are evaluated for assigned values of the horizontal Rayleigh number, of the Péclet number, and of the Gebhart number. The linear stability analysis is performed with reference either to transverse or to longitudinal roll disturbances. It is shown that generally the longitudinal rolls represent the preferred mode of instability.