Fluent鼓泡塔（欧拉欧拉+PBM）

[ZY_Feng]

圆柱鼓泡塔，几何模型和如图所示，Fluent三维计算

主要参考文献是 郭晓峰的博士论文（https://kns.cnki.net/kcms/detail/detail.aspx?dbcode=CDFD&dbname=CDFDLAST2018&filename=1018225706.nh&uniplatform=NZKPT&v=dxCmfTYfWep4o09HDEiFLBfLAiNXVIlckP4zPHTHnqC7e9nKS3hTYfdNbMs4EzKk）
网格结构化网格，离散格式空间上动量和体积分数采用QUICK格式，其他为二阶迎风，时间上一阶
边界条件如图所示，入口为气体边界0.01~0.03m/s
多相流模型为欧拉，气泡尺寸采用平均单一平均值或是PBM模型，相间相互作用考虑曳力、升力、壁面润滑力，表面张力
但是现在只考虑表面张力和曳力可以计算，但是径向含气率分布有问题，但添加升力和壁面润滑力就发散了，调整亚松弛因子和时间步长都尝试过了吗，各位老师怎么才能添加这些力？是我设置有什么问题吗？

[李东岳]

升力和壁面润华力确实会引起收敛性的问题。不过类似这种算例，如果不存在那种径向的峰值分布的话，只需要考虑液力就行。

针对你这个算例，只能多调试，没有太好的建议。

[ZY_Feng]

@李东岳 @谢谢东岳老师的回复。我这个模拟是比较关注径向含气率分布的（和实验数据对比），所以升力和壁面润滑力可能有较大影响。

[李东岳]

@zy_feng 壁面力这面主要是升力，有可能会过大导致气泡跑出去，OpenFOAM里面是增加了一个limiter，Fluent那面不知道怎么处理的。不过确实会存在这个问题。壁面润华力、湍流分散力问题都不大。你看下去掉升力会不会发散。

理论部分可以看这个文章的方程14 AIChE.pdf

[随风飘姚]

我记得有些鼓泡塔的文献提过升力会影响计算收敛，尤其是气泡聚集时（云图出现大面积红色）残差会变得很大。这里建议先检查网格和时间步长是否匹配，适当减小时间步长，然后在UDF中添加
if (vof > 0.99)
cl = 0;
祝收敛！