首页/文章/ 详情

Fluent一阶与二阶格式,差别不只是收敛速度

50分钟前浏览0

同一套网格、同一组边界条件,一阶迎风已经收敛,切换二阶后,压降或温度却发生变化。这种情况应先检查离散格式的影响。一阶通常更容易建立稳定解,但对局部梯度的描述可能偏平缓。判断结果能否用于工程分析,还要看目标指标对格式和网格是否敏感。

本文讨论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

来源:BB学长
Fluent湍流UM控制ANSYS
著作权归作者所有,欢迎分享,未经许可,不得转载
首次发布时间:2026-10-28
最近编辑:50分钟前
BB学长
硕士 | 研发工程师 公众号BB学长 知乎BB学长
获赞 167粉丝 317文章 266课程 1
点赞
收藏
作者推荐
未登录
还没有评论
课程
培训
服务
行家
VIP会员 学习计划 福利任务
下载APP
联系我们
帮助与反馈