CFD Online Discussion Forums

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   OpenFOAM Meshing & Mesh Conversion (https://www.cfd-online.com/Forums/openfoam-meshing/)
-   -   [blockMesh] blockMesh Problem (https://www.cfd-online.com/Forums/openfoam-meshing/128839-blockmesh-problem.html)

luigi21 January 22, 2014 09:53

blockMesh Problem
 
Hi all...

I'm trying to run a mesh using blockMesh but I get always the same error, if you guys know what the problem is, let me know it, I appreciate it so much

this is mesh file

vertices
(
( -8 0 0 )
( 6 0 0 )
( 17 0 0 )
( 32 0 0 )
( -8 0.3 0 )
( 6 0.3 0 )
( 12 0.3 0 )
( 14 0.3 0 )
( 17 0.3 0 )
( 32 0.3 0 )
( -8 0.5 0 )
( 6 0.5 0 )
( 12 0.5 0 )
( 14 0.5 0 )
( 17 0.5 0 )
( 32 0.5 0 )
( -8 0 1 )
( 6 0 1 )
( 17 0 1 )
( 32 0 1 )
( -8 0.3 1 )
( 6 0.3 1 )
( 12 0.3 1 )
( 14 0.3 1 )
( 17 0.3 1 )
( 32 0.3 1 )
( -8 0.5 1 )
( 6 0.5 1 )
( 12 0.5 1 )
( 14 0.5 1 )
( 17 0.5 1 )
( 32 0.5 1 )
);

blocks
(
hex (0 1 5 4 16 17 21 20) ( 700 38 1 ) simpleGrading (1 1 1)
hex (1 1 6 5 17 17 22 21) ( 300 38 1 ) simpleGrading (1 1 1)
hex (2 2 8 7 18 18 24 23) ( 150 38 1 ) simpleGrading (1 1 1)
hex (2 3 9 8 18 19 25 24) ( 750 38 1 ) simpleGrading (1 1 1)
hex (4 5 11 10 20 21 27 26) ( 700 25 1 ) simpleGrading (1 1 1)
hex (5 6 12 11 21 22 28 27) ( 300 25 1 ) simpleGrading (1 1 1)
hex (6 7 13 12 22 23 29 28) ( 100 25 1 ) simpleGrading (1 1 1)
hex (7 8 14 13 23 24 30 29) ( 150 25 1 ) simpleGrading (1 1 1)
hex (8 9 15 14 24 25 31 30) ( 750 25 1 ) simpleGrading (1 1 1)
);

edges
(
);

patches
(
patch inlet
(
(0 16 20 4)
(4 20 26 10)
)
wall bottom
(
(0 1 17 16)
(1 6 22 17)
(6 7 23 22)
(7 2 18 23)
(2 3 19 18)
)
patch outlet
(
(9 25 19 3)
(15 31 25 9)
)

patch atmosphere
(
(10 26 27 11)
(11 27 28 12)
(12 28 29 13)
(13 29 30 14)
(14 30 31 15)
)

empty frontBack
(
(0 4 5 1)
(1 5 6 6)
(7 8 2 2)
(2 8 9 3)
(4 10 11 5)
(5 11 12 6)
(6 12 13 7)
(7 13 14 8)
(8 14 15 9)
(16 17 21 20)
(22 21 17 17)
(18 24 23 23)
(18 19 25 24)
(20 21 27 26)
(21 22 28 27)
(22 23 29 28)
(23 24 30 29)
(24 25 31 30)
)
);

mergePatchPairs
(
);

and this is the error and the checkMesh what I always get. I've checked the mesh but I didn't catch the error.

Checking topology...
Boundary definition OK.
***Total number of faces on empty patches is not divisible by the number of cells in the mesh. Hence this mesh is not 1D or 2D.
Cell to face addressing OK.
Point usage OK.
Upper triangular ordering OK.
Face vertices OK.
Number of regions: 1 (OK).

Checking patch topology for multiply connected surfaces ...
Patch Faces Points Surface topology
inlet 63 128 ok (non-closed singly connected)
bottom 1626 3256 ok (non-closed singly connected)
outlet 63 128 ok (non-closed singly connected)
atmosphere 2000 4002 ok (non-closed singly connected)
frontBack 244400 247744 ok (non-closed singly connected)
defaultFaces 24 42 ok (non-closed singly connected)

Checking geometry...
Overall domain bounding box (-8 0 0) (32 0.5 1)
Mesh (non-empty, non-wedge) directions (1 1 0)
Mesh (non-empty) directions (1 1 0)
All edges aligned with or perpendicular to non-empty directions.
Boundary openness (-5.45052e-18 1.70473e-14 2.95858e-13) OK.
Max cell openness = 1.48256e-16 OK.
Max aspect ratio = 22.4667 OK.
***Zero or negative face area detected. Minimum area: 0
<<Writing 24 zero area faces to set zeroAreaFaces
Min volume = 2.07754e-06. Max volume = 0.00745783. Total volume = 18.1693. Cell volumes OK.
Mesh non-orthogonality Max: 100.102 average: 27.0394
*Number of severely non-orthogonal faces: 27801.
***Number of non-orthogonality errors: 7.
<<Writing 27808 non-orthogonal faces to set nonOrthoFaces
***Error in face pyramids: 16 faces are incorrectly oriented.
<<Writing 16 faces with incorrect orientation to set wrongOrientedFaces
***Max skewness = 12.977, 1 highly skew faces detected which may impair the quality of the results
<<Writing 1 skew faces to set skewFaces
Coupled point location match (average 0) OK.

Failed 4 mesh checks.


this is the blockMesh running

Create time

Creating block mesh from
"/home/miguel/OpenFoam/waves2Foam/tutorials/waveFoam/Sensibilidades/Caso_cluster/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 888
Found 2 undefined faces in mesh; adding to default patch.

Check topology

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

Creating block offsets
Creating merge list ..

Creating polyMesh from blockMesh
Creating patches
Creating cells
Creating points with scale 1

Writing polyMesh
----------------
Mesh Information
----------------
boundingBox: (-8 0 0) (32 0.5 1)
nPoints: 247744
nCells: 122200
nFaces: 490475
nInternalFaces: 242299
----------------
Patches
----------------
patch 0 (start: 242299 size: 63) name: inlet
patch 1 (start: 242362 size: 1626) name: bottom
patch 2 (start: 243988 size: 63) name: outlet
patch 3 (start: 244051 size: 2000) name: atmosphere
patch 4 (start: 246051 size: 244400) name: frontBack
patch 5 (start: 490451 size: 24) name: defaultFaces

End

alexeym January 23, 2014 07:53

Hi,

Maybe there's an error in the blockMesh you've posted as I got following error while creating a mesh with the text you've posted:

Code:

--> FOAM FATAL ERROR:
face 1 in patch 4 does not have neighbour cell face: 4(1 5 6 6)



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