雷诺（粘性）应力-fvc::div((nuEff*dev2(T(fvc::grad(U)))))项编译报错

[SuperBoy]

各位老师好， 我最近将动量方程中turbulence->divDevReff(U)这一项还原成-fvc::div((nuEff*dev2(T(fvc::grad(U)))))-fvm::laplacian(nuEff, U)，在编译时出现了报错的情况，如下图所示：

我不是很清楚报错的原因，请问有老师知道该怎么解决吗？

[SuperBoy]

@SuperBoy >divDevReff(U)的代码参考自linearViscousStress.C：

template<class BasicTurbulenceModel>
Foam::tmp<Foam::fvVectorMatrix>
Foam::linearViscousStress<BasicTurbulenceModel>::divDevRhoReff
(
    volVectorField& U
) const
{
    return
    (
      - fvc::div((this->alpha_*this->rho_*this->nuEff())*dev2(T(fvc::grad(U))))
      - fvm::laplacian(this->alpha_*this->rho_*this->nuEff(), U)
    );
}

[李东岳]

你的报错应该不是最终的，有一个最基本的，你的nuEff就没有定义。

[SuperBoy]

@李东岳 李老师好，我的nuEff定义在了构建动量方程之前：

	volScalarField nu = turbulence->nu();
	volScalarField nut = turbulence->nut();
	volScalarField nuEff = nu + nut;

而且，报错信息中没有指出nuEff未被声明，所以我感觉不是这个的原因

[SuperBoy]

@李东岳 李老师，我发现应该是dev2(T(……)）这一小段出的问题，删掉这个转置操作就能编译成功。因为我是在传热求解器的基础上进行修改的，T在createFields.H里声明了是温度场，所以是不是程序将T识别成了温度而进行不了转置操作了呢？