4. 网格质量诊断与优化

网格质量直接决定数值解的稳定性、收敛速度与物理可信度。在 Fluent 中，网格质量差通常不会立刻报错，而是表现为残差震荡、连续性方程难以收敛、局部出现非物理极值，甚至把本应光滑的梯度“抹平”。因此，网格质量诊断不是建模完成后的可选项，而是与边界条件、求解设置同等重要的前置工作。本节围绕三个核心指标——偏斜度（Skewness）、正交质量（Orthogonal Quality）、长宽比（Aspect Ratio）——说明如何检查、如何判断、如何优化，并进一步讨论网格无关性验证与自适应网格细化的工程做法。

一、三个核心质量指标的含义与判据

偏斜度衡量单元偏离理想形状的程度。Fluent 中偏斜度取值 0～1，0 表示等边/正方形单元，1 表示完全退化。工程上通常要求：最大偏斜度低于 0.85 可接受，低于 0.80 较好，超过 0.95 往往导致求解困难。偏斜度对对流项离散尤其敏感，因为高偏斜单元会使迎风格式中的面法向与单元中心连线严重偏离，引入数值扩散或非物理振荡。

正交质量衡量单元面法向与“面中心—相邻单元中心”连线的对齐程度，取值 0～1，越接近 1 越好。它与偏斜度互补：偏斜度关注单元自身形状，正交质量关注相邻单元之间的几何关系。一般要求最小正交质量大于 0.1，低于 0.01 时求解器可能直接拒绝计算。对于边界层网格，正交质量比偏斜度更能反映近壁面处的数值行为。

长宽比是单元最长边与最短边之比。在边界层内，长宽比达到 100～1000 是正常的，因为需要沿壁面法向加密以捕捉速度梯度；但在主流区，长宽比过大（例如超过 100）会加剧各向异性误差，使扩散项离散失真。因此长宽比不能孤立判断，必须结合单元所处区域：边界层允许高长宽比，远场与回流区应尽量接近 1～10。

二、在 Fluent 中执行质量检查

在 Fluent Meshing 或求解器中，可通过 Mesh → Check 或 TUI 命令获取质量统计。常用命令如下：

/mesh/check
/mesh/quality
/mesh/skewness
/mesh/orthogonal-quality
/mesh/aspect-ratio

执行后应重点关注三组数据：最小值、最大值和平均值。例如某散热器模型报告：最大偏斜度 0.91、最小正交质量 0.03、最大长宽比 320。此时不能只看平均值，因为求解稳定性由最差单元决定。应进一步用 Mesh → Display → Mesh 配合 Skewness 着色，把偏斜度大于 0.9 的单元高亮出来，通常集中在几何尖角、小圆角、薄壁交界或边界层与主流过渡区。

三、针对性优化策略

若最差单元出现在几何尖角处，优先做几何清理：把小于网格尺寸的圆角、倒角、小孔简化或删除，避免网格被迫在极小空间内生成退化单元。若出现在边界层与主流过渡区，可降低边界层增长率（Growth Rate 从 1.2 降到 1.1），或增加过渡层数，使单元尺寸平缓变化。若出现在回流区或尾流区，可在该区域添加局部尺寸控制（Body of Influence 或 Face Sizing），把目标尺寸缩小到原来的 1/2～1/3。

对于长宽比问题，若主流区出现长宽比大于 100 的单元，通常说明网格方向与流动方向严重错位。此时应改用扫掠网格（Sweep）或六面体主导（Hexcore），让单元长边沿流动方向排列，而不是横跨流动方向。若几何过于复杂无法扫掠，可考虑多区域划分，把规则区域单独切出做六面体，复杂区域保留四面体/多面体。

四、网格无关性验证

网格无关性验证的目标是证明计算结果不再随网格加密而显著变化。工程上通常做三套网格：粗、中、细，单元数比例约为 1:2:4 或 1:1.5:2.25。以某管道压降计算为例，三套网格单元数分别为 48 万、96 万、210 万，监控进出口压降 Δp：粗网格 1250 Pa，中网格 1180 Pa，细网格 1165 Pa。中网格到细网格变化约 1.3%，小于 2%，可认为中网格已达到无关性；后续优化或参数扫描可沿用中网格，以节省计算成本。

验证时需注意：加密方式应保持一致，例如统一按全局尺寸比例缩小，而不是只加密某一区域；监控量应选择对工程目标最敏感的量，如压降、升力、最高温度、流量分配比，而不是只看残差。若相邻网格结果变化超过 5%，说明网格仍偏粗，需继续加密或检查是否存在局部质量瓶颈。

五、自适应网格细化

自适应网格细化（Adaptive Mesh Refinement, AMR）根据计算得到的物理量梯度自动加密网格。Fluent 中可通过 Adapt → Gradient 或 Adapt → Region 实现。典型流程是：先完成一次初步收敛计算，再根据速度梯度、温度梯度或压力梯度标记需要加密的单元，执行 Adapt → Refine，然后继续迭代。例如在激波捕捉中，可按压力梯度阈值 5000 Pa/m 标记单元，加密后激波位置的分辨率明显提高，而远场网格保持不变，总单元数仅增加约 15%。

AMR 的代价是网格质量可能下降，因为加密后的过渡区容易出现悬挂节点和尺寸突变。因此每次自适应后都应重新检查偏斜度与正交质量，必要时配合光顺。对于瞬态问题，AMR 可动态跟随界面或涡结构，但需控制加密频率，避免网格数无限增长。

六、工程检查清单

把上述检查固化为流程，可以显著减少“算得动但算不对”的情况。网格质量诊断与无关性验证不是一次性动作，而是伴随几何修改、边界条件调整和物理模型切换反复执行的工程习惯。