流体仿真 — 从湍流模型选择到网格无关性验证的实战复盘

流体仿真（CFD）在工程中的应用极其广泛——汽车风阻优化、飞机气动设计、泵阀内部流动、建筑风载荷分析、电子设备散热，几乎每个工程领域都离不开它。做了大量CFD项目后，我最大的感触是：CFD的门槛看似不高（软件会自动求解），但做好CFD的门槛极高——选错一个湍流模型，结果可能差30%；边界层网格少一层，阻力预测可能偏20%。本文从实际项目出发，系统复盘流体仿真中湍流模型选择、网格策略、求解器配置和结果验证等核心环节的实战经验。

一、流体仿真的技术背景与求解体系

1.1 CFD的基本原理

流体仿真的核心是求解Navier-Stokes方程组——连续性方程、动量方程和能量方程。对于不可压缩牛顿流体：

– 连续性：∇·v = 0

– 动量：ρ(∂v/∂t + v·∇v) = -∇p + μ∇²v + f

– 能量：ρc_p(∂T/∂t + v·∇T) = k∇²T + Φ

其中v是速度场，p是压力，ρ是密度，μ是动力粘度，Φ是粘性耗散函数。这些方程是非线性的耦合偏微分方程，没有解析解（除了极少数简单情况），必须通过数值方法求解。

数值方法主要有三种：有限体积法（FVM）、有限元法（FEM）和有限差分法（FDM）。工程CFD中FVM是绝对主流——ANSYS Fluent、OpenFOAM、Star-CCM+都基于FVM。FVM的优势是天然满足局部守恒性，对守恒律方程有天然的适配性。

1.2 求解器类型

CFD求解器按处理方式分为两大类：

压力基求解器（Pressure-Based）：以压力为基本变量，通过SIMPLE/PISO算法耦合速度和压力。适用于不可压缩和低速可压缩流动。ANSYS Fluent的默认求解器。适合大多数工程应用（马赫数Ma < 0.3视为不可压缩）。

密度基求解器（Density-Based）：以密度为基本变量，直接求解连续性方程获得密度。适用于高速可压缩流动（Ma > 0.3），特别是超音速和高超声速流动。在密度变化大的流动中精度更高。

选型经验：

– 不可压缩流动（水、低速空气）：压力基

– 亚音速可压缩（Ma 0.3-0.8）：压力基（耦合求解器）

– 跨/超音速（Ma 0.8-5）：密度基

– 高超声速（Ma > 5）：密度基+真实气体效应

1.3 稳态vs瞬态

稳态计算：求解不含时间导数的方程。适用于达到稳定状态的流动（如管道稳态流、翼型定常绕流）。计算效率高，是设计迭代的默认选择。

瞬态计算：求解含时间导数的方程。适用于非定常流动（如旋翼流动、涡脱落、瞬态启动过程）。计算量大，时间步选择至关重要。

判断标准：如果流动的Strouhal数St=fL/U显著大于0（存在明显的特征频率），说明流动是非定常的，需要瞬态计算。例如圆柱绕流在Re>40时出现涡街（St≈0.2），必须用瞬态计算才能捕捉。

二、湍流模型选择与配置

2.1 湍流模型概述

湍流是CFD中最核心也是最具挑战性的问题。雷诺数Re=ρUL/μ超过临界值（管流约2300，外流约5×10^5）后，流动从层流转为湍流，出现多尺度涡结构。直接数值模拟（DNS）求解所有尺度的涡，计算量与Re^(9/4)成正比，工程应用中完全不现实。

大涡模拟（LES）直接求解大尺度涡，对小尺度涡用亚格子模型。计算量约为DNS的1/100，但仍需千万级网格和毫秒级时间步，目前仅用于科研。

工程CFD几乎都使用雷诺平均Navier-Stokes（RANS）方法，将瞬时速度分解为时均值和脉动值：v = V + v’。脉动项通过湍流模型封闭。RANS模型的计算量与层流计算相当，是工程实践的标准选择。

