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 the sphere layers + the last one for an open space around the sphere (put your boundary condition on it) and this should generate a coarse mesh with cell boundary sitting exactly on the layer boundary. It has already curved manifolds applied, so the refinement should be accurate...
среда, 21 сентября 2016 г., 18:15:05 UTC+3 пользователь Chenchen Liu написал: > > Hi all, > > > I have to use quadratic shape function (Q2 elements) to do my project > (2D). Previously the code works well with Q1 linear elements, when I want > to input the gird with Q2 elements generated by Gmsh, it gives the > following error: > > An error occurred in line <1418> of file <../source/grid/grid_in.cc> in > function > > void dealii::GridIn<2, 2>::read_msh(std::istream &) [dim = 2, spacedim > = 2] > > The violated condition was: > > false > > The name and call sequence of the exception was: > > ExcGmshUnsupportedGeometry(cell_type) > > Additional Information: > > The Element Identifier <8> is not supported in the deal.II library when > reading meshes in 2 dimensions. > > > And, it says the supported quad element should have 4 nodes and 4 edges. > Just want to check whether deal.ii support Q2 elements by Gmsh. > > My mesh looks like the following > > > <https://lh3.googleusercontent.com/-Lk96g0ZRgFg/V-KjPDZp3cI/AAAAAAAAAKs/NZpvWhx06Y4uL6x-6b1gLVfUAayRY5p2gCLcB/s1600/11.png> > > > > If deal.ii does not support higher order element by Gmsh, what other > format is recommended? Thanks! > > > Best, > > Chenchen > > > > -- 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 [email protected]. For more options, visit https://groups.google.com/d/optout.
