气液两相流PPT课件.ppt

《气液两相流PPT课件.ppt》由会员分享，可在线阅读，更多相关《气液两相流PPT课件.ppt（145页珍藏版）》请在三一文库上搜索。

1、1,气液两相流,主 讲： 谷 海 峰 学 院：核科学与技术学院 Email ： 2010.03,2,课程简介,学时：40 学时 学分：2.5 性质：专业基础课,反应堆热工水力、核动力设备,流体力学 传热学 工程热力学,两相流动 与传热,3,考核方式：,闭卷答题+平时成绩,4,参考教材：,两相流与沸腾传热，徐济鋆编，原子能出版社。 两相流与传热D.巴特沃思等编，陈学俊，陈听宽，曹柏林译，原子能出版社。 两相流与沸腾传热，鲁钟琪编，清华大学出版社。,5,课程目录,第一章 两相流基本参数及其计算方法（4学时） 第二章 两相流的流型和流型图（6学时） 第三章 两相流的基本方程（4学时） 第四章 截面含

2、气率的计算（8学时） 第五章 直管的两相流压降计算（10学时） 第六章 两相流局部压降计算（2学时） 第七章 两相临界流动（4学时） 第八章 两相流流动不稳定性（2学时）,6,1.1 基本概念 1.2 气相介质含量 1.3 两相流的流量和流速 1.4 两相介质密度及比容,第一章 两相流基本参数及其 计算方法,7,1.1 基本概念,1.物态：在某一条件下，物质存在的一种状态。常见的物态是气态、液态和固态。有时物态也称之为相，常见的物质三态也称为：气相、液相、固相。,判断：下列体系哪些是单相的？哪些是两相的？ 1、水和水蒸汽 4、有沙粒的小溪 2、盐水溶液 5、油水混合物 3、氢气和氮气 6、水和

3、冰,8,1.1 基本概念,2.相：通常指某一系统中具有相同成份且物理、化学性质完全均匀部分，各相之间有明显的界面。 3.辨别单相体系与两相体系 是否系统内各部分的性质均匀 是否存在明显的相间界面 4.两相流：任意两个相组合在一起、且具有相间界面的流动体系。,9,10,1.1.1 两相流的分类,（1）按化学成份分为 单组分两相流：同一种化学成分的物质的两种相 态混合在一起的流动，如水-水蒸汽、钠-钠蒸汽及冰 水混合物等。 双组分两相流：不同化学成分的两种物质处于同 一个系统内的流动，如空气-水，油-水，烟-气等。 （2）按流道是否存在热交换可分为 绝热两相流无相变，无相间质量交换。如汽水 分离。

4、 加热两相流有相变，有相间质量交换。如沸 腾，冷凝。,11,1.1.1 两相流的分类,（3）按两相物质所处的物态可分为 气液两相流：水和水蒸汽，水和空气 气固两相流：风沙，烟气 液固两相流：血液流动 液液两相流：两种不相溶的液体，油水混合物,12,1.2.1 定义 气相介质含量表示两相流中气相所占的份额。 1.2.2 几种表示方式 1.质量含气率x 单位时间内，流过通道某一截面的两相流体总质量M中气相所占的比例份额。 式中， 、 分别表示气相和液相的质量流量，kg/s。 那么，质量含液率(湿度)可以表示为,1.2 气相介质含量,13,又称为热平衡含汽率，它是由热平衡方程定义的含汽率。在有热量输

5、入的两相流系统中，可以根据输入的热量得到汽相的含量。 （1）热力学平衡（thermodynamic equilibrium) 系统内即无压力梯度和温度梯度，且该系统内所有共存相内也无化学势梯度时的状态。满足力学平衡、热平衡和化学平衡。 （2）热平衡方程 式中，i是流道某截面上两相流体的焓值，那么,2.热力学含汽率x （干度）,汽化潜热,饱和水焓,14,（3）讨论 当 ，两相流动处于欠饱和态， 过冷（欠热）沸腾或单相； 当 ，两相流动处于饱和 沸腾，这时热平衡含汽率等于质量含汽率； 当 ，过热蒸汽，单相汽体。,15,指单位时间内流过某一截面的两相总容积中，气 相所占的比例份额。 式中， ， 分别

