THE INF-SUP STABILITY OF THE LOWEST ORDER TAYLOR-HOOD PAIR ON ANISOTROPIC MESHES arxiv: v1 [math.na] 21 Oct 2017

Transcription

1 THE INF-SUP STABILITY OF THE LOWEST ORDER TAYLOR-HOOD PAIR ON ANISOTROPIC MESHES arxiv: v1 [mat.na] 21 Oct 2017 GABRIEL R. BARRENECHEA AND ANDREAS WACHTEL Abstract. Uniform LBB conditions are desirable to approximate te solution of Navier Stokes, Oseen, and Stokes equations on anisotropic meses and to enable anisotropic refinements. We prove suc conditions for te second order Taylor Hood pairs Q 2 Q 1 and P 2 P 1 on a class of anisotropic meses. Tese meses may contain refined edge and corner patces. To tis end, we generalise Verfürt s trick and recent results by some of te autors. 1. Introduction Te finite element metod for te Stokes problem is subject to te satisfaction of te discrete inf-sup condition. For an effective metod, te discrete velocity and pressure spaces sould be balanced correctly. Tis balance results in a discrete inf-sup constant tat is independent of mes properties, suc as te size and sape of te elements. Concerning te first requirement, tere are many finite element pairs wic ave been proved to be inf-sup stable on regular meses (see [BBF13] for an extensive recent review). Now, concerning te second requirement, for te vast majority of metods te discrete inf-sup constant as only been proved to be independent of te size of te elements, but it may depend on te aspect ratio of te elements of te partition. Tis work addresses te last point raised in te previous paragrap. Te stability of finite element pairs as been less studied in te anisotropic case, but some progress as been made for some specific pairs, especially using discontinuous pressures and, possibly, non-conforming elements for te velocity. For example, in [ANS01] te autors analyse te Crouzeix-Raviart element in anisotropic meses (see also [AD99]). In [SS98, SSS99] te autors consider quadrilateral and triangular elements for te p-fem, and analyse pairs like te Q 2 k+1 Q k 1 (k 1) wit continuous velocities and discontinuous pressures. Teir analysis sows tat tis family is uniformly inf-sup stable in edge patces. Later, in [AC00] te quadrilateral element Q 2 k+1 P k 1 was analysed and it was proved tat its inf-sup constant depends on te geometric properties of te quadrilaterals. As a matter of fact, in tat work it was sown tat, wile tis pair is uniformly inf-sup stable on edge patces, its stability constant degenerates wit te square root of te aspect ratio in corner patces. A solution to tis undesired beaviour as been recently proposed in [ABW15]. For a review and furter references on tis topic, see [AMR03]. Date: Version 9, July 22, Matematics Subject Classification. 65N12, 65N30, 65N50. Key words and prases. anisotropic mes, Taylor, Hood, LBB condition. Tis work as been partially supported by te Leverulme Trust under grant RPG AW wants to acknowledge te support given by te Asociación Mexicana de Cultura A.C.. 1

2 2 GABRIEL R. BARRENECHEA AND ANDREAS WACHTEL In tis work we study te stability of te lowest-order Taylor-Hood pair in anisotropic meses. Tis pair was originally proposed in [TH73], and was first analysed in [BP79] and [Ver84] for te lowest order on triangles. In te independent works [BF91] and [Ste90] te analysis was extended to iger-order families, bot in te triangular and quadrilateral cases, in two space dimensions. Later, in [Bof97], te tree-dimensional case was addressed, using tetraedral meses. As far as we are aware of, no general proof of stability is available for te Taylor- Hood pair in te anisotropic case. Te only exception is, up to our best knowledge, te work [BLR12], were te inf-sup condition is proved under te assumption tat no drastic cange of sizes of neigbouring cells in any direction may occur, see [BLR12, Assumption 1]. Moreover, for some anisotropic edge patces, negative results, in te form of numerical experiments, are given in [SSS99] and [AMR03]. In particular, in [AMR03] instabilities are reported for te P 2 2 P 1 pair wen te aspect ratio tends to zero in certain configurations. Te purpose of tis work is ten to give a rigorous proof of te inf-sup stability of te lowest order Taylor-Hood pair, bot in te quadrilateral and triangular cases. Our main result states tat if te partitions ave enoug internal degrees of freedom, and tese are located in te rigt locations, ten tis element is uniformly inf-sup stable. Tis is, te inf-sup constant is independent of te aspect ratio. Our proof will be valid bot for edge and corner patces, tus improving upon te results from [BLR12]. Finally, we will sow, by means of numerical experiments, tat our assumptions are optimal, wic complements te results from [AMR03]. Te rest of te paper is organised as follows. In 2 we present te problem of interest and give some notation. Our main results are ten stated in 3 and tested numerically in 4. Te proofs for te quadrilateral case are ten presented in 5, and for te triangular case in Preliminaries and notation LetΩ R 2 beanopen,bounded, connectedandpolygonaldomain. Trougout, we use standard notation for Sobolev spaces (see [GR86]), namely, for D Ω, L 2 (D) (resp., L 2 0 (D)) stands for te space of (generalised) functions wic are square integrable in D (resp., wic belong to L 2 (D) and ave zero mean value in D), H 1 (D) (H0 1(D)) are elements of L2 (D) wose first order derivatives belong to L 2 (D) (and wose trace is zero on D). Vector-valued spaces and functions will be denoted using bold-faced letters. Te inner product in L 2 (D) (or L 2 (D)) is denoted by (, ) D, wit associated norm 0,D, te norm (seminorm) in H 1 (D) is denoted by 1,D ( 1,D ). By virtue of te Poincaré inequality 1,D is a norm on H0(D). 1 Finally, we denote by χ D te caracteristic function of D. A classical result (see [GR86]) is te following inf-sup condition: Tere exists β Ω > 0, depending only on Ω, suc tat (1) inf q M sup v V (divv,q) Ω v 1,Ω q 0,Ω β Ω > 0, were we ave used V M := H 1 0 (Ω) L2 0 (Ω). Te purpose of tis work is to prove te discrete analogue of tis result for te lowest order Taylor Hood element on anisotropic meses.

