Imposing Dirichlet bc on the boundary of a hole embedded in a continuum body

[dokken]

As you have not supplied your mesh file, it is not very easy to help you. Please not that if you create your mesh with an external mesh creator (say GMSH) you can mark your boundaries in the creation step and import them, as covered in:

Transitioning from mesh.xml to mesh.xdmf, from dolfin-convert to meshio

Your msh file does not contain any physical markers (i.e Physical Line/Volume/Surface), thus you cannot read in cell_data. I would suggest using the following code for reading meshes (from Accessing and marking imported boundaries - #8 by dokken). Note that here I have defined a physical surface and physical curve. import meshio msh = meshio.read("cylinder.msh") def create_mesh(mesh, cell_type, prune_z=False): cells = mesh.get_cells_type(cell_type) cell_data = mesh.get_cell_data("g…