同一套网格、同一组边界条件,一阶迎风已经收敛,切换二阶后,压降或温度却发生变化。这种情况应先检查离散格式的影响。一阶通常更容易建立稳定解,但对局部梯度的描述可能偏平缓。判断结果能否用于工程分析,还要看目标指标对格式和网格是否敏感。
本文讨论Fluent中对流项的空间离散,重点是First Order Upwind与Second Order Upwind。瞬态计算的时间离散阶数、压力插值以及VOF界面格式,需要分别选择。
Fluent采用有限体积方法,把计算域划分为许多小单元,分别核算通过单元边界的输运量。温度等变量通常保存在单元中心,而计算对流输运时,需要知道单元面上的值。“迎风”表示根据穿过该面的流动方向,使用上游单元的信息。
一阶迎风把上游单元内的变量按常值处理,直接用该单元中心值估计面值。如果温度沿流动方向变化较快,这种估计就可能遗漏单元内部的变化。
二阶迎风在中心值之外,还利用空间梯度,即变量沿空间方向变化的快慢,向单元面作线性重构。Fluent会对梯度进行限制,以控制不合理的新极值。这里的“二阶”描述的是离散精度,并不表示结果误差会固定减半。[1]

图注:上方是相同的局部单元,下方示意变量从单元中心到单元面的估计方式。竖虚线标出目标单元面,横向为位置、纵向为变量值。原理示意,无实算数据。
热量、动量或组分随流动传递时,离散近似可能额外产生类似扩散的作用,称为数值扩散。它不同于介质真实的导热、黏性或分子扩散,不能直接理解为湍流混合增强。
例如,一股携带示踪物的窄流束斜穿网格。数值扩散较强时,计算出的示踪物分布可能过度展宽,中心值下降,边缘变化趋缓。射流速度梯度、尾迹和冷热气流交界,同样可能受到影响。[2]

图注:两排使用相同入口和网格作概念对比,灰度表示示踪物相对浓度。展宽程度仅用于说明数值扩散,不代表某个案例的一阶、二阶计算结果。
在传热分析中,温度梯度改变还可能影响壁面换热和热点判断。但不能据此断言一阶一定低估最高温度,或二阶一定使温度升高。最终变化还与流动结构、热源位置、边界条件和网格分辨率有关。
数值扩散与网格分辨率、流动相对网格的方向有关。对于简单、沿网格方向发展的层流,一阶与二阶的结果可能接近。流动斜穿单元、局部转向明显,或对流占主导时,则更需要关注格式影响。[2][3]

图注:两侧只比较输运方向与网格线的关系,流束线条不代表速度大小,也不展示计算误差。沿网格方向并不等于没有离散误差。
切换二阶不能补足缺失的空间分辨率。如果一个窄射流只有很少单元跨过其宽度,或近壁温度梯度落在过粗的网格上,仍需调整局部网格。一阶、二阶在同一网格上结果接近,也只能说明这次格式变化的影响较小,不能代替网格无关性验证。
对较难收敛的模型,可以先用一阶建立初始流场,再切换到适用的高阶格式继续计算。若模型本身容易稳定求解,也可以直接使用合适的二阶设置。先一阶、后二阶并非每个模型都必须经过的步骤。[3][4]
切换前保存case和data,记录压降、流量、温度等监测量。为判断格式本身的影响,应保持网格、物性、边界条件和物理模型一致。切换后继续迭代,等待监测量重新稳定,不能把刚切换时的波动作为最终结果。
在Solution Methods中,应分别核对动量、能量和相关湍流输运方程的设置。Pressure下的Second Order是压力插值选项,不能把它与Second Order Upwind当作同一个设置;瞬态时间离散也需要单独检查。[1][3]
如果二阶下持续振荡,需要检查网格质量、梯度分辨率、求解控制,以及工况能否合理地用稳态描述。必要时有针对性地调整求解控制并重新验证。一阶解中的数值耗散可能减弱波动,仅靠退回一阶获得平稳残差,仍不足以解释振荡的来源。
内流问题可比较压降、流量分配及局部速度剖面;传热问题可比较最高温度、平均温度、热点位置和壁面换热。含射流或尾迹时,还应在相同位置提取截面,观察分布宽度与梯度,而非只比较整场平均量。
云图必须采用相同截面和色标。分别自动缩放色标,可能把很小的差异放大,也可能掩盖实际变化。单点极值容易受局部网格影响,应结合平均量与分布判断;涉及流体进出口温度时,还要统一平均方式。
同时检查质量和能量收支。数值扩散改变局部输运,并不表示热量凭空消失;全局守恒良好,也不能证明局部温度分布已经足够准确。
没有适用于所有项目的统一格式差异合格线。若两种格式造成的指标变化足以改变设计判断,就应继续检查网格、格式和收敛程度,再确定可采用的结果。
1.明确当前调整的是空间对流格式,区分压力插值和时间离散。
2.格式对比保持工况一致,保存切换前后的计算文件。
3.切换后继续计算,确认目标监测量重新稳定。
4.用同一位置、色标和统计方式比较工程指标与局部分布。
5.结合局部网格检查和守恒核对,记录最终采用的设置。
需要学习从网格、求解设置到结果核对的完整流程,可结合相应案例教程。已有模型也可说明工况和异常现象,沟通模型排错与技术辅导。
[1] Ansys Fluent Theory Guide:Spatial Discretization
[2] Ansys Fluent User’s Guide:Mesh Topologies,Numerical Diffusion
[3] Ansys Fluent User’s Guide:Choosing the Spatial Discretization Scheme
[4] Ansys Innovation Courses:Turbulent Pipe Flow,Numerical Solution