ＬＥＳモデルでのfoam waraningについて

[shin o-nishi]

OpenFOAM Google ユーザーの皆様へ

いつもお世話になっております。

Smagorinskyモデルでpisofoamを実行すると以下の文が端末に出ます。

FOAM Warning : 

    From function cubeRootVolDelta::calcDelta()

    in file cubeRootVolDelta/cubeRootVolDelta.C at line 52

    Case is 2D, LES is not strictly applicable

解析領域は奥行き方向に長さを与えているのに

なぜこの警告が出るのかが分かりません。

どうかご存知の方

よろしくお願いします。

[zheng zhang]

奥行き方向にメッシュも数個切ってあるでしょうか?もし切っていなかったら、単に2次元となりますが。

[shin o-nishi]

お早い御回答ありがとうございます。

解析領域は

blockMesh

vertices

(

    (-1.2 -2.2 -0.25)

    ( 5   -2.2 -0.25)

    ( 5    2.2 -0.25)

    (-1.2  2.2 -0.25)

    (-1.2 -2.2  0.25)

    ( 5   -2.2  0.25)

    ( 5    2.2  0.25)

    (-1.2  2.2  0.25)

);

blocks

(

    hex (0 1 2 3 4 5 6 7) (72 48 10) simpleGrading (1 1 1)

);

となっているので、奥行きには１０分割しています。

parafoamでも確認しました。

2013年9月4日水曜日 16時51分35秒 UTC+9 zheng zhang:

奥行き方向にメッシュも数個切ってあるでしょうか?もし切っていなかったら、単に2次元となりますが。

[ONO Hiroki]

向かい合う2組のパッチ(おそらく奥行き方向)がemptyになっていませんか？
emptyは2D計算用の境界条件です。

2013年9月4日水曜日 15時14分22秒 UTC+9 shin o-nishi:

[shin o-nishi]

御回答ありがとうございます。

u,p,k,nuSgsの奥行き面の境界条件がemptyになっていました。

しかし、emptyから適当に条件を変えたのですが

同じ警告文が出るのですが計算は回ります

ＬＥＳのチュートリアルのpitzDailyでは奥行き面をemptyにしていたので

あまり神経質にならないでいいのでしょうか

また、自分なりに原因を調べてみます

ありがとうございました

2013年9月4日水曜日 22時18分10秒 UTC+9 ONO Hiroki:

[ONO Hiroki]

constant/polyMesh/boundaryにもemptyになっているpatchがあると思いますので，
そちらも適切なものに変更してください。


完全に2次元の世界(そもそも奥行という概念そのものが存在しない空間)の中で，LESをやるのであればいいのですが，
3次元空間での乱流は，たとえ平均流れが2次元とみなせる場合でも3次元的に絶えず変動するのでLESをやる場合は
3次元で計算領域を定義しないといけません。

ところが，おっしゃる通りチュートリアルのpitzDailyは確かにemptyを用いた2次元計算となっています。
そういった意味ではpisoFoam/pitzDailyのチュートリアルは，LESのチュートリアルとしてし不適切なものになってしまっていますね。
おそらく，RANS用のものをコピー＆ペーストして最低限LESで計算実行できるところまで変更しただけでリリースされているのでしょう。

2013年9月4日水曜日 23時00分23秒 UTC+9 shin o-nishi:

[shin o-nishi]

ありがとうございます。

確かにboundaryがemptyになっていました。

計算も警告文も出ることなく流れました。

最後に質問なのですが

出来れば奥行き方向の影響を小さくして

２次元に近い計算をしたいと思っています。

そこで、emptyに相当するメッシュとかはあるのでしょうか？

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

2013年9月5日木曜日 16時11分01秒 UTC+9 ONO Hiroki:

[ONO Hiroki]

3次元計算でemptyにもっとも対応するものはslipですが，これでは摩擦0の壁が存在しているのと同じことになり，
パッチを横切る方向の流れ成分は殺されてしまいます。
LESではcyclic，すなわち周期境界条件を用いるのがよいと思います。
cyclicにすることで，パッチを横切る方向の乱れをある程度殺さずに計算ができます。
ただし，妥当な計算を行おうと思えば，奥行き方向にもある程度の計算領域と格子分割が必要です。


2013年9月5日木曜日 19時53分19秒 UTC+9 shin o-nishi:

[shin o-nishi]

ありがとうございます。

やはり奥行き方向にある程度の領域を確保しなければならないのですね。

質問に答えて頂いてありがとうございました。

2013年9月5日木曜日 22時52分46秒 UTC+9 ONO Hiroki: