1 模型描述
本例的模型较为简单，如图1所示。计算域高1m，宽0.2m。顶部边界为压力出口，底部有一高温壁面hotwall，温度570K，其他壁面wall为绝热边界。计算域内初始充满0.9m深的水。划分网格如图2所示。
2 导入网格
打开fluent，导入上步生成的网格模型。Scale检查网格尺寸。如图3所示。
图3 scale计算域
确保计算域尺寸是我们所需要的。本例中x方向尺寸0~0.2m，y方向0~1m。
3 设置求解器
选择压力基(pressure-based)求解器，同时选择瞬态模拟。
由于水沸腾时水蒸气会在浮力作用下向出口运动，因此考虑重力。设置重力加速度为重力加速度为y方向，大小-9.81m/s2。
4 设置计算模型
添加多相流模型为mixture模型，勾选slip velocity及implicit body force，设置欧拉相数量为2。如
激活能量方程
此例为层流流动，不激活湍流模型。
添加材料water-vapor及water-liquid。修改材料属性。
Water-liquid属性：
Density：1000kg/m3
Viscosity：0.0009kg/m-s
Standard state enthalpy：0 j/kgmol
Reference temperature：298k
其他采用默认属性。
Water-vapor属性：
Cp：2014j/kg-k
Standard state enthalpy：2.99e7j/kgmol
Reference temperature：298.15k
其他属性采用默认设置。
6 相设置
设置液态水为主相，水蒸气为第二相。设置气泡直径0.0002m。如图7，图8所示。
图7 设置主相
图8 设置第二相
定义相间作用，即定义蒸发/冷凝模型。点击interaction按钮进入如图9所示对话框。选择from phase为liquid，选择to phase为vapor，即为由液态转化为气态，选择模型为evaproation-condensation。此时弹出模型定义对话框，如图10所示。这里采用默认设置即可。
7 设置边界条件
本例的边界条件较为简单。
Outlet：采用压力出口，设置出口蒸气含量100%，温度372K
Walls：设置为绝热边界。设置heat flux为0
Hot wall：设置temperature为570K
设置operating conditions如图11所示。
图11 operation condition
设置参考密度为气相密度0.5542。
8 求解方法
设置pressure为body force weighted，其他momentum、volume fraction、energy全采用quick算法，有利于提高计算精度。
在solution controls面板中设置亚松弛因子。设置pressurewie为0.5，momentum为0.2,volume fration为0.2。
9 初始化
初始化时设置temperature为372K。
需要patch两个区域：
1、与hotwall相邻网格节点。由于hotwall温度高达570K，超出水的沸点373.15K，因此需要adapt出hotwall边界相邻节点区域，设置其温度为373.15K。
2、Patch出初始水位。由于初始状态下计算域中有深0.9m的水，因此需要通过patch将其标记出来。
【Adapt】>【boundary】，弹出如图12所示对话框，在boundary zones中选择hotwall，设置number of cells为1，点击mark即可对此区域节点进行标记，下一步进行温度patch。
图12 标记临壁面节点
点击patch…按钮，弹出图13所示对话框。
图13 温度patch
下一步patch水位。
点击菜单【Adapt】>【region…】进入区域标记对话框。进行如图14所示设置。
点击mark进行标记。
图14 区域标记
进入patch对话框，进行如图15所示设置。点击patch完成设置。
图15 水位patch
此时可以查看云图，以检查patch是否正确。
10 其他设置
OK，该定义的都定义完了，剩下的就是一些监测量、动画什么的设置，这些都很基础，这里不再赘述。
11 设置计算
设置时间步长0.01s，时间步数1000步，迭代计算。