計算領域がうまく伝わらない

699 views
Skip to first unread message

ひげひげひげGN

unread,
May 1, 2021, 5:23:07 AM5/1/21
to OpenFOAM
お世話になります。

まだOpenFOAMは日が浅く、これが数値計算の最初です。
現在は、XSim を利用していろいろ試そうとしているところです。
いまは外部流れを計算しています。

--> FOAM FATAL IO ERROR: 
Unable to set reference cell for field p
    Please supply either pRefCell or pRefPoint

のエラーが出るのですが、

  出口で静圧ゼロ指定、

しているつもりのになぜかわからずいろいろと試していたところ、
計算領域がおかしいようであると気づきました。

XSim をもちいて、
計算領域を STL で指定して、
境界条件も形状名称を指定して設定し、
メッシュは3~4段の再分割設定しています。
XSimの形状表示では、各面が存在しているようで問題はなさそうでした。

ところが、エラー後、ParaVIEWでメッシュを確認したところ、
OpenFOAM で処理したときに、
一番大きな領域ではなく、
一番大きな再分割領域でしかメッシュが切れていません。
物体(複数試しました)の形状は、表示されていて問題はなさそうです。
そのため、四隅の角も丸くなまった領域になっています。
当然、境界条件を指定した各面は含まれていませんので、
そのために、前述のエラーが出ているような気がしました。

そもそも、
領域を指定しているのに、
それよりも小さな領域再分割領域でしかメッシュがきられていない原因は何でしょうか。
最初は、球や円周りの計算ができていましたので、
自覚なく、何か基本的な大きなミスをしているのでしょうか。
その原因を調べるために調べるべき点等はどこでしょうか。

とりあえず、
この問題を解決しなければ
圧力参照点の問題も解決しないように思っていますので、
なにかご助言、ご指導のほどよろしくお願いいたします。

質問、お読みいただいてありがとうございます。

haruka tsubota

unread,
May 6, 2021, 6:33:56 PM5/6/21
to OpenFOAM
質問の内容を読む限りでは メッシュ設定 の 「計算領域」を正しく設定できていない可能性が高いように思います。

1. XSimのメッシュ設定で「計算領域」はメッシュを切りたい閉空間(=流体が流れる領域)に正しく設定されているかを確認してください。

2. 「1.」が問題ない場合、再分割設定を全て外してメッシュを切った場合に正しくメッシュが切れるかを確認してください。
 これは 問題があるのが再分割設定か、それ以外の場所かの問題切り分けのためです。

3. 上記、「1.」「2.」を試して問題が解決しない場合、(可能であれば)XSimからダウンロードした、問題起きるZipファイルをここにアップロードしてください。
 確認してみます。


2021年5月1日土曜日 18:23:07 UTC+9 tomca...@gmail.com:

ひげひげひげGN

unread,
May 7, 2021, 12:20:58 PM5/7/21
to OpenFOAM
ご助言ありがとうございます。

1.につきましては、十分広い計算領域を設定できていると思います。
形状を追加、にて計算領域を設定した場合も、
計算領域を含めて stl でアップロードした場合も、
XSim の長方形空間の座標はモデルを包含しているようです。

ついでながら、XSim のメッシュの項の、計算領域、は何の大きさを表しているのでしょうか。

2.につきまして、計算領域を メッシュ作成後に確認したところ、
計算領域がメッシュ生成後は小さいことに気がつきました。
再分割領域のみがメッシュ精製されて、計算領域全体災害園部分がなくなっていました。
ParaVIEWで確認したところ、
二次元流れ場を計算しようとしていたので、
再分配領域の厚さが不足していたようで、
再分配領域でメッシュが一層しかありませんでした。
また、直方体の四隅が丸く欠けていました。

おそらく、計算領域の最も粗いメッシュが、
計算領域の厚さより大きかったために、無視されていたようです。

ご助言いただいたように、
再分割前の状態で領域が読み込まれるかを確認するのは、
一番大きなッシュの大きさが適切かどうかを確認するのによい方法とわかりました。
大変参考になりました。ありがとうございます。

圧力参照点が失われていたために、今回のエラーが出ていたということに納得ができました。
また、
おそらく同じ原因で、
計算はできましたが設定した境界条件とは全く異なり一様流が減衰していくおかしなケースもありましたが、
おそらく流入条件の領域が失われていたのかもしれず、
今回と同じ原因ではないかと納得できました。

XSim の画面ではなかなか厚さ方向のメッシュの確認が、
薄い計算領域の場合は難しいので、
最も粗いメッシュが適正稼働かをチェックするという切り分け方法として活用させていただきます。

初心者故の非常に初歩的なミスでした。
おかげさまで、先に進むことができそうです。
大変有用なご助言、大変ありがとうございました。
今後も何かありましたらよろしくお願い申し上げます。


2021年5月7日金曜日 7:33:56 UTC+9 haruka....@gmail.com:
Reply all
Reply to author
Forward
0 new messages