MPPICFoam/denseParticleFoam 对于重力的考虑

sxz0823

@李东岳东岳老师，您好，我还有个考虑，我的是气固流动颗粒在70um，可不可以将连续项里面的重力加速度项去掉 phicForces里面的 rAUcf*(g & mesh.Sf())

这样会影响到我粒子的重力计算吗，此时我的g仍然给（0 0 -9.81）

李东岳

不能去掉。去掉就不对了。

sxz0823

@李东岳那东岳老师，您在这篇帖子中提到的办法是指什么呢，1方法换成p_rgh，“直接将将pEqn.H改为p_rghEq.H，是无密度速度方程+prgh方程，虽然流场预测正确，但这种植入在密度不为1的时候耦合项不正确导致速度过低。（引用了您在这篇帖子的结论）”；所以我想着有没有办法去掉动量方程中的重力影响，但是保留颗粒相的重力（在不可压缩的前提下）
https://www.cfd-china.com/topic/394/如果用dpmfoam求解稀相流会怎么样-误差大么/59 ，

李东岳

这个已经过去太久了，我已经忘了怎么回事了.. 等我处理个bug之后再给你看看吧

sxz0823

@李东岳好的，谢谢东岳老师

sxz0823

@李东岳东岳老师，我把rAUcf*(g & mesh.Sf())注释掉以后，除了压力场为0，似乎没有造成其它影响，但是肯定不符合正常要求了，那还是要在MPPIC里面添加p_rghEqn.H吗？东岳老师可以简单指导一下吗

李东岳

我现在给你看一下你发的算例

李东岳

你这个发那个算例，我跑起来了，但是我这跑不动，有50万网格。能不能处理到5万。

sxz0823

@李东岳 tp_of10.zip 好了，东岳老师

李东岳

这个尺度正确么，x方向1公里，高度390米

sxz0823

@李东岳正确的

李东岳

http://dyfluid.com/openfoam.html#id2

跟这个算例一样，你把p的出口边界改成fixedMean就好了

sxz0823

@李东岳好的我试试，谢谢李老师，MPPICFoam也是一样的吗，不需要改成p_rgh吗

李东岳

不需要

sxz0823

@李东岳李老师，没有用过这个边界条件，总是给不对格式

李东岳

30     outlet
 31     {
 32         type            fixedMean;
 33         meanValue       uniform 0;
 34         value           uniform 0;

请忽略代码行号数字

sxz0823

@李东岳非常感谢李老师帮助

sxz0823

@李东岳李老师，最近在做实验，昨天又试验了一下，p场稳定了一些但是还在波动，但是U.air场变得混乱了，李老师有空可帮我看下吗，还有一个问题719e060c-5757-402d-b40a-0d978c2be0a7-tp_of10.zip

DPMFoam中的UcEqn.H:

fvVectorMatrix UcEqn
(
fvm::ddt(alphac, Uc) + fvm::div(alphaPhic, Uc)

fvm::Sp(fvc::ddt(alphac) + fvc::div(alphaPhic), Uc)

continuousPhaseTurbulence->divDevRhoReff(Uc)
==
(1.0/rhoc)*cloudSU
);

UcEqn.relax();

volScalarField rAUc(1.0/UcEqn.A());
surfaceScalarField rAUcf("Dp", fvc::interpolate(rAUc));

surfaceScalarField phicForces
(
fvc::flux(rAUccloudVolSUSu/rhoc) + rAUcf(g & mesh.Sf())
);

if (pimple.momentumPredictor())
{
solve
(
UcEqn
==
fvc::reconstruct
(
phicForces/rAUcf - fvc::snGrad(p)*mesh.magSf()
)
);
}

///////还有就是这一项代表了什么意思呢////////
surfaceScalarField phicForces
(
fvc::flux(rAUccloudVolSUSu/rhoc) + rAUcf(g & mesh.Sf())
);

疏影横斜水清浅

@sxz0823 请问，

///////还有就是这一项代表了什么意思呢////////
surfaceScalarField phicForces
(
    fvc::flux(rAUccloudVolSUSu/rhoc) + rAUcf(g & mesh.Sf())
);

这个有答案了吗

dxl

@sxz0823 请问我这把重力设置为（0 0 0 ）之后计算会报错是什么原因呢

CFD中文网

MPPICFoam/denseParticleFoam 对于重力的考虑