overpimpleDyMFoam +颗粒边界处理

lys

@星星星星晴你说的sourceTerm是颗粒对流体的作用力吗

李东岳

会有部分颗粒进入固体内部

具体没太细看啊，我插一句，是不是需要看一下你这个颗粒定义在哪套网格上了，如果定义在A网格上，壁面在B网格上，那是不是肯定会发生这个现象。

星星星星晴

@lys 比如ueq里面的 + parcels.SU(U) 这一项，其他的有加么

// Solve the Momentum equation

MRF.correctBoundaryVelocity(U);

tmp<fvVectorMatrix> tUEqn
(
    fvm::ddt(rho, U) + fvm::div(phi, U)
  + MRF.DDt(rho, U)
  + turbulence->divDevTau(U)
 ==
    rho()*g
  + parcels.SU(U)
  + fvOptions(rho, U)
);
fvVectorMatrix& UEqn = tUEqn.ref();

UEqn.relax();

fvOptions.constrain(UEqn);

if (pimple.momentumPredictor())
{
    solve(UEqn == -fvc::grad(p));

    fvOptions.correct(U);
    K = 0.5*magSqr(U);
}

星星星星晴

@李东岳正常情况下 cloud里面会有updateCellOccupancy这个方程，会自动更新当前cell中parcel的占用什么的，如果在solver中把parcels.evolve(); 放到mesh.changing之后就应该没事，但是很奇怪，他的parcel会进入空的地方。。

lys

@李东岳这是重叠网格，按道理讲是可以更新背景网格的流场信息，然后颗粒在更新后的流场信息作用下，是可以实现壁面效应捕捉的。但也有可能，颗粒只是感受到了流场，而没有实际发生壁面碰撞，所以导致颗粒穿透壁面，这是不是意味着overset 方法是不适用颗粒求解的，所以openFOAM-ESI版本下就没有现成的overset求解器+LPT

李东岳

你的A网格如果没有壁面，但是壁面定义在网格B上，那拉格朗日栗子同时也定义在A上，那A网格那面不就没有壁面了么。流场的话，重叠网格做了插值处理。栗子，是不是也要做B网格壁面-A网格壁面的映射，然后并反弹处理？

lys

@星星星星晴这项时肯定要加的

星星星星晴

@lys 嗯，就是确认一下别拉下什么东西没加就难受了。

所以说，那些particle 应该是在 $n - 1$ 步时还在网格里，然后网格发生变化，可能就不在网格里了，这部分parcel就没办法算了。

把parcel.evolve() 放到 bool changed = mesh.update() 前面，让parcel用上一步的旧网格试试看
另外如果你的Euler time step $Δ t_{E}$ 足够小，就是你每次更新网格的间隔减小，会不会也一定程度上减少这样的发生呢？
调整动网格或者在parcel部分，判断内部是否有parcel，有parcel就挪他们到wall表面。

lys

@星星星星晴

星星星星晴

@lys ？怎么有解决么？

dxl

@lys 你可以看下你的solver中有没有添加storeGlobalPositions这个函数，我对你使用的这个solver不是很了解，拉格朗日是你自己添加进去的吗

CFD中文网

overpimpleDyMFoam +颗粒边界处理