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

[crispy]

请教各位一个问题，我最近在计算圆柱绕流，雷诺数约为1.5e+4，使用的求解器是rhoPimpleFoam，使用的湍流模型是kOmegaSST，但是我换一个网格，换一个时间步长，都会得到差异不小的平均阻力系数，请问这是什么原因呢？
先附目前算得最好的网格


时间步长目前计算了1e-4和5e-6
1e-4的平均阻力系数1.197349，和实验结果1.2较为接近，St的误差在2%。
5e-6的平均阻力系数1.724599。
不是说时间步长越小，结果越精确吗？为什么会出现这样的问题呢？刚开始用openfoam，很多都不懂，谢谢大家了。

[coolhhh]

这是算二维的网格吗？之前算过三维方柱绕流，用ICEM的oblock，感觉这种网格质量很差。不知道是否可能网格质量引起的

[李东岳]

你可以把你的算例打包上传我看看

[crispy]

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

[crispy]

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

[李东岳]

是的

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

[crispy]

@李东岳 cylinder_case.7z 谢谢老师！

[crispy]

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

[李东岳]

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

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

[crispy]

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

[李东岳]

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

类似这种网格

http://dyfluid.com/openfoam.html#id5

[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不了这么长时间。