Boundary labelling of GMESH-generated meshes

mesh
1

Hello Everyone,
I was following this tutorial, where the blog author generates a three-dimensional mesh making use of the GMESH software, and studies a PDE on such a mesh via the Fenics software.

I tried to reproduce the steps in this tutorial using Fenics 2019.2.0.dev0 on Ubuntu 20.04 64 bit.

At a certain point, the blog author assigns a label to each boundary surface. The first utilised label index is 1.
I managed to successfully reproduce and convert the mesh by using the dolfin-converter.
However, when I go to examine the file geo_facet_region.xml, I see the string value="0" appearing.

I cannot understand why this happens, since the first used index in the boundary labelling process is 1.

Could anyone please explain to me what I did wrong?

Many thanks in advance and best regards.

2

The dolfin-convert function has not been maintained for a long time, and i would suggest you to use meshio to convert your mesh. See for instance: Transitioning from mesh.xml to mesh.xdmf, from dolfin-convert to meshio
And all the examples I’ve posted in this thread for intended usage.

1 Like