particleforce类内ErgunWenYuDrag的表达式感觉有错误

[upc_ngh]

@上级
您好，请问您上面的公式 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了吗？

[上级]

@upc_ngh
1.是的，$\mathbf{U}$代表的是KinematicParcel中的$\mathbf{U_c}$，也就是连续相的在颗粒所处位置的速度插值，具体的插值方式OF里面有cell，cellpoint那些，$\mathbf{U}_\rm{d}$表示的就是颗粒的速度。
2.是的，系数乘以速度差。