基于两相流EWF模型蒸发式空冷器数值模拟

求解两相蒸发和冷凝问题的气液相变模型孙东亮;徐进良;王丽【摘要】A vapor-liquid phase change model on the basis of the volume of fluid (VOF) method in FLUENT was proposed to solve the two-phase boiling or condensation problems. This model is suitable for the case where one is unsaturated phase and the other is saturated phase. In the model proposed here, the unsaturated-phase thermal conductivity and specific heat are actual physical parameters; the saturated-phase thermal conductivity is assumed to be zero and the saturated-phase specific heat is assumed to equal the unsaturated-phase specific heat; the interfacial mass-transfer rate is determined only by the unsaturated phase. In addition, our phase change model was verified by the analytical solutions of the one-dimensional Stefan problem and the accurate results of the two-dimensional boiling problem. The work will lay a foundation for the application of this model.%为了更精确简单地求解两相流中的蒸发和冷凝问题,基于FLUENT中的流体体积(VOF)方法提出了一种气液相变模型,该模型适应于两相中一相为非饱和相,即处于过热或过冷状态,另一相为饱和相,即处于饱和状态.该气液相变模型中:非饱和相导热系数和比热容为真实的物性参数；假设饱和相导热系数等于0,饱和相比热容等于非饱和相比热容；界面处相变率仅由非饱和相决定,最后,通过一维Stefan问题相界面位置的分析解和二维膜态沸腾气相体积比的精确解,验证了该相变模型的精确性和可行性.文中工作为该相变模型的推广应用奠定了基础.【期刊名称】《西安交通大学学报》【年(卷),期】2012(046)007【总页数】5页(P7-11)【关键词】两相流;气液相变模型;流体体积方法【作者】孙东亮;徐进良;王丽【作者单位】华北电力大学新能源电力系统国家重点实验室,102206,北京;华北电力大学低品位能源多相流与传热北京市重点实验室,102206,北京;华北电力大学低品位能源多相流与传热北京市重点实验室,102206,北京【正文语种】中文【中图分类】TK124模拟气液两相流动与传热的数值方法主要分为颗粒轨道模型、双流体模型和运动界面追踪技术，其中运动界面追踪技术最能准确、全面地反映气液两相间的信息.目前，FLUENT是最广泛使用的数值模拟软件，它采用流体体积（VOF）方法追踪相间的运动界面，但该方法仍缺乏准确的气液相变模型，无法准确模拟相间的蒸发和冷凝现象.在已有的相变模型中，Lee［1］提出的模型得到了最广泛的应用，其具体表达式为蒸发过程式中：˙mg 和˙ml 分别为气相和液相的相变率，kg／（m3·s）；r为控制相变强度的因子，该变量的设定依据为使包含相界面的计算网格的温度值维持在饱和温度，针对不同问题该变量具有不同的取值.De Schepper等［2］采用该模型对换热管内的沸腾过程进行了模拟，r设为0.1s－1；Alizadehdakhel等［3］对热管内的蒸发冷凝现象进行了计算，r也设为0.1 s－1；文献［4-5］中r均设为100s－1.在蒸发和冷凝问题中，相界面温度通常设定为饱和温度Tsat.根据傅里叶定律，可计算出相界面两侧热流密度的界跃值由于推导过程的假设不合理，导致在部分问题中该模型的计算结果与实际物理过程具有较大的偏差.例如在一维Stenfan问题中，界面位置随时间的变化规律与液体的导热系数λl无关，具体见其分析解表达式（28），而在式（4）中却包含了λl的信息，导致模拟结果与λl具有相关性.目前，许多作者［7－8］基于 VOF方法开发了求解相变问题的程序代码，这些方法的核心内容就是如何精确求解相界面两侧的热流密度，即式（3）中的－λl∂Τ／∂n｜l和－λg∂Τ／∂n｜g.虽然这些方法可以精确地模拟蒸发和冷凝问题，但求解过程较为复杂，需要获取相界面位置和相界面附近网格的信息，很难在FLUENT中推广.本文基于FLUENT中的VOF方法，提出了一种求解两相蒸发和冷凝问题的精确简单的气液相变模型，该模型适应于两相中一相为非饱和相，即处于过热或过冷状态，另一相为饱和相，即处于饱和状态.本文还通过两个经典的相变问题验证了该相变模型的精确性和可行性，表明该模型可应用于微通道内的蒸发和冷凝问题.VOF方法是Hirt等［9］在1981年首先提出的，它引入了流体体积分数α的概念，各相流体体积分数之和等于1，即该相变模型主要包括以下3部分.（1）非饱和相导热系数λuns和比热容Cp，uns为真实的物性参数.（2）假设饱和相导热系数λs＝0，比热容Cp，s＝Cp，uns.饱和相各点温度相等，不存在导热，因此饱和相导热系数的大小对结果没有影响.本模型从数值计算角度出发，强制设定λs＝0，从而可以强制保证饱和相内温度保持不变.同时，饱和相温度保持不变，Cp，s的大小对计算结果不产生影响，因此为了计算方便，本模型设定Cp，s＝Cp，uns.（3）求解界面处相变率˙m.该模型针对一相为非饱和相、另一相为饱和相的情况，故式（3）可写成式中：‖qs‖为相界面上热流密度界跃值，W／m2；n指向饱和相侧.与将表面张力由表面力转化为体积力的CSF模型相似，将表面上的‖qs‖转化为相界面两侧有限宽度内单位体积的热量值qv，即式中：L表示相界面两侧有限宽度值，m；qv表示单位体积的热量值，W／m3.假设相界面附近非饱和相温度呈线性分布，则可构建出通过式（21）～（26），即可验证式（20）的合理性.最后，基于式（20），可以得出界面处的相变率将本文提出的气液相变模型应用于一维Stefan问题和二维膜态沸腾问题，以验证模型的精确性和可行性，其中采用Youngs的PLIC方法捕捉相界面，采用SIMPLE算法对速度和压力进行耦合求解.一维Stenfan问题是带有相变的两相流动的经典问题，如图2所示.当给气体侧的固体壁面上施加一个较高的温度边界条件时，气体侧就会产生温度梯度，热量向液体部分传递，由于液体部分已经处于饱和状态，因此传递到液体部分的热量会使得液体部分开始气化，在相界面处产生相变.该问题相界面位置δ（t）的分析解［11］为该问题中非饱和相为过热气体，ρg＝0.001，0.01，0.1kg／m3，λg＝0.005W／（m·K），Cp，g＝200 J／（kg·K）；饱和相为饱和液体，ρl＝1kg／m3，根据本文提出的模型，假设λl＝0，Cp，l＝Cp，g＝200J／（kg·K）；hfg＝10kJ／kg，壁面温度与饱和温度的差值ΔT＝Tw－Tsat＝10K.通过网格独立性考核，选取宽度为0.001m的网格进行计算，时间步长为1×10－4 s.图3给出了一维Stefan问题中不同密度比率下气液界面位置随时间的变化曲线，从图中可以看出，本文相变模型的模拟结果与分析解相一致，验证了模型的精确性和可行性.图4给出了二维膜态沸腾示意图.该问题中非饱和相为过热气体，ρg＝5kg／m3，动力黏度μg＝0.005Pa·s，λg＝1W／（m·K），Cp，g＝200J／（kg·K）；饱和相为饱和液体，ρl＝200kg／m3，μl＝0.1Pa·s，根据本文提出的模型，假设λl＝0，Cp，l＝Cp，g＝200J／（kg·K）.hfg＝10kJ／kg；σ＝0.1N／m，g＝9.81m ／s2.图4中λ0为Taylor不稳定性波长由于流动结构的对称性，本文选取一半区域作为计算区域，宽度为λ0／2，高度为λ0.边界条件：底部壁面速度为无滑移边界条件，温度为定壁温边界条件，壁面温度比饱和温度高5K，即ΔT＝Tw－Tsat＝5K；左右均为对称边界条件；上部为压力出口边界条件.初始条件：速度初始条件为u＝v＝0 m／s；液体侧初始温度为饱和温度，气体侧初始温度分布为沿着气液界面到底部壁面呈线性分布；初始气液界面位置为通过网格独立性考核，计算区域网格数选为32×64，同时相界面附近进行局部加密，最小网格为初始网格的1／4，时间步长为1×10－5 s.图5显示了初始气泡形成过程中不同时刻的气液界面分布.图6给出了该过程中气相体积比（当前时刻的气相体积／初始时刻的气相体积）随时间变化的曲线.图6中精确解为Guo等［8］基于自己开发的程序代码得出的计算结果，从图中可以看出，本文相变模型的模拟结果与精确解完全吻合，进一步验证了模型的精确性和可行性.本文基于FLUENT中的VOF方法，提出了一种求解两相蒸发和冷凝问题的精确简单的气液相变模型，并通过一维Stefan问题和二维膜态沸腾问题验证了该模型的精确性和可行性.该气液相变模型主要针对两相中一相为非饱和相、另一相为饱和相的情况，主要包括以下3部分：（1）非饱和相导热系数λuns和比热容Cp，uns为真实的物性参数；（2）假设饱和相导热系数λs＝0，比热容Cp，s＝Cp，uns；（3）界面处相变率在该气液相变模型的基础上，作者正在开发可以应用于两相同时处于非饱和状态的相变模型，以进一步扩大模型的应用范围.【相关文献】［1］ LEE W H.A pressure iteration scheme for two-phase flow modeling［M］.Washington，USA：Hemisphere Publishing，1980.［2］ DE SCHEPPER S C K，HEYNDERICHX G J，MARIN G B.Modeling the evaporation ofa hydrocarbon feedstock in the convection section of a steam cracker［J］.Computersand Chemical Engineering，2009，33（1）：122-132.［3］ ALIZADEHDAKHEL A，RAHIMI M，ALSAIRAFI A A.CFD modeling of flow and heat transfer in a thermosyphon［J］.International Communications in Heat and Mass Transfer，2010，37（3）：312－318.［4］ YANG Z，PENG X F，YE P.Numerical and experimental investigation of two phase flow during boiling in a coiled tube［J］.International Journal of Heat and Mass Transfer，2008，51（5／6）：1003－1016［5］ FANG C，DAVID M，ROGACS A，et al.Volume of fluid simulation of boiling two－phase flow in a vaporventing microchannel［J］.Frontiers in Heat and Mass Transfer，2010，1（1）：1－11.［6］毛文彬.微流控系统中混合及汽液相变传热的数值模拟［D］.广州：中国科学院广州能源研究所，2009.［7］ WELCH S W J，WILSON J.A volume of fluid based method for fluid flows with phase change［J］.Journal of Computational Physics，2000，160（2）：662－682.［8］ GUO D Z，SUN D L，LI Z Y，et al.Phase change heat transfer simulation for boiling bubbles arising from a vapor film by VOSET method［J］.Numerical Heat Transfer：Part A，2011，59（11）：857－881.［9］ HIRT C W，NICHOLS B D.Volume of fluid（VOF）method for the dynamics of free boundary［J］.Journal of Computational Physics，1981，39（1）：201－225.［10］BRACKBILL J U，KOTHE D B，ZEMACH C.A continuum method for modelingsurface tension ［J］.Journal of Computational Physics，1992，100（2）：335－354. ［11］ ALEXIADES V，SOLOMON A D.Mathematical modeling of melting and freezing processes ［M］.Washington，USA：Hemisphere Publishing，1993.。

