LESモデルでのB,nuTildaの意味

1,532 views
Skip to first unread message

Sakuma

unread,
May 6, 2011, 8:09:48 PM5/6/11
to OpenFOAM
 OpenFoam初心者のsakumaです。

 LESの勉強としてチュートリアルにあるpisoFoam/les/pitzDailyを見ていますが,
良く理解できない設定があります。

 ①0ホルダにあるBの設定ファイルは何を意味しているのでしょうか。p,kと同じ次元を持つようですが,RASのチュートリア ルである
pisoFoam/ras/cavityを見ますとRがあり,これはレイノズル応力と理解していますが,これに対応するものでしょ うか。

 Rは対称テンソルになっているため,6個の独立した値を入れる事はわかりますが,Bは9個の独立した値を入れるようにな っています。Bはどのような
式で表されているのでしょうか。

 ②nuTildaはSpalart-Allmarasモデル用で,他のモデルを使う時は必要ないと考えてよいでしょうか。

よろしくお願いいたします。

Masashi Imano

unread,
May 7, 2011, 2:16:50 AM5/7/11
to open...@googlegroups.com
今野です。

2011/5/7 Sakuma <sakura....@gmail.com>:


>  OpenFoam初心者のsakumaです。
>
>  LESの勉強としてチュートリアルにあるpisoFoam/les/pitzDailyを見ていますが,
> 良く理解できない設定があります。
>
>  ①0ホルダにあるBの設定ファイルは何を意味しているのでしょうか。p,kと同じ次元を持つようですが,RASのチュートリア ルである
> pisoFoam/ras/cavityを見ますとRがあり,これはレイノズル応力と理解していますが,これに対応するものでしょ うか。
>
>  Rは対称テンソルになっているため,6個の独立した値を入れる事はわかりますが,Bは9個の独立した値を入れるようにな っています。Bはどのような
> 式で表されているのでしょうか。

これについては、第10回のOpenFOAM(R)勉強会 for beginnerでの小野さんの発表で質問が出ていました。

http://www.ustream.tv/recorded/14241568

の1:43:48辺りからご覧ください。

なお、${FOAM_SRC}/turbulenceModels/incompressible/LES/laminar/laminar.H

//- Return the sub-grid stress tensor B.
virtual tmp<volSymmTensorField> B() const;

とある通り、Bはサブグリッド応力テンソルで、対称テンソルです。

なぜかtutorialsのpisoFoam/les/pitzDailyの0/BはvolTensorFieldになってますが、volSymmTensorFieldの間違いです。
実際、Bを読み込むDeardorffDiffStressなどのLESモデルにするとエラーになります。
OF-1.3ではvolSymmTensoFieldrがなくて、volTesorFieldとして定義していたので、その時のBのファイルがそのまま
残っているのバグではないかと思います。

>  ②nuTildaはSpalart-Allmarasモデル用で,他のモデルを使う時は必要ないと考えてよいでしょうか。

grep -r nuTilda ${FOAM_SRC}/turbulenceModels/incompressible/LES/

等としてnuTildaが出てくるソースを調べると、SpalartAllmarasしかないことがわかると思います。

Sakuma

unread,
May 8, 2011, 4:05:38 AM5/8/11
to OpenFOAM
今野様 OpenFoam初心者のsakumaです。

①Bについては http://www.ustream.tv/recorded/14241568 を拝聴させて頂き,概要を理解しました。
②nuTildaは お教え頂きましたソースを見ました。

良く理解出来ていない部分が大分わかりました。ありがとうごさいました。



Reply all
Reply to author
Forward
0 new messages