2.2 RANS模型选型

主流RANS模型按复杂度递增：

Spalart-Allmaras（SA）一方程模型：只求解一个湍流运动粘度修正变量。计算速度快，对壁面流动（如航空外流场）表现好。弱点是对分离流动预测不准。适合航空航天外流场。

k-ε两方程模型：求解湍动能k和湍流耗散率ε。标准k-ε模型适用于远离壁面的高Re流动，但在壁面附近需要壁面函数。Realizable k-ε和RNG k-ε是改进版本，对分离流和旋流有更好表现。适合内部流动（管道、泵、阀门）。

k-ω两方程模型：求解湍动能k和比耗散率ω。在壁面附近不需要壁面函数（可以直接积分到壁面），对低Re和转捩流动表现好。弱点是对自由流来流条件敏感。

SST k-ω（Shear Stress Transport）模型：结合k-ε和k-ω的优点——壁面附近用k-ω，远离壁面用k-ε。通过混合函数实现平滑过渡。对逆压梯度下的分离预测显著优于其他两方程模型。这是工程CFD中最推荐的通用湍流模型。

RSM（Reynolds Stress Model）：直接求解雷诺应力张量的6个分量，不使用Boussinesq涡粘假设。对各向异性湍流（强旋流、强曲率流动）表现最好。计算量大（多求解5个方程），收敛困难。适合需要精确预测各向异性应力的场合。

选型决策表：

| 流动类型 | 推荐模型 | 理由 |

|———|———|——|

| 航空外流场 | SA或SST | 壁面流动主导，计算效率高 |

| 管道/泵/阀门 | Realizable k-ε | 内部流动，高Re |

| 翼型失速/分离 | SST | 逆压梯度分离预测好 |

| 强旋流（旋风分离器） | RSM | 各向异性显著 |

| 低Re流动 | SST（低Re修正） | 可以积分到壁面 |

| 热混合/自然对流 | SST或RNG k-ε | 浮力效应处理较好 |

2.3 壁面处理与y+策略

壁面附近的湍流行为对整体结果影响极大。壁面处理策略取决于y+值（无量纲壁面距离）：

壁面函数法（y+ > 30）：第一个网格节点位于对数律区域（30 < y+ < 300），用壁面函数近似壁面剪应力。计算效率高，适合高Re流动。不需要极细的边界层网格。Fluent默认使用Scalable Wall Functions。

壁面积分法（y+ < 1）：第一个网格节点位于粘性子层内，直接解析壁面附近的速度梯度。精度最高但需要极细的边界层网格（首层高度可能只有几微米）。适用于低Re流动、需要精确壁面热流的场合。

y+目标值配置：

– 外流场（航空/汽车）：y+ = 30-100，用壁面函数

– 内流场（泵/涡轮）：y+ = 30-50，用壁面函数

– 热传递分析：y+ < 1，壁面积分

– 分离流：y+ < 1，壁面积分（SST模型在y+<1时分离预测最准）

首层网格高度估算：

y₁ = y+ * μ / (ρ * U_τ)

其中摩擦速度U_τ = √(τ_w/ρ)，τ_w = 0.5 * C_f * ρ * U²。摩擦系数C_f的估算：平板湍流C_f ≈ 0.026/Re_x^(1/7)。

例如：空气（μ=1.8×10^-5 Pa·s，ρ=1.225 kg/m³），来流速度30m/s，平板长度1m，Re≈2×10^6：

– C_f ≈ 0.003

– τ_w ≈ 0.5 × 0.003 × 1.225 × 900 ≈ 1.65 Pa

– U_τ ≈ √(1.65/1.225) ≈ 1.16 m/s

– 目标y+=50时：y₁ = 50 × 1.8×10^-5 / (1.225 × 1.16) ≈ 0.63mm

2.4 实战案例：汽车外流场仿真