针对飞行器两相流问题的数值模拟研究飞行器两相流问题是飞行器研究中不可避免的一个难题。

随着飞行器速度的提高，乃至于超声速、高超声速等领域的研究，这一问题显得更加尖锐。

数值模拟是研究这一问题的重要手段，本文将讨论针对飞行器两相流问题的数值模拟研究。

一、飞行器两相流问题的研究背景随着航空航天工业的不断发展，以及人们对速度、效率等方面的日益追求，飞行器的设计、制造、使用等方面也发生了很大的变化。

尤其是在超声速、高超声速、喷气推进等领域，飞行器对气动力学、热力学等方面的要求更加苛刻，对设计和材料的要求更加高。

两相流问题是飞行器研究过程中的一个重要问题。

简而言之，就是指流动体系中存在两种物质时，其中一种物质的存在对另一种物质的流动形态产生影响，甚至可能引起凝结、相变等一系列物理过程。

由于飞行器研究中常常涉及到高速、高温、高压等条件，这一问题显得尤为重要和复杂。

二、数值模拟在飞行器两相流问题研究中的应用数值模拟是一种重要的研究手段。

通过建立数学模型，运用计算机等工具对物理问题进行数值求解，从而得到物理过程的定量描述和分析。

在飞行器两相流问题研究中，数值模拟发挥了重要作用。

