摘要

In this paper, we present a family of H(div)-compatible finite element spaces on strictly convex n-gons, whose construction makes use of generalized barycentric coordinates. In particular, for integers 0 <= k <= 2, we define finite element spaces with edge degrees of freedom that include polynomial vector fields of order k and whose vector fields have piecewise kth-order polynomial normal traces along the element boundary. These spaces suffer from the shortcoming that the image of the divergence operator includes nonpolynomial functions and, as such, their direct use in a mixed setting along with polynomial scalar fields may lead to unstable discretizations exhibiting degraded or no convergence. We present a general remedy for restoration of optimal convergence that involves "polynomial corrections" of vector fields and their divergence at the element level. These corrections are consistent with one another and require computation of suitable polynomial projection maps. In addition to the theoretical discussion, the performance of the discretization schemes based on the proposed spaces and the accompanying correction maps is numerically evaluated through their application to H(div) eigenvalue problems, mixed and least-squares approximations of linear and nonlinear porous media flow problems.

  • 出版日期2015