1. 课程导入与网格质量认知
在计算流体力学(CFD)、结构有限元以及多物理场耦合仿真中,网格是把连续的控制方程离散到计算机可求解代数系统的“载体”。工程实践中,几何模型往往来自CAD,边界形状复杂、曲率变化剧烈,结构化网格难以贴合,因此非结构网格成为主流选择:四面体、棱柱、六面体混合、多面体等单元可以灵活填充任意复杂空间。但灵活性是有代价的——非结构网格的单元形状、尺寸过渡和连接关系由算法自动生成,质量参差不齐,而质量直接决定了离散精度、求解稳定性和计算成本。
很多初学者把网格看作“前处理里点几下就能生成的中间产物”,直到计算发散、残差震荡、壁面热流出现非物理峰值,才回头怀疑网格。事实上,网格质量问题在仿真链条中具有放大效应:一个局部畸变单元可能只影响一个控制体积,但它会污染相邻单元的梯度重构,进而影响整个场的收敛路径。理解这一点,是本课程所有修复技术的认知起点。
非结构网格在仿真流程中的位置
从几何到结果,典型流程是:几何清理 → 网格生成 → 质量检查 → 求解设置 → 迭代计算 → 后处理。非结构网格处在“几何”与“物理”之间的枢纽位置:它既要忠实表达几何特征(如小圆角、薄壁、狭缝),又要满足数值方法对单元形态的约束。以四面体网格为例,求解器通常基于单元体积、面法向和梯度计算构建离散方程,若某个四面体出现极小体积或极大二面角,其面法向计算误差会被放大,导致通量计算失真。
因此,网格质量不是“好看不好看”的问题,而是数值误差的来源之一。工程上常把网格质量分为几何质量与拓扑质量两类:几何质量关注单元自身的形状指标,拓扑质量关注单元之间的连接合理性,例如是否存在孤立单元、重复面、非流形边。二者共同决定网格能否被求解器稳定使用。
质量问题的典型后果
为了建立直观认识,可以看一个教学示例(以下为合理假设的示例场景,非实测数据):某管道内流场采用四面体网格,全局尺寸设为 8 mm,弯头处局部加密到 2 mm。质量检查发现约 0.3% 的单元最小二面角低于 10°,集中在弯头内侧。直接求解时,连续性残差在前 200 步下降,随后出现周期性震荡,监测点速度波动约 15%。将这批单元局部重构、把最小二面角提升到 18° 以上后,残差单调下降,监测点波动收敛到 2% 以内。这个例子说明:少量低质量单元足以改变整个计算的收敛行为。
更严重的后果包括:负体积导致求解器初始化失败;高偏斜单元使对流项出现非物理过冲;边界层棱柱网格长宽比失控,导致壁面剪切应力预测偏差;多面体网格中相邻面夹角过小,引起矩阵条件数恶化。排障时若只调整松弛因子或时间步长,往往治标不治本,根源仍在网格。
课程目标与学习路径
本课程的目标不是罗列软件按钮,而是建立一套可迁移的判断与处理能力:看到质量报告,能判断哪些指标真正危险;面对低质量单元,能选择合适的修复手段而不是盲目全局重划;修复之后,能用同一套指标验证是否达标。围绕这一目标,课程按“认知—诊断—修复—验证”的主线展开。
第一阶段建立质量认知:理解各类质量指标的几何含义与数值影响,明确不同求解器、不同物理模型对指标的敏感度差异。例如结构分析更关注雅可比行列式与长宽比,而外流CFD更关注边界层棱柱的增长率与壁面法向尺寸。第二阶段进入诊断:学会阅读质量直方图、定位最差单元、区分“数量少但致命”和“数量多但温和”的问题。第三阶段是修复技术:光顺、拓扑修改、局部重划、边界层重构等,每种方法都有适用边界。第四阶段是验证与复盘:修复后重新检查指标分布,并通过短程试算确认收敛行为改善。
需要强调的是,网格质量没有统一的“万能阈值”。同一个偏斜度指标,在有限体积法中可能容忍度较高,在谱元法中却可能直接导致失败。因此课程中给出的数值(如二面角下限、长宽比上限)均作为教学参考区间,实际工程应以求解器文档和网格无关性验证为准。学习路径的终点,是形成“先看指标分布、再定位局部、后选择手段、最后闭环验证”的工作习惯,而不是记住某个固定数字。