第4章-管内气液两相流阻力计算

Chapter4.管内气液两相流的阻力计算

(Gas-liquidflowresistancecalculation)西安交通大学能源与动力工程学院王树众教授内容概要2引言摩擦压降计算加速压降局部阻力4.14.24.34.5重位压力降4.44.1引言3压力降计算是气液两相流研究中最重要的课题之一只有正确地进行压力降计算，才能使系统具有安全可靠和足够的压头，才能为动力设备的选型以及安全经济运行提供必要的依据。气液两相流的压力降包括四部分，即重位压力降、摩擦压力降、加速压力降和局部阻力压力降，亦即：

△PT=△Pg+△Pf+△Pa+△Pb式中△PT－总压力降△Pg－重位压力降-重力作用而引起的△Pf

－摩擦压力降-摩擦阻力引起的△Pa

－加速压力降-流体速度变化而引起的△Pb－局部阻力压力降-流动方向或管截面发生变化引起.4.2摩擦压降计算1-摩擦压降计算的经验方法2-依据流型的摩擦压降的计算3-影响气液两相流体摩擦阻力压力降的主要因素4.2摩擦压降计算1-摩擦压降计算的经验方法均相流模型

对单相：

（4-1）将两相看作均匀混合的介质：

（4-2）又混合相密度：

（x为平均干度）（4-3）

则上式成为：

（4-4）苏联计算方法均相模型（均相模型通常计算值偏低，有时差别还相当大）苏联50年代锅炉水循环计算法该式计算误差对水平均相为（>20%）4.2摩擦压降计算1-摩擦压降计算的经验方法均相流模型鉴于公式4-4误差较大1978年苏联修订后的锅炉水力计算标准状况为：

（4-5）x-为平均干度，考虑了校正系数

，其与x，p，m

有关，可查图求得。对于受热管上式中x用平均值

，

修正系数值可按照下式计算

（4-6）

出口、进口处的干度

苏联计算方法4.2摩擦压降计算1-摩擦压降计算的经验方法均相流模型欧美计算方法欧美采用的均相模型计算法和苏联采用的计算式方法略有不同采用均相模型计算

必须用一合适的摩擦阻力系数

苏联水动力计算标准采用单相流体的摩擦阻力系数，且认为流动工况已进入阻力平方区，

值已经与Re数无关，与相对粗糙度D/K有关其计算可按下列尼古拉兹计算式（4-7）D及K——分别为管子内直径及管壁粗糙度欧美则都采用勃拉休斯的光滑管计算式值，其公式为：（4-8）4.2摩擦压降计算1-摩擦压降计算的经验方法均相流模型欧美计算方法即

（4-9）按此法计算，当干度x=1时，不能正确反映全部为气体流过时的数值，因为在此计算法中

的函数不能转化为气体雷诺数的函数。为了避免这一不足，有些作者采用一个平均的两相动力黏度

，其值和干度x的关系应能满足：当x=0时，

当x=1时，

的条件。欧美在计算式4-8中的Re数时又有多种算法①Re数按液体的参数计算4.2摩擦压降计算1-摩擦压降计算的经验方法均相流模型欧美计算方法①Re数中的μ值采用平均的两相动力粘度计算不同作者采用的

值也不相同麦卡达姆（Mcadams）采用

希奇蒂（Cichitti）采用

杜克勒（Dukler）采用

班可夫（Bankoff）采用

如果使用麦卡达姆的式进行计算，则

值可按下式算得：

（4-10）应用式（4-10）算得的

值比按式（4-4）算得的低，较适宜用于高质量流速工况。气液两相完全分开流动不考虑两相间相互作用单位管长上的气、液两相的摩擦压力降相等，且等于两相管流的摩擦压力降104.2摩擦压降计算1-摩擦压降计算的经验方法分相流模型①

Lockhart—Martinelli和Martinelli—Nelson计算方法在1944-1947年间由Lockhart和Martinelli等提出，建立在分相模型流动基础上假设114.2摩擦压降计算1-摩擦压降计算的经验方法分相流模型①-Lockhart—Martnelli和Martinelli—Nelson计算方法推导-应用单相流动的Darcy公式由于

所以

（1）同理：

（2）其中，

、

分别代表混合物中的气相和混合物中的液相注：在英美文献中，定义摩擦阻力系数f为:摩擦切应力

，而

