CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM Native Meshers: blockMesh (http://www.cfd-online.com/Forums/openfoam-meshing-blockmesh/)
-   -   error message with modeling a cube with a hold at the center (http://www.cfd-online.com/Forums/openfoam-meshing-blockmesh/98563-error-message-modeling-cube-hold-center.html)

hsingtzu March 13, 2012 20:44

error message with modeling a cube with a hold at the center
 
my blockMeshDict code was like

Code:

/*--------------------------------*- C++ -*----------------------------------*\
| =========                |                                                |
| \\      /  F ield        | OpenFOAM: The Open Source CFD Toolbox          |
|  \\    /  O peration    | Version:  2.0.1                                |
|  \\  /    A nd          | Web:      www.OpenFOAM.com                      |
|    \\/    M anipulation  |                                                |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version    2.0;
    format      ascii;
    class      dictionary;
    object      blockMeshDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

convertToMeters 0.01;

vertices       
(
(0.45 0 0)  //21
(0.5 0 0)
(0.63 0 0)
(0.63 0.353553391 0)
(0.63 0.63 0)
(0.353553391 0.63 0)
(0.353553391 0.353553391 0)
(0.318198052 0.318198052 0)
(0 0.45 0)
(0 0.5 0)
(0 0.63 0)
(-0.353553391 0.63 0)
(-0.63 0.63 0)
(-0.63 0.353553391 0)
(-0.353553391 0.353553391 0)
(-0.318198052 0.318198052 0)
(-0.45 0 0)
(-0.5 0 0)
(-0.63 0 0)
(-0.63 -0.353553391 0)
(-0.63 -0.63 0)
(-0.353553391 -0.63 0)
(-0.353553391 -0.353553391 0)
(-0.318198052 -0.318198052 0)
(0 -0.45 0)
(0 -0.5 0)
(0 -0.63 0)
(0.353553391 -0.63 0)
(0.63 -0.63 0)
(0.63 -0.353553391 0)
(0.353553391 -0.353553391 0)
(0.318198052 -0.318198052 0)
(0.45 0 0.98)
(0.5 0 0.98)
(0.63 0 0.98)
(0.63 0.353553391 0.98)
(0.63 0.63 0.98)
(0.353553391 0.63 0.98)
(0.353553391 0.353553391 0.98)
(0.318198052 0.318198052 0.98)
(0 0.45 0.98)
(0 0.5 0.98)
(0 0.63 0.98)
(-0.353553391 0.63 0.98)
(-0.63 0.63 0.98)
(-0.63 0.353553391 0.98)
(-0.353553391 0.353553391 0.98)
(-0.318198052 0.318198052 0.98)
(-0.45 0 0.98)
(-0.5 0 0.98)
(-0.63 0 0.98)
(-0.63 -0.353553391 0.98)
(-0.63 -0.63 0.98)
(-0.353553391 -0.63 0.98)
(-0.353553391 -0.353553391 0.98)
(-0.318198052 -0.318198052 0.98)
(0 -0.45 0.98)
(0 -0.5 0.98)
(0 -0.63 0.98)
(0.353553391 -0.63 0.98)
(0.63 -0.63 0.98)
(0.63 -0.353553391 0.98)
(0.353553391 -0.353553391 0.98)
(0.318198052 -0.318198052 0.98)
   
);

blocks         
(
    hex (31 30 25 24 63 62 57 56) (10 10 5) simpleGrading (1 1 1)
    hex (0 1 30 31 32 33 62 63) (10 10 5) simpleGrading (1 1 1)
    hex (7 6 1 0 39 38 33 32) (10 10 5) simpleGrading (1 1 1)
    hex (8 9 6 7 40 41 38 39) (10 10 5) simpleGrading (1 1 1)
    hex (14 15 8 9 46 47 40 41) (10 10 5) simpleGrading (1 1 1)
    hex (17 16 15 14 49 48 47 46) (10 10 5) simpleGrading (1 1 1)
    hex (22 23 16 17 54 55 48 49) (10 10 5) simpleGrading (1 1 1)
    hex (25 24 23 22 57 56 55 54) (10 10 5) simpleGrading (1 1 1)
    hex (25 26 27 30 57 58 59 62) (10 10 5) simpleGrading (1 1 1)
    hex (27 28 29 30 59 60 61 62) (10 10 5) simpleGrading (1 1 1)
    hex (29 30 1 2 61 62 33 34) (10 10 5) simpleGrading (1 1 1)
    hex (1 2 3 6 33 34 35 38) (10 10 5) simpleGrading (1 1 1)
    hex (3 4 5 6 35 36 37 38) (10 10 5) simpleGrading (1 1 1)
    hex (5 6 9 10 37 38 41 42) (10 10 5) simpleGrading (1 1 1)
    hex (9 10 11 14 41 42 43 46) (10 10 5) simpleGrading (1 1 1)
    hex (11 12 13 14 43 44 45 46) (10 10 5) simpleGrading (1 1 1)
    hex (13 14 17 18 45 46 49 50) (10 10 5) simpleGrading (1 1 1)
    hex (17 18 19 22 49 50 51 54) (10 10 5) simpleGrading (1 1 1)
    hex (19 20 21 22 51 52 53 54) (10 10 5) simpleGrading (1 1 1)
    hex (21 22 25 26 53 54 57 58) (10 10 5) simpleGrading (1 1 1)
);

edges         
(
    arc 0 7 (0.4 0.2062 0)  //114
    arc 7 8 (0.2062 0.4 0)
    arc 8 15 (-0.2062 0.4 0)
    arc 15 16 (-0.4 0.2062 0)
    arc 16 23 (-0.4 -0.2062 0)
    arc 23 24 (-0.2062 -0.4 0)
    arc 24 31 (0.2062 -0.4 0)
    arc 31 0 (0.4 -0.2062 0)
    arc 32 39 (0.4 0.2062 0.98)
    arc 39 40 (0.2062 0.4 0.98)
    arc 40 47 (-0.2062 0.4 0.98)
    arc 47 48 (-0.4 0.2062 0.98)
    arc 48 55 (-0.4 -0.2062 0.98)
    arc 55 56 (-0.2062 -0.4 0.98)
    arc 56 63 (0.2062 -0.4 0.98)
    arc 63 32 (0.4 -0.2062 0.98)
arc 1 6 (0.4 0.3 0)
arc 6 9 (0.3 0.4 0)
arc 9 14 (-0.3 0.4 0)
arc 14 17 (-0.4 0.3 0)
arc 17 22 (-0.4 -0.3 0)
arc 22 25 (-0.3 -0.4 0)
arc 25 30 (0.3 -0.4 0)
arc 30 1 (0.4 -0.3 0)
arc 33 38 (0.4 0.3 0.98)
arc 38 41 (0.3 0.4 0.98)
arc 41 46 (-0.3 0.4 0.98)
arc 46 49 (-0.4 0.3 0.98)
arc 49 54 (-0.4 -0.3 0.98)
arc 54 57 (-0.3 -0.4 0.98)
arc 57 62 (0.3 -0.4 0.98)
arc 62 33 (0.4 -0.3 0.98)
);

boundary
(
    left
    {
//        type symmetryPlane;
        type cyclic;
        neighbourPatch right;
        faces
        (
            (12 13 45 44)
            (13 18 50 45)
            (18 19 51 50)
            (19 20 52 51)
        );
    }
    right
    {
//        type symmetryPlane;
        type cyclic;
        neighbourPatch left;
        faces
        (
            (3 4 36 35)
            (3 2 34 35)
            (2 29 61 34)
            (29 28 60 61)
        );
    }
    down
    {
//        type symmetryPlane;
        type cyclic;
        neighbourPatch up;
        faces
        (
            (20 21 53 52)
            (21 26 58 53)
            (26 27 59 58)
            (27 28 60 59)
        );
    }
    up
    {
//        type symmetryPlane;
                  type cyclic;
        neighbourPatch down;
        faces
        (
            (11 12 44 43)
            (11 10 42 43)
            (10 5 37 42)
            (4 5 37 36)
        );
    }
    holeDownLeft
    {
        type patch;
        faces
        (
            (16 23 55 48)
            (23 24 56 55)

        );
    }
    holeDownRight
    {
        type patch;
        faces
        (
            (24 31 63 56)
            (0 31 63 32)
        );
    }
    holeUpLeft
    {
        type patch;
        faces
        (
            (15 16 48 47)
            (8 15 47 40)
        );
    }
    holeUpRight
    {
        type patch;
        faces
        (
            (0 7 39 32)
            (7 8 40 39)
        );
    }
    front
    {
        type patch;
        faces
        (
            (24 25 30 31)
            (30 31 0 1)
            (0 1 6 7)
            (6 7 8 9)
            (8 9 14 15)
            (14 15 16 17)
            (16 17 22 23)
            (22 23 24 25)
            (25 26 27 30)
            (27 28 29 30)
            (29 30 1 2)
            (1 2 3 6)
            (3 4 5 6)
            (5 6 9 10)
            (9 10 11 14)
            (11 12 13 14)
            (13 14 17 18)
            (17 18 19 22)
            (19 20 21 22)
            (21 22 25 26)
        );
    }
    back
    {
        type patch;
        faces
        (
            (56 57 62 63)
            (62    63 32 33)
            (32    33 38 39)
(38    39    40    41)
(40    41    46    47)
(46    47    48    49)
(48    49    54    55)
(54    55    56    57)
(57    58    59    62)
(59    60    61    62)
(61    62    33    34)
(33    34    35    38)
(35    36    37    38)
(37    38    41    42)
(41    42    43    46)
(43    44    45    46)
(45    46    49    50)
(49    50    51    54)
(51    52    53    54)
(53 54 57 58)
        );
    }
);

mergePatchPairs
(
);

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

when I ran "blockMesh", the model could be compiled, but I got some warnings.
Code:

Creating block mesh topology
--> FOAM Warning :
    From function cellModel::mag(const labelList&, const pointField&)
    in file meshes/meshShapes/cellModel/cellModel.C at line 128
    zero or negative pyramid volume: -0.00259862 for face 0
--> FOAM Warning :
    From function cellModel::mag(const labelList&, const pointField&)
    in file meshes/meshShapes/cellModel/cellModel.C at line 128
    zero or negative pyramid volume: -0.00288735 for face 1
--> FOAM Warning :
    From function cellModel::mag(const labelList&, const pointField&)
    in file meshes/meshShapes/cellModel/cellModel.C at line 128
    zero or negative pyramid volume: -0.00274299 for face 2
--> FOAM Warning :
    From function cellModel::mag(const labelList&, const pointField&)
    in file meshes/meshShapes/cellModel/cellModel.C at line 128
    zero or negative pyramid volume: -0.00274299 for face 3
--> FOAM Warning :
    From function cellModel::mag(const labelList&, const pointField&)
    in file meshes/meshShapes/cellModel/cellModel.C at line 128
    zero or negative pyramid volume: -0.00274299 for face 4
--> FOAM Warning :
    From function cellModel::mag(const labelList&, const pointField&)
    in file meshes/meshShapes/cellModel/cellModel.C at line 128
    zero or negative pyramid volume: -0.00274299 for face 5
--> FOAM Warning :
    From function blockMesh::createTopology(IOdictionary&)
    in file blockMesh/blockMeshTopology.C at line 255
    negative volume block : 0, probably defined inside-out
--> FOAM Warning :
    From function cellModel::mag(const labelList&, const pointField&)
    in file meshes/meshShapes/cellModel/cellModel.C at line 128
    zero or negative pyramid volume: -0.00259862 for face 0
--> FOAM Warning :
    From function cellModel::mag(const labelList&, const pointField&)
    in file meshes/meshShapes/cellModel/cellModel.C at line 128
    zero or negative pyramid volume: -0.00288735 for face 1
--> FOAM Warning :
    From function cellModel::mag(const labelList&, const pointField&)
    in file meshes/meshShapes/cellModel/cellModel.C at line 128
    zero or negative pyramid volume: -0.00274299 for face 2
--> FOAM Warning :
    From function cellModel::mag(const labelList&, const pointField&)
    in file meshes/meshShapes/cellModel/cellModel.C at line 128
    zero or negative pyramid volume: -0.00274299 for face 3
--> FOAM Warning :
    From function cellModel::mag(const labelList&, const pointField&)
    in file meshes/meshShapes/cellModel/cellModel.C at line 128
    zero or negative pyramid volume: -0.00274299 for face 4
--> FOAM Warning :
    From function cellModel::mag(const labelList&, const pointField&)
    in file meshes/meshShapes/cellModel/cellModel.C at line 128
    zero or negative pyramid volume: -0.00274299 for face 5
--> FOAM Warning :
    From function blockMesh::createTopology(IOdictionary&)
    in file blockMesh/blockMeshTopology.C at line 255
    negative volume block : 1, probably defined inside-out
--> FOAM Warning :
    From function cellModel::mag(const labelList&, const pointField&)
    in file meshes/meshShapes/cellModel/cellModel.C at line 128
    zero or negative pyramid volume: -0.00259862 for face 0
--> FOAM Warning :
    From function cellModel::mag(const labelList&, const pointField&)
    in file meshes/meshShapes/cellModel/cellModel.C at line 128
    zero or negative pyramid volume: -0.00288735 for face 1
--> FOAM Warning :
    From function cellModel::mag(const labelList&, const pointField&)
    in file meshes/meshShapes/cellModel/cellModel.C at line 128
    zero or negative pyramid volume: -0.00274299 for face 2
--> FOAM Warning :
    From function cellModel::mag(const labelList&, const pointField&)
    in file meshes/meshShapes/cellModel/cellModel.C at line 128
    zero or negative pyramid volume: -0.00274299 for face 3
--> FOAM Warning :
    From function cellModel::mag(const labelList&, const pointField&)
    in file meshes/meshShapes/cellModel/cellModel.C at line 128
    zero or negative pyramid volume: -0.00274299 for face 4
--> FOAM Warning :
    From function cellModel::mag(const labelList&, const pointField&)
    in file meshes/meshShapes/cellModel/cellModel.C at line 128
    zero or negative pyramid volume: -0.00274299 for face 5
--> FOAM Warning :
    From function blockMesh::createTopology(IOdictionary&)
    in file blockMesh/blockMeshTopology.C at line 255
    negative volume block : 2, probably defined inside-out
--> FOAM Warning :
    From function cellModel::mag(const labelList&, const pointField&)
    in file meshes/meshShapes/cellModel/cellModel.C at line 128
    zero or negative pyramid volume: -0.00259862 for face 0
--> FOAM Warning :
    From function cellModel::mag(const labelList&, const pointField&)
    in file meshes/meshShapes/cellModel/cellModel.C at line 128
    zero or negative pyramid volume: -0.00288735 for face 1
--> FOAM Warning :
    From function cellModel::mag(const labelList&, const pointField&)
    in file meshes/meshShapes/cellModel/cellModel.C at line 128
    zero or negative pyramid volume: -0.00274299 for face 2
--> FOAM Warning :
    From function cellModel::mag(const labelList&, const pointField&)
    in file meshes/meshShapes/cellModel/cellModel.C at line 128
    zero or negative pyramid volume: -0.00274299 for face 3
--> FOAM Warning :
    From function cellModel::mag(const labelList&, const pointField&)
    in file meshes/meshShapes/cellModel/cellModel.C at line 128
    zero or negative pyramid volume: -0.00274299 for face 4
--> FOAM Warning :
    From function cellModel::mag(const labelList&, const pointField&)
    in file meshes/meshShapes/cellModel/cellModel.C at line 128
    zero or negative pyramid volume: -0.00274299 for face 5
--> FOAM Warning :
    From function blockMesh::createTopology(IOdictionary&)
    in file blockMesh/blockMeshTopology.C at line 255
    negative volume block : 3, probably defined inside-out
--> FOAM Warning :
    From function cellModel::mag(const labelList&, const pointField&)
    in file meshes/meshShapes/cellModel/cellModel.C at line 128
    zero or negative pyramid volume: -0.0117355 for face 0
--> FOAM Warning :
    From function cellModel::mag(const labelList&, const pointField&)
    in file meshes/meshShapes/cellModel/cellModel.C at line 128
    zero or negative pyramid volume: -0.0117355 for face 1
--> FOAM Warning :
    From function cellModel::mag(const labelList&, const pointField&)
    in file meshes/meshShapes/cellModel/cellModel.C at line 128
    zero or negative pyramid volume: -0.015964 for face 2
--> FOAM Warning :
    From function cellModel::mag(const labelList&, const pointField&)
    in file meshes/meshShapes/cellModel/cellModel.C at line 128
    zero or negative pyramid volume: -0.00750712 for face 3
--> FOAM Warning :
    From function cellModel::mag(const labelList&, const pointField&)
    in file meshes/meshShapes/cellModel/cellModel.C at line 128
    zero or negative pyramid volume: -0.0117355 for face 4
--> FOAM Warning :
    From function cellModel::mag(const labelList&, const pointField&)
    in file meshes/meshShapes/cellModel/cellModel.C at line 128
    zero or negative pyramid volume: -0.0117355 for face 5
--> FOAM Warning :
    From function blockMesh::createTopology(IOdictionary&)
    in file blockMesh/blockMeshTopology.C at line 255
    negative volume block : 10, probably defined inside-out
--> FOAM Warning :
    From function cellModel::mag(const labelList&, const pointField&)
    in file meshes/meshShapes/cellModel/cellModel.C at line 128
    zero or negative pyramid volume: -0.0117355 for face 0
--> FOAM Warning :
    From function cellModel::mag(const labelList&, const pointField&)
    in file meshes/meshShapes/cellModel/cellModel.C at line 128
    zero or negative pyramid volume: -0.0117355 for face 1
--> FOAM Warning :
    From function cellModel::mag(const labelList&, const pointField&)
    in file meshes/meshShapes/cellModel/cellModel.C at line 128
    zero or negative pyramid volume: -0.015964 for face 2
--> FOAM Warning :
    From function cellModel::mag(const labelList&, const pointField&)
    in file meshes/meshShapes/cellModel/cellModel.C at line 128
    zero or negative pyramid volume: -0.00750712 for face 3
--> FOAM Warning :
    From function cellModel::mag(const labelList&, const pointField&)
    in file meshes/meshShapes/cellModel/cellModel.C at line 128
    zero or negative pyramid volume: -0.0117355 for face 4
--> FOAM Warning :
    From function cellModel::mag(const labelList&, const pointField&)
    in file meshes/meshShapes/cellModel/cellModel.C at line 128
    zero or negative pyramid volume: -0.0117355 for face 5
--> FOAM Warning :
    From function blockMesh::createTopology(IOdictionary&)
    in file blockMesh/blockMeshTopology.C at line 255
    negative volume block : 13, probably defined inside-out
--> FOAM Warning :
    From function cellModel::mag(const labelList&, const pointField&)
    in file meshes/meshShapes/cellModel/cellModel.C at line 128
    zero or negative pyramid volume: -0.0117355 for face 0
--> FOAM Warning :
    From function cellModel::mag(const labelList&, const pointField&)
    in file meshes/meshShapes/cellModel/cellModel.C at line 128
    zero or negative pyramid volume: -0.0117355 for face 1
--> FOAM Warning :
    From function cellModel::mag(const labelList&, const pointField&)
    in file meshes/meshShapes/cellModel/cellModel.C at line 128
    zero or negative pyramid volume: -0.015964 for face 2
--> FOAM Warning :
    From function cellModel::mag(const labelList&, const pointField&)
    in file meshes/meshShapes/cellModel/cellModel.C at line 128
    zero or negative pyramid volume: -0.00750712 for face 3
--> FOAM Warning :
    From function cellModel::mag(const labelList&, const pointField&)
    in file meshes/meshShapes/cellModel/cellModel.C at line 128
    zero or negative pyramid volume: -0.0117355 for face 4
--> FOAM Warning :
    From function cellModel::mag(const labelList&, const pointField&)
    in file meshes/meshShapes/cellModel/cellModel.C at line 128
    zero or negative pyramid volume: -0.0117355 for face 5
--> FOAM Warning :
    From function blockMesh::createTopology(IOdictionary&)
    in file blockMesh/blockMeshTopology.C at line 255
    negative volume block : 16, probably defined inside-out
--> FOAM Warning :
    From function cellModel::mag(const labelList&, const pointField&)
    in file meshes/meshShapes/cellModel/cellModel.C at line 128
    zero or negative pyramid volume: -0.0117355 for face 0
--> FOAM Warning :
    From function cellModel::mag(const labelList&, const pointField&)
    in file meshes/meshShapes/cellModel/cellModel.C at line 128
    zero or negative pyramid volume: -0.0117355 for face 1
--> FOAM Warning :
    From function cellModel::mag(const labelList&, const pointField&)
    in file meshes/meshShapes/cellModel/cellModel.C at line 128
    zero or negative pyramid volume: -0.015964 for face 2
--> FOAM Warning :
    From function cellModel::mag(const labelList&, const pointField&)
    in file meshes/meshShapes/cellModel/cellModel.C at line 128
    zero or negative pyramid volume: -0.00750712 for face 3
--> FOAM Warning :
    From function cellModel::mag(const labelList&, const pointField&)
    in file meshes/meshShapes/cellModel/cellModel.C at line 128
    zero or negative pyramid volume: -0.0117355 for face 4
--> FOAM Warning :
    From function cellModel::mag(const labelList&, const pointField&)
    in file meshes/meshShapes/cellModel/cellModel.C at line 128
    zero or negative pyramid volume: -0.0117355 for face 5
--> FOAM Warning :
    From function blockMesh::createTopology(IOdictionary&)
    in file blockMesh/blockMeshTopology.C at line 255
    negative volume block : 19, probably defined inside-out

and when I used "paraFoam" to view the model, it gave me the following error message
Code:

--> FOAM FATAL ERROR:
More than six unsigned transforms detected:
6(((0.0126 -0.00248802 0) (1 0 0 0 1 0 0 0 1) 0) ((0.0126 -0.00221157 0.00196) (1 0 0 0 1 0 0 0 1) 0) ((0.0126 -0.00193513 0.00392) (1 0 0 0 1 0 0 0 1) 0) ((0.0126 -0.00165868 0.00588) (1 0 0 0 1 0 0 0 1) 0) ((0.0126 -0.00138223 0.00784) (1 0 0 0 1 0 0 0 1) 0) ((0.0126 -0.00082934 0) (1 0 0 0 1 0 0 0 1) 0))

    From function void Foam::globalIndexAndTransform::determineTransforms()
    in file primitives/globalIndexAndTransform/globalIndexAndTransform.C at line 185.

FOAM exiting

I would appreciate any suggestion.
Hsingtzu

hsingtzu March 14, 2012 09:54

1 Attachment(s)
If I replaced "cyclic" with "symmetryPlane", then I could run paraFoam even though I still got a series of warnings. I looked into the mesh, and found something strange (see the attachment). I have no idea what might cause this problem.:confused:

hsingtzu March 14, 2012 10:56

I got it. After reviewing the previous posts, I wondered if I should number the blocks in the same order. So I did it, and it worked. There was no more strange mesh even though there were some more warning messages. :)


All times are GMT -4. The time now is 11:26.