2.3 网格质量诊断与网格无关性验证

2.3 网格质量诊断与网格无关性验证

本节摘要：网格好不好要分开回答两个问题：网格本身质量是否过关（单元质量、偏斜度、长宽比），以及网格上算出的结果是否稳定（加密后关键量漂移是否收敛）。本节给出诊断指标的操作路径、0.95 这条偏斜度硬底线，以及一套四级网格无关性验证流程与判据。

两个不同的问题，别混在一起答

上一节末尾的统计表已经给出了平均单元质量 0.93 这样的数字，但"数字好看"与"结果可信"之间还隔着两道关卡。第一道是网格自身的几何品质：畸形单元会让局部方程病态，轻则精度损失，重则求解器直接罢工。第二道是离散误差：任何网格都是对连续空间的近似，网格不够细时，算出的最高温度、最大应力本身就带着系统性偏差。两道关卡都要过，且检查方法完全不同——前者看单个单元的形状指标，后者看结果随网格加密的走势。

质量指标：读什么、底线在哪

COMSOL 把每个单元的形状品质折算成一个零到一的标量，一代表完美，零代表退化到不可用。常用的统计口径有三项：平均单元质量反映整体水准，工程模型的合理区间在 0.5 以上；最小单元质量暴露最差的短板，低于 0.1 就要警惕，它通常与求解困难相伴；偏斜度（skewness）从相反方向度量畸变，零为理想、一为完全退化，业界通行的硬底线是最大偏斜度低于 0.95——超过这条线，意味着网格里存在接近塌缩的单元，任何求解器都可能在它身上栽跟头。

诊断操作固定成三步。第一步，右键"网格"节点→统计，弹出面板里逐项记录单元数、平均单元质量、最小单元质量与最大偏斜率，养成每次划网必看的习惯。第二步，把差单元找出来：在结果里新建一个"网格"绘图（三维绘图组→网格），颜色表达式选单元质量，色标设成从红到蓝，红色单元在图形窗口里无所遁形；放大红色集中区，十有八九在薄缝隙、小倒角或导入几何的碎面附近。第三步，对症下药：畸变源于几何就回 2.1 修几何或用虚拟操作压碎短边；源于尺寸突变得太猛就在局部加"尺寸"子节点平滑过渡，相邻单元尺寸增长控制在 1.5 倍以内；实在无解的个别单元，用"剖分单元"操作手工切开重划。

图：网格质量判定指标与处理路径

网格无关性验证：让结果不再随网格漂移

质量指标全绿，只说明网格"长得好"，不说明解已经稳定。离散误差随网格加密单调下降，验证的思路因此很直接：加密网格、重算、看关键工程量还动不动。什么时候必须做？出报告、写论文、做认证，一律必做；自己摸趋势、比方案的初步试算，可以豁免但心里要有数。

图：四级网格无关性验证流程

落地到散热器算例，流程走起来不超过半小时。第一步，选"芯片最高温"作关键量——它是第 3 章报告的主角，别拿"平均温度"这种不敏感的量糊弄自己。第二步，在用户控制网格的"尺寸"节点里只改一个全局缩放因子：1.4、1.0、0.7、0.5 四档，对应单元数约 2.6 万、5.2 万、10.8 万、21.4 万，单元数大致逐级翻倍。第三步，逐档重算，把单元数与最高温填进表格；注意除网格外什么都不许动，求解器设置一改，验证就失真了。第四步，判读：中到细变化 0.18%，细到超细变化 0.05%，远小于 2% 的判据，说明细网格上的最高温已经稳了，正式计算取细网格，验证表存档。若算例是应力集中或湍流这类对网格敏感的问题，变化率不会这么漂亮，可能要加密到第四、第五级才达标——判据不变，只是耐心要多一些。

本节要点回顾

• 两问分开答：网格质量看单元形状指标，结果可信度看加密后的漂移，两张卷子分开做。
• 三条指标线：平均单元质量大于 0.5、最大偏斜度低于 0.95、最小单元质量大于 0.1，配套增长率 1.5 以内。
• 诊断三步：统计面板读数 → 网格绘图按质量着色定位红块 → 按先几何后尺寸再手工的顺序修复。
• 无关性流程：选定关键量 → 四级缩放因子只动网格 → 逐级记录 → 相邻变化小于 2% 判稳，验证表进报告附录。
• 验证记录本身就是结果可信度的证据，第 9 章的验证与确认流程会把它归入证据链。

补一段网格无关性验证的完整操作示范，因为这是仿真可信度的基石。方法：从粗到密做三到四组网格，每组网格的自由度数至少差 50%，记录关键位置的温度值，画"温度 vs 自由度数"曲线。曲线应该先陡降后走平——走平的起点就是网格无关解，之后的网格再加密结果几乎不变。以散热器为例：粗网格（5000 单元）测得芯片温度 52 ℃，中网格（15000）测得 48 ℃，细网格（40000）测得 47.5 ℃，特细（80000）测得 47.4 ℃——47.5 与 47.4 差 0.1 ℃（<1%），说明 40000 单元已是网格无关解。在报告中只需要给出这张"温度 vs 自由度"表格加一句"后续仿真采用细网格设置"，审稿人或评审者就知道你做过验证——这一张表的成本是一小时的计算时间，买到的是结果的可信度认证。

补一个网格质量诊断的具体操作路径。COMSOL 的网格统计可以在"网格"节点右键→统计中查看：最小/平均/最大单元质量（0-1，1 最好）、最小角度、单元数与自由度数。质量低于 0.1 的单元会引起求解精度问题甚至不收敛——定位方法是右键网格→"按质量绘图"，颜色图上红色区域就是问题单元。修正低质量单元的手段按优先级：调整几何（删除过短边或过小面）→调整网格尺寸（在问题区域加 Size 子节点局部加密或放粗）→换网格类型（自由四面体换扫掠网格或映射网格）。修正后重新检查质量——这个迭代过程通常 2-3 轮就能达到所有单元质量大于 0.3 的合格线。记住：网格质量低于 0.1 的单元即使求解成功，其结果也是不可信的——一个烂单元就足以污染它周围的解。

网格无关性验证的报告方式也有讲究。在工程报告或学术论文中，标准的做法是附一张"网格收敛表"：列出每种网格设置的总单元数、自由度数、关键物理量（如最大温度、最大应力或流量）的数值，以及与最细网格的偏差百分比。偏差从粗到细逐级递减并在最细两级之间小于某个阈值（工程上常用 1-2%），说明已经达到网格无关。这张表的另一个用途是估算计算成本：从表中可以看出达到目标精度所需的最少自由度数，从而选择最经济的网格设置用于后续的参数化扫描或优化迭代——参数化扫描可能需要跑几十到几百次，每次省下一半计算量就是巨大的效率提升。网格无关性验证不是可选的学术仪式，而是仿真结果能否被信任的必要证据——没有做过验证的仿真，在工程评审中会被质疑甚至直接否决。