3 TAYLOR HOOD PAIRS IN ANISOTROPIC MESHES Finite element spaces. Te objective of te present work is to prove a discrete analogue of (1) for te lowest Taylor-Hood pair on anisotropic partitions. Te main assumptions on a partition P of Ω are te following: P is conforming, tat is, te non-empty intersection of two distinct cells is eiter a common point or a common edge e in te set of edges E P ; Pcanbesplitintoteunionofmacroelementsω Ω. Teseindividualsets ω are regular in sape, and te cells in P forming ω will be sape-regular, an edge patc (as in Figure 1(d-f)), or a corner patc (as in Figure 2(c-d) or as described below in Remark 4). In Section 6, we will analyse as well te Taylor-Hood space defined in simplicial triangulations. For tat case we will suppose tat P is a triangulation obtained by dividing eac quadrilateral of a partition satisfying te assumptions above into two triangles of equal area. We now define te discontinuous space Q l,p := { q L 2 (Ω): q K Q l (K) for K P } for l = 1,2, and for ω Ω te (locally) continuous spaces Q l,p (ω) := {q Q l,p : suppq ω} C 0 (ω). Let M P (ω) := Q 1,P (ω) L 2 0(ω); in te case of ω = Ω we write sortly M P. Let V P (ω) := [Q 2,P (ω)] 2 H 1 0(ω); in te case of ω = Ω we write sortly V P. Trougout tis manuscript, C (wit or witout subscript) will denote a positive constant, wic will be independent of te size, and aspect ratio of te elements of a given partition. Te value of suc a constant needs not to be equal wenever written in two different places. 3. Main results In tis section we state te main results of tis work. Te proofs of tese results are postponed to 5. To ease te readability, we ave restricted te presentation of te results to rectangular partitions, but te proofs can be extended, witout major difficulty, to parallelogram meses. It is wort mentioning tat, under tis assumption, te finite element spaces defined in te last section are mapped finite element spaces. Furtermore, we sketc in 6 ow te proofs can be extended to triangulated edge and corner patces. Below, Teorem 1 states te main result for edge patces, wile Teorem 3 covers te case of corner patces. Te results are presented for only one patc in eac case, but te proof of a global inf-sup condition (i.e., an inf-sup condition involving V P and M P ) can be obtained using a macro-element tecnique. Hypotesis 1. We will suppose tat every edge patc contains at least one division in te direction ortogonal to te long and tin elements. More precisely, every long and tin element is divided into two long and tin elements. Teorem 1. Let ω Ω be partitioned as an edge patc as depicted in Figure 1 (d,e,f). Ten, tere exists a constant β > 0, independent of te sape and size of te elements in P, suc tat, (2) inf sup q M P (ω) v V P (ω) (divv,q) ω v 1,ω q 0,ω β.

4 4 GABRIEL R. BARRENECHEA AND ANDREAS WACHTEL (a) (b) (c) (d) (e) (f) Figure 1. Edge patces. Te cases in te top row do not satisfy Hypotesis 1, but te cases in te bottom row do satisfy it. Remark 2. One clear consequence of te proof of Teorem 1 is tat te inf-sup constant β is te same one for te cases depicted in Figure 1(d) and Figure 1(e). Tis opens te door to study te case in wic an edge patc as been refined furter. More precisely, if we say tat te edge patc in Figure 1(d) is refined r = 0 times, and te patc in Figure 1(e) is refined r = 1 times, ten, continuing in te same way, i.e., bisecting te long and tin elements r times, we can define a patc tat as been refined r 2 times. Our conjecture, wic is supported by numerical evidence sown in Section 4, is tat te inf-sup constant is not affected by te value of r, but a formal proof of tis fact is lacking. Te next result concerns corner patces. A corner patc will be decomposed as ω = ω c ω E were ω c is te sape-regular small region (saded in Figure 2 for eac corner patc) and ω E is an overlapped edge patc as discussed in Teorem 1 (see Figure 1). For te proof of Teorem 3, we will need te following assumption. Hypotesis 2. Te Partition on ω c ω is suc tat te pair V P (ω c ) M P (ω c ) is (uniformly) inf-sup stable wit a constant β c. Teorem 3. Let ω Ω be partitioned into a corner patc, tat is, ω = ω c ω E wit ω c and ω E satisfying Hypotesis 2 and 1, respectively (e.g. as sown in Figure 2 (c-d)). Ten, tere exists β P > 0, depending only on β, from (2), and β c,from Hypotesis 2, suc tat inf sup q M P (ω) v V P (ω) (divv,q) ω v 1,ω q 0,ω β P. Remark 4. A comment on Hypotesis 2 is in place. In te case of a single corner patc, Hypotesis 2 requires ω c to be refined (as in Figure 2(c)). Now, if a corner patc is formed joining several patces of te type sown in Fig.2(a), see for instance Fig.2(d), ten te presence of enoug degrees of freedom inside of ω c ensures tat Hypotesis 2 is satisfied witout te need of furter refinement.