6、表示气相和液相介质的容积流量， m3/s。那么，容积含液率（1-） 由定义，可以推出质量含气率x与的关系,3.容积含气率,16,4.截面含气率,又称为空泡份额。指两相流某一截面上，气相 所占截面与总流道截面之比。 式中， ， 分别表示气相和液相所占的流道 截面积，m2。同理，截面含液率（1-),17,是容积流量比，表示流过通道的气相容积份额； 是体积比，表示存在于通道中的气相容积份额。,（1-12）,（1-13）,18,显然，对于单组分的两相流，,在两相流动系统中，两相之间会因物性不同(如密 度)存在不同程度的滑动。引入滑速比S的概念，则两相 之间滑动的大小可以表示为,19,1.3 两相流的流

7、量和流速,1.3.1 几种表达式形式 1.质量流量和质量流速 总质量流量M：单位时间内流过任一横截面的气液 混合物的总质量,kg/s. 总质量流速G：又称为质量流密度，指流道单位截 面通过的质量流量,kg/m2.s. 每一相的质量流速与总质量流速的关系,20,2.容积流量、相速度和折算速度,总容积流量：单位时间内，流过通道任一流通截面的气液混合物的总的容积，m3/s。 液相容积流量， m3/s ： 气相容积流量, m3/s ：,21,液相真实平均速度，m/s： 气相真实平均速度, m/s：,2.容积流量、相速度和折算速度,液相：,气相：,22,折算速度：又称容积流密度，又称为表观质量 流速（s

8、uperficial flow flux）, 定义为单位流道 截面上的两相流容积流量，m/s。它也表示两相流的平 均速度。 式中，Jg为气相折算速度，表示两相介质中气相单独流 过同一通道时的速度,m/s. Jf为液相折算速度，表示两相介质中液相单独流 过同一通道时的速度,m/s.,23,讨论,当气液两相无相对运动时 当气液两相存在相对运动，且 ，则 当气液两相存在相对运动，且 ，则,24,3.漂移速度和漂移通量,漂移速度：各相真实速度与两相混合平均速度J的差值。 气相漂移速度： 液相漂移速度： 漂移通量：各相相对于两相混合平均速度J运动的截面 所流过的体积通量。 气相漂移通量： 液相漂移通量：

9、,25,循环速度：与两相混合物总质量流量M相等的液相 介质流过同一截面的通道时的速度。 循环倍率：单位时间内，流过通道某一截面的两 相介质总质量与其中气相质量之比，也就是质量含气率 的倒数。,4.循环速度和循环倍率,26,自然循环锅炉原理,对于沸腾通道，循环速度在 数值上等于上升管入口处水 的速度 。 循环倍率是指循环水流量与 蒸汽产量的比值。表示上升 管中实际产生一公斤蒸汽需 要的循环水流量几公斤,27,对于自然循环蒸汽发生器，循环流速是指上升空间通道入口处饱合水的流速。 循环倍率是指循环水流量与蒸汽产量的比值。,立式自然循环蒸汽发生器,28,对于蒸汽发生器而言，其循环倍率是表征其二次侧流动

10、状 态的重要参数，它对于传热管的腐蚀、流动稳定、传热特性及 分离器工作等都有重要影响。一般认为在设计状态时应大于44.5，其主要考虑的因素如下： 1.传热要求：为了保证管壁润湿，特别是防止局部区域出现缺液或干涸。一般要求管束进出口处的蒸汽含量不超过20%25%，相当于循环倍率大于4。 2.流动稳定性：循环倍率低可能导致流动的不稳定，使流动产生振荡，这种流动振荡现象使传热效率下降，当流动振荡的幅度足够大时，就可能引起水和蒸汽流量的大幅度波动。 3.管材腐蚀：在局部滞流或低流速区，往往导致污垢沉积或浓缩。从防止腐蚀与流动的要求出发，应适当提高循环流速，以便在管板上表面及管束弯管区提高冲刷流速。,2

11、9,5.滑移速度和滑速比,由于气体和液体的流速不同，在两相之间存在滑动速度(滑移速度） 气体的速度与液体的速度之比称为滑速比S,竖直上升流：,竖直下降流：,30,1.4 两相介质密度及比容,1.两相介质密度表示方法 1）两相介质的流动密度 单位时间内流过流道某一截面的两相介质质量和 体积之比。 2）两相介质的真实密度 即单位体积内两相介质的质量，反映了存在于流 道中的两相介质的实际密度。,31,和 的关系：,1).当S=1时, 所以,，均质流动；,2). 当,时，,S1:,则,S1:,则,32,单位时间内流过流道某一截面的两相介质体积和质量之比 1) 截面平均比容vA 2) 动量平均比容vM,

