bunrs湍流耗散力

anubis

com版本和org版本的burns湍流耗散力好像不太一样

OpenFOAM8里的

Foam::turbulentDispersionModels::Burns::D() const
{
    const fvMesh& mesh(pair_.phase1().mesh());
    const dragModel& drag =
        mesh.lookupObject<dragModel>
        (
            IOobject::groupName(dragModel::typeName, pair_.name())
        );

    return
        drag.Ki()
       *continuousTurbulence().nut()
       /sigma_
       *pair_.dispersed()
       *sqr(pair_.dispersed() + pair_.continuous())
       /(
            max(pair_.dispersed(), pair_.dispersed().residualAlpha())
           *max(pair_.continuous(), pair_.continuous().residualAlpha())
        );
}

OpenFOAM2012里的

Foam::turbulentDispersionModels::Burns::D() const
{
    const fvMesh& mesh(pair_.phase1().mesh());
    const dragModel&
        drag
        (
            mesh.lookupObject<dragModel>
            (
                IOobject::groupName(dragModel::typeName, pair_.name())
            )
        );

    return
        0.75
       *drag.CdRe()
       *pair_.continuous().nu()
       *pair_.continuous().turbulence().nut()
       /(
            sigma_
           *sqr(pair_.dispersed().d())
        )
       *pair_.continuous().rho()
       *pair_.dispersed()
       *(
           1.0/max(pair_.dispersed(), residualAlpha_)
         + 1.0/max(pair_.continuous(), residualAlpha_)
        );
}

公式应该是

按照这个公式，感觉2012版本里的比较接近，但是还不太一样，相对速度Ur，在2012的代码里好像变成了连续相的粘度除以离散相直径，虽然量纲和速度是一样的，但是不太明白为什么会这样？

anubis

又仔细看了一下，两个版本应该是差不多的。
相对速度是因为CdRe不是曳力系数而是曳力系数×雷诺数所以可以用粘度、直径还有密度换算成相对速度。