在一个乘用车气动减阻项目中，需要将风阻系数Cd从0.32降到0.28。

模型配置：

– 求解器：ANSYS Fluent，压力基，稳态

– 湍流模型：Realizable k-ε（初始选型）→ SST k-ω（优化阶段）

– 计算域：车前3倍车长，车后5倍车长，上方3倍车高，侧面3倍车宽

– 边界条件：入口速度33.3m/s（120km/h），出口压力出口，地面移动壁面（模拟道路）

– 网格：Poly-Hexcore混合网格，约1500万网格

初始计算结果：Realizable k-ε模型预测Cd=0.35，实测0.32，偏差9.4%。切换到SST k-ω后Cd=0.328，偏差2.5%。

网格优化：

– y+分析：初始网格y+范围80-250，部分区域y+>300超出对数律区域

– 加密边界层：首层高度从0.5mm降到0.15mm，增长比1.2，15层→20层

– y+降到30-80，Cd=0.326，偏差1.9%

尾流优化：车尾分离区是阻力的主要来源。在车尾增加精细网格（加密到5mm），配合SST模型准确捕捉分离涡。最终Cd=0.283，偏差0.9%。

经验总结：汽车外流场的湍流模型选择对Cd预测精度影响极大。SST k-ω在分离流区域的优势是Realizable k-ε无法比拟的，但需要更细的边界层网格（y+<50）才能发挥优势。

三、网格策略与网格无关性验证

3.1 网格类型选择

CFD中常用网格类型：

结构化网格（Hexahedral）：网格排列规则，计算效率最高，数值扩散最小。适合简单几何。但复杂几何的六面体网格划分极其困难。

非结构化网格（Tetrahedral）：适应复杂几何，自动划分。但计算效率低于六面体（相同精度下网格数量多2-5倍），数值扩散较大。

多面体网格（Polyhedral）：ANSYS Fluent的特色网格类型。面数多（平均14面），相邻单元连接更丰富，计算效率比四面体高30-50%，精度也更好。适合复杂几何。

Poly-Hexcore（切割体网格）：Fluent Meshing的Watertight工作流输出的主流网格类型。远场用六面体（效率高），近壁和复杂区域用多面体（适应性好）。兼具两者的优势。

选型建议：优先使用Poly-Hexcore，兼顾效率和适应性。对于简单几何（如管道），可以用结构化六面体。纯四面体网格仅作为最后选择。

3.2 网格密度与局部加密

CFD网格密度的决定因素是流场梯度——梯度大的区域需要加密。关键加密区域：

壁面附近：边界层加密，按y+目标值设计首层高度和层数。

分离/尾流区：在预期分离点和尾流区域加密。分离区网格尺寸应小于涡特征尺寸的1/10。对汽车/飞机，尾流区网格尺寸建议为物体特征尺寸的1/50到1/100。

剪切层：射流、混合层等强剪切区域需要加密。网格应能分辨剪切层厚度（通常为物体尺寸的1-5%）。

激波区域：超音速流动中激波前后参数急剧变化，需要在激波位置加密。自适应网格加密（AMR）是处理激波的有效手段。

3.3 网格无关性验证

网格无关性验证是CFD结果可信度的基本保证。标准流程：

1. 生成三套网格：粗、中、细，网格数量比例约1:3:9。例如100万、300万、900万
2. 计算关键指标：选择对工程决策最重要的物理量（如Cd、压降、最大温度、质量流量）
3. Richardson外推：用Richardson外推法估计网格无关解和收敛阶

Richardson外推的核心公式：f_exact ≈ f_fine + (f_fine – f_medium)/(r^p – 1)，其中r是网格加密比（通常2），p是收敛阶。

网格收敛指数（GCI）：GCI = Fs * |(f_fine – f_medium)/(r^p – 1)| / f_fine，Fs为安全因子（1.25）。GCI < 5%表示网格基本无关，<1%表示高度无关。

