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

[李东岳]

是的

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

[crispy]

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

[crispy]

@李东岳

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

[李东岳]

这图做的好！

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

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

[李东岳]

@crispy 你把fvScheme里面的Euler，换成CN格式重新算一下1e-5的结果看看。你从第4秒开始算就可以

[crispy]

@李东岳 cylinder_1e-5_last.7z 好的老师，不好意思我昨天没看到消息，我先把包含最后一个时间步计算的文件夹发您，您看一下是这个意思吗？