MPPICFoam/denseParticleFoam 对于重力的考虑

[李东岳]

这个尺度正确么，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 ）之后计算会报错是什么原因呢

[sxz0823]

@dxl 给个联系方式，我加你。这种改法是错的，不考虑重力就不现实

[dxl]

@sxz0823 好的老师，我不是做这个大气流动方向的，但这个求解器相关的向您请教请教。我微信/QQ：710347000

[李东岳]

@疏影横斜水清浅 重力以及体积力

[dxl]

此回复已被删除！