@同学博 好的 十分感谢！

@队长别开枪 好的 十分感谢！

问一下各路大神，可不可以基于OpenFOAM开发出一个有界面的应用软件出来？

@西湖冷月 认同！CFD计算收敛与否不能仅仅看残差曲线，所有变量的残差都低于设定值不一定收敛，在不发散的情况下残差达不到预设的精度要求也不一定不收敛，需要查看进出口质量、能量的守恒情况，还可检测几何内关键点的计算结果变化情况。

@李东岳已经非常感谢了，我再继续努努力吧！:expressionless:

非常非常感谢，我的算法流程图跟这个基本类似，不过密度的更新我是的理想其他状态方程，因为马赫数很低，假定为不可压缩流体，我的程序也是用matlab写的，让我百思不得其解的是我的程序都是按照标准方程写的，也能启动运气，但随着计算的进行，速度方程会首先不满足boundedness条件，求解的速度值超大，然后就是整个程序的不收敛。我现在是怀疑我的密度的处理有问题，所以才有此问，不知道该往哪个方向debug:upset: !

@赵一铭 虽然是不可压缩，但其中涉及很多反应，气体组分会大大改变，所以密度也会变，当下的处理方程是用理想气体状态方程，rho=MP/RT 但结果并不理想，是否有其他更新密度手段？

SIMPLE 算法中连续方程被改为压力修正方程来求解，请问密度改怎么办？不可压缩流体

@赵一铭 整个程序我用的SIMPLE算法，所以连续性方程改成了压力修正方程来求解压力，不能直接求解出密度，怎么办？

在压力，速度，组分，温度都求解出来之后，对气体密度该如何更新？由于马赫数较低，假定为不可压缩流体，当下更新公式用的rho=MP/RT,但总感觉这个公式在温度升高较快时，导致密度降低太快，使得求解不收敛，请问该如何处理？:big_mouth: