初心者です。メッシュの品質とその管理方法を教えてください。

480 views
Skip to first unread message

potpot

unread,
Apr 10, 2013, 4:44:54 AM4/10/13
to open...@googlegroups.com

初心者です。

計算の前準備として、メッシュの品質を行うことは大事なステップであると認識しています。

DEXCS2012について、下記のことを教えていただければ幸いです。

①メッシュの品質の善し悪しはどのように定義されているのでしょうか?

②メッシュの品質を確認する為に下記が出来るのであれば幸いです。

 1)メッシュの全セルについて、品質の分布がヒストグラムで表示できる。

 2)メッシュの全セルについて、品質のランクごとの色分けを施し、当該モデルケースで

   表示できる。

 3)品質の最も悪いセルをモデルケース上で表示できる。

これまでの経緯をよくわかっていないので、重複している質問かもしれませんが、

よろしくお願いします。

                                          以上


E.Mogura

unread,
Apr 10, 2013, 6:47:47 AM4/10/13
to open...@googlegroups.com
E.Mogura です

DEXCS2012についての質問ということで回答しますが、DEXCSそのものは何らメッシュの品質チェック機能は有しておりません。
 ⇒ OpenFOAMの機能に関する質問ということで、回答させていただきます。
(DEXCSとOpenFOAMの関係を切り分けて質問することをお勧めします)

OpenFOAMには、checkMeshというユーティリティが付属しており、これを使えばある程度の品質チェックは出来ますが、
②-1)
②-2)
といったところまではやってくれません。
②-3)
については、上記ユーティリティによって、エラーと判定されたメッシュについて、paraFoamを使って確認することは出来ると思います。

DEXCS2012に搭載の標準モデル(DEXCSフォント回りの流れ計算)について、上記checkMeshユーティリティを使うと、以下のような出力がされるはずです。
ご質問の「品質の善し悪しの定義」は、この出力メッセージの最後のあたり、Checking geometry...以下で、ある程度判断できるとは思います。
正確な定義方法については。。。ソースコードを見るか、どこかに解説があるかもしれませんが、ご存知の方におまかせします。

et@et-FX ~/Desktop/test/exe $ checkMesh    ←(メッシュの存在するケースフォルダにて実行) 
/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.1.x                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
Build  : 2.1.x-6502e850de67
Exec   : checkMesh
Date   : Apr 10 2013
Time   : 19:16:18
Host   : "et-FX"
PID    : 15124
Case   : /home/et/Desktop/test/exe
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Disallowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create polyMesh for time = 0

Time = 0

Mesh stats
    points:           108953
    faces:            278371
    internal faces:   262270
    cells:            86099
    boundary patches: 4
    point zones:      0
    face zones:       0
    cell zones:       0

Overall number of cells of each type:
    hexahedra:     65624
    prisms:        2941
    wedges:        408
    pyramids:      0
    tet wedges:    636
    tetrahedra:    6
    polyhedra:     16484

Checking topology...
    Boundary definition OK.
    Cell to face addressing OK.
    Point usage OK.
    Upper triangular ordering OK.
    Face vertices OK.
    Number of regions: 1 (OK).

Checking patch topology for multiply connected surfaces ...
    Patch               Faces    Points   Surface topology                  
    inlet               196      225      ok (non-closed singly connected)  
    outlet              196      225      ok (non-closed singly connected)  
    wall                1624     1680     ok (non-closed singly connected)  
    Dexcs_ExportedfromBlender-2.64(sub0)14085    17005    ok (closed singly connected)      

Checking geometry...
    Overall domain bounding box (-7 -5.50001 -1.80002) (1 -1.5 2.20003)
    Mesh (non-empty, non-wedge) directions (1 1 1)
    Mesh (non-empty) directions (1 1 1)
    Boundary openness (-3.64703e-18 -4.38243e-17 -4.59883e-17) OK.
    Max cell openness = 3.69214e-16 OK.
    Max aspect ratio = 8.7879 OK.
    Minimum face area = 6.49788e-06. Maximum face area = 0.082251.  Face area magnitudes OK.
    Min volume = 7.70405e-08. Max volume = 0.0226516.  Total volume = 127.892.  Cell volumes OK.
    Mesh non-orthogonality Max: 63.7271 average: 12.8883
    Non-orthogonality check OK.
    Face pyramids OK.
    Max skewness = 3.21216 OK.
    Coupled point location match (average 0) OK.

Mesh OK.

End




2013年4月10日 17:44 potpot <atom...@s6.dion.ne.jp>:


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

potpot

unread,
Apr 10, 2013, 8:55:58 PM4/10/13
to open...@googlegroups.com

E.Mogura 様

お世話になります。

ご教示ありがとうございます。

勉強させていただきます。

                以上



2013年4月10日水曜日 17時44分54秒 UTC+9 potpot:
Reply all
Reply to author
Forward
0 new messages