openFOAM中颗粒压力梯度力于虚拟质量力的计算问题

液固两相流

压力梯度力的计算公式如下图1所示：
图片1.png 图1
虚拟质量力的计算公式如下图2所示：
图片2.png 图2
在openFOAM中的计算代码如下所示：
压力梯度力：
value.Su() = mass*td.rhoc()/p.rho()*DUcDt;
虚拟质量里：
value.Su() = Cvm*PressureGradient
OpenFOAM里面的计算公式与实际的计算公式有一定的差别，请教一下大家这种处理是基于什么考虑呢？

李东岳

我这面不怎么关注压力梯度力，给你分享个文献 https://www.jianguoyun.com/p/DSAdyIkQ9s3ZBhjMz6AD 你看看，希望你看完之后能分享一下你的理解，谢谢

星星星星晴

@液固两相流在 openFOAM中颗粒压力梯度力于虚拟质量力的计算问题中说：

ass*td.rhoc()/p.rho()*DUcDt

如果说压力梯度力这块你觉得没有乘up的话，在particleforces中的forceSuSp.H中有解释

 Helper container for force Su and Sp terms.

        F = Sp(U - Up) + Su

    Explicit contribution, Su specified as a force
    Implicit coefficient, Sp specified as force/velocity

然后PressureGradient这一部分是要去cloud中找的~