解析面にかかる圧力の算出について

[nakayama yuu]

お世話になっております。

OpenFOAM初心者のnakayamaともうします。

OpenFOAMで器の中を流れる流体の解析をしているのですが、

その器の表面にかかる圧力分布を調べたいと思っています。

器はstlファイルで定義したので三角形メッシュです。

その一つ一つの三角形にかかる圧力を求めれたらと考えています。

経験がある方、もしくは方法がわかる方、ご教授いただけたらと思います。

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

[Fumiya Nozaki]

器の表面が、流体計算の境界になっていると思いますので、

p ファイル内に出力されている器の境界の結果を見れば、

境界上のフェイスの中心点での圧力値が分かります。

ご質問の意図は、

境界フェイスの中心点で出力された圧力値を

stl データの三角パッチの頂点 or 重心等に補間した値を

求めたいということでしょうか？

野崎

2016年1月15日 21:51 nakayama yuu <taaaaa...@gmail.com>:

--
このメールは Google グループのグループ「OpenFOAM」に登録しているユーザーに送られています。
このグループから退会し、グループからのメールの配信を停止するには openfoam+u...@googlegroups.com にメールを送信してください。
このグループに投稿するには open...@googlegroups.com にメールを送信してください。
https://groups.google.com/group/openfoam からこのグループにアクセスしてください。
その他のオプションについては https://groups.google.com/d/optout にアクセスしてください。

[nakayama yuu]

返信ありがとうございます。

pファイルの中身は解析領域内の各cellの圧力値であると認識していたのすが、

それで正解でしょうか？

またそうであった場合、どのような順番で並んでいるのかがまだ理解できていません。

ドキュメント等に書かれている事柄なのでしょうか？

僕がやりたいと思っていることは、

stlデータの各三角パッチにかかる圧力を算出し、ファイル等でoutputできればと考えています。

そのために、三角パッチに対応しているcell(cell数はそれぞれ違うと思いますが)の圧力値の面平均などを

出すのかなと漠然と考えていました。

しかし、野崎さんがおっしゃるように重心に補間した値でも同様のことかなと思います。

そのような方法はあるのでしょうか？

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

2016年1月16日土曜日 10時02分36秒 UTC+9 nozaki:

[Fumiya Nozaki]

圧力 p (タイプ: volScalarField) ファイルには、

・internalField の部分に、セル中心値が、

・boundaryField の部分に、境界ごとにその境界上のフェイス中心値が

それぞれ出力されています。

ファイルを開いて確認してみてください。

セルとフェイスの並び順ですが、ParaView で視覚的に確認できます。

私のスライドの 61ページからをご覧ください。

http://www.slideshare.net/fumiyanozaki96/boundary-conditions-in-openfoam

また、sample ユーティリティを使用すれば、

指定した境界上のフェイス中心座標と圧力値をファイルに出力することも可能です。

***** 出力の例 *****

# p  FACE_DATA 223

#  x  y  z  p

-0.0198086 0.0254 0 -2.32339

-0.0182574 0.0254 0 -3.96789

-0.0167681 0.0254 0 -4.94179

-0.0153384 0.0254 0 -5.22516

-0.0139658 0.0254 0 -5.35275

sample ユーティリティの使用方法は、ネットにたくさんの情報がありますので、

検索してみてください。

野崎

2016年1月16日 12:13 nakayama yuu <taaaaa...@gmail.com>: