COMSOL三相流：界面粗糙算不准，初始步太大直接不一致收敛失败

石油开采、化工输送、多相分离中，油、水、气三相流极常见。两种不溶液体加一种气体，非线性强、相界面演化复杂，网格不精密界面就粗糙。这篇以二维三通注入模型为例，把 COMSOL 三相流相场设置的 5 个关键点讲清楚。

要点一：选三相流相场，研究用瞬态

模型向导选二维 → 流体流动 → 多相流 → 三相流相场（Three-Phase Flow, Phase Field），研究选瞬态。相界面靠相场变量捕捉，需较密网格和一定算力。

要点二：几何分区，用分割域做出两个入口

主矩形宽 2m、高 3m（空气域）；再建宽 1m、高 0.5m 的斜块，基准 X=-0.7、Y=2.8、旋转 ±45°，并用变换→镜像生成另一侧（镜像前保留输入对象）。然后用布尔与分割→分割域，按边切分入口区域并删除多余实体，最后形成联合体。

图 1　相场设置与三组表面张力配置

要点三：迁移率先用默认，材料配油、水、空气

旧版本需自定义迁移率变量（含体积分数 φA、φB、φC 与 Jac 算子）以改善收敛；新版本（6.x）对相场求解器已改进，迁移率保持默认即可，自定义反而可能报错。材料从库中添加空气、水，并选一种油（如玉米油）；一切以物理场里的材料分配为准。

流动勾选包含重力，湍流模型选 Algebraic（代数/ABL 类）。

要点四：表面张力按三对组合，相A油/B水/C空气

三种材料有三对表面张力：A-B 为液液（油-水），B-C、A-C 为液气（水-空气）。相场中分配：A=油、B=水、C=空气，并勾选「在动量方程中包含表面张力」——不勾时方程中对应项被去除。界面厚度控制取默认。

边界：左右入口速度默认 1.2m/s（收敛难时调小），加压力点约束。初始值分三个：右/左入口分别赋 B、A 相，上部空气域赋 C 相，与入口边界一一对应。

图 2　三相材料分配与初始相设置

要点五：初始不一致就降速度、用小步长起步

网格用「更细」。常见报错是「找不到初始一致值」「非整数次幂」。两个有效调整：①把两入口速度都减小（如降到 0.001m/s）；②变步长——初始时间步设 1×10⁻⁵，先算到 0.001s，再放大步长算到 0.5s，瞬态求解器勾「初始时间步」并与之一致。分离求解器下限若遇除零，可用 EPS 极小值代替零。

后处理：看速度、压力、各相体积分数与三相混合分布；颜色表切「离散」轮廓更鲜明；面箭头可加过滤器（如 B 相体积分数 >0.5）只看单相；还可做二维截线、二维拉伸三维体与等值面（水平 0.5）。

图 3　三相体积分数分布结果

小结

①选三相流相场+瞬态；②分割域做出两斜入口；③迁移率新版本用默认，材料配油/水/空气、含重力与代数湍流；④三对表面张力配齐，A油/B水/C空气并勾动量方程表面张力，三初始值对应三入口；⑤初始不一致就降入口速度、以 1×10⁻⁵ 小步长起步，除零用 EPS。

【系统学习 COMSOL 多物理场】

从单相到多相流、再到多物理场耦合，COMSOL 多物理场耦合仿真体系课程配合工程师能力认证与项目实训，带你独立完成复杂流场与耦合分析。想从流体方向拓展多物理场的工程师可重点了解。

【配套资料领取】

扫码添加小助理，私信回复关键词【COMSOL三相流】，领取本篇配套资料包：①三相流相场核心知识总结 ②表面张力与迁移率关键参数速查表 ③油水气三相流实战案例（分步） ④初始不一致与收敛易错清单。