arxiv: v1 [math.na] 28 Apr 2017

Transcription

1 THE SCOTT-VOGELIUS FINITE ELEMENTS REVISITED JOHNNY GUZMÁN AND L RIDGWAY SCOTT arxiv: v1 [matna] 28 Apr 2017 Abstract We prove tat te Scott-Vogelius finite elements are inf-sup stable on sape-regular meses for piecewise quartic velocity fields and iger (k 4) [2010]65N30, 65N12, 76D07, 65N85 1 Introduction In 1985 Scott and Vogelius [16] (see also [19]) presented a family of finite element spaces in two dimensions wic wen applied to te Stokes problem, produce velocity approximations tat are exactly divergence free In addition, tey proved tat te metod is stable by proving te pair of spaces satisfy te so-called inf-sup condition In teir inf-sup stability proof tey require te mes to be quasi-uniform In addition, te maximum mes size is assumed to be sufficiently small In tis paper we give an alternative proof of te inf-sup stability tat relaxes tese restrictions To be more precise, we prove tat te Scott-Vogelius finite element spaces for polynomial order k 4 are inf-sup stable assuming only tat te family of meses are nondegenerate (sape regular) One key aspect in te new proof is to use te stability of te P 2 P 0 (or te Bernardi-Raugel [2]) finite element spaces As a result te proof becomes significantly sorter In te last paragrap of [6], a modification of teir tecniques is sketced tat provides a proof of te inf-sup stability for te Scott-Vogelius elements tat is different from bot ours and te original proof Recently tere as been interest in developing finite element metods tat produce divergence free velocities or ave better mass conservation properties; see for example [1, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 18, 20, 21] In particular, te review paper [11] discusses in dept te effects of mass conservation in simulations Tere ave also been extensions of te Scott-Vogelius elements to tree dimensions [20, 21, 14], altoug a full general result is still out of reac One difficulty lies in generalizing te concept of singular (or non-singular) vertices to tree dimensions; see [14] To better describe te key differences between te proof of inf-sup stability in tis article compared to te original proof found in [16], we review te proof of [16] Rougly speaking, in [16] given a pressure function p from te discrete space, one wants to find a velocity vector field v from te discrete velocity space suc tat div v is close to p and v H 1 (Ω) C p L 2 (Ω) To do tis, te proof in [16] follows rougly tree steps In te first step a velocity field v 1 is found so tat p 1 = div v 1 p vanises at all vertices At tis step it would be desirable to find a vector field w so tat div w as te same average as p 1 on eac triangle and suc tat div w vanises at te vertices Tis can be done, owever, it is not clear ow to do tis in a stable way Terefore, alternatively in te second step, a continuous piecewise linear pressure function p is introduced 1

