overpimpleDyMFoam +颗粒边界处理

lys

@星星星星晴
两个版本都用了，of2206的overpimpleFoam的基础上，我给加了LPT，但是就是我觉得不理想

星星星星晴

@lys 哦哦，可能是你的其他模型不够吧，受力都有啥？particle turbulence dispersion model？1way 还是2way这些东西是有影响的，particle maxCo是多少？

lys

@星星星星晴受力主要阻力，重力，流体惯性力（压力梯度力），虚拟质量力，应该比较全了，2-way的，湍流脉动力确实没有考虑，颗粒Co=0.3，感觉跟这个关系不大，主要是边界没起到时时刻刻的不可穿透条件

星星星星晴

@lys 你可以试试Co小一点，这个也有影响的，判断你的拉格朗日dt. 正常情况下是需要考虑dispersion的，不然particle的RMS是错的（个人经验）。不过of自带的dispersion model 也不咋地，不过你可以试试看。也要看你用什么turbulence。 RANS 还是LES，RANS需要，LES怎么说呢，理论上也是是需要的，但是好多人不用。OF也没有默认的LES的turbulence dispersion。俩模型都$k-\varepsilon$的。你也没说你用的啥turbulence

我做过和你类似的，不过不是圆柱，不是你的重叠网格，是tablet周围被snappyhexmesh加密过，完全没有你说的情况，particle就是碰到壁面就没了，个人认为，是因为你的重叠网格，也不知道这个重叠网格怎么做的，另外如果你particle很小，不是dense的话，有些模型不用考虑，完全是浪费时间。shear-lift？不好意思，有些词我想不起来中文教材怎么叫了。。。

星星星星晴

另外请问你的source term在各个方程里面加了么？

@lys 在 overpimpleDyMFoam +颗粒边界处理中说：

@星星星星晴
两个版本都用了，of2206的overpimpleFoam的基础上，我给加了LPT，但是就是我觉得不理想

lys

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

李东岳

@lys 在 overpimpleDyMFoam +颗粒边界处理中说：

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

具体没太细看啊，我插一句，是不是需要看一下你这个颗粒定义在哪套网格上了，如果定义在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 $\Delta t_E$ 足够小，就是你每次更新网格的间隔减小，会不会也一定程度上减少这样的发生呢？
调整动网格或者在parcel部分，判断内部是否有parcel，有parcel就挪他们到wall表面。

lys

@星星星星晴

星星星星晴

@lys ？怎么有解决么？

dxl

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

CFD中文网

overpimpleDyMFoam +颗粒边界处理