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

[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版本算例