Commun. Comput. Phys., 14 (2013), pp. 1372-1414.


Approximation of H(div) with High-Order Optimal Finite Elements for Pyramids, Prisms and Hexahedra

Morgane Bergot 1, Marc Durufle 2*

1 CALVI project team, INRIA Nancy-Grand Est, Strasbourg, France.
2 BACCHUS project team, INRIA Bordeaux Sud-Ouest, Bordeaux, France.

Received 12 July 2012; Accepted (in revised version) 8 March 2013
Available online 5 July 2013
doi:10.4208/cicp.120712.080313a

Abstract

Classical facet elements do not provide an optimal rate of convergence of the numerical solution toward the solution of the exact problem in H(div)-norm for general unstructured meshes containing hexahedra and prisms. We propose two new families of high-order elements for hexahedra, triangular prisms and pyramids that recover the optimal convergence. These elements have compatible restrictions with each other, such that they can be used directly on general hybrid meshes. Moreover the H(div) proposed spaces are completing the De Rham diagram with optimal elements previously constructed for H^1 and H(curl) approximation. The obtained pyramidal elements are compared theoretically and numerically with other elements of the literature. Eventually, numerical results demonstrate the efficiency of the finite elements constructed.

AMS subject classifications: 65M60

Notice: Undefined variable: pac in /var/www/html/issue/abstract/readabs.php on line 164
Key words: Facet elements, high-order finite element, pyramids, H(div) approximation, De Rham diagram.

*Corresponding author.
Email: morgane.bergot@inria.fr (M. Bergot), marc.durufle@inria.fr (M. Durufle)
 

The Global Science Journal