回転翼周辺の相対速度について

[小寺 正毅]

始めて相談申し上げます。
東京大学工学部産業機械工学科学部4年荒川研究室所属の小寺と申します。

現在OpenFOAMを用いてウェルズタービンの解析を行っております。
ソルバーにはpimpleDyMFoamを用い計算結果もきちんと出ているのですが、
翼周辺の流れ場を表示するとき絶対速度でなく回転翼との相対速度を表示させたいのですが、
どのユーティリティを使えばよいのかわからず困っております。

foamCalcから計算させようと試みたりもしましたが、差となる翼の相対速度場の設定の仕方がわからずお手上げです。
何か良い方法や別のユーティリティなどないでしょうか。

[ohbuchi]

こんにちは。CFD-Onlineの下記記事でMRFSimpleFoam用のUrelユーティリティが

公開されています。

http://www.cfd-online.com/Forums/openfoam-solving/71277-how-derive-relative-velocity-mrfsimplefoam.html#post350891

MRFZonesを参照していますが、dynamicMeshDictと整合性のあるMRFZonesを作れば

使えるのではないかと思います。

もしダメであれば、回転領域のセルゾーンに対してdynamicMeshDictから回転軸、回転中心、

回転角速度の情報を取得し、回転速度成分を求め、静止系速度から差し引くユーティリティ

プログラムを作れば良い筈です。

2012年11月10日土曜日 5時11分39秒 UTC+9 小寺　正毅:

[E.Mogura]

E.Mogura　です。

自分は、ParaViewのCalculatorを使って相対速度場を作り、

これを使ってベクトル表示したり流線処理などやってます。

Z軸周りで回転速度Omega　であれば、相対速度場は、

iHat*(U_X-coordsY*Omega)+jHat*(U_Y+coordsX*Omega)+kHat*U_Z

iHat,jHat,kHat は単位ベクトル

U_*　は速度成分

coords＊　座標値成分

で、いちいち手入力する必要はなく、

CalculatorパネルのScalarsボタンからプルダウンで選択入力できます。

軸が傾いていたりすると、ちょっと面倒ですけど・・・

2012年11月10日土曜日 5時11分39秒 UTC+9 小寺　正毅:

始めて相談申し上げます。

[小寺正毅]

ohbuchi様E.Mogura様

返答ありがとうございます。
どちらも大変参考になりました。
今回は表示させるだけでよいのでparaviewを使った方法を試してみたいと思います。
重ねてありがとうございました。

--
----
東京大学工学部産業機械工学科
学部4年
荒川研究室所属
小寺正毅