rhoSimpleFoam计算翼型发散
-
理论上rhoSimpleFoam可以计算可压缩流动问题,官方给的tutorial中也有计算超音速流动的算例。但是我在计算一个三段翼的亚音速流动(U=125m/s)时,却无论如何都不能收敛,但是我把速度降到60之后就很容易收敛。不知道各位大佬有何高见?附件是计算的网格和设置文件,各位有兴趣可以试一试。rhoSimple.zip
-
-
@dzw05 抱歉网站上传附件有大小和格式限制。所以我把压缩文件分成三份。rhoSimple.z01和rhoSimple.z01文件,下载后请去掉zip后缀。rhoSimple.z02.zip
-
问题已经解决?下载这俩个文件可以?
-
@东岳 并没有解决,还发现了新问题。simpleC和simple算的升力系数大概差20%
-
非气动专业,尝试和大佬讨论,想问一下:
(1)div(phid,p) 的格式为什么不用 bounded Gauss upwind 呢?经验来看,bounded格式稳定很多;
(2)grad用最小二乘,以及interpolation用迎风也是为了数值稳定性吧?但貌似当前算例都用中心差分也可以,所以为了精度给改了。
按照以上两项修改试着跑了几万步,不过不是特别收敛,连续性方程不满足,升阻力系数还在震荡。是不是这个物理问题本身是非定常的,或者是网格太粗糙呢?
还有就是挺好奇这个网格使用sHM画的吗?像这类混合网格如何处理trailing edge这种特变尖锐的地方? -
@Jacobian 我也是初用OpenFOAM算气动问题。
(1)对于bounded格式,实际上看代码的话,会减去phid对应的场量的散度再乘以p。如果散度为0,那么没有问题。但是phid的散度实际上不为0,所以我倾向于不用bounded。
(2)最小二乘其实是为了保证精度。前面的cellLimited是为了数值稳定性。
我用Simple-C算没问题,但Simple一直发散,目前还没找到原因。
(3)网格其实实用Pointwise做的,SHM不太容易做。 -
下载rhoSimple.z02.zip rhoSimple.z01.zip 打不开呢
-
@dzw05 学习了,div(phid,p) 确实不用 bounded 为好
@东岳 把 z01~02 后面的zip去掉,然后对着 rhoSimple.zip 选择解压 -
zip去掉之后不久成了rhoSimple.z02和rhoSimple.z01了么?怎么解压
-
气动方向菜鸡选手,也想和大佬讨论。
- 大佬要不要试一下OpenFOAM里的SA模型?这个一方程模型我看他们在算气动时还挺常用的,可能会好收敛一些,不过kwsst精度好像会好一些。
- 网格可能也有影响。边界条件或许可以参考OpenFOAM里算翼型的case。
- 还有可能是这个翼型不一定适合这么高的流速,看到的一些文献里马赫数在0.2左右,所以收敛可能就是挺难的。
- 对了还想问楼主SIMPLEC和SIMPLE算法有什么区别呀,看您的讨论似乎会好收敛一些?
-
可算解压缩出来了。。。先看网格,大佬这网格画的很有意思
-
https://www.jianguoyun.com/p/DecMC2kQ9s3ZBhiw1P0EIAA 我粗糙的改了一下,目前可以收敛,不过总感觉网格好像有点稀,网格渐变很粗暴,
上U下p