实战经验：

– 不要只看全局指标（如Cd），还要检查局部量（如某点的压力系数Cp、某截面的速度剖面）。全局指标可能提前收敛，局部量还在变化

– 网格加密应该是系统性的——同时加密所有方向，不能只加密一个方向

– 对于瞬态计算，网格无关性验证要在稳态或代表性时间步上进行

3.4 边界层网格参数

边界层网格是CFD网格中最关键的部分。核心参数：

– 首层高度y₁：由y+目标值确定

– 增长比（Growth Ratio）：相邻层高度比，推荐1.15-1.2。>1.3会导致网格质量下降和数值误差增大

– 层数：壁面函数法10-15层，壁面积分法15-25层

– 总厚度：应覆盖边界层厚度δ的80%以上。湍流边界层δ ≈ 0.37L/Re^(1/5)

检查方法：在后处理中查看y+分布。y+的最大值应不超过目标范围的上限（壁面函数法<300，壁面积分法<1）。如果局部y+超限，需要在该区域调整边界层参数。

四、求解器收敛与瞬态计算

4.1 收敛判断标准

CFD计算的收敛判断比结构分析复杂得多。不能只看残差下降，需要多维度判断：

残差标准：

– 连续性残差：降到1E-4以下（工程精度）或1E-5以上（高精度）

– 动量残差：降到1E-4以下

– 能量残差：降到1E-6以下

– 湍流残差：降到1E-4以下

残差下降到目标值并不一定意味着收敛——特别是对于复杂流动，残差可能在某处振荡。需要结合以下指标判断：

监控量稳定性：设置关键物理量的监控点（如出口流量、某点压力、阻力系数），当这些量在连续100-500次迭代中变化小于0.1%时，可以认为收敛。

质量守恒检查：入口质量流量与出口质量流量之差应小于0.1%。如果不满足，说明收敛不充分或有泄漏面。

能量守恒检查：对于有热传递的流动，入口能量与出口能量之差应等于壁面热交换量。

4.2 不收敛的常见原因与对策

原因一：网格质量问题。Skewness > 0.95或Aspect Ratio > 1000的单元会导致局部发散。对策：在Fluent Mesh Quality中检查并改善低质量单元。

原因二：边界条件不当。出口设在回流区域会导致不收敛。对策：将出口移到远离回流区（至少10倍特征长度下游），或使用压力出口代替质量流量出口。

原因三：时间步过大（瞬态）。CFL条件不满足导致发散。对策：减小时间步，使CFL < 5（隐式）或CFL < 1（显式）。CFL = U*Δt/Δx。

原因四：欠松弛因子过大。对策：降低欠松弛因子。推荐值：压力0.3，动量0.5-0.7，湍流0.8，能量0.8-1.0。

原因五：物理模型不匹配。如用稳态求解器计算本质非定常的流动（如大分离区涡脱落）。对策：切换到瞬态求解。

4.3 瞬态计算的时间步选择

瞬态CFD的时间步选择直接影响精度和计算效率：

CFL条件：Δt < CFL * Δx / U_max。对于隐式求解器，CFL可以大于1但不建议超过10（精度下降）。对于显式求解器，CFL必须<1。

特征频率法：如果流动有主导频率f（如涡脱落频率），时间步应满足Δt < 1/(20f)，即一个周期至少20个时间步。

时间步无关性验证：类似于网格无关性，用三个时间步（比例1:2:4）计算，比较关键物理量。当时间步减半时结果变化<1%即可。

内迭代次数：每个时间步内的迭代次数通常设为20-50。如果达到最大迭代次数仍未收敛（残差未降到1E-4），说明时间步过大或网格有问题。

4.4 实战案例：离心泵瞬态性能预测

在一个离心泵（比转速n_s=120）的性能预测中，需要计算不同流量下的扬程和效率曲线。

模型配置：

