CdRe意思？

[youhaoyu]

最近在看颗粒受力模型，对于spheredragforce

最底下的曳力公式似乎是
这样的
但是这个量纲好像不是力的量纲，看了看其他力的量纲都是牛顿。
对比了一下其他论文里面的这个公式，后面的都有（Ul-Up），这样量纲就对了
是不是这个cdRe不是cd*Re（之前也看讨论过），应该是cdRe还有个绝对速度？？

template<class CloudType>
Foam::scalar Foam::SphereDragForce<CloudType>::CdRe(const scalar Re) const
{
    // (AOB:Eq. 35)
    if (Re > 1000.0)
    {
        return 0.424*Re;
    }
 
    return 24.0*(1.0 + (1.0/6.0)*pow(Re, 2.0/3.0));
}

上面这段代码也没看明白，始终没搞懂of是怎么编写这个曳力方程的

[李东岳]

http://dyfluid.com/reactingTwoPhaseEulerFoam.html

公式5 6 7 就是Cd乘以Re

[youhaoyu]

@李东岳
公式中不是有这项么，为什么我在代码里面没看出现呢，可能是我代码功力太弱了，我想的是不应该还
乘以"(U-p.U())"

[youhaoyu]

@youhaoyu
这行代码我没对应上

[lyc]

@youhaoyu 就是SchillerNaumann曳力系数模型Cd·Re

[youhaoyu]

@lyc 您这个是曳力系数表达式，我在想曳力公式那一块

return forceSuSp(Zero, mass*0.75*muc*CdRe(Re)/(p.rho()*sqr(p.d())));

主要是这个代码没看懂，曳力方程里面不是还有个速度之差么，这个代码没体现

[李东岳]

MomentumParcelI.H

template<class ParcelType>
inline Foam::scalar Foam::MomentumParcel<ParcelType>::Re
(
    const scalar rhoc,
    const vector& U,
    const vector& Uc,
    const scalar d,
    const scalar muc
)
{
    return rhoc*mag(U - Uc)*d/max(muc, rootVSmall);
}