I think I just find the problem. According to "problem with CylindricalManifold<3> after removing cells from a triangulation <https://groups.google.com/g/dealii/c/X1k5Beinw6o/m/GpPzJ5TbBAAJ>", using *set_all_boundary_ids* instead of *set_boundary_id* to set all of its children (line and vertex) to given boundary id.
Besides, you should pay extra attention to the intersection line of the bottom and the side boundary, which should be cylindrical manifold rather than flat manifold. 在2022年10月15日星期六 UTC+8 11:35:09<Chengjiang Yin> 写道： > Hi all, > > I'm trying to make a cylinder with a hole inside it. > > I use the *subdivided_cylinder* to generate a cylinder with 3 slice and > remove the central cell by *create_triangulation_with_removed_cells*. > However, after I setup the boundary and manifold, *the mesh refinement* *does > not seem to work correctly*. > > Here I show the mesh with two global refinements respectively.* It looks > like the manifold is now working for the first refinement* (like flat > manifold) and so the second refinement could not generate as a single > cylinder. > > [image: refine-1.png] > [image: refine-2.png] > > Attached is a demo to generate these two meshes. I'm not sure if this is a > bug of *CylindricalManifold* or I do not use it in the right way. > > Regards, > Chengjiang Yin > -- 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/f15dbdf2-4750-4b0b-940a-934c062afd6bn%40googlegroups.com.
