求助大佬，为了在compressibleInterFoam求解器中添加有效粘度muEff

surfaceScalarField muEff
    (
        "muEff",
         mixture.muf()      
      + fvc::interpolate(rho*turbulence->nut())
    );

需要定义muf()函数，mixture是twoPhaseMixtureThermo类型，在对应的H文件和C文件进行更改

virtual tmp<surfaceScalarField> muf() const;     //muf()声明
Foam::tmp<Foam::surfaceScalarField> Foam::twoPhaseMixtureThermo::muf() const   //muf()定义
{
    const surfaceScalarField alpha1f     
    (
        min(max(fvc::interpolate(alpha1_), scalar(0)), scalar(1))
        
    );

    return tmp<surfaceScalarField>    
    (
        new surfaceScalarField   
        (
            "muf",
            alpha1f*(thermo1_->rho())*fvc::interpolate(nuModel1_->nu())          
          + (scalar(1) - alpha1f)*(thermo2_->rho())*fvc::interpolate(nuModel2_->nu())
        )
    );
}

编译时主要报这些错误该如何解决？为什么类似的代码在incompressibleTwoPhaseMixture中可以编译？

/home/dell/OpenFOAM/OpenFOAM-4.x/src/OpenFOAM/lnInclude/GeometricField.H:99:40: error: invalid use of incomplete type ‘class Foam::surfaceMesh’
         typedef typename GeoMesh::Mesh Mesh;
/home/dell/OpenFOAM/OpenFOAM-4.x/src/OpenFOAM/lnInclude/GeometricField.H:103:48: error: invalid use of incomplete type ‘class Foam::surfaceMesh’
         typedef typename GeoMesh::BoundaryMesh BoundaryMesh;
twoPhaseMixtureThermo.C:83:17: error: ‘fvc’ has not been declared
         min(max(fvc::interpolate(alpha1_), scalar(0)), scalar(1))

CFD中文网

Miss Zhou_DUT