Skip to content

Possible bug in snappyHexMesh for large coordinates

Summary

For large coordinates, snappyHexMesh stops refinement after two levels.

Steps to reproduce

snappyHexMesh stops refinement in directions corresponding to large coordinates. This is a problem when working in geographic coordinates. In the attached casefile, will see that refinement is performed in the x and z directions, but stops after two levels in the y direction. I suspect that this is because the y coordinate is on the order of e+6, while the x direction is e+5, and z is close to zero. I don't know if this is a known issue, I have not found any information about it.

Example case

bad_mesh.tar.xz

I attached a very small casefile. Run blockMesh and snappyHexMesh.

What is the current bug behaviour?

No refinement in the y direction after two levels.

What is the expected correct behavior?

Refinement in all directions, regardless of coordinates.

Relevant logs and/or images

image

Environment information

  • OpenFOAM version : 2006
  • Operating system : CentOS 7
  • Hardware info : Dell Precision 3630
  • Compiler : gcc

Possible fixes

If the case is translated to the origin the problem goes away. This makes me think this is a bug, but I can't point to the part of the code that is causing this. snappyHexMesh finishes without errors.