计算流体力学仿真：网格无关性验证与湍流模型选择决策

计算流体力学仿真的第一课项目组是在管壳式换热器上学到的：用粗网格算出的壳侧压降与细网格结果差了30%。客户问”到底哪个对”，答案是细网格的更接近实验值。但粗网格不是”错”，是”没收敛”——网格密度不够时，数值扩散掩盖了物理扩散。从此项目组建立了一条铁律：任何CFD结果必须附网格无关性验证，否则不予交付。

网格无关性验证：四套网格的系统性对比

管壳式换热器的几何特征：壳径200 mm，管束含37根φ19管，正三角形排列，管间距25 mm。折流板间距100 mm，25%切割率。

项目组构建了四套网格：

监测指标：壳侧总压降、管侧总压降、总传热系数、壳侧最大流速。

M1到M3的壳侧压降变化显著（降低约30%），M3到M4的变化缩小到6%。M3和M4的传热系数差异仅1%——M3的320万网格已基本收敛。项目组选择M3作为生产网格，兼顾精度和计算效率。

y+的要求与边界层网格

y+（无量纲壁面距离）是CFD中最容易被忽视、却影响最大的参数。它决定了边界层内流场的解析质量：

• y+ < 1（粘性底层解析）：适用于k-ω SST和LES，需要极细的近壁网格
• 1 < y+ < 5（缓冲区）：应避免，湍流模型在此区域不确定
• 30 < y+ < 100（对数律区）：适用于k-ε模型配合壁面函数

项目组在M3网格中第一层高度0.05 mm，Re=12000（壳侧流速约0.8 m/s，水），计算y+≈1.2——刚好在粘性底层。这允许使用k-ω SST模型而无需壁面函数，精度最高。

但如果用k-ε模型，y+=1.2反而不好——k-ε的壁面函数要求y+>30。这时需要调整第一层网格高度到约1.2 mm（y+≈30），网格层数减少但每层更粗。项目组为两种湍流模型各构建了对应的网格，而非一个网格跑所有模型。

湍流模型对比

管壳式换热器的壳侧流动特征复杂：折流板间的交叉流、管束后的尾涡、泄漏流。项目组对比了三种湍流模型：

k-ε Realizable模型：

• 优点：收敛性好，计算量小，工业应用最广
• 缺点：对弯曲流线和分离流的预测精度差
• 壳侧压降：偏高约20%
• 流场特征：折流板后的分离区偏大，回流强度偏弱

k-ω SST模型：

• 优点：近壁处理精确，对分离流敏感
• 缺点：计算量比k-ε大50%，对入口湍流条件敏感
• 壳侧压降：偏高约6%
• 流场特征：分离区大小合理，尾涡频率与实验PIV数据吻合

LES（大涡模拟）：

• 优点：直接解析大尺度涡，无需湍流模型假设
• 缺点：计算量是RANS的100倍，需要非定常计算
• 壳侧压降：偏差约2%
• 流场特征：捕捉到了管束间的瞬态涡脱落，时均结果与SST接近

LES的精度最高但计算成本不可接受（单工况跑2周）。k-ω SST在精度和效率之间取得了最佳平衡——6%偏差在工程设计可接受范围内（通常±10%）。

边界条件与收敛判据

管壳式换热器CFD的边界条件设置：

• 壳侧入口：质量流量入口（0.8 kg/s，293K水）
• 壳侧出口：压力出口（0 Pa表压）
• 管侧入口：速度入口（1.5 m/s，353K热水）
• 管侧出口：outflow
• 管壁：耦合壁面（coupled，同时参与壳侧和管侧传热）
• 壳壁和折流板：绝热壁面

收敛判据的设置是另一个关键点。项目组不仅看残差，还监测关键物理量的稳定性：

• 残差标准：连续性<1E-4，动量<1E-5，能量<1E-7，k和ω<1E-4
• 监测点：壳侧出口平均压力、管侧出口平均温度、总传热率
• 收敛判定：残差达标 AND 监测点量在500步内变化<0.5%

如果只看残差，可能在物理量未稳定时就误判收敛——这在复杂几何（如管壳式换热器）中常见。项目组在M3网格上跑4000步后残差达到标准，但壳侧出口压力仍在缓慢下降（每100步降0.2%）。继续跑到8000步后压力稳定，才确认收敛。

后处理：数据的可视化与量化

CFD后处理的核心输出：

1. 流线图：展示壳侧流体的流动路径，识别短路流和死区
2. 压力云图：沿流动方向的压力分布，定位主要压损源
3. 温度场截面：管壁附近的温度梯度，计算局部传热系数
4. 速度矢量图：折流板转弯处的流速分布，识别高速冲刷区

量化输出包括：壳侧压降、管侧压降、总传热系数、壳侧传热系数、管侧传热系数、传热有效度ε。

反思：CFD不是”黑箱”

CFD仿真最大的风险是”调参到想要的结果”。项目组见过不少案例——为了让CFD结果与实验吻合，人为调整湍流模型参数或网格密度。这不是工程，是数字游戏。

正确的方法论是：先做网格无关性验证确定网格，再根据流动特征选择湍流模型，最后与实验数据对比验证。如果偏差>15%，先检查几何简化是否合理、边界条件是否准确，而不是调模型参数。更多CFD仿真实战经验，可以参考有限元仿真栏目，或返回科研学术网首页。