Nonstandard finite element de Rham complexes on cubical meshes arxiv: v1 [math.na] 12 Apr 2018

Transcription

1 Nonstandard inite element de Rham complexes on cubical meshes arxiv: v1 [math.na] 12 Apr 218 Andrew Gillette aibo Hu Shuo Zhang Abstract We propose two general operations on inite element dierential complexes on cubical meshes that can be used to construct and analyze sequences o nonstandard convergent methods. The irst operation, called DoF-transer, moves edge degrees o reedom to vertices in a way that reduces global degrees o reedom while increasing continuity order at vertices. The second operation, called serendipity, eliminates interior bubble unctions and degrees o reedom locally on each element without aecting edge degrees o reedom. These operations can be used independently or in tandem to create nonstandard complexes that incorporate Hermite, Adini and trimmed-adini elements. We show that the resulting elements provide convergent, non-conorming methods or problems requiring stronger regularity and satisy a discrete orn inequality. We discuss potential beneits o applying these elements to Stokes, biharmonic and elasticity problems. 1 Introduction To obtain compatible and stable numerical schemes in the ramework o discrete dierential orms and inite element exterior calculus [5, 6, 25], physical variables are discretized in discrete dierential complexes. There have been many successul inite element de Rham complexes, consisting o the Lagrange, 1st and 2nd Nédélec, Raviart-Thomas, Brezzi-Douglas-Marini (BDM) and discontinuous elements. On simplicial meshes, the construction o these elements is classical [33, 32, 37, 13]. On cubical meshes, the analogue o the Raviart-Thomas amily has a tensor product nature [9], while the analogue o the BDM amilies, denoted as S r Λ k, are developed in [2] as a resolution o the serendipity element (c.. [3]). These results are summarized in the Periodic Table o the Finite Elements [7]. The recently deined trimmed serendipity amily, denoted S r Λ k, provides a distinct resolution o the serendipity element as subcomplex o the de Rham sequence with some computational advantages [23]. While the de Rham complex describes regularity spaces appropriate or Hodge-Laplacian problems, it is not suited or low problems that involve the same dierential operators but require stronger regularities. A prominent example is the Stokes problem and the corresponding Stokes complex, Department o Mathematics, University o Arizona, 617 N. Santa Rita Ave., Tucson, Arizona, USA. agillette@math.arizona.edu Department o Mathematics, University o Oslo, PO Box 153 Blindern, NO 316 Oslo, Norway. kaibohu@math.uio.no LSEC, Institute o Computational Mathematics and Scientiic/Engineering Computing, Academy o Mathematics and System Sciences, and National Centre or Mathematics and Interdisciplinary Sciences, Chinese Academy o Sciences, Beijing 119, People s Republic o China szhang@lsec.cc.ac.cn 1

2 given by R H 2 (Ω) [ H 1 (Ω) ] 2 L 2, in two space dimensions. The importance o precisely preserving mass conservation (also called the ergence-ree condition) has been shown in [29]. Accordingly, inite dimensional subcomplexes o the Stokes complexes can be used to construct suitable elements. Falk and Neilan [21] construct a sequence starting with the C 1 Argyris element on 2D triangular meshes. Neilan and Sap [35] give an analogous construction on cubical meshes, which is a resolution o the Bogner-Fox-Schimt (BFS) C 1 element. Tensor product constructions based on splines in the context o the isogeometric analysis are proposed in [14]. We also reer to [34, 35] and the reerence therein or relevant work in 3D and reer to [43, 8, 16] or constructions on macroelements. To obtain these inite element complexes with higher regularity, one usually needs a high polynomial degree and extra derivative degrees o reedom, or certain mesh conditions. An attractive approach to constructing Stokes elements in a simple ashion while ensuring mass conservation is to use non-conorming elements. For instance, the discrete Stokes complex R Morley h [CrouzeixRaviart] 2 h DG, (1.1) consists o the Morley and the vector Crouzeix-Raviart elements and piecewise derivatives h and h. Although the Morley element is an elegant non-conorming choice or the biharmonic problem, it is not convergent or the Poisson equation (c.. [41, 36]). Correspondingly, the vector Crouzeix-Raviart element is not desirable as an H() element or the Hodge Laplacian problem. For the Darcy-Stokes- Brinkman problem, which is a perturbation o the low problem with an H() Hodge Laplacian, the numerical experiments in [31] demonstrate that the convergence with the Crouzeix-Raviart element deteriorates as the coeicient in ront o the Laplace operator tends to zero. As argued in [31, 39], a possible cure or this problem is to seek elements that are conorming in H() and non-conorming convergent in [ H 1] n, then urther construct discrete sequences that are conorming to the de Rham complex and non-conorming convergent as a Stokes complex. A similar complex is constructed on quadrilateral grids [44]. On a given mesh, any approximation space can be characterized by local unctions on each cell and the continuity across the boundaries o the cells, which we shall call the inter-element continuity in subsequent discussions. Finite elements are a special kind o approximation in which inter-element continuity can be ensured via constraints local to each element. Precisely, in Ciarlet s classical deinition o inite elements, a inite element is denoted as a triple (P,G,D) where P is the approximation space (in this case, a space o polynomial dierential orms), G is the geometry (in this case, an n-cube) and D is the set o unctionals on P that are associated to portions o G (e.g. evaluation at a speciic vertex o G, or integration against a speciic test unction along an edge o G). The inter-element continuity is imposed by requiring that these DoFs are single valued on neighboring elements. This is an extra locality condition and in general, approximation spaces, such as some spline spaces, may not ulill this condition. In the theory o inite element systems (FES) [17], another ramework o inite elements developed by Christiansen and collaborators, no degrees o reedom are explicitly used. The approximation space is viewed as a collection o spaces living on all subcells (o all dimensions) o the mesh. The inter-element continuity is imposed by requiring that all the traces o unctions on higher dimensional cells coincide with the unctions on lower dimensional cells (thereore the global unctions 2

3 are single valued). The original theory o the FES is tailored or inite element dierential orms [17], and is generalized to spaces with higher continuity in [16]. In general, there are two ways to reduce a piecewise-deined approximation space: reduce the size o local shape unction spaces or increase the inter-element continuity. These two operations can always be done. However, rom the perspective o inite elements, we require the irst operation (local reduction) to preserve certain local approximation properties(e.g. containment o polynomials o a certain degree) and require the second (global reduction) to retain the locality o the approximation (e.g. inite elements deined in Ciarlet s sense or in the FES sense). From a viewpoint o approximation theories, the locality o the approximation spaces plus the local approximation properties guarantees the global approximation properties. In light o the inite element exterior calculus, one urther demands that these two operations also preserve the exactness o a sequence o inite element spaces. In this paper, we provide a canonical procedure or modiying conorming derham subcomplexes into non-conorming convergent sequences o spaces o elements that can be used in problems requiring greater regularity. In 2D, our construction leads to simple elements or the Darcy-Stokes-Brinkman problem while in 3D our construction yields the Adini plate element [1] and an H() element that is conorming or Poisson/Stokes type problems. Instead o enriching an existing space with bubble unctions as in, e.g. [31, 39], we reduce the local and global degrees o reedom o standard de Rham inite element complexes to sequences that provide the desired properties. From the perspective o nonconorming methods, consistency is preserved in the space reduction. More precisely, i the consistency error in the Strang lemma converges to zero as h or a space V h, then the consistency error or any subspace o V h also tends to zero. The ideas presented here are related to recent investigations into nonstandard vector elements. Stenberg introduced a amily o nonstandard elements or the H() space on triangles and tetrahedra in [38]. These elements are subspaces o the classical BDM spaces that are continuous at vertices o the mesh. Enorcing this continuity constraint reduces the number o global degrees o reedom (DoFs). The canonical interpolations deined by the DoFs, however, do not commute with the dierential operations; instead, Stenberg uses macroelement techniques to show the in-sup condition holds [38]. In [15], Christiansen, Hu, and Hu derive complexes o nonstandard elements rom a dierent perspective. The elements are constructed as variations o copies o Lagrange and Hermite elements, and as a result, canonical nodal bases are available in some cases. The work is partly motivated by an element or linear elasticity [28, 26]. The extended periodic table in [15] includes the triangular and tetrahedral Stenberg elements as a special case or the H() space and several new H() inite elements. The results presented in this paper can be regarded as a continuation o the discussions in [15]. In both the simplicial and the cubical cases, the number o global DoFs is reduced. However, there are some prominent dierences. The triangular Hermite element does not lead to convergent schemes or biharmonic problems, although modiications on the shape unctions are possible to make it convergent [42]. In contrast, cubical elements with vertex derivative DoFs usually do lead to convergent schemes or the biharmonic equation due to the geometric symmetry. The Adini plate element serves as an illuminating example: it is globallyc with derivative DoFs at vertices and comparableto the Hermite element on triangular meshes (c.. [27] and the reerences therein). In this paper, we will construct the Adini element by a new procedure, and place it within an entire inite element complex on squares or cubes. Our approach to element construction is to provide two operations that can be applied to known 3

