一个关于 Axisymmetric Subsonic Jet 的问题

[yuky]

@King_RVM 方便问您一些具体的网格导入和修改边界条件的细节吗，因为我现在导入网格改边界条件总是失败

[King_RVM]

@yuky 可以参考这个文章https://www.bilibili.com/read/cv6292980/

[chenyu2016]

李老师,我最近也在算这个subsonic 这个算例 我发现,我把网格换成更加细的网格之后,结果差别很大,而且都不能收敛到实验解上, 我的openfoam算例就是您cfd中文网算的算例.

[李东岳]

@chenyu2016 https://cfd-china.com/topic/8018/免费帮助debug算例

可以的话，上传我看看

[chenyu2016]

非常感谢东岳老师, 算例我已经上传, 这两个后处理文件ipynb和py效果是一样的.
https://www.jianguoyun.com/p/Dfi7EvkQiLOUCBjtv4cGIAA (有效期至：2025-09-01 00:00)

[李东岳]

我明天给你看一下

[chenyu2016]

好的 李老师, 这个适用openfoam 2406 跑的

[李东岳]

我初步算了一下，不同网格下的nut的区别很大。这个应该是kOmegaSST湍流模型的原因。需要更详细的debug：OpenFOAM中植入的kOmegaSST与NASA的kOmegaSST的区别 你这个large网格算的比较慢，我今天要出去，你先从这个方向去debug，我持续跟进

https://turbmodels.larc.nasa.gov/openfoam_issues.html

[chenyu2016]

东岳老师 我发现了tolerance 从1e-9变回 1e-6之后这个问题就会就会好很多 , 但是变回去流场就会波动.

[李东岳]

我明天继续看一下

[李东岳]

找了半天可算找到了。没想到1e-6变成1e-9的差距这么大。2年过去了。这个是我给你们改的么，我自己都忘了..

东岳老师 我发现了tolerance 从1e-9变回 1e-6之后这个问题就会就会好很多 , 但是变回去流场就会波动.

有关这两个问题。我发现是omega的初始残差特别小。初始残差差不多在1e-7。如果设置为1e-6，则omega不会求解。最简单的办法就是设置minIter 1; 这样就可以解决了

[chenyu2016]

好的 非常感谢东岳老师, 我马上回去试一下.

[chenyu2016]

东岳 老师 , minIter 1 流场确实不波动,但是又回到了之前的问题. 不同的网格结果会差别很大.

[李东岳]

我初步给你调试了一下middle网格的算例，算的时间太长了，我尝试算完一次看看

ExecutionTime = 2264.0334 s  ClockTime = 2266 s

volFieldValue volFieldValue1 write:
    max() of U = (208.65074 3.361844e-09 59.389456)

Courant Number mean: 0.00018795879 max: 0.49110625
deltaT = 8.5609683e-07
Time = 0.00457517s

[李东岳]

初步结果出来了。湍流粘度还是不够大，导致速度场很尖锐。你的问题复现。我个人认为，继续debug的方向是湍流模型。但由于一个算例要算1小时，我这面无力debug了。

[chenyu2016]

好的 非常感谢 李老师. 我自己自己学学, 要是有结果的话我第一时间分享到网站上.

[李东岳]

我上传一个稳态的算例，middle网格的，稳态的算例依然要算1700秒。感兴趣的可以看看。

Jet_rhoSimpleFoam_3.tar.xz

[chenyu2016]

好的 李老师 我学习一下