12、2. 两相介质的比容,33,3）动能平均比容vE,对比上面四种比容表达式,当S=1时，则有,，,。,34,两相流特性参数的分类 (1) 流动参数： 可根据物质平衡或热平衡关系计算得到的参数,(2) 真实流动特性参数： 代表流体流动时真实流动特性的参数，这类参数与两相流体的相对运动有关，有,35,本章小结,1. 何谓两相流？ 两相流的分类？ 2. 有关含气率、速度、密度、比容定义及计算式。 3. 两相流特性参数的分类 4. 何谓热平衡？ 5. 滑速比与各参数的关系,作业：,1-4,36,不考虑两相间的相对滑移时,考虑两相间的相对滑移时,37,第五章 直管的两相流压降计算,？,38,序,两相流流动

13、压降计算的重要性 两相流流动压降计算的局限性 本章主要内容 1. 两相流摩擦压降的计算模型； 2. 影响两相流摩擦压降的主要因素； 3. 重位压降和加速压降的计算； 4. 环状流动的压降计算。,39,两相摩擦压降,单相摩擦压降,折算因子,计算模型,均相流模型,分相流模型,方法、 经验公式,40,重点 难点,41,5.2 均相流模型的摩擦压降计算,一.均相流模型计算法,两相摩擦压力梯度 对于圆管，控制体周界长度(m)： 通流面积(m2)： 流体与壁面的摩擦剪应力(N/m2)：,42,全液相摩擦压力梯度 与两相流总质量流量相同的液体流过通道时的压力梯度,全液相摩擦压降 梯度,单相水的摩阻系数一般按

14、布拉修斯（Blasius） 公式计算,43,全液相折算系数的求解,全液相折算系数,（1）令,44,两相平均粘度满足以下边界条件,（2）平均粘度法,把两相流体看作一种特殊的单相流体， 采用平均粘度来计算两相流体的等效摩阻系数,45,麦克达姆（Mecadam）计算式，应用最广,西克奇蒂（Chcchitti）计算式,德克勒（Dukler）计算式,代入（5-8）式,班可夫（Bankoff）计算式,46,计算 的步骤总结 1.求 (1) (2) 2.求 3.求,47,1.绝热、长度L的圆管内， 两相流动摩擦压降 2.非绝热 均匀加热 正弦加热,计算管长为L的一段圆管内的摩擦压降,48,例题：设油气混合物

15、在一内径为40mm的水平直管内流动，试按均相流模型计算两相摩擦压降梯度。 已知油气混合物的参数为：,解：（1）全液相摩擦压降梯度,全液相摩阻系数：,（2）全液相折算系数,（3）油气两相摩擦压降梯度,49,例：已知一均匀受热的试验段长为3.66m，管子内径为0.16mm，进口水温为204oC,压力为6.89MPa,试验段直立布置，管子进口水流量为0.108kg/s，加热功率为100KW，试用均相模型计算试验段内的摩擦阻力压降。 已知：压力为6.89MPa下 饱和水焓为1.26MJ/kg, 入口水焓值为0.872MJ/kg. 入口水比容为1.16510-3m3/kg, 饱和水的比容为1.3510-

16、3m3/kg， 蒸汽比容为2.7810-2m3/kg. 入口水的动力粘度为1.3510-4N.s/m2, 饱和水的动力粘度为0.97210-4N.s/m2, 蒸汽动力粘度为0.1910-4N.s/m2。,50,解：由于入口水的焓值低于饱和水的焓值，因此入口为欠饱和水。则试验段内的摩 擦阻力压降应由单相水的摩擦阻力压降和气液两相摩擦阻力压降两部分组成。,1.单相水的摩擦阻力压降,（1）根据热平衡方程，定出预热段长度，即单相水段的长度,由于沿管长均匀加热，因此,51,（2）单相水摩擦阻力系数,（3）单相水摩擦阻力压降,52,2. 汽水混合物管段的摩擦阻力压降,（1）全液相摩擦压降梯度,（2）全液相

