一个关于 Axisymmetric Subsonic Jet 的问题

[cccrrryyy]

@King_RVM 不是这个，是PIMPLE里面的nOuterCorrectors，本来设的是多少？

[King_RVM]

@东岳 是的李老师，这个是计算出来的，公式如下：


其中k取的是0.001，这个值比较小是因为当入口湍动能值较大时，得到的速度值偏大，改小后结果有所改善（但是相比大的k值，小的k值残差振荡比较严重）。turbulent viscosity 取得是0.001，这个是参考的NASA官网提供的值。

[King_RVM]

@东岳 李老师，我做了一个rhoSimpleFoam 10w步的结果与您提供的13000的结果的速度云图对比，因为没监测Ux（轴向速度），目前只能提供一个瞬时的速度Umagnitude的对比结果，如下图：

13000这个结果有所改善，但同样在大概2倍的出口半径处仍然存在小的波动。李老师可否提供postProcessing中的监测点数据，这样我可以画个沿轴向的速度曲线来对比。下图是之前我算rhoSimpleFoam的速度曲线，也可以看出在射流出口处速度波动较大，和云图表现出的结果是一致（我们猜测会不会是射流喷嘴处壁面的影响）。

[King_RVM]

@cccrrryyy 不好意思，我刚开始学还没弄明白这些值。pimple里这个的设置如下：

SIMPLE
{
    residualControl
    {
        p               1e-4;
        U               1e-4;
        "(k|omega|e)"   1e-4;
    }

    nNonOrthogonalCorrectors 0;
    pMinFactor      0.1;
    pMaxFactor      2;
}

PIMPLE
{
    nCorrectors              2;
    nNonOrthogonalCorrectors 1;
    nOuterCorrectors         1;
    pMinFactor      0.1;
    pMaxFactor      2;
}

relaxationFactors
{
    equations
    {
        ".*"     1;
    }
}

nOuterCorrectors默认的是1，这个算例是我根据OpenFoam自带的tutorials-rhoPimpleFoam 中的aerofoilNACA0012的算例修改的，因此这个值就没变。

[李东岳]

1. 你可以直接把13000的数据放到你的算例下面，后处理数据再出图
2. https://www.jianguoyun.com/p/DTTBLxIQ9s3ZBhjxi6cD 或者你下载我配置后的rhoSimpleFoam的算例文件 出个结果图我看看 这几天我更新在rhoSimpleFoam解析里面
3. rhoPimpleFoam那个我感觉是湍流模型的问题，算的太慢了不好调试，你试试调整下湍流边界条件，看看有无改善

[李东岳]

这个是rhoPimpleFoam算的，结果只出到0.02秒，不过看起来也行了。我更新在rhoPimpleFoam解析，然后你去下载一下玩玩

等你的plot图啊

[King_RVM]

@东岳 太感谢了李老师， 我用您的配置算例，然后监测了singleGraph一条线上的速度值，U的15000步迭代后的结果如图1和2所示，射流喷嘴出口处仍存在小的波动，但整体计算结果相比NASA的数值结果更靠近实验数据（NASA的实验数据是在不同半径处沿360度周向取了20个监测点得到的速度均值），这个在速度云图2中也可以看出来：


经过时间光滑后Umean的结果如图3和4所示，出口处速度波动会消失，只是远离喷嘴处结果和实验偏离的有些大：

[King_RVM]

@东岳 这个rhoPimpleFoam的结果也太好了，是修改了湍流模型吗？真期待您的计算结果。
还想问一下李老师，rhoSimpleFoam里面您修改了outlet的边界条件，并且增加了Ux，Uy，和Uz的边界值和初始值，主要是outlet边界条件对结果有重要影响吗？

[King_RVM]

这个是计算结果文件：https://www.jianguoyun.com/p/DQHmutkQ1IPKCBjH46cD

[李东岳]

https://www.jianguoyun.com/p/Dbf3umcQ9s3ZBhjhk6cD 这个是瞬态的文件，你对比一下吧，我改了一些内容，记不太清哪些有影响了，欢迎你的理解与分享！

感谢你的plot

[King_RVM]

@东岳 好的，我学习一下，也非常感谢东岳大哥的指导与交流。

[King_RVM]

@东岳 经过对比验证，我们发现fvSolution中绝对残差tolerance对结果的影响非常大，默认值是1e-6，调整成1e-9后结果有很大改善，改变这个值可以提高计算精度。

[李东岳]

@King_RVM 感谢分享！！

[李东岳]

@king_rvm 在 一个关于 Axisymmetric Subsonic Jet 的问题 中说：

我们发现fvSolution中绝对残差tolerance对结果的影响非常大，默认值是1e-6，调整成1e-9后结果有很大改善，改变这个值可以提高计算精度。

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

[King_RVM]

@李东岳 是的李老师，之后就是用您调的结果算的

[King_RVM]

@李东岳 对，我们也没想到影响会这么大，有时候这类原因都想不到

[yuky]

求问大佬如何绘制楔形模型的，是在icem里画完了导入openfoam的嘛，再改的wedge的边界条件嘛，我是要算一个轴对称圆柱燃烧室

[King_RVM]

@yuky 是的，用画网格软件画的，这个网格最初是来源于NASA官网

[yuky]

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

[King_RVM]

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