📚 课程目录第2章网格基础

网格基础：网格类型与生成工具

网格是什么？说白了，就是把连续的空间切成小块。这些小块用来近似物理场。我最早接触这个概念时，觉得很简单，不就是画格子嘛。但后来发现，格子的画法直接影响仿真准不准、算得快不快。嗯，这就是咱们这门课的基础。

网格类型决定了数值离散的方式。网格生成工具决定了你能不能把思路变成实际的网格。这两个东西，是参数敏感度分析的起点。

什么是网格

你想想看，一个连续的物体，比如机翼，空气流过它时压力是连续变化的。但计算机只能处理离散的数字。所以，我们得把这机翼分成很多小单元，每个单元内用一个简单的函数近似物理量。这就是网格——空间的离散化。

网格分两大类：结构化与非结构化。我在项目中看到过不少新手，上来就用非结构化，说它方便。但其实，结构化网格在某些场景下优势明显。不理解这两者的区别，后面的敏感度分析容易跑偏方向。

提示：网格参数敏感度分析的核心，就是看参数变化（比如网格尺寸、拓扑）对结果的影响。不理解网格类型，你都不知道该动哪里。

结构化网格

结构化网格，说白了就是每个节点都有规则的拓扑连接。像棋盘一样，每个点都能通过(i,j,k)索引定位。典型的单元形状是六面体（3D）或四边形（2D）。

优点：

内存占用低。只需要存储节点坐标，连接关系隐含在索引里。
计算效率高。因为数据访问是连续的，CPU缓存友好。
数值精度好。尤其是对于流动、传热等梯度大的问题，六面体比四面体更准。

缺点：

几何适应能力差。复杂的模型，比如带倒角、孔洞的零件，很难生成纯结构化网格。
手工工作量大。我最早接触结构化网格是在做气动分析的时候，光分块就花了两周。

避坑指南：我曾经在划分结构化网格时吃过亏——某个区域网格歪斜严重，导致计算发散。后来才发现，分块接口处的节点没对齐。结构化网格对拓扑一致性要求极高。

生成结构化网格，可以用 blockMesh（OpenFOAM）、ICEM CFD 里的结构化分块功能，或者 Gambit 的分块设计。我最常用的逻辑是用映射法：先画二维拓扑，再拉伸成三维。

# 一个简单的结构化网格生成逻辑（Python 伪代码示意）
nx, ny, nz = 20, 10, 30
dx, dy, dz = 0.1, 0.05, 0.2
for i in range(nx):
    for j in range(ny):
        for k in range(nz):
            x = i * dx
            y = j * dy
            z = k * dz
            # 生成节点 (x, y, z)

非结构化网格

非结构化网格，节点之间的连接没有固定模式。单元可以是四面体、三棱柱、金字塔，甚至混合使用。说白了，就是更自由，但代价是存储连接关系需要额外数组。

优点：

几何适应能力强。复杂模型，比如带有曲面、细小特征，很容易自动生成网格。
自动化程度高。很多工具一键生成，不需要手工分块。

缺点：

内存占用高。因为要显式存储每个单元的节点编号。
计算效率稍低。数据访问不连续，缓存命中率下降。
数值精度可能会差一点。尤其是四面体，在某些流动问题上会产生伪扩散。

我第一次用非结构化网格，是被逼出来的。那是分析一个带有复杂散热齿的元件，用结构化网格分块，我得疯了。Gmsh 一键生成四面体，省了大力气。但后来发现，网格质量差的地方，温度结果跳来跳去。

避坑指南：我曾经在非结构化网格上踩过一个坑——某个区域单元尺寸变化太剧烈（从1mm直接跳到0.01mm），导致相邻单元体积差了几个数量级，求解器直接报错。非结构化网格对尺寸过渡要求很敏感。

# 使用 Gmsh 生成非结构化网格的命令行示例
gmsh -3 -algo front2d -clscale 0.1 my_model.step -o my_mesh.msh

网格生成工具

工具是干活的家伙。我用过不少，从商业软件到开源都有。个人习惯，简单列几个常见的：

工具	类型	主要特点	适用场景
Gmsh	开源	轻量，支持多种算法，有 Python API	快速原型、教学、小模型
Salome	开源	集成 CAD 修复、网格划分，功能全面	复杂几何、多物理场耦合
cfMesh / snappyHexMesh	开源 (OpenFOAM)	基于八叉树迭代切割，自动化程度高	大规模流体网格、汽车外流场
Ansys Meshing	商业	强大的自动网格生成 + 手动控制	工业级仿真，质量稳定

为什么会有这种差异？说白了，工具的背后是算法不同。Gmsh 对 Delaunay 算法支持好，snappyHexMesh 走的是迭代切割路线。你选工具，先得看你的几何类型和网格需求。

个人经验：我建议你，用 Gmsh 做技术验证，用 snappyHexMesh 做大规模计算。Salome 适合需要 CAD 修复的场景，但学习曲线陡。不要贪多，一个工具吃透就够了。

我的选择建议

网格类型和工具的选择，其实是个权衡。我一般按这个逻辑走：

几何简单？用结构化网格，精度和效率都高。
几何复杂但区域规整？考虑混合网格：核心区域结构化，外部非结构化。
完全复杂的曲面模型？老老实实上非结构化网格，用 Delaunay 或前沿推进法。
工具选择？新手从 Gmsh 或 Ansys Meshing 起步，老手再碰开源框架。

我曾经见过有人用非结构化网格算一个简单的管道流动，结果四面体太多，内存爆了。其实用结构化网格，一个 block 就解决。所以，别盲目跟风。