5 TAYLOR HOOD PAIRS IN ANISOTROPIC MESHES 5 (a) (b) (c) (d) Figure 2. Corner patces. Case(a) does not satisfy Hypotesis 2, te oters do. However, te assumptions of Teorem 3 also exclude case (b) since te overlapped edge patc does not satisfy Hypotesis Numerical confirmation In tis section we report te discrete inf-sup constants for te Taylor-Hood pairs in different configurations of corner and edge patces. Our aim is to confirm te validity of te results from last section. Ten, we introduce te parameter in suc a way tat te sort side of an anisotropic element is of size, and te long one is O(1)(see Figure 3 for details). For te refinement on te direction ortogonal to te long and tin element required by Hypotesis 1, we ave introduced an extra edge at te midpoint of te long and tin edges (for example, in Figure 1(e) te extra vertical edge is located at x = 1/2, and in Figure 1(f) is located at x = (1+)/2) Figure 3. Parametrised geometries for te coarse partitions of te edge and corner patces sown in Figures 1 and Results on edge patces. Our first experiments (in Table 1) aim at sowing tat Hypotesis 1 imposes sufficient (and necessary) requirements on edge patces suc tat te pair Q 2 2 Q 1 is uniformly stable on tem. Te results reported in te

6 6 GABRIEL R. BARRENECHEA AND ANDREAS WACHTEL last two columns oftable 1 sow tat, once te requirements from Hypotesis 1 are fulfilled, te inf-sup constant remains bounded below by a constant independent of. Now, to assess te necessity of tis restriction, we also report, in te first two columns of Table 1 te results for partitions tat do not satisfy tis ypotesis, were we see tat te inf-sup constants degenerate wit. Table 1. Discrete inf-sup constants of te pair Q 2 2 Q 1 on edge patces sown in Figure 1 Fig. 1(a) Fig. 1(c) Fig. 1(d) Fig. 1(f) Results on corner patces. In Table 2 we report te results obtained for different configurations of corner patces. We confirm te results of Teorem 3 in te sense tat wenever te ypoteses imposed on te partitions are satisfied, te inf-sup constant remains bounded below by a constant independent of (as it can be seen in te last two columns). On te contrary, te first column of Table 2 sows tat if te ypoteses are violated, ten te inf-sup constant decays wit. Table 2. Discrete inf-sup constants of te pair Q 2 2 Q 1 on corner patces sown in Figure 2 Fig. 2(b) Fig. 2(c) Fig. 2(d) Te stability on an L-saped domain. Te previous examples were restricted to one single edge, or corner, patc. Next, we sow an example ( in wic edge and corner patces are combined. For tis, we coose Ω = int [ 2,1] [0,1] ) [0,1] [ 1,0], and partition it into an edge patc and a corner patc, as depicted in Figure 4. Once again te inf-sup constant stays bounded below by a constant independent of te value of Refined patces. In tis section we verify te claims made in Remark 2 numerically. For tis, we consider te refined edge and corner patces depicted in Figure 5. In Table 3 we report te values of te inf-sup constants. We observe tat, as was conjectured in Remark 2, for te case of a refined edge patc, te inf-sup constants reported in te first two columns remain independent of te number of

7 TAYLOR HOOD PAIRS IN ANISOTROPIC MESHES 7 1 β P Figure 4. Discrete inf-sup constants for te Taylor-Hood element in te L-sape domain. refinements, tus confirming tat conjecture, at least numerically. Now, for te case of a refined corner patc, we see tat, as mentioned in Remark 11 (after te proof of Teorem 3), te value of te inf-sup constants (reported on te tird column of Table 3) decreases wit te number of refinements, tus confirming te sarpness of te proof of Teorem 3. Figure 5. Tese patces are refinements of level r = 2 of edgepatces and a corner patc. Table 3. We fix te parameter = Te columns below contain te LBB constants on (r-times) refined patces r Fig. 5(left) Fig. 5(centre) Fig. 5(rigt) Proofs on rectangular meses We start wit a preliminary result. Tis sligt generalisation of[jo16, Teorem 3.89] will be te main tool we will base our proof of stability on.