1. 建立数学模型在进行数值模拟之前，必须建立一个合适的数学模型。

建立数学模型是将实际问题抽象成为符号和公式的过程。

对于飞行器两相流问题，建立数学模型可以将问题分解为宏观平均场和微观实体物质之间的相互作用。

通过求解数学模型中的方程组，可以得到两相流过程中的速度、密度、温度等重要参数。

2. 选择数值方法在建立数学模型之后，需要选择合适的数值方法对方程组进行求解。

数值方法是将微积分中的连续问题抽象成离散问题，然后通过计算机等工具对离散化问题进行求解的一种方法。

在针对飞行器两相流问题的数值模拟研究中，常用的数值方法包括欧拉法、拉格朗日法、伽辽金法等。

每种数值方法都有其适用范围和优劣点，需要根据具体情况进行选择。

3. 求解计算通过建立数学模型和选取数值方法，就可以开始进行数值模拟。

基于CFD两相流模型的涡旋压缩机抗液击性能预测方法林江波;郑康;王柳
【期刊名称】《制冷与空调》
【年(卷),期】2024(24)1
【摘要】当涡旋压缩机在某些特殊工况运行时,液态制冷剂会进入压缩机,进而部分液态制冷剂进入涡旋腔,对涡旋壁和十字环产生很大的冲击力。

