reactingTwoPhaseEulerFoam添加空化模型

[冠竹]

各位大佬，我尝试在reactingTwoPhaseEulerFoam中添加空化模型，参考的是
www.tfd.chalmers.se/~hani/kurser/OS_CFD_2017/SuryaKaundinyaOruganti/Final_presentation_Surya_Oruganti.pdf
中的方式，但是在修改phaseChangeModel.C文件时产生了报错。
报错如下

这部分我的代码与参考资料中给的一致：

namespace Foam
{
    defineTypeNameAndDebug(phaseChangeModel, 0);
    defineRunTimeSelectionTable(phaseChangeModel, components);
}

// * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //

Foam::phaseChangeModel::phaseChangeModel
(
    const word& type,
    const dictionary& interfaceDict,
    const phaseModel& phase1,
    const phaseModel& phase2

)
:
    interfaceDict_(interfaceDict),
    phase1_(phase1),
    phase2_(phase2),
    phaseChangeModelCoeffs_(interfaceDict.optionalSubDict(type + "Coeffs")),
    pSat_("pSat", dimPressure, interfaceDict.lookup("pSat"))

{}


// * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * * //

Foam::Pair<Foam::tmp<Foam::volScalarField>>
Foam::phaseChangeModel::vDotAlphal() const
{

    volScalarField alphalCoeff(1.0/phase1_.rho() - phase1_*(1.0/phase1_.rho() - 1.0/phase2_.rho()));
    Pair<tmp<volScalarField>> mDotAlphal = this->mDotAlphal();

    return Pair<tmp<volScalarField>>
    (
        alphalCoeff*mDotAlphal[0],
        alphalCoeff*mDotAlphal[1]
    );
}

Foam::Pair<Foam::tmp<Foam::volScalarField>>
Foam::phaseChangeModel::vDotP() const
{
    dimensionedScalar pCoeff(1.0/phase1_.rho() - 1.0/phase2_.rho());
    Pair<tmp<volScalarField>> mDotP = this->mDotP();

    return Pair<tmp<volScalarField>>(pCoeff*mDotP[0], pCoeff*mDotP[1]);
}

这种问题是什么原因导致的呢？我刚开始接触这种模型编译，不太懂

[李东岳]

openfoam这面要注意版本问题。你这个链接里面我看起码是2017年之前的文章。然后你用的openfoam2206。这里面版本不对。按照之前的弄，100%会有问题。