（其中p为管子周长）。根据受力平衡知上式等于

而

，显然

。、

分别为两相流中的气相和液相与管壁的摩擦阻力系数（指在其流动截面AG、AL中流速时）124.2摩擦压降计算1-摩擦压降计算的经验方法分相流模型①-Lockhart—Martnelli和Martinelli—Nelson计算方法定义分相折算系数（摩擦阻力乘积因子）ϕ，通过它将两相混合物的摩擦压降梯度同气相或液相单独在管内流动时的摩擦压降梯度联系起来）分气相折算系数：

（3）

为假定气相单独流过流动截面A时的摩阻。分液相折算系数：

（4）

为假定液相单独流过流动截面A时的摩阻。134.2摩擦压降计算1-摩擦压降计算的经验方法分相流模型①-Lockhart—Martnelli和Martinelli—Nelson计算方法其中：（5）（因为

,所以

）（6）144.2摩擦压降计算1-摩擦压降计算的经验方法分相流模型①-Lockhart—Martnelli和Martinelli—Nelson计算方法将（5）、（1）代入（3）式中得：

（7）（较精确的，没有不合理假设）

都与各自的Re数有关，并按通用的勃拉修斯公式计算有：

（8）

（注：n值取决于流态）其中:Re=2000-105时，n=0.25，

（水力光滑管）；或Re=5000-200000时，n=0.2，

）

（9）（光滑管区的Blasius公式为：

）154.2摩擦压降计算1-摩擦压降计算的经验方法分相流模型①-Lockhart—Martnelli和Martinelli—Nelson计算方法又因

则（

所以，

（10）（其n存在的前提是假定气相流过AG时的流态和其单独流过管道A时的流态相同，因而才认为式（8），（9）中的t，n相等）同理可得：

（11）又定义：

（12）X2—称为Lockhart-Martinelli参数

所以，

（13）

（14）16再将式（14）代入式（10）、（11）中得

（15）

（16）由上述结果表明，两相流的分相折算系数可以用参数X加以整理，可通过查图确定。而且对于一定的n值（主要取决于流态），

（或

）

只是参数X的函数171-摩擦压降计算的经验方法分相流模型①-Lockhart—Martnelli和Martinelli—Nelson计算方法4.2摩擦压降计算流态液相Re气相Re紊流（液）—紊流（气）（tt）>1000>1000紊流—粘性（tv）>1000≦1000粘性—紊流（vt）≦1000>1000粘性—粘性（vv）≦1000≦1000注:把Re数等于1000作为层流、紊流的界限，是因为上述Re中的速度是折算速度，因此一种相的Re数有效值会因为另一相的存在而增大，而实际的各相雷诺数都要大。184.2摩擦压降计算1-摩擦压降计算的经验方法分相流模型①-Lockhart—Martnelli和Martinelli—Nelson计算方法与X的关系如图所示194.2摩擦压降计算1-摩擦压降计算的经验方法分相流模型①-Lockhart—MartInelli和Martinelli—Nelson计算方法计算方法和适用原则在Lockhart—Martinelli方法计算

时，先计算出

，求得X值。然后再在线算图上（或计算）求得

或者

，最后再由

或

求得

。该方法适用于双组分的气液两相流在低压时摩擦阻力计算，因其计算数据是建立在低压的气液流动基础之上的。为了适用于汽—水混合物的摩擦阻力压力降的计算，Martinelli-Nelson对此方法进行了改进。204.2摩擦压降计算1-摩擦压降计算的经验方法分相流模型①-Lockhart—Martnelli和Martinelli—Nelson计算方法maRtinelli-nelson改进方法假设汽水两相分开流动时都呈紊流状态，同时利用常压下的空气—水混合物试验数据和高压汽—水混合物的试验数据建立了

的关系曲线。定义

，全液相折算系数－两相压降与假设汽水混合物全部为液相时的摩擦压力降（以总的质量流向相等为原则）之比

为假设两相混合物全部为液相时的摩擦阻力系数214.2摩擦压降计算1-摩擦压降计算的经验方法分相流模型①-Lockhart—Martnelli和Martinelli—Nelson计算方法mastinelli-nelson改进与x，P的关系如下图A图B224.2摩擦压降计算1-摩擦压降计算的经验方法分相流模型①-Lockhart—Martnelli和Martinelli—Nelson计算方法mastinelli-nelson改进按图A根据干度x及压力P查出

