opencascade and 2d problems

45 views
Skip to first unread message

Praveen C

unread,
Jan 2, 2016, 5:56:41 AM1/2/16
to Deal.II Googlegroup
Dear all

I am trying to use opencascade for 2-d flow problems, e.g., flow over an airfoil. Learning from step-54, I created an iges file for rae2822 airfoil, a Q1 mesh in gmsh and I can now refine it. The code is here

https://bitbucket.org/cpraveen/deal_ii/src/master/occ/refine/?at=master

The files are in occ/rae2822.

This approach requires using Triangulation<dim,3> since occ works in 3d. If my problem is only in 2d (dim=2), is there a way to use occ without introducing a spacedim template everywhere in my code ? Can I set all the manifold info on a Triangulation<dim,3> object, then get a Triangulation<dim> from that, and use it in my code ?

Thanks
praveen


Luca Heltai

unread,
Jan 2, 2016, 8:39:46 AM1/2/16
to dea...@googlegroups.com
Dear Praveen,

If you do not intend to refine the grid afterwards, then you can use Gridtools::flatten to extract a flattened version of the triangulation, I.e., a two dimensional one after you have done all the refinement you need. 

Best,

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

Praveen C

unread,
Jan 2, 2016, 8:47:16 AM1/2/16
to Deal.II Googlegroup
Dear Luca

Will the manifold info be available for high order mappings after doing flatten ?

Thanks
praveen

luca.heltai

unread,
Jan 2, 2016, 10:05:04 AM1/2/16
to Deal.II Users
No. You have to make all refinement on the <2,3> grid, and then flatten it. Alternatively, you can still use the 2,3 grid as if it was a two dimensional grid… nothing really changes, except the fact that gradients will be three dimensional instead of two dimensional, with third component equal to zero…

Best,

Luca.
Reply all
Reply to author
Forward
0 new messages