17、折算系数,（3）气液两相摩擦阻力压降,53,3. 按均相流模型法计算得到的试验段内摩擦阻力压降,54,同样，也可以以气相摩擦压降计算为基础，计 算两相摩擦压力梯度。 全气相折算系数可定义为,全气相摩擦压降 梯度,55,二.前苏联锅炉水动力计算标准方法,K-绝对粗糙度；碳钢、珠光体钢K=0.08mm; 不锈钢K=0.01mm。,56,1950年,水平管：偏差20%，垂直管：偏差较大,用于计算气液两相流,1978年,1961年,57,(2)均匀加热 1)入口饱和水，出口xe1; 2)入口x1，出口x2。,查图5-1(a),绝热流动 x=const, 查图5-1(b),58,3) 的查图确定法 查图

18、的已知条件: x,P,G. 计算： 判别: 当P17.64MPa时，按右行路线在图右边中间 下边横坐标上查 值.,59,图5-1(a) 加热管,60,图5-1(b) 不加热管,61,四.我国电站锅炉水动力计算方法,摩擦压降修正系数，,62,摩阻系数修正系数 按以下进行计算,当 时： (2)当 时：,(3)当 时：,使用条件：,1)非绝热：,2)绝热两相流：x=const,63,5.3 分相流模型的摩擦压降计算,洛克哈特马蒂内里(L-M)关系式 1.L-M法的基本假设 (1)两相之间无相互作用，气相压降等于液相压降，且沿 管子径向不存在静压差； (2)液相所占管道体积与气相所占管道体积之和等于管

19、道 的总体积。 根据以上假设，各相的压降梯度彼此相等，且等于两相 流摩擦压降梯度。,液相部分压降梯度,气相部分压降梯度,64,液相部分的压降梯度可表示为 式中： 表示液相所占截面的当量直径，m； 表示液相部分摩阻系数。 气相部分的压降梯度可表示为 式中： 表示气相所占截面的当量直径，m； 表示气相部分摩阻系数。,65,2.摩擦压降梯度关系式,分液相折算系数 分气相折算系数,分液相摩擦压降梯度,分气相摩擦压降梯度,66,和 分别表示液相和气相单独流过同一 管道时的摩阻压降梯度。 表达式分别为 式中， 分液相摩阻系数; 分气相摩阻系数.,67,(1) 和 计算式的推导,同理可得:,68,洛克哈特-

20、马蒂内里参数,69,(2).分相流态的判别,分组原则是看各相单独流过相同管径管子时是 层流还是紊流而定。,层流-层流（ll）,紊流-层流（tl）,层流-紊流（lt）,紊流-紊流（tt）,70,(3).不同流态下Martinelli参数的计算,71,对于t-t工况：,72,对于l-l工况：,同理，可以求出,73,二.奇斯霍姆关系式,奇斯霍姆（Chisholm）把管内两相流摩擦压力 梯度表示成两相混合物“实际动压头”有关的表达式,74,其中,奇斯霍姆推荐的系数C值列在下表中。,75,奇斯霍姆法计算摩擦压降梯度的步骤：,计算分相Re数，判断组合工况，确定C值; （ll，lt, tl, tt ) 计算

21、分液相摩阻梯度和分气相摩阻梯度； 3. 计算马蒂内里参数,76,4.计算分液相折算系数或分气相折算系数 5.计算两相摩阻压降梯度,或,或,77,L-M方法和奇斯霍姆方法的应用条件,1、实验压力：0.110.35MPa 2、实验介质：空气和水、空气和油,78,上节内容回顾,L-M方法：,Chisholm方法：,79,三.马蒂内里-纳尔逊关系式（M-N）,1.基本假设 (1).流动工况：湍流-湍流（t-t）型； (2).流场中任意点均处于热力学平衡； (3).L-M法基本思想也适用于单组分两相流(沸腾工况) 2.摩擦压力梯度的计算 (1). 和 的计算,80,对于水工质，在临界压力Pcr=22.1

