Working with Imported Meshes

[dokken]

Instead of reading in the 1D and 2D meshes as meshes, they should be considered as MeshValueCollections. See for instance

Integration over subdomain

The boundary with tag 30 in your mesh is not an interior surface. It is an external surface, as there are only cells on one side of the facets. Thus you can use the ds Measure, without the “+” and “-” restriction to integrate over the surface. For instance consider the minimal code (using your mesh): from dolfin import * import meshio # Convert mesh mesh_from_file = meshio.read("sim_2.msh") def create_mesh(mesh, cell_type, prune_z=False): cells = mesh.get_cells_type(cell_type) cell…