FLUENT管内沸腾及蒸汽继续加热案例

    之前我们用自带的Evaporation&Condensation模型模拟一下水的沸腾问题，并提出一些可能在其他案例没有的观点供参考（可以点击文末“阅读原文”查看）。今天我们继续做一个管内水沸腾制蒸汽的案例。

    本次案例模型来自于FLUENT 2019R3官方验证案例VMFL067，只是在原来的基础上稍加修改。原案例是欠热水的临界热流密度（CHF）模拟，本案例将入口水的温度设置为饱和温度，这样壁面的热功率全部用于制蒸汽，因为压力不变，水的饱和温度不变，饱和状态的水吸热即沸腾蒸发生成水蒸气。

    原案例重要参数：管径10mm，管长7m，入口水的质量流率1495kg/m2s（单相），饱和压力7.0176MPa，饱和温度286℃，入口水温282℃（欠热），壁面热流密度797000W/m2。本案例部分参数修改为：入口水温286℃，壁面热流密度1060000 W/m2（总热功率233106.17W），其余参数不变。

    本案例需进行如下设置，否则可能无法得到案例的结果。

    首先，选择稳态求解，开启重力选项。

    选择欧拉多相流模型。

    选择水为主相，水蒸气为次相，水蒸气的气泡直径为boiling-dia。

    开启沸腾模型，饱和温度559K。

    设置水和水蒸气的物性参数，这里重点需要注意：水蒸气和水的标准状态焓值之差等于蒸发潜热。我们可以通过IF97数据库查对应温度下的饱和水以及饱和蒸汽焓值，水的标准状态焓为2.2849e+07，水蒸气的标准状态焓4.9945e+07j/kmol。IF97数据库查到水和水蒸气的焓值（286℃）分别为1268.342637kj/kg和2772.339798kj/kg，FLUENT标准状态焓的单位为j/kmol，之前说过我们只需要将查到的相变潜热乘以水的摩尔质量18.0152kg/kmol即可。

    沸腾模型只能在湍流中可用，我们设置流态为湍流，选择RNG k-e湍流模型。

    求解器进行如下设置。

    亚松驰因子进行如下设置，亚松驰因子的设置对计算结果的收敛影响极大，需要多次进行尝试。

    我们设置两个监视器，一个监测入口的水流量，一个监测出口蒸汽流量。

    从入口进行初始化。

    入口的水流量变化趋势如下图。

    出口的蒸汽流量变化趋势如下图。

    我们读取出入口的总流量如下图。

    从以上可以看出，饱和水被加热后沸腾蒸发成水蒸气，达到稳态时全部蒸发为水蒸气了，这一点也可以从下图所示管内蒸汽体积分数看出。

    下面我们进行能量守恒等核算。

    出口水蒸气的流量为0.11878518kg/s（入口水的流量0.1174kg/s，不平衡率1.2%），根据IF97数据，286℃的蒸发潜热为1503997.161j/kg，因此沸腾蒸发吸收的热量为178652.5735W，壁面传热量为233106.17W，相变消耗的热量占总热量的76.64%，其余的热量应该是继续加热蒸汽了。我们读取出口的蒸汽温度为364.8℃，此时的蒸汽处于过热状态了。我们从IF97读取7.0176MPa，364.8℃过热蒸汽的焓值3060.65kj/kg，以此值减去水的焓值之差乘以流量积为介质吸收的总热量为212899.8985W，与壁面传热量不平衡率为8.7%。

    以上是用IF97的数据进行核算，接下来我们用FLUENT计算结果数据进行核算。进出口的能量之差为236688.99W，与壁面传热量不平衡率1.5%。

    出口蒸汽的焓值为4667567.1j/kg，如前所述286℃蒸汽的焓值为2772359.56j/kg，焓差为1895207.544j/kg，因此蒸汽总吸收热量为225122.5693W（焓差乘以蒸汽流量），与壁面传热量不平衡率为3.4%，与IF97数据计算不平衡率为5.7%。

    根据以上核算，本次案例结果还是相对较好的，但是沸腾传热模拟相当难收敛，稍微设置不合适就发散，建议大家在算完之后务必进行守恒核算。