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.
