rhoSimpleFoam求解超音速横向射流问题

[OFtao]

问题描述：
我想在方管左侧给一个高马赫数的超音速流动，同时在底部的小喷嘴给一个亚音速的射流然后捕捉横向射流中的激波现象。

仿真设置：
OpenFOAM版本：OpenFoam-9
求解器：rhoSimpleFoam
湍流模型：KOmegaSST模型
热物理模型：

边界条件和初场：
左侧流速：3Ma,喷嘴入口流速：50m/s
alphat

k

nut

omega

p

U

T

运行报错

[李东岳]

你有没有试试rhoCentralFoam
rhoSimpleFoam稳态容易发散，要加很大的松弛

[OFtao]

@李东岳 还没有，因为我老师想要先看看稳态情况，所以我就选了rhoSimpleFoam求解器求解

[李东岳]

稳态求解器一般来说比瞬态要不稳定一些。你可以先试试瞬态，跑通了再跑稳态。

[OFtao]

@李东岳 好的，我试试，谢谢李老师

[OFtao]

@李东岳
李老师，我将求解器改成rhoCentralFoam之后，修改了一下controlDict文件和fvSchemes文件，运行之后还是同样的报错。

我感觉求解器好像都还没开始求解，是不是我上面配置文件中设置有问题呢。

[李东岳]

你这个算发散了。

如果你的算例可以公开，可以把网格处理小后，我给你调一下。

[OFtao]

李老师，我已经将网格尺寸调小了。链接: https://pan.baidu.com/s/15MF64vuhqHbjcxpmMOyQUQ 提取码: 7zxz

[李东岳]

百度网盘我下不下来，你把网格调小，算例清空，发到论坛附件试试

[OFtao]

supersonic_injet.rar

[李东岳]

supersonic_injet-stead.tar.xz

更新在此： http://dyfluid.com/openfoam.html

对了，我没加湍流，湍流你自己加一下

[OFtao]

@李东岳 好的，感谢李老师

[OFtao]

@李东岳 李老师，我在您上面的基础上修改了一下几何模型。

其他设置不变的情况下，仿真只跑了两步就发散了。

我试着修改了边界条件和初始条件，发现仿真最多跑了七步就发散了。
请问对于这种喷嘴口在外侧的情况，可能是哪里的设置导致仿真发散了呢。

[李东岳]

把你网格blockMesh发我我看看

[OFtao]

@李东岳 log.mesh

[李东岳]

blockMeshDict文件

[OFtao]

@李东岳 blockMeshDict

[李东岳]

我很久之前算过这种，你如果这个东西突出来，网格应该变成这样类似，并且内部要加密。你试一下。

[OFtao]

@李东岳 好的，谢谢李老师，我试一试

[wangpc]

@OFtao 请问您最后解决了吗，可以加个联系方式交流一下吗