Smagorinsky模型不输出k_sgs场

[李东岳]

可以做如下改动进行输出，至少适用于OpenFOAM-8

1. 打开Smagorinsky.H，添加k_:

    // Protected data

        dimensionedScalar Ck_;
        
        volScalarField k_;//这是添加的一行


    // Protected Member Functions

2. 打开Smagorinsky.C，添加下面代码:

:
    LESeddyViscosity<BasicMomentumTransportModel>
    (
        type,
        alpha,
        rho,
        U,
        alphaRhoPhi,
        phi,
        transport
    ),

    Ck_
    (
        dimensioned<scalar>::lookupOrAddToDict
        (
            "Ck",
            this->coeffDict_,
            0.094
        )
    ),//这里加个逗号
    
    //下面是添加的
    k_
    (
        IOobject
        (
            IOobject::groupName("k", this->alphaRhoPhi_.group()),
            this->runTime_.timeName(),
            this->mesh_,
            IOobject::MUST_READ,
            IOobject::AUTO_WRITE
        ),
        this->mesh_
    )
{
    if (type == typeName)
    {
        this->printCoeffs(type);
    }
}

correctnut()函数通过下面代码替换：

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给我趴下]

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

[李东岳]

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你就得按照主贴的方式试一下了