– 求解器：ANSYS Fluent，压力基，瞬态

– 湍流模型：SST k-ω

– 网格：约500万，结构化六面体，叶轮和蜗壳分别网格化，用滑移网格界面连接

– 转速：2900 rpm

– 时间步：Δt = 5.77×10^-5 s（叶轮每转1°的时间），每步30次内迭代

计算策略：

1. 先做稳态计算（MRF方法，多重参考系），获得初始流场
2. 切换到瞬态滑移网格，计算5个完整旋转周期（1725步），确保周期性收敛
3. 在第3-5个周期内时间平均，获得扬程和效率

结果对比：

– 设计流量下扬程：仿真32.1m，实测31.5m，偏差1.9%

– 最高效率：仿真76.2%，实测74.8%，偏差1.9%

– 偏离设计流量（0.5Q_design）时扬程：仿真38.5m，实测36.2m，偏差6.4%

经验总结：

1. 稳态MRF方法在设计流量附近精度可接受（偏差<3%），但在偏离工况偏差大
2. 瞬态滑移网格在所有工况下精度都更好，但计算量是稳态的20-50倍
3. 蜗壳隔舌附近的网格密度对结果影响很大——该区域流动极度复杂（冲击、分离、二次流），网格尺寸应小于1mm
4. 时间步选择1°/步是叶轮机械的标准配置，2°/步精度损失约2%，0.5°/步精度提升不明显但计算量翻倍

五、结果验证与精度控制

5.1 CFD验证的层次

CFD结果的验证需要多层次对标：

基准实验验证：选择标准算例（如NACA 0012翼型、突扩管道、方柱绕流）验证模型设置的正确性。这些算例有大量公开实验数据，偏差应控制在5%以内。

缩比实验验证：用缩比模型的风洞/水洞实验数据验证。注意雷诺数效应——缩比模型的Re可能与全尺寸不同，需要修正。

全尺寸验证：用实际产品的测试数据验证。这是最终验证，但数据有限且条件不可控。

5.2 常见精度问题与对策

数值扩散（Numerical Diffusion）：网格与流动方向不对齐时，截断误差产生额外的数值粘性，导致涡量被耗散。对策：对准网格与流向（结构化网格优势）；使用高阶离散格式（Second-Order Upwind或MUSCL）。

离散格式选择：

– First-Order Upwind：最稳定但精度最低，数值扩散严重。仅用于初始计算

– Second-Order Upwind：标准选择，精度和稳定性平衡好

– QUICK：三阶格式，对结构化网格精度更高

– MUSCL：三阶格式，适用于非结构化网格

推荐策略：先用一阶格式计算1000步获得收敛初始场，再切换到二阶格式继续计算。能量方程始终用二阶。

入口湍流参数：入口边界需要指定湍流强度I和湍流长度尺度L：

– 湍流强度：I = 0.16 * Re^(-1/8)（管内充分发展湍流）

– 长度尺度：L = 0.07D（管径D），或L = 0.1δ（边界层厚度δ）

入口湍流参数不当会影响边界层转捩和分离预测，特别是在翼型前缘附近。

5.3 CFD报告的关键验证指标

一份合格的CFD报告应包含以下验证信息：

1. 网格信息：网格类型、总网格数、最小/最大网格尺寸、y+分布直方图
2. 网格无关性：三套网格的对比结果和GCI值
3. 收敛历史：残差曲线和监控量曲线
4. 守恒检查：质量/能量平衡的数值
5. 与实验对比：关键物理量的对比表和误差分析
6. 敏感性分析：湍流模型对比、离散格式对比（如有）

六、专业流体仿真服务

需要流体仿真服务？

科研学术网提供专业的流体仿真服务：

✅ 博士级工程师团队，一对一技术支持

✅ 计算结果可靠，可提供详细的技术报告

✅ 周期灵活，加急项目最快3天交付

✅ 价格透明，无隐形费用

立即咨询报价 →