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

李东岳

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

类似这种网格

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 好的老师，不好意思我昨天没看到消息，我先把包含最后一个时间步计算的文件夹发您，您看一下是这个意思吗？

crispy

@李东岳老师我现在正从第4s开始使用CN格式计算，系数设置的是0.9。目前还不知道会不会发散，我之前使用backward是算了一会儿就发散了。等计算完成后我再整理结果发上来。

crispy

@李东岳老师，对于dT=1e-5，4s开始使用CN计算得到的平均阻力系数约为1.637893e+000，Eluer计算得到的平均阻力系数约为1.636890e+000。

李东岳

@crispy 我给你重新写一个求解器看一下，得下周。

crispy

@李东岳好的谢谢老师！

李东岳

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

你把1e-4的结果也发上来，我从1e-4的结果上开始算，从1e-5开始算要算好久

CFD中文网