2. Fluent多相流建模的几何与网格策略

多相流仿真的精度上限，往往在网格生成阶段就已经被决定了。与单相流不同，多相流计算域中同时存在连续相与离散相，相间界面（自由液面、气泡表面、颗粒边界、液膜界面）会随时间迁移、变形甚至破碎，网格必须既能解析连续相的边界层与剪切层，又要在界面经过的位置保持足够的分辨率。本节围绕计算域简化、多区域网格划分、边界层处理、网格质量对界面捕捉的影响四条主线展开，并给出一个工业级网格案例的完整配置思路。

一、计算域简化：先做减法，再做网格

工业几何往往包含大量与流动无关的细节：法兰螺栓、仪表接管、外部支撑、倒角、小孔、铭牌凹槽。这些特征会迫使网格在这些位置加密，导致单元数暴涨而物理收益为零。简化的基本原则是：凡是不影响相分布、不改变主流通道截面积、不产生显著局部阻力的特征，都可以删除或等效替代。

具体操作上，可按以下顺序处理：

• 对称性利用：若几何与边界条件关于某平面近似对称，且多相流本身不破坏该对称性（如竖直管内均匀气液两相流），可只建一半或四分之一模型，用 symmetry 边界替代。但若存在重力驱动的分层、旋流或非对称进料，则不可简化。
• 小特征删除：直径小于主流通道特征尺寸 1/20 的接管、支管、测压孔，若仅用于监测而不影响流动，可直接删除；若必须保留，则用等效阻力元件或一维模型替代。
• 入口段与出口段延长：多相流入口处相分布往往不均匀，出口处可能存在回流。建议入口延长 3～5 倍水力直径，出口延长 5～10 倍水力直径，避免边界条件直接干扰界面形态。
• 三维到二维的降维判断：对于长直管道、轴对称喷嘴、平面液膜等，若相分布沿某一方向基本不变，可降为二维轴对称或二维平面模型，计算量可降低一到两个数量级。但涉及气泡三维摆动、颗粒旋转、液滴破碎等本质三维现象时，不可降维。

二、多区域网格划分：按物理需求分区

多相流计算域内不同区域的物理特征差异很大：壁面附近需要边界层网格解析速度梯度，自由液面附近需要各向同性细网格捕捉界面变形，主流核心区可用较粗网格降低计算量。因此，分区划分、分区控制尺寸是多相流网格策略的核心。

典型的分区方式包括：

• 近壁区：采用棱柱层（prism layer）或六面体边界层网格，第一层高度由目标 y+ 反算，增长率控制在 1.1～1.2，层数 10～20 层。
• 界面活动区：自由液面可能经过的区域，采用各向同性六面体或切割体网格，尺寸取气泡/液滴直径的 1/10～1/20。若使用 VOF 模型，界面附近网格尺寸直接决定界面厚度和虚假扩散程度。
• 主流核心区：尺寸可放宽至界面区的 2～4 倍，但相邻区域过渡比不超过 1.2，避免数值扩散突增。
• 离散相轨迹区：若使用 DPM 或欧拉-拉格朗日方法，颗粒/液滴可能经过的路径需保证网格尺寸大于颗粒直径的 2～3 倍，否则颗粒会被“卡”在网格中，导致轨迹失真。

分区之间的交界面应尽量与流动方向垂直或平行，避免斜交导致插值误差。在 Fluent 中，可通过 Face Meshing、Body Sizing、Inflation 等局部控制实现分区，也可借助 Fluent Meshing 的 Watertight Geometry 工作流中的“Local Sizing”和“Boundary Layer”任务完成。

三、边界层处理：多相流中的特殊考量

边界层网格的核心参数是第一层高度 Δy 和层数 N。第一层高度由目标 y+ 决定：

Δy = y+ · ν / uτ
uτ = sqrt(τw / ρ)

其中 ν 为运动黏度，uτ 为摩擦速度，τw 为壁面剪应力，ρ 为密度。对于 k-ε 模型，通常要求 y+ 在 30～300 之间，使用壁面函数；对于 k-ω SST 或雷诺应力模型，要求 y+ ≈ 1，直接解析黏性底层。

多相流中边界层处理有三个特殊问题：

• 相变与壁面沸腾：若壁面存在沸腾、冷凝或蒸发，边界层内温度梯度极大，需额外加密温度边界层，第一层高度按热边界层厚度估算，通常比速度边界层更薄。
• 壁面润滑与液膜：使用欧拉壁膜模型时，壁膜厚度可能只有几十微米，边界层网格需能解析液膜厚度，否则液膜速度、厚度预测严重失真。
• 颗粒与壁面碰撞：DPM 中颗粒碰撞壁面时，边界层网格过粗会导致碰撞点位置和反弹速度计算误差增大，建议碰撞区域局部加密。

四、网格质量对相间界面捕捉的影响

相间界面的捕捉精度直接受网格质量支配。以下指标对多相流尤为关键：

以 VOF 模型为例，界面厚度通常被数值扩散控制在 2～4 个网格单元内。若网格尺寸为 1 mm，则界面物理厚度约 2～4 mm；若实际液膜厚度只有 0.5 mm，则网格根本无法分辨，液膜会被虚假扩散到多个单元中，导致液膜速度、传热传质计算完全失真。此时必须将界面区网格加密到 0.1 mm 以下，或改用欧拉壁膜模型专门处理薄液膜。

五、工业级网格案例：气液分离器

以一个典型卧式气液分离器为例，直径 1.2 m，长度 3.5 m，入口为顶部切向进料，气相从顶部出口排出，液相从底部出口排出。操作压力 0.8 MPa，气相表观速度 0.15 m/s，液相表观速度 0.02 m/s。

网格策略如下：

• 计算域：保留完整筒体与进出口接管，删除外部支撑、法兰螺栓、仪表接口。入口延长 0.5 m，出口延长 1.0 m。
• 分区：将筒体分为入口射流区、液池区、气相分离区、出口区四个区域。入口射流区尺寸 8 mm，液池区 10 mm，气相分离区 15 mm，出口区 12 mm。
• 边界层：壁面第一层高度 0.3 mm，增长率 1.15，层数 12 层，目标 y+ ≈ 50，配合 Realizable k-ε 模型。
• 界面加密：液池自由液面附近 0.3 m 高度范围内，网格尺寸加密至 5 mm，确保液面波动和气泡夹带可分辨。
• 网格类型：主体采用六面体切割体网格（cutcell），边界层采用棱柱层，总单元数约 320 万。

质量检查结果：最低正交质量 0.32，最大纵横比 85，最大偏斜度 0.78，满足多相流计算要求。若将液池区网格粗化到 20 mm，则自由液面波动被抹平，气相出口带液量预测值偏低约 40%，说明界面区加密不可省略。

该案例说明，多相流网格不是“越密越好”，而是“该密的地方必须密，该粗的地方可以粗”。入口射流区、自由液面区、边界层区是三个必须重点控制的区域，其余区域可在保证过渡比的前提下适当放宽。最终网格质量指标应作为仿真报告的一部分，与结果一同交付。