openFoam与Fluent模拟简单的Laval蒸汽喷管

[Yolo]

@李东岳 东岳老师，这个用rhoCentralFoam的case可以发上来参考一下吗

[wangxiaowei]

@李东岳 东岳老师，这个用rhoCentralFoam的case可以发上来参考一下吗

[李东岳]

http://dyfluid.com/openfoam.html#cdv

这个里面有个类似的

[wangxiaowei]

好的，感谢李老师解惑

[wangxiaowei]

@李东岳 李老师，我参考您给网址里面的案例设置了边界，发现出口压力边界用fixedValue，就会发散，而用waveTransmissive边界，计算出的出口处压力远远大于我预期设定数值，而且整个的流场感觉也不对。请问李老师，您还记得您之前用rhoCentralFoam实现的case中具体用的是哪种类型的边界条件？

[李东岳]

你的流动是超音速、跨音速、还是低音速

[wangxiaowei]

@李东岳 李老师，出口是超音速的。

[李东岳]

出口超音速边界条件不能用fixedValue

[wangxiaowei]

是的，算不了一会就发散了。那还是用waveTransmissive边界吗？感觉用这个边界算出来的结果不太对。

[shanyeyun]

@李东岳 李老师，您的这个算例，不知道 有没有关注壁面温度沿流向的变化？
我最近在用openfoam 7 ， 计算一个收敛扩张喷管的算例。壁面给 的绝热壁面 ，流体从亚音速至出口 的超音速，属于under expand flow. 核心流马赫数，温度，压力及壁面压力，壁面剪切应力 都与fluent吻合特别好，只有壁面 的温度与fluent差别太大。
请问李老师 ，帮忙解决下哈

[李东岳]

@wangxiaowei 对，waveTransmissive

@shanyeyun 没有关注过。你那个openfoam跟fluent的压力看起来是完全重合了？但是T有区别？

[shanyeyun]

@李东岳 是的，核心流的温度，压力，马赫数都重合了。壁面压力和shear stress也重合了，只有壁面温度没重合

[李东岳]

@shanyeyun 你这个如果能把算例处理成allrun的形式，其中包含一个gnuplot的脚本一键allrun后能直接生成T的那个图，如果不涉密的话，可以发上来我看看。

如果不能发上来，那就不知道了。很奇怪。流体变量都是耦合到一起的，竟然压力能温和的这么好，然后T区别那么大。

[shanyeyun]

我在怀疑，turbulence->计算alphaEff是，调用的this->alpha_为零，在哪里可以看到你的邮箱，我发给你看看

[李东岳]

网格大么？网格太大我也debug不了。网格小就发在这里。把算例处理干净。

[李东岳]

@李东岳

上面这个图是OpenFOAM-9算的。使用OpenFOAM-10计算的话，直接发散掉了。经对比应该是pEqn.H的方程，OpenFOAM-10里面这一段代码有变动：

const fvScalarMatrix divPhidp(fvm::div(phid, p));
    phiHbyA -= divPhidp.flux();

我把他改回去后，OpenFOAM-10可以收敛。但是结果异常。

然后测试OpenFOAM-12，经发现OpenFOAM-12与修改后的OpenFOAM-10均结果异常，如下图：

主要是出口的位置。经对比，出口边界条件的代码无差异。具体是哪一个代码块引起的差异目前还需要更深层次的debug。

[李东岳]

我在找bug的时候，在《无痛笔记》更新了一下openfoam无反射边界条件，目前来看植入的没有问题。bug应该不在这里。

[李东岳]

laval-of12.zip

OpenFOAM-12版本算例

[李东岳]

@李东岳 在 openFoam与Fluent模拟简单的Laval蒸汽喷管 中说：

上面这个图是OpenFOAM-9算的。使用OpenFOAM-10计算的话，直接发散掉了。经对比应该是pEqn.H的方程，OpenFOAM-10里面这一段代码有变动：

const fvScalarMatrix divPhidp(fvm::div(phid, p));
phiHbyA -= divPhidp.flux();
我把他改回去后，OpenFOAM-10可以收敛。但是结果异常。

针对这部分内容，今天碰到一个bug处理的帖子，确实是这部分原因导致的。在新版OpenFOAM-10中已经删掉

https://github.com/OpenFOAM/OpenFOAM-10/commit/910a10bf395da16de85ad12e88bc5cb9abb4f7e0