CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM Meshing & Mesh Conversion (http://www.cfd-online.com/Forums/openfoam-meshing/)
-   -   Mesh not 1D or 2D (http://www.cfd-online.com/Forums/openfoam-meshing/107967-mesh-not-1d-2d.html)

Sagun October 10, 2012 12:34

Mesh not 1D or 2D
 
Hello,

I found this online tutorial of transient filling of water in a bottle initially filled with air. I copied the blockMeshDict file mentioned as it is in the tutorial, but whenever I do checkMesh I get the following output:

Mesh stats
points: 18712
internal points: 0
faces: 36805
internal faces: 18095
cells: 9150
boundary patches: 5
point zones: 0
face zones: 0
cell zones: 0

Overall number of cells of each type:
hexahedra: 9150
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: 1 (OK).

Checking patch topology for multiply connected surfaces ...
Patch Faces Points Surface topology
inlet 10 22 ok (non-closed singly connected)
inletWall 30 64 ok (non-closed singly connected)
bottleWall 290 582 ok (non-closed singly connected)
atmosphere 18300 18712 ok (non-closed singly connected)
defaultFaces 80 164 ok (non-closed singly connected)

Checking geometry...
Overall domain bounding box (-0.05 0 0) (0.05 0.23 0.001)
Mesh (non-empty, non-wedge) directions (1 0 1)
Mesh (non-empty) directions (1 0 1)
***Number of edges not aligned with or perpendicular to non-empty directions: 100
<<Writing 120 points on non-aligned edges to set nonAlignedEdges
Boundary openness (-2.52649e-18 2.85249e-19 5.08717e-15) OK.
Max cell openness = 1.32349e-16 OK.
Max aspect ratio = 2 OK.
Minumum face area = 1e-06. Maximum face area = 4e-06. Face area magnitudes OK.
Min volume = 2e-09. Max volume = 4e-09. Total volume = 2.045e-05. Cell volumes OK.
Mesh non-orthogonality Max: 24.2277 average: 1.05611
Non-orthogonality check OK.
Face pyramids OK.
Max skewness = 0.456717 OK.
Coupled point location match (average 0) OK.

Failed 1 mesh checks.

Also, I named the empty patch as frontAndBack (as in the tutorial) but after running blockMesh I see that the empty patch gets renamed to defaultFaces. I don't seem to be able to understand why this is happening.

Here is the blockMeshDict file:

convertToMeters 0.001;

vertices
(
(-10 230 0)
(10 230 0)
(10 220 0)
(-10 220 0)
(5 210 0)
(-5 210 0)
(5 200 0)
(-5 200 0)
(5 0 0)
(-5 0 0)
(50 200 0)
(50 0 0)
(-50 200 0)
(-50 0 0)

(-10 230 1)
(10 230 1)
(10 220 1)
(-10 220 1)
(5 210 1)
(-5 210 1)
(5 200 1)
(-5 200 1)
(5 0 1)
(-5 0 1)
(50 200 1)
(50 0 1)
(-50 200 1)
(-50 0 1)
);

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

edges
(
);

boundary
(
inlet
{
type patch;
faces
(
(0 1 15 14)
);
}

inletWall
{
type wall;
faces
(
(0 3 17 14)
(3 5 19 17)
(5 7 21 19)
(1 2 16 15)
(2 4 18 16)
(4 6 20 18)
);
}

bottleWall
{
type wall;
faces
(
(10 11 25 24)
(11 8 22 25)
(8 9 23 22)
(9 13 27 23)
(13 12 26 27)
);
}

atmosphere
{
type patch;
faces
(
(6 10 24 20)
(12 7 21 26)
):
}

frontAndBack
{
type empty;
faces
(
(0 1 2 3)
(3 2 4 5)
(4 6 7 5)
(6 10 11 8)
(6 8 9 7)
(7 9 13 12)

(14 15 16 17)
(17 16 18 19)
(19 18 20 21)
(20 24 25 22)
(20 22 23 21)
(21 23 27 26)
);
}
);

mergePatchPairs
(
);

Please help.

Regards,
Sagun

wyldckat October 13, 2012 05:18

Greetings Sagun,

You can use ParaView to help you diagnose what's wrong:
Code:

paraFoam -block
Best regards,
Bruno

Sagun October 15, 2012 07:34

Thank you. :)


All times are GMT -4. The time now is 09:28.