rhoSimpleFoam计算翼型发散

[dzw05]

理论上rhoSimpleFoam可以计算可压缩流动问题，官方给的tutorial中也有计算超音速流动的算例。但是我在计算一个三段翼的亚音速流动（U=125m/s）时，却无论如何都不能收敛，但是我把速度降到60之后就很容易收敛。不知道各位大佬有何高见？附件是计算的网格和设置文件，各位有兴趣可以试一试。rhoSimple.zip

[dzw05]

@dzw05 rhoSimple.z01.zip

[dzw05]

@dzw05 抱歉网站上传附件有大小和格式限制。所以我把压缩文件分成三份。rhoSimple.z01和rhoSimple.z01文件，下载后请去掉zip后缀。rhoSimple.z02.zip

[李东岳]

问题已经解决？下载这俩个文件可以？

[dzw05]

@东岳 并没有解决，还发现了新问题。simpleC和simple算的升力系数大概差20%

[Jacobian]

非气动专业，尝试和大佬讨论，想问一下：
（1）div(phid,p) 的格式为什么不用 bounded Gauss upwind 呢？经验来看，bounded格式稳定很多；
（2）grad用最小二乘，以及interpolation用迎风也是为了数值稳定性吧？但貌似当前算例都用中心差分也可以，所以为了精度给改了。
按照以上两项修改试着跑了几万步，不过不是特别收敛，连续性方程不满足，升阻力系数还在震荡。是不是这个物理问题本身是非定常的，或者是网格太粗糙呢？
还有就是挺好奇这个网格使用sHM画的吗？像这类混合网格如何处理trailing edge这种特变尖锐的地方？

[dzw05]

@Jacobian 我也是初用OpenFOAM算气动问题。
（1）对于bounded格式，实际上看代码的话，会减去phid对应的场量的散度再乘以p。如果散度为0，那么没有问题。但是phid的散度实际上不为0，所以我倾向于不用bounded。
（2）最小二乘其实是为了保证精度。前面的cellLimited是为了数值稳定性。
我用Simple-C算没问题，但Simple一直发散，目前还没找到原因。
（3）网格其实实用Pointwise做的，SHM不太容易做。

[李东岳]

下载rhoSimple.z02.zip rhoSimple.z01.zip 打不开呢

[Jacobian]

@dzw05 学习了，div(phid,p) 确实不用 bounded 为好
@东岳 把 z01~02 后面的zip去掉，然后对着 rhoSimple.zip 选择解压

[李东岳]

zip去掉之后不久成了rhoSimple.z02和rhoSimple.z01了么？怎么解压

[shampoo]

气动方向菜鸡选手，也想和大佬讨论。

1. 大佬要不要试一下OpenFOAM里的SA模型？这个一方程模型我看他们在算气动时还挺常用的，可能会好收敛一些，不过kwsst精度好像会好一些。
2. 网格可能也有影响。边界条件或许可以参考OpenFOAM里算翼型的case。
3. 还有可能是这个翼型不一定适合这么高的流速，看到的一些文献里马赫数在0.2左右，所以收敛可能就是挺难的。
4. 对了还想问楼主SIMPLEC和SIMPLE算法有什么区别呀，看您的讨论似乎会好收敛一些？
   

[李东岳]

可算解压缩出来了。。。先看网格，大佬这网格画的很有意思

[李东岳]

https://www.jianguoyun.com/p/DecMC2kQ9s3ZBhiw1P0EIAA 我粗糙的改了一下，目前可以收敛，不过总感觉网格好像有点稀，网格渐变很粗暴，

上U下p

[AppleKiller]

@李东岳 这些算例应该也可以使用rhoCentralFoam求解器计算吧？

[李东岳]

@AppleKiller 忘了多少音速了，如果是跨音速是可以的。

[AppleKiller]

@李东岳 老师好，rhoCentralFoam对于超音速的情况不能算吗？

[李东岳]

@AppleKiller 不可压缩不可以

[三吉]

@李东岳 在 rhoSimpleFoam计算翼型发散 中说：

https://www.jianguoyun.com/p/DecMC2kQ9s3ZBhiw1P0EIAA 我粗糙的改了一下，目前可以收敛，不过总感觉网格好像有点稀，网格渐变很粗暴，

上U下p

您好，为什么我打开这个网页下载数据结果是喷管的结果，不是翼型的

[李东岳]

http://dyfluid.com/openfoam.html#id7 看看这个对不对