本案例利用Fluent中的Evaporation-Condensation模拟模拟闪蒸过程。

1 问题描述

计算模型如下图所示。高温高压的液态水从入口流入计算区域，入口温度400 K，压力0.3 MPa。容器内初始压力1 atm，初始温度293 K，液相进入容器后由于其温度高于饱和温度而发生闪蒸。

2 Fluent数值

• 以2D、Double Precision 方式启动Fluent
• 利用菜单File → Read → Mesh… 读取网格文件flashing.msh

2.1 General设置

• 进入General 面板，选择选项Transient及Axisymmetric，并指定重力加速度为X轴-9.81 m/s

• 点击Scale… 按钮打开网格缩放对话框，如下图所示，将计算域尺寸缩放到0.642*0.075 m

2.2 Models设置

• 打开多相流设置对话框
• 选择Mixture模型
• 激活选项Implicit Body Force
• 其他参数保持默认设置

• 激活能量方程

• 选择采用Realizable k-epsilon湍流模型

2.3 Materials设置

• 从材料数据库中添加液态水与水蒸气

• 修改液态水物性参数，如下图所示

• 修改水蒸气物性参数，如下图所示

注：利用REFPROP计算得到1 atm条件下的水的汽化热为2256.5 kJ/kg，将该值与分子量相乘换算为标准状态焓。

2.4 设置相

• 设置液态水为主相

• 设置水蒸气为次相

• 选择相间传质模型为evaporation-condensation

• 设置Lee模型参数，如下图所示

注：From Phase Frequency与To Phase Frequency分别为蒸发频率因子与冷凝频率因子，这两个参数通常用于调整蒸发量与冷凝量与试验值进行对标用的。若没有试验值，对于同时存在蒸发与冷凝的问题（如热管等），通常将冷凝因子调大，冷凝因子与蒸发因子的比率大约为液相与气相密度的比率。本案例主要为蒸发，为了更容易观察到蒸发现象，可以将蒸发频率设置一个较大的值。千万注意，没有经过试验数据校准的Lee模型，是不可用用于定量研究的。

2.5 设置边界条件

1、velocity-inlet_inlet边界

• 将边界velocity-inlet_inlet 设置为Pressure-Inlet类型
• 指定总压为0.3e6 Pa

• 设置入口温度400 K

• 设置vapor相的体积分数为0

2、output边界

• 将output边界类型指定为pressure-outlet
• 设置出口静压为0 Pa

• 指定出口回流温度为293 K

• 指定出口汽相回流体积分数为0

2.6 Methods设置

• 如下图所示指定计算方法

2.7 Controls设置

• 如下图所示设置求解控制参数

2.8 Monitor设置

• 将连续方程的收敛残差设置为1e-8

2.9 初始化

• 如下图所示设置初始值

2.10 进行计算

• 设置时间步数10000，时间步长为0.0001 s进行计算

3 计算结果

计算时间很长，这里只计算了0.03 s。

• 速度变化
• 汽相分布

免责声明：本文系网络转载或改编，未找到原创作者，版权归原作者所有。如涉及版权，请联系删

武汉格发信息技术有限公司，格发许可优化管理系统可以帮你评估贵公司软件许可的真实需求，再低成本合规性管理软件许可,帮助贵司提高软件投资回报率，为软件采购、使用提供科学决策依据。支持的软件有: CAD,CAE,PDM,PLM,Catia,Ugnx, AutoCAD, Pro/E, Solidworks ,Hyperworks, Protel,CAXA,OpenWorks LandMark,MATLAB,Enovia,Winchill,TeamCenter,MathCAD,Ansys, Abaqus,ls-dyna, Fluent, MSC,Bentley,License,UG,ug,catia,Dassault Systèmes,AutoDesk,Altair,autocad,PTC,SolidWorks,Ansys,Siemens PLM Software,Paradigm,Mathworks,Borland,AVEVA,ESRI,hP,Solibri,Progman,Leica,Cadence,IBM,SIMULIA,Citrix,Sybase,Schlumberger,MSC Products...