buoyantBoussinesqPimpleFoam动量方程是否有误？

[ctt-twophase]

这是solver内的定义。
这是我的推导过程。
请教路过的朋友看看哪里出了问题，是不是rhok定义有误？链接文本

[李东岳]

问题解决了么？
对这个比较感兴趣，能否用latex代码推一下公式？...

[noringname]

Hi你好，如果你问的是分母里的rho缺失的话，我大概能回答一下，你方程的推导没问题，但是，你理解的p_rgh和openfoam里面的p_rgh不是一样的，你可以查看buoyantBoussinesqPimpleFoam算例hotroom里面p和p_rgh的维度，对比一下buoyantPimpleFoam算例hotroom里面p和p_rgh的维度，会发现，不可压流里面的p啊p_rgh啊早已除掉了rho，因为真正matter的是这一项的gradient，这一项的绝对值没有任何意义，但是在可压流里面，p是要用来算状态方程的。很多的发表的文章里面方程右侧的1/rho*gradP 都直接写成gradP，并不是他们写错了，而是此p非彼p，可以称为modified kinematic pressure了。

[ctt-twophase]

@noringname 看了下量纲，发现是除过rho的了，谢谢你的帮助。