Dear Riku, 

this is not (yet) supported by deal.II. You will have to attach manually a 
CylindricalManifold after extrusion. The library does not know yet how to 
compute the extrusion of a general manifold object, therefore it removes all 
manifold objects from the extruded triangulation (as they would not be valid 
objects, being of type Manifold<2>, while you’d require a Manifold<3> object).

L.

> On 20 Sep 2018, at 20:24, Riku Suzuki <rikusuzuki.to...@gmail.com> wrote:
> 
> Hi all
> 
> I've been using the following chunk of code to generate mesh for 2D disk or 
> 3D sphere:
> 
>   template <int dim>
>   void sphere(Triangulation<dim> &tria,
>                       const Point<dim> &center,
>                       double radius)
>   {
>     SphericalManifold<dim> spherical_manifold(center);
>     TransfiniteInterpolationManifold<dim> inner_manifold;
>     GridGenerator::hyper_ball(tria, center, radius);
>     tria.set_all_manifold_ids(1);
>     tria.set_all_manifold_ids_on_boundary(0);
>     tria.set_manifold(0, spherical_manifold);
>     inner_manifold.initialize(tria);
>     tria.set_manifold(1, inner_manifold);
>   }
> 
> The attached file p1.png shows a nice mesh in 2D generated from it (after 
> refinement).
> 
> Now I am trying to extrude the 2D mesh to get a 3D cylinder using 
> extrude_triangulation as follows:
> 
>   void cylinder(Triangulation<3> &tria,
>                         const double      radius,
>                         const double      length)
>   {
>     Triangulation<2> tria2d;
>     Point<2> center(0, 0);
>     sphere(tria2d, center, radius);
>     GridGenerator::extrude_triangulation(tria2d, 10, length, tria);
>   }
> 
> However the SphericalManifold is gone after extrusion, resulting in a mesh 
> shown in p2.png.
> 
> What did I do wrong?
> 
> Thank you
> Riku
> 
> 
> -- 
> 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 dealii+unsubscr...@googlegroups.com.
> For more options, visit https://groups.google.com/d/optout.
> <p1.png><p2.png>

-- 
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 dealii+unsubscr...@googlegroups.com.
For more options, visit https://groups.google.com/d/optout.

Reply via email to