，再乘上算出的

值，即可得出不受热管中的

值。对于进口处干度x=0，出口处x=xe的受热管，可按出口干度xe及P值在图B中查出自x=0到x=xe的

平均值

，然后按同法求得

值。一般认为此方法适用于低质量流速工况。。234.2摩擦压降计算1-摩擦压降计算的经验方法分相流模型②-奇斯霍姆（Chisholm）计算法奇斯霍姆对两相流动摩擦阻力压降进行了许多研究，得出了

与X值的关系为或者

其中

式中C为系数，可按下式确定：

式中

—系数C2—系数

、

—气体及液体的比容，m3/kg.

系数C是压力P和干度x以及质量流速的函数244.2摩擦压降计算1-摩擦压降计算的经验方法分相流模型Chisholm推荐C值如下②-奇斯霍姆（Chisholm）计算法LiquidGasSubscriptCTurbulentViscousTurbulentViscousTurbulentTurbulentViscousViscousttvttvvv2012105对压力P>3MPa的汽水混合物，C值的计算A:当质量流速

时对光滑管：

C值按P31页公式计算，此时

。对粗糙管：

此时

。254.2摩擦压降计算1-摩擦压降计算的经验方法分相流模型②-奇斯霍姆（Chisholm）计算法对压力P>3MPa的汽水混合物，C值的计算B:当质量光流速

时对光滑管（

和粗糙管

式中

，

，

值可按下法计算；

C值仍然按照P31公式计算。对于粗糙管

；对光滑管

。对于粗糙管n=0，对于光滑管n=0.25.由上所述可见，奇斯霍姆计算法的计算过程较麻烦。但是计算结果和试验偏差较小。264.2摩擦压降计算1-摩擦压降计算的经验方法分相流模型③-Friedel（1979）计算方法（分相模型）Friedel把一个有25000个点的数据库用他本人的关系式以及当时存在的其他公式进行了比较，得出全液相折算系数其中

—分别为在总质量流速相同的情况下，具有气相或液相物性的单相流体在管内流动下的摩擦阻力系数，即全气相、全液相时的摩擦阻力系数）

适用于垂直向上流动与水平流动（应用于垂直向下流动的关系式稍有不同，其对单组分流动的标准偏差约为30%，对双组分流动标准偏差约为40—50%，可见误差仍然相当大）其中假设两相混合物的总质量全部以液体流过时的压降274.2摩擦压降计算1-摩擦压降计算的经验方法分相流模型当

时，应用Friedel关系式当

时，应用Chisholm关系式当

时，应采用Martinelli关系式在“传热和流动服务中心”的专利数据库基础上得到的计算结果对以上几个分相模型公式的应用建议：（通过试验比较给出的建议）284.2摩擦压降计算1-摩擦压降计算的经验方法我国摩擦压力降的计算方法①林宗虎—林宗振计算法作者等应用能谱理论和混合长度理论相结合的方法导得一种计算管内气液两相流动的摩擦阻力压力降的方法设有一绝热的均匀气液两相紊流沿水平管道流动。气液两相流沿管道中心线方向的速度为u，在垂直中心线方向的速度为v。令u’为速度u的脉动分量，v’为速度v的脉动分量，ρ’为密度ρ的脉动分量。则可得到气液两相流的摩擦切应力如下： （4-1）式中

，

—均为统计平均值；y为从管子内壁算起的距离。由混合长度理论可知：

（4-2）

（4-3）式中

，

分别为速度分布及密度分布的混合长度。294.2摩擦压降计算1-摩擦压降计算的经验方法我国摩擦压力降的计算方法①林宗虎—林宗振计算法

根据能谱理论，前述两种速度脉动分量之间的关系为：

（4-4）式中k—紊流流动的外尺寸；

—无穷大

—谱密度谱密度可用下式表示：

（4-5）式中s——任意值；b—常数；

—总能量；v—运动黏度。将式（4-5）代入（4-4），可得

（4-6）304.2摩擦压降计算1-摩擦压降计算的经验方法我国摩擦压力降的计算方法①林宗虎—林宗振计算法式中令

。将公式（4-6）代入（4-1），并令

，

可得

（4-7）速度u和密度ρ之间的关系可用下式表明：

（4-8）

上式中的A’和B’均为常数。应用下式边界条件：y=0处（管壁上）

