Tomiyama升力模型bug

[李东岳]

Tomiyama升力模型可针对离散相不同的粒径预测不同的曳力系数，在直径分布区间较大的情况下，Tomiyama升力模型预测的升力系数甚至是反向的。

OpenFOAM中的Tomiyama升力模型有一个小bug，此bug导致对于大直径离散相，其升力系数依旧是正值。结果是相分数分布反向。做如下修正：

Foam::tmp<Foam::volScalarField> Foam::liftModels::TomiyamaLift::Cl() const
{
    // volScalarField EoH(pair_.EoH2());
    volScalarField EoH(pair_.EoH1());

    volScalarField f
    (
        0.0010422*pow3(EoH) - 0.0159*sqr(EoH) - 0.0204*EoH + 0.474
    );

    return
        neg(EoH - scalar(4))*min(0.288*tanh(0.121*pair_.Re()), f)
      + pos0(EoH - scalar(4))*neg(EoH - scalar(10.7))*f
      + pos0(EoH - scalar(10.7))*(-0.288);
}

未更正结果（相分数靠近壁面较大，不合理）

更正后结果（相分数靠近管道中心，合理）

黑线为CFD值，点为实验值

[李东岳]

这个bug在OpenFOAM-6里面还存在

[金石为开]

李老师 这是什么bug呢 仅仅改变了EoH1() 我以前计算用的原来的升力模型 难道以前计算的错了？？

[李东岳]

5mm一下直径没关系，5mm以上直径的颗粒有bug

[金石为开]

@东岳 看了看只有EoH1和EoH2的区别 里面cbrt是什么函数 没看懂

[李东岳]

1/3次幂咯