6.3 相变与特殊多相流

6.3 相变与特殊多相流

多相地图的最后一块：当两相之间发生质量交换——水闪蒸成汽、蒸汽凝回水、低压区液体拉出空泡——界面不再是几何边界而是反应前沿。本节梳理蒸发、冷凝、空化三类相变在 Fluent 里的建模入口，顺带清点多孔介质与活塞流（后文运动章节的近亲）这两类"特殊形态"，把多相模型关的清单收尾。

相变建模的账本结构

相变的核心是界面质量源项：单位体积单位时间从一相转移到另一相的质量。它同时挂进三个账本——质量账本（一相减、一相加，总量守恒）、能量账本（潜热吸收或释放，量级巨大：水的汽化潜热约为同质量升温一百度的五倍多）、动量账本（界面动量交换随之调整）。Fluent 的相变模型本质都是在回答同一个问题：界面处的传质速率由什么驱动、怎么算。

Lee 模型（蒸发-冷凝机理模型）。给界面传质速率安一个与偏离饱和程度成正比的系数式：过热度驱动蒸发、过冷度驱动冷凝。两个传质强度系数（蒸发系数与冷凝系数）是模型的心脏，量级通常在零点一到一百每秒之间，需要按经验或敏感性分析标定。它挂在 VOF 或 Eulerian 框架下，是蒸发器、冷凝器、热管仿真的常用配置。

空化模型族。空化是局部压力跌破饱和蒸汽压时液体的暴力汽化，出现在泵进口、螺旋桨叶面、喷嘴内部。Fluent 提供从简化到精细的一族：Rayleigh-Plesset 模型从单气泡动力学导出传质率，Schnerr-Sauer 与 Zwart-Gerber-Belamri 是它的两个工程化参数化版本，差别在气泡数密度与成核假设。空化模型还有一个隐藏开关——不凝结气体质量分数，它控制完全空化区的密度下限，防止方程组在纯汽区奇异。

图 6-3：相变模型的驱动机制与挂载框架

案例：离心泵进口空化的性能预测

背景。某离心泵在低吸入压力工况出现性能跌落与噪声，厂方需要评估空化初生裕度并定位空化区。试验台架可测扬程曲线但看不到叶轮内部。

操作。单叶轮流道周期模型（第八章会讲旋转域的 MRF 与滑移网格处理，此处用 MRF 初算定常性能），湍流 SST，两相用 Mixture 框架挂 Schnerr-Sauer 空化模型，饱和压力按进口温度查物性表，气泡数密度按经验区间做敏感性扫描。先算无空化基准扬程对标试验，再逐级下调进口压力触发空化。

结果。空化区首先出现在叶片吸力面前缘背面，与高速摄影的白色区域位置一致；扬程跌落拐点的进口压力预测偏差在百分之五以内。扫描发现气泡数密度取值影响空化区面积三成以上，取中值并记录敏感性区间进报告。

解读。空化仿真的可信度建立在与试验的"两处对表"：空化位置对（定性验证）加扬程拐点对（定量验证）。只对扬程不对位置，可能蒙对结果错认机理。案例还示范了多相仿真的不确定度声明写法：敏感参数给出扫描区间与中值，结论限定在区间内成立。

变式。喷嘴内部空化与射流雾化的耦合需要 VOF 加空化的高保真组合，网格在喷孔内大幅加密；冷凝类任务（蒸汽轮机湿蒸汽）用湿蒸汽模型或组分法，物性与成核模型绑定选型。

相变配置问答与异常对照

Lee 模型的两个传质系数取多少。它们是速率常数不是物理常数：文献里常用量级在零点一到一百每秒之间，不同任务差几个数量级都正常——正确做法是中值起步加敏感性扫描，把结论对系数区间的依赖写进报告。

空化算例怎么会算出全流场空化。多半是进口压力边界给得太低或参考压力错位，空化模型诚实地把低压区全拉了汽。先核对压力参考系与饱和压力的量级，再怀疑模型。

不凝结气体质量分数的作用。它是空化区的密度下限保护：取零会让纯汽区方程奇异，取过大又压制空化发展。工程常用小分数起步，随敏感性分析微调。

相变算例异常对照：

关于多孔介质的补充说明（很多任务在此与多相流狭路相逢）：多孔介质把固体骨架的阻力效果压缩成动量方程里的黏性损失与惯性损失两项，参数来自压降-流速实验曲线的拟合。它与多相模型的组合有两个高频坑：其一，浸泡相含率设置错误会让"孔隙里的流体量"与实际不符，直接污染停留时间类结论；其二，多孔区网格粗化与相界面捕捉的冲突——界面穿过多孔区时分辨率骤变，界面表现会失真，关键相变过程避开布置在多孔区内是布局层面的解法。多孔与相变的联合任务，先在局部小模型上验证两者的兼容性再上整机。

把相变任务的验收思路做个收拢：相变仿真比普通多相多出两份必须签字的账。第一份是强度账：相变速率的合理性——蒸发率换算成界面后退速度或产汽量级，与热流密度、潜热的收支自洽（单位面积蒸发掉的水需要的功率与当地热流对得上）。第二份是分布账：相变的落点合理——蒸发集中在过热最强的位置、冷凝集中在过冷最强的位置，分布与温度场互相印证。两份账都过，相变模型才算真正"挂对"，而不是只挂了个名字。做相变任务久了你会发现，绝大多数诡异结果都源于这两份账里的一份没对上，而查账的成本远低于盲调参数的成本。

补一个跨章联动提醒：相变任务的瞬态化往往不可避免。稳态假设在强相变场景（闪蒸、深度空化）里很难自洽，因为相变界面本身就是非定常结构——此时的处理思路参考第八章：界面形态用瞬态统计量描述（时均界面位置、脉动幅值），而不是期待一张"稳态界面图"。瞬态相变的时间步受界面推进速度与潜热释放速率双重约束，比普通瞬态更严。把这些提前写进计划，比算到一半发现稳态解不下来再改要体面得多——相变仿真的项目排期，应该按"必然瞬态"来打预算，稳态跑通算意外之喜。

本节验收点

• 相变模型选择与驱动机制匹配：温差驱动用 Lee，压力驱动用空化族。
• 传质强度系数或气泡数密度做过敏感性扫描，区间写入报告。
• 空化算例完成位置与性能拐点两处对表验证。
• 三本账（质量、能量、动量）平衡核对通过，潜热源项已激活。