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

[液固两相流]

压力梯度力的计算公式如下图1所示：
图1
虚拟质量力的计算公式如下图2所示：
图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中找的~

[AppleKiller]

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

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

大佬，您最后是怎么处理的呢？我也意识到这个问题，如果要计算虚拟质量力的话，我发现粒子速度的全导数不易计算（我模拟的是气泡在水中的上浮）

[youhaoyu]

@AppleKiller 您好，请问您做气泡上浮使用什么求解器