多相流求解器中默认的曳力模型与文献不一致？

[四季之中的隐星]

如题，检查曳力模型代码
位于dragModel.C

Foam::tmp<Foam::volScalarField> Foam::dragModel::Ki() const
{
    return
        0.75
       *CdRe()
       *swarmCorrection_->Cs()
       *pair_.continuous().rho()
       *pair_.continuous().thermo().nu()
       /sqr(pair_.dispersed().d());
}

两个问题，一是为什么有粘度，文献中基本是不考虑粘度的把，第二个问题是为什么液滴直径要开根号？
这里以李东岳老师的教程为例，文献形式一般为

检查Cd的计算，与文献中一致。检查到方程的部分，也是直接乘以的速度。
前后都没有问题，所以这里为什么会有粘度和被开了根号的直径呢？

[李东岳]

OpenFOAM里面的植入，植入的是CdRe，所以后面还有Re的部分，这样就一直了

[四季之中的隐星]

@李东岳 谢谢李老师，但是还是有一点问题。
Re的计算有：phasePair.C

Foam::tmp<Foam::volScalarField> Foam::phasePair::Re() const
{
    return magUr()*dispersed().d()/continuous().thermo().nu();
}

CdRe的计算有
SchillerNaumann.C

Foam::tmp<Foam::volScalarField> Foam::dragModels::SchillerNaumann::CdRe() const
{
    volScalarField Re(pair_.Re());

    return
        neg(Re - 1000)*24*(1.0 + 0.15*pow(Re, 0.687))
      + pos0(Re - 1000)*0.44*max(Re, residualRe_);
}

最后，K的计算有
dragModel.C

Foam::tmp<Foam::volScalarField> Foam::dragModel::Ki() const
{
    return
        0.75
       *CdRe()
       *swarmCorrection_->Cs()
       *pair_.continuous().rho()
       *pair_.continuous().thermo().nu()
       /sqr(pair_.dispersed().d());
}


Foam::tmp<Foam::volScalarField> Foam::dragModel::K() const
{
    return max(pair_.dispersed(), pair_.dispersed().residualAlpha())*Ki();
}

综上，Kf的表达式为：

比较您的文章中的

可以发现液滴直径的参与方式上不同，这个是什么情况，您了解吗

[四季之中的隐星]

@李东岳 谢谢李老师，我把sqr函数看成sqrt了，现在我明白了，是能对上的