|Reference : Efficient evaluation of the geometrical validity of curvilinear finite elements|
|Scientific congresses and symposiums : Unpublished conference/Abstract|
|Engineering, computing & technology : Multidisciplinary, general & others|
|Efficient evaluation of the geometrical validity of curvilinear finite elements|
|Johnen, Amaury [Université de Liège - ULg > Dép. d'électric., électron. et informat. (Inst.Montefiore) > Applied and Computational Electromagnetics (ACE) >]|
|Remacle, Jean-François [Université Catholique de Louvain - UCL > Institute of Mechanics, Materials and Civil Engineering (iMMC) > Applied mechanics and mathematics (MEMA) > >]|
|Geuzaine, Christophe [Université de Liège - ULg > Dép. d'électric., électron. et informat. (Inst.Montefiore) > Applied and Computational Electromagnetics (ACE) >]|
|Fifth International Conference on Advanced COmputational Methods in ENgineering (ACOMEN 2011)|
|du 14 novembre 2011 au 17 novembre 2011|
|Université de Liège|
|[en] The development of high-order numerical techniques on unstructured grids has been underway for many years. The accuracy of these methods strongly depends of the accuracy of the geometrical discretization, and thus depends on the availability of quality curvilinear meshes.
The usual way of building such curvilinear meshes is to first generate a straight sided mesh. Then, mesh entities that are classified on the curved boundaries of the domain are curved accordingly. Some internal mesh entities may be curved as well. If we assume that the straight sided mesh is composed of well shaped elements, curving elements introduces a kind of "shape distortion" that should be controlled so that the final curvilinear mesh is also composed of well shaped elements.
In this work we propose a method to analyze curvilinear meshes in terms of their elementary jacobians. The method does not deal with the actual generation of the high order mesh. Instead, it provides an efficient way to guarantee that a curvilinear element is geometrically valid, i.e., that its jacobian is strictly positive in all its reference domain. It also provides a way to measure the distortion of the curvilinear element. The key feature of the method is to adaptively expand the elementary jacobians in a polynomial basis, built using Bézier functions, that has both properties of boundedness and positivity.
The algorithm has been implemented in the open-source mesh generator Gmsh, and allows to control the geometrical validity of curvilinear meshes made of triangles, quadrangles, tetrahedra, hexahedra and prisms of any order.
|File(s) associated to this reference|
All documents in ORBi are protected by a user license.