ガス拡散の解析のためのpimpleFoamコード改良

[Yota]

皆様、

こんにちは、お世話になっております。

現在コードを改良中でして、皆様のお知恵をお借りしたいと思い投稿いたしました。

LES解析に汚染ガス（パッシブスカラ）の拡散を盛り込もうと思い、

pimpleFoam.CにCEqn.Hを#includeさせるという形で輸送方程式を盛り込みました。

CEqn.Hは以下の通りです。

----------------------------------------------------------------------------------------------------------------

{

    volScalarField DCEff("DCEff", turbulence->nuEff()/Sct);

    fvScalarMatrix CEqn

    (

        fvm::ddt(C)

      + fvm::div(phi, C)

      - fvm::laplacian(DCEff, C)

    );

    CEqn.relax();

    CEqn.solve();

}

-----------------------------------------------------------------------------------------------------------------

質問したいところは1行目の「turbulence->nuEff()」についてです。

私は以前まで、「turbulence->nu()/Sct + turbulence->nuSgs()/Sct」

としていたのですが、コンパイル時にエラーで「nuSgs()」は引っ張って

これないと言われました。

そこでコードをたどっていたら、nuEffを見つけ、説明にも「nuEff = nu + nuSgs」と

書かれていたので、最終的に上記のようなコードになったわけですが、

これで大丈夫でしょうか？

もしご意見があればいただけると大変うれしいです。

どうぞよろしくお願いいたします。

[ohbuchi]

こんにちは。

この表記で問題ありません。

pimpleFoamは、LES、RASどちらも扱える様に、turbulenceModelというアブストラクト・クラスを介して

有効粘性係数にアクセスしています。

実行時にRASモデルを選べば渦粘性＋動粘性で、LESを選べばSGS粘性＋動粘性で評価されます。

2014年4月24日木曜日 8時25分54秒 UTC+2 Yota:

[Yota]

Ohbuchi様

教えていただきどうもありがとうございます！

C++には不慣れなもので少し心配していました。

RANSでもLESでも使えるようにするための表記だったのですね。

これで心置きなく解析が回せます。

どうもありがとうございました。

2014年4月30日水曜日 17時38分02秒 UTC+9 ohbuchi:

[ONO Hiroki]

こんにちは。

>>私は以前まで、「turbulence->nu()/Sct + turbulence->nuSgs()/Sct」

とありますが，シュミット数と乱流シュミット数は同値なのでしょうか？それともtypo？同値なら上記で問題ありませんが，違う場合は問題になります。

確かv2.0ぐらいからだったと思いますが，LESの場合でもturbulence->nut()でnuSgsにアクセスできるようになっています。

2014年5月7日水曜日 20時59分14秒 UTC+9 Yota:

[Yota]

ONO様、

返信遅くなり申し訳ございません。

ご意見ありがとうございます！

シュミット数と乱流シュミット数ですが、現在は同値で行っているので、そのまま解析を進めているのですが、

今後、それぞれ変更する可能性もあるので、できればばらしておきたいところです。

>> 確かv2.0ぐらいからだったと思いますが，LESの場合でもturbulence->nut()でnuSgsにアクセスできるようになっています。

このプログラムはv2.1で改造しまして、最近v2.3にバージョンアップしたのでもう一回確かめてみようと思います！

貴重なご意見ありがとうございました。

Yota

2014年6月3日火曜日 17時28分33秒 UTC+9 ONO Hiroki: