Steffen Weißer

Clemens Hofreither, Ulrich Langer, Steffen Weißer

We present a new discretization method for homogeneous convection-diffusion-reaction boundary value problems in 3D that is a non-standard finite element method with PDE-harmonic shape functions on polyhedral elements. The element stiffness matrices are constructed by means of local boundary element techniques. Our method, which we refer to as a BEM-based FEM, can therefore be considered a local Trefftz method with element-wise (locally) PDE-harmonic shape functions. The Dirichlet boundary data for these shape functions is chosen according to a convection-adapted procedure which solves projections of the PDE onto the edges and faces of the elements. This improves the stability of the discretization method for convection-dominated problems both when compared to a standard FEM and to previous BEM-based FEM approaches, as we demonstrate in several numerical experiments.

Steffen Weißer

Only a few numerical methods can treat boundary value problems on polygonal and polyhedral meshes. The BEM-based Finite Element Method is one of the new discretization strategies, which make use of and benefits from the flexibility of these general meshes that incorporate hanging nodes naturally. The article in hand addresses quasi-interpolation operators for the approximation space over polygonal meshes. To prove interpolation estimates the Poincaré constant is bounded uniformly for patches of star-shaped elements. These results give rise to the residual based error estimate for high order BEM-based FEM and its reliability as well as its efficiency are proven. Such a posteriori error estimates can be used to gauge the approximation quality and to implement adaptive FEM strategies. Numerical experiments show optimal rates of convergence for meshes with non-convex elements on uniformly as well as on adaptively refined meshes.

Steffen Weißer

New interpolation and quasi-interpolation operators of Clément- and Scott-Zhang-type are analyzed on anisotropic polygonal and polyhedral meshes. Since no reference element is available, an appropriate linear mapping to a reference configuration plays a crucial role. A priori error estimates are derived respecting the anisotropy of the discretization. Finally, the found estimates are employed to propose an adaptive mesh refinement based on bisection which leads to highly anisotropic and adapted discretizations with general element shapes in two- and three-dimensions.