@人来人往 您好，请问您是石油大学华东这边的吗，最近也准备在集群上安装OF，完全没有头绪，可以请教下吗？谢谢

@东岳 好的，我去了解一下，谢谢各位老师的指点

@东岳 不好意思东岳老师，刚看到消息，我没看过相关的文章就是自己用超算时候发现的

@东岳 好的，谢谢东岳老师回复

DPM并行计算的时候，发现好像仅把求解域分块，单独求解流体的U、P，但是对于粒子相好像还是在一个process里面单独求解，并不是想象中的也按照分块计算，粒子数目多了计算变得非常慢，请问各位老师有什么好的办法吗？谢谢各位老师

@东岳 好的，谢谢东岳老师的回答，另外想问一下老师，DPM中颗粒碰撞力的PairMOdel有没有什么参考资料？？谢谢

在使用DPMFoam计算管道颗粒流时，单核计算比较慢，可以使用把管道分块的并行计算方法吗？这种方法提高计算速度可行吗？谢谢各位老师

@yuanlee2011 您好，请问您在DPMFoam中添加粒子标量温度T实现了吗，最近在做这方面内容，希望得到一些指点

@上级
您好，请问您上面的公式 F=Sp(U−Ud)+Su用来计算颗粒受力的，这里面的U和Ud分别表示的是什么东西？和KinematicParcel中的Uc_以及U_是对应关系吗

const parcelType& p = static_cast<const parcelType&>(*this);
 const forceSuSp Fcp = forces.calcCoupled(p, dt, mass, Re, mu);
    const forceSuSp Fncp = forces.calcNonCoupled(p, dt, mass, Re, mu);
    const forceSuSp Feff = Fcp + Fncp;
    const scalar massEff = forces.massEff(p, mass);
    // Update velocity - treat as 3-D
    const vector abp = (Feff.Sp()*Uc_ + (Feff.Su() + Su))/massEff;
    const scalar bp = Feff.Sp()/massEff;

    Spu = dt*Feff.Sp();

    IntegrationScheme<vector>::integrationResult Ures =
        td.cloud().UIntegrator().integrate(U_, dt, abp, bp);

    vector Unew = Ures.value();

    // note: Feff.Sp() and Fc.Sp() must be the same
    dUTrans += dt*(Feff.Sp()*(Ures.average() - Uc_) - Fcp.Su());

    // Apply correction to velocity and dUTrans for reduced-D cases
    const polyMesh& mesh = td.cloud().pMesh();
    meshTools::constrainDirection(mesh, mesh.solutionD(), Unew);
    meshTools::constrainDirection(mesh, mesh.solutionD(), dUTrans);

    return Unew;
}

如果我想要计算颗粒所受的ErgunWenYuDrag这个力的时候是不是就等于

(mass/p.rho() *(150.0*(1.0 - alphac)/alphac +1.75*Re)*muc/(alphac*sqr(p.d()))*（U-Ud）//这样就OK了吗？

在DPM中粒子所受力例如重力、虚拟质量力等都是通过forceSuSp中的vector Su以及scalar Sp来决定的，最后通过KinematicParcel.C来计算Fcp以及Fncp(这个地方是每个力单独耦合求解速度还是把所有的力都考虑进来一起耦合),不知道理解的对不对，现在我想要把每个粒子的单独受力比如重力、虚拟质量力等都单独写出来，不知道各位老师有没有什么办法去实现？Ps:我看另一个帖子中写到F=Su*(U-Ud)+Sp.这个是在哪里出现的我也没找到，Ud是啥，可以直接用吗？

这是单独一个力的以重力为例：

template<class CloudType>
Foam::forceSuSp Foam::GravityForce<CloudType>::calcNonCoupled
(
    const typename CloudType::parcelType& p,
    const scalar dt,
    const scalar mass,
    const scalar Re,
    const scalar muc
) const
{
    forceSuSp value(Zero, 0.0);

    value.Su() = mass*g_*(1.0 - p.rhoc()/p.rho());//这里求出来Su

    return value;
}

这是KinematicParcel.C出现通过力算速度的：

const forceType& forces = td.cloud().forces();
//这是最后粒子受力耦合求解速度
    // Momentum source due to particle forces
    const parcelType& p = static_cast<const parcelType&>(*this);
//这个地方是每个力单独耦合求解速度还是把所有的力都考虑进来一起耦合？？
    const forceSuSp Fcp = forces.calcCoupled(p, dt, mass, Re, mu);
    const forceSuSp Fncp = forces.calcNonCoupled(p, dt, mass, Re, mu);
    const forceSuSp Feff = Fcp + Fncp;
    const scalar massEff = forces.massEff(p, mass);
    // New particle velocity
    //~~~~~~~~~~~~~~~~~~~~~~
    // Update velocity - treat as 3-D
    const vector abp = (Feff.Sp()*Uc_ + (Feff.Su() + Su))/massEff;
    const scalar bp = Feff.Sp()/massEff;

    Spu = dt*Feff.Sp();

    IntegrationScheme<vector>::integrationResult Ures =
        td.cloud().UIntegrator().integrate(U_, dt, abp, bp);

    vector Unew = Ures.value();