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/)
-   -   FOAM FATAL ERROR: Inconsistent number of faces blockMesh::createMergeList() line 193 (http://www.cfd-online.com/Forums/openfoam-meshing-blockmesh/80098-foam-fatal-error-inconsistent-number-faces-blockmesh-createmergelist-line-193-a.html)

Hengel September 15, 2010 08:30

FOAM FATAL ERROR: Inconsistent number of faces blockMesh::createMergeList() line 193
 
Hallo!
I have got some troubles with blockMesh:


Here my problem
Code:franz@franz-desktop:~/OpenFOAM/franz-1.7.1/DA_aktuell/2D_FreeConvection/020_FreeConvection_SKE_buoyantSimpleFoam$ blockMesh
/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 1.7.1 |
| \\ / A nd | Web: www.OpenFOAM.com |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Build : 1.7.1-03e7e056c215
Exec : blockMesh
Date : Sep 15 2010
Time : 13:12:38
Host : franz-desktop
PID : 4753
Case : /home/franz/OpenFOAM/franz-1.7.1/DA_aktuell/2D_FreeConvection/020_FreeConvection_SKE_buoyantSimpleFoam
nProcs : 1
SigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time


Creating block mesh from
"/home/franz/OpenFOAM/franz-1.7.1/DA_aktuell/2D_FreeConvection/020_FreeConvection_SKE_buoyantSimpleFoam/constant/polyMesh/blockMeshDict"


Creating blockCorners

Creating curved edges
Creating blocks
Creating patches
Creating block mesh topology
Default patch type set to empty
Check block mesh topology
Basic statistics
Number of internal faces : 12
Number of boundary faces : 30
Number of defined boundary faces : 30
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 3

From function blockMesh::createMergeList()
in file createMergeList.C at line 193.

FOAM exiting



You can see the blockMeshDict-File below:
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 1.6 |
| \\ / A nd | Web: http://www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
object blockMeshDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
convertToMeters 1;
vertices
(
(0 0 0) //Punkt 0
(0.00762 0 0) //Punkt 1
(0.06858 0 0) //Punkt 2
(0.0762 0 0) //Punkt 3

(0 0.00762 0) //Punkt 4
(0.00762 0.00762 0) //Punkt 5
(0.06858 0.00762 0) //Punkt 6
(0.0762 0.00762 0) //Punkt 7

(0 2.17238 0) //Punkt 8
(0.00762 2.17238 0) //Punkt 9
(0.06858 2.17238 0) //Punkt 10
(0.0762 2.17238 0) //Punkt 11

(0 2.18 0) //Punkt 12
(0.00762 2.18 0) //Punkt 13
(0.06858 2.18 0) //Punkt 14
(0.0762 2.18 0) //Punkt 15

(0 0 0.01) //Punkt 16
(0.00762 0 0.01) //Punkt 17
(0.06858 0 0.01) //Punkt 18
(0.0762 0 0.01) //Punkt 19

(0 0.00762 0.01) //Punkt 20
(0.00762 0.00762 0.01) //Punkt 21
(0.06858 0.00762 0.01) //Punkt 22
(0.0762 0.00762 0.01) //Punkt 23

(0 2.17238 0.01) //Punkt 24
(0.00762 2.17238 0.01) //Punkt 25
(0.06858 2.17238 0.01) //Punkt 26
(0.0762 2.17238 0.01) //Punkt 27

(0 2.18 0.01) //Punkt 28
(0.00762 2.18 0.01) //Punkt 29
(0.06858 2.18 0.01) //Punkt 30
(0.0762 2.18 0.01) //Punkt 31

);
blocks
(
hex (0 1 5 4 16 17 21 20) (10 1 1) simpleGrading (1 1 1) // Block 0
hex (1 2 6 5 17 18 22 21) (1 1 1) simpleGrading (1 1 1) // Block I
hex (2 3 7 6 18 19 23 22) (1 1 1) simpleGrading (1 1 1) // Block II
hex (4 5 9 8 20 21 25 24) (1 1 1) simpleGrading (1 1 1) // Block III
hex (5 6 10 9 21 22 26 25) (1 1 1) simpleGrading (1 1 1) // Block IV
hex (6 7 11 10 22 23 27 26) (1 1 1) simpleGrading (1 1 1) // Block V
hex (8 9 13 12 24 25 29 28) (1 1 1) simpleGrading (1 1 1) // Block VI
hex (9 10 14 13 25 26 30 29) (1 1 1) simpleGrading (1 1 1) // Block VII
hex (10 11 15 14 26 27 31 30) (1 1 1) simpleGrading (1 1 1) // Block VIII

);
edges
(
);

patches
(
wall hot
(
(19 23 7 3)
(23 27 11 7)
(27 31 15 11)
)
wall cold
(
(0 16 20 4)
(4 20 24 8)
(8 24 28 12)
)
wall topAndBottom
(
(0 1 17 16)
(1 2 18 17)
(2 3 19 18)

(12 28 29 13)
(13 29 30 14)
(14 30 31 15)
)
empty frontAndBackPlanes
(
(0 4 5 1)
(1 5 6 2)
(2 6 7 3)

(4 8 9 5)
(5 9 10 6)
(6 10 11 7)

(8 12 13 9)
(9 13 14 10)
(10 14 15 11)


(16 20 21 17)
(17 21 22 18)
(18 22 23 19)

(20 24 25 21)
(21 25 26 22)
(22 26 27 23)

(24 28 29 25)
(25 29 30 26)
(26 30 31 27)
)
);

mergePatchPairs
(
);

// ************************************************** *********************** //
If I replace the values for the meshrefinement from (10 1 1) to (1 1 1), see below, then blockMesh will work
Code: hex (0 1 5 4 16 17 21 20) (1 1 1) simpleGrading (1 1 1) // Block 0


Can anyone help me to solve this problem?


I guess there can be a problem with the block IV coz it is positioned in the middle of the geometry and around there are all other blocks. The block IV does not have any contacts to the walls



Thanks for your support


best regards


Franz


marval September 17, 2010 04:27

Quote:

Originally Posted by Hengel (Post 275290)
If I replace the values for the meshrefinement from (10 1 1) to (1 1 1), see below, then blockMesh will work
Code: hex (0 1 5 4 16 17 21 20) (1 1 1) simpleGrading (1 1 1) // Block

Hello Franz!

I don't see where the problem is, you seem to have figured it out by yourself! :)
The change you did had nothing to do with the refinement (grading), it was the number of elements in each direction (x y z).

Regards
Marco

Hengel September 17, 2010 05:07

Hallo marval!

Thank you for your attention!

I solved the problem.

I think the problem is the amout of cells just in the first block. The neighbor cells requires the same amount of cells. If you do not do this then you will get an error like this:

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

From function blockMesh::createMergeList()
in file createMergeList.C at line 193.

FOAM exiting


Solution:

I set the some number of cells to the neighbor cells and then it works fine.

Btw., can blockMesh work with hanging nodes?

Thx

best regards

Franz

marval September 17, 2010 05:36

Quote:

Originally Posted by Hengel (Post 275567)

I think the problem is the amout of cells just in the first block. The neighbor cells requires the same amount of cells. If you do not do this then you will get an error like this:

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

From function blockMesh::createMergeList()
in file createMergeList.C at line 193.

FOAM exiting


Solution:

I set the some number of cells to the neighbor cells and then it works fine.


Exactly! You only have one cell in each direction for each block?

Quote:

Originally Posted by Hengel (Post 275567)
Btw., can blockMesh work with hanging nodes?


I'm not sure what that means, is it a defined point that's not used to construct the geometry?

Regards
Marco

Hengel September 17, 2010 06:01

Quote:

Quote:
Originally Posted by Hengel http://www.cfd-online.com/Forums/ima...s/viewpost.gif

I think the problem is the amout of cells just in the first block. The neighbor cells requires the same amount of cells. If you do not do this then you will get an error like this:

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

From function blockMesh::createMergeList()
in file createMergeList.C at line 193.

FOAM exiting


Solution:

I set the some number of cells to the neighbor cells and then it works fine.



Exactly! You only have one cell in each direction for each block?

Yes secondly.

First of all I had 10 cells in the x direction in the first block, so I had to change the upper and the lower neighbor. Both have 10 cells in the x-direction now.


Quote:

Quote:
Originally Posted by Hengel http://www.cfd-online.com/Forums/ima...s/viewpost.gif
Btw., can blockMesh work with hanging nodes?



I'm not sure what that means, is it a defined point that's not used to construct the geometry?
Hanging-node splits e.g. cells into two pieces but do not change the primary size of the cell. You can compare it with fluent at Yplus/Ystar Adaption. Here you can adapt the y+ value with hanging nodes.


Now I think, If the hanging-node-methode is working in OpenFoam the main problem could be the aspect ratio. The max. of aspect ratio is 10:1.

best

Franz

Ray092 December 30, 2015 10:42

help me with this error please
/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.3.x |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
/* Windows port by CFD support (www.cfdsupport.com) [based on Symscape] *\
\*---------------------------------------------------------------------------*/
Build : 2.3.x-819030ed51bd
Exec : D:\openfoamwindows\OpenFOAM\cygwin64\opt\OpenFOAM\ OpenFOAM-2.3.x\platfo rms\cygwin64mingw-w64DPOpt\bin\blockMesh.exe
Date : Dec 30 2015
Time : 10:26:23
Host : "SATYAKI-PC"
PID : 4004
Case : D:/openfoamwindows/OpenFOAM/cygwin64/home/satyaki/FOAM_RUN/tutorials/mu ltiphase/multiphaseEulerFoam/bubbleColumnMod1
nProcs : 1
fileModificationChecking : Monitoring run-time modified files using timeStampMas ter
allowSystemOperations : Disallowing user-supplied system call operations

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

Creating block mesh from
"D:/openfoamwindows/OpenFOAM/cygwin64/home/satyaki/FOAM_RUN/tutorials/multip hase/multiphaseEulerFoam/bubbleColumnMod1/constant/polyMesh/blockMeshDict"
Creating curved edges
Creating topology blocks
Creating topology patches

Reading patches section

Creating block mesh topology

Reading physicalType from existing boundary file

Default patch type set to empty
--> FOAM Warning :
From function polyMesh:: polyMesh(... construct from shapes...)
in file meshes/polyMesh/polyMeshFromShapeMesh.C at line 903
Found 14 undefined faces in mesh; adding to default patch.

Check topology

Basic statistics
Number of internal faces : 4
Number of boundary faces : 22
Number of defined boundary faces : 22
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



here's my blockmeshdict file

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

convertToMeters 0.001;

vertices
(
(0 0 0)
(90 0 0)
(90 200 0)
(0 200 0)
(99 0 0)
(99 200 0)
(101 200 0)
(101 0 0)
(110 0 0)
(110 200 0)
(200 200 0)
(200 0 0)
(0 0 100)
(90 0 100)
(90 200 100)
(0 200 100)
(99 0 100)
(99 200 100)
(101 200 100)
(101 0 100)
(110 0 100)
(110 200 100)
(200 200 100)
(200 0 100)
);

blocks
(
hex (0 1 2 3 12 13 14 15) (10 200 1) simpleGrading (1 1 1)
hex (2 1 4 5 14 13 16 17) (10 200 1) simplegrading (1 1 1)
hex (5 4 7 6 17 16 19 18) (10 200 1) simplegrading (1 1 1)
hex (6 7 8 9 18 19 20 21) (10 200 1) simplegrading (1 1 1)
hex (9 8 11 10 21 20 23 22) (10 200 1) simplegrading (1 1 1)
);

edges
(
);

patches
(
patch inlet
(
(4 16 19 7)
)
patch outlet
(
(3 15 14 2)
(2 14 17 5)
(5 17 18 6)
(6 18 21 9)
(9 21 22 10)
)
wall walls
(
(0 12 15 3)
(11 23 22 10)
)
);

mergePatchPairs
(
);

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

Persia December 30, 2015 22:21

The first block should be defined like this to be consistent with the other blocks:

hex (3 0 1 2 15 12 13 14) (10 200 1) simpleGrading (1 1 1)


All times are GMT -4. The time now is 21:01.