[求助]天然气泄漏问题

bruceyao

我现在做的是天然气泄漏问题，具体模型就是：
6mp的天然气，经直径为10mm的小孔喷射出来，然后就模拟喷射出来的流动，及其在空气中扩散情况，现在做的是自由扩散。
问题是经6Mp的压力出来，空口处的气体速度应该很大，马赫数远大于0.3，应该作为可压缩气体进行模拟，但是在稍远离空口流动速度又降下来了，又是不可压缩流动。
我现在把它当作可压缩算的化，算到中途都不能收敛，提示发散
当作不可压缩算的话，收敛后，空口处的压力显示很低，只有几个大气压而已，跟现有的文献相比较，明显过低。然而出口处的速度倒是很大，有好几千米每秒，与现有的文献相比，明显过大。
我自己认为，是当作不可压缩处理，气体压力膨胀作功都转化为动能，因而使得气体速度很大，事实上，是应该有一部分转化为气体内能，使得气体温度上升，气体密度变化。
但是我在把它当作不可压缩算收敛后，再当可压缩算的话仍然无法收敛。
我用的边界条件是压力入口，压力出口。
还有一个问题是在扩散边界，也就是压力出口边界的回流比较明显。