VOF法 — CAE用语解说

体积流体法（Volume of Fluid）是一种用来追踪多种互不混溶流体（例如水和空气）界面位置和形状的手法。在每个单元中定义"体积率"这样一个标量量。如果水占据单元体积的60%，体积率就是0.6；空气占40%，就是0.4。这个值在0到1之间连续变化，从而以模糊的形式表现界面。

很好的观察。确实界面在数个单元中会变模糊，但这样做计算非常稳定。为了保持界面清晰，在移流计算后进行"界面重构"处理。例如，PLIC（分段线性界面计算）法通过体积率分布推算界面的法向方向和位置，在每个单元内重构直线界面。这样界面的厚度实际上被保持在约1个单元的程度，相对清晰。

完全正确。连续方程和动量方程在整个流体中是相同的，但密度和粘度根据体积率按混合规律在每个单元确定。体积率α的输运方程假设不可压缩，形式如下：

你的担忧是对的。一阶迎风法的数值扩散很大，界面会崩溃。因此，开发了几何VOF法。例如，OpenFOAM的`interFoam`求解器使用的MULES（多维通用限制器显式解法）通过限制移流通量来保证体积率的有界性（0≤α≤1），同时实现高分辨率。另一方面，Ansys Fluent的几何重构格式通过几何计算流体通过单元界面的体积，是一种更精确的方法。

首先，从当前的体积率分布，用PLIC法求出每个单元内界面的方向（法向量）和位置。然后，根据流速和时间步长，将流向相邻单元的"流体形状"作为多面体进行几何计算。这个多面体的体积本身就是被移流的体积率。计算成本较高，但可以非常精确地保持界面的形状。根据Fluent手册，这个方法对于液滴合并或分裂等剧烈的界面变化也相对稳健。

这是VOF法的一个大课题。表面张力与界面曲率κ成正比，但直接从模糊的α场计算κ会产生很大的噪声。因此使用CSF（连续表面力）模型。这个模型将表面张力作为界面附近的体积力，光滑分布。曲率κ通过先对α场进行光滑化，然后计算其梯度得到的单位法向量

首先，在"多相模型"的选择中选VOF。接下来，准确输入各相的物性值。水的密度998.2 kg/m³，动粘度1.003e-6 m²/s（20°C）。空气1.225 kg/m³，1.4607e-5 m²/s。为了捕捉界面的尺度，界面经过的区域至少需要5～10层细网格，单元尺寸以界面曲率半径的1/10以下为目标。例如，分析直径5mm的气泡，界面附近的网格尺寸应该在0.5mm以下。

根据CFL（Courant-Friedrichs-Lewy）条件确定。特别是为了保证界面移流的稳定性，界面附近的CFL数建议控制在0.25以下。在Ansys Fluent中，设置"自适应时间步进"，将最大CFL数设为0.25，求解器会自动调整时间步长。另外，当表面张力是主导现象时，基于毛细波时间尺度的"表面张力时间步长限制"也可能起作用。

完全同意，沿网格的阶跃状初始界面是不理想的。大多数软件都有"补丁"功能。Fluent中可用"自适应区域"选择水面高度的区域，在其中将水的体积率设为1。更精确的方法是初始化后执行几次"界面光滑化"。或者，将从CAD数据导入的水面形状定义为"壁面"，然后在其一侧填充水。

核心差异在于界面追踪算法和用户界面。

**Autodesk CFD** 也装备了VOF法。重点放在易用性和CAD连接，可以相对容易地开始自由表面流动分析。但对于学术论文基准化的极端精度或复杂现象（激烈的界面破裂等），Fluent或OpenFOAM的算法选项更丰富。
**FreeCAD** 的CfdOF工作台实质上是OpenFOAM的前端，可以从GUI利用`interFoam`的功能。作为免费进行本格VOF分析的优势最大。

取决于算法和并行化效率。几何重构法（Fluent）精度高，但比代数式VOF（OpenFOAM的MULES）每个时间步的计算成本高。但因为能稳定地使用较大的时间步，总计算时间取决于具体情况。而且，Fluent和COMSOL是高度优化的商业求解器，对核数增加的并行化效率非常好。OpenFOAM也可以并行计算，但在大规模计算集群上的设置需要专业知识。

那是"虚假电流"或"寄生流"，是VOF法的经典问题。主要原因是离散化误差导致表面张力平衡被破坏。界面曲率计算有误差时，会产生不应该的虚假压力梯度，界面附近会出现微小的流动。这个流速通常低于毛细速度

有几个对策。
1. **细化网格**：使界面曲率能准确表示。特别是将界面附近设为均匀的直交网格很有效。
2. **增强曲率光滑化**：Fluent中增加"曲率精化"的迭代次数（默认2次，可增加到5次左右）。
3. **更改界面张力模型**：用比CSF模型更高级的SSF（光滑表面力）模型，或重新审视Brackbill模型的"壁面粘附"设置。
4. **强化压力-速度耦合**：使用PISO等对非定常计算强健的耦合解法。同时，将压力修正的松弛系数设得小（约0.1）。

这是个严重的问题。原因和对策如下。
**原因1：时间步长过大** → 严格控制CFL数在0.25以下。
**原因2：压力-速度收敛标准太松** → 将连续方程的残差从1e-4严格控制到1e-5左右。在Fluent中通过"残差监视"确认连续方程的收敛。
**原因3：界面移流格式问题** → 尽可能改用保存性更高的几何重构格式。OpenFOAM中，`interFoam`的MULES相对保存性较好，也可考虑`multiphaseEulerFoam`。
建立定期监视领域内水体积的习惯，确认从初值的变化在数个百分点内。

这是"数值破碎"现象。物理上应该合并的微小界面因数值扩散和移流误差而人工分裂。对策如下：
1. **人工调整界面厚度**：某些求解器有调整界面光滑系数的选项，可增大实际界面厚度。但缺点是界面形状本身会变钝。
2. **使用界面锐化算法**：试用Fluent 2023R1以后"锐利界面建模"这样的、有意保持界面清晰的算法。
3. **根本改变算法**：如果界面分裂合并是主要现象，可能该考虑从VOF法切换到基于粒子的SPH法或界面追踪型Level Set法。