，

处（管子中心），

，并将上述条件代入式（4-8），可以得到下列密度和速度之间的关系式：

（4-9）

上式中

及

分别表示液体密度及管子中心处气液两相流密度。314.2摩擦压降计算1-摩擦压降计算的经验方法我国摩擦压力降的计算方法①林宗虎—林宗振计算法令

（4-10）

则上式（4-9）可改写成为下列形式：

（4-11）

（4-12）设两相流动混合长度和单相流动的相同，则可应用下式计算混合长度

（4-13）式中K1及K2均为无因次混合长度常数。联立解式（4-7）、（4-11）、（4-12）及（4-13），并令

（4-14）

可得

（4-15）

式中

并

令

324.2摩擦压降计算1-摩擦压降计算的经验方法我国摩擦压力降的计算方法①林宗虎—林宗振计算法进一步推导可得

（4-16）式（4-16）的边界条件为

。应用这一边界条件，且在a，K1，

和

给定时，通过逐步积分可得出

和y之间的函数关系。对于管子半径为R的圆管，平均混合物物密度

可由下式求得：

（4-17）

式（4-17）也可以改写成

（4-18）

式中

（4-19）式中ReL—按总流量和液体物性计算的雷诺数334.2摩擦压降计算1-摩擦压降计算的经验方法我国摩擦压力降的计算方法①林宗虎—林宗振计算法

（4-20）质量流速

（4-21）速度分布

（4-22）

将式（4-22）及前述的

等计算式代入式（4-21），可得

（4-23）或改写成

（4-24）当总流量为液体时，摩擦切应力可按下式计算：

（4-25）两相摩擦阻力压力降

和单相液体摩擦阻力压力降

之比可写作为

（4-26）

或者

（4-27）当常数K1，

和

已知时，应用式（4-16）、（4-18）（4-19）和式（4-27）可得出两相摩擦阻力压力降

的解。莱维（Levy）的研究表明：K1=0.4，

=26，

。因而根据上述4个方程式可解得

和

的值如下图所示344.2摩擦压降计算1-摩擦压降计算的经验方法我国摩擦压力降的计算方法①林宗虎—林宗振计算法354.2摩擦压降计算1-摩擦压降计算的经验方法我国摩擦压力降的计算方法①林宗虎—林宗振计算法曾以此法算得的气液两相摩擦阻力压力降值和多种不同压力下得出的实验值进行比较，比较表明，计算值和实验值相符很好，特别是在高压及超高压时。若知

可按

值在图上查得

值。

可按照多列查尔（Dolezal）计算法

算得，

为单相液体的摩擦阻力系数，可按式（1—8）计算,可求出

。

，其中

按林宗虎《气液两相流和沸腾传热》书中公式（4-28）计算，因而

（4-28）式中的截面含气率

值可按照按林宗虎《气液两相流和沸腾传热》书中第四章中式（4-4）计算.364.2摩擦压降计算1-摩擦压降计算的经验方法我国摩擦压力降的计算方法②陈立勋计算法陈立勋计算采用的型式和苏联61年锅炉水力计算标准方法建议使用的公式一样；

（4-29）只是校正系数ψ的值计算方法有所不同。陈立勋以林宗虎《气液两相流和沸腾传热》书中图3-2的试验曲线为基础得出的这些试验曲线校正系数的近似数学表达式如下：当

≦1500kg/m2·s时，

值按照下式计算

（4-30）

≧1500kg/m2·s时，

值按照下式计算

（4-31）式中x—蒸汽干度。按照此法计算结果与前述奇斯霍姆计算式的计算结果相近，但是计算过程要简便得多。374.2摩擦压降计算1-摩擦压降计算的经验方法我国摩擦压力降的计算方法③我国电站锅炉水动力计算方法我国电站锅炉水动力计算方法中计算汽—水混合物两相摩擦阻力压力降的计算式的计算型式和式（4-29）型式相同，但是校正系数值的计算方法主要采用（4-30）和式（4-31）的型式。根据林宗虎的研究，如果将式（4-30）和式（4-31）所适用的质量流速分界值

改写成

