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);
    }
    

    未更正结果(相分数靠近壁面较大,不合理)

    0_1542100859169_捕获.JPG

    更正后结果(相分数靠近管道中心,合理)

    0_1542100864999_捕获2.JPG

    黑线为CFD值,点为实验值



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



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



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



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



  • 1/3次幂咯


Log in to reply