wallLubrication模型中有bug

[李东岳]

https://github.com/OpenFOAM/OpenFOAM-5.x/blob/master/applications/solvers/multiphase/reactingEulerFoam/interfacialModels/wallLubricationModels/Antal/Antal.C
第73行，力的方向搞反了，应该改为

Foam::tmp<Foam::volVectorField> Foam::wallLubricationModels::Antal::Fi() const
{
    volVectorField Ur(pair_.Ur());

    const volVectorField& n(-nWall());//负号

    return zeroGradWalls
    (
        max
        (
            dimensionedScalar("zero", dimless/dimLength, 0),
            Cw1_/pair_.dispersed().d() + Cw2_/yWall()
        )
       *pair_.continuous().rho()
       *magSqr(Ur - (Ur & n)*n)
       *n
    );
}

[李东岳]

看错了看错了，符号在pEqn.H里面处理了