[LasseJacobsen]

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 (https://transfer.pcloud.com/download...eQCkhTlSpen1Nk).

Meshing is done using the bash script "run_mesh_single.sh".

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!

[LasseJacobsen]

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.