MPPICInterFoam中的bug？

李东岳

在ESI的版本中，MPPICInterFoam中的pEqn，其对应双向耦合的连续相方程

fvScalarMatrix p_rghEqn
(
    fvm::laplacian(Dp, p_rgh)
   ==
   fvc::ddt(alphac) + fvc::div(alphacf*phiHbyA)
);

经测试其中的fvc::ddt(alphac)为0。 Jasak也表示：

Thus, if you try to evaluate the ddt the way it is written, it will recognise that the field does not have the oldTime() level stored and will use the current value to initialise it. Basically, fvc::ddt(magSqr(U)) will always give you zero, which is wrong.

在fvc::ddt(alphac)为0的情况下将导致忽略相分数的影响，难道是个bug？用算例简单测试，将fvc::ddt(alphac)去掉之后求解结果也没变化。

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"

李东岳

http://www.cfd-china.com/topic/1501/openfoam中动量方程扩散项的问题/2

看这个回答

winsway_zero

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

linhan.ge

@winsway_zero 是官网上的doxygen，https://www.openfoam.com/documentation/cpp-guide/html/index.html

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).

CFD中文网

MPPICInterFoam中的bug？

但这绝对是一个bug