8 8 GABRIEL R. BARRENECHEA AND ANDREAS WACHTEL Lemma 5. Let us suppose tere exists a space B L 2 (Ω), suc tat (3) sup v V P (q,divv) Ω v 1,Ω β 1 q 0,Ω for all q B, and a projection Π B : M P B suc tat (4) sup v V P (q,divv) Ω v 1,Ω β 2 q Π B q 0,Ω for all q M P, were β 1,β 2 are positive constants. Ten V P M P is inf-sup stable wit inf-sup constant β 0 := β1β2 1+β 1+β 2, tis is (q,divv) sup Ω β 0 q v V P v 0,Ω for all q M P. 1,Ω Proof. Let q M P. Ten, using (3) we get (q,divv) sup Ω v V P v 1,Ω = sup v V P { (Π B q,divv) Ω v 1,Ω β 1 Π B q 0,Ω q Π B q 0,Ω. + (q Π Bq,divv) Ω v 1,Ω Multiplying (4) by (β 1 +1)/β 2 and adding it to te last inequality yields ( 1+ β 1 +1 ) (q,divv) sup Ω ( ) β 1 ΠB q β 0,Ω + q Π B q 0,Ω β1 q 0,Ω, 2 v V P v 1,Ω } as required Proof of Teorem 1. Our aim is to prove te assumptions of Lemma 5. For simplicity of te presentation we will consider a local coordinate system, and will consider ω = [ H,H] [, + 2H] to be partitioned into an edge patc as sown in Figure 6. We will consider ω = M M, were te bottom M = [ H,H] [,] is divided into 2 anisotropic elements, as depicted in Figure 7 (left), or 4 anisotropic elements, as in Figure 7 (rigt), respectively. On te oter and, te top M = [ H,H] [,+2H] is divided into 2 sape-regular elements. Figure 6. Edge patces in te local coordinate system.

9 TAYLOR HOOD PAIRS IN ANISOTROPIC MESHES 9 (5) For furter use we define te following linearly independent functions φ 0 := χ M M M χ M, φ 1,M (x,y) := x ( H χ M (x,y) and φ 2,M(x,y) := 1 2 x H ) χ M (x,y). Te pressure space M P (ω) is included in M P := M P (M ) M P (M) span{φ 0 }. We will prove te inf-sup condition between tis larger pressure space and V P (ω). As a consequence, te Taylor Hood pair will be uniformly stable on ω. As a first step, in te next result we state a decomposition of te space M P (M). Figure 7. Flat parts of edge patces. Lemma 6. Let M = [ H,H] [,]be partitioned as in Figure 7. Ten, M P (M) can be decomposed as M P (M) = B M G M, were and B M := {q M P (M): y q = 0 in M} = span{φ 1,M, φ 2,M }, G M := {q M P (M): (q,q s ) M = 0 for all q s B M }. Moreover, for all v H 1 0(M) te following olds (6) ( y v,q) M = 0 for all q B M, and, for every q G M tere exists (0,w) V P (M) H 1 0 (M) suc tat (7) ( y w,q) M = q 2 0,M and w 1,M C 1 q 0,M, were C 1 is independent of te size and sape of te elements in P. Proof. As it is clear from te context, we will omit te subscript M in tis proof. Te ortogonality (6) follows using integration by parts. Hence, we only need to prove (7). We split te remainder of te proof in two cases. Case 1 : M = [ H,H] [,] is split at x = 0 into two anisotropic cells belonging to P, see Figure 7(left). First, we extend te set {φ 1,φ 2 } to a basis of M P (M) = span{φ 1,φ 2,φ 3,φ 4,φ 5 } were φ 1,φ 2 are defined as in (5), and φ 3 (x,y) := y, φ 4 := 3φ 1 φ 3 and φ 5 := 3φ 2 φ 3. A direct computation, using Fubini s Teorem and te fact tat eac of te functions {φ 1,φ 2,φ 3 } as zero average, yields te ortogonalities (φ i,φ j ) M = δ ij φ i 2 0,M = δ M ij for i,j {1,...,5}, 3 and ten G M = span{φ 3,φ 4,φ 5 }, as φ 1,φ 2 B M. To prove (7), let g 1 G M, tat is, g 1 = 5 q iφ i, forsomerealcoefficients q 3,q 4,q 5, and g 1 2 0,M = 1 3 M 5 q2 i. Let p 0,p 1,p 2 M be te locations of te degrees of freedom (dof) of V P (ω) in M (see Figure 7, left), and let b 0,b 1,b 2 H0 1(M) be te piecewise Q 2 bubble functions, suc tat b i (p j ) = δ ij and (0,b i ) V P (ω). Let, in addition, v 3 := b 2 +b 0,

