简介：在CFD仿真中，网格质量直接决定计算能否收敛以及结果是否可信。本文围绕正交质量、偏斜度、长宽比等核心指标展开，给出可落地的判断标准、常见问题表现与优化策略，帮助工程师在求解前识别风险网格、在求解中定位发散原因，从而减少无效迭代与结果偏差。
在进行计算流体力学（CFD）仿真时，一个常被低估的事实是：求解器的稳定性与结果精度，很大程度上在网格划分阶段就已经被决定了。很多工程人员习惯把不收敛、残差震荡、流场异常等问题首先归因于湍流模型选择不当或边界条件设置错误，但实际排查中，网格质量往往才是更深层的原因。
从数值求解的角度看，CFD本质上是将偏微分方程离散到有限体积单元上进行迭代求解。网格单元的形状、大小、过渡方式以及相邻单元之间的关系，都会直接影响离散方程的系数矩阵特性。当网格质量较差时，系数矩阵可能出现病态，导致迭代误差无法有效衰减，表现为残差曲线长期高位震荡或突然发散。更隐蔽的情况是计算最终“收敛”了，但收敛到的是一个由网格缺陷引入的伪解，与实验值或物理直觉严重不符。
因此，建立一套系统的网格质量判断与优化流程，对提升仿真效率和可信度具有直接工程价值。
工程中常用的网格质量评价指标有多种，但最需要重点关注的通常是以下三类：正交质量、偏斜度和长宽比。它们从不同侧面描述网格单元与“理想单元”的偏离程度，各自对求解过程的影响机制也不完全相同。
正交质量描述的是网格单元面法向量与相邻单元中心连线方向之间的夹角关系。理想情况下，对于六面体或四边形单元，相邻单元中心之间的连线应当垂直于共享面，此时通量计算最准确，离散误差最小。
当正交质量下降时，共享面上的通量插值会产生较大误差，尤其是对于扩散项和压力项。实践表明，对于大多数稳态RANS求解，正交质量最低值应不低于0.1～0.15，平均值应高于0.7。如果最低正交质量低于0.05，求解器在迭代中即使使用较保守的松弛因子，也可能难以建立稳定的压力-速度耦合关系。
判断方法上，可以通过网格统计报告查看正交质量的最小值和分布。如果最小值过低但占比极小，可以尝试局部修复；如果大量单元集中在中低区间，则需要重新评估网格划分策略。
偏斜度衡量的是实际网格单元与理想规则单元的偏差程度。对三角形或四面体而言，偏斜度反映其形状偏离正三角形或正四面体的程度；对四边形或六面体，则反映其偏离矩形或正六面体的程度。偏斜度越高，单元越“扭曲”。
高偏斜度单元对求解的影响主要体现在两个方面：其一，它会使梯度重构精度下降，因为从相邻单元中心到当前单元中心的几何关系变得不再规则；其二，它会降低局部矩阵的对角占优特性，从而削弱迭代的稳定性。一般来说，偏斜度低于0.85～0.9是可接受范围，超过0.95的单元应尽量避免出现在关键流动区域。对于边界层等对精度要求高的区域，建议将偏斜度控制在0.8以下。
长宽比是网格单元最长边与最短边的比值。在边界层网格中，为了用可接受的网格数量分辨近壁面流动，必然存在较大的长宽比，这本身并非错误。但长宽比过高且过渡不均时，会导致相邻单元的尺寸差异过大，使离散方程中各系数的量级差异显著，从而影响迭代收敛速度。
对于边界层区域，长宽比达到数百甚至上千是常见做法，但关键是保证过渡平滑。对于主流区域，建议长宽比控制在10～50以内。判断是否合理的一个重要原则是：流动参数梯度大的方向上，单元尺寸变化应尽可能缓，而不是单纯追求某一个指标绝对小。
网格质量问题并不总会在求解前被网格统计报告完全暴露，很多情况下它是在迭代过程中通过异常现象表现出来的。熟悉这些表现与网格质量之间的关联，有助于快速定位问题。
残差曲线无法下降或持续震荡：如果初始残差下降非常有限，或者下降到某一量级后就不再降低并在一个区间内反复震荡，同时伴随局部物理量出现非物理波动，应优先检查网格偏斜度和正交质量在问题区域的分布。特别是当某区域的偏斜度超过0.95且单元数量不少时，局部迭代误差很难被有效抑制。
特定位置出现明显异常值：当压力场、速度场或温度场在某个局部区域出现尖刺状或棋盘状分布，而这些区域恰好存在质量较差的网格单元时，通常可以判断为离散误差被局部放大。调整网格后，这些异常往往会明显改善甚至消失。
不同网格数量下结果差异巨大：如果网格加密后求解结果发生显著变化，且加密前后均未达到网格无关性，说明粗网格中可能存在较严重的质量缺陷，导致离散解尚未进入渐近收敛区间。此时应优先改善网格质量而非继续盲目增加数量。
面对一个不收敛或结果异常的仿真案例，建议按照以下顺序进行排查，而不是直接调整求解参数。
在将网格导入求解器之前，应完成一次面向质量的全面检查。除了查看统计报告中的最小值、最大值和平均值外，更应关注质量指标的空间分布。一个有效的方法是：将偏斜度或正交质量以云图方式显示，观察差质量单元是否集中在关键区域，如边界层、流动分离区、曲率剧烈变化处或几何狭窄通道中。
对于发现的局部差质量区域，应判断其是否处于高敏感流动区。如果位于远离关键区域的远场边界附近，其影响可能有限；如果位于壁面附近、激波区域或分离区内，则必须优先修复。
如果求解过程中出现不收敛，可以先降低松弛因子以稳定迭代。但更重要的是在较低的Courant数或较保守的松弛条件下，观察残差异常的位置。部分求解器支持输出残差的空间分布或场变量的残差信息，这能帮助定位问题单元。另一种实用方法是，在疑似问题区域局部加密或重新划分网格后重新求解，如果残差行为明显改善，即可确认网格质量是主要矛盾。
即使计算收敛，也应进行网格无关性验证。建议至少准备三套网格，且在加密时注意不只是增加单元数量，还要关注质量指标的同步改善。若加密后结果变化仍超出工程允许范围，说明当前网格方案尚未达到可接受的分辨率与质量水平。
改善网格质量通常需要从划分策略和局部修复两个层面同时入手。
整体策略层面，优先采用结构化或混合网格方案。对于形状规则的计算域，结构化网格能提供更好的正交性和可控的长宽比。对于复杂几何，可以采用近壁面结构化或棱柱层、核心区非结构化的混合方式，兼顾边界层精度和几何适应性。边界层网格的层数、第一层高度和增长率应基于目标y+值提前估算，避免在求解后反复返工。
局部修复层面，对于几何尖角、狭缝、曲率突变等区域，应通过几何简化和局部加密来降低网格生成难度。例如，对不参与流动的关键小特征进行合理简化，或对尖角区域设置更小的局部尺寸。若网格生成后仍存在少量差质量单元，可以使用专门的网格修复功能进行局部光顺或节点移动，但需注意修复范围不宜过大，以免引入新的过渡问题。
参数设置层面，在网格质量无法完全理想的情况下，适当采用更稳健的离散格式。例如，在迭代初期使用一阶迎风格式以获得稳定的初始流场，收敛到一定程度后再切换到二阶格式提高精度。但需要明确，这种策略只能缓解求解稳定性问题，无法从根本上弥补网格缺陷带来的精度损失，最终仍应以合格的网格质量为前提。
综合实践经验，可以将网格质量的工程判断准则归纳为以下几条可执行规则：
第一，偏斜度最大不超过0.95，关键流动区域控制在0.85以下；第二，正交质量最小值不低于0.1，对于精度敏感的边界层区域应尽量高于0.3；第三，长宽比在主流区控制在50以内，边界层内允许较高但过渡必须平滑；第四，网格无关性验证中，相邻两套网格结果的关键物理量变化应小于工程允许偏差。
需要强调的是，这些数值并非刚性阈值，而是基于大多数工程问题的经验参考。实际判断时，应结合具体问题的物理特性、求解目标以及可接受误差范围灵活调整。但无论具体数值如何变化，核心原则是不变的：在关键流动区域保持足够好的网格质量，比在非关键区域追求均匀的网格分布更有价值。
随着工程问题复杂程度不断提升，网格质量对仿真结果的影响只会更加显著。建立一套从划分前规划、划分中控制、求解前检查到求解后验证的完整闭环，是每个CFD工程师提升工作效率和结果可信度的必经之路。