，并将式（4-30）和式（4-31）中的1500改写成1000，则计算所得的摩擦阻力压力降值和其他实验值以及运行数据吻合的更好。压力为4.5-10.5MPa，质量流速为500-2700Kg/m2·s，热流密度为0-0.57MW/m2,进口干度为0-0.81的试验条件所进行的水平管中汽—水混合物摩擦阻力压力降试验证实了这一点。384.2摩擦压降计算1-摩擦压降计算的经验方法我国摩擦压力降的计算方法③我国电站锅炉水动力计算方法

在我国电站锅炉水动力计算方法中，汽—水混合物在水平、垂直以及倾斜管中流动时，摩擦阻力压力降按下式计算：

（4-32）

式中

-汽-水混合物的质量流速，kg/m2·s，也可用

代替，因为在等直径直管中各截面的质量流速相等；X—管内平均质量含气率；

λ

—摩擦阻力系数394.2摩擦压降计算1-摩擦压降计算的经验方法我国摩擦压力降的计算方法③我国电站锅炉水动力计算方法摩擦压降校正系数按以下方法计算：

（1）

；（2）

时，

（4-33）（3）

时，

（4-34）404.2摩擦压降计算1-摩擦压降计算的经验方法我国摩擦压力降的计算方法③我国电站锅炉水动力计算方法ψ值也可以按照线算图查得。对于不受热管，ψ

值按照林宗虎《气液两相流和沸腾传热》书中图3-12查得；对于受热管

ψ值按林一书上图3-13查得。对于管子入口处工质已经是汽—水混合物的受热管，当出口干度和进口干度之差xe-xi

≧0.1时，值按下式求得：

（4-35）

式中

及

分别按xe及xi在林宗虎《气液两相流和沸腾传热》书中图3-13查出。414.2摩擦压降计算1-摩擦压降计算的经验方法我国摩擦压力降的计算方法③我国电站锅炉水动力计算方法此法经国内各主要锅炉厂对各种类型锅炉的水力工况试算后，和实际工况接近。此法宜用于压力大于1MPa的锅炉，用于低压时，误差增大。当

xe-xi

≦0.1时，

值按下式计算：

（4-35）式中，

值为按

查出的

值，

及

则分别按xe及xi在林宗虎《气液两相流和沸腾传热》书中图3-12中查出。424.2摩擦压降计算2-依据流型的摩擦压降的计算对气液两相流摩擦阻力压力降的计算至今已做了不少研究工作，提出了许多纯经验的或半理论半经验的计算式或计算方法。但大多数都存在着适用范围较窄和考虑影响因素不够全面的缺点。对基于某种流型的具体特点，建立在理论分析基础上的摩擦阻力压力降计算方法，也曾进行了一些研究。Mandhane曾经收集了近万个气液两相流体在水平管中流动的试验数据，并和16个摩擦阻力压力降计算式的结果进行了比较，表明根据管中流型应用相应计算式的计算结果最精确。本节就将根据不同流型的具体特点，分别建立相应的摩擦阻力压力降计算式。4.2摩擦压降计算在泡状流中气相弥散在连续的液相中，这时摩擦阻力的产生主要是由以速度v运动的混合物中的液相造成的，因而泡状流的摩擦阻力压力降可表示为：

式中fM为两相混合的摩擦阻力系数，ρM为两相混合物密度，G为混合物的质量流速.

①

泡状流计算2-依据流型的摩擦压降的计算

Beattie&Whalley通过与大量的水平或垂直管中的压降数据的分析比较认为，对泡状流如fM、ρM按以下方法计算，将获得较好的结果：

式中ε/D为管壁相对糙度，ReM为两相混合物的雷诺数。

4.2摩擦压降计算①

泡状流计算2-依据流型的摩擦压降的计算混合物密度ρM按均相模型计算如下;

两相混合物的动力粘度μM可由下式计算;

4.2摩擦压降计算①泡状流计算2-依据流型的摩擦压降的计算间歇流中，每一个流动单元的摩擦阻力压降由两部分组成：液弹段摩擦阻力压力降和Taylor气泡段摩擦阻力压力降，即：

4.2摩擦压降计算②间歇流计算2-依据流型的摩擦压降的计算式中fs，fG分别为液弹的摩擦阻力系数和Taylor气泡段的摩擦阻力系数DG为Taylor气泡的当量直径fs，fG可分别根据ReMS，由Chen式求得

其中，液弹段混合物密度液弹段混合物粘度 4.2摩擦压降计算②间歇流计算2-依据流型的摩擦压降的计算Taylor气泡的当量直径DG对垂直管为：

