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

[snappyHexMesh] Making surface layers extend to end patches without collapsing

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

LinkBack Thread Tools Search this Thread Display Modes
Old   August 30, 2022, 06:51
Default Making surface layers extend to end patches without collapsing
New Member
Lasse Jacobsen
Join Date: May 2022
Posts: 5
Rep Power: 3
LasseJacobsen is on a distinguished road
I am currently working on meshing a case with internal pipe flow in a simple cylinder.

The mesh is to be used in conjunction with buoyantSimpleFoam, for compressible turbulent flow with buoyancy effects and heat transfer, and later on chtMultiRegionSimpleFoam for two fluid domains and a solid interface.

The mesh is created using snappyHexMesh from the "water.stl" file. See "current_mesh.jpg".

My issue here is making the surface layers extend to the end of the domains (I know that I can do a 2D mesh and simply extrude, but I need my procedure to generalize to more complex geometries in the future).

What I have tried so far:
  • Split the meshing procedure into two steps: (castellate + snap) and (addLayers). The castellation and snapping processes seem to work fine, and the problem I think lies in addLayers phase.
  • Use featureEdgeRefinement with 90 degree extraction to further refine the intersection of the cylinder itself and the end patches.
  • Used refinementSurfaces on the end patches but to no avail.
  • Messed around with minMedialAxisAngle - here I found 110 gives the best results so far.
  • Varying the minThicknessToMedialRatio with constant minMedialAxisAngle, but to no avail.
  • Changing the featureAngle and slipFeatureAngle to no avail.

In general it seems that further refinement of the end patches, in any way, only adds to the issue. Again, what I want is the surface mesh to extend throughout the entire way to the end patches.

I have attached the case file as a link (since I was unable to upload .zip somehow) if anyone would be so kind as to have a look or try and replicate my results (

Meshing is done using the bash script "".

By the way I run OpenFOAM v 2112 (.com distro) on linux subsystems for windows.

Any help would be much appreciated, as I am running out of ideas on what to try. Thank you in advance!
Attached Images
File Type: jpg current_mesh.jpg (79.8 KB, 21 views)
LasseJacobsen is offline   Reply With Quote

Old   August 30, 2022, 08:03
Default Solved
New Member
Lasse Jacobsen
Join Date: May 2022
Posts: 5
Rep Power: 3
LasseJacobsen is on a distinguished road
Solved the issue by changing the background blockMesh, such that and odd number of cells were placed in the radial directions of the cylinder (with the center of the box coinciding with the cylinder's center axis).

Somehow, this worked.
LasseJacobsen is offline   Reply With Quote


addlayers, meshing 3d, meshing cylinder, snappyhexmesh, surface layer

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
Table bounds warnings at: END OF TIME STEP CFXer CFX 4 July 16, 2020 23:44
mass flow in is not equal to mass flow out saii CFX 12 March 19, 2018 05:21
Low Mixing time Problem Mavier CFX 5 April 29, 2013 00:00
Error finding variable "THERMX" sunilpatil CFX 8 April 26, 2013 07:00
Constant velocity of the material Sas CFX 15 July 13, 2010 08:56

All times are GMT -4. The time now is 17:38.