CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM (http://www.cfd-online.com/Forums/openfoam/)
-   -   BlockMeshDict for a Pipe (http://www.cfd-online.com/Forums/openfoam/98791-blockmeshdict-pipe.html)

sen.1986 March 19, 2012 11:34

BlockMeshDict for a Pipe
 
Hi!

I am trying to create a pipe, with inlet and outlet in x direction, so that I can avail the flow along x-direction.

I am getting the following error:

Quote:

/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.1.0 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Build : 2.1.0-0bc225064152
Exec : blockMesh
Date : Mar 19 2012
Time : 16:30:54
Host : "fwss116"
PID : 30357
Case : /home/sen04/OpenFOAM/sen04-2.1.0/testMesh
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

Creating block mesh from
"/home/sen04/OpenFOAM/sen04-2.1.0/testMesh/constant/polyMesh/blockMeshDict"
Creating curved edges
Creating topology blocks
Creating topology patches

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: -30 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: -30 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: -10 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: -50 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: -30 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: -30 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: -30 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: -30 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: -50 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: -10 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: -30 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: -30 for face 5
--> FOAM Warning :
From function blockMesh::createTopology(IOdictionary&)
in file blockMesh/blockMeshTopology.C at line 255
negative volume block : 4, probably defined inside-out

Check topology

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

Creating block offsets
Creating merge list

--> FOAM FATAL ERROR:
Inconsistent number of faces between block pair 0 and 1

From function blockMesh::calcMergeInfo()
in file blockMesh/blockMeshMerge.C at line 221.

FOAM exiting

Could someone tell me my possible mistake? I am also attaching my BlockMeshDict:

Quote:

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

convertToMeters 0.004;

vertices
(
/*0*/(0 2 2)
/*1*/(0 0 0)
/*2*/(0 0 5)
/*3*/(0 2 3)
/*4*/(0 3 3)
/*5*/(0 5 5)
/*6*/(0 5 0)
/*7*/(0 3 2)
/*8*/(30 2 2)
/*9*/(30 0 0)
/*10*/(30 0 5)
/*11*/(30 2 3)
/*12*/(30 3 3)
/*13*/(30 5 5)
/*14*/(30 5 0)
/*15*/(30 3 2)

);

blocks
(
hex (0 7 4 3 8 15 12 11) (80 11 11) simpleGrading (5 1 1)
hex (3 4 5 2 11 12 13 10) (80 11 51) simpleGrading (5 1 1)
hex (4 7 6 5 12 15 14 13) (80 51 11) simpleGrading (5 1 1)
hex (0 7 6 1 8 15 14 9) (80 11 51) simpleGrading (5 1 1)
hex (1 2 3 0 9 10 11 8) (80 51 11) simpleGrading (5 1 1)
);

edges
(
arc 5 6 (0 6.53553 2.5)
arc 6 1 (0 2.5 -6.53553)
arc 1 2 (0 -6.53553 2.5)
arc 2 5 (0 2.5 6.53553)
arc 5 6 (30 6.53553 2.5)
arc 6 1 (30 2.5 -6.53553)
arc 1 2 (30 -6.53553 2.5)
arc 2 5 (30 2.5 6.53553)
);

boundary
(
inlet
{
type cyclic;
neighbourPatch outlet;
faces
(
(1 2 3 0)
(3 2 5 4)
(4 5 6 7)
(7 6 1 0)
(0 3 4 7)
);
}

outlet
{
type cyclic;
neighbourPatch inlet;
faces
(
(9 8 11 10)
(11 12 13 10)
(15 14 13 12)
(8 9 14 15)
(15 12 11 8)
);
}

wall1
{
type wall;
faces (

(9 1 6 14)
(9 10 2 1)
(10 13 5 2)
(14 6 5 13)
);

}

);


mergePatchPairs
(
);

anon_a March 19, 2012 16:41

Code:

negative volume block
A simple search in the forum would give you

http://www.cfd-online.com/Forums/ope...-vertexes.html
http://www.cfd-online.com/Forums/ope...m-warning.html

for example. Your geometry is defined inside out. This can be alleviated by rearranging the sequence of vertices for the problematic hexes in the blocks section.

Code:

Inconsistent number of faces between block pair 0 and 1
You are probably trying to match in one edge 51 cells from one side with 11 cells from the other. You have probably misplaced one of the 51's. Always start with something simple uniform and then move on to what you want.


All times are GMT -4. The time now is 10:18.