[thiagomarinho]

I've been a follower of this forum for a long time and found many answers to several issues here.

I've been through many threads with snapping problems and tried a bunch of different 'solutions' with no success. I've ended up with a much simplified case in order to iterate more quickly, still without success.

I am using openfoam 4 on ubuntu 16.04. The attached picture shows how the mesh is failing to snap to the feature edges of a simple box.

Under snappy log there's mention to the feature edges:

Reading features.
Read edgeMesh box.eMesh
    points      : 8
    edges       : 12
    boundingBox : (0 0 0) (0.5 0.5 0.5)

But during the morph section when its supposed to snap this is the status it provides is:

Feature analysis : total master points:638 attraction to :
    feature point   : 0
    feature edge    : 0
    nearest surface : 0
    rest            : 638

Along with the error/warning:

--> FOAM Warning : Displacement (-0.0001886980013 0.0003330792974 -5.659949559e-05) at mesh point 45330 coord (0.03065552061 0.002054266242 0.01850754925) points through the surrounding patch faces

I've observed this warning in the tutorial run as well but unlike my case the mesh snapped to the edges.

So far I've tried:

• Increasing the snapControls tolerance up to 10 in 0.5 to 2 increments
• implicit feature recognition
• explicit feature recognition
• disabling castellating by having the surface mesh refinement at level 0
• converting the eMesh feature file to .obj format before running snap
• upgrading to ubuntu 18.04 and openfoam 6
• Increasing the refinement level on the curve up to level 7
• refining the base mesh as a whole before running snap
• providing a much coarser base mesh and refining only on the body really stressing the castellation
• simplified geometry (as the box on the picture) to provide reduced number of feature lines
• extracting the feature line file from the gz format before running SHM
• increasing snap iterations up to 200
• using the snapEdge utility post SHM

All of these without success.

Right now it seems I'm on a dead end where the only 'solution' is to really refine the edge mesh so that the absolute deviation of the mesh cells are much smaller, but this is going to be time consuming and not guaranteed to provide sucessfull meshing of the layers.

Im sure that I'm missing something obvious here, could you please provide some input?

I've attached the Log and dictionary files for SHM.

kind regards, T