Efficient assembly of H(div) and H(curl) conforming finite elements
Artikel i vetenskaplig tidskrift, 2009

In this paper, we discuss how to efficiently evaluate and assemble general finite element variational forms on H(div) and H(curl). The proposed strategy relies on a decomposition of the element tensor into a precomputable reference tensor and a mesh-dependent geometry tensor. Two key points must then be considered: the appropriate mapping of basis functions from a reference element, and the orientation of geometrical entities. To address these issues, we extend here a previously presented representation theorem for affinely mapped elements to Piola-mapped elements. We also discuss a simple numbering strategy that removes the need to contend with directions of facet normals and tangents. The result is an automated, efficient, and easy-to-use implementation that allows a user to specify finite element variational forms on H(div) and H(curl) in close to mathematical notation. © 2009 Society for Industrial and Applied Mathematics.

Mixed finite element

Variational form compiler

Piola

Författare

M.E. Rognes

R.C. Kirby

Anders Logg

Chalmers, Matematiska vetenskaper, Matematik

Göteborgs universitet

SIAM Journal of Scientific Computing

1064-8275 (ISSN) 1095-7197 (eISSN)

Vol. 31 6 4130-4151

Ämneskategorier

Matematik

Beräkningsmatematik

DOI

10.1137/08073901X