@液固两相流 在 openFOAM中颗粒压力梯度力于虚拟质量力的计算问题 中说:
压力梯度力的计算公式如下图1所示:
图片1.png 图1
虚拟质量力的计算公式如下图2所示:
图片2.png 图2
在openFOAM中的计算代码如下所示:
压力梯度力:
value.Su() = mass*td.rhoc()/p.rho()*DUcDt;
虚拟质量里:
value.Su() = Cvm*PressureGradient
OpenFOAM里面的计算公式与实际的计算公式有一定的差别,请教一下大家这种处理是基于什么考虑呢?
大佬,您最后是怎么处理的呢?我也意识到这个问题,如果要计算虚拟质量力的话,我发现粒子速度的全导数不易计算(我模拟的是气泡在水中的上浮)