blockMeshでの軸対称モデル作成時のエラー

[yug]

あけましておめでとうございます。
OpenFOAM初心者のyugです。

件名の通りblockMeshで軸対称モデルの作製を検討しているのですが、
checkMeshを行った際にエラーが出てしまいました。

エラーは軸対称モデルの軸部(axis)で起こっているようなのですが、
このエラーの対処法についてアドバイス等頂ければと思います。

添付致しました画像のように、モデルは円柱を軸対称にしたシンプルなものです。

blockMeshDictの中身、blockMesh時のlog、checkMesh時のlogを以下に載せます。

●blockMeshDict

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
convertToMeters 1;

vertices
(
    (-0.503 0 0)//0
    (-0.493 0 0)//1
    (-0.503 0.0029971446647521 -0.000130858161946576)//2
    (-0.493 0.0029971446647521 -0.000130858161946576)//3
    (-0.503 0.0029971446647521 0.000130858161946576)//4
    (-0.493 0.0029971446647521 0.000130858161946576)//5
);

blocks
(
    hex (0 1 3 2 0 1 5 4) (100 30 1) simpleGrading (1 1 1)
);

edges
(
);

boundary
(
    wall
    {
        type wall;
        faces
        (
        (2 4 5 3)
        );
    }
    outlet
    {
        type patch;
        faces
        (
        (1 3 5 1)
        );
    }
    axis
    {
        type symmetryPlane;
        faces
        (
        (0 1 1 0)
        );
    }
    inlet
    {
        type patch;
        faces
        (
        (0 0 4 2)
        );
    }
    back
    {
        type wedge;
        faces
        (
        (0 2 3 1)
        );
    }
    front
    {
        type wedge;
        faces
        (
        (0 1 5 4)
        );
    }
);

mergePatchPairs
(
);

// ************************************************************************* //



●blockMeshのlog

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

Creating block mesh from
    "/home/(ユーザ名)/OpenFOAM/(ユーザ名)-2.2.2/run/tutorials/incompressible/pimpleFoam/wedge-mizo7/constant/polyMesh/blockMeshDict"
Creating curved edges
Creating topology blocks
Creating topology patches

Creating block mesh topology

Check topology

    Basic statistics
        Number of internal faces : 0
        Number of boundary faces : 6
        Number of defined boundary faces : 6
        Number of undefined boundary faces : 0
    Checking patch -> block consistency

Creating block offsets
Creating merge list .

Creating polyMesh from blockMesh
Creating patches
Creating cells
Creating points with scale 1

Writing polyMesh
----------------
Mesh Information
----------------
  boundingBox: (-0.503 0 -0.000130858) (-0.493 0.00299714 0.000130858)
  nPoints: 6181
  nCells: 3000
  nFaces: 12065
  nInternalFaces: 5870
----------------
Patches
----------------
  patch 0 (start: 5870 size: 100) name: wall
  patch 1 (start: 5970 size: 30) name: outlet
  patch 2 (start: 6000 size: 35) name: axis
  patch 3 (start: 6035 size: 30) name: inlet
  patch 4 (start: 6065 size: 3000) name: back
  patch 5 (start: 9065 size: 3000) name: front

End


●checkMeshのlog

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

Create polyMesh for time = 0

Time = 0

Mesh stats
    points:           6181
    internal points:  0
    faces:            12065
    internal faces:   5870
    cells:            3000
    faces per cell:   5.97833
    boundary patches: 6
    point zones:      0
    face zones:       0
    cell zones:       0

Overall number of cells of each type:
    hexahedra:     2905
    prisms:        65
    wedges:        30
    pyramids:      0
    tet wedges:    0
    tetrahedra:    0
    polyhedra:     0

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                 
    wall                100      202      ok (non-closed singly connected) 
    outlet              30       61       ok (non-closed singly connected) 
    axis                35       68       ok (non-closed singly connected) 
    inlet               30       61       ok (non-closed singly connected) 
    back                3000     3131     ok (non-closed singly connected) 
    front               3000     3131     ok (non-closed singly connected) 

Checking geometry...
    Overall domain bounding box (-0.503 0 -0.000130858) (-0.493 0.00299714 0.000130858)
    Mesh (non-empty, non-wedge) directions (1 1 0)
    Mesh (non-empty) directions (1 1 1)
    Wedge back with angle 2.5 degrees
    Wedge front with angle 2.5 degrees
    All edges aligned with or perpendicular to non-empty directions.
    Boundary openness (1.37724e-18 1.78173e-15 2.29178e-16) OK.
    Max cell openness = 1.65594e-16 OK.
    Max aspect ratio = 2.00191 OK.
 ***Zero or negative face area detected.  Minimum area: 0
  <<Writing 35 zero area faces to set zeroAreaFaces
    Min volume = 4.35779e-14. Max volume = 2.57109e-12.  Total volume = 3.92201e-09.  Cell volumes OK.
    Mesh non-orthogonality Max: 0 average: 0
    Non-orthogonality check OK.
    Face pyramids OK.
 ***Max skewness = 1.59823e+12, 35 highly skew faces detected which may impair the quality of the results
  <<Writing 35 skew faces to set skewFaces
    Coupled point location match (average 0) OK.

Failed 2 mesh checks.

End



blockMeshDictにて分割数を(100 30 1)としているところを(10 30 1)と粗くすると2つのエラーはなくなるのですが、
もう少し細かいメッシュにしたいため困っております。

またblockMeshのlogを確認すると、patchでaxisのsizeが35となっており、
本来ここは0であるべきではないのかと思うのですが、
どう修正すべきか悩んでおります。

因みにOpenFOAMのバージョンは2.2.2で、
Ubuntu13.10にて使用しております。

皆様お忙しい中大変恐縮ですが、アドバイス等頂けたらと思います。

[ohbuchi]

こんにちは。
OpenFOAM-2.2以降で何か問題がある様です。ちなみにOpenFOAM-2.1では正常にメッシュが作成でき、
checkMeshのエラーもありませんでした。
Mantisを見ても該当するバグレポートはありません。
CFD-Onlineに同様の報告がありますが、結論に至っていません。
http://www.cfd-online.com/Forums/openfoam/117476-blockmesh-axisymmetry-error.html

2Dメッシュから軸対称メッシュを作成するツールmakeAxialMeshというものがあります。
https://openfoamwiki.net/index.php/Contrib/MakeAxialMesh
こちらを試したらどうでしょうか？

また、OpenFOAM Foundationにバグ報告をお願いします。
http://www.openfoam.org/bugs/



2015年1月5日月曜日 15時56分59秒 UTC+9 yug:

[ohbuchi]

補足です。OpenFOAM-2.2以降でもチュートリアルmovingConeなどwedgeを含むメッシュを正常に
作成できています。今回のものと何が違うのか、何度も見比べましたが違いを見つけることができませんでした。
したがって、単純なバグではない様です。なんとも謎な現象です。


2015年1月6日火曜日 10時55分27秒 UTC+9 ohbuchi:

[yug]

ohbuchi様

お世話になっております。
yugです。
早速のご返信、誠にありがとうございます。

OpenFOAMのバージョンやmakeAxialMesh等の貴重なご意見、
是非参考にさせて頂きたいと思います。

進展がありましたら再度こちらに投稿させて頂きたいと思います。