CFD应用科普｜离散误差估计

离散误差的估计通常依赖于细化网格解的对比分析，具体方法可参见第 3.11.1.2 节与第 3.9 节。如前所述，此类误差根源于方程及边界条件中引入的近似处理。针对具有光滑解的问题，近似质量可通过“阶数”这一指标来表征，该指标揭示了截断误差与网格间距幂次之间的内在联系：若空间导数的截断误差正比于 ，则称该近似为  阶。

需要强调的是，阶数并非误差绝对值的直接度量，而是反映了误差随网格间距变化的渐近规律。在特定网格上，相同阶数的近似方案其误差量级可能相差显著；甚至在某些情况下，低阶近似的误差反而小于高阶近似。然而，随着网格间距的不断减小，高阶近似终将展现出更优越的精度特性。

泰勒级数展开法是确定近似解阶数的有效工具。然而，当方程中不同项采用各异的近似方法时，整体求解阶数往往难以直观判断，通常受限于关键项中精度最低的近似阶数。此外，代码实现中的潜在缺陷也可能导致实际收敛阶偏离理论预期。因此，通过数值实验验证各类问题的实际求解阶数显得尤为重要。

在结构化网格中，分析离散化误差的理想策略是将各方向间距减半，但这一方法在实际应用中常受限制。尤其在三维情形下，网格间距减半会导致节点数量增至八倍；若进行第三次细化，节点总数将达到初始网格的64倍，计算成本往往令人望而却步。另一方面，鉴于误差分布通常具有非均匀性，全局网格细化并无必要。对于采用任意控制体积或单元的无结构网格而言，由于缺乏明确的局部坐标方向，其单元细化机制也与结构化网格存在显著差异。

网格细化的核心在于确保其具备实质性与系统性。例如，仅将某一方向的节点数从 112 增至 128，除非针对误差分布均匀且网格本身均匀的纯学术问题，否则在实际应用中意义有限。有效的网格细化应遵循严格标准：各方向节点数至少比原始网格增加 50%，或网格间距缩小至原来的 1/1.5。同时，系统性细化要求在所有网格层级上保持拓扑结构及节点相对空间密度的一致性。值得注意的是，即便节点总数相同，节点分布的差异也可能导致离散化误差发生显著变化。如图 8.10 所示，在近壁面加密的非均匀网格上计算所得  结果，其精度比同节点数的均匀网格高出一个数量级。尽管两种方案均以二阶收敛于同一网格无关解，但在节点数相同的前提下，误差量级差异仍可达 10 倍甚至更多。图 6.4、6.6 及 8.17 中的案例也印证了这一结论。

上述分析凸显了优化网格设计的重要性。在实际工程中，网格生成往往最为耗时，获得合格网格已非易事，高质量网格更是难求。优质网格应尽可能接近正交，但需注意“正交性”在不同数值方法中定义各异；在 FV 方法中，单元面法线与相邻单元中心连线的夹角最为关键，四面体网格在此意义上亦可视为正交。更为重要的是，在预期截断误差较大的区域应适当加密网格，这要求设计者预先掌握解的部分特征信息（详见第 12.1.1 节）。采用局部细化的非结构网格最能满足这一准则。此外，网格平滑度、纵横比及扩展比等其他质量标准，则需根据具体数值方法确定。关于网格质量度量的更多细节，将在第 12.3 节中详述。

估计离散化误差最简便的方法基于理查森外推法（Richardson Extrapolation）。该方法的核心前提是：在足够精细的网格上进行计算时，结果应呈现单调收敛特性；若无法满足这一条件，实际误差往往超出预期。因此，只有当两个最精细网格均达到足够精度，且误差收敛阶数已知时，该方法才能保证准确性。进一步地，若三个网格均满足精细度要求，则可依据下式计算收敛阶（详见 Roache 1994；Ferziger 和 Peri ́c 1996）：

式中， 为网格加密系数（例如当网格间距减半时，取 ）， 表示平均间距为  的网格上的数值解。相应的离散化误差估计方法如下（详见第 3.9 节）：

因此，当网格间距减半时，对于二阶方法，当前网格解的误差约为该网格与前一级网格解之差的三分之一；而对于一阶方法，其误差则直接等于该差值。

