一个关于 Axisymmetric Subsonic Jet 的问题

[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]

好的 李老师 我学习一下