VOF法（流体体积法）多相流分析

正是。海浪仿真、储罐晃动（液体晃动）、液滴下落、溃坝——这些有"自由界面"或"多流体界面"的问题都是VOF法的主战场。VOF（Volume of Fluid，流体体积法）用体积分数 $\alpha$ 表示每个网格单元中各相流体的占比，隐式追踪界面。1981年由Hirt和Nichols提出，至今仍是多相流CFD的标准方法之一。

对，穿越界面的网格单元里 $0 < \alpha < 1$。如何让这些"界面单元"保持清晰锐利是精度的关键，这就是"界面重构格式"要解决的问题。Geo-Reconstruct（PLIC）方式在每个单元内把界面重构为一条直线，能保持非常清晰的界面。

关键看气泡尺寸与网格尺寸的比值。VOF要求气泡直径至少覆盖3～5个网格单元，才能分辨界面。大气泡（直径>几mm）、Taylor气泡、液弹流（Slug Flow）都适合VOF。微米级弥散气泡（烟气洗涤、浮选槽）则不适合，应用欧拉双流体模型（Euler-Euler）或Mixture模型。工程里经常看到混用：大气泡用VOF，弥散小气泡用Mixture Model。

VOF的核心是体积分数 $\alpha$ 的输运方程：

$\alpha = 1$ 代表流体1（如水），$\alpha = 0$ 代表流体2（如空气）。混合密度和粘度用 $\alpha$ 线性插值：

N-S方程只解一套（共享速度场）。OpenFOAM的interFoam求解器在 $\alpha$ 输运方程中加入压缩项 $\nabla \cdot (\mathbf{u}_r \alpha(1-\alpha))$，防止界面扩散变模糊。

CSF（连续表面力，Continuum Surface Force）模型把表面张力作为体积力处理：

$\sigma$ 是表面张力系数，$\kappa = -\nabla \cdot \hat{n}$ 是界面曲率，$\hat{n} = \nabla\alpha / |\nabla\alpha|$ 是界面法向量。CSF的弱点是"寄生流（Spurious Currents）"——界面曲率的离散化误差会在静止液滴周围产生非物理涡流。加密网格或用Height Function法高精度计算曲率可以缓解。

这两个无量纲数帮你判断什么力主导：

$We \ll 1$：表面张力主导，液滴保持球形（毛细现象）；$We \gg 1$：惯性主导，液滴破碎变形。$Bo \ll 1$：表面张力大于重力，液滴呈球形（毛细管中的小液滴）；$Bo \gg 1$：重力主导，液体为平坦自由液面（储罐、海浪）。这些参数直接指导你是否需要精细处理表面张力。

各主要界面重构格式的对比：

必须控制界面Courant数——一个时间步内界面移动不能超过一个网格单元，否则 $\alpha$ 输运会出错：

实务中推荐用自适应时间步长，让界面Courant数始终保持在0.5以内。表面张力很强时（小液滴、毛细现象），还需满足毛细Courant数限制 $Co_{cap} = \Delta t / \sqrt{\rho \Delta x^3 / (2\pi\sigma)} \le 1$，时间步要更小。

VOF的质量守恒验证是基本质量控制：

• 总体积分率守恒：监控 $\int \alpha \, dV$ 的时间变化，变化率应小于初始值的0.1%。OpenFOAM仿真结束后检查 log 文件中的 Phase-1 Volume fraction。
• 界面锐利度：后处理中查看 $\alpha$ 等值面（等值面0.5），界面宽度（$\alpha$ 从0.1到0.9的法向距离）不应超过2～3个网格单元。
• 寄生流检查：静止液滴案例中，最大速度应小于物理速度尺度的1%。否则需要改善曲率计算方法。

液体晃动是VOF法的经典应用，几个关键点：

1. 初始液位设置：用patch命令或场函数精确初始化 $\alpha$ 分布，避免初始界面模糊。
2. 网格加密：在预期液面波动范围内（约±20%液位高度）加密网格，Geo-Reconstruct法要求界面处至少3～5层网格。
3. 自适应时间步长：设 $Co_{max} = 0.5$，让求解器自动调整 $\Delta t$。
4. 壁面润湿性：接触角（Contact Angle）对毛细现象有影响，纯水对金属约60°～80°，可通过UDF自定义。
5. 结果验证：与晃动固有频率（$f = \frac{1}{2\pi}\sqrt{g\pi/L \cdot \tanh(\pi h/L)}$）的理论值对比。

液滴撞击问题涉及高速动态接触线，是VOF的高难度应用：

• 超细网格：液滴直径至少20～30个网格，接触线附近还要进一步加密（AMR）。
• 动态接触角：速度相关的动态接触角模型（Kistler, Cox-Voinov）对铺展/回缩精度至关重要。
• 表面张力精度：Weber数低（$We < 10$）时必须用Height Function法计算曲率，否则寄生流会主导流动。
• 飞溅预测：$We > 80$（近似临界值）时液滴会飞溅，需要超细网格（液滴直径50+网格）才能准确预测。

来看主要工具的VOF功能对比：

VOF的优势是质量守恒性好，Level Set在界面法向量和曲率计算上更精确。目前热门的研究方向包括：

• CLSVOF（耦合Level Set与VOF）：兼顾两者优点，界面清晰且质量守恒。
• 自适应网格细化（AMR）：界面附近自动加密，其他区域用粗网格，大幅节省计算量。
• 机器学习辅助界面重构：用神经网络预测界面位置，比传统PLIC更高效。
• Phase-Field方法融合：适合纳米尺度润湿与毛细现象的精细模拟。

激光粉末床熔融（LPBF）的熔池VOF仿真是当前研究热点。挑战在于：时间尺度跨度极大（激光扫描微秒→冷却毫秒），界面要同时处理液-气界面（表面张力）和液-固界面（凝固）。通常做法：VOF追踪液-气界面 + 焓法处理固液相变 + Marangoni热毛细对流。这类多物理场耦合问题在COMSOL和Ansys Fluent都有专门的增材制造模块。