MPPICInterFoam中的bug？

[李东岳]

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

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

[linhan.ge]

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

[李东岳]

现在手头没有DPM算例，要不你测试一下

[李东岳]

简单测试了一下，alphac.oldTime() = alphac不能解决，留给其他感兴趣的朋友们，谁有兴趣修复一下。

但是DPMFoam中没有这个问题。

[溯光]

@linhan-ge 请问这是哪篇文献中的控制方程

[linhan.ge]

@溯光 Jing, L., Kwok, C.Y., Leung, Y.F., Sobral, Y.D., 2016. Extended CFD-DEM for free-surface flow with multi-size granules. International Journal for Numerical and Analytical Methods in Geomechanics 40, 62-79.