DG＝D－2δL

δL为Taylor气泡段液膜平均厚度，

因此 4.2摩擦压降计算②间歇流计算2-依据流型的摩擦压降的计算对水平或微倾斜管;

SG,Si可由εLF按以下方法求得;由于 于是 可由该式中求得θ值则无量纲液位高度

4.2摩擦压降计算②间歇流计算2-依据流型的摩擦压降的计算

在式中， 间歇流的平均摩擦阻力压力降梯度为：

在以下计算中，所有的间歇流特性参数，如εLS、εLF、LS、LU、VT均等可由第四章中相应的间歇流特行参数的封闭方程组中求得。4.2摩擦压降计算②间歇流计算2-依据流型的摩擦压降的计算倾角为β的气芯和液膜的动量方程为：

4.2摩擦压降计算③环状流计算2-依据流型的摩擦压降的计算以上两式相加消去τiSi，并整理得：

又因AL/A=1-α且ρM=αρG+(1+α)ρL,于是上式可写成:

总的压力降由两部分组成组成，即摩擦阻力压力降（右式第一项）和重位压力降（右边第二项），加速压力降在上面的模型中被忽略。于是环状流的摩擦阻力压力降为：

τB为液膜与管壁的摩擦切应力，SB为液膜湿周,A为管道流通截面积。4.2摩擦压降计算③环状流计算2-依据流型的摩擦压降的计算

将以上三式代入（5-30）中得：

液膜摩擦阻力系数 式中DL为液膜的当量直径，由第（3.3.3）节中的分析知：

于是

当ReLF<2100时，CL=16,n=1当ReLF≥2100时，CL＝0.0791，n＝0.254.2摩擦压降计算③环状流计算2-依据流型的摩擦压降的计算在水平或微倾斜管中，由于重力的影响，在两相流过程中，常常会出现分层流，对分层流的流动结构分析，已有叙述，液相和气相的动量方程分别为：

4.2摩擦压降计算④分层流计算2-依据流型的摩擦压降的计算以上两式相加，消去τiSi项，并整理得：

式中，气相和液相的摩擦切应力τG,,τL

可分别计算如下：

fL、fG可分别按ReL、ReG求得。4.2摩擦压降计算④分层流计算2-依据流型的摩擦压降的计算又AG/A=1-α,AL/A=1-α,ρM=αρG+(1+α)ρL,SL=SL/D,SG=SG/D,连同式（5-41），（5-42）一起代入式（5-40）中可得：

上式中第一项就是分层流的摩擦阻力压力降，即

4.2摩擦压降计算④分层流计算2-依据流型的摩擦压降的计算无量纲湿周SG,SL可由截面含气率α计算如下：先由 通过迭代求得湿角θ（弧度），然后有：SG=θ/2,SL=π-θ/24.2摩擦压降计算④分层流计算2-依据流型的摩擦压降的计算4.2摩擦压降计算3-影响气液两相流体摩擦阻力压力降的主要因素压力P的影响可参见林宗虎《气液两相流和沸腾传热》中的3-15.由于压力增高，气液密度比增大，因而

下降。当压力为临界压力时，

等于1.0。干度x的影响也可以参见林宗虎《气液两相流和沸腾传热》书中图3—15.一般的规律是在管子不受热、当x增加时，

先增大后减小。在受热管中，随着x的增加，

先增大，后减小，然后增大。这主要是由于干度变化和受热时流型变化造成的。①压力P及干度x的影响现在的有关摩擦阻力压力降的试验资料大多是对管子内径小于50mm的工况的，在此范围内，管径对摩擦阻力压力降的影响不大，可以略去不计,但对于大管径用这些方法误差就很大。目前对这方面的研究方兴未艾（由于海上油气开采的进行），建立依据流型的压力降计算方法并同时考虑管径等因素的影响（实际上式管径的不同造成两相流型转变的变化，从而引起摩阻计算公式的不适用）。②管径的影响4.2摩擦压降计算3-影响气液两相流体摩擦阻力压力降的主要因素热流密度q对

的影响已经有不少研究者进行过试验，但由于试验时的热流密度值不够高，一般低于85kw/m2,所以q对

的影响显示不出来。

③管子热流密度的影响泰勒苏娃的试验证明：若其他条件相同，在干度x小于某值时，受热管的