22、2Mpa时 两相流可作为单相流处理 此式为单组分两相混合物在Pcr下必须满足的条件。,(2). 压力效应的考虑,81,3. M-N法关系曲线,（1）绝热流动 a. P=Pcr时，取n=0.25,b. P=0.1MPa,取n=0.25，用L-M法（5-45）式,c. 中间压力，采用内插法决定。,82,83,（2）均匀加热,入口饱和水，出口xe1,对于汽水系统，其积分结果示于图5-4。,84,85,植田辰洋为M-N法提供拟合关系式,当 时，上式可用下列近似公式代替：,86,5.4 影响两相流摩擦压降的主要因素,一.压力P的影响 对于水，Pcr=22.12Mpa P=Pcr , PPcr ,87,二

23、.质量含气率x的影响,1.不加热管 2.加热管 出现马鞍形。,原因：在沸腾工况下，因为相变导致x沿途发生变化， 造成流型改变。比如：泡-弹-环。,88,三.质量流速的影响 1. 工程上在计算摩擦压降时认为 G2000 kg/m2.s , 采用均相流模型（如前苏联 水动力标准方法),巴罗塞方法：,89,90,91,2. 巴罗塞方法 优点: (1).适用于液态金属，如Na,K,F-22,汞等； (2).计算精度高。 缺点： 只有曲线，使用很不方便。 3. 巴罗塞奇斯霍姆拟合关系式 奇斯霍姆在研究巴罗塞曲线的基础上，给出了有关的拟合关系式。即,K值适用于 当G2000kg/m2.s,92,5.5 重

24、位压降计算,一.均相流模型的重位压降,1. 绝热流动,93,2.均匀加热,均相流模型主要用于低质量含气率x、高质量流速G的情 况。一些文献建议，只要符合下列条件之一，便可考虑采用均 相模型，即 但是，当液相粘度 时，建议不采用均相 模型。,X 沿途呈线性分布，当入口：饱和水；出口：Xe1,94,二. 分相流模型的重位压降计算,真实密度 与含气量和加热方式有关。 1.绝热流动 2.均匀加热,95,均相模型 1.绝热流动 2.加热流动 入口为饱和水 出口Xe1 入口为X1，出口X2,5.6 加速压降计算,96,1.入口X1,出口X2 2.入口饱和水，出口Xe1,二.分相流模型的加速压降,积分,97

25、,5.7 环状流动的压降计算,一. 气流中没有夹带液滴的环状流动 假设： (1) 两相完全分开 (2) 气液界面光滑或波动 1.连续性方程（ ） 气相： 液相: 气液两相总的容积流量：,98,2.动量方程 推导过程： (1)对于一微段，如果忽略进出口的动量变化，且为 一元流动 (2)沿流动方向，分别对气相和液相列出力平衡方程 气相： 液相： (3)简化、整理得出液膜内切应力分布方程,99,在管壁处：Y=0,管壁对流体的切应力为：,液膜内流体的切应力方程：,100,3.液膜内流速分布和液膜流量 (1) 液膜内速度分布方程 (2) 液膜流量,若,101,(3)液膜厚度 (4)截面含气率,102,4

26、.环状流摩阻压降梯度 根据界面状况可分为界面光滑和界面波动两种情况 (1) 环状流摩阻压降梯度推导过程 1)以气流核心为研究对象，给出气流核心摩阻梯度 2)根据L-M法计算两相摩阻梯度公式 3)给出分气相折算系数表达式,为环状流分气相折算系数与截面含气率的一般关系式，也适用于夹带环状流,103,的计算 气液界面上的剪应力来源于气流与液膜的速度差： 分气相剪应力： 通常，,104,的计算 a.光滑环状流 气液界面比较光滑时，则 适用于分液相折算系数 ，当液膜流速 较高时，气液界面出现波动，i将显著高于按光滑管估 计的数值，因此，按照上式计算得到的 偏低。,105,b.波动环状流 当液相速度较高时

27、，由于液膜的紊流流动引起了 液膜脉动振动，从而出现波动界面，此时 Wallis建议,106,二.雾环状流动 1.假设： (1)气流核心为气相和液滴的均匀混合物 (2)液膜速度呈线性，界面上切应力来源于气流与液 膜的速度差。 2.分析雾环状流 与 的关系,107,引入窜流比E概念 窜流比（夹带率）=气流中夹带的液相流量/液相总流量,夹带液相流量：,液膜流量：,108,气流核心流量：,气流核心混合物密度：,低压下，,109,110,3. 窜流比E的确定 (1)影响E的因素 液膜流速 通常用Rel来描述液相流量对液膜波动的影响 气相流速 通常用无因次量 来表征气流对液膜卷吸作用 (2)wallis-

28、steen试验结果 通过实验研究认为 Rel200时， 较低，即使气相流速很高，也几 乎没有携带，E=0; 当200Rel3000时，液膜中开始紊流流动，则,111, 当Rel3000 或无因次液相速度 时， 液膜处于完全紊流中 Wallis通过试验给出了E与 的关系 其中 a.当 时， b.当 时，,112,4. 补充方程,5. 计算程序 例题 已知：气-水混合物向上流过内径为30mm垂直管，,求：,113,计算过程 1.计算迭代初值 (1).求jg、jf (2).求Rel、Reg、判别流态 (3).求 (4).求 (5).求 (6).求 2.求窜流比E 3.求液膜的分液相摩阻梯度 (1).

29、液膜流量 (2).液膜的分液相雷诺数,114,(3).液膜流量的分液相摩阻系数 (4).液膜摩阻梯度 4.液膜流量的分液相折算系数 5.利用补充方程求 和 6.求 和 7.返回到第4步，进行迭代计算 8.判别：最后的压降梯度值与前一次的压降梯度差值小于容许误差，迭代结束。,115,本章小结,1.基本概念 2.直管中压降计算(摩擦、重位、加速) 3.摩擦压降计算（均相、分相模型） 4.重位压降和加速压降计算 5.什么叫窜流比，影响E大小的因素有哪些，如何影响？ 6.雾环状流 的计算程序,116,第七章 两相临界流动,主要内容: 1.两相临界流特征及判断临界流动的判据； 2.掌握运用平衡均相模型，

30、福斯克模型求 解临界质量流速的方法； 3. 短孔道内亚稳态流动。,117,二. 单相临界流,上游条件：Po,to 下游条件：Pb,(1)Po=Pb ,没有流动；,(2)PbPo ,开始流动； 随着Pb的降低，压力 梯度不断增大，而流量 在不断的增加,(3)Po不变，Pb降低到某 一值时，流量不再随着 Pb的减小而增大，此时 流量达到最大值。,118,临界流动的判据,1.临界流动的定义： 系统的某一部分中的流动(流量)仅受在上游条件(Po,to) 的限制，而不受一定范围内变化的下游条件(Pb)影响的流动。 2.当地音速 临界流动时孔口处的流速称为当地音速。,或,或,119,研究两相临界流动的目的

31、 1.为安全注射系统提供设计依据，根据冷却 剂系统、主管道尺寸等，算出临界流动下的流量， 为事故保护系统设计参考。 2.估计产生事故的危害性，这与反应堆的安装 场所有关。根据具体情况，从安全的角度提出有效 的对策。,120,假设: (1)一元水平稳定流动，且dA=0; (2)等熵过程，忽略摩阻，与外界无热交换； (3)对外不做功。 推导： (1)动量方程 (2)连续性方程 (3)临界流动判据,3.单相临界流Gc方程的推导,121,7.3 两相临界流的平衡均相模型,一.两相临界流的特征 1.发生相变； 2.流型变化; 3.两相间存在相对运动，S1; 4.两相间存在热力学不平衡现象。 二.平衡均相

32、模型 1.假设条件 (1)两相处于热平衡，流道足够长； (2)两相间无相对运动，S=1，均质单相流； (3)等熵过程； (4)一元水平稳定流动，对外不做功。,122,2.公式推导 根据假设条件，把两相流当作均质单相流,3.讨论 图7.3中横坐标：Xe 孔口处含气量 纵坐标： (1) (2),即,，计算结果偏低。,说明平衡均相模型不能用于短管内临界流量的计算。,123,计算（按等焓条件处理）,等焓过程：,热平衡含气率：,整理得到：,式中：,与 的计算,124,3.讨论 图7.3中横坐标：Xe 孔口处含气量 纵坐标： (1),即,，计算结果偏低。,说明平衡均相模型不能用于短管内临界流量的计算。,(

33、2),125,上节内容回顾：,两相临界流动定义 两相临界流动的判据 第一判据 或 两相临界流量的计算模型 平衡均相模型,126,7.4 长孔道内两相临界流,一.Fauske（福斯克）模型 1.假设条件 (1)两相处于热平衡； (2)两相存在相对运动，S1;且按分相流处理； (3)等熵过程； (4)一元水平流动，对外不做功。,127,2.Gc计算式,128,(1).S和 的计算 孔口处：压力梯度的绝对值达到有限最大值，即,=（有限）最大值,临界流的第二个判据,Fauske(1962),129,(1).S和 的计算,临界流动时，必有,根据动量方程,130,在临界流动条件下，孔口处 最大，S应达到最