4 inite element complexes: a DoF-transer operation and a serendipity operation. The DoF-transer operation moves some edge DoFs to vertex DoFs. The local shape unction spaces and bubbles do not change and the DoF-transer preserves unisolvence and exactness. The serendipity operation eliminates some interior bubbles and DoFs at the same time in a way that preserves unisolvence and exactness. The idea o the serendipity operation bears some similarity to the serendipity reduction process in virtual element methods as described in [18, 19], however, virtual element techniques are chiely concerned with applicability to many-sided polygonal and generic polyhedral element geometries. Combining these two operations, we reduce the local and global spaces o various known discrete complexes on cubical meshes [2, 4, 23] to obtain Hermite, Adini and trimmed-adini type amilies. The H() element in the 2D Hermite complex coincides with the nonstandard element introduced in Stenberg [38]. The rest o this paper will be organized as ollows. In Section 2, we introduce notation and background material. In Section 3, we deine the DoF-transer and serendipity operations. In Section 4, we deine the Hermite, Adini, and trimmed-adini elements and complexes. Section 5 is devoted to the approximation and convergence properties o the new inite elements, including the discrete orn inequality and the convergence as non-conorming elements. In Section 6, we urther use the idea o the serendipity operation to construct a complex o minimal elements that are conorming as a de Rham complex and non-conorming as a Stokes complex. We conclude in Section 7 with some remarks and uture directions. 2 Notation and background Let Ω R n be a polyhedral domain and T h be a mesh on Ω consisting o n-cubes. Sometimes we use Th 2 and T h 3 to denote 2D and 3D meshes respectively. A single square element is denoted by 2 and a single cube element by 3. Given a mesh, we use V, E, F and to denote the sets o vertices, edges, aces and three dimensional cells (zero to three dimensional cells) respectively. We use ν and τ to denote the unit normal and tangential vectors o a cell, respectively, when applicable. We use P r,s (), or simply P r,s, to denote the polynomial space o degree r in x and degree s in y on a domain, i.e. P r,s := P r (x) P s (y). Similarly, in three dimensions, we deine P r,s,t := P r (x) P s (y) P t (z), When the orders are equal, we get the tensor product polynomial spaces, denoted Q r := P r,r, and Q r := P r,r,r. Following [2], we deine the superlinear degree o a monomial to be the total degree o actors that enter with quadratic or higher degrees. For example, the superlinear degree o x 2 y 3 z is 5, where we count the degrees o x and y, but not z. We use S r to denote the polynomials with superlinear degree less than or equal to r. There is a simple nesting o spaces given by P r S r Q r. We review some basic acts rom homological algebra; urther details can be ound, or instance, in [1, 2, 5]. A dierential complex is a sequence o spaces V i and operators d i such that V1 d 1 V2 d 2 d n1 Vn d n, (2.1) 4

5 satisying d i+1 d i = or i = 1,2,,n1; typically, this condition is denoted as d 2 =. Let ker(d i ) be the kernel space o the operator d i in V i, and Im(d i ) be the image o the operator d i in V i+1. Due to the complex property d 2 =, we have ker(d i ) Im(d i1 ) or each i 2. Furthermore, i ker(d i ) = Im(d i1 ), we say that the complex (2.1) is exact at V i. At the two ends o the sequence, the complex is exact at V 1 i d 1 is injective (with trivial kernel), and is exact at V n i d n is surjective (with trivial cokernel). The complex (2.1) is called exact i it is exact at all the spaces V i. I each space in (2.1) has inite dimension, then a necessary (but not suicient) condition or the exactness o (2.1) is the ollowing dimension condition: n (1) i dim(v i ) =. i=1 For inite element spaces, the exactness o the local shape unction space on each mesh element is called local exactness, while that o the whole space, with the inter-element continuity taken into account, is called global exactness. In subsequent discussions, the local exactness ollows rom classical results. Thereore we mean global exactness as a deault. Another complex W1 d 1 W2 d 2 d n1 Wn is called a subcomplex o (2.1) i or each i = 1,2,,n, W i is a subspace o V i and the dierential operators are the same. Deine H(d,Ω) := { u L 2 (Ω) : du L 2 (Ω) }, where d is a dierential operator. We use u s,ω to denote the H s norm o u on Ω. When there is no possible conusion, we also use s instead o s,ω and use instead o the L 2 norm. There are two de Rham sequences with the above type o Sobolev regularity in 2D: and R H(,Ω) R H(grad,Ω) H(,Ω) grad H(rot,Ω) L 2 (Ω), (2.2) rot L 2 (Ω). (2.3) Here, the 2D operator maps a scalar unction to a vector valued unction, deined by u := ( 2 u, 1 u) T. (2.4) One can obtain (2.3) by rotating the spaces in (2.2) by π/2. Accordingly, in the remainder o this paper, we only consider (2.2) in 2D. There is only one 3D version o the sequence, which reads: R H(grad,Ω) grad H(,Ω) H(,Ω) L 2 (Ω). (2.5) Here, the 3D has its usual meaning, mapping a vector valued unction to a vector valued unction: u := ( 2 u 3 3 u 2, 3 u 1 1 u 3, 1 u 2 2 u 1 ) T. (2.6) When there is no possible conusion, we will use to denote both (2.4) and (2.6), as determined by the spatial dimension. The sequences (2.2), (2.3) and (2.5) are all exact on contractible domains. We 5

