Skip to content

BUG: leak closure in snappyHexMesh does not run in parallel

Summary

Leak closure produces very different results in parallel.

Steps to reproduce

Run shm (e.g. with debugFlags (mesh)) non-parallel and parallel side by side. At some point the cell count starts to diverge, e.g.

Removing mesh beyond surface intersections
------------------------------------------
..
Keeping all cells containing inside points
Selected for keeping : 1876028 cells.

Example case

Not easy. Make tiny case with small gaps. Use e.g. hierarchical. Run parallel with debug flags. Run

mpirun -np 9 postProcess -func processorField -parallel

and look at where the outside connects to the inside.

What is the current bug behaviour?

Inside mesh leaks outside but not at gap but at interface between two different processors.

Environment information

  • OpenFOAM version : develop