SIMPLE 算法中密度如何求解?



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


  • 管理员

    如果是不可压缩,为什么要求解密度?



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



  • 0_1492182483521_草图.png

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



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



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

    只能祝好运:big_mouth:



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


登录后回复
 

与 CFD 中国 的连接断开,我们正在尝试重连,请耐心等待