Hello Daniel,
I have fixed the problem using Your proposed algorithm of mesh creation.
The strange behaviour of velocity-pressure field was really caused by the
bad quality of the mesh (possibly doubled vertices at the boundary between
merged entities).
I want to ask a question, which could be interesting for other users of
dealii-
do we have some bulletproof procedure, how to produce (not so simple)
geometries,
meshes in gmsh and then merge them "safely" in the dealii? Safely meaning-
without
the creation of the doubled vertices at the boundary between merged
geoetries, meshes?
Is there some procedure in dealii to cope with this problem?
The function merge_triangulation requieres "matching" triangulation, could
it be possible
to define a function, which would somehow merge vertices at the boundary of
two not-so-nice
matching trianguations (possibly within some toleration, i.e. if two
vertices are less then \eps
away from each other, than merge them....)
Does my idea make sense?

Thank You

Marek

2016-09-26 15:29 GMT+02:00 Daniel Arndt <d.ar...@math.uni-goettingen.de>:

> Marek,
>
> the attached file should be close to what you want to do.
> Basically, you need to create two geometries: One for the cylinder shell
> and one for the cylinder.
> Here, we take advantage of the fact that GridGenerator::cylinder always
> creates four subdivisions w.r.t height.
> Furthermore, the cylinder needs to be refined once to match the vertices
> of the cylinder shell. Since we can only merge
> coarse triangulations, we need to flatten the resulting triangulation.
> If you want to create a gemoetry with a different height, you can delete
> cells in the cylinder triangulation using GridGenerator::create_
> triangulation_with_removed_cells
> or add more cells by merging another cylinder.
>
> 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 dealii+unsubscr...@googlegroups.com.
> For more options, visit https://groups.google.com/d/optout.
>

-- 
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