CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Native Meshers: blockMesh

Is it possible to have a mesh to represent a domain with a hole

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

Reply
 
LinkBack Thread Tools Display Modes
Old   May 3, 2006, 05:09
Default Hi, (I had originaly posted
  #1
New Member
 
kmurphy
Join Date: Mar 2009
Posts: 10
Rep Power: 8
kmurphy is on a distinguished road
Hi,

(I had originaly posted this message in Meshing / Mesh conversion: gmsh sub group in error, sorry)

I am trying to construct a mesh to represent a rectangular region with a single circular hole. I started with the mesh used in tutorials/potentalFoam/cylinder (which matches half of my required mesh) and tried to extend this. When I do not split the blocks, i.e., using

hex (# # # # # # # #) (1 1 1) simpleGrading(1 1 1)

for each block, everything seems OK. Program checkBlock does not generate any warnings and the resulting mesh looks like .

The output from checkMesh is in checkMesh_2.log and the blockMeshDict file is blockMeshDict_1 .

Adding grading to most of the blocks works as I would have expected, to produce the following mesh .

The output from checkMesh is in checkMesh_2.log and
the blockMeshDict file is blockMeshDict_2 and everything appears OK.

However, when I try to add grading to any of the other blocks I do not obtain the expected (to me) mesh. In the following diagram and files I have only
added grading to two of the problem blocks to keep the mesh as simple as possible and the problem seem to be independent of the grading used in the other blocks.

.

The output from checkMesh is checkMesh_3.log and the blockMeshDict file is blockMeshDict_3 .

In the log of checkMesh there were reports of

High aspect ratio for cell 10: 4.44057e+199
High aspect ratio for cell 11: 4.44058e+199

and

Severe non-orthogonality detected for face ...

which doesn't appear to be too healthy.

I have tried to switch the order of the vertices without success. Any advice would be great,

Thanks, Kieran
kmurphy is offline   Reply With Quote

Old   May 3, 2006, 05:24
Default Looks like you've messed up th
  #2
Senior Member
 
Hrvoje Jasak
Join Date: Mar 2009
Location: London, England
Posts: 1,758
Rep Power: 21
hjasak will become famous soon enough
Looks like you've messed up the mesh definition - this kind of mesh should be pretty easy to do. As an encouragement, have one of mine :-) This is a mesh for vertex shedding behind a cylinder.

From the pictures it seems to me that the grading is giving you trouble: remember, the grading coefficients are defined in the coordinate system of the block and not the global one.

Good luck,

Hrv

blockMeshDict
__________________
Hrvoje Jasak
Providing commercial FOAM/OpenFOAM and CFD Consulting: http://wikki.co.uk
hjasak is offline   Reply With Quote

Old   May 3, 2006, 17:10
Default Hi Hrvoje, Thanks for your bl
  #3
New Member
 
kmurphy
Join Date: Mar 2009
Posts: 10
Rep Power: 8
kmurphy is on a distinguished road
Hi Hrvoje,
Thanks for your blockMeshDict. I should able to use it to reconstruct mine (without the error this time). Regards, Kieran.
kmurphy is offline   Reply With Quote

Old   October 13, 2014, 17:23
Default
  #4
New Member
 
Join Date: Aug 2012
Posts: 24
Rep Power: 4
cesarjets is on a distinguished road
Quote:
Originally Posted by hjasak View Post
Looks like you've messed up the mesh definition - this kind of mesh should be pretty easy to do. As an encouragement, have one of mine :-) This is a mesh for vertex shedding behind a cylinder.

From the pictures it seems to me that the grading is giving you trouble: remember, the grading coefficients are defined in the coordinate system of the block and not the global one.

Good luck,

Hrv

blockMeshDict

could you share the whole file of this case? I mean I would like to check how you set up the U/P files?
thanks
cesarjets 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
How to represent the divergence of a tensor times a velocity titio OpenFOAM Running, Solving & CFD 0 January 15, 2009 17:51
How to find hole in mesh prateek jain CFX 5 May 24, 2007 07:06
How to represent two phase flow in Gambit K.Baker FLUENT 4 May 22, 2006 08:03
Is it possible to have a mesh to represent a domain with a hole kmurphy Open Source Meshers: Gmsh, Netgen, CGNS, ... 1 May 3, 2006 04:45
What does the variables R and nuTilda represent in the icofoam ztdep OpenFOAM Running, Solving & CFD 2 August 18, 2005 22:21


All times are GMT -4. The time now is 21:11.