MPPICFoam/denseParticleFoam 对于重力的考虑
-
(1)如果在求解时加入g(0 0 -9.81)则流场outlet处会出现不合适的速度同时有压差出现:
(2)改变g(0 0 0),流场速度正常,没有出现其他变化,
以下附上p,U.air文件
pinlet { type zeroGradient; } outlet { type fixedValue; value uniform 0; } walls { type symmetry; value uniform 0; } top { type symmetry; value uniform 0; } bottom { type zeroGradient; }
U.air文件
inlet { type fixedValue; value uniform (11 0 0); } outlet { type zeroGradient; } walls { type symmetry; } top { type symmetry; } bottom { type fixedValue; value uniform (0 0 0); }
但是在我没有添加kinematicCloudProperties以及相关颗粒离散格式、求解方式设置时这个边界条件是正常的,如果我在 kinematicCloudProperties 文件中打开了 gravity,在constant文件中g为(0 0 0),那么例子所受重力还存在吗(这个我还没有测试)。同时我在查阅以往的记录里面发现了相似的问题(链接如下,不过时间有点遥远),有没有相关大佬遇到过这方面问题,可以指点一下吗,谢谢?
https://www.cfd-china.com/topic/394/如果用dpmfoam求解稀相流会怎么样-误差大么 -
@李东岳 TP_of10.zip 东岳老师,这是of10能运行的文件
https://www.cfd-china.com/topic/394/如果用dpmfoam求解稀相流会怎么样-误差大么/47
和这个一样的情况,在不改变g情况下应该怎么修改呢,因为我的粒子需要考虑重力 -
@李东岳 李老师,最近在做实验,昨天又试验了一下,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())
);