rhoSimpleFoam计算翼型发散


  • OpenFOAM教授

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


  • OpenFOAM教授


  • OpenFOAM教授

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



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


  • OpenFOAM教授

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



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


  • OpenFOAM教授

    @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了么?怎么解压



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

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


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

    替代文字



  • 感觉湍流边界的设置可以改善下,速度是125m/s,为什么k是10?

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

    上U下p
    捕获.JPG


Log in to reply
 

CFD中文网 2016 - 2020 | 京ICP备15017992号-2