151、Fluent VOF模型显式和隐式格式的对比

Fluent中的VOF模型可以用来计算气-液、液-液这类不互溶流体。

VOF模型存在两种格式，分别是显式Explicit和隐式Implicit格式，这两种方式决定了VOF界面的计算精度、稳定性以及适用场景。

这篇文章就来详细对比VOF模型的显式和隐式格式，我们先从概念上进行对比，然后在Fluent中实际操作计算，看看两者的计算区别

VOF模型的核心思想很简单：引入一个标量场，体积分数αq表示第q相在某个网格单元内所占的体积比例。

αq=1：单元内全是q相αq=0：单元内完全没有q相0<αq<1：单元内存在两相交界面αq

体积分数输运方程：

这个方程形式上很简单，但真正困难的地方在于：如何保证界面在输运过程中保持尖锐，不发生数值扩散。这正是显式与隐式格式最大的区别所在

2.1 显式格式

在每个时间步开始时，利用上一时间步已知的α(n)和速度u，直接计算出新的 α(n+1)。离散形式可以写成：

特点：

• 按照界面几何形状搬运流体，更加清晰

• 强依赖时间步长，也就是库朗数要小于1.0

• Fluent中通常配合Geo-Reconstruct（几何重构）

2.2 隐式格式

和显式方法不同，隐式方法使用未来时刻的未知量来计算体积分数：

特点：

• 对库朗数没有限制，类似时间隐式格式，可参考文章：150、Fluent瞬态计算的库朗数必须小于1吗？库朗数和时间积分格式的关系？

• 通常采用代数格式（如HRIC）

• 用数值扩散方式传播体积分数，界面不够清晰

3.1 优势

显式格式之所以能实现极高的界面精度，核心原因在于它可以配合Geo-Reconstruct格式。

Geo-Reconstruct 做的事情非常直观：在每个时间步开始时，根据每个混合单元内的α值，实际地计算出界面的位置和方向，然后利用这个几何信息精确计算通过网格面的流量。

这样做的好处是：

几何重构保留了界面的物理形状，而不是用代数函数去插值，所以界面更加清晰。

数值扩散被降至最低。因为流量计算基于真实的几何界面，而不是基于梯度或插值，所以人为的“抹平”效应几乎为零。

类比一下：显式+Geo-Reconstruct 就像用钢笔描边，线条锐利；隐式格式更像用毛笔晕染，边界模糊。

3.2 局限

显式格式计算出的界面不基于梯度或插值，因此界面比较清晰。但这也使得数值收敛性较差。存在以下局限性：

• 显式格式要求 CFL ≤ 1

• 每个时间步内还要进行子循环，总计算量往往比隐式高出数倍甚至一个数量级。

• 仅支持瞬态计算：显式格式无法用于稳态问题。

• 时间精度受限：只能配合一阶时间格式

4.1 优势

隐式格式虽然界面不够清晰，但收敛性更好。

• 库朗数数可以设为 2、5 甚至更高，Fluent 依然能稳定求解。

• 由于体积分数与压力-速度场紧密耦合，收敛行为非常好，不易发散。

• 隐式格式支持稳态计算

4.2 局限

隐式VOF通常配合HRIC或CICSAM等代数压缩格式。这些格式本质上是通过限制器来压缩界面，但无法完全消除数值扩散，界面总是比较模糊。

对于液滴变形、气泡聚并等表面张力主导的流动，界面模糊会直接导致受力计算失真。

我们使用完全相同的参数，仅仅更改显式和隐式的格式。两者时间步长都是0.001s，体积分数空间离散格式：显式格式选择几何重构方法；隐式格式选择Modified HRIC格式。其他的保持二阶格式，如下图

结果如下：明显显式格式的界面更加清晰，而隐式格式界面模糊。一般如果研究涉及液滴动力学和表面张力效应这类问题，优先使用显式+Geo-Reconstruct