自由表面流
自由表面流的理论基础
概述
老师，自由表面流是什么？它和VOF法不同吗？
自由表面流特有的物理是什么？
表面张力和接触角在小尺度上变得重要。在大尺度上，重力和惯性的比值，即Froude数 $Fr$ 是支配参数。
表面张力的CSF模型
表面张力如何建模？
Brackbill等人（1992）的CSF（连续表面力）模型应用最广泛。从界面曲率 $\kappa$ 计算为体积力。
我听说曲率计算精度是个问题。
CSF的最大问题是寄生流（寄生电流）。来自界面离散化误差的非物理速度场会产生。在毛细管数 $Ca = \mu U / \sigma$ 很小的流动（表面张力主导）中特别严重。
界面追踪法的比较
| 方法 | 界面表示 | 质量守恒 | 界面锐度 | 计算代价 |
|---|---|---|---|---|
| VOF法 | 体积分率 $\alpha$ | 严格 | 依赖方案 | 低-中 |
| Level Set法 | 符号距离函数 $\phi$ | 非守恒（需重新初始化） | 锐利 | 中 |
| CLSVOF | VOF + Level Set | 良好 | 锐利 | 中-高 |
| Phase-Field | 有序参数 | 守恒 | 扩散 | 高 |
| SPH | 粒子 | 守恒 | 粒子分辨率依存 | 高 |
应该选择哪种方法？
在工业应用中，VOF法因成本和精度的平衡而应用最广泛。当界面曲率精度很重要时（微流体等），CLSVOF或Phase-Field法更有优势。在涉及大变形和液体飞散的情况，SPH法也是选项。
自由表面的力学——Bernoulli与毛细管相遇的地方
自由表面（气液界面）是流体力学中物理现象最丰富的地方。重力、表面张力、粘度、惯性全部竞争，它们的比值所表示的无量纲数（Bond数、Weber数、Capillary数、Froude数）的组合决定了支配物理。当Bond数 Bo = ρgL^2/σ 超过1时，重力压过表面张力；反之则毛细管现象占主导。在微流体装置中，Bo << 1，自由表面收缩为球面；在大型船舶波浪中，Bo >> 1，表面张力可以忽略。这种跨越20多个数量级的Bo数多样性，使自由表面CFD成为所有工业中都需要的技术。
自由表面流的数值计算方法
数值求解的详细
请教我自由表面流的数值要点。
界面的Courant数管理和网格设计最重要。在VOF法中，界面Courant数 $Co_\alpha < 0.25$（最大0.5）需要保持。
开放渠道流
河流和水渠流中，Fluent或CFX有开放渠道流的专用设置。入口、出口可以直接指定水位，水面初始化容易。
开放渠道流的设置要点是什么？
波的生成与吸收
在海洋工程的波浪模拟中，需要造波和消波的边界条件。
| 方法 | 概述 | 工具支持 |
|---|---|---|
| Stokes波理论 | 在入口指定波的速度、水位 | Fluent、STAR-CCM+ |
| 五阶Stokes波 | 高阶非线性波 | STAR-CCM+ |
| 消波区 | 用阻尼吸收波 | 全部工具（UDF/Field Function） |
| 数值波浪池 | 造波板 + 消波区 | OpenFOAM (waves2Foam) |
OpenFOAM的waves2Foam很有名吧。
这是Jacobsen等人（丹麦技术大学，2012）开发的库，实现了波的生成、吸收的边界条件和弛豫区。在海洋结构物的波力分析中广泛使用。
寄生流对策
怎样才能抑制寄生流？
| 对策 | 效果 | 备注 |
|---|---|---|
| 网格细分 | 高 | 界面附近的单元要足够小 |
| 高度函数法 | 非常高 | 适合结构网格，提高曲率精度 |
| 锐利表面力 | 高 | 可在STAR-CCM+中使用 |
| CLSVOF | 高 | 用Level Set改进曲率计算 |
| cAlpha调节 | 中 | OpenFOAM的界面压缩参数 |
Level Set法与VOF法的融合——CLSVOF的出现
Level Set法虽然能精确计算界面的几何信息（法线、曲率），但质量守恒不完美。VOF在质量守恒上优秀，但界面形状描述精度较低。2000年代，Sussman等人提出的耦合Level Set-VOF（CLSVOF）法是混合方法，用Level Set函数计算曲率，同时用VOF保证质量守恒。在飞机燃料罐晃荡分析中，CLSVOF比单纯VOF在相同精度下少用了25%的网格，有对标基准的成果。
自由表面流的实务应用
实践指南
请教我自由表面流的分析步骤。
以油罐晃荡分析为例来说明。
1. 形状定义：油罐形状（矩形或圆筒）、液面高度
2. 网格：界面附近细分（液面高度方向至少20个单元）
3. 初始条件：setFields（OpenFOAM）或Patch（Fluent）设定液面
4. 物性值：液体和气体的密度、粘度、表面张力
5. 外力：用Dynamic Mesh或Moving Reference Frame设定油罐振动
6. VOF设置：显式 + Geo-Reconstruct（Fluent）、MULES（OpenFOAM）
7. 时间步：界面Co < 0.25
AMR（适应性网格细化）的应用
我想只在界面附近细分网格。
AMR（自适应网格细化）非常有效。根据体积分率梯度设置细化标准，在界面追踪中动态细化网格。
| 工具 | AMR支持 | 细化标准 |
|---|---|---|
| Fluent | 体积分率的梯度 | 自动 |
| STAR-CCM+ | 基于表格的AMR | 可定制 |
| OpenFOAM | dynamicRefineFvMesh | 基于场 |
| Basilisk | Octree AMR | 原生 |
Basilisk之前在VOF法的文章中也提过。
它是Stephane Popinet（巴黎第六大学）开发的多相流专用开源求解器，集成了octree型AMR。在学术领域评价很高，在自由表面流高精度模拟中有众多成果。
接触角设置
壁面和液面的接触角怎样设置？
静态接触角作为VOF法的壁面边界条件指定。动态接触角模型（前进、后退接触角的差异、速度相关性）也可用。在微流体装置或喷墨等润湿性重要的系统中，接触角设置对结果影响很大。
舰船阻力、耐航性分析——CFD改变了造船
舰船造波阻力（Wave-Making Resistance）占总阻力的30-50%，其预测精度决定了设计的优劣。2000年后，Reynolds平均Navier-Stokes方程+VOF法的CFD成为大型商船设计的必备工具，造波系数预测精度已达5%以内。现代大型LNG油轮设计中，通过CFD优化船首形状（球型船首鼓包），燃油消耗可改善2-4%，年度燃料费差可达数亿元。船级协会已规范了认可的分析条件和验证方法，CFD已从"试验补充"变成"设计主角"。
自由表面流的软件比较
商业工具比较
请比较支持自由表面流的工具。
Flow-3D是什么？
Flow Science公司开发的自由表面流专用CFD求解器。TruVOF法在结构网格上实现高精度界面追踪，在充填、凝固、波浪等自由表面问题中有悠久的业绩。是铸造模拟的业界标准之一。
应用场景推荐
| 应用 | 推荐工具 | 理由 |
|---|---|---|
| 舰船造波阻力 | STAR-CCM+、Fluent | 海洋工程业绩 |
| 油罐晃荡 | Fluent、Flow-3D | 动网格 + VOF |
| 溃坝、洪水 | OpenFOAM、SPH | 大变形流动 |
| 铸造充填 | Flow-3D、Fluent | 充填 + 凝固耦合 |
| 微流体 | COMSOL、OpenFOAM | 表面张力精度 |
| 海洋结构物波力 | OpenFOAM (waves2Foam)、STAR-CCM+ | 波生成、吸收 |
SPH法在什么场景使用？
在津波、碎浪、绿水（船甲板上的海水打击）等界面剧烈变形、液体飞散的问题中，网格型方法困难的领域很有效。DualSPHysics是开源的GPU兼容SPH代码，在港口工程研究中广泛使用。
OpenFOAM vs STAR-CCM+——自由表面分析的实务选择
自由表面流CFD工具选择由成本和功能的权衡决定。OpenFOAM的interFoam免费，拥有世界最大用户群，新算法（isoAdvector等）最先实现。STAR-CCM+（Siemens）强在一体化GUI和CAD联动、自动网格生成，在造船、海上工程的交付周期缩短方面被评价。ANSYS Fluent在通用性和耦合分析（FSI、热传导）的一体化上强。年度许可成本STAR-CCM+约400-600万人民币，Fluent约500-700万，而OpenFOAM为零，但内部培训、技术支持成本不可忽视。
自由表面流的前沿研究
前沿技术与研究动向
自由表面流的最新研究有什么？
我们来看几个方向。
isoAdvector法
Roenby等人（丹麦技术大学，2016）提出的几何VOF法，即使在非结构网格上也能准确重构界面为iso-surface。OpenFOAM v1806后已实装，与传统MULES法相比，界面锐度和质量守恒都有改进。
混合Euler-Lagrange法
连续液体区用VOF法求解，飞散液滴自动转换为DPM（Lagrangian粒子）的混合方法在发展。Fluent 2020后的VOF-to-DPM转换功能是代表。
一次分裂用VOF捕捉，二次分裂用DPM模拟，对吧。
正是如此。反过来，DPM液滴撞击壁面形成液膜的DPM-to-VOF转换也有。可以高效处理多尺度多相流。
GPU计算
大规模自由表面流CFD的GPU加速研究在进行。Fluent 2024后开始支持GPU原生求解器中的VOF法。SPH法因算法特性，与GPU亲和性高，DualSPHysics用一个GPU可实现数百万粒子的实时计算。
机器学习波浪预测
AI也有应用吗？
碎浪的直接模拟——乱流和自由表面的边界
碎浪（Breaking Wave）是自由表面CFD的最难问题之一，波浪前端乱流化并伴随空气卷吸（Air Entrainment）现象。RANS模型对碎浪乱流能量常常高估，需要LES，但自由表面+乱流+气泡生成三者同时求解的计算负荷成倍增加。东京大学研究组用并行VOF-LES直接计算波高1 m级的碎浪，空气卷吸率在实验值±15%以内。海面CO₂交换速度在碎浪中增加2-3倍，对地球科学来说碎浪CFD精度提高是重要课题。
自由表面流的故障排除
故障排除
自由表面流常见的故障有哪些？
一个一个来看。
1. 界面模糊
症状：液面扩散，水和空气的边界看不清。
对策：
- 界面重构方案从HRIC改为Geo-Reconstruct（Fluent）或isoAdvector（OpenFOAM）
- 界面附近网格细化（AMR推荐）
- OpenFOAM中调整 cAlpha 为1.0-1.5
- 界面Courant数控制在0.25以下
2. 寄生流产生
本应静止的液滴周围产生涡…
对策：
- 细分网格（液滴直径至少20单元）
- 使用高度函数法的曲率计算（结构网格）
- Fluent中启用Implicit Body Force
- STAR-CCM+考虑Sharp Interface Method
3. 液面无法初始化
对策：
- Fluent中用Patch操作或Open Channel Flow Initialization
- OpenFOAM中用setFields工具
- 初始压力场用静水压初始化（$p = \rho g h$）
4. Open Channel Flow液面振荡
对策：
- 入口和出口水位保持一致
- 出口附近设置阻尼区
- 时间步足够小
- 启用Implicit Body Force
5. 工具特定的注意事项
| 工具 | 注意点 |
|---|---|
| Fluent | Operating Pressure设置要注意。VOF法推荐设为0（相对压力基准） |
| STAR-CCM+ | VOF的Sharpening因子设置影响界面品质 |
| OpenFOAM | interFoam的alphaEqnSubCycles和nAlphaSubCycles调节 |
| Flow-3D | 仅结构网格。复杂形状用FAVOR法表示，需注意形状分辨率 |
波消失——数值扩散杀死自由表面
自由表面CFD最常见的品质问题是"波传播时振幅衰减"的数值扩散。如果移流方案采用一阶迎风差分，波传播只需数倍波长就会振幅半减。对策是使用高阶精度方案（MUSCL、CICSAM等）和确保充分的网格分辨率（波长方向至少80格子为目标）。不在出口附近设置数值波消波器会导致反射波重新进入计算域污染结果。碰到"波消失"时，首先要按系统进行时间步和网格分辨率的收敛性检验，这是铁则。
帮助
细节
错误