第三讲:网格生成技术基础
各位同学好,我是老张。今天咱们聊聊网格生成。说实话,网格这东西,看着不起眼,但往往是决定CFD计算成败的关键。我见过太多人,模型建得漂亮,边界条件设得合理,结果算出来一塌糊涂——十有八九是网格出了问题。
网格生成,说白了就是把连续的空间切成一个个小单元。你想想看,计算机没法直接处理连续的流体,它只能处理离散的点。网格就是这些点的集合。网格质量好不好,直接决定了你的计算结果靠不靠谱。
核心观点:网格是CFD计算的骨架。骨架歪了,肉长得再好看也没用。
3.1 网格类型:结构网格 vs 非结构网格 vs 混合网格
网格类型的选择,是我每次做项目时第一个要纠结的问题。三种主流类型各有千秋,我一个个说。
3.1.1 结构网格
结构网格,就是那种整整齐齐、像棋盘一样的网格。每个节点都有规则的编号,邻居关系一目了然。
- 优点:计算效率高,内存占用少,边界层处理效果好
- 缺点:对复杂几何适应性差,生成耗时
- 适用场景:简单几何、翼型、管道等
我个人习惯,能用结构网格的地方绝不用别的。为什么?因为结构网格的计算精度和收敛速度,非结构网格真比不了。我记得有一次做二维翼型分析,用结构网格跑了500步就收敛了,换成非结构网格,同样的设置跑了2000步还在振荡。
3.1.2 非结构网格
非结构网格就灵活多了。三角形、四面体、多面体,想怎么切就怎么切。
- 优点:几何适应性强,自动化程度高
- 缺点:计算效率低,内存占用大,边界层处理困难
- 适用场景:复杂几何、多部件装配体
我建议,遇到复杂的几何,比如风机叶片带扭曲、带弯掠的,老老实实用非结构网格。别硬着头皮画结构网格,那是在跟自己过不去。
3.1.3 混合网格
混合网格,说白了就是取长补短。近壁面用结构网格,远离壁面用非结构网格。
嗯,这里要注意:混合网格的交接面处理是个技术活。我曾经在一个项目中,因为交接面没处理好,导致计算结果出现虚假的涡结构,排查了整整三天才找到原因。
我的建议:对于风机翼型分析,推荐使用混合网格。近壁面用O型结构网格捕捉边界层,远场用非结构网格减少网格量。这样既保证了精度,又控制了计算成本。
3.2 网格质量评价指标
网格画完了,怎么判断好不好?三个核心指标:正交性、长宽比、偏斜度。我一个个讲。
3.2.1 正交性
正交性,衡量的是网格线之间的夹角。理想情况下,网格线应该互相垂直,夹角90度。
为什么会这样?因为正交性差的网格,会导致数值扩散。说白了,就是计算出来的结果被"抹平"了,该有的细节没了。
| 正交性范围 | 质量等级 | 建议 |
|---|---|---|
| 0.9 - 1.0 | 优秀 | 理想状态 |
| 0.7 - 0.9 | 良好 | 可接受 |
| 0.5 - 0.7 | 一般 | 需要改进 |
| < 0.5 | 差 | 必须重新生成 |
警告:我曾经在一个项目中,正交性最低只有0.3,结果计算出来的升力系数比实验值高了15%。重新画网格后,误差降到了3%以内。正交性差,计算结果就是废纸。
3.2.2 长宽比
长宽比,就是网格单元最长边和最短边的比值。对于边界层网格,长宽比可以很大,比如100:1甚至更高。但对于主流区,长宽比最好控制在5:1以内。
你想想看,一个网格单元被拉得特别长,就像一根面条。流体在这个单元里流动,数值误差会被放大。我建议,主流区的长宽比控制在3:1以内,边界层可以放宽到100:1。
3.2.3 偏斜度
偏斜度,衡量的是网格单元偏离理想形状的程度。理想情况下,偏斜度为0。偏斜度越大,网格质量越差。
- 0 - 0.25:优秀
- 0.25 - 0.5:良好
- 0.5 - 0.75:一般
- 0.75 - 1.0:差
嗯,这里有个坑要注意:偏斜度高的网格,会导致收敛困难,甚至发散。我见过有人算了一个星期,结果因为几个偏斜度0.9的网格,残差死活降不下去。
3.3 网格独立性验证
网格独立性验证,说白了就是回答一个问题:你的计算结果跟网格数量有没有关系?
为什么会这样?因为理论上,网格越密,计算结果越接近真实解。但网格太密,计算成本又太高。我们需要找到一个平衡点。
我个人的做法是这样的:
- 粗网格:先画一套粗网格,跑一遍,看看趋势
- 中等网格:加密一倍,再跑一遍
- 细网格:再加密一倍,跑第三遍
- 对比:比较三个结果,如果中等网格和细网格的差异小于1%,就认为网格独立了
技巧:做网格独立性验证时,重点关注你关心的物理量。比如做翼型分析,就看升力系数和阻力系数。这两个量对网格最敏感。
我记得有一次做风机叶片分析,粗网格算出来升力系数1.2,中等网格1.35,细网格1.38。粗网格和细网格差了15%,但中等网格和细网格只差了2.2%。最后我选了中等网格,既保证了精度,又控制了计算时间。
网格独立性验证,不是走过场。我见过太多人,随便画一套网格就开始算,结果论文被审稿人质疑。你想想看,审稿人一问"你的网格独立了吗?",你拿什么回答?
好了,这一讲的内容就到这里。网格生成是CFD的基础,也是决定成败的关键。下一讲,我们会进入求解器设置,到时候再聊。
公众号:蓝海资料掘金营,微信deep3321