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

Meshes made of several blocks

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

Reply
 
LinkBack Thread Tools Display Modes
Old   April 9, 2015, 04:51
Default Meshes made of several blocks
  #1
New Member
 
Philippe Gambron
Join Date: Apr 2015
Location: Oxford
Posts: 5
Rep Power: 4
zaphod is on a distinguished road
Hi!

I am trying to make a mesh made of several blocks in order to have a geometry containing obstacles but I keep obtaining error messages about faces that don't have neighbour cell faces.

I can reproduce the problem by trying to make a mesh made of 2 blocks inside a cube.

In blockMeshDict, I create a cube and I add 4 points in the middle, along the y axis, to slice it in half.
Quote:
vertices
(
( 0 0 0)
(10 0 0)
(10 10 0)
( 0 10 0)
( 0 0 10)
(10 0 10)
(10 10 10)
( 0 10 10)

( 0 5 0)
( 0 10 0)
(10 5 0)
(10 5 10)
);
Then I attempt to make 2 blocks.
Quote:
blocks
(
hex ( 0 8 9 4 1 10 11 5) (10 5 10) simpleGrading (1 1 1)
hex ( 3 7 9 8 2 6 11 10) (10 5 10) simpleGrading (1 1 1)
);
But I keep obtaining the following error message.
Quote:
--> FOAM FATAL ERROR:
face 0 in patch 0 does not have neighbour cell face: 4(0 3 2 1)
What am I doing wrong?

Thank you!

Cheers,

Philippe
zaphod is offline   Reply With Quote

Old   April 9, 2015, 05:05
Default
  #2
Senior Member
 
Alexey Matveichev
Join Date: Aug 2011
Location: Nancy, France
Posts: 1,760
Rep Power: 29
alexeym will become famous soon enoughalexeym will become famous soon enough
Send a message via Skype™ to alexeym
Hi,

In general the error means you have messed up definition of patches (boundary section). Can you post your whole blockMeshDict as an attachment?
alexeym is offline   Reply With Quote

Old   April 9, 2015, 05:09
Default
  #3
New Member
 
Philippe Gambron
Join Date: Apr 2015
Location: Oxford
Posts: 5
Rep Power: 4
zaphod is on a distinguished road
Hello!

Thank you for your reply.

Here is the boundary section.

Quote:
boundary
(
walls
{
type wall;
faces
(
(0 3 2 1)
(0 1 5 4)
(1 2 6 5)
(2 3 7 6)
(0 4 7 3)
(4 5 6 7)
);
}
);
Should the boundaries be cut in the same way as the cube? I had thought it was independent. I am going to try.

Cheers,

Philippe
zaphod is offline   Reply With Quote

Old   April 9, 2015, 05:18
Default
  #4
Senior Member
 
Alexey Matveichev
Join Date: Aug 2011
Location: Nancy, France
Posts: 1,760
Rep Power: 29
alexeym will become famous soon enoughalexeym will become famous soon enough
Send a message via Skype™ to alexeym
Hi,

Patches in general composed of boundaries of the blocks. If you cut cube, you cut boundaries. So what was described before cutting as a single face should be described as several faces.
alexeym is offline   Reply With Quote

Old   April 9, 2015, 05:42
Default
  #5
New Member
 
Philippe Gambron
Join Date: Apr 2015
Location: Oxford
Posts: 5
Rep Power: 4
zaphod is on a distinguished road
I have sliced the boundaries in 2.

Quote:
boundary
(
walls
{
type wall;
faces
(
( 0 4 9 8)
( 1 10 11 5)
( 3 8 9 7)
( 2 6 11 10)
( 0 8 10 1)
( 4 5 11 9)
( 2 10 8 3)
( 6 7 9 11)
( 0 1 5 4)
( 3 7 6 2)
);
}
);
But, now, I obtain the following error message:
Quote:
--> FOAM FATAL ERROR:
Inconsistent number of faces between block pair 0 and 1
However, t seems to me that there is the same number of cells on both sides of the slice.

Thanks.

Philippe
zaphod is offline   Reply With Quote

Old   April 9, 2015, 06:58
Default
  #6
Senior Member
 
Alexey Matveichev
Join Date: Aug 2011
Location: Nancy, France
Posts: 1,760
Rep Power: 29
alexeym will become famous soon enoughalexeym will become famous soon enough
Send a message via Skype™ to alexeym
Hi,

Start with visualizing your blocks, to check if you have made mistake there. I.e. you comment out patches definitions, run blockMesh, open mesh in paraview. Your vertices look rather suspicious, among four additional points for some reason ( 0 10 0) appeared. Also I am not quite sure in correctness of your vertices numbering in blocks definitions.

I can try to reconstruct the case from the pieces you have posted and try to look for error but I would rather wait till you post archive with the case files.
alexeym is offline   Reply With Quote

Reply

Tags
block, blockmesh

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] Advice on genetrating 2D blocks required user_1 ANSYS Meshing & Geometry 0 December 3, 2014 06:46
[ICEM] Two boxes connected by a cylinder, blocks merge where I don't want them to merge. Polarbear ANSYS Meshing & Geometry 4 April 30, 2014 15:21
Hex and Tet meshes - simplefoam comparison danvica OpenFOAM Running, Solving & CFD 10 January 4, 2013 02:18
Getting prism to inflate into mixed tet-hex meshes Joe CFX 16 October 10, 2011 07:06
Large 3D tetrahedral meshes Aldo Bonfiglioli Main CFD Forum 4 August 27, 1999 03:33


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