Hi Jie, One needs to move some of the vertices around. I have some code that does generate this geometry in one of my projects:
https://github.com/kronbichler/adaflo/blob/master/tests/flow_past_cylinder.cc (check the create_triangulation function) The mesh that comes out from this code is accurate but not optimal yet it is not fine enough around the cylinder as compared to be bulk of the domain, at least for some of the benchmark test cases. You can find a picture which gives better convergence results in our recent publication (Fig 11): http://www.sciencedirect.com/science/article/pii/S0021999117306915 However, I do not have that code online right now and it takes too much time to look it up. But it should be easy enough to compose it of several triangulation, adding the transfinite interpolation manifold to the transition from polar->straight elements that we recently introduced in the deal.II and that is mentioned in that paper. Best, Martin On 20.10.2017 01:15, Jie Cheng wrote: > Hi Lucas > > > I am trying to replicate the results here > > <http://www.featflow.de/en/benchmarks/cfdbenchmarking/flow/dfg_benchmark1_re20.html> > of the benchmark problem of flow around a cylinder (circle in 2D). > Here is my code so far: > > > It's been a long time since you posted this question. But if you are > still reading the threads, could you please tell me how you managed to > generate a non-symmetric triangulation in the benchmark case? I tried > to use one hypercube_with_cylindrical_hole and two hyper_rectangles to > construct it, but the result looks ugly upon global refinement. > > Thank you. > > Jie > -- > 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 > <mailto: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.
