Hello, Instead of creating a mesh with holes in it. You can use FE_Nothing for the cells that are in a hole. The advantage is that you can use a much simpler mesh. The disadvantage is that the hole is inside the domain. It's not on the boundary of the domain and so it might be a little more complicated to deal with boundary condition at the interface with the hole.
Best, Bruno On Tuesday, March 8, 2022 at 7:46:16 PM UTC-5 [email protected] wrote: > Hi, > I am trying to create a voxel-based mesh by reading TIFF files obtained > from the CT scanning of the specimen. The geometry has a lot of pores, and > of course, these elements should be removed from the triangulation object. > I tried to create first a structured grid using subdivided hyper-rectangle > and then removed the cells which corresponded to pores using the > create_triangulation_with_removed_cells. I managed to do this for one > single image by using parallel::distributed::triangulation. However, this > strategy won't work when I try to read a lot of images and create a 3d > mesh. Because in this case, the coarse mesh itself is quite large. And also > reading from an externally created mesh file is taking forever in this case. > I want to ask if there is a simple way to create a voxel-based mesh for > the parallel::fullydistributed::triangulations and easily remove the cells > which correspond to air? > -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/forum/dealii?hl=en --- You received this message because you are subscribed to the Google Groups "deal.II User Group" group. To unsubscribe from this group and stop receiving emails from it, send an email to [email protected]. To view this discussion on the web visit https://groups.google.com/d/msgid/dealii/400575ab-f685-4567-8772-ca5f2922cc42n%40googlegroups.com.
