CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Meshing & Mesh Conversion

Mesh not 1D or 2D

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Display Modes
Old   October 10, 2012, 12:34
Default Mesh not 1D or 2D
  #1
Member
 
Sagun Tripathi
Join Date: Aug 2012
Location: Amherst, USA
Posts: 78
Rep Power: 5
Sagun is on a distinguished road
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
Sagun is offline   Reply With Quote

Old   October 13, 2012, 05:18
Default
  #2
Super Moderator
 
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 8,507
Blog Entries: 34
Rep Power: 86
wyldckat is just really nicewyldckat is just really nicewyldckat is just really nicewyldckat is just really nice
Greetings Sagun,

You can use ParaView to help you diagnose what's wrong:
Code:
paraFoam -block
Best regards,
Bruno
wyldckat is offline   Reply With Quote

Old   October 15, 2012, 07:34
Default
  #3
Member
 
Sagun Tripathi
Join Date: Aug 2012
Location: Amherst, USA
Posts: 78
Rep Power: 5
Sagun is on a distinguished road
Thank you.
Sagun is offline   Reply With Quote

Reply

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
[ICEM] surface mesh merging problem everest ANSYS Meshing & Geometry 39 June 5, 2013 19:02
[ICEM] Problem making structural mesh on a surface froztbear ANSYS Meshing & Geometry 1 November 10, 2011 09:52
snappyHexMesh won't work - zeros everywhere! sc298 OpenFOAM Native Meshers: snappyHexMesh and Others 2 March 27, 2011 21:11
Meshing aifoil in ICEM student123a ANSYS Meshing & Geometry 13 December 8, 2010 11:40
2d irregular grid Remy Main CFD Forum 1 December 22, 2008 05:49


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