Collection of abstracts 14th GAMM-Seminar Kiel on Concepts of Numerical Software January 23rd to 25th, 1998.

Friday, January 23rd, 1998

A new flexible implementation og general hp-FEM in 2D based on triangular and quadrilateral elements

A general 2D-hp-adaptive Finite Element (FE) implementation in Fortran 90 is described. The implementation is based on an abstract data structure, which allows to incorporate the full hp-adaptivity of triangular and quadrilateral finite elements.

The h-refinement strategies are based on h2-refinement of quadrilaterals and h4-refinement of triangles. For p-refinement we allow the approximation order to vary within any element. The mesh refinement algorithms are restricted to 1-irregular meshes. Anisotropic and geometric refinement of quadrilateral meshes is made possible by additionally allowing double constrained nodes in rectangles.

The capabilities of this hp-adaptive FE package are demonstrated on various test problems.