10 10 GABRIEL R. BARRENECHEA AND ANDREAS WACHTEL v 4 := b 2 b 0 and v 5 := b (b 0+b 2 ). A direct calculation using Fubini s teorem gives (8) ( y v 3,φ 3 ) M = 16H 9, ( y v 3,φ 4 ) M ( y v 3,φ 5 ) M ( y v 4,φ 3 ) M ( y v 4,φ 4 ) M = 8H 3 3, ( y v 4,φ 5 ) M ( y v 5,φ 3 ) M ( y v 5,φ 4 ) M ( y v 5,φ 5 ) M = 4H 3 3. Hence, defining v = α 3 q 3 v 3 + α 4 q 4 v 4 + α 5 q 5 v 5 wit α 3 = 3 4,α 4 = 3 2 and α 5 = 3, we obtain 5 (9) ( y v,g 1 ) M = q i φ i 2 0,M = g 1 2 0,M. Finally, since α i 3 and v i 2 1,M C 2 M /3, we get 5 (10) v 2 1,M 3 α i q i v i 2 1,M 3 C 5 qi 2 2( 2 M ) = 9 C 5 q i φ i 2 0,M, wic, since (0,v ) V P (M), finises te first case wit C 1 = 3 C. Case 2 : We now consider te case were M := [ H,H] [,] is divided into four rectangles K P obtained by splitting M along te lines x = 0 and y as depicted in Figure 7 (rigt). We extend te set {φ 1,φ 2,φ 3,φ 4,φ 5 }, introduced in te proof of te first case, to a basis of M P (M) by adding te following functions { φ 3 (x,2y ), if y > 0, φ 6 (x,y) := φ 3 (x,2y +), if y 0., φ 7 := 3φ 1 φ 6, φ 8 := 3φ 2 φ 6. Tese functions are even in y, and ten, proceeding as before we get (φ i,φ j ) M = δ ij φ i 2 0,M = δ M ij for i,j {1,...,8}. 3 Let now g = g 1 +g 2 G M, were g 1 = 5 q iφ i, and g 2 = 8 i=6 q iφ i, for some real coefficients q 3,...,q 8. Tanks to te definition of te function φ 6, we observe tat g 2 can be expressed as follows g 2 (x,y) = { (q 6 φ 3 +q 7 φ 4 +q 8 φ 5 )(x,2y ), if y > 0, (q 6 φ 3 +q 7 φ 4 +q 8 φ 5 )(x,2y +), if y 0. Ten, using te function v from te previous case, we consider te function ṽ := v +v2, were { v2(x,y) (α 6 q 6 v 3 +α 7 q 7 v 4 +α 8 q 8 v 5 )(x,2y ), if y > 0, := (α 6 q 6 v 3 +α 7 q 7 v 4 +α 8 q 8 v 5 )(x,2y +), if y 0, were α i,i = 6,7,8 are geometric constants to be determined. We note tat (0,v2) V P (M). Now, we prove ( y (v +v2),g 1 +g 2 ) M = g 2 0,M. To tis end, we recall (9) and since g 2 and y v2 are even in y, and g 1 and y v are odd in y, we obtain te ortogonalities ( y v 2,g 1 ) M = 0 and ( y v,g 2 ) M = 0.

11 TAYLOR HOOD PAIRS IN ANISOTROPIC MESHES 11 In addition, integrating by substitution and applying (8) gives ( y v 2,g 2 ) M = 2 H 0 H y v 2g 2 dxdy = 2 wic on coosing α i := α i 3 /2 (i = 6,7,8) gives ( y (v +v 2), g) M = 8 qi 2 α i ( y v i 3,φ i 3 ) M, i=6 5 8 q i φ i 2 0,M + q i φ i 3 2 0,M = g 2 0,M. Finally, we prove v 1,M C 1 q 0,M. First, applying Caucy s inequality and considering te scaling (w.r.t. y) inside v2 as well as α i = α i 3 /2, we get 8 ) ( 8 ) v2 2 1,M (4 3 α i q i v i 3 2 1,M = 3 qi α 2 i 3 v i 3 2 1,M. i=6 Now, using α i v i 2 1,M C M = 3 C φ i 2 0,M v 2 2 1,M 3 C i=6 i=6 8 q i φ i 3 2 0,M = 9 C g 2 2 0,M, i=6 (as in Case 1) we get were C is te same constant as in (10). Finally, recalling tat v is even in y, and v 2 is odd in y, we get ( v, v 2) M = 0 and arrive at (11) v +v 2 2 1,M = v 2 1,M + v 2 2 1,M 9 C ( g 1 2 0,M + g 2 2 0,M ) = 9 C g 2 0,M, wic finises te proof in tis case, again wit te same constant C 1 = 3 C. Now we rewrite M P = M P (M ) M P (M) span{φ 0 } = G ω B ω were (12) G ω := M P (M ) G M and B ω := B M span{φ 0 }, wit G M and B M defined as in Lemma 6. Let (13) X := V P (M ) {(0,v) V P (M)} V P (ω). Ten, using Lemma 6 and te stability of te Taylor Hood pair on M (wic is due to M being partitioned in a sape regular way), we conclude te uniform inf-sup condition (divv,q) (14) inf sup ω β 2, q G ω v X q 0,ω v 1,ω were β 2 is independent of te quotient H 1. From (14) it follows tat G ω is controlled by a subset of V P (ω) tat vanises on γ := M M. To control te remaining part of M P we need te bubble functions connecting M and M. Tis is stated in te next result. Lemma 7. Let B ω be defined as in (12). Ten, tere exists a constant β 1 > 0 independent of and H, suc tat (divv,q) sup ω β 1 q v V P (ω) v 0,ω for all q B ω. 1,ω

12 12 GABRIEL R. BARRENECHEA AND ANDREAS WACHTEL Proof. Let n 0,n 1,n 2 γ be te locations of te degrees of freedom (dof) of V P (ω) on γ (see Figure 6), and let f 0,f 1,f 2 H0(ω) 1 be te unique piecewise Q 2 functions suc tat f i (n j ) = δ ij, (0,f i ) V P (ω), and f i are piecewise affine in y. Let q B ω be arbitrary, tat is, q = q 0 φ 0 +q 1 φ 1,M +q 2 φ 2,M for some q 0,q 1,q 2 R. Let w := 2 i=0 α iw i were α i are coefficients to be cosen and w 0 := f 0 + f 2, w 1 := f 2 f 0 and w 2 := f (f 0+f 2 ). Using te fact tat γ consists of two edges of equal lengts, a direct computation gives te ortogonalities (w 0,φ 1,M ) γ (w 0,φ 2,M ) γ Ten, integration by parts gives In addition (w 1, φ 0 ) γ (w 1,φ 2,M ) γ ( y w,q) ω = (w, q ) γ = α 0 q 0 (w 0, φ 0 ) γ + (w 2, φ 0 ) γ (w 2,φ 1,M ) γ = 0. 2 α i q i (w i,φ i,m ) γ. (w 1,φ 1,M ) γ = 2H 3, (w 2,φ 2,M ) γ = (f 1,φ 2,M ) γ = H 3, and (w 0, φ 0 ) γ = 4H 3 φ 0 γ. i=1 Since M φ 0 = φ 0 2 0,ω and φ i,m 2 0,ω = φ i 2 0,M = (2H/3)2, we set α 0 := 3q 0 α 1 = 2q 1 and α 2 = 4q 2 to obtain ( y w,q) ω = q 2 0,ω. Finally, since α i 4q i and w i 2 1,ω 3 2 j=0 f j 2 1,ω C( 2 M +H 2 M ), we get 2 2 ( w 2 1,ω 3 α i w i 2 1 1,ω C qi ) 2 M + H 2 M i=0 C i=0 2 ( qi M 2 1+ M ) M C i=0 wic finises te proof since (0,w ) V P (ω). 2 q i φ i 2 0,ω i=0 =: β 2 1 q 2 0,ω, Te proof of Teorem 1 appears ten as an application of Lemma 5. In fact, recall te decomposition M P (ω) M P = G ω B ω, were G ω,b ω are defined in (12). Let q G ω B ω, and recall te definition (13) of X. Let also Π B q be te L 2 (ω) projection of q onto B ω. Ten, from (6) and te fact tat Π B (q) is constant in M, we get (Π B q,divv) ω = 0 for all v X. Hence (q,divv) sup ω v V P v 1,ω (q,divv) sup ω v X v 1,ω (q Π B q,divv) = sup ω v X v 1,ω β 2 q Π B q 0,ω, were we applied q Π B q G ω and (14). Finally Lemma 7 gives (Π B q,divv) sup ω v V P v 1,ω β 1 Π B q 0,ω,

13 TAYLOR HOOD PAIRS IN ANISOTROPIC MESHES 13 wic finises te proof upon application of Lemma 5 wit a constant β = β1β2 1+β 1+β 2 independent of size and aspect ratios of te edge patc. Remark8. Te essential part of te independence of te constant β (in Teorem 1) of te aspect ratio is given by te proof of Lemma 6. We observed in tat proof tat, for bot cases, te constant C 1 is given by 3 C, and ten, it was unaffected by te additional refinement of te edge patc (see te inequalities (10) and (11)). Tis confirms wat was claimed in Remark 2, and is affirmed by numerical experiments. Remark 9. It is wort noticing tat te proof of Teorem 1 provides te inf-sup stability of a family of elements wic contains te lowest order Taylor Hood pair. In fact, te pressure space is allowed to be discontinuous across γ, wic someow generalises te results known so far. Remark 10. A closer look at te proofs of Lemmas 6, 7 and Teorem 1 sows tat te tecnique can be applied to overlapped edge patces, as sown for instance in Figure 1(f) or Figure 8(centre). In tis case, Lemma 6 as to be applied twice, tat is, tere are two flat parts, say M 1 and M 2, were M 2 contains te vertically aligned long and tin cells. Ten, M P (M 2 ) as to be decomposed into B M2 = {q M P (M 2 ): x q = 0 in M 2 } and G M2. Hence, using te same arguments tere exists (v,0) V P (M 2 ) suc tat ( x v,g) M2 = g 2 0,M 2 and v 1,M2 C 1 g 0,M2 for all g G M2. Ten, similar to (12) and te arguments tereafter, te pair consisting of te velocity space X := {(0,v) V P (M 1 )} V P (M ) {(v,0) V P (M 2 )} and te pressure space G ω := G M1 M P (M ) G M2 is uniformly inf-sup stable. Ten, an adapted proof of Lemma 7 sows tat B ω = B M1 B M2 span{φ 0,φ 0}, were φ 0 is anoter piecewise constant function, is be controlled by te six bubbles on te edges connecting M 1,M and M 2,M. Finally, te proof of Teorem 1 follows analogously considering tese modified definitions of G ω,b ω and X. Figure 8. A corner patc decomposed into ω c and ω E satisfying Hypotesis 2 and 1, respectively Proof of Teorem 3. In tis section we use notations ω c,ω E,ω as suggested in Figure 8 for te complete and parts of te (refined) corner patc. As done in te previous section we see tat M P (ω) M P := M P (ω E ) M P (ω c ) span{φ c }, were φ c is defined as { 1 if (x,y) ωc, φ c (x,y) = ωc ω E oterwise.

14 14 GABRIEL R. BARRENECHEA AND ANDREAS WACHTEL We apply te tecnique developed in [ABW15] to prove Teorem 3. Tanks to Teorem 1, we know tat te pair V P (ω E ) M P (ω E ) is uniformly inf-sup stable, tat is, its stability constant does not depend on H 1. In addition, if we recall te Hypotesis 2, ten te Taylor Hood pair is inf-sup stable on ω c. Let now q M P (ω), tat is, q = q +Π c q, were q M P (ω c ) M P (ω E ) and Π c q span{φ c } is constant on ω c and ω E. Hence, te inf-sup conditions on ω c and ω E (mentioned above) imply (q,divv) sup ω (q,divv) sup ω min{β,β c } q v V P (ω) v 1,ω v V P (ω E) V P (ω c) v 0,ω, 1,ω were β is from Teorem 1 and β c is from Hypotesis 2. We will ten finis te proof by sowing tere exists a constant C > 0 suc tat q 0,ω C Π c q 0,ω. To tis end, we note tat te projection Π c q is given by Π c q = q ωc χ ωc + q ωe χ ωe, and since 0 = q ω, we ave q ωe = ω c ω E 1 q ωc. Hence ( ) (15) Π c q 2 0,ω = ω E q 2 ω E + ω c q 2 ωc ω c = ω E +1 ω c q 2 ω c = C 3 ω c q 2 ω c, were C 3 = ω ω E 1 1. Now let Γ c := ω c ω E, and let e Γ c be an arbitrary edge satisfying e = K e K e wit K e ω c. Ten, from q ωc q ωe = Π c q e, and since q is continuous in ω, we get C 3 q ωc = Π c q e = q e and q 2 ω c = C 2 3 q 2 e. Ten, te inequality of Caucy and Scwarz and [ABW15, Lemma 2.1] yield q 2 ω c = 1 q 2 ω Γ c c Γ c 1 = C3 2 Γ q 2 e c e Γ e c 1 1 C3 2 Γ c e Γ e e q 2 0,e c 1 ( e 1 C3 2 Γ c K e + 1 ) K e Γ c e 4 q 2 0,K e K e. Now, using K e K e, as required. e ω c Γ c Π c q 2 0,ω = C 3 ω c q 2 ω c 1 e C 3 Γ c e Γ c K e = 1 and 1/C 3 1 we get from (15) ω c K e 8 q 2 0,K e K e 8 q 2 0,ω, Remark 11. In tis proof we ave not considered te possibility of refining te corner patces. If we were to consider tat case, ten in te last step of te above

15 TAYLOR HOOD PAIRS IN ANISOTROPIC MESHES 15 proof we would be led to use te fact tat, for a corner patc tat as been refined r times, we ave e ω c Γ c K e = 2r 1. Tis would give as a result te inf-sup constant β P C2 r 1 2 wic sows a dependency of te inf-sup constants on te refinement of te partition. Tis dependency is not an artefact of te proof, as te numerical results in Section 4.4 sow. 6. Te triangulated case In tis section we prove Teorem 1 and 3 for triangulated edge- and corner patces. A single cange is required to prove Teorem 3, i.e., we ave to replace te used discrete trace estimate by one tat is valid on anisotropic triangles, see e.g. [WH03, Teorem 3]. Ten it remains to prove Teorem 1. To tis end, we note tat te proofs of Lemma 5, Lemma 7 remain valid and we only need to prove Lemma 6 for triangulated anisotropic (flat) macro elements. Proof. We perform a direct computation wic is simplified by using te following tree-point quadrature formula f dx = K 3 f(m i ), 3 K i=1 were f P 2 (K) and m i are te mid-points of te edges of triangle K. Now, we let M = [ H,H] [,] and notice tat owever we triangulate M Te space B M (i.e., φ 1,φ 2 ) defined in Lemma 6 is part of M P (M). Next we prove (7) on te triangulated macro element. Case 1: M is split at x = 0 and ten triangulated. Again, we ave dimg M = 3 and for a basis of G M we coose {ϕ i } 5 to be te nodal interpolants of {φ i} 5 in Lemma 6. Tese functions take te values +1/ 1/0 at te signs +/ /0 sown in Figure 9. For simplicity (but witout loss of generality) we cose a triangulation on wic ϕ 3,ϕ 5 remain even in x, wile ϕ 4 is odd in x. Additionally te product ϕ 3 ϕ 5 vanises in eac ofte quadraturepoints, and terefore{ϕ i } i is an ortogonalbasis of G M. Hence, for g G M wit g = 5 q iϕ i, te quadrature formula yields 5 g 2 0,M = q i ϕ i 2 0,M = M ( 2q q4 2 +2q5 2 ). ϕ 3 ϕ 4 ϕ 5 Figure 9. Basis of G M wit sign-patterns of ϕ i and y ϕ i.

16 16 GABRIEL R. BARRENECHEA AND ANDREAS WACHTEL Now, let b 0,b 1,b 2 be te quadratic bubble functions, eac supported on two triangles and taking te value 1 at te midpoint of te sared edge. Ten defining v 3,v 4,v 5 similar to Lemma 6 by v 3 := b 2 +b 0, v 4 := b 2 b 0 and v 5 := 2b 1 (b 2 +b 0 ) we get ( y v i,g) M = (v i, y g) M = q i (v i, y ϕ i ) M for i = 3,4,5, and ten v := 5 α iv i satisfies ( y v,g) M = (v, y g) M = 5 α i q i (v i, y ϕ i ) M. To calculate tese products we use te quadrature formula and te values of ( y ϕ i ) K, given inside te triangles sown in Figure 9. We obtain (v 3, y ϕ 3 ) M = M /(3), (v 4, y ϕ 4 ) M = M /(6), (v 5, y ϕ 5 ) M = (2b 1, y ϕ 5 ) M = M /(3), and coosing α i = q i we obtain 5 ( y v,g) M = q i ϕ i 2 0,M = g 2 0,M. Finally, since α i q i and v i 2 1,M C 2 M we get 5 5 v 2 1,M C α i v i 2 1,M C 2 qi( 2 2 M ) = C g 2 0,M, wic finises case 1 as (0,v ) V P (M). Case 2: On appropriate triangulations of a 2-by-2 macro element one can construct ortogonal bases of te velocity and pressure spaces, again even and odd on y. Ten, te proof follows similar to te direct computation above A numerical confirmation. Te results presented in Table 4 sow tat te Taylor Hood pair is also uniformly stable on triangulated edge and corner patces. Furtermore, we see tat te difference between te inf-sup constants in column 1 and 2 is marginal, tat is, te additional refinement does not affect te inf-sup constant. 1/2 1/2 Figure 10. Triangulated edge and corner patces satisfying te ypoteses of Teorems 1 and 3. In all cases, Ω = (0,1) 2.

17 TAYLOR HOOD PAIRS IN ANISOTROPIC MESHES 17 Table 4. Discrete inf-sup constants of te P 2 2 P 1 pair on te partitions sown in Fig. 10 Fig. 10(left) Fig. 10(centre) Fig. 10(rigt) References [ABW15] M. Ainswort, G. R. Barrenecea, and A. Wactel, Stabilization of ig aspect ratio mixed finite elements for incompressible flow, SIAM Journal on Numerical Analysis 53 (2015), no. 2, [AC00] M. Ainswort and P. Coggins, Te stability of mixed p-finite element metods for Stokes flow on ig aspect ratio elements, SIAM J. Numer. Anal. 38 (2000), no. 5, (electronic). [AD99] G. Acosta and R. G. Durán, Te maximum angle condition for mixed and nonconforming elements: application to te Stokes equations, SIAM J. Numer. Anal. 37 (1999), no. 1, (electronic). [AMR03] T. Apel and H. Maaravo Randrianarivony, Stability of discretizations of te Stokes problem on anisotropic meses, Mat. Comput. Simulation 61 (2003), no. 3-6, [ANS01] T. Apel, S. Nicaise, and J. Scöberl, Crouzeix-Raviart type finite elements on anisotropic meses, Numer. Mat. 89 (2001), no. 2, [BBF13] D. Boffi, F. Brezzi, and M. Fortin, Mixed finite element metods and applications, Springer Series in Computational Matematics, vol. 44, Springer, Heidelberg, [BF91] Franco Brezzi and Ricard S. Falk, Stability of iger-order ood-taylor metods, SIAM Journal on Numerical Analysis 28 (1991), no. 3, [BLR12] M. Braack, G. Lube, and L. Röe, Divergence preserving interpolation on anisotropic quadrilateral meses, Comput. Metods Appl. Mat. 12 (2012), no. 2, [Bof97] D. Boffi, Tree-dimensional finite element metods for te stokes problem, SIAM Journal on Numerical Analysis 34 (1997), no. 2, [BP79] M. Bercovier and O. Pironneau, Error estimates for finite element metod solution of te stokes problem in te primitive variables, Numerisce Matematik 33 (1979), no. 2, [GR86] V. Girault and P.-A. Raviart, Finite element metods for Navier Stokes equations, Springer, Berlin, [Jo16] V. Jon, Finite element metods for incompressible flow problems, Springer, Berlin, [SS98] D. Scötzau and C. Scwab, Mixed p-fem on anisotropic meses, Mat. Models Metods Appl. Sci. 8 (1998), no. 5, [SSS99] D. Scötzau, C. Scwab, and R. Stenberg, Mixed p-fem on anisotropic meses. II. Hanging nodes and tensor products of boundary layer meses, Numer. Mat. 83 (1999), no. 4, [Ste90] Rolf Stenberg, Error analysis of some finite element metods for te Stokes problem, Matematics of Computation 54 (1990), no. 190, [TH73] C. Taylor and P. Hood, A numerical solution of te Navier-Stokes equations using te finite element tecnique, Computers & Fluids 1 (1973), no. 1, [Ver84] R. Verfürt, Error estimates for a mixed finite element approximation of te Stokes equations, RAIRO Anal. Numér. 18 (1984), no. 2, [WH03] T. Warburton and J. S. Hestaven, On te constants in p-finite element trace inverse inequalities, Comput. Metods Appl. Mec. Engrg. 192 (2003), no. 25,

18 18 GABRIEL R. BARRENECHEA AND ANDREAS WACHTEL Department of Matematics and Statistics, University of Stratclyde, 26 Ricmond Street, Glasgow G1 1XH, Scotland address: Departamento Académico de Matemáticas, ITAM, Río Hondo 1, CDMX, México address: