下载文档原格式
基于CFD技术的超音速喷嘴两相流破碎机制研究朱玲玲;吴建军;刘明翔;隋大山;崔振山【摘要】采用CFD(computational fluid dynamics,计算流体力学)软件系统研究超音速气雾化喷嘴两相流的雾化过程.利用VOF(volume of fluid,流体体积)函数两相流模型模拟验证金属液不同质量流率下的2种初级破碎模式,并研究雾化压力和液体表面张力对金属液初级破碎过程的影响.模拟结果表明:金属液质量流率较小(0.053 kg/s)时,初级破碎模式为液膜破碎,金属液质量流率较大(0.265 kg/s)时,初级破碎模式为"微型喷泉"破碎;随雾化压力从0.5 MPa增加到1.5 MPa,初级破碎程度加剧,但雾化压力过高反而会削弱雾化效果;将金属液表面张力由1.2 N/m降至0.4 N/m,初级破碎时能够获得尺寸更细小的液滴,通过随后的二次破碎形成更加均匀细小的液滴,从而获得高质量的沉积锭.%The atomization process of gas-liquid two-phase flows in ultrasonic gas atomizer was analyzed using computational fluid dynamics software.The two kinds of primary-breakup models of liquid metal at different mass flow rates were verified by volume of fluid model,and the effects of atomizing pressure and liquid surface tension on the primary-breakup process were also studied.The simulation results indicate that the transformation of primary-breakup model will occur with decreasing the liquid mass flow rate.When the liquid mass flow rate is 0.053 kg/s,the primary breakup model is melt sheet,but when it is 0.265 kg/s,the primary breakup model is "fountain".When the atomizing pressure increases from 0.5 MPa to 1.5 MPa,the primary-breakup degree aggravates.However,excessive atomizing pressure will weaken the atomization effect.The smaller-size drop can be formed in the primary-breakup process when the surface tension of liquid metal decreases from 1.2 N/m to 0.4 N/m.Through the following secondary-breakup,more uniform and finer powder particles and high-quality as-spray ingot will be obtained.【期刊名称】《粉末冶金材料科学与工程》【年(卷),期】2018(023)003【总页数】9页(P229-237)【关键词】喷射成形;CFD;超音速喷嘴;两相流;破碎机制【作者】朱玲玲;吴建军;刘明翔;隋大山;崔振山【作者单位】上海交通大学材料科学与工程学院塑性成形系,上海 200030;上海交通大学材料科学与工程学院塑性成形系,上海 200030;上海交通大学材料科学与工程学院塑性成形系,上海 200030;上海交通大学材料科学与工程学院塑性成形系,上海 200030;上海交通大学材料科学与工程学院塑性成形系,上海 200030【正文语种】中文【中图分类】TF123.112喷射成形工艺中，金属液在高速气流的冲击作用下发生雾化并破碎成细小液滴，这是气液两相机械作用和传热的复杂物理过程，涉及两相之间剧烈的能量、动量和热量交换[1]。
第３９卷 第４期Ｖｏｌ．３９ Ｎｏ．４袁寿其基于ＭＵＳＩＧ模型的气液两相流离心泵内部流动数值模拟袁寿其，何文婷，司乔瑞 ，袁建平，张皓阳，崔强磊（江苏大学国家水泵及系统工程技术研究中心，江苏镇江２１２０１３）收稿日期：２０２０－０４－１０；修回日期：２０２０－０５－０２；网络出版时间：２０２１－０４－０７网络出版地址：ｈｔｔｐｓ：／／ｋｎｓ．ｃｎｋｉ．ｎｅｔ／ｋｃｍｓ／ｄｅｔａｉｌ／３２．１８１４．ｔｈ．２０２１０４０６．１１３９．０２６．ｈｔｍｌ基金项目：国家自然科学基金资助项目（５１７７９１０７）第一作者简介：袁寿其（１９６３—），男，上海金山人，研究员，博士生导师（ｓｈｏｕｑｉｙ＠ｕｊｓ．ｅｄｕ．ｃｎ），主要从事流体机械及排灌机械研究．通信作者简介：司乔瑞（１９８６—），男，河南开封人，研究员，博士生导师（ｓｉｑｉａｏｒｕｉ＠ｕｊｓ．ｅｄｕ．ｃｎ），主要从事泵内部流动机理及特性研究．摘要：为了解决气液两相流泵内部流动数值模拟中所采用的欧拉－欧拉双流体非均相流模型无法考虑气泡离散相粒子直径变化以及气相之间的聚合作用与破碎作用，导致在高含气量时的模拟结果与试验存在一定差距的问题，文中将一种新型的欧拉－欧拉双流体拓展模型，即ＭＵＳＩＧ模型用于气液两相流泵内部流动的数值模拟，通过与气液两相流工况外特性试验数据对比发现，入口含气率在５％左右时，ＭＵＳＩＧ模型计算得到的外特性曲线与试验结果整体趋势吻合较好；普通的欧拉－欧拉两相流模型在大含气率下与试验相差较大．基于ＭＵＳＩＧ模型，分析入口含气率对内部流动特性的影响，发现入口含气率的增加会引起内流失稳和流线紊乱，气相逐渐聚集在前盖板与流道中间部位，最后引起能量损失，叶轮出口压力下降．这些现象会随着入口含气率增加而逐渐加剧，最终扬程与效率均会随着含气率增加而下降．含气率小于３％时，内流较稳定；当入口含气率为５％时，扬程下降至３２ｍ，效率下降至５５％，推测此时流道内气相聚合，生成气囊．关键词：离心泵；气液两相流；ＭＵＳＩＧ模型；欧拉－欧拉双流体模型；外特性曲线中图分类号：Ｓ２７７．９ 文献标志码：Ａ 文章编号：１６７４－８５３０（２０２１）０４－０３２５－０６Ｄｏｉ：１０．３９６９／ｊ．ｉｓｓｎ．１６７４－８５３０．２０．０１０４ 袁寿其，何文婷，司乔瑞，等．基于ＭＵＳＩＧ模型的气液两相流离心泵内部流动数值模拟［Ｊ］．排灌机械工程学报，２０２１，３９（４）：３２５－３３０，３３７． ＹＵＡＮＳｈｏｕｑｉ，ＨＥＷｅｎｔｉｎｇ，ＳＩＱｉａｏｒｕｉ，ｅｔａｌ．Ｎｕｍｅｒｉｃａｌｓｉｍｕｌａｔｉｏｎｏｎｇａｓ－ｌｉｑｕｉｄｔｗｏ ｐｈａｓｅｆｌｏｗｉｎｃｅｎｔｒｉｆｕｇａｌｐｕｍｐｂａｓｅｄｏｎＭＵＳＩＧｍｏｄｅｌ［Ｊ］．Ｊｏｕｒｎａｌｏｆｄｒａｉｎａｇｅａｎｄｉｒｒｉｇａｔｉｏｎｍａｃｈｉｎｅｒｙｅｎｇｉｎｅｅｒｉｎｇ（ＪＤＩＭＥ），２０２１，３９（４）：３２５－３３０，３３７．（ｉｎＣｈｉｎｅｓｅ）Ｎｕｍｅｒｉｃａｌｓｉｍｕｌａｔｉｏｎｏｎｇａｓ－ｌｉｑｕｉｄｔｗｏ ｐｈａｓｅｆｌｏｗｉｎｃｅｎｔｒｉｆｕｇａｌｐｕｍｐｂａｓｅｄｏｎＭＵＳＩＧｍｏｄｅｌＹＵＡＮＳｈｏｕｑｉ，ＨＥＷｅｎｔｉｎｇ，ＳＩＱｉａｏｒｕｉ，ＹＵＡＮＪｉａｎｐｉｎｇ，ＺＨＡＮＧＨａｏｙａｎｇ，ＣＵＩＱｉａｎｇｌｅｉ（ＮａｔｉｏｎａｌＲｅｓｅａｒｃｈＣｅｎｔｅｒｏｆＰｕｍｐｓ，ＪｉａｎｇｓｕＵｎｉｖｅｒｓｉｔｙ，Ｚｈｅｎｊｉａｎｇ，Ｊｉａｎｇｓｕ２１２０１３，Ｃｈｉｎａ）Ａｂｓｔｒａｃｔ：ＴｈｅＥｕｌｅｒｉａｎ－Ｅｕｌｅｒｉａｎｔｗｏ ｐｈａｓｅｆｌｏｗｍｏｄｅｌｕｓｅｄｉｎｔｈｅｓｉｍｕｌａｔｉｏｎｏｆｔｈｅｔｈｅｇａｓ－ｌｉｑｕｉｄｔｗｏ ｐｈａｓｅｆｌｏｗｐｕｍｐｃａｎｎｏｔｃｏｎｓｉｄｅｒｔｈｅａｉｒｂｕｂｂｌｅｓｄｉａｍｅｔｅｒｃｈａｎｇｉｎｇ．Ｔｈｅｒｅｉｓａｃｅｒｔａｉｎｇａｐｂｅｔｗｅｅｎｔｈｅｓｉｍｕｌａｔｉｏｎｒｅｓｕｌｔａｎｄｔｈｅｔｅｓｔｗｈｅｎｔｈｅｇａｓｃｏｎｔｅｎｔｉｓｒｅｌａｔｉｖｅｌｙｈｉｇｈ．ＴｈｅＭＵＳＩＧｍｏｄｅｌ，ｗｈｉｃｈｉｓａｎｅｗｔｙｐｅｏｆｅｘｐａｎｓｉｏｎｍｏｄｅｌｂａｓｅｄｏｎｏｒｉｇｉｎａｌＥｕｌｅｒｉａｎ－Ｅｕｌｅｒｉａｎｍｏｄｅｌ，ｃａｎｃａｌｃｕｌａｔｅｍｏｒｅａｃｃｕｒａｔｅｌｙ．Ｃｏｍｐａｒｉｓｏｎｗｉｔｈｔｈｅｅｘｐｅｒｉｍｅｎｔａｌｄａｔａｓｈｏｗｓｔｈａｔｗｈｅｎｔｈｅｇａｓｃｏｎｔｅｎｔｉｓａｂｏｕｔ５％，ｔｈｅｈｙｄｒａｕｌｉｃｃｈａｒａｃｔｅｒｉｓｔｉｃｃｕｒｖｅｃａｌｃｕｌａｔｅｄｂｙｔｈｅＭＵＳＩＧｍｏｄｅｌｈａｓａｓｌｉｇｈｔｄｅｖｉａｔｉｏｎｆｒｏｍｔｈｅｅｘｐｅｒｉ ｍｅｎｔａｌｒｅｓｕｌｔｓ，ｂｕｔｔｈｅｏｖｅｒａｌｌｔｒｅｎｄｉｓｉｎｇｏｏｄａｇｒｅｅｍｅｎｔ．Ｂａｓｅｄｏｎｔｈｅｏｒｉｇｉｎａｌｍｏｄｅｌ，ｔｈｅｆｌｏｗｃｈａ ｒａｃｔｅｒｉｓｔｉｃｓｉｎｓｉｄｅｔｈｅｐｕｍｐｆｏｒｄｉｆｆｅｒｅｎｔＩＧＶＦ（ｉｎｌｅｔｇａｓｖｏｌｕｍｅｆｒａｃｔｉｏｎ）ｓｈｏｗｔｈａｔｔｈｅｉｎｃｒｅａｓｅｏｆｔｈｅｉｎｌｅｔｇａｓｃｏｎｔｅｎｔｗｉｌｌｃａｕｓｅｅｎｔｉｒｅｆｌｏｗｉｎｓｔａｂｉｌｉｔｙａｎｄｓｔｒｅａｍｌｉｎｅｔｕｒｂｕｌｅｎｃｅ；ｔｈｅｇａｓｐｈａｓｅｗｉｌｌｇｒａｄｕａｌ ｌｙａｃｃｕｍｕｌａｔｅｉｎｔｈｅｍｉｄｄｌｅｏｆｔｈｅｓｈｏｕｌｄａｎｄｔｈｅｆｌｏｗｃｈａｎｎｅｌ；ｔｈｅｐｒｅｓｓｕｒｅｏｆｔｈｅｉｍｐｅｌｌｅｒｄｒｏｐｓ．ＴｈｅｓｅｐｈｅｎｏｍｅｎａｗｉｌｌｇｒａｄｕａｌｌｙｉｎｔｅｎｓｉｆｙｗｉｔｈｔｈｅＩＧＶＦｉｎｃｒｅａｓｉｎｇ，ａｎｄｈｅａｄａｎｄｅｆｆｉｃｉｅｎｃｙｗｉｌｌｄｅ ｃｒｅａｓｅｏｂｖｉｏｕｓｌｙ．Ｗｈｅｎｔｈｅｇａｓｃｏｎｔｅｎｔｉｓｌｅｓｓｔｈａｎ３％，ｔｈｅｉｎｔｅｒｎａｌｆｌｏｗｉｓｒｅｌａｔｉｖｅｌｙｓｔａｂｌｅ；ｗｈｅｎｔｈｅｉｎｌｅｔｇａｓｃｏｎｔｅｎｔｉｓ５％，ｔｈｅｐｅｒｆｏｒｍａｎｃｅｗｉｌｌｓｕｄｄｅｎｌｙｄｒｏｐ，ａｎｄｔｈｅｇａｓａｃｃｕｍｕｌａｔｉｏｎｏｃｃｕｒｓｉｎｆｌｏｗｃｈａｎｎｅｌ．Ｋｅｙｗｏｒｄｓ：ｃｅｎｔｒｉｆｕｇａｌｐｕｍｐ；ｇａｓ－ｌｉｑｕｉｄｔｗｏ ｐｈａｓｅｆｌｏｗ；ＭＵＳＩＧｍｏｄｅｌ；Ｅｕｌｅｒｉａｎ－Ｅｕｌｅｒｉａｎｔｗｏ ｐｈａｓｅｍｏｄｅｌ；ｐｕｍｐｐｅｒｆｏｒｍａｎｃｅｃｕｒｖｅ 离心泵广泛应用于石油、化工及天然气输送、核电站等领域，与国民生计关系紧密［１－２］．由于泵的工作环境千差万别，常常出现各种特殊情况，例如多相流、空化、汽蚀等现象的发生．气液两相流是较为常见的一种多相流，指纯水流体中混入一定含量的空气．诸多研究发现，离心泵的性能与内部流动都会受到入口气体含量的影响．袁建平等［３］认为含气率达到１０％时，会出现气液分离现象，造成离心泵输水性能急剧下降；含气率逐渐增加的过程中，叶轮出口压力逐渐降低，说明气液两相流会造成离心泵一定程度的扬程损失．唐苑峰等［４］认为离心泵外特性会随着含气率增加而下降．ＳＨＡＯ等［５］通过可视化试验，证实了随含气率的增加，内部流动会呈现４种气液两相流流态，即泡状流、聚合泡状流、气囊流和气液分离流；结合离心泵外特性变化规律分析发现，两相流流型与外特性存在紧密联系，例如扬程与效率出现骤降时，同时也会伴随气囊流的出现．当前，随着纯水工况下离心泵内部流动数值模拟研究的深入，已经可以准确预测泵的性能．而对气液两相流泵的数值模拟工作开展相对较少，ＳＩ等［６］基于欧拉－欧拉双流体非均相模型，模拟分析了不同含气率下离心泵性能以及内流变化规律，但模拟中并未考虑气泡形变和气泡碰撞合并现象．当前，基于欧拉－欧拉双流体模型，将气体粒子的直径变化规律、聚合以及分裂等加入计算，延伸出ＰＢＭ－ＣＦＤ耦合模型与ＭＵＳＩＧ模型［７］．戈振国［８］基于ＰＢＭ－ＣＦＤ耦合模型，模拟分析了气泡直径和流型转变的规律对离心泵性能的影响．文中进一步探究另一种基于欧拉－欧拉双流体的拓展模型———ＭＵＳＩＧ模型用于气液两相流泵内部流动数值模拟的可靠性，进而分析含气率对内流场以及外特性的影响．１ 数值计算１．１ 模型泵参数模型泵采用单极单吸离心泵，基本参数中，比转数ｎｓ＝８８．６，叶片数Ｚ＝６，纯水工况下额定流量Ｑｄ＝５０ｍ３·ｈ，额定扬程Ｈ＝３４ｍ，额定转速ｎ＝２９００ｒ／ｍｉｎ．叶轮进口直径Ｄ１＝７４ｍｍ，叶轮出口直径Ｄ２＝１７４ｍｍ，叶片出口宽度ｂ２＝１２ｍｍ，泵进口直径Ｄｓ＝６５ｍｍ，泵出口直径Ｄｄ＝６５ｍｍ．１．２ 网格划分及边界条件设置１．２．１ 网格划分与无关性分析图１，２分别为模型泵部件网格划分与总装配图．采用结构化网格，水体结构分为进口管、口环、叶轮、泵腔、蜗壳与出口管６个部分，网格总数为４００万．图１ 模型泵水利部件网格示意图Ｆｉｇ．１ Ｇｒｉｄｓｏｆｐｕｍｐ图２ 模型泵三维总装配图Ｆｉｇ．２ Ｔｈｒｅｅ ｄｉｍｅｎｓｉｏｎａｌｇｅｎｅｒａｌａｓｓｅｍｂｌｙｄｒａｗｉｎｇｏｆｍｏｄｅｌｐｕｍｐ为了确定适用于计算的网格数，进行无关性分３２６析，发现网格数达到４００万时，扬程趋于平稳．最终确定网格总数为４００万，并且最小角度未小于１８°，可以保证计算精度．１．２．２ 边界条件设置基于定常计算离心泵的扬程与效率，基于非定常分析内流规律．将２５℃的纯水以及空气混合作为流体介质，入口边界设定气相体积分数，入口压力为１０１．３２５ｋＰａ，出口设置为质量流量．ＭＵＳＩＧ模型设置：气相离散流体，入口气泡组分组设置为１２组，最小粒径为０．１ｍｍ，最大粒径为１．０ｍｍ．粒子分组过少，导致计算不符合实际情况，分组过多，计算耗时较长且不易收敛［９］．欧拉－欧拉非均相流模型设置：气相为离散流体，平均粒径设置为０．１ｍｍ．非定常设置：非定常计算时设为ＴｒａｎｓｉｅｎｔＲｏｔｏｒＳｔａｔｏｒ，并且以定常计算结果作为非定常计算的初始条件，取叶轮每转３°计算１次，则时间步长Δｔ＝１．７２４×１０－４ｓ，总时间为０．２０６９ｓ，对最后稳定的５周数据进行处理．１．３ ＭＵＳＩＧ模型气液两相流泵内部两相流模型是选择了基于欧拉－欧拉双流体模型的一种新型拓展模型，即ＭＵＳＩＧ模型．该模型与欧拉－欧拉双流体非均相流模型不同的是，考虑离散相粒子直径的变化，以及气相之间的聚合作用与破碎作用．其基本原理是，将所有粒子划分为粒径不同的Ｎ组，每一组粒子拥有各自独立的连续性方程，但却拥有同一套动量方程，最终去求解Ｎ＋１组方程．连续性方程为ｔ（αｋρｋ）＋ ·（αｋρｋωｋ）＝０，（１）动量方程为ｔ（αｋρｋωｋ）＋ ·（αｋρｋωｋ ωｋ）＝－αｋ ｐｋ＋ ·｛αｋμｋ［ ωｋ＋（ ωｋ）Ｔ］｝＋Ｍｋ＋αｋρｋｆｋ，（２）式中：ｋ为任意相；ρｋ为ｋ相密度；ｐｋ为ｋ相压力；αｋ为ｋ相体积分数；μｋ为ｋ相动力黏度；ωｋ为ｋ相流体相对速度；Ｍｋ为ｋ相所受相间作用力；ｆｋ为与叶轮旋转有关的质量力．ＭＵＳＩＧ模型中采用的气泡粒子破碎模型来源于ＬＵＯ等［１０］针对湍流场中的粒子破碎行为的研究，气体粒子的合并模型来源于ＰＲＩＮＣＥ等［１１］的研究．２ 外特性试验验证２．１ 气液两相流泵实验台及试验过程介绍试验采用开式实验台，由开式储水箱和进出口管路等组成纯水试验回路，纯水试验完成后，加入压缩机、气液混合器和输气管等组成气体输送回路［１２］，实验台示意图如图３所示．图３ 试验示意图Ｆｉｇ．３ Ｅｘｐｅｒｉｍｅｎｔａｌｓｃｈｅｍａｔｉｃ２．２ 气液两相流数值模拟与试验验证图４为额定工况下，普通两相流模型与ＭＵＳＩＧ模型的扬程与效率随入口含气率（ＩＧＶＦ）变化曲线对比．额定工况下，纯水扬程为３５ｍ，随着入口含气率增加扬程有明显的下降趋势．其中，入口含气率未超过３％时，扬程无明显波动，２种模型计算结果与试验结果都较为接近，说明这时内流相对较为稳定，气相并未造成明显的水力损失．当入口含气率到达５％时，试验扬程出现了骤降，试验过程中，通过调节气体阀门增加入口含气率，同时调节液体阀门使得液体质量流量不变，这个过程入口质量流量出现波动，但是最终维持稳定状态．试验结果显示扬程从３５ｍ降至２６ｍ左右，ＭＵＳＩＧ模型结果也出现了骤降，从３５ｍ降至３１ｍ左右，较试验略有偏高，普通模型结果只出现略微下降趋势．入口含气率继续增加，ＭＵＳＩＧ模拟结果与试验结果都出现了更剧烈的骤降．当含气率到达１０％时，模拟扬程与试验扬程都下降至２０ｍ以下．整体比较发现，普通气液两相流模型在含气率超过３％时，与实际扬程误差增大，不能模拟出试验中扬程骤降，而ＭＵＳＩＧ模型能够较好地与试验扬程相匹配．３２７图４ 额定流量下的扬程与效率随入口含气率（ＩＧＶＦ）变化曲线Ｆｉｇ．４ Ｈｅａｄａｎｄｅｆｆｉｃｉｅｎｃｙｃｈａｎｇｅｗｉｔｈｉｎｌｅｔｇａｓｃｏｎｔｅｎｔ（ＩＧＶＦ）ｕｎｄｅｒｔｈｅｇａｓ－ｌｉｑｕｉｄｔｗｏ ｐｈａｓｅｆｌｏｗｃｏｎｄｉｔｉｏｎｆｏｒＱ＝Ｑｄ随着含气率增加，离心泵的效率也逐步下降，含气率为５％时，试验效率从７９％骤降至４８％，效率骤降说明含气率到达５％时，输水能力急剧下降．最终效率在含气率为１０％时下降至３５％左右，离心泵输水能力损失较严重．当含气率在５％左右时，ＭＵＳＩＧ模型计算得到的外特性曲线与试验结果有少量偏差，扬程相差３ｍ左右，效率相差８％左右，而在含气率较小时则均有非常好的一致性．因为少量的气体并不能形成严重的水利损失，流场稳定程度接近于纯水工况，但是ＭＵＳＩＧ模型整体趋势的符合度非常好．入口含气率超过５％时，普通两相流模型没有将气相粒子的形变考虑在计算过程中，计算结果与实际情况差距较大，入口含气率增加至７％时，误差达到最大值．３ 内流分析３．１ 流线分布图５为纯水工况下，叶轮内液体流线分布（额定工况）．流体流过高速旋转的叶轮时，流线排列越顺畅，造成的水力损失就越小．额定流量，纯水工况下，液体流线均匀有序地充满整个流道，并没有出现剧烈的速度波动，也没有出现由流场不稳定引起的流动行为，例如回流等．所以，纯水工况离心泵几乎可以达到稳定状态，流场中没有多余的能量损耗以及水力损失，扬程和效率都较高．图５ 额定流量下入口含气率为０时流线分布Ｆｉｇ．５ ＳｔｒｅａｍｌｉｎｅｄｉｓｔｒｉｂｕｔｉｏｎｆｏｒＱ＝Ｑｄ，ＩＧＶＦ＝０图６为ＩＧＶＦ＝３％时，离心泵内液体流线分布．入口含气率为３％时，流线沿着流道排列顺畅，未发现回流等流动行为．３％的含气率并没有对液体流线造成明显干扰，和纯水工况相比，流线排列与速度大小无明显波动．速度流线显示，叶轮出口到蜗壳部位速度较高，蜗壳部位整体速度略高于叶轮内的流体速度，蜗壳出口部位速度有所减小．图６ 额定流量下入口含气率为３％时流线分布Ｆｉｇ．６ ＳｔｒｅａｍｌｉｎｅｄｉｓｔｒｉｂｕｔｉｏｎｆｏｒＱ＝Ｑｄ，ＩＧＶＦ＝３％图７为ＩＧＶＦ＝５％时，离心泵内液体流线分布．入口含气率为５％时，在叶轮出口处，紧贴叶片吸力面一侧的流线速度有明显下降，大约从１７ｍ下降至２ｍ左右，并且这一部分液体无法顺利流至流道出口．液体流速下降剧烈甚至有向低压区倒流的趋势，说明水体的内能被消耗，导致流体没有充足的能量流至出口．含气率为５％时，相比较含气率为３％，水力损失加剧，所以入口的含气率与形成的水力损失成正比，这一点与许多研究成果相符．图７ 额定流量下入口含气率为５％时流线分布Ｆｉｇ．７ ＳｔｒｅａｍｌｉｎｅｄｉｓｔｒｉｂｕｔｉｏｎｆｏｒＱ＝Ｑｄ，ＩＧＶＦ＝５％３２８图８为ＩＧＶＦ＝７％时，在额定流量下离心泵内液体流线分布．入口含气率为７％时会造成较为严重的水力损失，因为流道中２／３的流体流线速度下降至２ｍ左右，随着液体速度骤降，输水能力随之减弱．个别流道中部位置出现明显回流，回流区域堵塞了一部分流道，使得流道的利用率降低．同时回流造成流体能量损耗，速度波动形成新的回流，加剧了水力损失．图８ 额定流量下入口含气率为７％时流线分布Ｆｉｇ．８ ＳｔｒｅａｍｌｉｎｅｄｉｓｔｒｉｂｕｔｉｏｎｆｏｒＱ＝Ｑｄ，ＩＧＶＦ＝７％整体看来，气液两相流必然会影响离心泵的内流与外特性，并且随着含气率增加，这种不良影响会加剧，扬程与效率会出现骤降，流线速度骤降并且出现回流等流动行为．３．２ 局部含气率分布图９为随入口含气率增加，叶轮内部局部含气率α分布情况．局部含气率可以直接反映气相在流道中的聚集情况．可视化试验研究发现气相聚集现象会引发气囊流流型（气相在流道中逐渐合并，形成气袋），由于气囊占用一部分流道，消耗流体内能，从而输水性能骤降．图９ 不同入口含气率下叶轮内局部含气率分布Ｆｉｇ．９ ＬｏｃａｌｇａｓｆｒａｃｔｉｏｎｄｉｓｔｒｉｂｕｔｉｏｎｓｉｎｉｍｐｅｌｌｅｒｗｉｔｈｄｉｆｆｅｒｅｎｔＩＧＶＦ入口含气率从３％依次递增至７％，前盖板气相聚集越来越明显．最初含气率为３％时，流道整体气体含率在１０％以下，因为气体流动状态整齐有序，粒子分布松散，相互聚合碰撞的概率很小．含气率达到５％之后，流道中部局部含气率超过８０％，并且紧贴叶片吸力面的含气率有所上升；最终可以在入口含气率为７％时明显观察到气囊，并可发现个别流道很大一部分面积局部含气率超过８０％．观察后盖板，局部的气体含量较低，从侧面说明了气囊更易发生在贴近前盖板的位置．入口含气率达到７％时，后盖板贴近叶轮的极小部位含气率有所上升．３种入口含气率条件下，后盖板的局部含气率均维持在１０％以下，说明含气率的增加对后盖板附近部位的流场并没有明显影响．３．３ 压力分布图１０所示为随入口含气率增加叶轮内部压力分布情况．由于扬程代表了离心泵的加压能力，如果离心泵内部的水力损失造成压降，会直接引起扬程的骤降．由于叶轮为主要的加压部位，所以流道的压力分布变化可以说明气液两相流对叶轮加压能力的影响．图１０ 不同入口含气率下叶轮压力分布Ｆｉｇ．１０ ＰｒｅｓｓｕｒｅｄｉｓｔｒｉｂｕｔｉｏｎｓｉｎｉｍｐｅｌｌｅｒｗｉｔｈｄｉｆｆｅｒｅｎｔＩＧＶＦ 叶轮是主要离心部件，通过将水体加压赋予水体能量，流体在出口与进口的压差越大，离心泵的扬程就越高．纯水工况下，流体正常加压，离心泵运转正常．流体在流道中部压力迅速升高，叶片工作面中部与叶片吸力面尾部压力急剧升高．当含气率逐步增加时，叶轮的加压能力显著减弱，由于叶片工作面与背面的高压区依次减少．含气率为３％时，相比较纯水工况，压力有微小损失，叶片压力面中部偏下部位，压力达到３００ｋＰａ，到叶片末尾压力逐步增加．含气率超过５％时，叶片压力面的压降十分明显，直到叶片尾部压力才升至３００ｋＰａ，叶片背面尾部的最大压力下降至２７０ｋＰａ左右．３２９４ 气相形变规律图１１为额定流量下，气泡直径ｄｂ随入口含气率的变化规律．ＭＵＳＩＧ模型与普通的欧拉－欧拉双流体模型相比较，引入了气相气泡粒子聚合与破碎的模型．图１１ 额定流量下不同含气率下气泡直径分布Ｆｉｇ．１１ ＢｕｂｂｌｅｓｉｚｅｄｉｓｔｒｉｂｕｔｉｏｎｓｉｎｉｍｐｅｌｌｅｒｗｉｔｈｄｉｆｆｅｒｅｎｔＩＧＶＦｆｏｒＱ＝Ｑｄ入口含气率为３％时，扬程与效率没有剧烈波动，因此推测此时流态为均匀泡状流，整个流道中的粒子直径几乎均大于０．３ｍｍ，叶轮出口边缘部位气泡粒子直径小于０．２ｍｍ，入口部位至流道中部的气泡直径为０．８ｍｍ左右．含气率增加至５％，这时出现性能骤降，推断入口部位有气囊初生．流道中的气体粒子直径跨度较明显，叶片吸力面靠近出口处为直径小于０．２ｍｍ的粒子，入口部位紧贴压力面的小部分流道粒子直径超过了０．８ｍｍ，说明较大的气泡在入口处合并为气囊，只有较小的气泡流至出口部位；入口含气率增加至７％，整体流道中的气泡明显变大，直径增加至０．５ｍｍ左右，流道入口沿着压力面的部位气相聚集成为气囊，所以显示气泡直径超过１ｍｍ．通过多次计算发现，持续增加最大气泡粒子直径对最终气泡直径分布没有明显影响，因为粒子聚合与分裂模型并不仅仅作用于最大直径的粒子．５ 结 论１）外特性结果显示整体试验与ＭＵＳＩＧ模拟结果匹配度较高，尤其当入口含气率小于３％时．入口含气率大于５％时，普通两相流模型计算准确性低于ＭＵＳＩＧ模型．２）通过ＭＵＳＩＧ模型非定常计算，分别研究了含气率为３％，５％，７％这３种情况时的内流规律．随着含气率增大，液体流线紊乱程度与回流程度加剧；紧贴前盖板位置的气相聚集逐渐严峻（但是后盖板部位没有气体聚集），逐渐从泡状流转变为气囊流；叶轮的加压能力急剧减弱．３）ＭＵＳＩＧ模型非定常计算显示，含气率的增大同时引起流体回流与局部含气率的升高，气囊的出现会占据一部分流道空间，使得叶轮出现一定程度的空转；流体能量被损耗和叶轮加压能力减弱是导致扬程和效率下降的主要原因．ＭＵＳＩＧ模型考虑了气相粒子之间的聚合以及破碎作用，较传统的欧拉－欧拉两相流模型更接近实际观测结果，说明性能骤降与气相的聚合与破碎有密切关系．参考文献（Ｒｅｆｅｒｅｎｃｅｓ）［１］ 王东伟，刘在伦，曾继来．离心泵非定常空化流场及空泡特征分析［Ｊ］．流体机械，２０２０，４８（１２）：２８－３５．ＷＡＮＧＤｏｎｇｗｅｉ，ＬＩＵＺａｉｌｕｎ，ＺＥＮＧＪｉｌａｉ．Ａｎａｌｙｓｉｓｏｆｕｎｓｔｅａｄｙｃａｖｉｔａｔｉｏｎｆｌｏｗｆｉｅｌｄａｎｄｃａｖｉｔａｔｉｏｎｂｕｂｂｌｅｃｈａｒａｃｔｅｒｉｓｔｉｃｓｆｏｒａｃｅｎｔｒｉｆｕｇａｌｐｕｍｐ［Ｊ］．Ｆｌｕｉｄｍａｃｈｉｎｅｒｙ，２０２０，４８（１２）：２８－３５．（ｉｎＣｈｉｎｅｓｅ）［２］ 刘宇宁，王秀勇，刘志远，等．多级离心泵水力性能数值模拟精度影响因素研究［Ｊ］．流体机械，２０２０，４８（１０）：４１－４７．ＬＩＵＹｕｎｉｎｇ，ＷＡＮＧＸｉｕｙｏｎｇ，ＬＩＵＺｈｉｙｕａｎ，ｅｔａｌ．Ｓｔｕｄｙｏｎｔｈｅｉｎｆｌｕｅｎｃｉｎｇｆａｃｔｏｒｓｏｆｎｕｍｅｒｉｃａｌｓｉｍｕｌａｔｉｏｎａｃｃｕｒａｃｙｆｏｒｔｈｅｐｅｒｆｏｒｍａｎｃｅｐｒｅｄｉｃｔｉｏｎｏｆｍｕｌｔｉｓｔａｇｅｃｅｎｔｒｉｆｕｇａｌｐｕｍｐ［Ｊ］．Ｆｌｕｉｄｍａｃｈｉｎｅｒｙ，２０２０，４８（１０）：４１－４７．（ｉｎＣｈｉｎｅｓｅ）［３］ 袁建平，张克玉，司乔瑞，等．基于非均相流模型的离心泵气液两相流动数值研究［Ｊ］．农业机械学报，２０１７，４８（１）：８９－９５．ＹＵＡＮＪｉａｎｐｉｎｇ，ＺＨＡＮＧＫｅｙｕ，ＳＩＱｉａｏｒｕｉ，ｅｔａｌ．Ｎｕｍｅｒｉｃａｌｉｎｖｅｓｔｉｇａｔｉｏｎｏｆｇａｓ－ｌｉｑｕｉｄｔｗｏ ｐｈａｓｅｆｌｏｗｉｎｃｅｎｔｒｉｆｕｇａｌｐｕｍｐｓｂａｓｅｄｏｎｉｎｈｏｍｏｇｅｎｅｏｕｓｍｏｄｅｌ［Ｊ］．ＴｒａｎｓａｃｔｉｏｎｓｏｆｔｈｅＣＳＡＭ，２０１７，４８（１）：８９－９５．（ｉｎＣｈｉｎｅｓｅ）［４］ 唐苑峰，袁建平，司乔瑞，等．基于Ｅｕｌｅｒｉａｎ－Ｅｕｌｅｒｉａｎ模型的轴流泵气液两相流动数值分析［Ｊ］．排灌机械工程学报，２０１８，３６（６）：４７２－４７８．ＴＡＮＧＹｕａｎｆｅｎｇ，ＹＵＡＮＪｉａｎｐｉｎｇ，ＳＩＱｉａｏｒｕｉ，ｅｔａｌ．Ｎｕｍｅｒｉｃａｌａｎａｌｙｓｉｓｏｆｔｈｅｔｗｏ ｐｈａｓｅｆｌｏｗ（ｌｉｑｕｉｄ／ｇａｓ）ｉｎａｘｉａｌｐｕｍｐｂａｓｅｄｏｎＥｕｌｅｒｉａｎ－Ｅｕｌｅｒｉａｎｆｌｏｗｍｏｄｅｌ［Ｊ］．Ｊｏｕｒｎａｌｏｆｄｒａｉｎａｇｅａｎｄｉｒｒｉｇａｔｉｏｎｍａｃｈｉｎｅｒｙｅｎｇｉｎｅｅｒｉｎｇ，２０１８，３６（６）：４７２－４７８．（ｉｎＣｈｉｎｅｓｅ）［５］ ＳＨＡＯＣＬ，ＬＩＣＱ，ＺＨＯＵＪＦ，ｅｔａｌ．Ｅｘｐｅｒｉｍｅｎｔａｌｉｎ ｖｅｓｔｉｇａｔｉｏｎｏｆｆｌｏｗｐａｔｔｅｒｎｓａｎｄｅｘｔｅｒｎａｌｐｅｒｆｏｒｍａｎｃｅｏｆａｃｅｎｔｒｉｆｕｇａｌｐｕｍｐｔｈａｔｔｒａｎｓｐｏｒｔｓｇａｓ－ｌｉｑｕｉｄｔｗｏ ｐｈａｓｅｍｉｘｔｕｒｅｓ［Ｊ］．Ｉｎｔｅｒｎａｔｉｏｎａｌｊｏｕｒｎａｌｏｆｈｅａｔ＆ｆｌｕｉｄｆｌｏｗ，２０１８，７１：４６０－４６９．（下转第３３７页）３３０ｄｏｕｂｌｅｓｕｃｔｉｏｎｃｅｎｔｒｉｆｕｇａｌｐｕｍｐ［Ｊ］．Ｃｈｉｎａｒｕｒａｌｗａｔｅｒｃｏｎｓｅｒｖａｎｃｙａｎｄｈｙｄｒｏｐｏｗｅｒ，２０１７（１０）：２１８－２２２．（ｉｎＣｈｉｎｅｓｅ）［１０］ 周晓红，刘国营，胡永金，等．叶片交错布置对水泵蜗壳内部流动影响［Ｊ］．水电与抽水蓄能，２０１８，４（２）：８２－８５．ＺＨＯＵＸｉａｏｈｏｎｇ，ＬＩＵＧｕｏｙｉｎｇ，ＨＵＹｏｎｇｊｉｎ，ｅｔａｌ．Ｅｆｆｅｃｔｏｆｓｔａｇｇｅｒｅｄｂｌａｄｅａｒｒａｎｇｅｍｅｎｔｏｎｉｎｔｅｒｎａｌｆｌｏｗｏｆｐｕｍｐｖｏｌｕｔｅ［Ｊ］．Ｈｙｄｒｏｐｏｗｅｒａｎｄｐｕｍｐｅｄｓｔｏｒａｇｅ，２０１８，４（２）：８２－８５．（ｉｎＣｈｉｎｅｓｅ）［１１］ 张金凤，黄茜，袁寿其，等．基于ＰＩＶ的低比转数离心泵网格无关性［Ｊ］．排灌机械工程学报，２０１６，３４（７）：５６７－５７２．ＺＨＡＮＧＪｉｎｆｅｎｇ，ＨＵＡＮＧＸｉ，ＹＵＡＮＳｈｏｕｑｉ，ｅｔａｌ．Ｇｒｉｄ ｉｎｄｅｐｅｎｄｅｎｔｏｆｌｏｗｓｐｅｃｉｆｉｃｓｐｅｅｄｃｅｎｔｒｉｆｕｇａｌｐｕｍｐｂａｓｅｄｏｎＰＩＶ［Ｊ］．Ｊｏｕｒｎａｌｏｆｄｒａｉｎａｇｅａｎｄｉｒｒｉｇａｔｉｏｎｍａｃｈｉｎｅｒｙｅｎｇｉｎｅｅｒｉｎｇ，２０１６，３４（７）：５６７－５７２．（ｉｎＣｈｉｎｅｓｅ）［１２］ 李晓俊，袁寿其，潘中永，等．离心泵边界层网格的实现及应用评价［Ｊ］．农业工程学报，２０１２，２８（２０）：６７－７２．ＬＩＸｉａｏｊｕｎ，ＹＵＡＮＳｈｏｕｑｉ，ＰＡＮＺｈｏｎｇｙｏｎｇ，ｅｔａｌ．Ｉｍｐｌｅｍｅｎｔａｔｉｏｎａｎｄａｐｐｌｉｃａｔｉｏｎｅｖａｌｕａｔｉｏｎｏｆｂｏｕｎｄａｒｙｌａｙｅｒｇｒｉｄｏｆｃｅｎｔｒｉｆｕｇａｌｐｕｍｐ［Ｊ］．ＴｒａｎｓａｃｔｉｏｎｓｏｆｔｈｅＣＳＡＥ，２０１２，２８（２０）：６７－７２．（ｉｎＣｈｉｎｅｓｅ）［１３］ 丁成伟．离心泵与轴流泵［Ｍ］．北京：机械工业出版社，１９８１．［１４］ 赵万勇，张亮，雒军．双吸离心泵径向力数值分析［Ｊ］．排灌机械，２００９，２７（４）：２０５－２０９．ＺＨＡＯＷａｎｙｏｎｇ，ＺＨＡＮＧＬｉａｎｇ，ＬＵＯＪｕｎ．Ｎｕｍｅｒｉｃａｌａｎａｌｙｓｉｓｏｆｒａｄｉａｌｆｏｒｃｅｏｆｄｏｕｂｌｅｓｕｃｔｉｏｎｃｅｎｔｒｉｆｕｇａｌｐｕｍｐ［Ｊ］．Ｄｒａｉｎａｇｅａｎｄｉｒｒｉｇａｔｉｏｎｍａｃｈｉｎｅｒｙ，２００９，２７（４）：２０５－２０９．（ｉｎＣｈｉｎｅｓｅ）［１５］ ＧＯＮＺ ＬＥＳＪ，ＳＡＮＴＯＬＡＲＩＡＣ，ＰＡＲＲＯＮＤＯＪＬ，ｅｔａｌ．Ｕｎｓｔｅａｄｙｒａｄｉａｌｆｏｒｃｅｓｏｎｔｈｅｉｍｐｅｌｌｅｒｏｆａｃｅｎｔｒｉｆｕｇａｌｐｕｍｐｗｉｔｈｒａｄｉａｌｇａｐｖａｒｉａｔｉｏｎ［Ｃ］／／ＡＳＭＥ／ＪＳＭＥ２００３４ｔｈＪｏｉｎｔＦｌｕｉｄｓＳｕｍｍｅｒＥｎｇｉｎｅｅｒｉｎｇＣｏｎｆｅｒｅｎｃｅ，２００３：１１７３－１１８１． ［１６］ 张霞，袁寿其，张金凤，等．不同叶片包角对螺旋形单蜗壳离心泵叶轮径向力特性的影响［Ｊ］．流体机械，２０１７，４５（１）：４３－４７．ＺＨＡＮＧＸｉａ，ＹＵＡＮＳｈｏｕｑｉ，ＺＨＡＮＧＪｉｎｆｅｎｇ，ｅｔａｌ．Ｅｆｆｅｃｔｏｆｄｉｆｆｅｒｅｎｔｂｌａｄｅｗｒａｐａｎｇｌｅｓｏｎｒａｄｉａｌｆｏｒｃｅｃｈａｒａｃｔｅｒｉｓｔｉｃｓｏｆｓｐｉｒａｌｓｉｎｇｌｅｖｏｌｕｔｅｃｅｎｔｒｉｆｕｇａｌｐｕｍｐｉｍｐｅｌｌｅｒ［Ｊ］．Ｆｌｕｉｄｍａｃｈｉｎｅｒｙ，２０１７，４５（１）：４３－４７．（ｉｎＣｈｉｎｅｓｅ）（责任编辑 盛杰）檨檨檨檨檨檨檨檨檨檨檨檨檨檨檨檨檨檨檨檨檨檨檨檨檨檨檨檨檨檨檨檨檨檨檨檨檨檨檨檨檨檨檨檨檨檨（上接第３３０页）［６］ ＳＩＱ，ＣＵＩＱ，ＺＨＡＮＧＫ，ｅｔａｌ．Ｉｎｖｅｓｔｉｇａｔｉｏｎｏｎｃｅｎｔｒｉ ｆｕｇａｌｐｕｍｐｐｅｒｆｏｒｍａｎｃｅｄｅｇｒａｄａｔｉｏｎｕｎｄｅｒａｉｒ－ｗａｔｅｒｉｎｌｅｔｔｗｏ ｐｈａｓｅｆｌｏｗｃｏｎｄｉｔｉｏｎｓ［Ｊ］．Ｈｏｕｉｌｌｅｂｌａｎｃｈｅ，２０１８（３）：４１－４８．［７］ 韩路长．多相反应器内流体颗粒传质、破裂分散和聚并行为的研究［Ｄ］．湘潭：湘潭大学，２０１０．［８］ 戈振国．基于ＣＦＰ－ＰＢＭ耦合模型的离心泵气液两相流动特性研究［Ｄ］．西安：西安理工大学，２０１７．［９］ 文键，王斯民，厉彦忠．基于ＭＵＳＩＧ模型的低温流体过冷沸腾数值模拟［Ｊ］．化学工程，２０１０，３８（１１）：２９－３３．ＷＥＮＪｉａｎ，ＷＡＮＧＳｉｍｉｎ，ＬＩＹａｎｚｈｏｎｇ．ＮｕｍｅｒｉｃａｌｓｉｍｕｌａｔｉｏｎｆｏｒｓｕｂｃｏｏｌｅｄｂｏｉｌｉｎｇｐｒｏｃｅｓｓｏｆｌｏｗｔｅｍｐｅｒａｔｕｒｅｆｌｕｉｄｗｉｔｈＭＵＳＩＧｍｏｄｅｌ［Ｊ］．Ｃｈｅｍｉｃａｌｅｎｇｉｎｅｅｒｉｎｇ（Ｃｈｉｎａ），２０１０，３８（１１）：２９－３３．（ｉｎＣｈｉｎｅｓｅ）［１０］ ＬＵＯＨ，ＳＶＥＮＤＳＥＮＨＦ．Ｔｈｅｏｒｅｔｉｃａｌｍｏｄｅｌｆｏｒｄｒｏｐａｎｄｂｕｂｂｌｅｂｒｅａｋｕｐｉｎｔｕｒｂｕｌｅｎｔｄｉｓｐｅｒｓｉｏｎｓ［Ｊ］．Ａｉｃｈｅｊｏｕｒｎａｌ，１９９６，４２（５）：１２２５－１２３３．［１１］ ＰＲＩＮＣＥＭＪ，ＢＬＡＮＣＨＨＷ．Ｂｕｂｂｌｅｃｏａｌｅｓｃｅｎｃｅａｎｄｂｒｅａｋ ｕｐｉｎａｉｒ ｓｐａｒｇｅｄｂｕｂｂｌｅｃｏｌｕｍｎｓ［Ｊ］．ＡＩＣｈＥｊｏｕｒｎａｌ，１９９０，３６（１０）：１４８５－１４９９．［１２］ 崔强磊．入流含气条件下离心泵内部不稳定流动及其诱导特性［Ｄ］．镇江：江苏大学，２０１９．（责任编辑 盛杰）３３７。
液—固水力旋流器两相流动数值模拟研究进展液—固水力旋流器广泛应用于各个行业，如石油化工、选矿、造纸、医药卫生、环境保护和食品等。
近年来，随着各工业领域的不断发展，对液固分离技术与装备提出了新的挑战和更高的要求同时，为了适应当今不断高涨的降低能耗的要求，目前迫切需要开发出高分离速度、高脱水度、高分离精度的高性能液—固分离技术然而，由于实验条件的限制，单纯通过实验来研究旋流器的性能不仅周期长而且费用高，如果辅助以理论分析计算和流场模拟等方法来研究旋流器内部流体流动规律，以及结构尺寸变化对分离性能和压力特性的影响等，则可缩短研究周期和实验费用，具有重要的理论研究和工程应用价值。
液—固水力旋流器结构及工作原理图1所示为典型的液—固旋流器，主要由进料口、溢流口、底流口和器壁组成。
其工作原理为，当液体高速旋转受离心力作用时，轻相向轴心迁移，从溢流口排出，重相向壁面运动，由底流口排出，从而实现轻相( 液) 和重相( 固) 分离的目的。
图1 典型液—固旋流器两相流场数值模拟研究进展液—固水力旋流器中液固分离的问题在数值模拟中定义为液体和固体不相容的两相流问题。
由于液固两相之间的相互作用和每一相的运动，传热传质和反应等的影响，颗粒相的模拟基本分为两类：一类是Euler方法，该方法除将流体作为连续介质外，把颗粒群也作为拟连续介质或拟流体，设其在空间有连续的速度和温度分布及等价的输运性质( 粘性扩散导热等) ；另一类是Lagrange方法，该方法把流体作为连续介质，而将颗粒群看作离散体系，并以此来探讨颗粒动力学颗粒轨道等基于这两种方法，研究者采用了不同的模型对旋流器内的两相分离过程进行了模拟研究。
K T Hsien和R K Rajamani( 1991)根据颗粒的受力平衡，用代数逼近法求出固体颗粒的滑移速度和轨迹，P He，M Salcudean和I S Gartshore(1997 )分别用二维和三维模型计算了旋流器的分离效率。
气液两相流模型及相应控制方程1. 分相流模型一般把两相看成分开的两股流体流动，分别按单相流处理，同时计及相间作用，然后把两相方程相加，得出两相流的基本方程。
以下为分相流动模型下的一元两相流动基本控制方程。
(1) 连续性方程： 气相：g g g g g f A f w Am tzρρδ∂∂+=∂∂液相：(1)(1)l g g g l f A l f w Am t z ρρδ∂-∂-+=-∂∂ 两相：0A GAt zρ∂∂+=∂∂ 其中：m δ为控制体单位长度的质量交换率，即考虑相变。
真实密度：(1)g g l g f f ρρρ=+- 质量流速: (1)g g g g l l G f w f w ρρ=+- (2) 动量方程：因为相间有摩擦，相界面上作用有摩擦力i τ，相间有质量交换，则有动量交换，若相界面周长为i P ，界面流速为i w ，则各自的动量方程为： 气相：21sin g og g gg g g g g g g i f w f w A p m f P gf w tA z z A Aρτρρθ∂∂∂∂+=---+∂∂∂ 液相：2(1)(1)1(1)(1)sin ol l g l l g l g l l g if w f w A p m f Pg f w t A z z A Aτρρρθ∂-∂-∂∂+=-----+∂∂∂两相：22(1)1sin l g l g g g o A f w f w P G P g t A z z Aρρτρθ⎡⎤∂-+∂∂∂⎣⎦+=---∂∂∂其中：o og g ol l P P p τττ=+ 且总压等于分压之和。
（3）能量方程：依照单相流动方程的方法的推法，即在控制体上应用热力学第一定律，再考虑相间的相互作用，即相间摩擦耗功，相变递能和界面传热，得各自的能量方程为： 气相能量方程：22()(())g g g g g g g g g Af U w w Af U w dz dz tzρρ⎡⎤⎡⎤∂+∂+⎣⎦⎣⎦+=∂∂21sin 2g g g g g g i i i i i i pAf w dQ dz Af w gdz Pw dz q Pdz mw dzZρθτδ⎡⎤∂⎣⎦----+∂/*注意与单相流能量方程的区别，热力学能第一定律只关注于内能，机械能与热能数量上的守恒，动量定理则能清楚的看出各种力对流动特性的影响，比如在能量方程中你看不出粘性耗散对流动的影响*/g dQ 单相流中推导的气体与外界，不包括液体之间的热量交换。
两相流：通常把含有大量固体或液体颗粒的气体或液体流动称为两相流；其中含有多种尺寸组颗粒群为一个“相”，气体或液体为另一“相”，由此就有气—液，气—固，液—固等两相流之分。
两相流的研究：对两相流的研究有两种不同的观点：一是把流体作为连续介质，而把颗粒群作为离散体系；而另一是除了把流体作为连续介质外，还把颗粒群当作拟连续介质或拟流体。
引入两种坐标系：即拉格朗日坐标和欧拉坐标，以变形前的初始坐标为自变量称为拉格朗日Langrangian 坐标或物质坐标；以变形后瞬时坐标为自变量称为欧拉Eulerian 坐标或空间坐标。
离散相模型FLUENT在求解连续相的输运方程的同时，在拉格朗日坐标下模拟流场中离散相的第二相；←离散相模型解决的问题：煤粉燃烧、颗粒分离、喷雾干燥、液体燃料的燃烧等；←应用范围：FLUENT中的离散相模型假定第二相体积分数一般说来要小于10-12%（但颗粒质量承载率可以大于10-12%，即可模拟离散相质量流率等/大于连续相的流动）；不适用于模拟在连续相中无限期悬浮的颗粒流问题，包括：搅拌釜、流化床等；←颗粒-颗粒之间的相互作用、颗粒体积分数对连续相的影响未考虑；湍流中颗粒处理的两种模型：Stochastic←Tracking，应用随机方法来考虑瞬时湍流速度对颗粒轨道的影响；Cloud Tracking，运用统计方法来跟踪颗粒围绕某一平均轨道的湍流扩散。
通过计算颗粒的系统平均运动方程得到颗粒的某个“平均轨道”多相流模型FLUENT中提供的模型：VOF模型(Volume of Fluid Model)←混合模型(Mixture Model)←←欧拉模型(Eulerian Model)VOF模型(Volume of Fluid Model)← VOF模型用来处理没有相互穿插的多相流问题，在处理两相流中，假设计算的每个控制容积中第一相的体积含量为α1，如果α1=0，表示该控制容积中不含第一相，如果α1=1，则表示该控制容积中只含有第一相，如果0<α1<1，表示该控制容积中有两相交界面；← VOF方法是用体积率函数表示流体自由面的位置和流体所占的体积，其方法占内存小，是一种简单而有效的方法。
1. 引言在科学研究和工程领域，液滴在流场中的破碎现象一直备受关注。
这个主题涉及到流体力学、计算物理学以及工程应用等多个领域，涉及的问题也非常复杂。
近年来，随着计算机模拟技术的不断发展，光滑粒子法（Smoothed Particle Hydrodynamics，SPH）作为一种数值模拟方法，在研究液滴破碎现象中发挥了越来越重要的作用。
本文将围绕液滴在流场中破碎的光滑粒子法数值模拟展开深入探讨。
2. 液滴破碎的现象和意义液滴在流场中破碎是一种普遍存在的现象，例如在空气动力学、生物工程和化工等领域都有着广泛的应用。
研究液滴破碎的机理和规律，有助于我们更好地理解复杂的流体行为，优化工程设计和预测事故发生。
进行液滴破碎的数值模拟研究具有重要的理论和实际意义。
3. 光滑粒子法的基本原理光滑粒子法是一种基于流体微元的数值模拟方法，其基本原理是将流体划分为离散的粒子，通过对粒子之间的相互作用力进行数值计算，来模拟流体的运动。
SPH方法由于其适用于连续介质的性质、处理自由表面流体和多相流等优点而在流体力学领域得到了广泛应用，特别是在液滴破碎的数值模拟中显示出了独特的优势。
在液滴破碎的光滑粒子法数值模拟中，首先需要建立流场的数值模型。
通过对流场的一系列参数进行离散化，可以得到一系列粒子，并利用连续介质力学和粒子内外相互作用力来模拟液滴在流场中的破碎过程。
在模拟过程中，需要考虑表面张力、粘性、重力和碰撞等因素的影响，并结合数值计算方法对其进行全面的模拟分析。
5. 液滴破碎的数值模拟结果和分析通过光滑粒子法进行液滴破碎的数值模拟，可以得到液滴形状、体积、速度、压力等参数随时间的变化规律。
通过对模拟结果的分析，可以深入理解液滴在复杂流场中的运动规律和破碎机理，为液滴破碎现象提供更加深刻的理论解释和工程应用指导。
6. 总结和展望通过本文对液滴在流场中破碎的光滑粒子法数值模拟的探讨，可以更好地理解这一复杂现象的物理本质和数值模拟方法的原理。
气液两相段塞流中双流体模型的分析赵志勇 王树立 彭 杰抚顺石油学院机械工程分院储运工程系 辽宁省抚顺市 113001【摘要】采用一维瞬态双流体模型描述了倾斜管道中气液两相段塞流的流动情况 。
提出了段塞单元的简化水力模型 。
认为段塞流是由一系列段塞单元构成 ,该段塞单元由液相段塞区和膜区 两部分组成 ;通过对段塞单元的分析 ,采用新方法来计算段塞频率 、段塞单元及液相段塞的长度 ;根 据其流动特点推导出表征气液两相间相互作用及气液两相与管壁间相作用的关联式 。
本文在处理 相间曳力系数时摒弃以往采用经验常数的方法 ,而是在稳态段塞流动量守恒方程的基础上推导出 新型的曳力系数表达式从而达到封闭模型的目的 。
关键词 : 气液两相 瞬态 双流体模型 段塞流A Two 2Fl u id Model U sed In G a s an d L i quid Two 2Pha s e Sl u g Fl o wZha o Z hiyong Wang Shul i P eng JieDepart m ent of Mechanical engineering , Fushun Pet r oleum Instit u te , Fushun Liao n ing 113001Abstract : A o ne d imensio n t ransient t w o 2f luid mo del was emplo yed to describe t he slu g f lo w in a tilting t ube which carries t he gas 2liqu id t w o p hase mixt ure. It is su ppo sed t hat slug flo w be co nstit ut e of series of slu g unit s which co ntain t w o p art s : t h e liquid slug zo n e and film zo n e. t h e co n stit u tive relatio n shi p s for int eractio n bet w een gas and liquid ,and also t h at of t w o p h ase and t h e t u be wall were derived f r o m it ’s mo vement p r op ert y o n t h e basis of mo m ent u m balance. A co n serva 2 tive equatio n ,but no t f r o m t h e empirical dat a .K ey Words : G as 2Liquid t w o p h ase ; Transient ; T w o 2Flu id mo del ; S lug f lo w .〔2〕无压波模型( N FM ) 。
气液两相流理论分析模型及实验平台介绍摘要：介绍了气液两相流的理论分析模型，均相流模型、分相流模型、漂移流模型、双流体模型，各模型的特点及适用条件，为实验研究提供理论基础。
并介绍了国家化工设备质量监督检验中心换热实验平台，为气液两相流实验研究提供支撑。
关键词：气液两相流；两相流模型；实验平台伴随着管壳式换热器的普遍应用，两相流动及沸腾换热特性的研究越来越得到各国学者们的重视。
两相流动及沸腾广泛存在于各个工业领域中，如工业锅炉设备中，水经过软化后通过电加热或者化石燃料燃烧加热沸腾，气化为水汽两相状态，过热的饱和蒸汽经过汽包，送至各类用户；产热电站及动力站也是通过各类热源的作用，使得原始介质发生相变产生两相或者多相介质直接或者间接作为动力及发电源动力；汽轮机、核反应堆堆芯、蒸汽发生器等等同样伴随着各类气液两相或者多相流动及沸腾换热现象。
一、气液两相流理论分析模型（一）均相流模型均相流模型是最简单的物理模型，结构示意图如图1。
将气液两相流混合物看做是一种均匀介质，混合物的物性参数即两相介质各自参数的平均值。
运用一般流体力学的研究方法，对气液混合物的平均流动特性进行确定和描述。
该模型基于以下几项基础的假设：1.气液两相的流速相等，故滑动比，滑差，流动密度都相等，把气液两相流体看作均匀流。
2.气液两相混合物处于热力学平衡状态。
即气液两相的温度是相等的，且不存在热量传递。
3.两相流的摩擦系数基于单相流摩擦系数，且气液两相混合物物理特性遵循等熵变化。
这种均相模型一般适用于高速高压条件下的泡状流和雾状流。
这是由于在一般实际情况下，两相流速是不可能相同的。
只有在高速高压下两相流质混合才更加均匀。
由于假设两相之间速度是相同的，其与实验计算值误差较大，随着质量流量的减小偏差逐渐增大。
所以准确确定气液两相混合物的平均物理特性对于选择模型十分重要。
（二）分相流模型分相流模型的结构示意图如图2。
分相流模型将气液两相流的两相分别看作连续的有自身平均流速的单独流体，即气液介质有各自的物性参数，当气液两相流速相等时就成为了均相流模型。
1.多相流动模式我们可以根据下面的原则对多相流分成四类：•气-液或者液-液两相流：o 气泡流动：连续流体中的气泡或者液泡。
o 液滴流动：连续气体中的离散流体液滴。
o 活塞流动: 在连续流体中的大的气泡o 分层自由面流动：由明显的分界面隔开的非混合流体流动。
•气-固两相流：o 充满粒子的流动：连续气体流动中有离散的固体粒子。
o 气动输运：流动模式依赖诸如固体载荷、雷诺数和粒子属性等因素。
最典型的模式有沙子的流动，泥浆流，填充床，以及各向同性流。
o 流化床：由一个盛有粒子的竖直圆筒构成，气体从一个分散器导入筒内。
从床底不断充入的气体使得颗粒得以悬浮。
改变气体的流量，就会有气泡不断的出现并穿过整个容器，从而使得颗粒在床内得到充分混合。
•液-固两相流o 泥浆流：流体中的颗粒输运。
液-固两相流的基本特征不同于液体中固体颗粒的流动。
在泥浆流中，Stokes 数通常小于1。
当Stokes数大于1 时，流动成为流化（fluidization）了的液-固流动。
o 水力运输: 在连续流体中密布着固体颗粒o 沉降运动: 在有一定高度的成有液体的容器内，初始时刻均匀散布着颗粒物质。
随后，流体将会分层，在容器底部因为颗粒的不断沉降并堆积形成了淤积层，在顶部出现了澄清层，里面没有颗粒物质，在中间则是沉降层，那里的粒子仍然在沉降。
在澄清层和沉降层中间，是一个清晰可辨的交界面。
•三相流(上面各种情况的组合)各流动模式对应的例子如下：•气泡流例子：抽吸，通风，空气泵，气穴，蒸发，浮选，洗刷•液滴流例子：抽吸，喷雾，燃烧室，低温泵，干燥机，蒸发，气冷，刷洗•活塞流例子：管道或容器内有大尺度气泡的流动•分层自由面流动例子：分离器中的晃动，核反应装置中的沸腾和冷凝•粒子负载流动例子：旋风分离器，空气分类器，洗尘器，环境尘埃流动•风力输运例子：水泥、谷粒和金属粉末的输运•流化床例子：流化床反应器，循环流化床•泥浆流例子: 泥浆输运，矿物处理•水力输运例子：矿物处理，生物医学及物理化学中的流体系统•沉降例子：矿物处理2. 多相流模型FLUENT中描述两相流的两种方法:欧拉一欧拉法和欧拉一拉格朗日法，后面分别简称欧拉法和拉格朗日法。
液液两相流流体破碎模型
液一液两相中的液滴变形和破碎现象不论在自然界还是在生产实践中都广泛存在，如石油行业中，会遇到油和水两相混合流动的现象。
在油水两相流动过程中，常常会形成油。
水分散体系，其中分散相粒径大小和分布对于油水混合物的输运、检测以及分离等过程都会产生重要的影响，而分散相的粒径大小和分布又与液滴变形破碎过程密切相关。
此外，液滴的变形破碎过程也广泛存在于化工、环保等领域，如液液萃取、多相反应、悬浮聚合及乳状液的制备等化工过程，这些过程速率取决于二相间的相际面积，研究液滴变形破碎对于增加相际面积，加快反应速率提高效率等具有举足轻重的作用。
由此可见，认识和掌握较大液滴的变形、破碎过程机理及其规律无疑是非常必要的。
要准确预测系统中液滴群的运动和传质过程，最基本的出发点就是对单个液滴的流体力学行为的预测。
同时，对单个液滴流体力学行为的准确把握为我们提供了一个理解更为复杂的实际多相流系统的基础。
液滴运动规律的研究越来越受到国内外的关注，虽然理论方法和实验手段都取得了一定的进展，但仍有其本身的局限性。
随着CFD的发展，数值模拟成为探索液滴运动规律的重要手段。
Rallison 和Acrivos[1]首先将边界积分方法应用于液滴变形数值模拟，该方法的主要优点是使用势函数将二维问题降为一维问题，精度比较高，但由于其数值稳定性较差，只能用来模拟变形不大的液滴运动。
由Hirt和Nichols[2]提出的VOF方法被广泛用于两相流的数值模拟中，流体体积分数概念的引入大大简化了捕捉界面的计算，然而由VOF只能得到控制单元中的流体体积分数，要得到物质界面还需要进行界面重构，不同的重构方法会有不同的效果。
Osher和Sethia[3]提出的水平集(Levelset)方法是目前处理、追踪物质界面效果较好的一种方法，通过引入水平集函数的概念自动捕捉界面的拓扑变化，在处理复杂结构变化方面优势明显。
尤学一，刘伟[4]采用VOF法追踪了重力流液液和气液相界面的迁移，发现VOF法可很好地追踪强非线性、大形变的相界面随时间的变化。
重点讨论了不同相密度比、黏性比条件下，液液相界面和气液相界面的运动形式，验证了液液相界面主要受Kelvin-Helmholtz不稳定性控制。
结果显示：液相密度比越大，液液相界面随时间变化越快，气液交界面不再保持水平位置，而以摆动的形式存在；在黏性比小于5时，黏性比不太影响液液相界面随时间变化。
熊燃华[5]在数值求解两相混合的质量和动量守恒方程基础上、运用VOF(volume of fluid )模型求解每个网格单元的体积分数离散方程、然后用界面重构技术追踪液滴与环境流体界面,从而模拟了液滴整个演变过程。
数值计算结果虽然在后期破碎过程与实验结果有所区别，但也反映出这几种模式的演变特
点；根据数值得到的流场结构，分析了实验推理的液滴演变机理的合理性，而且进一步深化了液滴演变机理解释。
高一娟[6]基于多相流理论，应用Level-Set方法，建立了微尺度两相流相界面跟踪模型，实现了微槽道中液液两相流的相界面追踪的动态模拟。
建立的两相流微流动数值模型能很好的实现界面的动态跟踪和流场模拟，为微槽道内多相流流动模拟打下了好的基础。
董贺飞，张德良，杨国伟[7]采用PFV算法，并结合LevelSet方法，对化工系统中液液两相系统问题进行数值模拟，精确地捕捉两相流动中的界面信息。
并以典型的液滴的坠落过程为例，得到了液滴的坠落过程中界面变形和分裂的全过程。
数值模拟和已发表论文结果十分吻合，证实了采用摄动有限体积新算法数值模拟液液两相系统中液滴的坠落过程是成功的，为两相流动的数值模拟提供了一种新途径。
国内相关报道较多的是VOF模型，PFV算法是比较不错的模型。
参考文献
[1]GaoZhi，Y angGuowei．Perturbation finite volume method for convective—difusion integral euqation[J]．Acta Mechanica Sinica，2004，20(6)：580—590．
[2]Rallison J M,Acrivos A.A numericla study of the deformation and burst of a viscous drop in an extensional flow[J]．JfluidMech，1978，89：191—200．
[3]Hitr C W,Nicobls B D．V olume of fluid(VOF)method for the dynamics of free boundaries [J]．JComputPhys，1981，39：201—225．
[4]尤学一，刘伟．两相流相界面迁移的数值模拟[J]．水动力学研究与进展，2006，21(6)：724-729。
[5]熊燃华．液-液两相介质中液滴在瞬时来流作用的演变过程研究[D]．中国科学技术大学，2010。
[6]高一娟．微通道内两相流界面追踪的数值模拟[D]．天津大学，2009。
[7]董贺飞，张德良，杨国伟．摄动有限体积算法及其在两相流中应用．计算物理，2009，26(6)：857—864。