6 will work with three known inite element subcomplexes o the de Rham complex. In 2D these are: R Q r Λ (T 2 h ) Q r Λ 1 (T 2 h ) Q r Λ 2 (T 2 h ), R Sr+2 Λ (Th 2 ) Sr+1 Λ 1 (Th 2 ) Sr Λ 2 (Th 2 ), (2.7) R S r Λ (T 2 h ) S r Λ 1 (T 2 h ) S r Λ 2 (T 2 h ), namely, the 2D Raviart-Thomas (tensor product) type inite element complex [4], the 2D BDM type complex [2], and the 2D trimmed serendipity complex [23], respectively, or r 1. The 3D amilies can be similarly denoted by R Q r Λ (T 3 h ) grad Q r Λ 1 (T 3 h ) Q r Λ 2 (T 3 h ) Q r Λ 3 (T 3 h ), R Sr+3 Λ (Th 3 ) grad Sr+2 Λ 1 (Th 3 ) Sr+1 Λ 2 (Th 3 ) Sr Λ 3 (Th 3 ), R S r Λ (Th 3 S r Λ 1 (Th 3 S r Λ 2 (Th 3 S r Λ 3 (Th 3. (2.8) Some problems involve de Rham complexes with enhanced smoothness. For example, or Stokes lows, the velocity belongs to [ H 1 ] n while the pressure is in L 2 with vanishing integration, and : [ H 1 ] n L 2 (Ω) is onto. In this case, we consider the 2D complex and the 3D version R H 2 (Ω) R H 2 (Ω) grad H 1 (,Ω) [ H 1 (Ω) ] 2 [ H 1 (Ω) ] 3 L 2 (Ω), (2.9) L 2 (Ω), (2.1) where H 1 (,Ω) := { u [ H 1 (Ω) ] 3 : u [ H 1 (Ω) ] 3 }. The subsequent discussions involve another 3D smooth de Rham complex R H 2 (Ω) grad [ H 1 (Ω) ] 3 H(,Ω) The sequence (2.11) is exact on any contractible domain Ω. L 2 (Ω). (2.11) 3 Serendipity and DoF-Transer Operations In this section, we propose two operations or reducing piecewise-deined approximation spaces on cubical elements: a serendipity operation, which reduces local bubble unctions, and a DoF-transer operation, which enhances inter-element continuity by moving degrees o reedom. All the abovementioned properties are achieved, including preservation o the exactness. The serendipity operation, denoted by Q, eliminates some interior DoFs and corresponding shape unctions on inidual elements. Figure 1 shows an example o the serendipity operation applied to the tensor product element Q 3 Λ ( 2 ), resulting in the serendipity element S 3 Λ ( 2 ). On -orms, 6

7 Q Figure 1: An instance o the serendipity operation, taking Q 3 Λ ( 2 ) to S 3 Λ ( 2 ). T T Figure 2: An instance o the DoF-transer operation, taking Q 3 Λ ( 2 ) to G3 Λ ( 2 ) and taking Q 3 Λ1 ( 2 ) to G3 Λ1 ( 2 ). there is only one kind o serendipity operation, but on k-orms with k 1, there are serendipity operations to both the regular and trimmed serendipity sequences; we denote both by Q as it will be clear rom context which one is being used. The DoF-transer operation, denoted by Q, moves two DoFs per edge rom edge-association to vertex-association, interpolating either a directional derivative at the vertex (in the case o -orms) or a coordinate value o the vector ield (in the case o 1-orms). Figure 2 shows an example o each type. The serendipity operation reduces local DoFs while the DoF-transer operation reduces global DoFs and, as we will show, T and Q both preserve unisolvence and exactness. Our subsequent constructions o inite element sequences are summarized in diagrams (3.1) and (3.2). They employ the inite element sequences or Q Λ, SΛ, and S Λ given in (2.7). Q r+2 Λ Q r+2 Λ1 Q r+2 Λ2 T T Gr+2 Λ Q Gr+2 Λ1 Q Gr+2 Λ2 T Q Q S r+2 Λ Q S r+1 Λ 1 Q S r Λ 2 (3.1) T T T A r+2 Λ A r+1 Λ 1 A r Λ 2 In diagram (3.1), we start rom the standard tensor product complex Q Λ (top, back), with the dierential operators as shown. By a serendipity operation Q, we obtain the BDM-type complex SΛ 7

8 (bottom, back). By the DoF-transer operation T, we get the Hermite complex G Λ (top, ront). Applying the serendipity operation to G Λ or DoF-transerring SΛ, we obtain the Adini complex AΛ (bottom, ront). We will show that the diagram (3.1) commutes, i.e. that the operations Q, T and the dierential operators commute with each other. The extension o this diagram to complexes in 3D is straightorward. A variant o this construction is described in diagram (3.2). Here, we start again with Q Λ but use a serendipity operation that results in the trimmed serendipity sequence S Λ. Applying both Q and T gives a trimmed-adini complex, denoted A Λ. As indicated in the diagram, there are some equalities between trimmed and non-trimmed-adini elements when k = or n, based on identities or serendipity-type elements, but or < k < n, the Adini and trimmed-adini elements are truly distinct. Q r+2 Λ Q r+2 Λ1 Q r+2 Λ2 T T Gr+2 Λ Q Gr+2 Λ1 Q Gr+2 Λ2 T Q Q Sr+2 Λ Q Sr+2 Λ1 Q Sr+2 Λ2 (3.2) T T T A r+2 Λ A r+2 Λ1 A r+2 Λ2 = = A r+2 Λ A r+1 Λ 2 We now ormalize the operations T and Q. Degrees o reedom on an n-cube n when P is a space o polynomial dierential k-orms are typically described by u (tr u) q, q I, (3.3) where n is a sub-ace o the cube and I is a space o polynomial dierential (nk)-orms. We call I the index space associated to and there is a canonical association D n d=k n dim = d I (3.4) For instance, DoFs associated to Q r Λk ( n ) have I := Q r1 Λdk () or each n [4, equation (13)]. For inite elements conorming with respect to the standard derham sequence, i.e. (2.2) and (2.5), DoFs are associated to aces o dimension only in the case o -orms (k = ). In those cases, dimi = 1 and (3.3) is interpreted as evaluation o the -orm at the vertex. The non-standard elements considered here have two additional kinds o DoFs. For a -orm u and vertex x, we allow partial derivative evaluation at vertices, i.e. u i u(x) (3.5) 8

9 where i indicates the direction o an edge e i incident to x. For a 1-orm v and vertex x, we allow evaluation o the 1-orm at the vertex, i.e. v v(x), (3.6) where v(x) is a vector in R 2 or R 3. While these DoFs are associated to a vertex o the geometry in the global setting, we observe that they require both a vertex and an incident edge direction to be characterized locally, e.g. the component o the vector v(x) in the 1-orm case. Thereore, we continue to treat these as edge DoFs to simpliy the upcoming ormalism. Let n = 2 or 3. The serendipity operation Q : (P 1, n,d 1 ) (P 2, n,d 2 ) can be thought o as acting component-wise on P 1 and D 1 and should satisy the ollowing conditions: (P 1, n,d 1 ) and (P 2, n,d 2 ) are inite elements, and Q is onto. P 2 P 1 and Q : P 1 P 2 is a projection. Q is the identity on DoFs associated to vertices or edges o n. For max(k,2) d n and each n with dim = d, let I,1 and I,2 denote the index spaces rom (3.4) rom D 1 and D 2, respectively. Then I,2 I,1 and Q : I,1 I,2 is a projection. Q : D 1 D 2 is the linear morphism that is determined by these conditions via (3.4). The DoF-transerring operation T : (P 1, n,d 1 ) (P 2, n,d 2 ) can be thought o as acting component-wise on P 1 and D 1, and should satisy the ollowing conditions: (P 1, n,d 1 ) and (P 2, n,d 2 ) are inite elements, and Q is onto. P 1 = P 2 and T : P 1 P 2 is the identity map. D 1 only has DoFs o type (3.3). T : D 1 D 2 only aects DoFs associated to edges e n with dimi e 2, and leaves all other DoFs ixed. On each such edge e, T changes two DoFs rom type (3.3) to type (3.5) or (3.6), according to whether P 1 is a space o -orms or 1-orms, respectively. Lemma 1. I Q leads to the same local shape unctions or both the DoF-transerred and the original elements, then diagrams (3.1) and (3.2) are commutative. Proo. The result holds because T Q and Q T yield elements with same local spaces and same inter-element continuity. The condition or Q in Lemma 1 is ulilled in all the examples below or the Hermite, Adini and trimmed-adini amilies. We now veriy some properties o the two operations. Lemma 2. The DoF-transer operation preserves unisolvence, i.e. i (P, n,d) deines a unisolvent inite element, then T (P, n,d) is also unisolvent. Proo. The operation T only changes the type o some edge DoFs, rom type (3.3) to type (3.5) or (3.6), but preserves association o DoFs to edges, as discussed ater (3.6). The resulting set o DoFs is still linearly independent on each edge, so unisolvence o the element is preserved by the operation. 9

10 Lemma 3. The DoF-transer operation preserves exactness, i.e. i (P, n,d ) deines an exact sequence o inite elements, then T (P, n,d ) is also exact. Proo. Since the DoF-transer only changes some - and 1-orms, the exactness at T (P k, n,d k ), k 3, holds due to the exactness o the original sequence (P, n,d ). For the exactness at k = 1, we observe that i u T (P 1, n,d 1 ) and du =, then there exists φ (P, n,d ) such that gradφ = u, by the exactness o (P, n,d ) at index 1 and the act that T (P, n,d ) are subspaces o (P, n,d ). Since u is C at vertices, we conclude that φ is correspondingly C 1 at vertices. This implies that φ T (P, n,d ) and shows the exactness o T (P, n,d ) at index 1. It only remains to check the exactness at index 2. This ollows by a dimension count. In n dimensions, on each element, n DoFs are added per vertex and n DoFs are removed per edge or -orms and likewise or 1- orms. Thus, the alternating sum o dimension counts o spaces in the sequence is unaected, meaning exactness is preserved. Thereore the exactness o the whole sequence T (P, n,d ) holds. Remark 1. Starting with 2-orms, the DoF-transerred sequence branches into the standard inite element sequence. A direct check o the exactness at the space o two orms is not trivial. For example, the 2D 1-orms (H() element) with vertex continuity and the 2-orms o piecewise polynomials have been analyzed in Stenberg [38] by a macroelement technique. From the point o view presented in this paper and [15], the exactness seems more clear and natural. In subsequent discussions, we use the term bubble unction to mean a shape unction on n or which all DoFs vanish on all l-dimensional acets o n or l n1. Lemma 4. The DoF-transer operation does not change the local space o bubble unctions on acets o dimension l 2. Proo. Since the DoF-transerred sequence T (P, n,d ) has stronger inter-element continuity than (P, n,d ),weconcludethatthebubblespaceso T (P, n,d )arecontainedinthoseo(p, n,d ). On the other hand, the DoF-transer operation does not change the interior DoFs associated to acets with dimension l 2. So, the spaces o bubble unctions o T (P, n,d ) and (P, n,d ) have the same dimension. Thereore, the spaces o bubble unctions or these two amilies are the same. We now examine properties o the serendipity operation. First, we address the issue o unisolvence. Degrees o reedom o type (3.3) are unisolvent i and only i on each k dimensional ace m,k m n, the trace space with vanishing boundary conditions on (the space o bubble unctions on ) coincides with the local space I used in the deinition o degrees o reedom in (3.3), i.e., tr (P k H Λ k ) = I. Here, the Hodge star is used since in the language o dierential orms, the inner product is deined by (u,v) := u vdx. In terms o vector proxies, the test space I is just the same as the space o bubbles on. We assume that the serendipity operation only aects degrees o reedom o type (3.3). I the image o a the serendipity operation Q : (P 1, n,d 1 ) (P 2, n,d 2 ) satisies the corresponding relation tr (P 2 H Λ ) = I,2, then unisolvence will be preserved. 1

11 For exactness, consider the diagram (3.7): d ker(q ) d ker(q1 ) d ker(qn ) ι ι ι ι R V d V1 d d Vn (3.7) Q Q 1 Q 2 Q 3 R Ṽ d Ṽ1 d d Ṽn The irst row o (3.7) consists o the kernel spaces o the serendipity operation, indicating the subspaces that have been eliminated. The second row is the inite element sequence that we start with and the last row is the reduced sequence. By deinition, each column ker(qk ) ι Vk Q k Ṽk, (3.8) is a short exact sequence, where ι is the inclusion map. Assume that the row in the middle, i.e. R V d V1 d d Vn, (3.9) is exact. From a general algebraic result, the exactness o either the irst or last row o (3.7) implies exactness o the other (c.. [2] or the general result and [17, Proposition 5.16, Proposition 5.17] or an application in the context o inite element systems). In particular, the exactness o the kernel sequence ker(q ) d ker(q1 ) d d ker(qn ), (3.1) implies that o the reduced sequence R Ṽ d Ṽ1 d d Ṽn, (3.11) and vice versa. The ollowing lemma summarizes this result. Lemma 5. In any setting o the orm (3.7), i any two rows o the diagram are exact then the third is exact as well. In particular, exactness o the irst and second rows implies that the serendipity operation preserves exactness. 11

12 Figure3: The lowestorderhermite complex, with regularityh() H() L 2. Thelocal shape unction spaces are Q 3 P 3,2 P 2,3 Q 2, or, equivalently, Q 3 Λ ( 2 ) Q 3 Λ 1 ( 2 ) Q 3 Λ 2 ( 2 ). Remark 2. From the perspective o inite element systems, the serendipity operation reduces local (bubble) spaces on various dimensions and the DoF-transer only changes the spaces carried by zero dimensional cells (vertex spaces) and the trace/restriction rom one dimensional cells to their boundaries. Although a rigorous theory has not been studied in the literature, we observe that the structures o the FES are not changed or two and higher dimensions. The DoF-transer deined in this section changes the system associated with one dimensional cells rom Lagrange to Hermite type. This partly explains the act that ace and interior bubble unctions are not aected by the DoF-transer. 4 Hermite, Adini and trimmed-adini complexes In this section we introduce the Hermite, Adini and trimmed-adini complexes obtained rom standard complexes by the serendipity and DoF-transer operations. We describe the unction spaces and degrees o reedom deinitions in two and three space dimensions. 4.1 Hermite complex Two dimensions The 2D Hermite sequence is denoted by R G r Λ (T h ) G r Λ 1 (T h ) G r Λ 2 (T h ). (4.1) The shape unction spaces coincide with those o the tensor product amily on squares, i.e. Q Λ k, or Raviart-Thomas elements. However, the Hermite amily is only deined when r 3, as smaller r values do not provide enough degrees o reedom to ensure the requisite continuity at vertices or k = and k = 1. The lowest order case is shown in Figure 3. Space G r Λ ( 2 ). The shape unction space o G r Λ ( 2 ) is Q r. The DoFs can be given by unction evaluation and irst order derivatives at each vertex x V: u(x), i u(x), i = 1,2, moments on each edge e E: e uq dx, q P r4 (e), 12

13 interior DoFs in each F: up dx, p Q r2 (). Space G r Λ 1 ( 2 ). The shape unction space o G r Λ 1 ( 2 ) is P r,r1 P r1,r, the same as Q r Λ 1 ( 2 ). The DoFs can be given by unction evaluation at each vertex x V: v(x), moments o the normal components on each edge e E: v ν e q dx, q P r3 (e), interior degrees o reedom in each F: v ψ dx, ψ P r1,r2 () P r2,r1 (). e Space G r Λ2 ( 2 ). The space G r Λ2 ( 2 ) consists o piecewise Q r1 polynomials on each element Three dimensions The 3D Hermite sequence is denoted by R G r Λ (T h ) grad G r Λ 1 (T h ) G r Λ 2 (T h ) G r Λ 3 (T h ). (4.2) Space G r Λ ( 3 ). The shape unction space o G r Λ ( 3 ) is Q r. The DoFs can be given by unction evaluation and irst order derivatives at each vertex x V: u(x), i u(x), i = 1,2, moments on each edge e E: moments on each ace F: interior DoFs in each : uq dx, q P r4 (e), e up dx, p Q r2 (), uw dx, w Q r2 (). 13

14 Space G r Λ 1 ( 3 ). The discrete H() space o the 3D Hermite complex has the shape unction space P r,r+1,r+1 P r+1,r,r+1 P r+1,r+1,r. DoFs or a unction v G r Λ1 ( 3 ) can be given by unction evaluation at each vertex x V: v(x), moments o the tangential components on each edge e E: v τ e q dx, q P r3 (e), e moments o the tangential components on each ace F: (v ν ) p dx, p P r1,r2 () P r2,r1 (), where v ν is treated as a 2-vector in the plane o the ace, interior DoFs in each : v s dx, s P r1,r2,r2 () P r2,r1,r2 () P r2,r2,r1 (). Spaces G r Λ 2 ( 3 ) and G r Λ 3 ( 3 ). The last two spaces in the sequence coincide with the standard Raviart-Thomas element o degree r and the piecewise polynomial tensor product element Q r1, i.e. Q r Λ 2 ( 3 ) and Q r Λ 3 ( 3 ). The DoFs o G r Λ 2 ( 3 ), can be given by (c.. [9]) moments on each ace F: u ν q dx, q Q r1 (), interior DoFs or each : u s dx, s P r2,r1,r1 () P r1,r2,r1 () P r1,r1,r2 (). 4.2 Adini complexes The Adini element or the plate problem [1] can be regarded as a serendipity version o the cubical Hermite element. The shape unction space o the 2D Adini element in the lowest order reads Q 1 +span { x 2 iq : 1 i 2, q Q 1 }, which coincides with the shape unction space o the serendipity element S 3 Λ ( 2 ). The degrees o reedom in this case can be given as the unction evaluation and the irst order derivatives at each vertex x V. We now explain how this deinition can be expanded to an entire exact seqeuence o inite elements. 14

15 +2 +3 Figure 4: The lowest order Adini complex, with regularity H() H() L 2. The local shape unction spaces are S 3 Λ ( 2 ) S 2 Λ 1 ( 2 ) S 1 Λ 2 ( 2 ) Two dimensions The 2D Adini sequence is denoted by R Ar+2 Λ (T h ) Ar+1 Λ 1 (T h ) Ar Λ 2 (T h ), (4.3) The space A r+2 Λ ( 2 ) represents the Adini element o degree r+2 and A r Λ 2 ( 2 ) = S r Λ 2 ( 2 ) = P r, the space o piecewise polynomials o total degree at most r. The lowest order case o the sequence is shown in Figure 4. Space A r Λ ( 2 ). The shape unction space is S r Λ ( 2 ) = P r + span{x r y,xy r }. The DoFs are given by unction evaluation and irst order derivatives at each vertex x V: u(x), i u(x), i = 1,2, moments on each edge e E: interior DoFs in each F: uq dx, q P r4 (e), e up dx, p P r4 (). Space A r Λ 1 ( 2 ). The shape unction space is S r Λ 1 ( 2 ) = P r +span { ( x r+1 y ), ( xy r+1)}, (4.4) which coincides with the BDM element o degree r. The degrees o reedom can be given by unction evaluation at each vertex x V: u i (x), i = 1,2, moments on each edge e E: interior DoFs in each : u ν e q dx, q P r2 (e), e u p dx, p [P r2 ()] 2. 15

16 4.2.2 Three dimensions The 3D Adini sequence is denoted by R Ar+3 Λ (T h ) grad Ar+2 Λ 1 (T h ) Ar+1 Λ 2 (T h ) Ar Λ 3 (T h ). (4.5) Space A r Λ ( 3 ). The shape unction space is the same as that o the serendipity element S r Λ, i.e. polynomials with superlinear degree at most r. For r = 3, the shape unction space can be equivalently represented as Q 1 +span { x 2 i q : 1 i 3, q Q 1}, and A r Λ (T h ) coincides with the 3D Adini element. The degrees o reedom can be given by unction evaluation and irst order derivatives at each vertex x V: u(x), i u(x), i = 1,2,3, moments on each edge e E: moments on each ace F: interior DoFs in each : uq dx, q P r4 (e), e uv dx, v P r4 (), uw dx, w P r6 (). Space A r Λ 1 ( 3 ). The shape unction space coincides with that o the H() inite element given in [2], i.e. yz(w 2 (x,z)w 3 (x,y)) S r Λ 1 () := [P r ()] 3 +span zx(w 3 (x,y)w 1 (y,z)) +gradq(x,y,z), xy(w 1 (y,z)w 2 (x,z)) where w i P r () and q is a polynomial on with superlinear degree at most k+1. The DoFs or a vector unction u can be given by unction evaluation at each vertex x V: moments on each edge e E: moments on each ace F: where u ν is considered as a 2D vector, e u(x), u τ e q dx, q P r2 (e), (u ν ) p dx, p [P r2 ()] 2, 16

17 +5 +6 Figure 5: The lowest order trimmed-adini complex, with regularity H() H() L 2. The local shape unction spaces are S 3 Λ S 3 Λ1 S 3 Λ2. interior DoFs in : u s dx, s [P r4 ()] 3. Thelasttwospacesinthesequencearethesameasthe BDMtypeH() elementandthe piecewise polynomial space in [2]. For completeness, we include the deinitions o these spaces. Space A r Λ 2 ( 3 ). The shape unction space o the cubical BDM space is given by yz(w 2 (x,z)w 3 (x,y)) BDM r := [P r ()] 3 +span zx(w 3 (x,y)w 1 (y,z)). yz(w 1 (y,z)w 2 (x,z)) The DoFs or a vector unction v can be taken as moments on each ace F: interior DoFs in : v ν q dx, q P r (), v p dx, p [P r2 ()] 3. Space A r Λ 3 ( 3 ). Piecewise polynomials o degree r in 3 variables. 4.3 Trimmed-Adini complex The Adini element can also be treated as the irst element in a dierent inite element sequence, which we call the trimmed-adini complex. The shape unction spaces are those o the trimmed serendipity amily, deined in [23] as S r Λ k := S r1 Λ k +κs r1 Λ k+1. The operator κ is the oszul operator, which is discussed in detail in a inite element context in [5]. As shown in [23], the trimmed serendipity spaces nest in between the regular serendipity spaces via the inclusions S r Λ k S r+1 Λk S r+1 Λ k and satisy identities at orm order and n, given by S r Λ = S r Λ and S r Λ n = S r1 Λ n. The associated complexes in 2D and 3D are stated in (2.7) and (2.8). Figure 5 shows the lowest order trimmed-adini complex in two dimensions. 17

18 4.3.1 Two dimensions The 2D trimmed Adini sequence is denoted by R A r Λ (T h ) A r Λ 1 (T h ) A r Λ 2 (T h ). (4.6) Given the identities on the serendipity and trimmed serendipity spaces or -orms and n-orms just mentioned, the sequence could also be written as A r Λ A r Λ 1 A r1 Λ 2. Thus, only the space A r Λ1 ( 2 ) is distinct rom a 2D Adini-type element already deined. Space A r Λ 1 ( 2 ). The degrees o reedom can be given by unction evaluation at each vertex x V: u i (x), i = 1,2, moments on each edge e E: e u n e q dx, q P r3 (e), interior DoFs in : u p dx, p [P r3 ()] 2 H r1 Λ (), where H r1 Λ () denotes the space o (scalar) homogeneous polynomials o degree r1 on Three dimensions The 3D trimmed Adini complex is denoted by R A r Λ (T h ) grad A r Λ 1 (T h ) A r Λ 2 (T h ) A r Λ 3 (T h ). (4.7) Using the identities or -orms and n-orms, we can re-write this sequence as R Ar Λ (T h ) grad A r Λ 1 (T h ) A r Λ 2 (T h ) Ar1 Λ 3 (T h ). (4.8) Thus, we will only describe the 1-orm and 2-orm spaces in 3D. Space A r Λ1 ( 3 ). The DoFs or a vector valued unction u can be given by value at each vertex x V: u(x), moments on each edge e E: u τ e q dx, q P r3 (e), e moments on each ace F: (u ν ) p dx, p [P r3 ()] 2 gradh r1 Λ (), where u ν is considered as a 2D vector, interior DoFs in : u s dx, s [P r5 ()] 3 H r3 Λ 1 (). 18

19 Space A r Λ 2 ( 3 ). The DoFs or a vector unction v can be taken as moments on each ace F: v ν q dx, q P r1 (), interior DoFs in : 4.4 Unisolvence and exactness v p dx, p [P r3 ()] 3 gradh r1 Λ (). The unisolvence o the elements and the exactness o the various complexes introduced in this section ollow rom the properties o the serendipity and the DoF-transer operations shown in Lemma 2 Lemma 5. For instance, since the Hermite and Adini complexes are the images o the T applied to the tensor product and serendipity complexes, respectively (recall the diagram in (3.1)), and since T preserves unisolvence by Lemma 2 and exactness by Lemma 3, it ollows that the Hermite and Adini complexes are are unisolvent and exact. It is possible to prove exacness without applying the DoF-transer to the S Λ sequence. Consider the diagram shown in (4.9), which does not assume any relation between S Λ and A Λ. ker(d,q) ker(d 1,Q) ker(d 2,Q) = = ker(d,g) ι ker(d,g) ι ker(d,g) = ι ι Q r+2 Λ ι Q r+2 Λ1 ι Q r+1 Λ2 T T Gr+2 Λ Q Gr+2 Λ1 Q Gr+1 Λ2 T Q Q S r+2 Λ Q S r+1 Λ 1 Q S r Λ 2 A r+2 Λ A r+1 Λ 1 A r Λ 2 (4.9) By Lemma 2 and Lemma 3, the Hermite amily G Λ is unisolvent and exact. Since the serendipity amily S Λ is known to be exact on contractible domains [2], we know that the kernel sequence ker(d,q) ker(d1,q) ker(d2,q), is exact by Lemma 5. Since the DoF-transer does not change bubble spaces (by Lemma 4) the kernel sequences are in act equal, i.e. ker(d,q) = ker(d,g). Using Lemma 5 again, we conclucde the Adini sequence A Λ is also exact. 19

20 5 Approximation and convergence properties The DoF-transer and serendipity operations aid in establishing approximation and convergence properties. The DoF-transer operation ensures a discrete orn inequality and non-conorming convergence. The serendipity operation preserves these properties since it is a local space reduction that does not change inter-element continuity. 5.1 A discrete orn inequality The orn inequality is an indispensable tool in linear elasticity, asserting that i a vector ield is orthogonal to the rigid body motion, then its H 1 norm can be controlled by its symmetric gradient. Speciically, the orn inequality reads: i u and v are vector valued unctions in a domain Ω, then there exists a positive constant C > such that and u 2 1 C ( ǫ(u) 2 + u 2), u [ H 1 (Ω) ] n, v 2 1 C ( ǫ(v) 2 + n i=1 Ω v i ds ) v 2 1 C ǫ(v) 2, v [ H 1 (Ω) ] n., v [ H 1 (Ω) ] n, Here ǫ is the symmetric gradient mapping a vector to a symmetric matrix, deined by [ǫ(u)] ij := 1/2( i u j + j u i ). For non-conorming inite element spaces, a discrete version o orn s inequalities is desirable. In this case, the derivatives appearing in the original orn s inequality are replaced by the piecewise derivatives and the broken Sobolev norms ( ) 1/2 u m,h :=. T T h u 2 m,t We irst recall the results in [4]. Following the general results, we will see that the DoF-transer operation leads to inite elements satisying the discrete orn inequality. The ollowing two conditions or a inite dimensional space V h are crucial or the discrete orn inequality in [4]. (H1) There exists an integer r 1 such that or any v V h, v T P r (T), T T h. (H2) For any v V h, T T h, let F be an arbitrary (n1) dimensional ace o T. Then v is continuous on a set consisting o at least n points that are not on a common n2 dimensional hyperplane, and the set is aine invariant. Theorem 1 ([4]). Assume that the assumptions (H1) and (H2) are true. Then there exists a constant C independent o the mesh size h, such that the ollowing discrete orn inequalities hold: ( ) v 2 1,h C ǫ(v) 2,T + v 2,Ω, v V h, T T h v 2 1,h C ( T T h ǫ(v) 2,T + 2 n i=1 Ω v i ds ), v V h,

21 and v 2 1,h C T T h ǫ(v) 2,T, v V h, where V h is the subspace o V h with vanishing degrees o reedom on the boundary. In some cases, the condition(h2) is a strong assumption. For example, generally the shape unctions o the Crouzeix-Raviart non-conorming element [11] are only continuous at n 1 points on a n 1 dimensional ace. However, with the DoF-transer operation, the resulting inite element unctions are continuous at each vertex. On an n1 dimensional ace, this amounts to at least n points. Thereore the assumption (H2) is ulilled and we conclude the ollowing result. Theorem 2. The - and 1-orms in the Hermite, Adini and trimmed-adini complexes satisy the discrete orn inequalities in Theorem Convergence as non-conorming elements The DoF-transer operation ensures extra continuity at vertices compared to standard elements, namely, C 1 continuity or -orms and C continuity or 1-orms. This order o continuity allows provable convergence results or problems requiring more regularity in the solution. Speciically, the resulting -orms can be regarded as convergent non-conorming elements or the biharmonic equations [1, 27] while the 1-orms are convergent non-conorming elements or Poisson type equations involving the scalar Laplacian operator. First, recall the Strang lemma (c.. [12]). Lemma 6 (Strang Lemma). Let V and V h be Hilbert spaces and dimv h <. Let a h (, ) be a symmetric positive-deinite bilinear orm on V +V h that reduces to a(, ) on V. Let u V solve where F V V h. Let u h V h solve a(u,v) = F(v), v V, a h (u h,v) = F(v), v V h. Then uu h h in uv h h + v h V h sup w h V h \{} a h (uu h,w h ) w h h. (5.1) The irst term on the right hand side o (5.1) is the local approximation error, which is determined by the local shape unction spaces. The second term, which does not appear in conorming methods, is the consistency error. The consistency error is generally determined by the inter-element continuity, although not strong enough to guarantee conormity. As we shall see, or the new elements proposed in this paper, the consistency error is controlled also because o the proper local shape unction spaces and the geometric symmetry. In particular, we note that i the consistency error o a space V h tends to zero asymptotically as h, then any o its subspaces inherit this property. This means that the space reductions, either by the serendipity or by the DoF-transer operations, do not break the convergence o non-conorming methods. 21

22 In this section, we prove the convergence o the DoF-transerred 1-orms or Poisson equations. Since one oten seeks simple nonconorming elements and one cannot improve the consistency error by increasing local polynomial degrees, we ocus on the lowest order Adini 1-orm A 2 Λ 1. The space A 2 Λ 1 is not a tensor product o copies o scalar elements, meaning that one cannot collect the shape unctions and degrees o reedom o one component to get a unisolvent scalar element. However, one can obtain a scalar element with the same inter-element continuity by adding some interior degrees o reedom. We now show that each component o A 2 Λ 1 (in the above sense) yields a convergent scalar element or the Poisson problem. Then the convergence o each component implies that o the element A 2 Λ 1 or vector Poisson problems. We start the analysis with the ollowing explicit orm o the bases. Lemma 7. The shape unction space o the Adini 1-orm A 2 Λ 1 has a direct sum decomposition A 2 Λ 1 = P 2 Λ 1 J 2 Λ 1 dj 3 Λ, where in three space dimensions xyz J 2 Λ 1 = span x 2 z, y 2 z xyz, xyz xy 2, xz 2 xyz, xyz x 2 y, yz 2 xyz, and dj 3 Λ = span 3x 2 y x 3 y 2 z 2xyz xy 2,, 3x 2 z x 3 y 3 z 3xy 2 z xy 3,, 2xyz x 2 z x 2 y z 3 3yz 2,, 3x 2 yz x 3 z x 3 y z 3 3xz 2,, 3y 2 z y 3 yz 2 xz 2 2xyz,, y 3 3xy 2 yz 3 xz 3 3xyz 2,. Without loss o generality, we ocus on the details o the x-component o the Adini 1-orms in 3D on a reerence cubic = (1,1) 3. We use E x to denote the edges o parallel to the x-axis, and Fx or the aces o parallel to the x-axis. Collecting the shape unctions o the x-component, we deine Equivalently, P x := P 2 ()+span{xyz,x 2 y,x 2 z,yz 2,y 2 z,y 3,z 3,zy 3,z 3 y,x 2 yz}. P x = span{1,x,x 2 }{1,y}{1,z}+span{y 2,z 2,y 2 z,yz 2,y 3,z 3,y 3 z,yz 3 }, where the three sets o y braces indicate one choice is to be made rom each set o braces beore multiplying together. The ollowing degrees o reedom deine a unisolvent inite element with local shape unction space P x: unction evaluation at each vertex x V: moments on edges parallel to the x-axis: e v(x), vdx, e E x, 22

23 moments on aces parallel to the x-axis: interior degrees o reedom on : vdx, F x, vqdx, q P x H 1 (). We note that the above degrees o reedom on aces (including vertices and edges) are those or the x-component o A 2 Λ 1, which are linearly independent or P x and have dimension 16. On the other hand, the space P x has dimension 2. This dierence is due to the non-tensor-product nature o A 2 Λ 1 as explained above. Thereore we also include the interior degrees o reedom with dimension 216 = 4. Lemma 8. The above degrees o reedom are unisolvent or A 2 Λ 1. The local shape unction space P x and the above degrees o reedom give a inite element space on a mesh T h on Ω, which we denote as V h. Deine V h as the corresponding space with homogeneous boundary conditions: { } V h := v V h : v(a) =, a V b ; v =, e E xb, v =, F xb, e where V b is the set o boundary vertices, E xb the boundary edges parallel to the x-axis and F xb the boundary aces parallel to the x-axis. Consider the variational orm o the Poisson equation: given L 2 (Ω), ind u H 1 (Ω), such that ( u, v) = (,v), v H 1 (Ω); (5.2) and its non-conorming element discretization: ind u h V h, such that ( h u h, h v h ) = (,v h ), v h V h. (5.3) Here h is the piecewise gradient deined on T h. To show the convergence o (5.3), we introduce some auxiliary spaces, again using the convention that a list inside y braces indicates one o the options should be chosen: B x = span{ (x 2 1) [ (y 2 1){1,x}+{1,y}(z 2 1) ]}, Q x = span{1,x,x 2 }{1,y}{1,z}; and A x = Q x +B x. The auxiliary space A x is deined so that it is unisolvent with the vertex, edge and aces DoFs o V h, and also orms a C -conorming global space. Denote by Π Q Π B and Π A the interpolations such that Π Q v Q x ; (Π Q v)(a) = v(a), a V(); (Π Q v) = v, e E; x e e 23

24 and Π B v B; x (Π B v) = e e v, F x, Π A v A x ; (Π A v)(a) = v(a); (Π A v) = v, e E; x It is straightorward to check that the interpolations are all well-deined. Lemma 9. For any u C (), it holds that (Π A Π B )(I Π Q )u =. F (Π A v) = v, F x. Proo. For any vertex a o, [(I Π Q )u](a) =, and thereore Π A (I Π Q )u B x ; actually, {v A : v(a) = } = B x. The assertion ollows. Lemma 1. For any v C (), it holds that (v Π A v)dx =, v P x. Proo. Note that Π A Π Q = Π Q, and we have by Lemma 9 that uπ A u = (uπ Q u)π A (uπ Q u) = (uπ Q u)π B (uπ Q u) = (I Π B )(I Π Q )u. By direct calculations, we get Table 1. u y 2 z 2 y 3 y 2 z yz 2 z 3 y 3 z yz 3 Π Q u 1 1 y z y z yz yz uπ Q u y 2 1 z 2 1 y 3 y y 2 z z yz 2 y z 3 z y 3 z yz yz 3 yz Table 1: Interpolations o shape unctions. Denote aces by (Let, Right, Out, Interior, Bottom, Top) F L (1,y,z); F R (1,y,z); (y,z) (1,1) 2 F O (x,1,z); F I (x,1,z); (x,z) (1,1) 2 F B (x,y,1); F T (x,y,1); (x,y) (1,1) 2. Here, or instance, F L (1,y,z) means the ace {x = 1,y,z (1,1)}. For any v P x, it can be veriied that uπ A u = uπ A u, F L F R uπ A u = uπ A u, F I F O uπ A u = uπ A u. (5.4) F T F B Actually, 1. i v Q x, then (I Π A)v = (I Π B )(I Π Q )v = ; 2. i v span{y 2,z 2,y 3,y 2 z,yz 2,z 3,y 3 z,yz 3 } (1) on {F I,F O,F T,F B }, by deinition o Π B, (I Π B )((I Π Q )v) = ; F I,F O,F T,F B (2) on {F L,F R }, (I Π Q )v = (I Π Q )v, or v span{y 2,z 2,y 3,y 2 z,yz 2,z 3,y 3 z,yz 3 } F L F R and w = or w B, thus (I Π B )((I Π Q )v) = (I Π B )((I Π Q )v). F L,F R F L F R 24

25 The assertion o the lemma then ollows. Deine an auxiliary space A h with the local shape unction space A x and the same vertex, edge and ace degrees o reedom as V h. Deine the nodal interpolation Π A h with respect to A h. We use u to denote the jump o u on a ace, i.e. i = T 1 T 2, then u := u T1 u T2. When the context is clear, we omit the subscript and just write u instead. Theorem 3. Let u H 1 (Ω) H 2 (Ω) and u h be the solutions o (5.2) and (5.3), respectively. Then uu h 1,h Ch u 2,Ω. (5.5) Proo. By the Strang lemma, we can estimate the consistency error. We have ( u, h v h )(,v h ) = n uv h = n uv h Π A hv h = n u(v h Π A hv h ) F F T h = 3 xi u(v h Π A h v h)n i = 3 ( xi uc i )(v h Π A h v h)n i, T h i=1 T h i=1 where the second equality is due to the C conormity o A h and the last is by Lemma 1. Here C i is an arbitrary constant. This way, ( u, h v h )(,v h ) 3 in ( xi uc i (v h Π A hv h ) ) Ch u 2,Ω v h 1,h, C i R#(T h ) 3 T h i=1 where #(T h ) is the number o elements in T h. The consistency error is controlled, and the assertion ollows. Now we have proved that the space V h, which consists o the local space P x and the same interelement continuity asthe x-component oa 2 Λ 1, convergesas a non-conormingelement or the Poisson equation. This means that each component (in a proper sense) o A 2 Λ 1 converges as a nonconorming element. However, since we have added interior degrees o reedom to V h, there is one more step to conclude with the convergence o A 2 Λ 1 as a non-conorming element or the vector Poisson equation. Consider the vector Poisson equation: given [L 2 (Ω)] 3, ind u [ H(Ω) ] 1 3, such that ( u, v) = (,v), v [ H 1 (Ω) ] 3 ; (5.6) and its non-conorming element discretization: ind u h Å2Λ 1, such that ( h u h, h v h ) = (,v h ), v h Å2Λ 1. (5.7) Here Å2Λ 1 is the subspace A 2 Λ 1 with vanishing boundary degrees o reedom. Theorem 4. Let u [H 1 (Ω) H2 (Ω)] 3 and u h be the solutions o (5.6) and (5.7), respectively. Then uu h 1,h Ch u 2,Ω. (5.8) Proo. The theorem ollows rom the Strang lemma. In act, the local approximation error satisies the above estimate since piecewise constants are contained in A 2 Λ 1. Moreover, the consistency error also satisies the estimate as we have shown or each component o A 2 Λ 1 with the same local spaces and inter-element continuity. Remark 3. The two dimensional results ollow the same way. 25