当进入的液态制冷剂量超过涡旋腔所能承受的极限,对涡旋和十字环造成损坏并导致压缩机失效。

在新压缩机产品开发过程中,传统的方式是通过试验方法去验证压缩机的抗液击性能,但这种方法需要组装样机并在实验室进行测试,时间长且成本高。

文章提出了一种基于两相流模型预测涡旋压缩机抗液击性能的方法,通过该方法,在新产品设计的早期阶段无需试验就可以预测涡旋压缩机的抗液击性能,能够有效地降低项目成本且加快产品上市时间。

【总页数】6页(P69-74)
【作者】林江波;郑康;王柳
【作者单位】丹佛斯(天津)有限公司
【正文语种】中文
【中图分类】TH4
【相关文献】
1.基于CFD的液固两相流冲刷腐蚀预测研究
2.基于均质两相流的滚动转子压缩机径向泄漏预测模型
3.基于LabVIEW的涡旋压缩机液击电流检测方法与实现
4.基于
DSMC-CFD方法的气固两相流冲蚀预测研究5.基于CFD-DPM耦合方法的多级串联弯管内液固两相流冲蚀特性研究
因版权原因，仅展示原文概要，查看原文内容请购买。

y型微通道气液两相流的数值模拟简
微通道中的气液两相流是一种复杂的流动现象，其数值模拟可以通过多相流模型来求解。

下面是一种简单的数值模拟方法：
1. 网格划分：根据微通道的几何形状，将计算域划分为若干个小单元，形成网格。

可以使用结构化网格或非结构化网格，根据实际情况选择合适的网格类型。

2. 方程建立：建立气液两相流动的守恒方程，包括质量守恒、动量守恒和能量守恒方程。

对于气液两相流动，需要考虑两相间的相互作用，使用数学模型描述气液两相的物理过程。

