首页/文章/ 详情

01 基础知识 | 网格无关性验证

7月前浏览844

    最近可能是进入大论文重要节点了,经常有人问我网格无关性验证的问题,公 众号前期都是从实用工程的角度出发,并没有对Verification and validation 这个环节进行介绍,今天就浅显讲一下网格无关性验证的问题。

1 网格无关性验证是什么

网格无关性验证也叫做网格收敛性/独立性验证,其本质是:当网格加密到一定程度后,继续细化网格,计算得到的关键物理量(如速度、压力、阻力系数、换热系数等)的变化量小于预先设定的误差阈值,此时可以认为模拟结果不再依赖网格数量,达到了网格无关解。

2 为什么要进行此类验证

众所周知: CFD 常用的方法是将连续的流场离散为有限个网格单元进行求解,网格过粗会导致离散误差过大,结果失真;网格过细则会大幅增加计算资源消耗和时间成本,且无实际必要,在实际的仿真计算中,受限于计算资源与时间,我们需要一套设置合理的网格进行仿真计算。

而且 Verification and validation 也是CFD仿真领域发表高水平学术论文的必备环节之一,如果缺少该验证,论文大概率会被相关审稿人质疑。这也是为什么大部分导师让大家进行网格无关性验证的原因。

isope

Ocean Engineering

Journal of Engineering and Applied Science

3 网格无关性验证该怎么做

3.1 网格无关性验证案例

此处以常用标模SUBOFF为例展开说明:

    1. 通过对相关文献的阅读,首先确定一套基准网格进行计算,以阻力系数、压力系数等作为物理量,以阻力作为计算方案验证的物理量(这个后续推文会展开计算方案验证的物理量)

    2. 以基准网格为Medium,分别细化Fine和粗化Coarse,网格数量增加比例可以在1.2-3之间

    3. 绘制观察物理量随网格数量的变化曲线/表格,并观察收敛状况。

    4. 若无法达到收敛,则重复步骤2,在Fine的基础上进一步按照1.2-3的增长率细化网格

    5. 重复步骤3,以下图为网格无关性验证结果,可以发现当网格持续细化,相关的物理量已经趋于收敛,因此可以选择合适的网格量进行后续的验证(个人感觉这篇文献还可以再细化一下网格

3.2 注意事项

  •  在流体计算中经常进行多工况的仿真计算,但网格无关性验证一般只需要考虑其中一种工况,该工况并不是随意选择,而是应该选择速度最快、模型最大等最应该去细化的工况展开验证。

  • 所选择的检测物理量,优先为后续计算所需的物理量,在仿真计算中,不同的湍流模型和边界条件会导致监测物理量数值存在差异,尽量选择对仿真结果影响最大的物理量进行检测。

  • 在进行网格加密的时候,可以无需对整个流体域进行加密,重点对需要精细化处理的位置进行加密,例如在外流场的计算中,可以不对外围网格进行加密,重点对计算体周边展开加密。

PS

  • 个人其实感觉这个无关性验证还有一点问题,例如

    • 在保证检测到的物理量与标模一致的情况下,可以通过虚假的网格数量提升(加密不重要区域的网格),使检测到的物理量几乎不变,从而到达网格无关性验证的目的。

    • 受限于计算资源的限制,面对复杂模型无法做到完全收敛的时候该怎么办,是否会有人采用上节提到的方案人为控制收敛。:

  • 想要深入了解Verification and validation的读者可以去看一下《Guide for the Verification and alidation of Computational Fluid Dynamics Simulations



来源:CFD仿真库
湍流UM控制
著作权归作者所有,欢迎分享,未经许可,不得转载
首次发布时间:2026-01-13
最近编辑:7月前
CFD仿真库
硕士 分享知识,互帮互助,共同进步。
获赞 21粉丝 105文章 142课程 0
点赞
收藏
作者推荐

Actran 基础案例 | 点声源计算

本案例为声学计算基础案例,为点声源在自由场中声传播的计算,主要涉及Actran软件自建模、网格划分、计算设置与后处理结果查看。可以覆盖该软件的大部分操作,建议大家仔细阅读,为大家后续学习该软件打下坚实基础。1Actran前处理设置1.1几何设置首先点击Mesh选择Ellipsoid,建立球体域,相关尺寸设置如下图所示:1.2网格划分声学网格可以通过最大分析频率和声波波长确定,最大的声学网格要小于波长的1/6,波长可以由声速确定,即:其中,c是声速。本文计算的频率范围为500Hz,声速为340,因此可以得到声学网格的尺寸,并进行网格划分首先采用meshonmesh对上节生成的2D网格进行重新绘制,具体的绘制方法如下图所示:其次对3D区域展开网格划分,选择volumemesh对3D域进行网格填充,网格尺寸选择0.06,具体划分方法如下图所示:1.3计算域划分设置首先采用Autocreatedomains对现有的计算域进行自动划分,并对各部分进行重新命名,具体设置如下图所示:1.4材料设置此处对材料进行了设置,相关设置如下图所示,流体材料选择空气,相关设置如下图所示:2Actran设置2.1分析方法设置右击Aanlysis,选择DirectFrequencyResponse,并对求解的频率范围进行设置,相关设置如下图所示:2.2组分设置此处对有限流体区域进行了设置,相关设置如下图所示,流体材料选择空气,相关设置如下图所示:其次对无限元边界进行设置(无限元边界是actran所采用的一种特殊边界条件),相关设置如下:2.3边界条件设置对点声源进行设置,该设置为边界条件设置,具体操作如下图,幅值设置为1:2.4后处理设置2.4.1监测点设置首先建立相对应的监测点,具体设置情况如下图所示:2.4.2云图设置其次对云图展开设置,具体设置如下图所示:2.5运行设置将计算文件导出,统一选择ACTRAN格式导出,具体设置如下图:打开launch程序,选择计算文件,点击RUN进行计算:3Actran后处理结果3.1PLTviewer结果首先打开PLT文件,对监测点的频谱图结果进行读取,相关结果如下图:3.2Maps结果其次对声学云图展开可视化处理,导入NFF文件,具体数据导入如下图所示,将声压结果进行导入,不同频率下的云图结果如图所示:来源:CFD仿真库

未登录
还没有评论
课程
培训
服务
行家
VIP会员 学习计划 福利任务
下载APP
联系我们
帮助与反馈