MPPICInterFoam中的bug？

linhan.ge

@李东岳有另外一个疑问，在MPPICInterFoam的UEqn中的第7行，为什么没有alphac项呢，

，
即是否应为：

fvm::Sp(fvc::ddt(alphacRho))+fvc::div(alphaRhoPhic)

linhan.ge

@linhan-ge 抱歉上面的公式有误，应为：fvm::Sp(fvc::ddt(alphacRho)+fvc::div(alphaRhoPhic),U)?

李东岳

UEqn中的第7行目前已经在湍流模型里面包含进去了

linhan.ge

@东岳能否请大神详细解释下，您所指的包含在湍流模型，是否是指这里的line99~102

还想请教下，为什么MPPICinterfoam里的湍流模型要用phaseCompressibleTurbulenceModel这个模板类呢？

李东岳

是的，你看他已经包含进去了

#include "PhaseCompressibleTurbulenceModel.H"

李东岳

看这个回答

winsway_zero

@linhan-ge 您好，您这个使用的IDE吗？

linhan.ge

winsway_zero

@linhan-ge 哦哦，明白了

linhan.ge

@东岳还是没有理解，包括第8行，不应该也有流体项的体积分数吗？即 fvm::div(alphacRhoPhi,U),湍流模型里面看起来是没问题的,包含了体积分数我参考的方程是这样的：

李东岳

rhoPhi中已经在createFields.H中包含了alpha，您看看

linhan.ge

@东岳这里的rho是除颗粒相外，流体的体积分数吧？您说的是这里吗？

难道是这里的phi里包含了？

李东岳

再往上看，这个rho是混合密度。rho=alpha1 * rho1+alpha2 * rho2

linhan.ge

ρf = α1ρ1 + α2ρ2 前面应该还有个 αf = α1 + α2 ？

李东岳

这个下标f应该是fluid的意思。这一行代码：

rho = alpha1*rho1 + alpha2*rho2

乍一看应该是

rho = alpha1*rho1

这样就跟公式符合了。

目前只能跟进到这里，你再看看？

李东岳

一楼中提到的bug，fvc::ddt(alpha)的值取决于时间格式，将其换为backward就可以了

更新下：问题在于alphac.oldTime()。

在Euler中，计算的是

(alphac - alphac.oldTime())/delta T

Euler不起作用的原因在于alphac.oldTime()=0。

backward是可以解决问题。但是在backward中调用的是三层时间步。即vf、vf.oldTime()和vf.oldTime().oldTime()。起作用的是vf.oldTime().oldTime()。alphac.oldTime()还依然是0。

但是，问题在于alphac，其他的场oldTime()正确。

李东岳

在我自己的求解器中，我验证了双向耦合中界面交换相、相分数的耦合作用大小。在我测试的算例中，相分数的耦合作用较小。耦合主要体现在彼此之间的力交互。

但这绝对是一个bug

linhan.ge

@东岳当颗粒位于界面处时，MPPICInterFoam是不是没考虑颗粒对界面的影响？因为MPPICInterFOAM中的相方程还是跟interFOAM中一样的吧？MULES那块不怎么看得懂。

linhan.ge

@东岳这里有个问题，mppicInterFOAM中，alpha1的表达式是什么？是占整个流体的体积分数alpha1=V1/( V1+V2)。还是占三相的？ alpha1 = V1/(V1+V2+Vparticle).

linhan.ge

@东岳这个bug能否通过手动赋值解决，alphac.oldTime() = alphac ?