snappyMultiRegionHeaterImplicit does not constrain decomposition

Summary

The heatTransfer/chtMultiRegionFoam/snappyMultiRegionHeaterImplicit does not use faceZones to constrain the decomposition. Hence it might lead to errors in the implicit solution (.. lduAssembly267 .. incorrect number of faces)

Steps to reproduce

Seems to depend on the scotch level. My version does not showcase problem.

Example case

See tutorial case.

Environment information

• OpenFOAM version : v2206

Possible fixes

Use the faceZones generated by snappyHexMesh directly or have the topoSet step generate faceZones instead of currently faceSets.