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

[snappyHexMesh] Using pentagonal cells for refinement

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

Like Tree2Likes
  • 1 Post By Yann
  • 1 Post By katiuskas

LinkBack Thread Tools Search this Thread Display Modes
Old   April 14, 2022, 05:51
Default Using pentagonal cells for refinement
New Member
Join Date: Mar 2022
Location: Mallorca
Posts: 16
Rep Power: 3
katiuskas is on a distinguished road
Dear all,

I have been trying to refine my (2D) mesh in a certain region, let's say with a refinement level of 1. I have done so in two different ways: using blockMesh and then mergePatchPairs, and through snappyHexMesh. Both methods yield similar results: the "transition zones" between the different regions of the mesh use tetra elements (see first figure). In spannyHexMesh, I know you can control the "smoothness" of this transition with the nCellsBetweenLevels entry. So far so good...

My doubt comes from having seen meshes where no transition zone exists at all by using snappyHexMesh (see the second figure, taken from the wiki). Note that this is neither a non-conformal mesh nor a "hanging node" situation. It's just that some of the cells are now pentagons instead of squares (if you allow me to use the 2D jargon) and thus can have 5 neighbours instead of 4. Now, my question can be summarized into the following:
  1. How can I create such a mesh? Every time I try, I always get tetra elements.
  2. Would this "penta" mesh behave better than mine in terms of mesh quality and solver?

Thank you very much!
Attached Images
File Type: png my_mesh.PNG (122.2 KB, 17 views)
File Type: jpg penta_mesh.jpg (83.7 KB, 19 views)
katiuskas is offline   Reply With Quote

Old   April 14, 2022, 10:36
Senior Member
Join Date: Apr 2012
Location: France
Posts: 993
Rep Power: 24
Yann will become famous soon enough
Hi Pere,

This is a matter of visualization in ParaView.
When loading your case in ParaView, untick the "Decompose polyhedra" option and you should get the same visualization as in the wiki.

katiuskas likes this.
Yann is offline   Reply With Quote

Old   April 15, 2022, 06:27
New Member
Join Date: Mar 2022
Location: Mallorca
Posts: 16
Rep Power: 3
katiuskas is on a distinguished road
Hi Yann,

Thank you very much for the rapid response! I would have never guessed it was a visualization thing. I did not know that some ParaView filters are not good at handling polyhedra and that's why it decomposes them into tetra and pyramids. Good to know!
Yann likes this.
katiuskas is offline   Reply With Quote


mesh, pentagon, refinement

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
problem with parallel calculations OF8 mpirun otaolafr OpenFOAM Programming & Development 0 January 5, 2021 08:33
[snappyHexMesh] Error snappyhexmesh - Multiple outside loops avinashjagdale OpenFOAM Meshing & Mesh Conversion 53 March 8, 2019 10:42
[snappyHexMesh] sHM layer process keeps getting killed MBttR OpenFOAM Meshing & Mesh Conversion 4 August 15, 2016 04:21
[snappyHexMesh] No layers in a small gap bobburnquist OpenFOAM Meshing & Mesh Conversion 6 August 26, 2015 10:38
snappyhexmesh remove blockmesh geometry philipp1 OpenFOAM Running, Solving & CFD 2 December 12, 2014 11:58

All times are GMT -4. The time now is 12:02.