3、网格划分问题：网格质量差导致发散，局部加密策略与网格无关性验证
网格这东西，说白了就是仿真计算的「骨架」。骨架歪了，肉长得再好看也没用。我见过太多新手一上来就猛画网格，结果算到一半直接发散，气得拍桌子。其实很多发散问题，根源就在网格质量上。
3.1 网格质量差为什么会发散？
你想想看，CFD求解器本质上是在解偏微分方程。网格质量差，意味着单元的形状扭曲、长宽比过大、偏斜度太高。这些畸形的网格会导致数值误差急剧放大，就像用一把歪尺子量东西，越量越离谱。
我个人习惯把网格质量指标分成三个等级：
| 指标 | 优秀 | 可接受 | 危险区 |
|---|---|---|---|
| 偏斜度 (Skewness) | < 0.5 | 0.5 - 0.85 | > 0.85 |
| 正交质量 (Orthogonal Quality) | > 0.7 | 0.3 - 0.7 | < 0.3 |
| 长宽比 (Aspect Ratio) | < 5 | 5 - 20 | > 20 |
我在项目中遇到过一台服务器散热仿真，网格偏斜度到了0.92，结果残差曲线直接飞上天。后来把那个区域的网格重新画了一遍，偏斜度降到0.6以下，收敛就稳了。嗯，这里要注意：偏斜度超过0.95的网格，基本就是定时炸弹。
⚠️ 避坑指南： 我曾经在IGBT模块仿真中，因为忽略了边界层网格的正交质量，导致结温计算结果偏差了15%。后来花了整整两天排查，才发现是近壁面网格太歪了。从那以后，我每次画完网格第一件事就是检查正交质量。
3.2 局部加密策略：把钱花在刀刃上
网格不是越密越好。全模型加密，计算量翻倍，精度提升却有限。说白了，局部加密才是性价比最高的方案。
我个人常用的加密策略有这几种：
- 温度梯度加密：热源附近、散热器翅片根部、芯片结区，这些地方温度变化剧烈，必须加密。
- 几何特征加密：小圆角、窄缝隙、薄壁结构，这些地方网格容易扭曲，需要局部细化。
- 边界层加密：壁面附近速度梯度大，第一层网格高度要满足y+要求。我一般用这个公式估算：
第一层网格高度 = (y+ × 特征长度) / (0.199 × Re^(7/8))举个例子，一个风冷散热器仿真，我通常这样设置加密区域：
- 芯片表面：全局尺寸的1/5
- 翅片间隙：至少3层网格
- 入口区域：保持均匀网格，避免数值耗散
- 出口区域：适当稀疏，节省计算资源
💡 小技巧： 我习惯用「曲率加密」功能，它能自动识别几何曲率大的地方进行加密。但要注意，曲率加密的阈值别设太低，否则会在平面区域也生成大量无用网格。
3.3 网格无关性验证：别被假收敛骗了
网格无关性验证，说白了就是回答一个问题：你的结果还依赖网格吗？
我记得有一次做LED灯具热仿真，用粗网格算出来结温85°C，加密一次变成92°C，再加密一次变成88°C。当时我就懵了——这波动也太大了。后来一直加密到第四套网格，结果才稳定在89.5°C左右。
网格无关性验证的标准流程，我总结为三步：
- 生成3-5套网格：从粗到细，网格数量按1.5-2倍递增。比如10万、20万、40万、80万。
- 监控关键指标：温度、压力降、流量这些工程关心的量。我一般选3-5个监测点。
- 判断收敛标准：当相邻两套网格的结果偏差小于1%-2%，就认为网格无关了。
这里有个坑：不要只看全局量。我曾经只看总热阻，觉得网格无关了，结果局部温度分布还是有问题。所以我现在都会同时检查几个关键点的温度值。
✅ 网格无关性验证的实用建议：
- 用GCI（网格收敛指数）定量评估，GCI < 5% 算合格
- 关注梯度大的区域，比如热源正上方
- 如果计算资源有限，至少做两套网格对比
- 记录每套网格的计算时间和内存占用
3.4 本章知识体系
下面这张图是我自己整理的网格问题排查逻辑，你可以照着这个思路走：
3.5 实战中的几个小建议
最后分享几个我踩过的坑，希望能帮你少走弯路：
- 别迷信自动网格：自动网格生成器确实方便，但复杂几何下经常出问题。我习惯先用自动网格跑一遍，然后手动检查关键区域。
- 网格过渡要平滑：粗细网格之间体积比不要超过1.3，否则会在过渡区产生数值振荡。
- 保存网格版本：每次修改网格都保存一个版本，方便回溯。我曾经改了一版网格后结果变差，还好有备份能回退。
- 多核并行时注意：有些求解器在并行计算时对网格分区有要求，分区不好会导致负载不均，反而更慢。
网格这东西，做多了就有感觉了。刚开始可能觉得麻烦，但等你吃过几次发散的亏，就会明白——花在网格上的时间，永远是最值得的投资。