[deal.II] Re: Input Q2 element generated by Gmsh using Gridin

2016-09-22 Thread Konstantin Ladutenko
Chenchen, If you just need to simulate a multlayer sphere (or a part of it) and the bounding box should not have to be a box (e.g. a free space problem) you can also take a look at my attempt for this here https://groups.google.com/d/msg/dealii/pvqNQUM3eyM/ftqQK1ypBQAJ Just set all raddii of t

Re: [deal.II] Re: Input Q2 element generated by Gmsh using Gridin

2016-09-22 Thread Bruno Turcksin
Chenchen, 2016-09-22 10:54 GMT-04:00 Chenchen Liu : > Thank you for your sharing. To my mind, after we generate Q2 mesh in Gmsh, > we have all the information about vertices and cells. We can read these > information to deal.ii as point list and cell list, and then generate mesh > like step-14 and

Re: [deal.II] Re: Input Q2 element generated by Gmsh using Gridin

2016-09-22 Thread Chenchen Liu
Thank you for your sharing. To my mind, after we generate Q2 mesh in Gmsh, we have all the information about vertices and cells. We can read these information to deal.ii as point list and cell list, and then generate mesh like step-14 and 49. But still, I am not sure whether deal.ii supports hi

[deal.II] Re: Input Q2 element generated by Gmsh using Gridin

2016-09-22 Thread Chenchen Liu
I see your point. Thank you! But it seems that hanging nodes is really annoying :) I will have a try. 在 2016年9月22日星期四 UTC-4上午7:38:35,Konstantin Ladutenko写道: > > BTW, the element in bottom left corner is too triangle for good FEM. To > solve this providing an additional element inside the curved

[deal.II] Re: Input Q2 element generated by Gmsh using Gridin

2016-09-22 Thread Chenchen Liu
Thank you! I almost done with this. 在 2016年9月21日星期三 UTC-4下午12:04:19,Bruno Turcksin写道: > > Chenchen, > > On Wednesday, September 21, 2016 at 11:56:31 AM UTC-4, Daniel Arndt wrote: >> >> We only support to read in meshes that are not refined. Translated to >> your situation this means that you are

[deal.II] Re: Input Q2 element generated by Gmsh using Gridin

2016-09-22 Thread Konstantin Ladutenko
BTW, the element in bottom left corner is too triangle for good FEM. To solve this providing an additional element inside the curved boundary looks to be a common approach, you can find more details here http://dealii.org/developer/doxygen/deal.II/step_6.html#Abettermesh среда, 21 сентября 201

Re: [deal.II] Re: Input Q2 element generated by Gmsh using Gridin

2016-09-21 Thread Praveen C
Having curved boundaries alone is not sufficient sometimes. One would like to use curved elements in the interior also. It would be very useful if one could read high order meshes, e.g., from gmsh. The extra nodes can be used to create an euler vector that can be used in a MappingFEField. There ha

[deal.II] Re: Input Q2 element generated by Gmsh using Gridin

2016-09-21 Thread Daniel Arndt
Chenchen, We only support to read in meshes that are not refined. Translated to your situation this means that you are only able to use Q1 meshes as initial grid. Of course, you can then use triangulation.set_manifold to define the correct behavior when you refine your mesh within dealii along

[deal.II] Re: Input Q2 element generated by Gmsh using Gridin

2016-09-21 Thread Bruno Turcksin
Chenchen, On Wednesday, September 21, 2016 at 11:56:31 AM UTC-4, Daniel Arndt wrote: > > We only support to read in meshes that are not refined. Translated to your > situation this means that you are only able to use Q1 meshes as initial > grid. > Of course, you can then use triangulation.set_ma