CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Native Meshers: blockMesh

error message with modeling a cube with a hold at the center

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Display Modes
Old   March 13, 2012, 20:44
Default error message with modeling a cube with a hold at the center
  #1
New Member
 
Hsingtzu Wu
Join Date: Feb 2011
Posts: 21
Rep Power: 6
hsingtzu is on a distinguished road
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 is offline   Reply With Quote

Old   March 14, 2012, 09:54
Default
  #2
New Member
 
Hsingtzu Wu
Join Date: Feb 2011
Posts: 21
Rep Power: 6
hsingtzu is on a distinguished road
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.
Attached Images
File Type: jpg problem.jpg (63.6 KB, 110 views)
hsingtzu is offline   Reply With Quote

Old   March 14, 2012, 10:56
Default
  #3
New Member
 
Hsingtzu Wu
Join Date: Feb 2011
Posts: 21
Rep Power: 6
hsingtzu is on a distinguished road
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.
hsingtzu is offline   Reply With Quote

Reply

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Explosion inside a cube - blockMesh shangzung OpenFOAM 3 March 6, 2013 04:59
Gambit help: Cube inside cube Jack Martinez FLUENT 13 August 11, 2010 06:29
Modeling Flow/Saturation/Absorption in Fibers Gene Dougherty Main CFD Forum 0 June 6, 2003 14:49


All times are GMT -4. The time now is 15:07.