CFD中文网

如题，除去遍历网格中心坐标的方法外，还有没有其他方法可以确定离某个坐标点的距离为d的所有网格的编号。

@李东岳 后面我测试了粒径更小的颗粒，发现同等条件下更小粒径的颗粒不会出现上述情况，应当是MPPICInterFoam要求网格体积和颗粒体积的比值比DPMFoam的更高

@李东岳 没有液相时1、2号颗粒也会在初始状态下产生向上运动的趋势

@李东岳 先是颗粒替换成了碰撞颗粒，发现颗粒运动异常，然后使用未修改的MPPICInterFoam尝试，颗粒运动还是出现类似情况，去除液相后的DPMFoam预测的基本符合预期。

我想对MPPICInterFoam做一些修改，故想用颗粒沉降实验测试一下MPPICInterFoam的求解结果，但是如图所示，颗粒1和颗粒2在其他颗粒下降时发生了向上的位移，而颗粒3在其他颗粒下降到中间位置时有一个突然的加速，网格上显示的时流体的速度。颗粒的注射方式采用的时manualInjection，颗粒的受力采用的都是默认的设置

particleForces
    {
        WenYuDrag
        {
            alphac      alphac;
        }
        gravity;
        interface
        {
            C            -10;
            alpha        alpha.water;
        }
    }

反应器上部分是空气下部分是水如下图所示， transportProperties采用的也是默认设置

phases          (water air);

water
{
    transportModel  Newtonian;
    nu              1e-06;
    rho             1000;
}

air
{
    transportModel  Newtonian;
    nu              1.48e-05;
    rho             1;
}

sigma           0.07;

单颗粒沉降时颗粒还算正常，但是颗粒一多就会出现这种情况，想问下各位同行有遇到过类似情况吗？

最近在看DPMFoam的代码，但是发现动量方程好像和其数学表达式对不上

fvVectorMatrix UcEqn
(
    fvm::ddt(alphac, Uc) + fvm::div(alphaPhic, Uc)
  - fvm::Sp(fvc::ddt(alphac) + fvc::div(alphaPhic), Uc)
  + continuousPhaseTurbulence->divDevTau(Uc)
 ==
    (1.0/rhoc)*cloudSU
);

UcEqn.relax();

fvOptions.constrain(UcEqn);

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

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

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

    fvOptions.correct(Uc);
}

此处通过fvc::reconstrct返回的

fvc::snGrad(p)*mesh.magSf()

是否需要乘上alpha