34、大值。 故孔口处,131,(2). 计算（按等焓条件处理）,等焓过程：,热平衡含气率：,整理得到：,式中：,132,1)滞止参数 流动系统中流速为零状态下相应参数称为滞止参数。 2)公式推导 假设长孔道中，流体即不做功，又没有热交换。则能 量方程为,(3).X计算,气液两相可分别写出,两相混合物滞止焓,133,(4).临界压力Pc的确定,临界压力比,当 长孔道。,时，属于,134,3. Fausko模型计算Gc的步骤,建立方 程组联 立求解 Gc,135,计算步骤： (1)求Pc (2)有Pc查表： 用差值比计算： (3)由Po,to查表： (4)用上述方程组进行迭代求解： 迭代初值取,136

35、,二. Moody模型 1.基本假设 (1) 两相处于热平衡； (2) 两相存在相对运动，S1;且按分相流处理； (3) 等熵过程； (4) 一元水平流动，对外不做功。 与Fauske模型相同。 2.根据分相流能量方程，导出,137,3.S的计算,临界流动时滑速比,4. Gc的计算,138,例题：7-1,分析步骤：,（1）确定求解临界流速的公式 （2）根据出口压力，确定出口状态参数 （3）查表确定出口参数 （4）确定计算公式中的未知量,139,单相流：,以平衡均相模型为例：,140,7.5 短孔道内两相临界流,短孔道内的亚稳态流动 1.亚稳态流动 在临界流动情况下，液体的相变过程(扩容膨胀)

36、滞后于压力降低的变化而使液体过热的现象。 2.特点 短孔道内亚稳态流动时出口含气率大大低于假定 热力学平衡时的计算含气率。 3.产生亚稳态流动的原因 (1)气化核心数目不足； (2)表面张力阻碍气泡的生成； (3)相间传热困难。,141,二.短孔道内Gc的计算,1. 孔板，锐边进口,2. 喷嘴，锐边进口,式中，po上游压力，Pa; pb 背压，Pa; 水的密度，由Po,to查表。,式中，pc孔口处临界压力，Pa;,Po,to=const,142,3. ：短管,EL.Wakil建议按照图7-10查图确定 横坐标：Po , MPa 纵坐标： , kg/m2.s 或： 对于饱和流动：,143,(1)

37、直角锐边进口 收缩流离壁面较远，壁面条件对Gc没有影响。 (2)圆角进口 受壁面粗糙度的影响，壁面与液芯射流之间的蒸汽泡 或空气泡与液体之间产生一定的搅混，因而汽化有两种 方式： a.液芯射流表面汽化 b.汽环中蒸汽泡、空气泡起到汽化核心的作用,汽化增多， 阻塞提前。,因此，对于圆角进口，受壁面效应的影响，汽化核心增多， 汽化加剧，导致临界压力增加，临界质量流速减小；而且粗糙度 越大，临界质量流速越低。,三.入口条件及壁面条件对Gc的影响,144,四.伯内尔(Burnell)方法,式中，po上游压力，Pa; 由po,to查表，kg/m3; pS上游温度to对应的饱和压力，Pa; C 考虑表面张力影响的经验系数，查图7-11 或,Ps对应的表面张力，N/m; Ps=1.2MPa时的,适用条件：直角锐边节流孔、孔洞或管嘴 。,145,本 章 小 结,1.何谓临界流动,长管内两相临界流的判据？ 2.两相临界流动的特征？ 3.平衡均相模型、Fanske模型基本假设？ 4.何谓亚稳态流动？产生亚稳态流动的原因？ 5.短孔道内入口条件及壁面条件对Gc有何影响。 6.长孔道内Gc计算（运用平衡均相模型）。 8.短管内Gc计算。,