周期境界の境界条件の設定

[Y.Mura]

初めてこちらに書き込ませていただきましたY.Mと申します。
大学で流体の数値解析をするために数ヶ月前よりopenfoamを使わせてもらっています。

現在行っている計算が大変時間がかかるため出口の条件を入り口側に持っていく周期境界を設定しようと考え、その前段階として液単相の円管で周期境界を試
してみようとしているところで問題が起こったのでこちらに答えていただきたく投稿しました。

計算はpisoFoamで回しており、乱流モデルには標準k-epsilonモデルをもちいています。
blockMeshDictのpatchesの部分は以下のようになっています。

patches
(
wall
wall
(
(@ @ @ @)
(@ @ @ @)
)
cyclic
inout
(
(@ @ @ @)
(@ @ @ @)
)
);

また、0ディレクトリの各パラメータのinoutの境界条件のtypeはすべてcyclicになっています。
速度に関しては全体に入り口で流入する予定の流速を設定し、
圧力に関しては全体が0になっています。

この状態で計算を回すとだんだんと流速が遅くなってしまいました。
いろいろ調べたり聞いたりした結果圧力の境界条件の問題ではないかと考えているのですが、
どのようにすれば解決できるかがわかりませんでした。
こちらの解決方法がわかる方がいらっしゃいましたらご教授お願いいたします。

[Youhei Takagi]

周期境界条件で円管内乱流やチャネル乱流で計算するときは、

流量固定か圧力勾配一定の流入・流出条件になるかと思います。

自分の要望にあった境界条件を選択し、境界条件の設定をお願い

します。

2010年12月20日20:12 Y.Mura <dfhw...@kcc.zaq.ne.jp>:

--
このメールは Google グループのグループ「OpenFOAM」の登録者に送られています。
このグループに投稿するには、open...@googlegroups.com にメールを送信してください。
このグループから退会するには、openfoam+u...@googlegroups.com にメールを送信してください。
詳細については、http://groups.google.com/group/openfoam?hl=ja からこのグループにアクセスしてください。

[ohbuchi]

jumpCyclicという境界条件を使うと圧力降下を指定できそうです。

また、imcompressible/pisoFoam/les/pitzDailyDirectMappedチュートリアル
の様に速度をdirectMappdedで指定したオフセット下流の値をコピーすれば
圧力は通常通り出口で規定し、入口はZeroGradientでよいのではないかと
思います。

ご参考まで。

[Y.Mura]

Youhei Takagi さん、ohbuchi さんご回答ありがとうございました。

patch を cyclic　にすると　patchField　に圧力勾配一定などが設定できないためうまくいっていませんでしたが、
ohbuchi さんが書かれていた
imcompressible/pisoFoam/les/pitzDailyDirectMappedチュートリアル
を参考に設定したところ流速が下がることなく計算を回すことができました。
すばやいご返答ありがとうございました。