流体仿真服务 — 气液两相流VOF模型的网格与自适应时间步耦合策略

发布时间：2026-07-17 来源：科研学术网

小中大

字号：

为什么VOF仿真要么”假收敛”要么”永远不收敛”

VOF（Volume of Fluid）是CFD中模拟两相流最常用的方法——通过求解体积分数α的输运方程来追踪气液界面的位置和形状。然而在流体仿真服务实践中，VOF模型存在两种非正常的极端行为：

一种是”假收敛”——残差正常下降，看似一切顺利，但监测的物理量（如液膜厚度、液滴直径）完全没有达到稳态。原因是界面区域的Courant数远远高于主流区，时间步太大导致界面被”抹平”——物理上的尖锐界面变成数值上的模糊过渡带。

另一种是”永远不收敛”——施加了小时间步和严格收敛判据后，求解器在高界面曲率区域反复振荡，残差曲线像锯齿一样波动，迭代永远无法收敛到设定的容差。原因是界面的毛细波（capillary wave）在极小的时间和空间尺度上不断产生和耗散，求解器在追逐这些物理上真实的、但数值上难以收敛的小尺度流动特征。

本文复盘一个压力旋流喷嘴（pressure-swirl atomizer）的雾化仿真——燃油以8 bar压力注入，在喷嘴出口形成锥形液膜，随后液膜破碎成液滴。重点分享在ANSYS Fluent中绕过上述两大陷阱的实战策略。

VOF的Courant数控制：远比你以为的复杂

Courant数（Co = u·Δt/Δx）是显式VOF计算中稳定性的黄金法则。在Fluent中，使用显式VOF格式时，界面区域的Courant数必须小于1（推荐<0.25）才能保证界面捕捉的精度。

但问题出在”界面区域”的定义。对于高速雾化流，界面附近的局部速度可以轻松达到100 m/s——远高于主流区速度。而VOF中Courant数的计算只基于主流速度而非界面局部速度，导致在界面附近局部Courant数远超1——你的全局Courant数是0.2，但界面附近实际上是5-10。这就是”假收敛”的根本原因。

我们的解决策略是使用Fluent中的自适应时间步功能：将时间步与全局Courant数绑定（设置Target Global Courant Number = 0.5），让求解器自动在界面区域变密时缩小时间步。同时在界面区域使用High Resolution Interface Capturing（HRIC）格式而非标准的几何重建——HRIC在Courant数略大于1时仍能保持界面的尖锐性，提供了额外的数值鲁棒性。

对于极端界面变形区域（如液膜破碎点），即使HRIC+自适应时间步也不足以精确解析——这些区域天然的物理不稳定性（Rayleigh-Plateau和Kelvin-Helmholtz）要求亚微米级的网格和亚纳秒级的时间步，在工程计算中不可能达到。此时应该接受现实：用亚网格模型（如波动破碎模型）来参数化小尺度雾化过程，而不是试图用VOF直接解析它们。

界面张力与壁面粘附：两个被低估的物理效应

对于毫米级以下的流动，表面张力是主导力。在VOF模拟中，表面张力的计算通过CSF（Continuum Surface Force）模型实现——将表面张力转化为体积力。CSF模型的精度高度依赖界面曲率的计算准确度——而曲率又依赖界面法向量梯度的精度。

这里有一个数值上的恶性循环：界面区域的网格不够细→界面法向量计算不准→曲率计算有误差→表面张力不准确→界面变形预测错误→更差的计算精度。打破这个循环的唯一方法是保证界面区域有足够的网格分辨率——通常要求至少3-5个单元覆盖液膜厚度。对于我们喷嘴案例中约100微米的液膜厚度，这意味着界面区域需要20-30微米的网格——在三维计算中导致数百万的单元数。

壁面粘附（wall adhesion）是另一个容易被忽视的效应。通过指定接触角，VOF可以模拟液体在固体表面的润湿行为。但在高速流动中，动态接触角可能与静态值相差很大——小流量区域的静态接触角和喷嘴出口高流速区域的动态接触角可能差异30°-50°。使用单一静态接触角会导致喷嘴出口处的液膜形状预测偏移。在Fluent中，可以通过UDF（用户自定义函数）实现接触角的动态修正（如Kistler模型），使接触角随界面速度变化——这个看似微小的改进对液膜锥角的预测精度有显著影响。