devReff()通过devRhoReff()实现其功能（rho =1）？（OF4.x）

Yongbo

IncompressibleTurbulenceModel.C 中的代码：
首先给出devReff() 和devRhoReff()的代码，如下：

template<class TransportModel>
Foam::tmp<Foam::volSymmTensorField>
Foam::IncompressibleTurbulenceModel<TransportModel>::devReff() const
{
    return devRhoReff();
}

...

template<class TransportModel>
Foam::tmp<Foam::volSymmTensorField>
Foam::IncompressibleTurbulenceModel<TransportModel>::
devRhoReff() const
{
    NotImplemented;

    return devReff();
}

可以看到，devReff()调用了 devRhoReff()。然而在IncompressibleTurbulenceModel.C中，并没有找到这两个函数更进一步的细节

接着，在linearViscousStress.C 中的的代码中，
发现devRhoReff()的定义

template<class BasicTurbulenceModel>
Foam::tmp<Foam::volSymmTensorField>
Foam::linearViscousStress<BasicTurbulenceModel>::devRhoReff() const
{
    return tmp<volSymmTensorField>
    (
        new volSymmTensorField
        (
            IOobject
            (
                IOobject::groupName("devRhoReff", this->U_.group()),
                this->runTime_.timeName(),
                this->mesh_,
                IOobject::NO_READ,
                IOobject::NO_WRITE
            ),
            (-(this->alpha_*this->rho_*this->nuEff()))
           *dev(twoSymm(fvc::grad(this->U_)))
        )
    );
}

先将问题总结如下：
1）NotImplemented是什么意思？为什么不可压中devRhoReff()返回的是devReff()？
是否可以理解为，不可压缩流里面的devReff()执行的是其它代码（linearViscousStress.C）中的devRhoReff()？

template<class TransportModel>
Foam::tmp<Foam::volSymmTensorField>
Foam::IncompressibleTurbulenceModel<TransportModel>::
devRhoReff() const
{
    NotImplemented;

    return devReff();
}

2）若devReff() 实际上执行的是devRhoReff()（linearViscousStress.C）

并且在CFDonline一个帖子里，链接文本，有人说devReff()实际上是devRhoReff()中rho为1时的情况
''it seems like it's best to simply search for "devRhoReff()", because "devReff()" is just when "rho" is simply 1.0.''。
若是如此，这里面rho_=1是在哪里定义的，我一直没找到相应代码

望各位不吝赐教：）！

李东岳

手头没有OpenFOAM。简单来说，notimplementsd是一个宏，返回到一个error函数。ρ=1是因为定义了geometricOneField。你顺着这个思路看看。

Samuel-Tu

我记得构造哪个class的时候就是默认了1

devReff()通过devRhoReff()实现其功能 （rho =1）？（OF4.x）

devReff()通过devRhoReff()实现其功能（rho =1）？（OF4.x）