Your browser does not seem to support JavaScript. As a result, your viewing experience will be diminished, and you have been placed in read-only mode.
Please download a browser that supports JavaScript, or enable it if it's disabled (i.e. NoScript).
@李东岳 老师,这个是我做过网格无关性验证的网格,因为我的雷诺数约为1.5e+4,所以第一层网格高度比较小,导致网格量可能稍微有点大。处理到3万会不会网格质量不好了。
采用非均一网格的方式,你把圆柱壁面那一片区域加密就行,其他的不用那么密,你这个是瞬态还需要跑好久
类似这种网格
http://dyfluid.com/openfoam.html#id5
@李东岳 cylinder_case.7z 好的老师,麻烦您有空看一下这个算例,谢谢!
你发给我这个算例,默认跑了一下,目前时间不长1e-4,跑到0.89秒的时候,Cd=0.89431,是正常的么
或者说你这个需要跑到多少秒看结果
@李东岳 老师我在自己这里又跑了一下,我在0.89s的时候Cd是1.089155
你这个需要跑到多少秒看结果。然后你的Cd是看瞬时的某个时间步的结果还是要做个时间平均
@李东岳 老师我这个跑到5s,算的是4-5s的时间平均结果。
我目前跑到0.29秒,用了1873 s...
目前时间不长1e-4,跑到0.89秒的时候,Cd=0.89431,
0.089秒,说错了
@李东岳 老师我现在跑了30min左右,算到1.45s了
@李东岳 老师我这边算完了,最后1s平均的Cd是1.204409e+000。
你缩小时间步10倍,再试一下。不行我这debug不了这么长时间。
@李东岳 好的老师,我算完了把结果告诉您,我这边算1e-5比较慢。
@李东岳
老师,麻烦您看下这个结果呢
这图做的好!
你把1e-5最后一个时间步的U以及p压缩一下发上来。我直接跑。
有关时间不长对结果的影响,之前我见过这方面讨论,但是没有仔细研究也没有好算例去弄。我一直感觉这玩意就像是理论与实际的区别。类似隐性格式照样会发散一样。你这个算例看起来是一个很好的验证算力 https://cfd-china.com/topic/746
@crispy 你把fvScheme里面的Euler,换成CN格式重新算一下1e-5的结果看看。你从第4秒开始算就可以
@李东岳 cylinder_1e-5_last.7z 好的老师,不好意思我昨天没看到消息,我先把包含最后一个时间步计算的文件夹发您,您看一下是这个意思吗?
@李东岳 老师我现在正从第4s开始使用CN格式计算,系数设置的是0.9。目前还不知道会不会发散,我之前使用backward是算了一会儿就发散了。等计算完成后我再整理结果发上来。
@李东岳 老师,对于dT=1e-5,4s开始使用CN计算得到的平均阻力系数约为1.637893e+000,Eluer计算得到的平均阻力系数约为1.636890e+000。
@crispy 我给你重新写一个求解器看一下,得下周。
@李东岳 好的谢谢老师!