Smagorinsky模型不输出k_sgs场

[ZZ给我趴下]

@李东岳 李老师，请问为什么做了这些改变之后会编译错误呢？

[李东岳]

156-167行删掉，用下面的替换：

void Smagorinsky<BasicMomentumTransportModel>::correctNut()
{
    k_ = (this->k(fvc::grad(this->U_)));

    this->nut_ = Ck_*this->delta()*sqrt(k_);
    this->nut_.correctBoundaryConditions();
    fv::options::New(this->mesh_).correct(this->nut_);
}

[ZZ给我趴下]

@李东岳这里提示重复定义，在59行有定义，而且和这个函数一样。

[李东岳]

更简单的方法就是不用动求解器了，在算例文件下面的controlDict下面把这个放进去：

cacheTemporaryObjects
(
    k
);

functions
{
    #includeFunc writeObjects("k")
}

然后直接跑就行

[疏影横斜水清浅]

@李东岳epsilon()函数要不要改

template<class BasicTurbulenceModel>
tmp<volScalarField> tshbSmagorinsky<BasicTurbulenceModel>::epsilon() const
{
   k_ = (this->k(fvc::grad(this->U_)));

   return volScalarField::New
   (
       IOobject::groupName("epsilon", this->alphaRhoPhi_.group()),
       this->Ce_*k_*sqrt(k_)/this->delta()
   );
}

[疏影横斜水清浅]

@李东岳 我在controlDict中写入

cacheTemporaryObjects
(
    k
);

functions
{
    #includeFunc writeObjects("k")
}

结果出现错误，请问李老师，能知道是什么原因吗

[李东岳]

你这个是openfoam7么，openfoam7可能不能这么用，你有没有openfoam10，你试一下

[疏影横斜水清浅]

@李东岳 是的，我使用的OpenFOAM7，我安装个10试试。

[李东岳]

openfoam7你就得按照主贴的方式试一下了

[78star]

@ZZ给我趴下 在 Smagorinsky模型不输出k_sgs场 中说：

@李东岳这里提示重复定义，在59行有定义，而且和这个函数一样。

老帖子了，为了方便后面的人看，这里指出下问题。原始的代码159行这里是调用correctNut，而它的实现是在上面的59行。李老师的意思是换掉59行位置correctNut的函数实现代码里面的k_，不是把此处调用correctNut换成另一个新的实现，因此会报重复定义错误