LES smagorinskyモデルの設定について

[kure]

お世話になっております、
僕は大学院で流体数値解析の研究をしている学生でkureと申します。


すみません、

OpneFOAMのLESモデルでsmagorinskyモデルの設定について
お聞きしたいのですが、

いろいろウェブなどで検索しながらsmagorinskyモデルを設定しようとしているのですが、
（理由は比較する解析でこのモデルを使用していたからです。）
いろいろ試行錯誤した結果、turbulencePrpopertiesで下のような設定をしました。

これってsmagorinskyモデルなのでしょうか？
OpenFOAMは4.0バージョンを使用しています。
結構前のバージョンでは明確にLESModelにsmagorinskyと入れていたように見えますが
バージョンが変わってこのようになったのか。。
それとも自分が間違っているのか教えて頂けないでしょうか？

また下の乱流設定では0.origでk（乱流エネルギー）ファイルを設定しなくては
いけないようですが、、smagorinskyモデルでkを使用するのはあっているのでしょうか？
乱流エネルギーや乱流喪失率は使わずメッシュサイズとCsで渦粘性係数を決めると
知っていたのでさらに確信が持てない状態です。
（Smagorinsky.Cを開いて見るとkが計算されているように見えましたが、epsilonも見えたので
さらに混乱しています。。）

すみません、どなたかご存知の方教えて頂ければ幸いです。

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  4.0                                   |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    location    "constant";
    object      turbulenceProperties;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

simulationType  LES;

LES
{
    LESModel        kEqn;

    turbulence      on;

    printCoeffs     on;

    delta           cubeRootVol;

    cubeRootVolCoeffs
    {
        deltaCoeff      1;
    }

    PrandtlCoeffs
    {
        delta           cubeRootVol;
        cubeRootVolCoeffs
        {
            deltaCoeff      1;
        }

        smoothCoeffs
        {
            delta           cubeRootVol;
            cubeRootVolCoeffs
            {
                deltaCoeff      1;
            }

            maxDeltaRatio   1.1;
        }

        Cdelta          0.158;
    }

    vanDriestCoeffs
    {
        delta           cubeRootVol;
        cubeRootVolCoeffs
        {
            deltaCoeff      1;
        }

        smoothCoeffs
        {
            delta           cubeRootVol;
            cubeRootVolCoeffs
            {
                deltaCoeff      1;
            }

            maxDeltaRatio   1.1;
        }

        Aplus           26;
        Cdelta          0.158;
    }

    smoothCoeffs
    {
        delta           cubeRootVol;
        cubeRootVolCoeffs
        {
            deltaCoeff      1;
        }

        maxDeltaRatio   1.1;
    }
}

// ************************************************************************* //

[すもも]

マニュアルにあるとおり、kEqn は Smagorinsky とは別のモデルのようですが。

https://cfd.direct/openfoam/user-guide/turbulence/

2017年6月17日土曜日 15時12分17秒 UTC+9 kure: