Dear Chucui, I would like to extend Daniel's answer slightly.
The source code for channel_with_cylinder is is at about line 2180 of https://www.dealii.org/developer/doxygen/deal.II/grid__generator_8cc_source.html If you want to generate just the grid around the cylinder (i.e., without the channel) then look at how cylinder_tria is set up in lines 2250 - 2305: you can get the grid you posted by adding more shells (and increasing the size of the shell region). The development version of deal.II can copy manifold ids, but 9.0 cannot. Thanks, David On Wednesday, December 19, 2018 at 12:07:16 PM UTC-5, Daniel Arndt wrote: > > Chucui, > > the mesh you are showing degenerates close to the midpoint. This isvery > likely not what you wan't. > For the inner part something like shown in > https://www.dealii.org/developer/doxygen/deal.II/classTransfiniteInterpolationManifold.html > would be a good idea. > Then for the middle part GridGenerator::hyper_cube_with_cylindrical_hole ( > https://www.dealii.org/9.0.0/doxygen/deal.II/namespaceGridGenerator.html#add14cab546d033c1eaacc9234c64ebcd > ) > should work, you just have to make sure that the vertices align when > merging the two meshes via > GridGenerator::merge_triangulations ( > https://www.dealii.org/developer/doxygen/deal.II/namespaceGridGenerator.html#a7cd88e7eacd46697dee80ad2b8438d54 > ). > Make sure too choose copy_manifold_ids if you are using a developer > version or set the manifold ids after creating and merging the whole mesh. > Finally, for the outer part you would just merge some > GridGenerator::subdivided_hyper_rectangle ( > https://www.dealii.org/developer/doxygen/deal.II/namespaceGridGenerator.html#ac76417d7404b75cf53c732f456e6e971 > ). > Again make sure that the verices align and make sure the manifold ids are > correct after merging all the meshes. > > Best, > Daniel > -- 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]. For more options, visit https://groups.google.com/d/optout.
