圆柱计算中无关性验证的问题

crispy

@coolhhh 是的是用的二维计算的，我现在换了O型网格，在时间步长为1e-4的情况下做了网格无关性验证。使用新的网格仍然会出现减小时间步长后结果并没有随着时间步长的减小而收敛的问题。

crispy

@李东岳好的老师，直接使用回复框的这个文件上传吗？麻烦您了！

李东岳

是的

1）算例可公开 2）网格数少 3）确保可以跑通

crispy

crispy

@李东岳
网格我是存在grid里的，运行Allclean，Allrun后，会自动复制已生成好的网格到constant中，麻烦老师有空帮我看一下算例，万分感谢！

李东岳

这个算例网格太大跑不动，你处理到3万个网格我试一下

1）算例可公开 2）网格数少 3）确保可以跑通

crispy

@李东岳老师，这个是我做过网格无关性验证的网格，因为我的雷诺数约为1.5e+4，所以第一层网格高度比较小，导致网格量可能稍微有点大。处理到3万会不会网格质量不好了。

李东岳

采用非均一网格的方式，你把圆柱壁面那一片区域加密就行，其他的不用那么密，你这个是瞬态还需要跑好久

类似这种网格

crispy

@李东岳 cylinder_case.7z 好的老师，麻烦您有空看一下这个算例，谢谢！

李东岳

你发给我这个算例，默认跑了一下，目前时间不长1e-4，跑到0.89秒的时候，Cd=0.89431，是正常的么

或者说你这个需要跑到多少秒看结果

crispy

@李东岳老师我在自己这里又跑了一下，我在0.89s的时候Cd是1.089155

李东岳

你这个需要跑到多少秒看结果。然后你的Cd是看瞬时的某个时间步的结果还是要做个时间平均

crispy

@李东岳老师我这个跑到5s，算的是4-5s的时间平均结果。

李东岳

我目前跑到0.29秒，用了1873 s...

目前时间不长1e-4，跑到0.89秒的时候，Cd=0.89431，

0.089秒，说错了

crispy

@李东岳老师我现在跑了30min左右，算到1.45s了

crispy

@李东岳老师我这边算完了，最后1s平均的Cd是1.204409e+000。

李东岳

你缩小时间步10倍，再试一下。不行我这debug不了这么长时间。

crispy

@李东岳好的老师，我算完了把结果告诉您，我这边算1e-5比较慢。

crispy

老师，麻烦您看下这个结果呢

李东岳

这图做的好！

你把1e-5最后一个时间步的U以及p压缩一下发上来。我直接跑。

有关时间不长对结果的影响，之前我见过这方面讨论，但是没有仔细研究也没有好算例去弄。我一直感觉这玩意就像是理论与实际的区别。类似隐性格式照样会发散一样。你这个算例看起来是一个很好的验证算力 https://cfd-china.com/topic/746

CFD中文网