值比不受热管的大，两者的差别随x的增大先增加后减小，当x=100%时，两者差别消失。x较小时，比如是泡状流时，对受热管由于其在截面也产生气泡，而导致对中心流的阻力增加且使得截面上的含气率增加，因贴壁上的气泡流动较慢，在总x不变时，α就增加，受热管的摩阻比不受热管的大；但是，当干度较高时，如果是受热管就会产生更多的汽相，使得截面上的含气率相对于不受热管增大（保持x一定，而贴壁的汽体流动较慢，所以α）增大，两相混合物就更近乎单相汽时的流动，所以这是受热管的摩阻大而不受热管的小。4.2摩擦压降计算3-影响气液两相流体摩擦阻力压力降的主要因素质量流速对摩擦阻力压降均有影响（无论是什么布置方式）在某些参数范围内，这种影响是相当大的。比如Osbin（等）的试验结果。⑤质量流速的影响因流动方向不同时，在相同流动条件下的流型也不一样，因而流动时压力降也不一样，这时实际上是流型变化导致的影响。其他一些试验研究证实：（1）M—N法比较适用于低质量流速范围。（2）在高质量流速范围（G>2000—2500kg/m2·s）采用均相模型较为适宜。（两相湍动度大，两相分布更为均匀）⑥流动方向β4.2摩擦压降计算3-影响气液两相流体摩擦阻力压力降的主要因素

研究较少，通常管壁粗糙物突出贴壁液膜较多时，两相流体流动沿粗糙管流动时的摩擦阻力系数约比沿光滑管流动时大一倍。但当液膜能盖住管壁粗糙度的凸尖时，则光滑管和粗糙管的摩擦阻力压力降相同。

不同的流动型态造成两相流动规律的不同，因此根据不同的流动形态的具体特点建立相应的数学物理模型，从而求得相应的摩擦阻力压力降，在本质上能接近真实的流动状态。6.管内粗糙度

4.2摩擦压降计算3-影响气液两相流体摩擦阻力压力降的主要因素该方法也是目前正在广泛研究的课题之一。无非是

等的影响，归结起来很大程度上都是由于流型的变化，而使得摩擦阻力发生变化。统一的经验的方法是需要综合考虑各种因素。显然是力不从心。建立依据流型的压力降计算模型和准确确立管内的流型的流型的转换界限是紧密相联的。目前这方法研究主要有：Bendiksen（SlugFlow），陈宣政（垂直上升管），王树众(各种布置状况下)、Kokal和Stanislav（1989）水平管等。4.2摩擦压降计算作业直径D=5.08cm管子，P=180bar，进口流量M=2.14kg/s,进口为饱和水，

粗糙管

，出口干度

，管长100m，求

。

分别用M—N法、Chisholm方法（经验的C公式）、苏联78年计算标准、我国水动力计算方法。（参考：陈之航、曹柏林、赵在三，“气液双相流流动和传热”，机械工业出版社）§4.3加速压降

在实际工程应用中，在多相管流的压降梯度中，加速压降所占的份额很低。某些研究者导出多相管流压降关系式时，考虑了加速压降，有的则忽略加速压降。本章对加速压降的计算和特点进行了讨论。气液两相流体在管路中流动时，加速压力降的产生通常由以下两个原因造成：（1）由于加热、冷却以及压力变化，使两相混合物的组成或流速发生变化；（2）由于管路流通截面积发生变化，从而引起两相流速发生变化。当两相混合物在等截面直管中流动时，第二部分的加速压力降就消失了。在通常的气液两相流操作速度范围内，由于两相流速较低，加速压力降常忽略不计。当两相混合物流速较高或两相混合物由于压力、温度等的变化而凝析或气化，造成两相流速发生较大变化时，就可能会导致较大的加速压力降。这时，必须计及加速压力降的影响。对加速压力降的计算基本上有两种方法:1、基于气液两相完全分开流动的分相模型进行；2、根据均相模型进行推导。由于后者不考虑两相间的相对滑动，因而误差较大。本节将采用分相模型计算方法来导得气液两相流的加速压力降。由两相流的动量方程知稳定流动的气液两相流的加速压力降为:

