📚 课程目录第4章网格收敛性验证

4. 网格收敛性验证

网格收敛性验证是判断数值解是否“足够接近真实解”的第一道可信度关口。它的核心逻辑并不复杂：在保持物理模型、边界条件、材料参数和求解器设置不变的前提下，只改变网格密度或网格类型，观察目标量（如最大温度、最大应力、流量、电容、谐振频率等）是否随网格加密而趋于稳定。如果目标量在连续两次或多次加密后变化小于预设容差，就可以认为该网格下的解已进入收敛区间；反之，则说明当前网格仍不足以解析关键梯度，需要继续细化或调整网格策略。

网格类型选择：先看物理场，再看几何特征

COMSOL 中常用的网格类型包括自由四面体、自由六面体（映射/扫掠）、棱柱/金字塔过渡层以及边界层网格。选择时不要只凭习惯，而应同时考虑物理场类型和几何拓扑。结构力学中，四面体适合复杂曲面和倒角，但一阶四面体偏“刚”，容易高估刚度，通常建议使用二阶单元；六面体在规则拉伸体、薄板、长杆中效率更高，相同自由度下精度往往更好。传热与流体中，边界层网格决定壁面热流和剪切的解析能力，第一层厚度需与热边界层或黏性底层尺度匹配。电磁波问题中，网格尺寸通常按最高频率对应的波长控制，一般要求每波长至少 5～10 个单元，且曲面上需保证足够采样。

多物理场模型还要注意“网格服务于谁”。例如焦耳热问题中，电流密度在电极边缘集中，温度梯度也在此处最大，网格应优先加密电极边缘和导体拐角；流固耦合中，流体侧边界层和固体侧应力集中区往往不在同一位置，需要分别设置局部细化，而不是全局均匀加密。

局部细化策略：把自由度花在梯度大的地方

全局加密的代价是自由度按三维几何近似三次方增长，而局部细化可以把计算量集中在真正需要的地方。COMSOL 中常用的局部细化手段有：在几何上添加“细化域”或“网格控制域”；对特定边界设置“尺寸”节点并指定最大单元尺寸；对边或点设置更小尺寸；使用“边界层”节点在壁面附近生成各向异性单元；对扫掠网格设置分布单元数。实际操作中，建议先做一次粗网格求解，查看通量、应力、温度梯度或电流密度的分布，再回到网格节点对高梯度区域加密。这样比一开始就盲目全局细化更高效。

局部细化还要避免“孤岛式”突变。相邻区域单元尺寸比不宜过大，一般控制在 2～3 倍以内，否则会在过渡区引入新的数值误差。对于薄层结构，如涂层、薄膜或接触层，应使用扫掠或边界层网格保证厚度方向有足够单元，而不是让四面体在厚度方向只落一层。

网格收敛指数计算：用 GCI 给出量化判据

网格收敛指数（Grid Convergence Index, GCI）是 CFD 和传热领域常用的收敛性量化方法，也可迁移到其他物理场。其基本流程是：选取至少三套系统性加密的网格，计算目标量 f，用细网格结果估计离散误差。设三套网格的单元尺寸分别为 h₁ < h₂ < h₃，对应目标量为 f₁、f₂、f₃，加密比 r = h₂/h₁。先计算收敛阶 p：

p = ln[(f3 - f2)/(f2 - f1)] / ln(r)

GCI_fine = Fs * |(f2 - f1)/f1| / (r^p - 1)

其中 Fs 为安全因子，三套网格比较时通常取 1.25。

若 f₂−f₁ 与 f₃−f₂ 同号且量级合理，说明解处于渐近收敛区间；若出现振荡或反号，则说明网格尚未进入收敛区，需要继续加密或检查局部奇异性。工程上常把 GCI 小于 5% 作为可接受判据，但具体阈值应结合目标量用途确定：用于趋势判断可放宽到 10%，用于定量设计则建议控制在 2%～3%。

COMSOL 网格统计工具：看什么、怎么判断

COMSOL 的“网格统计”窗口会给出单元数、最小单元质量、平均单元质量、体积/面积比、偏度等指标。单元质量接近 1 表示形状规则，接近 0 表示畸变严重。四面体网格中，若最小质量低于 0.1，通常意味着存在狭长或扁平单元，可能导致刚度矩阵病态和收敛困难。此时应检查几何中的小面、尖角、短边，或调整网格尺寸和增长速率。另一个重要工具是“网格收敛图”：在“研究”中启用“网格收敛”研究步骤，COMSOL 会自动生成多套网格并绘制目标量随自由度变化的曲线，可直接观察是否趋于平台。对于多物理场，建议分别监控每个场的代表性目标量，而不是只看一个总残差。

多场耦合下的网格匹配与收敛判据

多物理场模型最容易出现“单场收敛、耦合不收敛”的情况。原因通常有两个：一是不同物理场对网格尺度要求不同，共用一套网格时某一侧解析不足；二是耦合变量在界面上的传递对网格匹配敏感。例如流固耦合中，流体压力传递到固体表面，若固体侧网格远粗于流体侧，压力插值会引入额外误差；焦耳热中，电流密度和温度梯度都集中在电极边缘，若只加密温度场而忽略电流场，焦耳热源本身就不准。

因此，多场模型的收敛判据应分层设置。第一层是各物理场自身的目标量收敛，如最高温度、最大应力、压降、电容等，分别做网格收敛检查；第二层是耦合量收敛，如流固耦合中的界面位移、焦耳热中的总焦耳热功率、压电中的电荷输出等。只有当各场目标量和耦合量都进入稳定区间，才能认为网格收敛性验证通过。若某一耦合量始终随网格加密单调漂移，往往提示存在几何奇异性、接触非线性或边界条件不匹配，需要回到模型设置排查，而不是继续盲目加密。

实用建议：在 COMSOL 中可先建立“网格收敛”研究步骤，选择 3～4 套网格，指定需要监控的全局量（如 max(T)、max(solid.mises)、comp1.force 等），求解后直接查看收敛曲线。若曲线已趋于水平，再取最细网格结果作为参考解；若曲线仍明显下降，则需继续加密或改用更高阶单元。

常见误区：把“网格质量高”等同于“网格收敛”。单元质量只反映单元形状是否规则，不反映网格密度是否足够。一个质量很好但整体过粗的网格，仍然可能给出严重偏离的结果。收敛性验证必须通过系统性加密和结果对比来完成，不能只看统计窗口中的质量指标。

最后需要强调，网格收敛性验证不是一次性的收尾工作，而应贯穿建模过程。建议在模型定型前就确定目标量和容差，例如“最高温度变化小于 0.5 K”或“最大应力变化小于 2%”，然后按此标准选择最终网格。这样既能避免过度加密带来的计算浪费，也能防止网格不足导致结论不可信。对于多场耦合模型，还应把各场目标量和耦合量列成检查清单，逐项确认收敛，再进入参数扫描或优化设计阶段。