在 CFD模拟实践中，存在一类不易直观察觉但能显著影响预测精度的误差源——数值伪扩散。该误差不同于残差振荡或网格发散等易于识别的不稳定现象，其作用方式更为隐蔽：通过对输运变量（速度、温度、浓度及湍流参数）的空间分布实施非物理平滑，致使流动界面、大梯度区域及典型流动结构（如涡旋、剪切层）的原始形态趋于模糊。深入认识并有效抑制此类误差，是准确评估CFD模拟置信度、合理选用离散格式及规划网格拓扑结构的基本前提。
1. 物理扩散的数学物理内涵
从连续介质力学视角审视，物理扩散表征了某一守恒量在分子或湍流涡团作用下由高势区向低势区的净迁移过程。典型实例包括：固体壁面附近的热量经分子热运动传导至低温区；污染物质量因浓度梯度驱动而展宽分布范围；流体微团间动量经黏性剪切力作用而逐层传递。上述过程均受确定的物理规律支配，并可通过物性参数（如热导率、质量扩散系数、动力黏度）在输运方程中定量刻画。
以通用标量输运方程为例，其描述某一标量场 \phi 的演化时，通常涵盖瞬态项（时间变化）、对流项（流体运动携带）与扩散项（分子或湍流混合导致的均化效应）。其中扩散通量与标量梯度呈正比，其强度完全取决于实际材料的物理属性，而与计算采用的网格离散方式或数值算法无关，是真实物理过程的客观映射。
2. 数值伪扩散的定义与数学成因
数值伪扩散特指在偏微分方程离散化求解过程中，由数值近似操作（包括对流项离散、变量插值及梯度重构）人为引入的、非物理的附加平滑效应。该效应并非流体真实黏性、导热或质量传递能力的体现，而是离散误差在特定形式下呈现出的拟扩散行为，其强度随数值格式与网格配置而变化。
在对流输运占主导地位的问题中，该现象尤为突出。例如，一高浓度污染物团块沿主流方向迁移，若物理分子扩散微弱，其浓度锋面理应保持陡峭形态。然而，采用一阶迎风格式等耗散型格式计算时，锋面附近的浓度过渡区会被人为展宽，形成类似于高扩散系数作用的结果。实际上，这种展宽完全源于数值截断误差中的二阶导数项，其表现形式恰与物理扩散项同构。
从数学上可这样理解：对对流项的一阶迎风离散，其截断误差的首项包含一个数值黏性系数，该系数与网格尺度、流速及流动方向与网格线的夹角有关。这个人为附加的黏性项即构成了伪扩散的数学根源。
3. 数值伪扩散的主要成因
伪扩散的强度受多重数值因素耦合影响，主要归纳为以下三类：
（1）离散格式的内禀耗散性：低阶插值格式（尤以一阶迎风格式为代表）为保证有界性及迭代稳定性，在重构界面通量时严重依赖上游节点信息，其截断误差中包含显著的二阶导数项，等效于在原输运方程上叠加了额外的扩散项。格式精度越低，此项耗散效应越强。
（2）流动方向与网格拓扑的失配：当主流方向与网格线平行时，对流信息传递的数值误差相对较小。但当流动倾斜穿越正交网格时，变量需依赖多个方向的插值来估算界面通量，导致误差沿非主流方向积累并在横向产生非物理“抹平”效应。此即多维问题中伪扩散效应显著加剧的核心原因，常被称为“交叉扩散”。
（3）空间解析度的不足：在存在局部强梯度区域（如边界层、激波、混合层）时，若网格尺度过于粗放，则流场变化无法被足够多的单元解析。数值格式在此情形下被迫在少数网格间距内完成变量的剧烈过渡，其结果等效于人为降低了梯度幅值，增大了平滑作用。
4. 数值伪扩散与物理扩散的本质差异
两类扩散虽在输运方程中表现形式相似，但其内在属性存在根本不同，对比总结表所示。
数值伪扩散与物理扩散的核心特征对比
| 特征维度 | 物理扩散 | 数值伪扩散 |
| 物理根源 | 来源于分子热运动、浓度梯度或速度梯度等真实物理机制 | 来源于数值离散格式、网格布局、插值方式等非物理的数值处理 |
| 参数依赖性 | 由材料物性参数（黏度、热导率、扩散系数）及流动状态决定，不随网格改变 | 由CFD设置（格式阶次、网格密度、时间步长）决定，随数值参数变化而显著改变 |
| 方向特性 | 通常沿物理梯度方向（如温度梯度、浓度梯度），各向同性或具有明确的物理各向异性 | 具有显著的“网格依赖性”，沿网格线方向或离散模板方向表现更突出 |
| 对结果解释的影响 | 代表真实的混合、传热或动量传递过程，是物理规律的体现 | 代表计算误差，可能高估混合效率、低估峰值，导致错误的物理结论 |
5. 数值伪扩散的有效控制策略
针对伪扩散的成因，通常可从离散格式与网格策略两个维度进行抑制：
优选高阶离散格式：对于对流占优问题，应避免将一阶迎风格式作为最终解的输出格式。该格式宜用作迭代初期的稳定化处理。推荐采用二阶迎风格式、QUICK格式、TVD（格式或MUSCL格式等高分辨率格式，它们在保持解的有界性的同时能大幅降低数值耗散。
自适应加密关键区域网格：对于存在强剪切、剧烈热交换、浓度锋面或涡核区的位置，应布置足够细密的网格。
优化网格拓扑以匹配主流方向：在工程实践中虽难以完全实现，但对于具有明确主导流向的问题（如管道流、喷管流、边界层流动），采用结构化网格或贴体坐标使网格线尽可能顺沿流向，可有效抑制横向伪扩散。
开展网格无关性与格式敏感性研究：通过在相同工况下系统加密网格或更换离散格式，若关键指标（如温度峰值、阻力系数、换热系数、回流区长度）产生显著漂移，则说明原结果受伪扩散污染严重，需改进计算策略。
6. 对数值伪扩散的客观评价
需强调的是，数值伪扩散虽本质上是一种误差，但在工程实践中并非“百害而无一利”。适度的数值耗散有时可作为稳定化机制，抑制高雷诺数或强间断问题中非物理振荡。特别是在激波捕捉等场景中，完全无耗散的格式往往难以收敛。因此，CFD 工作者的目标并非“彻底消除”伪扩散（这在数学上也难以实现），而是确保其量级控制在工程允许的误差范围内，保证关键物理结论（如涡结构演化、热斑位置、分离点）由真实的物理模型主导，而非数值副作用。
数值伪扩散是CFD离散误差的一种特殊表现形式，在对流占优、网格较粗、流向倾斜及采用低阶格式时尤为显著。其与真实物理扩散的区别在于：后者源于物性参数，描述客观物理规律；前者源于数值处理，反映计算方法的局限性。一个收敛残差小、分布云图光滑的计算结果，并不能自动等同于真实物理状态。只有经过严格的网格无关性检验、格式敏感性分析与物理合理性论证的数值解，方能被认为是可靠的计算结论。