2 2 J GUZMÁN AND R SCOTT so tat p 2 = p 1 p as average zero on non-overlapping patces of rougly size K were K is a sufficiently large constant Ten one finds a vector field v 2 so tat div v 2 = p 2 As a result div(v 1 + v 2 ) = p p on Ω Te final step is to find a vector field v 3 so tat div v 3 p, were ere quasi-uniformity and sufficiently small mes size is used Ten one arrives at div(v 1 + v 2 + v 3 ) p on Ω In contrast, we reverse te order of te steps First we find a piecewise quadratic v 1 suc tat p 1 = div v 1 p as average zero on eac triangle Tis can be done in a stable way, using te Bernardi-Raguel finite elements [2] (or te P 2 P 0 finite element space) Ten one finds v 2 tat is piecewise quartic so tat p 2 = div v 2 p 1 vanises at te vertices and as average zero at eac triangle, in wic case we require tat div v 2 as average zero on eac triangle To construct v 2 we combine two basis functions tat are used implicitly in [16] Finally, following [19] a local argument will find v 3 so tat div v 3 = p 2 Hence, div(v 1 + v 2 + v 3 ) = p Te only reason tat we are restricted to k 4 is tat te v 2 constructed above is piecewise quartic In a subsequent paper we sow tat if we impose furter mild restrictions on nonsingular vertices ten we can coose v 2 to be piecewise cubic and, ence, prove inf-sup stability in te case k = 3 Te paper is organized as follows: In te following section we begin wit Preliminaries In Section 3 we prove te inf-sup stability for k 4 2 Preliminaries We assume Ω is a polygonal domain in two dimensions We let {T } be a family of nondegenerate (sape regular) triangulations of Ω; see [3] Te set of vertices and te set of internal edges are denoted by S ={x : x is a vertex of T }, E ={e : e is an edge of T and e Ω} For every z S te function ψ z is te continuous, piecewise linear Lagrange basis function corresponding to te vertex z Tat is, for every y S { 1 if y = z, (21) ψ z (y) = 0 if y z We also define te internal edges and triangles tat ave z S as a vertex: Finally, we define te patc E (z) ={e E : z is a vertex of e}, T (z) ={T T : z is a vertex of T } Ω (z) = T T (z) Te diameter of tis patc is denoted by z = diam(ω (z)) T

3 INF-SUP STABILITY OF UNFITTED STOKES ELEMENTS 3 Figure 1 Example of singular vertices z Dased edges denote boundary edges In order to define te pressure space we need to define singular and non-singular vertices Let z S and suppose tat T (z) = {T 1, T 2, T N } If z is a boundary vertex ten we enumerate te triangles suc tat T 1 and T N ave a boundary edge Moreover, we enumerate tem so tat T j, T j+1 sare an edge for j = 1, N 1 and T N and T 1 sare an edge in te case z is an interior vertex Let θ j denote te angle between te edges of T j originating from z We define { max{ θ 1 + θ 2 π,, θ N 1 + θ N π, θ N + θ 1 π } if z is an interior vertex Γ(z) = max{ θ 1 + θ 2 π,, θ N 1 + θ N π } if z is a boundary vertex Later, we will also need te following definition: { max{ sin(θ 1 + θ 2 ),, sin(θ N 1 + θ N ), sin(θ N + θ 1 ) } if z is an interior vertex Θ(z) = max{ sin(θ 1 + θ 2 ),, sin(θ N 1 + θ N ) } if z is a boundary vertex Definition 21 A vertex z S is a singular vertex if Γ(z) = 0 It is non-singular if Γ(z) > 0

4 4 J GUZMÁN AND R SCOTT We denote all te non-singular vertices by S 1 = {x S : x is non-singular }, and all singular vertices by S 2 = S \S 1 Let q be a function suc tat q T C(T ) 1 for all T T For eac vertex z S 2 define N A z (q) = ( 1) N j q Tj (z) Now we are ready to define te Scott-Vogelius finite element spaces for k 1: j=1 V k ={v [C 0(Ω)] 2 : v T [P k (T )] 2 T T }, Q k 1 ={q L 2 0(Ω) : q T P k 1 (T ) T T, A z (q) = 0 z S2 } Here P k (T ) is te space of polynomials of degree less tan or equal to k defined on T Also, L 2 0 (Ω) denotes te subspace of L2 of functions tat ave average zero on Ω Note tat if all te vertices are non-singular, ten Q k 1 is te space of discontinuous piecewise polynomials of degree k 1 wit average zero on Ω If singular vertices exist ten te pressure space is constrained on tose vertices Finally, te only way an interior vertex z is singular is if z as four edges coming out of it and tey lie on two lines Te definition of Q k 1 is natural as te following result wic was proved in [17] sows Lemma 1 For k 1, it olds div V k Qk 1 Tis lemma is a consequence of te following result Lemma 2 Let v be a vector field suc tat v T C 1 (T ) for every T T In addition, assume v C(Ω) and v vanises on Ω Ten, A z (div v) = 0 for every singular vertex z S2 We prove tis result for completeness in te appendix following te argument given in [17] Te goal of tis article is to prove te inf-sup stability of te pair V k, Qk 1 for k 4 We recall te definition of inf-sup stability Definition 22 Te pair of spaces V k, Qk 1 are inf-sup stable on a family of triangulations {T } if tere exists β > 0 suc tat for all Ω q div v dx (22) β q L 2 (Ω) sup q Q k 1 v V k,v 0 v H 1 (Ω) Let z S and e E (z) were e = {z, y} Define t z e = (y z) e tangent to e It is clear tat (23) t z e ψ y = 1 e on e to be te unit vector tat is 1 We define C k (T ) to be te subset of C k (T ) consisting of functions wit continuous limits on T, k 0

5 INF-SUP STABILITY OF UNFITTED STOKES ELEMENTS 5 Figure 2 Illustration of geometric quantities in one triangle We will also need to compute te derivatives of ψ y in all directions Suppose tat T T (y) and let g be te edge of T tat is opposite to y If we let n y T be te unit normal vector to g tat points out of T ten (24) ψ y T = 1 y n y T, T were y T is te distance of y to te line defined by te edge g If z is anoter vertex of T and we denote te edge e = {z, y} ten a simple calculation gives (25) y T = sin(θ) e were θ is te angle between te edges of T originating from z See Figure 2 for an illustration

6 6 J GUZMÁN AND R SCOTT 3 Establising te inf-sup condition 31 Preliminary stability results In order to prove te inf-sup stability of te Scott-Vogelius finite element spaces we need two well known results Te first follows from te stability of te Bernardi-Raugel [2] finite element space Proposition 1 Let k 1 Tere exists a constant α 1 suc tat for every p Q k 1 tere exists a v V 2 suc tat div v dx = p dx for all T T, and T T v H 1 (Ω) α 1 p L 2 (Ω) Te constant α 1 is independent of p and only depends on te sape regularity of te mes and Ω Te second result we need is contained in Lemma 25 of [19] A tree dimensional analog is given in [14] We provide te proof ere for completeness Proposition 2 Let k 1 Tere exists a constant α 2 suc tat for every T T and any p T P k 1 (T ) wit p vanising on te vertices of T and T p T dx = 0 tere exists a v T [P k (T )] 2 wit v T = 0 on T suc tat div v T = p T on T, and v T H 1 (T ) α 2 p T L 2 (T ) Here te constant α 2 depends only on te sape regularity of te mes and k Proof Consider te spaces M k 1 (T ) ={q P k 1 (T ) : q vanises at te vertices of T and and B k (T ) ={v [P k (T )] 2 : v = 0 on T } = {b T v : v [P k 3 (T )] 2 } Here b T is te cubic bubble of T tat vanises on T It is easy to sow tat te { 1 dim M k 1 (T )= 2k(k + 1) 4 for k 3, 0 for k 2, and dim B k (T )= { (k 2)(k 1) for k 3, 0 for k 2 T q(x) dx = 0} Moreover, it is clear tat div B k (T ) M k 1 (T ) Finally, we set Z k (T ) = {v B k (T ) : div v 0 on T } ten a simple argument gives Hence, Z k (T ) = { curl (ψb 2 T ) : ψ P k 5 (T )} dim Z k (T )= { 1 2 (k 4)(k 3) for k 5, 0 for k 4

7 INF-SUP STABILITY OF UNFITTED STOKES ELEMENTS 7 We claim tat M k 1 (T ) = div B k (T ) To sow tis, we will use a dimension count First note tat dim [B k (T )] = dim [div B k (T )] + dim [Z k (T )] For k 4 we know tat Z k (T ) is empty so tat dim [div B k (T )] = dim [B k (T )] = dim [M k 1 (T )] For k 4 dim [div B k (T )] = dim [B k (T )] dim [Z k (T )] = (k 2)(k 1) 1 (k 4)(k 3) 2 = 1 2 k2 3k (7/2)k 6 = 1 2 k k 4 = dim[m k 1 (T )] Tus we ave proved tat for any p T M k 1 (T ) tere exists a v T suc tat div v T = p T on T Te bound v T H 1 (T ) α 2 p T L 2 (T ) follows from a scaling argument mapping to te reference element wit te Piola transform wic preserves te divergence Summing up te result for every T T we can prove te following lemma Lemma 3 Let k 1 Let α 2 be te constant from te previous proposition For every p Q k 1 suc tat p(z) = 0 for all z S and T p dx = 0 for all T T tere exists v V k suc tat and div v = p on Ω, v H 1 (Ω) α 2 p L 2 (Ω) 32 Interpolating vertex values: Fundamental Vector Fields We will need to define vector fields tat will elp in interpolating pressure vertex values To do tis, we first need to define te following functions For every z S and e E (z) wit e = {z, y} we define te two functions η z e =ψ 2 zψ y γ z e =η z e 5 2 ψ2 zψ 2 y, were ψ z is defined in (21) Let T 1 and T 2 be te two triangles tat ave e as an edge Ten we can easily verify te following: (31a) (31b) (31c) support η z e T 1 T 2, support γ z e T 1 T 2, (ηe)(σ) z = 0 = (γe z )(σ) for σ S and σ z, γe z ds = 0 e

8 8 J GUZMÁN AND R SCOTT For example, to prove te last equation we used tat (32) ψzψ 2 y ds = e e 12, (33) ψzψ 2 y 2 ds = e e 30, wic can be done by transforming te edge e to te unit interval Te first vector field we define is (34) w z e = e t z eη z e We are going to need to consider w z e and vector fields of te form cγ z e Te following lemmas collect properties of tese functions Lemma 4 Let z S and e E (z) wit e = {z, y} and denote te two triangles tat ave e as an edge as T 1 and T 2 Let v = c e γ z e were c is a constant vector and let w z e be given by (34) It olds, (35a) (35b) (35c) (35d) (35e) (35f) (35g) w z e V 3, v V 4 support we z T 1 T 2, support v T 1 T 2 div we z dx = 0 = div v dx for all K T, K K div v(σ) = 0 = div w z e(σ) for σ S and σ z, div v Ts (z) = e c ψ y Ts = e y T s c n y T s for s = 1, 2 div w z e Ts (z) = 1 for s = 1, 2 w z e L 2 (T 1 T 2 ) C z, and v L 2 (T 1 T 2 ) C z c Te constant C only depends on te sape regularity Proof Te results (35a) and (35b) follow directly from te definitions of ηe z and γe z To prove (35c) note tat it trivially olds for K T if K is not T 1 or T 2 Terefore, let K = T s (for s = 1, 2) ten, using (35b) and integration by parts we get div v dx = v n ds K were n is te unit vector normal to e pointing out of K Terefore, by (31c) div v dx = e (c n) γe z ds = 0 Since t z e n = 0 we can also sow Te equations (35d) follow from (31b) K K e e div w z e dx = 0

9 INF-SUP STABILITY OF UNFITTED STOKES ELEMENTS 9 To prove (35e), we use tat ψ y (z) = 0 and ψ z (z) = 1 to get div v Ts (z) = e ψ 2 z(z) ( 1 5ψ y (z) ) c ψ y Ts + e 2ψ z (z) ( ψ y (z) 5 2 ψ y(z) 2) c ψ z Ts = e c ψ y Ts Te result follows from (24) Similarly, div w z e Ts (z) = e t z e ψ y Ts = e t z e ψ y e = 1, were we used tat ψ y is continuous along e and terefore t z e ψ y T1 = t z e ψ y T2 Finally, we used (23) To prove (35g), we use te Caucy-Scwarz inequality, and an inverse estimate v L 2 (T 1 T 2 ) C z v L (T 1 T 2 ) C v L (T 1 T 2 ) C c e γ z e L (T 1 T 2 ) C z c Here we used te sape regularity of te mes Te bound for w z e is similar We define te oter fundamental vector field in te following lemma Lemma 5 For every z S 1 and T T (z) tere exists a v z T V 4 wit te following properties (36a) (36b) (36c) (36d) div v z T (σ) = 0 for all σ z, div v z T T (z) = 1, and div v z T K (z) = 0 for all K T (z) and K T support vt z Ω (z), div vt z dx = 0 for all K T K Te following bound olds (37) v z T L 2 (Ω (z)) C z ( 1 Θ(z) + 1 ) Te constant C is independent of and only depends on te sape regularity and k Proof We adopt te notation given in te preliminary section wic we recall ere For z S we enumerate te triangles tat ave z as a vertex: T (z) = {T 1, T 2, T N } If z is a boundary vertex ten we enumerate te triangles suc tat T 1 and T N eac ave a boundary edge Moreover, we enumerate tem so tat T j, T j+1 sare an edge e j, for j = 1, N 1 and T N and T 1 sare an edge e N in te case z is an interior vertex Let θ j denote te angle between te edges of T j originating from z We let 1 s N 1 be suc tat sin(θ s + θ s+1 ) = Θ(z) Witout loss of generality we can assume tat s N; tis is immediate for a boundary vertex, and for an interior vertex we can enumerate te triangles accordingly Let e s = {z, y}, ten recall tat n y T s and n y T s+1 are te unit normal vectors pointing out of T s, T s+1, respectively, at te edges opposite to y Let t s+1 be te tangent vector to e s+1 pointing away from z ortogonal n y T s+1 (ie t s+1 n y T s+1 = 0) See Figure 3 for an illustration We need to define vt z j for 1 j N We start by defining vt z s = e sin(θ s) t s+1 ηe z Θ(z) s

10 10 J GUZMÁN AND R SCOTT Figure 3 Example, on non-singular vertex z, N = 5, s = 1 Ten, for every 1 l N s we define Also, for 1 l s 1 v z T s+l = ( 1) l 1( v z T s w z e s + w z e s ( 1) l 1 w z e s+l 1 ) v z T s l = ( 1) l 1( v z T s w z e s 1 + w z e s ( 1) l 1 w z e s (l 1) ) Wit tese definitions, (36a), (36c) and (36d) clearly follow from (35d),(35b) and (35c), respectively We are left to prove (36b) and (37) Let us first prove tis for T = T s By te definition of η z e s it as support in T s T s+1 and, terefore, we only need to consider K = T s and K = T s+1 First, by (35e) (38) div v z T s Ts+1 (z) = sin(θ s) Θ(z) On te oter and, (39) div v z T s Ts (z) = sin(θ s) Θ(z) e y t s+1 n y T s+1 = 0 T s+1 e y t s+1 n y T s = t s+1 n y T s = 1, T s Θ(z) were we used tat t s+1 n y T s = cos(θ s + θ s+1 + π 2 ) = sin(θ s + θ s+1 ) = Θ(z) Te inequality (37) follows from (35g) to get (310) v z T s L 2 (Ω (z)) C e s Θ(z) C z Θ(z) For T = T s+l, 1 l N s and T = T s l for 1 l s 1 we can use (38), (39) and (35f) to prove (36b) Te bound (37) follows from (310) and (35g), using N is bounded depending only on te sape regularity Using tese vector fields we can prove a crucial lemma

11 INF-SUP STABILITY OF UNFITTED STOKES ELEMENTS 11 Lemma 6 For every p Q k 1 and z S tere exists a v suc tat te following properties old: (311a) (311b) (311c) (311d) div v(σ) = 0 for all σ z, div v T (z) = p T (z) for all T T (z), support v Ω (z), div v dx = 0 for all K T K If z is a singular vertex ten v V 3 wit te bound (312) v L 2 (Ω (z)) C p L 2 (Ω (z)) If z is a non-singular vertex ten v V 4 (313) v L 2 (Ω (z)) wit te bound C Θ(z) + 1 p L 2 (Ω (z)) Te constant C is independent of p and and only depends on te sape regularity and k Proof We adopt te notation from te proof of te previous lemma We set a j = p Tj (z) First, suppose tat z is a singular vertex Ten, we know by te definition of Q k 1 N (314) ( 1) N j a j = 0 We define were we set j=1 v = N 1 j=1 b j w z e j ( j ) b j = ( 1) j l a l for all 1 j N 1 l=1 We immediately see tat (311a), (311c) and (311d) from (35d), (35b), (35c), respectively Moreover, (35a) gives tat v V 3 Using (35f) we ave For j = 1 we ave Also, div v Tj (z) = b j 1 + b j = a j for 2 j N 1 div v T1 (z) = b 1 = a 1 N 1 div v TN (z) = b N 1 = ( 1) N 1 l a l = a N, were we used (314) Terefore, we ave sown (311b) We see from (35g) tat ( N 1 ) 1/2 N 1 N N v L 2 (Ω (z)) C z C z b j C z N a j C z N p L (T j ), j=1 b 2 j j=1 l=1 j=1 j=1

12 12 J GUZMÁN AND R SCOTT were te equivalence of norms uses tat N is bounded depending only on te sape regularity Tis bound on N also implies inequality (312) after applying te inverse estimate p L (T j ) C Tj p L 2 (T j ) Next, we assume z is a non-singular singular vertex In tis case, we define v = N a j vt z j j=1 Clearly, (311a), (311b), (311c), (311d) follow from (36a), (36b), (36c), (36d) Using (37) we get ( 1 ) v L 2 (Ω (z)) C z Θ(z) + 1 p L (Ω (z)) Te inequality (313) follows after applying inverse estimates We can use te previous lemma to prove a global result First, we define Lemma 7 For every p Q k 1 (315a) (315b) Θ min = min Θ(z) z S 1 tere exists v V 4 (div v p)(z) = 0 for all z S, div v dx = 0 for all K T and ( 1 ) (316) v L 2 (Ω) C Θ min K p L 2 (Ω), were te constant C 1 is independent of p and depends only on te sape regularity of te meses and k Proof Let p Q k 1 be given Given z S let v z denote te vector field satisfying te properties of te previous lemma Ten, we set v = z S v z Clearly, from te previous lemma, (315a) olds Finally, since only tree v z s are non-zero on eac given triangle T we can easily sow tat v 2 L 2 (Ω) = v 2 L 2 (T ) = v z 2 L 2 (T ) = v z 2 L 2 (Ω (z)) T T T T z T z S C p 2 L 2 (Ω (z)) + C p 2 L 2 (Ω (z)) Θmin 2 z S 1 z S \S 1 Te inequality (316) now easily follows

13 INF-SUP STABILITY OF UNFITTED STOKES ELEMENTS Te Final Step We can now combine all te above results to prove te inf-sup condition Teorem 1 Suppose tat our family of meses {T } is non-degenerate (sape regular) Ten, Q k 1, V k satisfy te inf-sup condition (22) for k 4 were te constant β depends on Θ min and k, but is independent of Proof Let p Q k 1 Let v 1 V 2 be from Proposition 1 and let p 1 = p div v 1 We ave tat T p 1 dx = 0 for all T T By Lemma 1 we ave tat p 1 Q k 1 Given p 1 let v 2 V 4 be te corresponding vector field from Lemma 7 Ten, p 2 = p 1 div v 2 satisfies T p 2 dx = 0 for all T T and p 2 vanises at all te vertices We can, terefore, apply Lemma 3 and ave a v 3 V k so tat div v 3 = p 2 on Ω Setting v = v 1 + v 2 + v 3 V k we ave and v L 2 (Ω) (α 2 p 2 L 2 (Ω) + C 1 ( 1 Θ min + 1 div v = p on Ω, ) ( 1 ) p 1 L 2 (Ω) + α 1 p L 2 (Ω)) C Θ min p L 2 (Ω), were C 2 depends only on sape regularity, k and Ω Terefore, using Poincare s inequality p 2 L 2 (Ω) = Ω p div w dx p div v dx v H 1 (Ω) sup Ω w V k,w 0 w H 1 (Ω) Hence, ( 1 ) p L 2 (Ω) C Θ min sup w V k,w 0 p div w dx Ω w H 1 (Ω) Te result now follows by letting β = 1/ ( C 2 ( 1 Θ min + 1 )) 4 Appendix: Proof of Lemma 2 Let z S 2 We first assume tat y is an interior vertex In tis case, Ω (y) = 4 i=1 T i; see Figure 4 Our ypotesis tells us tat v Ti C 1 (T i ) for eac i and tat v C(Ω (y)) Let F (ˆx) = M ˆx + z were M = [t 1, t 2 ] is a 2 2 matrix Define ˆv(ˆx) = M 1 v(f (ˆx)) Note tat tis is te Piola transformation Let R 1,, R 4 be te four standard quadrants, and let L 1,, L 4 be te four semi-finite lines; see Figure 5 If we let B be a ball wit center at te origin and wit small enoug radius ten we know tat ˆv C(B) and ˆv B Ri C 1 (B R i ) for eac i A standard calculation sows tat div v(f (ˆx)) = div ˆv(ˆx) for ˆx B Terefore, 4 A y (div v) = i=1 ( 1) 4 i div v Ti (z) = 4 ( 1) 4 i div ˆv Ri (0) i=1

14 14 J GUZMÁN AND R SCOTT Figure 4 Illustration of Ω (z) If we let ˆv i = ˆv Ti ten we see tat 4 ( 1) 4 i div ˆv Ri (0) = div ˆv 1 (0) + div ˆv 2 (0) div ˆv 3 (0) + div ˆv 4 (0) i=1 = ( 1ˆv ˆv 2)(0) 1 + ( 1ˆv 1(0) 2 + 2ˆv 2(0)) 2 ( 1ˆv ˆv 2)(0) 3 + ( 1ˆv ˆv 2)(0) 2 = 1 (ˆv 1 1 ˆv 1)(0) 4 1 (ˆv 1 3 ˆv 1)(0) 2 2 (ˆv 2 1 ˆv 2)(0) 2 2 (ˆv 2 3 ˆv 2)(0) 4 =0 In te last step we used tat since ˆv is continuous ˆv 1 1 ˆv4 1 is identically zero on L 1 B Similarly, ˆv 2 1 ˆv2 2 vanises on L 2 B and so on Tis proves tat A y (div v) = 0 If, y S2 is a boundary vertex ten we extend v by zero to Ω c and apply te previous result References [1] D N Arnold and J Qin Quadratic velocity/linear pressure stokes elements Advances in computer metods for partial differential equations, 7:28 34, 1992

15 INF-SUP STABILITY OF UNFITTED STOKES ELEMENTS 15 Figure 5 Illustration of four quadrants [2] C Bernardi and G Raugel Analysis of some finite elements for te stokes problem Matematics of Computation, pages 71 79, 1985 [3] S Brenner and R Scott Te matematical teory of finite element metods, volume 15 Springer Science & Business Media, 2007 [4] B Cockburn, G Kanscat, and D Scötzau A locally conservative ldg metod for te incompressible navier-stokes equations Matematics of Computation, 74(251): , 2005 [5] J A Evans and T J Huges Isogeometric divergence-conforming b-splines for te unsteady navier stokes equations Journal of Computational Pysics, 241: , 2013 [6] R S Falk and M Neilan Stokes complexes and te construction of stable finite elements wit pointwise mass conservation SIAM Journal on Numerical Analysis, 51(2): , 2013 [7] J Guzmán and M Neilan A family of nonconforming elements for te brinkman problem IMA Journal of Numerical Analysis, page drr040, 2012 [8] J Guzmán and M Neilan Conforming and divergence-free stokes elements in tree dimensions IMA Journal of Numerical Analysis, 34(4): , 2014 [9] J Guzmán and M Neilan Conforming and divergence-free stokes elements on general triangular meses Matematics of Computation, 83(285):15 36, 2014 [10] S Harald Cristiansen and K Hu Generalized Finite Element Systems for smoot differential forms and Stokes problem ArXiv e-prints, May 2016 [11] V Jon, A Linke, C Merdon, M Neilan, and L G Rebolz On te divergence constraint in mixed finite element metods for incompressible flows SIAM Review, 2016 [12] A Linke Collision in a cross-saped domain a steady 2d navier stokes example demonstrating te importance of mass conservation in cfd Computer Metods in Applied Mecanics and Engineering, 198(41): , 2009 [13] C Lovadina and G Vacca Divergence free virtual elements for te stokes problem on polygonal meses ESAIM: Matematical Modelling and Numerical Analysis, 2015 [14] M Neilan Discrete and conforming smoot de ram complexes in tree dimensions Matematics of Computation, 84(295): , 2015

16 16 J GUZMÁN AND R SCOTT [15] J Qin On te convergence of some low order mixed finite elements for incompressible fluids PD tesis, Penn State, 1994 [16] L Scott and M Vogelius Norm estimates for a maximal rigt inverse of te divergence operator in spaces of piecewise polynomials RAIRO-Modélisation matématique et analyse numérique, 19(1): , 1985 [17] L R Scott and M Vogelius Conforming finite element metods for incompressible and nearly incompressible continua Tecnical report, DTIC Document, 1984 [18] X-C Tai and R Winter A discrete de ram complex wit enanced smootness Calcolo, 43(4): , 2006 [19] M Vogelius A rigt-inverse for te divergence operator in spaces of piecewise polynomials Numerisce Matematik, 41(1):19 37, 1983 [20] S Zang On te family of divergence-free finite elements on tetraedral grids for te stokes equations Preprint University of Delaware, 2007 [21] S Zang Divergence-free finite elements on tetraedral grids for k 6 Matematics of Computation, 80(274): , 2011 Division of Applied Matematics, Brown University, Providence, RI 02912, USA address: jonny guzman@brownedu Departments of Computer Science and Matematics, Committee on Computational and Applied Matematics, University of Cicago, Cicago IL 60637, USA address: ridg@ucicagoedu