式中，G为两相混合物的质量流速，流动过程中质量守恒。在上式中:第一项为由于两相混合物组成（x）或压力变化（P）而引起的压力降;第二项为由于管路流通截面积发生变化而造成的加速压力降。若在等截面直管中，气液两相混合物从位置Z1流到Z2，则上式中第二项为零，于是加速度压力降梯度为：

§4.4重位压力降管路内气液两相流的重位压力降为

其中，β为流动方向和水平线的夹角，向上流动时为正。a为相应流型下的截面含气率。在间歇流中，由于每个流动单元总是由气泡段和液弹段两部分组成，环绕Taylor气泡的是贴壁液膜，而Taylor气泡内压力沿长度方向基本上是恒定的，因而Taylor气泡造成的重位压力降可以忽略不计，这样段塞流的重位压力梯度可表示为：式中，ρMS为液弹中两相混合物的密度，可按下式计算：

试验表明，重位压力降是管内气液两相流竖直流动时总压降的重要组成部分，因此准确计算重位压力降，对提高管路内两相流压力降的计算精度具有重要意义。而重位压力降计算的关键就在于对截面含气率（或持液率）的准确计算。§4.5局部阻力目前，对气液两相流流过局部阻力件时，局部阻力压力降的研究还很不够，还主要是依据试验确定的经验的方法。通常通过人工查表求得各种经验系数，然后再进行计算。下面重点介绍弯头、节流式测量元件、突扩、突缩接头、三通、阀门的局部阻力压力降。§4.5.1弯头的局部阻力压力降汽液两相流通过弯头时的阻力可分为两部分：在弯头内部产生的，是由于两相流体通过弯头时出现二次流和流场变化而引起的；由于汽液两相滑动比发生了变化，在下游又要恢复到直管内的滑动比所引起的动量变化。

Chisholm提出了一种半理论半经验的方法来计算流过弯头的局部阻力压力降。对第一部分压力降按均相模型有：

PBL0为假设两相流体全部为液体时，流过弯头的摩擦压力降：

式中,λ－单相流体摩擦阻力系数L/D－弯头当量长度G－两相混合物的质量流速，kg/(m2·s)

对单相流体的摩擦阻力系数，Churchill推荐在整个流动范围内（在全Re数范围内）可按下式计算：

其中

B=(37530/Re)16弯头当量长度L／D，对于90°弯头可根据弯头相对弯曲半径R／D求得。通过数据拟合L／D和R／D的关系如下：对第二部分压力降按分相模型计算：

式中MFL0为假定两相混合物全部为液体时单位面积的动量MFL0

=G2/ρL系数B反映了流过弯头过程中，出口处和进口处的滑动比的变化

于是两相流体流过弯头时的局部阻力为：

经整理得

式中系数B′对90°弯头有：

当弯曲角度小于90o时，取B′＝B90′当弯曲角度大于90o时，其中——单相液体流过90度弯头的阻力系数

——单相液体时，弯曲角度为的弯头阻力系数

其中ξ0－弯头的原始阻力系数，与其形状和相对弯曲半径有关；k－考虑管壁粗糙度的影响系数；

(kξ0)对缓转弯与焊接弯头可按下式给出（经数据拟合）：B－是与弯曲角度θ有关的修正系数，经数据拟合可按下式计算：

式中θ的单位为度°

应该注意的是，上式仅用于质量含气率x<0.04的情形。因为在x<0.04时，该式是有一定的试验基础的。对x>0.04时弯头的局部阻力计算如下方法:

式中,PL0－假定两相全部为液体时，流过弯头的局部阻力压力降,

-单相液体流过弯曲角度为θ

的弯头时的局部阻力系数；α－相对弯曲半径R/D的修正系数，可由下式查出经多项式数据拟合可按下式计算R/D≤34567α00.71.31.82.0§4.5.2节流式测量元件的局部阻力压力降

孔板、文丘利管和喷嘴是节流式测量元件，可以测流量和干度。因为它们结构简单、无转动部件、使用简便、运行可靠且有足够多的精度，所以主要用作两相测量设备。它们的工作原理相同，下面以孔板为例。1）对单相流体：(kg/s)

Y—流体膨胀系数，若认为流体不可压缩，则=1

—孔板的膨胀系数

—孔口直径和管子内径之比，

A—孔口截面积，

c—流出系数，dischargecoefficient，其值与取压方式（角接取压、法兰取压、径距取压等）、Re数等有关

—孔板前的流体密度

阻力损失

—孔板前后压力差

动能