3. 物理模型：选择合适的两相流模型。

常见的两相流模型包括欧拉模型、VOF模型、多流体模型等。

根据实际情况选择适应的模型。

4. 边界条件：根据实际情况设置合适的边界条件，包括压力、速度和温度等。

边界条件的选择对模拟结果有较大影响，需要根据实际情况进行合理设定。

5. 数值求解：使用数值方法对守恒方程进行离散化，得到离散方程。

常用的数值方法有有限差分法、有限体积法和有限元法等。

通过迭代求解离散方程，得到气液两相流场的数值解。

6. 边界处理：处理流场的边界，使得计算结果满足物理约束条件。

边界处理包括边界设置、边界条件的施加和边界修正等。

7. 结果分析：对计算结果进行分析与评估，包括流速分布、温度分布、压力分布等。

根据模拟结果对气液两相流动进行分析，并与实验数据进行对比。

以上是一种简单的数值模拟方法，当涉及到更复杂的问题时，可能需要使用更高级的模型和方法来进行数值模拟。

此外，数值模拟的准确性还受到边界条件和网格划分的影响，需要进行适当的验证和调整。

液体蒸发的数值模拟
蔡玉强;李亚丛
【期刊名称】《华北理工大学学报：自然科学版》
【年(卷),期】2017(039)003
【摘要】基于fluent多相流混合模型,结合自定义函数（UDF）对水箱内部水的蒸发过程进行模拟分析,得出不同加热时刻的温度场、压力场、速度场和蒸汽体积分布,从而清楚地了解容器内的相变程度;分析当加热时间相同时,不同加热温度对沸腾情况的影响,便于更加经济有效地控制加热温度。

结果表明,该方法能够快速准确地
模拟液体的相变过程,并可反映相变程度,从而缩短相关容器设计、性能改进的周期。

【总页数】6页(P99-104)
【作者】蔡玉强;李亚丛
【作者单位】华北理工大学机械工程学院,河北唐山063210;华北理工大学机械工
程学院,河北唐山063210
【正文语种】中文
【中图分类】TQ026.4
【相关文献】
1.水平管降膜蒸发器管外液体流动数值模拟 [J], 宋小曼;杨丽;王伟洁
2.隧道内甲醇液体蒸发及蒸气扩散规律数值模拟分析 [J], 陈长坤;王玮玉;刘晅亚
3.液体蒸发的数值模拟 [J], 蔡玉强;李亚丛
4.水平滴形管降膜蒸发器管外液体流动数值模拟 [J], 王伟洁;杨丽;冀哲
5.隧道内障碍物对可燃液体蒸发特性影响的数值模拟研究 [J], 刘连冲
因版权原因，仅展示原文概要，查看原文内容请购买。

颗粒拟流体模型中的颗粒拟流体的静压力、颗粒相的切应力的各种表达方法、物理基础、优缺点1•颗粒拟流体模型的简介:颗粒拟流体模型又叫多流体(双流体)模型。

该模型将弥散颗粒相与连续流体相看作是连续介质，对颗粒相的处理方法与对连续介质相的处理方法类似，认为颗粒相是欧拉坐标系中与连续相流体相互渗透的一种“假想”流体，称为拟流体。

因此，这种模型又叫做颗粒拟流体模型。

该模型不仅考虑连续流体相与颗粒相之间存在的显著速度滑移和温度滑移，并且认为这种滑移与颗粒相的扩散是两种完全不同的作用，而且颗粒相的扩散是独立于流体相扩散之外的另一种运动特性。

该模型还引入了颗粒相粘性、扩散和导热系数这些与连续流体类似的物理性质。

颗粒拟流体模型的基本假设包括：(1)在流场中弥散颗粒相与连续流体相共存并且相互渗透，连续流体相和弥散相在计算区域中的任何一点共存，占据同一空间，但分别具有各自的速度、浓度、温度和体积分数等，而且在每个计算单元内只有一个值；若是将颗粒相按尺寸分组，则每个尺寸组的颗粒具有相同的速度和温度。

