snappyHexMeshでのエラーについて

888 views
Skip to first unread message

Sakuma

unread,
Jun 12, 2010, 8:31:06 PM6/12/10
to OpenFOAM
いつも色々とお世話になっております。初心者のsakumaです。

snappyHexMeshでモデルを作成していますが,メッシュ作成後にcheckMeshを実施するとエラーが検出されます。
Time=0,1はMesh OKとなりますが,
Time=2で,
Checking topologyは全部OK
Checking patch topology for multiply connected surfacesは全部OK
Checking geometryは次のよう診断されました。
Overall domain bounding box (-0.13 -0.26 -0.13) (0.13 0.26 0.13)
Mesh (non-empty, non-wedge) directions (0 0 0)
Mesh (non-empty) directions (0 0 0)
***Number of edges not aligned with or perpendicular to non-empty
directions: 888720
<<Writing 530065 points on non-aligned edges to set nonAlignedEdges
Boundary openness (-3.331998e-16 -1.437487e-16 -1.435692e-15) OK.
Max cell openness = 3.1175e-16 OK.
Max aspect ratio = 9.780169 OK.
Minumum face area = 1.28114e-08. Maximum face area = 0.0004. Face
area magnitudes OK.
Min volume = 4.936347e-12. Max volume = 8e-06. Total volume =
0.03507671. Cell volumes OK.
Mesh non-orthogonality Max: 64.74152 average: 10.07508
Non-orthogonality check OK.
Face pyramids OK.
Max skewness = 3.64128 OK.
Failed 1 mesh checks.

Time=3も,
Checking topologyは全部OK
Checking patch topology for multiply connected surfacesは全部OK
Checking geometryはTime=2と同じようになります。
Checking geometry...
Overall domain bounding box (-0.1328339 -0.26 -0.1328357)
(0.1328411 0.26 0.132838)
Mesh (non-empty, non-wedge) directions (0 0 0)
Mesh (non-empty) directions (0 0 0)
***Number of edges not aligned with or perpendicular to non-empty
directions: 2800727
<<Writing 1141262 points on non-aligned edges to set nonAlignedEdges
Boundary openness (-5.588483e-16 -4.82978e-16 -1.532081e-15) OK.
Max cell openness = 9.149178e-16 OK.
Max aspect ratio = 33.75798 OK.
Minumum face area = 2.398806e-09. Maximum face area =
0.0004012592. Face area magnitudes OK.
Min volume = 4.247982e-13. Max volume = 8.001721e-06. Total
volume = 0.03593525. Cell volumes OK.
Mesh non-orthogonality Max: 64.99829 average: 10.02507
Non-orthogonality check OK.
Face pyramids OK.
Max skewness = 2.324375 OK.
Failed 1 mesh checks.

グーグルで”Number of edges not aligned with or perpendicular to non-empty
directions”と打つとcfd-onlineのスレッドに当たりますが,良くわかりません。何方か似たようなエラーが出て経験のある方がおられ
ましたら,対応方法を教えてください。よろしくお願いします。

E.Mogura

unread,
Jun 15, 2010, 5:27:39 PM6/15/10
to OpenFOAM
E.Mogura です

どなたからもコメントがないようなので、口火を切らせていただきます。

Sakumaさんの経験されたエラー
***Number of edges not aligned with or perpendicular to non-empty
であったかどうか記憶は定かでありませんが、
snappyHexMeshで作成したメッシュをcheckMeshにかけてエラーになる、
という経験は、現実の製品CADなど複雑形状を使った場合、毎度のことです。

MaxSkewのエラーなんかは、Dictパラメタで調整でき、これを小さくして
収束しやすくする(そのかわりメッシュの作成時間は長くなる)
といった経験はありますが、

それ以外は、せいぜいエラー箇所を表示確認する程度。
よほどたくさんの箇所で、主流領域にも存在。
それが原因で計算が動かない(発散する)可能性が大
ということでもあればともかく、計算が動いてくれればそのまま流しています。

経験的には、この種のエラーが原因で発散した例は非常に少ないです。
初めてのケースで発散した時点では、まずこれを疑うのですが、
ほとんどの場合、境界条件や乱流条件に問題ありの場合が多いです。

以上、反例(論)もあるかとは思いますが・・・

Takuya OSHIMA

unread,
Jun 16, 2010, 7:42:41 AM6/16/10
to OpenFOAM
大嶋です。

今回のエラーに関しては、エラーメッセージを意訳すると
「emptyパッチを含む場合は、メッシュは2次元メッシュの押出しでなければならない」
と言っているのではないかと思います。
ただ、non-empty directionsが(0 0 0) となっていたりする辺り、謎です。
私はsHMの経験はあまり無いので、違っていたらご容赦を。

Sakuma

unread,
Jun 17, 2010, 9:12:43 AM6/17/10
to OpenFOAM
 いつも色々とありがとうごさいます。sakumaです。

 元のCADモデルを調べたり,surfaceCheckでstlファイルを調べたり,cfd-onlineのスレッドを見たり,
頂いたアドバイスを参考とさせて頂いたりして原因を調査しました。でも決定的なエラーがどこにあるか
上手く見つけられませんでしたが,snappyHexMeshDictのrefinementSurfaceのpatch名のつづりが
一文字ミスタイプという全く間抜けなミスをしていました。このpatch名を直したら,checkMeshは全て
Mesh OKになりました。

 Ubuntu 8.04 64ビットの上にOpenFOAM1.6.xを入れてコマンドベースで使っていますが,やはりこの
ようななポカミスを無くするには,やはりDEXCSのようなランチャーが必要と実感しています。

 DEXCS 64ビット版が公開されるのは,いつごろでしょうか。
Reply all
Reply to author
Forward
0 new messages