Suppose I have a domain as below:
I want to subtract the domain tagged 8 from the whole structure but retain the mesh of the surface 8 and of the rest of the structure as well.
So what I want are these two structures:
- Here the mesh should be the same as above (the mesh is coarser than the above since this was created with different scripts in gmsh)
and
- the innermost domain with the same mesh as the original figure at the beginning.
In FEniCS, I can subtract domains before meshing and mesh them separately. But as soon as I am trying to subtract subdomains after meshing, it’s returning error.
Here is a rough idea what I tried in FEniCS to do this:
- I extracted the the innermost submesh and named it mesh1.
- The whole mesh is named mesh
- mesh2 = mesh - mesh1
This is returning the error that ‘-’ operation is not supported.
Is there a way to do this? Also, please let me know if you have any questions.