(2)在做体积平均后，每一尺寸组的颗粒相在空间中具有连续的速度分布、温度分布和容积分数的分布。

(3)每一个尺寸组的颗粒相除了与连续流体相具有质量、动量和能量间的相互作用之外，还具有自身的湍流脉动，并由此造成颗粒相自身的质量、动量和能量的湍流运输，因而具有其自身的湍流粘性、扩散和导热等湍流输运性质；对于稠密颗粒悬浮体，颗粒相之间的碰撞还会引起附加的颗粒粘性、扩散和热传导；因此，颗粒相具有类似于连续流体相的“拟”物理性质。

(4)弥散颗粒相可按初始尺寸分布分为不同的群组。

(5)连续流体相和颗粒相都在欧拉坐标系内描述，因此该模型也称为“双流体模型”，也叫“欧拉-欧拉模型”。

颗粒拟流体模型的基本方程组包括下述一些方程。

连续流体相的连续方程dp d / 、v1示+鬲Mr $—2/皿弥散颗粒相的连续方程连续流体相动量方程弥散颗粒相动量方程d d , 、无(n k v ki ) + 苏(n k v kj v ki J = n k (Vi 一 v ki )/T rk + n k gi连续流体相能量方程专（PCpT ）+右 3jCpT ）dxj 匕丁幽弥散颗粒相能量方程 —(n k c k T k ) + ^(n k v kj c k T k ) = n k (Q h 一 Qk - Qrk)/m kd流体的组分方程 爲（P 。

航空发动机滚动轴承内油气两相流动与温度场的数值模拟徐让书;李骏;牛玲;常柱宇【摘要】基于VOF方法和滑移网格技术，研究喷射润滑下主轴承内部压力、温度场及油气体积比随转子转速及滑油流量的变化规律。

给出滑油在轴承内的相界面分布，描述轴承内油膜分布及运动。

结果表明：轴承内压力和油气两相在轴承内分布不均匀；油液体积分数随转速的升高而降低，随供油量的增加而增加，呈非线性关系；轴承内最高工作温度随主轴转速的增加而增加；轴承润滑效果与供油量有直接的关系，不同的转速下对应不同的最佳供油量。

%Based on the volume of fluid(VOF)method and sliding mesh technique,the effect of rotor speed and the oil inlet flow rate on the pressure,temperature field and oil volume fraction in the rolling bearing with spay lubrication were in-vestigated.The film distribution and movement in oil and air interface were showed.The results show that the pressure and the oil-air two-phase flow distribution is not homogenized in the bearing.The oil volume fractionis decreased with the in-crease of revolution speed,and is increased with speeding up the lubrication flow rate nonlinearly.The highest temperature of the bearing is increased with the increase of revolution speed.The effect of the bearing lubrication is directly related with the lubrication flow rate,and an optimum oil feeding rate exists at different rotor speed.【期刊名称】《润滑与密封》【年(卷),期】2014(000)003【总页数】5页(P91-95)【关键词】航空发动机;滚动轴承;气液两相流;温度分布;VOF模型【作者】徐让书;李骏;牛玲;常柱宇【作者单位】沈阳航空航天大学航空航天学部辽宁沈阳 110136;沈阳航空航天大学航空航天学部辽宁沈阳 110136;沈阳航空航天大学航空航天学部辽宁沈阳110136;沈阳航空航天大学航空航天学部辽宁沈阳 110136【正文语种】中文【中图分类】V228.2航空发动机主要采用滚动轴承，高DN值下轴承将产生大量的摩擦热量并引起较高的温升[1]。

2006年第5期 总第153期低 温 工 程CRYOGENICSNo.5 2006Sum No.153垂直上升管内气液两相流动特性的数值模拟何鸿辉 刘国青 刘波涛（北京卫星环境工程研究所 北京 100029）摘 要：某空间模拟器的液氮制冷系统采用了结构简单、可靠性高、最节能的重力自循环系统。

基于气液两相流理论，建立了垂直上升管内液氮气液两相流动的一维均相流模型；数值模拟不同入口压力下的管内压力、质量流量、管壁温度和质量含气率等参量，研究其变化规律。

计算结果表明，建立的数学模型能够详细描述重力自循环系统垂直上升管内液氮流动的流体动力特性，为空间模拟器液氮制冷系统的研制提供了重要的基础数据和理论参考依据。

关键词：气液两相流 数值模拟 液氮 重力自循环 空间模拟器中图分类号：TB663，V573 文献标识码：A 文章编号：1000-6516（2006）05-0022-05收稿日期：2006-06-12；修订日期：2006-07-13Numerical simulation of gas-liguid two phase flowin the vertical tubeHe Honghui Liu Guoging Liu Botao（Beijing Institute of Spacecraft Environment Engineering ，Beijing 100029，China ）Abstract ：The naturai circuiation system with simpie construction ，high reiiabiiity and efficient ener-gy saving was appiied in the iiguid nitrogen refrigeration system of space simuiator.One-dimensionai hom-ogeneous modei based on the gas-iiguid two-phase fiow theory was estabiished.The pressure ，the mass fiux ，the waii temperature and the porosity parameters were simuiated numericaiiy.The reguiarity varied with the time-space was studied.The resuits showed the mathematicai modei couid describe the hydrody-namic characteristics of iiguid nitrogen in the verticai tube and provided the important originai dada and theoretic reference for the deveiogoment of the iiguid nitrogen refrigeration system of space simuiator.Key words ：gas-iiguid two-phase fiow ；numericai simuiation ；iiguid nitrogen ；naturai circuiation ；space simuiator1 引 言某真空热试验设备，主要由热真空系统（包括真空系统、液氮系统、控制系统等）和太阳模拟器系统组成，可模拟太空的真空、冷黑及太阳辐照环境。

2012年4月Water Conservancy Science and Technology and Economy Apr.，2012［收稿日期］2011－12－19［基金项目］国家自然科学基金(51076041)［作者简介］冯俊(1988－)，男，湖南长沙人，硕士研究生，从事流体机械工程的研究;郑源(1964－)，男，山东日照人，教授，博士生导师，从事流体机械工程及水利水电工程的研究.基于气液两相流模型的水轮机数值模拟研究进展冯俊，郑源(河海大学能源与电气学院，南京211100)［摘要］随着水电站机组容量不断增大和计算机技术的飞速发展，传统的数值模拟手段已不能满足人们对深入分析水轮机性能的要求，因此更贴近实际、结果更准确的基于气液两相流模型的数值模拟越来越受到人们的关注。

介绍了气液两相流动的主要计算模型，概述近年来国内外水轮机气液两相流数值模拟进展，重点阐述在其应用于水轮机空化性能、尾水管压力脉动等方面的研究成果和应用前景，并对需要进一步深入研究的问题及研究中应遵循的一些原则作出初步探讨。

［关键词］气液两相流;水轮机;数值模拟;空化;压力脉动［中图分类号］TV74［文献标识码］A［文章编号］1006－7175(2012)04－0004－04Study of Water Turbine's Gas －Liquid Two －PhaseFlow Based Numerical SimulationFENG Jun ，ZHENG Yuan(The College of Energy and Electrical Engineering HoHai University ，Nanjing 211100，China )Abstract :With the development of water turbines'capacity and computer technology ，people start to focus on the need for accurate numerical simulation to analyse turbines'internal flow.The Gas －Liq-uid two －phase flow model which is closer to actual is now drawing more and more attention from people.This article will be followed by a brief instruction of main Gas －Liquid two －phase flow mod-els and an overview of latest international progress of water turbine's Gas －Liquid flow based numeri-cal simulation.Details on applications in water turbines'cavitation and draft tube's pressure pulsation are described in later sections.Finally ，it shows a discussion on problems and methods to solve them.Key words :gas －liquid two －phase flow ;water turbine ;numerical simulation ;cavitation ;pressure pulsation效率、空化性能和稳定性是影响水力机械机组性能的3项指标。