I am using a MATLAB toolbox to generate my mesh for my research in fluid flow through porous rocks. From a stack of binary images I generate a volumetric mesh of the rock.
My problem comes from the exterior surface facets that enclose my volumetric Mesh. The toolbox that I use has the limitation of not being able to generate perfectly flat boundaries for my cuboid mesh. Because of that, I can not mark the boundaries with the usual workflow of:
class Face1(SubDomain): def inside(self, x, on_boundary): return on_boundary and near(x, li) Face1().mark(boundaries, 1)
Can I explicitly identify the Facets which compose each of the six faces of my volumetric mesh of the rock with some other workflow?
I know that I can compute the exterior boundary mesh with the function
but until now I was not able to mark all six faces.
Thanks in advance,