CFD Online Logo CFD Online URL
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Meshing & Mesh Conversion

[blockMesh] Meshing inside a cylinder so that the edges contain more cells than the center

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

LinkBack Thread Tools Search this Thread Display Modes
Old   April 11, 2010, 04:36
Default Meshing inside a cylinder so that the edges contain more cells than the center
New Member
Bertrand Delvaux
Join Date: Apr 2010
Posts: 29
Rep Power: 16
Bertrand is on a distinguished road

I am pretty new to OF, so I apologize if I may say trivial or stupid things...

Here is my problem : I gotta mesh a variable cross-sectional cylinder (axi-symmetrical geometry, but not axi-symmetrical solution).

I have tried several ways, but I would like some input on which one I should investigate further :

1. "Classical cylinder meshing", i.e. 5 blocks-based, with a cuboid in the centre and 4 blocks to fill the cylinder.
Problem : doesn't allow to have more cells on the edges of the cylinder, where limit layer imposes it, since the number of cells on an arc of the cylinder is dictated by the cuboid number of cells on one side. Except if one uses mergePatch, which makes it harder.

2. Small cylinder in the center, covered by 4 blocks that cover the rest of the shape.
Problem : using mergePatch doesn't seem to work on this kind of geometry, OF argues that one face would have less than 3 points, I guess this circle in the middle is not really OF-friendly

3. "Analytical function", I still didn't try this possibility, but I was wondering wether it could work to specify a function in the conventional cylindrical coordinates, like r(z), with the radius depending on the z-dimension only, with the number of cells depending somehow on the radius?

I hope it was "clear"...

Thank you for your time!


Bertrand is offline   Reply With Quote

Old   April 13, 2010, 04:46
Senior Member
romant's Avatar
Roman Thiele
Join Date: Aug 2009
Location: Eindhoven, NL
Posts: 374
Rep Power: 20
romant is on a distinguished road

have you tried looking into the possibility to use snappyHexMesh instead of blockMesh to create your mesh? Creating cylinders with this is quite easy and straightforward.

In linux there is salome and blender which can easily create these cylinders and output an stl file for use with snappyHexMesh.

otherwise, for your second approach I would suggest reading the tutorial on stress analysis, which could help.
romant is offline   Reply With Quote


Thread Tools Search this Thread
Search this Thread:

Advanced Search
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 Off
Pingbacks are On
Refbacks are On

Similar Threads
Thread Thread Starter Forum Replies Last Post
Foam::error::PrintStack almir OpenFOAM Running, Solving & CFD 91 December 21, 2022 04:50
[snappyHexMesh] Error snappyhexmesh - Multiple outside loops avinashjagdale OpenFOAM Meshing & Mesh Conversion 53 March 8, 2019 09:42
[snappyHexMesh] sHM layer process keeps getting killed MBttR OpenFOAM Meshing & Mesh Conversion 4 August 15, 2016 03:21
snappyHexMesh in parallel - FOAM Fatal IO Error mturcios777 OpenFOAM Running, Solving & CFD 4 August 10, 2012 19:18
physical boundary error!! kris Siemens 2 August 3, 2005 00:32

All times are GMT -4. The time now is 16:24.