[upuli]

Dear all



I want to mesh the following shape given in the attachment. Following is the block Mesh entry used. When the commented lines are included a check Mesh error occurs. I will be very much thankful if someone can help me to overcome the problem.


convertToMeters 1.0;

vertices
(
(-0.06 -0.17 0.06)
(0.06 -0.17 0.06)
(0.06 -0.17 -0.06)
(-0.06 -0.17 -0.06)
(-0.06 0.0 0.06)
(0.06 0.0 0.06)
(0.06 0.0 -0.06)
(-0.06 0.0 -0.06)
(-0.06 0.1 -0.06)
(-0.06 0.1 0.06)
(0.06 0.1 0.06)
(0.06 0.1 -0.06)
(0.06 0.27 -0.06)
(-0.06 0.27 -0.06)
(-0.06 0.27 0.06)
(0.06 0.27 0.06)
(0.06 0.27 0.06)
/*(-0.03 -0.17 0.03)//16
(0.03 -0.17 0.03)//17
(0.03 -0.17 -0.03)
(-0.03 -0.17 -0.03)
(-0.03 0 0.03)//20
(0.03 0 0.03)//21
(0.03 0 -0.03)
(-0.03 0 -0.03)*/
);

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

patches
(
wall Walls
(
(0 3 7 4)
(4 7 8 9)
(9 8 13 14)
(1 5 6 2)
(5 10 11 6)
(10 15 12 11)
//(19 23 20 16)
//(17 21 22 18)
//(16 20 21 17)
//(18 22 23 19)
//(0 1 2 3)
)
patch Inlet
(
// (16 17 18 19)
(0 1 2 3)
)
patch Outlet
(
(12 15 14 13)
)
empty back
(
(2 6 7 3)
(6 11 8 7)
( 11 12 13 8)
)
empty front
(
(4 9 10 5)
(0 4 5 1)
(9 14 15 10)
)
);

edges
(
);

mergePatchPairs
(
);

error log




Create time

Create polyMesh for time = 0

Time = 0

Mesh stats
points: 1422
internal points: 0
faces: 2599
internal faces: 1181
cells: 630
faces per cell: 6
boundary patches: 6
point zones: 0
face zones: 0
cell zones: 0

Overall number of cells of each type:
hexahedra: 630
prisms: 0
wedges: 0
pyramids: 0
tet wedges: 0
tetrahedra: 0
polyhedra: 0

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: 2
The mesh has multiple regions which are not connected by any face.
<<Writing region information to "0/cellToRegion"
<<Writing region 0 with 528 cells to cellSet region0
<<Writing region 1 with 102 cells to cellSet region1

Checking patch topology for multiply connected surfaces...
Patch Faces Points Surface topology
Walls 338 454 ok (non-closed singly connected)
Inlet 6 14 ok (non-closed singly connected)
Outlet 12 26 ok (non-closed singly connected)
back 528 585 ok (non-closed singly connected)
front 528 585 ok (non-closed singly connected)
defaultFaces 6 14 ok (non-closed singly connected)

Checking geometry...
Overall domain bounding box (-0.06 -0.17 -0.06) (0.06 0.27 0.06)
Mesh has 0 geometric (non-empty/wedge) directions (0 0 0)
Mesh has 0 solution (non-empty) directions (0 0 0)
***Number of edges not aligned with or perpendicular to non-empty directions: 318
<<Writing 234 points on non-aligned edges to set nonAlignedEdges
Boundary openness (-1.76899e-16 1.11873e-17 2.44721e-18) OK.
Max cell openness = 7.39533e-16 OK.
Max aspect ratio = -1 OK.
Minimum face area = 1.56671e-05. Maximum face area = 0.00464364. Face area magnitudes OK.
***Zero or negative cell volume detected. Minimum negative volume: -9.70098e-06, Number of negative volume cells: 54
<<Writing 54 zero volume cells to set zeroVolumeCells
Mesh non-orthogonality Max: 164.519 average: 33.1723
*Number of severely non-orthogonal (> 70 degrees) faces: 31.
***Number of non-orthogonality errors: 105.
<<Writing 136 non-orthogonal faces to set nonOrthoFaces
***Error in face pyramids: 306 faces are incorrectly oriented.
<<Writing 215 faces with incorrect orientation to set wrongOrientedFaces
***Max skewness = 40.6103, 85 highly skew faces detected which may impair the quality of the results
<<Writing 85 skew faces to set skewFaces
Coupled point location match (average 0) OK.

Failed 5 mesh checks.

End

[senthilathiban]

/*(-0.03 -0.17 0.03)//17
(0.03 -0.17 0.03)//18
(0.03 -0.17 -0.03)
(-0.03 -0.17 -0.03)
(-0.03 0 0.03)//21
(0.03 0 0.03)//22
(0.03 0 -0.03)
(-0.03 0 -0.03)*/