snappyHexMeshでの形状再現性向上
435 views
Skip to first unread message
村井和彦
Sep 11, 2015, 10:50:16 PM9/11/15
to OpenFOAM
村井です。 こんにちは。
もう十分議論されてきた内容のように思えますが、
snappyHexMeshで形状が欠けてしまいます。
今回、DEXCSに搭載されているTreeFoamを使ってのメッシュ作成を
行いましたが、何故か私が作成すると、ボックスのような簡単な形状でも
欠けるようなメッシュができてしまいます。
何か間違っているかと思いますが、回避方法を教えて頂けないでしょうか？
メッシュ作成した際のcaseファイル（polyMeshは削除）と、
どんな不具合が起きているかの説明資料を送ります。
よろしくお願いします。
Masashi Imano
Sep 12, 2015, 2:05:59 AM9/12/15
to OpenFOAM
今野です．
applications/utilities/mesh/generation/snappyHexMesh/snappyHexMeshDict
constant/polyMesh/blockMeshDict が無いので，blockMesh が実行できず，
不具合の再現ができませんが， locationInMeshの座標値が切りの良い数字なので，
恐らく格子の界面か頂点に厳密に一致しており，内外判定に失敗しているのだと思います．
system/snappyHexMeshDict
locationInMesh (0.07 0 0);
に書いてある以下のコメントも参考にしてください．
// NOTE: This point should never be on a face, always inside a cell, even
// after refinement.
locationInMesh (5 0.28 0.43);
// NOTE: This point should never be on a face, always inside a cell, even
// after refinement.
locationInMesh (5 0.28 0.43);
https://github.com/OpenFOAM/OpenFOAM-2.4.x/blob/master/applications/utilities/mesh/generation/snappyHexMesh/snappyHexMeshDict
注)内部点はベースメッシュの界面や辺上に あってはならないので，形状ファイルの領 域の角の少し内側を指定するのが良い
2015年9月12日土曜日 11時50分16秒 UTC+9 村井和彦:
以下の講習会資料のpp.56も参考にしてください．
まお，ユーティリティのソースのディレクトリにあるDictファイルは，
通常コメントが詳しく書いてありますので，ユーティリティで問題が起きた場合には，
まずソースのディレクトリにあるDictファイルにあるコメントに目を通したほうが良いです．
それでも解決できない場合には，ソースファイル(*.C, *.H)のコメントに目を通す必要が
生じます．
2015年9月12日土曜日 11時50分16秒 UTC+9 村井和彦:
MURAI Kazuhiko
Sep 13, 2015, 1:29:23 AM9/13/15
to open...@googlegroups.com
村井です。 こんにちは。
アドバイス通り、locationInMesh を色々変えてやってみましたが、うまく行きません。
端っこにしてもそうですし、あまり端過ぎると、うまく切れませんでした。
（外側の部分だけが切れてしまう。）
ソースファイルも見てますが・・・決め手に欠いております。
もう少し見てみます。
--
このメールは Google グループのグループ「OpenFOAM」に登録しているユーザーに送られています。
このグループから退会し、グループからのメールの配信を停止するには openfoam+u...@googlegroups.com にメールを送信してください。
このグループに投稿するには open...@googlegroups.com にメールを送信してください。
http://groups.google.com/group/openfoam からこのグループにアクセスしてください。
その他のオプションについては https://groups.google.com/d/optout にアクセスしてください。
Masashi Imano
Sep 13, 2015, 3:55:10 AM9/13/15
to OpenFOAM
今野です．
constant/polyMesh/blockMeshDict を送ってください．
2015年9月13日日曜日 14時29分23秒 UTC+9 村井和彦:
村井です。 こんにちは。アドバイス通り、locationInMesh を色々変えてやってみましたが、うまく行きません。端っこにしてもそうですし、あまり端過ぎると、うまく切れませんでした。（外側の部分だけが切れてしまう。）ソースファイルも見てますが・・・決め手に欠いております。
もう少し見てみます。
このグループから退会し、グループからのメールの配信を停止するには openfoam+unsubscribe@googlegroups.com にメールを送信してください。
村井和彦
Sep 14, 2015, 6:41:49 AM9/14/15
to OpenFOAM
村井です。 こんばんは。
2015年9月13日日曜日 16時55分10秒 UTC+9 Masashi Imano:
遅くなりましたが、blockMeshDictを送付します。
あと、ご参考までにsnappyHexMeshDictとsurfaceFeatureExtractDictも送付します。
よろしくお願いします。
2015年9月13日日曜日 16時55分10秒 UTC+9 Masashi Imano:
もう少し見てみます。
このグループから退会し、グループからのメールの配信を停止するには openfoam+u...@googlegroups.com にメールを送信してください。
0 new messages