以图 8.10 为例，将理查森外推法分别应用于均匀网格与非均匀网格时，尽管两者的解误差相差一个数量级，但所得网格无关解的估计值在五位有效数字内仍保持高度一致。

需要指出的是，误差评估既可针对阻力、升力等积分量进行，也可基于物理场值展开，但这些物理量的收敛阶次往往存在差异。对于层流等具有光滑解的问题，实际收敛阶次通常与理论阶次（如二阶）相符；然而，当涉及湍流、燃烧、两相流等复杂模型，或采用包含开关、限制器的数值格式时，准确确定收敛阶次往往较为困难。在此情形下，精确计算阶次 p 或 Roache (1994) 提出的 网格收敛指标 （grid convergence index）并非绝对必要。通常仅需展示计算量在一系列网格（建议三个）上的变化趋势即可：若变化呈现单调性且差异随网格细化而减小，便可直观估算出网格无关解的大致位置。当然，在条件允许的情况下，仍应优先采用理查森外推法来获取更为精确的网格无关解估计。

此外，网格细化并非必须在整个计算域内均匀实施。当误差评估表明特定区域的误差显著低于其他区域时，采用局部细化策略往往更为合理。这一方法在绕体流动模拟中尤为适用，因为高分辨率通常仅需集中在物体表面及尾流区域。相较于全域或整块网格的均匀细化，局部细化策略在计算效率上更具优势。对于非结构化网格而言，即便是最粗层级的网格，往往也已内置了局部细化机制。然而，应用该策略时需格外审慎：若未能对主要误差源（即截断误差较大的区域）进行有效细化，受对流和扩散等传输机制的影响，未细化区域的误差可能向周边传播，从而导致整体细化效果大打折扣。

在评估离散化误差时，有两类难点需特别关注。首先是湍流计算中的壁函数问题（详见 10.3.5.5 节）。在此类方法中，壁面并未施加唯一的边界条件，而是通过建立壁面剪切应力与紧邻壁面网格中心流速之间的关联，并假定该位置处于边界层对数律适用范围内。然而，当流动受复杂壁面几何形状影响，导致法向与切向存在显著梯度时，对数律假设往往不再严格成立。随着网格细化，近壁网格中心的位置随之移动，这实质上改变了动量方程的边界条件。此类变动常引起积分量的异常波动，且其变化规律不符合所用离散化格式的预期收敛行为。当网格细化致使近壁网格中心落入缓冲层（）时，该问题尤为突出。若必须采用壁函数，建议采取以下策略之一：确保所有网格的近壁计算点均位于对数层内（）；或者在所有网格序列中保持近壁棱柱层第一层厚度固定，使近壁单元的  维持不变，仅在切向方向对壁面附近网格进行细化。

当网格分辨率足以解析边界层时（即采用“低-Re”方法，近壁单元中心满足 ），动量方程可直接施加无滑移条件，从而确保边界条件的唯一性。在此情形下，针对力、力矩及热传递等积分量的变化，通常更适合采用理查森外推法进行误差估计。需要特别强调的是，应用该方法时必须严格保证所有近壁单元的 。若  落入中间过渡区域，虽然“全 y+”壁函数的表现优于“高 Re”壁函数，但最佳实践策略仍是确保在使用壁函数时，计算域内绝大部分壁面的  维持在 30 以上。

网格对几何特征（尤其是壁面曲率）的解析能力不足，是数值模拟中另一类常见问题。涡轮机械叶片前缘便是典型示例，涵盖船用螺旋桨、风扇、水轮机、燃气轮机以及风力或潮汐涡轮机叶片等部件。此类区域通常具有极大的曲率，且压力极值点间距极短，因此要求在壁面的法向与切向上均布置精细网格。然而，商业软件的自动网格工具若未在生成过程中被特别配置（例如指定前缘网格间距，或要求特定节点以解析曲率），往往难以充分捕捉前缘的几何细节。这可能导致在粗网格条件下，原本圆滑的前缘呈现出尖锐状。此时，网格细化过程实质上改变了求解域的几何形状，致使部分物理量的变化偏离预期规律，进而增加了应用理查森外推法的难度。因此，确保初始网格能够精确表征求解域的几何形状至关重要。图 12.7 展示了一个相关示例。

注：本文内容取自《Computational Methods for Fluid Dynamics》 12.2.2节。

”

（完）