Fluent 多相流仿真：竖直管内气液两相弹状流特性分析

本项目对竖直上升管内的气液两相弹状流（slug flow）进行 CFD 仿真，目标是预测 Taylor 气泡结构、弹单元长度与压降梯度，为石油、化工及核能领域的两相流设备设计提供依据。管道内径 40 mm，介质为空气-水，在 20 °C、常压下运行。仿真在 ANSYS Fluent 中完成，采用 VOF（Volume of Fluid）模型捕捉气液界面，湍流模型选用 SST k-ω。

几何模型为 2D 轴对称竖直管，管长 2 m，以保证充分发展。边界条件：底部为速度进口，液相表观速度 0.8 m/s，气相表观速度 1.6 m/s；顶部为压力出口。表面张力取 0.072 N/m，壁面接触角 60°。网格在管壁与气液界面处自适应加密，时间步长 0.001 s，总计算时间 1.6 s。分析输出：截面含气率时序、Taylor 气泡频率、弹单元长度与压降梯度。

案例分析

在表观液速 0.8 m/s、表观气速 1.6 m/s 的设计点，仿真得到 Taylor 气泡弹状频率 1.9 Hz，单个弹单元长度 285 mm，管道截面平均含气率 34%。压降梯度为 3.96 kPa/m，主要由液相重力与界面摩擦贡献。滑移比（气相实际速度与液相实际速度之比）为 1.42，说明气相速度高于液相，符合弹状流物理特征。

与 Lockhart-Martinelli 经验关联式对比，在设计点 L-M 关联式预测压降梯度比 VOF 仿真高约 18%，原因是经验关联未充分考虑弹状流的周期性加速度压降。流场动画显示，Taylor 气泡尾部存在明显液膜回流，尾部夹带的小气泡数量随气速增加而增多。若气速提高到 2.5 m/s，弹状流会逐渐向搅混流转变，频率升至 3.2 Hz，弹单元长度缩短至 180 mm。

总结教训

本项目的第一个教训是 VOF 模型对网格与时间步长敏感。若轴向网格尺寸大于 5 mm，Taylor 气泡界面会被抹平，频率与长度都会失真；时间步长超过 0.002 s 则会出现数值弥散。第二个教训是入口边界条件必须能够形成稳定的弹状流，采用均匀速度进口需要较长的入口发展段，否则入口效应会污染结果。第三个教训是表面张力系数必须准确。20 °C 空气-水表面张力 0.072 N/m，若取 0.06 N/m，气泡形状与频率偏差明显。

项目经验

1. 竖直管弹状流仿真建议采用 Fluent VOF 模型。
2. 管道内径 40 mm，空气-水，液相表观速度 0.8 m/s，气相 1.6 m/s。
3. Taylor 气泡弹状频率 1.9 Hz，弹单元长度 285 mm，截面含气率 34%。
4. 设计点压降梯度 3.96 kPa/m，滑移比 1.42。
5. Lockhart-Martinelli 关联式在设计点高估压降约 18%。
6. 轴向网格建议小于 5 mm，时间步长小于 0.002 s。
7. 入口发展段需足够长（≥40 倍管径），避免入口效应。
8. 表面张力取 0.072 N/m，接触角 60°，对气泡形态影响大。
9. Taylor 气泡尾部存在液膜回流与小气泡夹带。
10. 气速升高到 2.5 m/s 时，弹状流向搅混流过渡。
11. 压降梯度由重力、摩擦与周期性加速度压降共同组成。
12. 本项目方法可用于油气输送、核反应堆冷却剂通道等两相流分析。

补充经验

1. 水平管内弹状流受重力影响，气泡偏向管顶部，与竖直管相比频率降低、长度增加。
2. 环状流出现在高气速区，液膜沿壁面稳定流动，压降预测需采用液膜厚度经验关联。
3. 油气分离器设计中，入口挡板与重力沉降段长度应根据弹单元长度 285 mm 的整数倍设计。

复现参数速查（案例图硬参数）

更多 Fluent 多相流仿真经验，可参考本站的 Fluent 仿真案例栏目。