3.3 网格质量评估指标
3.3 网格质量评估指标
网格造出来只是半成品，本节做质检：用偏斜度、正交质量、长宽比与 y+ 把网格放进验收单的第一栏。这些数字不是官僚手续——每一个都对应第二章讲过的离散误差机制，质检不合格的网格会让最好的格式与模型一起失效。本节还负责一件贯穿全书的事：网格无关性验证的标准流程，它是第四关确认验收的核心证据。
指标红线的物理出处
偏斜度（Skewness）。衡量单元偏离理想形状的程度，零是完美，一是最差。它杀伤求解器的机制在面通量：歪斜单元的面法向与单元中心连线夹角过大，梯度重构失真，矩阵对角占优被破坏，轻则收敛慢重则发散。工程红线：最大值不超过零点九五（必须修），平均值最好低于零点三三。
正交质量（Orthogonal Quality）。单元中心到面的法向向量与面中心向量的夹角余弦类指标，一为完美、零最差。它对应扩散通量计算的干净程度——Fluent 的扩散通量默认用正交修正两步法，正交质量低于红线时修正量膨胀，数值噪声随之而来。工程红线：最小值不低于零点一，边界层网格可适度放宽。
长宽比（Aspect Ratio）。最长边与最短边之比。边界层网格故意拉长（比值上百属正常），但主流区大长宽比单元意味着方向分辨率失衡。判断原则：只警惕"出现在不该出现位置"的大长宽比。
y+ 与边界层分辨率。湍流近壁面第一层网格的无量纲高度，是湍流模型近壁策略与网格的接口：壁面函数路线要求第一层节点落在对数律区（y+ 约三十到三百），低雷诺数解析路线要求 y+ 约等于一。它错配的后果在第四章展开，这里先记住：y+ 是网格与湍流模型之间的合同条款，违约一方是网格。
| 指标 | 合格 | 警戒 | 红线 |
|---|---|---|---|
| 偏斜度最大值 | 低于 0.9 | 0.9 到 0.95 | 大于 0.95 必须修 |
| 偏斜度平均值 | 低于 0.33 | — | — |
| 正交质量最小值 | 高于 0.15 | 0.1 到 0.15 | 低于 0.1 必须修 |
| 主流区长宽比 | 低于 100 | 按区域评估 | 局部评估 |
| y+（壁面函数） | 30 到 300 | 越界复核 | 与模型错配 |
图 3-3：网格质检看板
y+ 反推第一层高度：一个完整算例
背景。某换热器管束外流任务，选用增强壁面处理的 SST 模型，要求 y+ 约等于一。来流速度每秒三米，特征长度取管径二十五毫米，空气运动黏度取每秒一点五乘十的负五次方米平方。
操作。先算摩擦速度的估算链：雷诺数约五千，查平板经验关联式得到表面摩擦系数约零点零零五三，壁面切应力等于系数乘二分之一密度乘速度平方，摩擦速度等于切应力开根号，约每秒零点一二米。再算黏性长度：运动黏度除以摩擦速度，约零点一二五毫米。首层高度等于目标 y+ 乘黏性长度，取一点二左右，得零点一五毫米。换壁面函数路线（y+ 五十），首层放宽到六毫米上下。
结果与解读。两条路线的首层高度差四十倍——这就是为什么 y+ 目标必须在划网格之前定下来：它由湍流模型决定，而湍流模型由物理判断决定，链条倒过来走必然返工。计算完成后用 Fluent 的壁面 y+ 后处理核对实际分布，若分布越界，调整 inflation 首层高度重新生成网格，而不是假装没看见。
变式。旋转机械壁面切应力分布变化剧烈，y+ 沿叶展不均匀，按最严处校核；共轭传热问题对 y+ 敏感度更高，因为壁面传热系数直接由近壁网格插值而来。
网格无关性验证：第四关的标准动作
流程固定四步。第一步定基准：用当前网格算出关键工程量（阻力、压降、换热系数）。第二步加密：按三套网格的惯例取细化比约一点三到一点五（网格量约翻倍到三倍）。第三步比较：关键量的相对变化随加密收窄。第四步定量：用 Richardson 外推与网格收敛指数（GCI）给出离散误差带，工程惯例 GCI 低于百分之一到五视为网格无关。这套动作在第九章验证与确认一节还有一次正式登场，那里会给它配上完整算例——本节你只需要记住：验收单第一栏的"签字"签的就是这次核对。
指标不合格的修法对照
指标读数的价值在于处置：每一个超线数字都对应一种修法，先修法后重算，别带着红线硬跑。
| 指标读数 | 常见成因 | 对应修法 |
|---|---|---|
| 偏斜度零点九五以上成片 | 尺寸场与几何冲突 | 局部加密或回几何简化特征 |
| 正交质量低于零点一 | 畸形单元聚集在狭缝 | 清狭缝或换 Mosaic 拓扑 |
| 边界层层间正交质量差 | 首层过厚或增长率过大 | 调层数与增长率组合 |
| 主流区大长宽比 | 过渡尺寸跳跃太大 | 平滑尺寸场降低梯度 |
| y+ 全壁面偏低 | 首层按旧工况设定 | 按 3.3 算例公式反推重设 |
三个高频疑问。平均值达标但最大值超标能放行吗：不能一刀切——最大值超标单元若远离关注区且数量极少，可书面豁免；落在关注区（分离区、传热面）必须修，位置比数量更决定成败。负体积单元在求解前要处理吗：必须清零，负体积在任何求解器里都是数学非法。指标与结果矛盾时信谁：指标是代理变量，极端工况下代理会失灵——用第九章的网格无关性做最终裁决，指标负责日常。
关于指标还有一层进阶认知：质量指标是离散误差的代理，不是误差本身。同样偏斜度零点九的单元，放在压力梯度平缓的主流区与放在激波附近，杀伤力差一个量级——所以成熟团队的网格报告不只看全局直方图，还会把"最差百分之一个单元"的空间位置标出来做位置审查。位置比数值更有信息量：最差单元扎堆在关注区，整体再好看也要返工；散落在影响圈之外，书面豁免即可。下次审网格时，把直方图与空间分布图并排打开，你看到的会完全是另一个世界。
补一个与 9.2 节衔接的伏笔式提醒：本节的所有指标与无关性验证，解决的都是"离散误差可不可信"；而"模型误差"是另一本账，由湍流与物理模型的假设决定，网格再好也治不了模型病。两本账在验收单上分栏登记、分开签字，评审时的归因讨论才有章法——网格栏签了字而结论仍与实验偏差大，问题自动指向模型栏，这是流程给工程师的免费定位服务。
本节验收点
- 偏斜度、正交质量报告已归档，红线项全部清零或书面豁免。
- y+ 目标在划网格前已由湍流模型决定，算后核对分布无越界。
- 网格无关性验证已执行，GCI 或等效误差带写进报告。
- 网格统计（总量、边界层层数、最小尺寸）作为报告附件交付。