[blockMesh] error message with modeling a cube with a hold at the center

hsingtzu · March 13, 2012, 20:44

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

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.

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.

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
Power placement question on a cube	Octagon1	FloEFD, FloWorks & FloTHERM	1	December 24, 2016 02:30
Data Center Modeling	Aaron Chen	STAR-CCM+	0	February 17, 2016 18:18
Center of a mesh cell	pyroknife	Main CFD Forum	1	July 10, 2015 06:25
Cooling cube in a room	John_Major	System Analysis	0	July 4, 2015 04:24
Gambit help: Cube inside cube	Jack Martinez	FLUENT	13	August 11, 2010 07:29

March 14, 2012, 10:56		#3
hsingtzu New Member Hsingtzu Wu Join Date: Feb 2011 Posts: 24 Rep Power: 15	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.