乱流二次元角柱の解析について

[tiyotsuru]

OpenFOAM初心者のtiyotsuruというものです。
大学の研究でOpenFOAMを使っているのですが、inletとoutletの乱流(標準のk-ε)の設定がわかりません。
ある程度は調べてみたのですが、

epsilon
R
K
nut(←この意味は？）
nuTilda

の値について

inlet
outlet
wall

の設定をどうすればいいのか情報をお願いします。

流速
up down
{ type fixedValue;
value uniform (10 0 0);
}
inlet
{
type fixedValue;
　 value uniform (10 0 0);
}
outlet
{
type zeroGradient;
}

kakuchuu
{
type fixedValue;
value uniform (0 0 0);
}

圧力
up down
{
type zeroGradient;
}

inlet
{type fixedValue ;
　　　 　value uniform 0 ;
}
outlet
{ type fixedValue ;
　 　　 value uniform 0.0;
}
kakuchuu
{
type zeroGradient;
}

[ohb...@amber.plala.or.jp]

こんにちは。
simpleFoamのチュートリアルと同じ様に設定すれば良いと思います。
ただ、R（レイノルズ応力）、nut,nuTilda(渦粘性係数）は0ディレクトリに含まれてますが
削除しても計算できる様に、OpenFOAMが勝手に作ってくれるものです。
kとepsilonだけを正しく設定すれば良い筈です。

入口で、値をfixedValueで定義し、出口でZeroGradientにするのが一番簡単な方法です。
turbulentIntensityKineticEnergyInletを使うと、kをIntensityで
turbulentMixingLengthDissipationRateInletを使うと、epsilonをmixingLengthで
与えることもできます。

以上、ご参考まで。

[Sakuma]

tiyotsuruさん，同じくOpenFOAM初心者のsakumaです。

k-εですが，2010年5月14日に開催されました第1回オープンCAE講習会資料
「DEXCS-OpenFOAMの基礎」のスライド61が参考になるかと思います。この
資料はオープンCAE学会のホームページからダウンロードできるようになって
います。

また，ohbuchiさんがご紹介されたturbulentIntensityKineticEnergyInletと
turbulentMixingLengthDissipationRateInlet
の書き方は，下記のチュートリアルに事例があるようです。
incompressible/pimpleFoam/t-junction/0のkとepsilon

以上，ご参考になればと思います。

[nakagawa]

すでにいくつかのアドバイスが投稿されていますが，追加で投稿します。

> nut(←この意味は？）
νt　渦粘性係数　です。

標準のk-eモデルでは，nuTildaは不要だと思います。

壁面では壁関数を使っているはずです。
下記マニュアルにも，関連する記述があります。
http://www.openfoam.com/docs/user/turbulence.php#x38-2050007.2
0ディレクトリにnutというファイルがあり，使用する壁関数が記載されていますよね？

inletでは，kとεを与える必要があります。
OpenFOAMのProgrammersGuideには，基本的な乱流の例題が説明されています。入り口での乱れ強さから，kはeを見積もっていま
す。
インストールしたOpenFOAMのドキュメントがあるディレクトリに，ProgrammersGuideが入っていると思います。
その中の「3.2　Steady turbulent flow over a backward-facing step」をごらんください。

outletでは，zerogradientでよいのでは？

[tiyotsuru]

お忙しい中、アドバイスありがとうございます。
ぜひとも参考にさせていただきます.