求助： SonicFoam求解二维叶栅时无法收敛（求解文件可下载）

[宝丁]

计算虽然可以收敛，但是出口总压却一直不能固定，而且进口马赫数一旦达到0.9，整个计算就会崩掉，这难道是求解器的问题？

[宝丁]

@宝丁 出口静压不能固定

[宝丁]

这应该与我的出口边界条件有关

        /*type           fixedValue;
        //meanValue      $pOut;
        value          uniform 1e5;*/

        type            waveTransmissive;
        field           p;
        psi             thermo:psi;
        gamma           1.4;
        fieldInf        120000;
        lInf            0.2;
        value           uniform 110000;

原先想用fixedMean这个边界条件，但是一用出口边界的压力值是固定了，但是出口附近的压力会异常的高，好像出口把流道给堵住了一样；
但是用 waveTransmissive这个边界条件，边界的压力值就会一直下降，进口的压力也下降，速度升高，Ma数达到0.9左右整个计算就会崩掉。

[宝丁]

目前来看这个问题是不是已经要修改ＯＦ的一些底层代码
不知道这个难度如何，应该要放弃这个算例了

[李东岳]

哪个求解器？后来神么错误

[宝丁]

@东岳 还是用的sonicFoam， 就是出口背压一种下降，虽然可以收敛一会，但是紧接着进口Ma 数超音以后整个计算就崩掉了

[宝丁]

感觉就差最后一点，但是就是做不出来

[宝丁]

@东岳 老师，您能否给些建议，我现在应该从哪方面去想办法，现在我觉得已经试了很多，但是依然找不到方向。

[李东岳]

最近事情实在是太多了非常抱歉，找我的所有合作全都是调节算例的问题。你这个进行到目前，应该也不是我下载下来算例调节几个条件就能回发给你的，按照目前的进度，5月下旬可能会抽出一点时间给你看看你的算例。
问题挺有意思，不过非常抱歉目前时间太紧事情太多

[宝丁]

@东岳 好的，谢谢老师的关注，原先您能给的指导我已经很感动了，希望您的工作顺利

[宝丁]

这个链接里的算例是我现在迭代的最新的算例设置了，我再整理一下问题，希望浏览的人如果遇到过类似的情况的人能给些建议，我只想知道我错在哪里了。
算例可以稳定的收敛一段时间，但是出口背压不能固定，出口边界是 waveTransmissive，出口背压不断减小，进口静压也是减小，进口总压一定，进口速度也会不断增大，当进口Ma数大于1的时候整个计算就会崩掉。
如果强行换成fixedMean或者fixedValue出口背压就会不断增大，出口就像被堵住了一样。

链接：https://pan.baidu.com/s/1hzOJcuYQw0GOYWyNK_TqOw
提取码：sr9c

[Foamer24]

@宝丁 您好！请问下出错的提示是什么？

[宝丁]

@Foamer24 温度超限或者参数错误

[Foamer24]

@宝丁 温度可以在system/fvOptions文件中设置最小值和最大值。参数错误就不清楚了。

[宝丁]

@Foamer24 这个做法只是辅助性的，应该不会影响本质的算法不稳定

[Foamer24]

@宝丁 最近也在调试非定常计算，感觉openfoam里面算非定常也要像fluent一样才行，先算一个定常的初始场，然后再转非定常。不过不知道这个操作如何在openfoam里面实现。

[宝丁]

@Foamer24 有一个mapField，可以试一下，另一种也是我现在使用直接修改计算收敛值最后一个时刻的边界条件（比如总压），作为初始流场计算

[Foamer24]

@宝丁 mapfield相当于向用套粗糙的网格算一下，然后插值到密网格做初始场

[宝丁]

@Foamer24 哈哈，既然粗网格可以插值，那相同的网格不是更没有问题吗

[Foamer24]

@宝丁 细网格主要问题是存在边界层，要求时间步长较小。算的确实慢，而且很容易发散。有初始的化会好一些