The Approximation of Solutions of Elliptic Boundary-Value Problems via the p-Version of the Finite Element Method

Abstract

The approximation theory developed in [11] is used to determine the piecewise polynomial approximability of solutions of elliptic problems on polygonal domains in $\mathbb{R}^2 $ and polyhedra in $\mathbb{R}^3 $. From these estimates, convergence orders for the p-version of the finite element method applied to such problems are readily obtained. The critical issue is the approximation of the singularities which occur at the nonsmooth parts of the domain boundaries. It is seen that the estimates of [11] involving the weighted Sobolev spaces $Z_l^s $ are well suited for treating such singular functions, yielding directly the optimal approximation degree, up to an arbitrarily small $\varepsilon $.Numerical results for two problems from two-dimensional linear elasticity are also presented. The computations show that the predicted order of convergence is achieved even for low values of p. Moreover, in contrast to the usual h-version of the finite element method, the point at which the p-version enters the asym...