容器注水模拟，欧拉气液两相流VOF模型，气液两相流算例设置流程。

看着啤酒杯注水时气泡翻涌的画面，突然想用CFD复现这个场景。今天就以OpenFOAM为例，聊聊怎么用VOF模型模拟这种气液两相流现象。别被专业名词吓到，咱们直接上硬菜。

1. 几何建模偷懒指南

先画个二维矩形域模拟水箱（长3m高1m），用blockMesh生成结构化网格。这里有个坑：在气液交界面区域需要加密网格。试过用snappyHexMesh处理复杂几何？记得加个边界层——表面张力对结果影响很大。

// system/blockMeshDict
convertToMeters 0.1;
vertices (...);
blocks ( hex [0 1 2 3 4 5 6 7] (60 20 1) simpleGrading (1 1 1) );

2. 求解器选择玄学

选interFoam就对了。这个求解器自带VOF模型，采用PIMPLE算法处理瞬态问题。注意在phaseProperties里设置水和空气的属性时，粘度参数别手滑写反——我干过把水粘度设为1e-5的蠢事，结果模拟出"神仙水"。

// constant/phaseProperties
phases (water air);
water {
    transportModel Newtonian;
    nu 1e-06;
    rho 1000;
}
air { ... }

3. 初始化骚操作

用setFields定义初始水位，比GUI操作快多了。想要不规则液面？试试UDF自定义初始化。这里有个骚操作：在0目录下同时定义alpha.water和p_rgh，别漏了动量守恒方程需要的初始场。

// system/setFieldsDict
defaultFieldValues ( volScalarFieldValue alpha.water 0 );
regions (... box (0 0 -1) (3 0.2 1) ) { alpha.water 1 }

4. 边界条件里的魔鬼细节

入口用速度入口（velocityInlet），出口建议压力出口（pressureOutlet）。重点注意alpha的边界条件设置——顶部设为"inletOutlet"比固定值更稳定。曾经因为壁面接触角没设置好，模拟出液体爬墙的灵异现象。

// 0/alpha.water
boundaryField {
    inlet { type surfaceNormalFixedValue; }
    outlet { type inletOutlet; inletValue 0; }
    walls { type constantAlphaContactAngle; theta0 70; }
}

5. 时间步长翻车现场

用自适应时间步长能救命。Courant数建议控制在0.5以内，特别是表面张力大的情况。见过因为时间步长太大导致界面破碎吗？像极了倒啤酒时泡沫突然炸开的样子。

// system/controlDict
adjustTimeStep yes;
maxCo 0.5;
maxAlphaCo 0.5;

6. 后处理装逼技巧

用paraView看界面演化时，把alpha=0.5的等值线提取出来最直观。想渲染水面波动效果？试试用alpha场做透明贴图，再用streamLine显示速度场，截图发朋友圈绝对收获一堆"大佬666"。

模拟完成后对比实验数据时可能会发现：实际注水产生更多涡旋。这时候别慌，检查是否忽略了湍流模型。不过对于低速注水，层流模型其实更合适——别为了炫技乱加k-epsilon模型。

避坑总结：

• 表面张力系数别照搬教科书，不同液体差异很大
• 加密网格要适度，否则算到地老天荒
• 初始化液位高度比实际稍高，防止数值耗散导致提前溢出
• 输出频率别太密，否则硬盘比CPU先罢工

最后说个冷知识：用VOF模拟倒啤酒时，适当调低粘度能看到更明显的"挂杯"效果。不过千万别让老板发现你在用超算模拟这个——除非你是啤酒厂CFD工程师。