cartesian座標系で解析した結果を円筒座標系に変換する方法

S.kita

unread,

Nov 9, 2010, 2:46:22 AM11/9/10

to OpenFOAM

S.kitaといいます。どなたか、ご存知でしたら教えて下さい。

２D軸対称でモデリングした計算結果を円筒座標系に変換して後処理したいのですが、
OpenFOAM上のツールで解析結果を変換可能でしょうか？
また、paraviewに関してもよくわかっていないのですが、座標系を変換するやり方を
ご存知の方いらっしゃいましたら、教えて頂けると大変助かります。

宜しくお願いします。

ohb...@amber.plala.or.jp

unread,

Nov 9, 2010, 5:06:36 AM11/9/10

to OpenFOAM

どんな変換をお望みでしょうか？
例えばｘ軸を円筒中心軸にとり、y軸を半径方向にとる場合なら、
　Ux→Uz、　Uy→Ur
と読み替えるだけだと思います。

S.kita

unread,

Nov 9, 2010, 9:17:00 AM11/9/10

to OpenFOAM

ohbuchi様

ご回答ありがとうございます。
２D軸対称ならば、読み替えで十分だったことは恥ずかしながら
気づきませんでした。勉強になりました。

もしご存知でしたら、３次元の場合の円筒座標系変換の
方法に関してもご教授頂けれると幸いです。
（例えば、Xが円筒軸方向、Y・Zが軸直角面というような場合）

標準のUtilitityにはどうもなさそうで、paraviewのコマンドを
探してみようかと思いますが。

On 11月9日, 午後7:06, "ohbu...@amber.plala.or.jp"

ohb...@amber.plala.or.jp

unread,

Nov 9, 2010, 4:07:59 PM11/9/10

to OpenFOAM

おはようございます。
ｚ軸を円筒中心軸に合わせた場合、速度成分の変換は
　　Vr=Vx*cosθ＋Vy*sinθ
　　Vθ=Vy*cosθ－Vx*sinθ
　　Vz=Vz
とすれば良いと思います。

S.kita

unread,

Nov 9, 2010, 6:38:50 PM11/9/10

to OpenFOAM

ご回答ありがとうございました。

参考にして、変換処理をやってみようと思います。

On 11月10日, 午前6:07, "ohbu...@amber.plala.or.jp"

propeller

unread,

May 2, 2014, 1:37:09 AM5/2/14

to open...@googlegroups.com

円筒座標系に変換するする際には、以下の式からとあります。

ｚ軸を円筒中心軸に合わせた場合、速度成分の変換は
　　Vr=Vx*cosθ＋Vy*sinθ
　　Vθ=Vy*cosθ−Vx*sinθ
　　Vz=Vz

ParaViewで、上記のように径方向，周方向，軸方向の円筒座標系を表示するためには

どのように行えば良いのでしょうか？

ohbuchi

unread,

May 4, 2014, 7:21:36 AM5/4/14

to open...@googlegroups.com

こんにちは。

paraviewのCalculatorフィルタ機能を使って、半径方向速度と周方向速度を

定義すればよいと思います。

半径方向速度

　 U_r = (U_X * coordsX + U_Y * coordsY) / sqrt(coordsX^2+coordsY^2)

周方向速度

　 U_t = (U_Y * coordsX - U_X * coordsY) / sqrt(coordsX^2+coordsY^2)

2014年5月2日金曜日 7時37分09秒 UTC+2 propeller:

Reply all

Reply to author

Forward