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换成另一个新的实现，因此会报重复定义错误

CFD中文网

Smagorinsky模型不输出k_sgs场