Praveen, we support high order elements also with another technique. You could use MappingQEulerian. On top of a standard (bilinear) triangulation, you define a dim dimensional FESystem of FEQ elements (of the order you like), and generate a vector of FEQ^n displacements of the original mesh.
After you have done this, you can control very finely your Q^n mesh, by simply changing the displacement vector. It requires a little work, but I don't think it would be too difficult to assign the displacements according to a high order mesh file... Luca. -- Luca Heltai <[email protected]> http://people.sissa.it/~heltai/ Scuola Internazionale Superiore di Studi Avanzati Phone: +39 040 3787 449, Office: 732 -- There are no answers, only cross references On Mar 26, 2011, at 4:57 AM, Praveen C wrote: > Thanks. It does look like deal.II cannot read higher order elements. Is there > any plan to add this in future ? > > For cylinder, I can use the boundary descriptions available inside deal.II > but this would become impractical for more complex situations. > > praveen > > On Fri, Mar 25, 2011 at 10:18 PM, Michael Rapson <[email protected]> wrote: > Hi, > > I don't know about support for reading in second order elements as I > understand it deal.II treats the mesh as the 'geometry' without > specifying the FE order, which is done later. It is recommended that > you read in a fairly coarse mesh and then refine a few times in > deal.II since this is faster. I presume your concern is that you need > the boundary of your cylinder to be represented well. The way to do > this is to use a Boundary object to describe the shape you want, then > when you refine in deal.II new nodes are placed correctly on the > boundary. The second mesh in Step-1 shows how this can be done. > > Hope this helps. > Cheers, > Michael > > On Fri, Mar 25, 2011 at 6:16 PM, Praveen C <[email protected]> wrote: > > Hello all, > > I am trying inviscid flow over a cylinder using DG scheme. For good results, > > the boundary must be represented by second order elements. I generate second > > order elements in gmsh. But it looks like deal.ii does not accept second > > order elements in gmsh format. Can you confirm this ? > > If this is so, what other grid options do I have to represent curved > > geometries ? Which grid format should I use ? > > Thanks > > praveen > > _______________________________________________ > > dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii > > > > > > _______________________________________________ > dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii _______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
