SIMPLE 算法中密度如何求解？

[谷柏辰]

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

[赵一铭]

如果是不可压缩，为什么要求解密度？

[谷柏辰]

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

[CFD中文网]

是否是这个图的意思？The Finite Volume Method in Computational Fluid Dynamics An Advanced Introduction with OpenFOAM® and Matlab®，page 666。一般这种用来求解大马赫数问题，或者那种那个所谓的全速流动。

[谷柏辰]

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

[李东岳]

@赵一铭 不可压缩一般指的是密度针对压力的不可压缩，但是可能会因为温度的变化而改变，比如HVAC中的浮力驱动流。 @谷柏辰 很遗憾至今你的问题还没有解决，不过错位网格+matlab编程+数值不准的问题很难判断到底是哪里错了。 @cfd-china 那个书中的方法应该是错位网格所用的方法。

只能祝好运:big_mouth: ...

[谷柏辰]

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