arxiv: v1 [math.na] 20 Jul 2009

Transcription

1 STABILITY OF LAGRANGE ELEMENTS FOR THE MIXED LAPLACIAN DOUGLAS N. ARNOLD AND MARIE E. ROGNES arxiv: v1 [mat.na] 20 Jul 2009 Abstract. Te stability properties of simple element coices for te mixed formulation of te Laplacian are investigated numerically. Te element coices studied use vector Lagrange elements, i.e., te space of continuous piecewise polynomial vector fields of degree at most r, for te vector variable, and te divergence of tis space, wic consists of discontinuous piecewise polynomials of one degree lower, for te scalar variable. For polynomial degrees r equal 2 or 3, tis pair of spaces was found to be stable for all mes families tested. In particular, it is stable on diagonal mes families, in contrast to its beaviour for te Stokes equations. For degree r equal 1, stability olds for some meses, but not for oters. Additionally, convergence was observed precisely for te metods tat were observed to be stable. However, it seems tat optimal order L 2 estimates for te vector variable, known to old for r > 3, do not old for lower degrees. 1. Introduction In tis note, we consider approximations of te mixed Laplace equations wit Diriclet boundary conditions: Given a source g, find te velocity u and te pressure p suc tat u grad p = 0, div u = g in Ω, p = 0 on Ω, for a domain Ω R 2 wit boundary Ω. Te equations offer te classical weak formulation: Find a square integrable vector field wit square integrable divergence u H(div, Ω) and a square integrable function p L 2 (Ω) suc tat (1.1) u v + q div u + p div v = g q Ω for all v H(div, Ω) and q L 2 (Ω). Te above formulation can be discretized using a pair of finite dimensional spaces V H(div, Ω), Q L 2 (Ω), yielding discrete approximations u V and p Q satisfying (1.1) for all v V and q Q. As is well-known, te spaces V and Q must satisfy certain stability, or compatibility, conditions for te discretization to be well-beaved [7]. More precisely, tere must exist positive constants α and β suc tat for any, (1.2a) (1.2b) 0 < α α = inf 0 < β < β div u Z = inf q Q Ω u, v, v Z u div v div div v, q. v V v div q Matematics Subject Classification. Primary: 65N30. Key words and prases. mixed finite elements, Lagrange finite elements, stability. Te work of te first autor was ported in part by NSF grant DMS Te work of te second autor was ported by a Center of Excellence grant from te Norwegian Researc Council to Centre of Matematics for Applications at te University of Oslo and by a Center of Excellence grant from te Norwegian Researc Council to Center for Biomedical Computing at Simula Researc Laboratory. 1

2 2 DOUGLAS N. ARNOLD AND MARIE E. ROGNES Here, div and 0 denote te norms on H(div, Ω) and L 2 (Ω),, is te L 2 (Ω) inner product and (1.3) Z = {v V div v, q = 0 q Q }. Te two conditions will be referred to as te Brezzi coercivity and te Brezzi inf- condition for te mixed Laplacian. Te classical conforming discretizations of (1.1) rely on te finite element families of Raviart and Tomas [18] or Brezzi, Douglas and Marini [8] for te space V H(div) in order to satisfy tese conditions. In tis note, we sall consider te Lagrange vector element spaces, tat is, continuous piecewise polynomial vector fields defined relative to a triangulation T, for te space V. Tis is motivated by te following reasons. First, tese spaces are fairly inexpensive, simple to implement and post-process and in frequent use for oter purposes. Second, suc pairs would allow continuous approximations of te velocity variable, or wen viewed in connection wit linear elasticity, lay te ground for continuous approximations of te stress tensor. Moreover, in te recent years, tere as been an interest in mixed finite element discretizations tat are bot stable for (1.1) and for te Stokes equations: (1.4) grad u : grad v + q div u + p div v = f v Ω for all v H 1 (Ω; V) suc tat Ω v = 0 and all q L2 (Ω). Te searc for conforming suc discretizations is complicated by te fact tat te existing, stable discretizations of (1.1) are suc tat V H 1 (Ω; V). On te oter and, te existing stable discretizations of (1.4) are typically unstable for (1.1) [14]. Te existence of stable discretizations V Q of (1.1) suc tat V H 1 (Ω; V) becomes a natural separate question. Unfortunately, tere are no known suc finite element discretizations tat are stable for any admissible triangulation family {T }. In tis note, we aim to numerically examine cases were a reduced stability property may be identified. In tis sense, te investigations ere are in te spirit of te work of Capelle and Bate [10] and Qin [17]. For a family of conforming discretizations {V Q } of (1.1) suc tat div V Q for eac, te condition (1.2a) is trivial. Te stability conditions tus reduce to te condition (1.2b), namely te question of bounded Brezzi inf- constant β div. On te oter and, recall tat for te Stokes formulation (1.4), te corresponding Brezzi coercivity condition is trivial by te Poincaré inequality. Hence, for any family of conforming discretizations, te stability conditions for Stokes reduce to tat of a uniform bound for te Brezzi inf- constant β 1. Here, (1.5) β 1 = inf q Q div v, q. v V v 1 q 0 wen V H 1 (Ω; V), Q L 2 (Ω) and 1 denotes te norm on H 1 (Ω). Furter, suc a bound immediately gives (1.2b) since β 1 βdiv by definition. Hence, if div V Q, stability for Stokes immediately gives stability for te mixed Laplacian. Te conditions div V Q and V H 1 (Ω; V) are clearly satisfied by te element pairs consisting of continuous piecewise polynomial vector fields of degree less tan or equal to r and discontinuous piecewise polynomials of degree r 1, for r = 1, 2,.... Tis family could be viewed as an attractive family of elements for bot te Stokes equations and te mixed Laplacian. However, te Brezzi inf constant(s) will not be bounded for all r. For r 4, Scott and Vogelius demonstrated tat tese finite element spaces will be stable for te Stokes equations on triangulations tat ave no nearly singular vertices, tat is, triangulations tat are not singular in te appropriate sense [19]. Te lower order cases, 1 r Ω

3 STABILITY OF LAGRANGE ELEMENTS FOR THE MIXED LAPLACIAN 3 3, were studied carefully by Qin, concluding tat te elements are not stable in general [17]. However, tey are stable for some specific families of triangulations, and can be stabilized by removal of spurious pressure modes on some oter classes of triangulations. (Te space N of spurious pressure modes is defined in (2.10) below.) In general, te stability of finite element spaces for te Stokes equations as been extensively investigated. In addition to te previous references, surveys are presented in [3, 9]. However, to our knowledge, a careful study of te lower order cases as not been conducted for te mixed Laplacian. As te stability for te mixed Laplacian is a weaker requirement wen div V Q, tere may be a greater class of triangulations for wic te elements form a stable discretization. In fact, tis is known to be true. One example is provided by te pairing of continuous piecewise linear vector fields V and te subspace of discontinuous piecewise constants Q suc tat Q = div V on crisscross triangulations of te unit square. Qin proved tat tere does not exist a β > 0 suc tat β 1 > β for any [17, Lemma 7.3.2]. On te oter and, Boffi et al. proved tat suc a bound does exists for β div [4]. We sall present numerical evidence suggesting tat tere is a range of triangulations for wic tis penomenon occurs. Te main results are summarized below. Spurious modes: For r = 2, 3 and for all triangulations tested, te dimension of te space N of spurious modes is equal to te number of interior singular vertices σ. However, for r = 1 and one of te triangulation families studied ( Flipped, wic is defined in Figure 1), dim N is strictly greater tan σ. Stability: For all triangulations we ave tested, te metod seems at least reduced stable (i.e., stable after removal of spurious modes, if any), for r = 2, 3. Tis is in contrast to te situation for te Stokes equations, were for some triangulations, suc as te diagonal triangulation, te metod is not reduced stable for r = 2, 3, wile for oter triangulations, it is. For r = 1, reduced stability olds for some triangulations, but fails for oters, including te diagonal triangulation. Convergence: We also studied convergence of te metod on diagonal triangulations. For suc meses, te metod was observed to be stable for r > 1, but unstable for r = 1. Teory predicts optimal convergence of p in L 2 (Ω) and u in H(div, Ω) for a stable metod and tis is in fact wat was observed. Suc optimal convergence olds for r = 2, 3, 4, but not in te apparently unstable case r = 1. In te case r 4, it is known tat u converges at one order iger in L 2 (Ω) tan in H(div, Ω). No suc increase of order was observed for r < 4. Te note is organized as follows. We introduce furter notation and summarize some key points of te teory of mixed finite element metods in Section 2. Furter, we derive some eigenvalue problems associated wit te stability conditions and give a caracterization of te Brezzi inf- constant for te mixed Laplacian β div in Section 3. Tese eigenvalue problems applied to te Stokes equations were also stated by Malkus [13] and (in part) by Qin [17] and provide a foundation for numerical investigations of te Brezzi stability conditions. Based on tese general results, Section 4 is devoted to te study of continuous piecewise polynomials in two dimensions for te velocity and discontinuous piecewise polynomials for te pressure. 2. Notation and preliminaries Te notion of reduced stability of families of mixed finite element spaces is a key point in tis note. In order to make tis notion precise, tis preliminary section

4 4 DOUGLAS N. ARNOLD AND MARIE E. ROGNES aims to introduce notation and summarize te stability notions for finite element discretizations of abstract saddle point problems. If V is an inner product space, we denote te dual space by V, te inner product on V by, V and te induced norm by V. Let Ω be an open and bounded domain in R d wit boundary Ω. We let H m (Ω), for m = 0, 1,..., denote te standard Sobolev spaces of square integrable functions wit m weak derivatives and denote teir norm by m. Accordingly, H 0 (Ω) = L 2 (Ω). Te space of polynomials of degree r on Ω is denoted P r (Ω). Te space of vectors in R d is denoted V and in general, X(Ω; V) denotes te space of vector fields on Ω for wic eac component is in X(Ω). For brevity owever, te space of vector fields in L 2 (Ω; V) wit square integrable divergence is written H(div, Ω) wit norm div and semi-norm div = div 0. Te subscripts and te reference to te domain Ω will be omitted wen considered erfluous. Let T denote an admissible simplicial tessellation of Ω, measuring te mes size of te tessellation. We sall frequently refer to spaces of piecewise polynomials defined relative to suc, and label te spaces of continuous, and discontinuous, piecewise polynomials of degree less tan or equal to r as follows. P c r = P c r(t ) = {p H 1 (Ω) p K P r (K) K T } r = 1, 2,..., P r = P r (T ) = {p L 2 (Ω) p K P r (K) K T } r = 0, 1,.... Te classical abstract saddle point problem reads as follows [7, 9]: for given Hilbert spaces V and Q and data (f, g) V Q, find (u, p) V Q satisfying (2.1) a(u, v) + b(v, p) + b(u, q) = f, v + g, q (v, q) V Q, were a and b are assumed to be continuous, bilinear forms on V V and V Q, respectively. We sall assume ere and trougout tat a is symmetric. Following [2], tere exists a unique solution (u, p) of (2.1), if and only if te continuous Babuška inf- constant (2.2) γ = inf 0 (u,p) 0 (v,q) a(u, v) + b(v, p) + b(u, q) (u, p) V Q (v, q) V Q is positive. By [7], tis olds if and only if te continuous Brezzi coercivity and Brezzi inf- constants are positive. Tese are defined as (2.3) (2.4) α = β = inf 0 u Z 0 v Z inf 0 q Q 0 v V a(u, v) u V v V, b(v, q) v V q Q, respectively, were Z = {v V b(v, q) = 0 q Q}. Given finite dimensional spaces V V and Q Q, defined relative to a tessellation T of Ω, te Galerkin discretization of (2.1) takes te form: Find (u, p ) V Q satisfying (2.5) a(u, v) + b(v, p ) + b(u, q) = f, v + g, q (v, q) V Q.

5 STABILITY OF LAGRANGE ELEMENTS FOR THE MIXED LAPLACIAN 5 On te discrete level, te Babuška inf-, Brezzi coercivity and Brezzi inf- constants are defined as (2.6) (2.7) (2.8) were γ := α := β := inf 0 (u,p) V Q inf 0 u Z inf 0 q Q a(u, v) + b(v, p) + b(u, q), 0 (v,q) V Q (u, p) V Q (v, q) V Q a(u, v), 0 v Z u V v V b(v, q), 0 v V v V q Q (2.9) Z = {v V b(v, q) = 0 q Q }. For given V Q, tere exists a unique solution of (2.5) if and only if α and β (or equivalently γ ) are positive. Furtermore, for a family of discretization spaces V Q parameterized by, if α and β are uniformly bounded from below, ten one obtains te quasi-optimal approximation estimate [7]: ( ) u u V + p p Q C inf u v V + inf p q Q, v V q Q wit C depending only on te bounds for α and β and bounds on te bilinear forms a and b. Te uniform boundedness condition motivates te notion of stability for pairs of finite element spaces. Definition 2.1 (Stable discretization). A family of finite element discretizations {V Q } is stable in V Q if te Brezzi coercivity and inf- constants α and β (or equivalently te Babuška inf- constant γ ) are bounded from below by a positive constant independent of. In accordance wit standard terminology, we say tat {V Q } satisfies te Brezzi coercivity or inf- conditions if α or β, respectively, are uniformly bounded from below. Tere are families of discretizations tat are not stable in te sense defined above, but ave a reduced stability property. More precisely, for a pair V Q consider te space of spurious modes N Q : (2.10) N = {q Q b(v, q) = 0 v V }. For a stable discretization, N contains only te zero element. Indeed, β = 0 if and only if N contains non-zero elements. If N is non-trivial, it is natural to consider te reduced space N, te ortogonal complement of N in Q, in place of Q. Tis motivates te definition of te reduced Brezzi inf- constant, relating to te stability of V N : (2.11) β = inf 0 q N b(v, q), 0 v V v V q Q and te following definition of reduced stable. By definition, β 0. Definition 2.2 (Reduced stable discretization). A family of discretizations {V Q } is reduced stable in V Q if te Brezzi coercivity constant α and te reduced Brezzi inf- constant β, defined by (2.7) and (2.11), are bounded from below by a positive constant independent of.

6 6 DOUGLAS N. ARNOLD AND MARIE E. ROGNES 3. Eigenvalue problems related to te Babuška-Brezzi constants For a given set of discrete spaces, te Babuška and Brezzi constants defined by (2.6) (2.8) can be computed by means of eigenvalue problems. Te form and properties of te eigenvalue problem associated wit te Brezzi inf- constant for te Stokes equations were discussed by Qin in [17]. Since also te Brezzi coercivity constant plays a role for te mixed Laplacian, we begin tis section by deriving ow te Brezzi coercivity constant can be computed by similar eigenvalue problems. Actually, in our application in Section 4, te Brezzi coercivity condition will be automatic, owever, we discuss it ere in te abstract case, for te sake of completeness. Tese eigenvalue problems were also stated, and carefully analysed from an algebraic view-point, by Malkus [13] in connection wit te displacementpressure formulation of te linear elasticity equations. We continue by observing tat te continuous Brezzi inf- constant can be naturally associated wit te smallest eigenvalue of te Laplacian itself Eigenvalue problems for te discrete Babuška-Brezzi constants. Let V V and Q Q be given finite dimensional spaces as before. It follows easily from te definition tat te Babuška inf- constant γ = λ min wen λ min is te smallest (in modulus) eigenvalue of te following generalized eigenvalue problem: Find λ R, 0 (u, p) V Q satisfying (3.1) a(u, v) + b(v, p) + b(u, q) = λ ( u, v V + p, q Q ) (v, q) V Q. Te following lemma identifies an eigenvalue problem associated wit te Brezzi inf- constant. Lemma 3.1 (Qin [17, Lemma ]). Let λ min be te smallest eigenvalue of te following generalized eigenvalue problem: Find λ R, 0 (u, p) V Q satisfying (3.2) u, v V + b(v, p) + b(u, q) = λ p, q Q, (v, q) V Q. Ten, λ 0 and for β defined by (2.8), β = λ min. It can also be sown tat te reduced Brezzi inf- constant β equals te squareroot of te smallest non-zero eigenvalue of (3.2) [17, Teorem 5.1.1]. Te Babuška and Brezzi inf- constants are tus easily computed, given bases for te spaces V and Q. As for (3.1), it is easily seen tat te Brezzi coercivity constant α = λ min were λ min is te smallest (in modulus) eigenvalue of te eigenvalue problem: Find λ R and 0 u Z suc tat (3.3) a(u, v) = λ u, v V v Z. However, a basis for Z is usually not readily available, tus indering te actual computation of te eigenvalues of (3.3). Instead, te above eigenvalue problem over Z can be extended to a generalized eigenvalue problem over V Q : Find λ R and 0 (u, p) V Q suc tat (3.4) a(u, v) + b(v, p) + b(u, q) = λ u, v V (v, q) V Q. Te following lemma establises te equivalence between (3.3) and (3.4). Lemma 3.2. If (λ, u) is an eigenpair of (3.3), tere exists a p Q suc tat (λ, (u, p)) is an eigenpair of (3.4). Conversely, if (λ, (u, p)) is an eigenpair of (3.4) and u 0, ten u Z and (λ, u) is an eigenpair of (3.3). For 0 p N and any scalar λ, (λ, (0, p)) is an eigenpair of (3.4), and tese are te only eigenpairs of (3.4) wit u = 0.

7 STABILITY OF LAGRANGE ELEMENTS FOR THE MIXED LAPLACIAN 7 Proof. Let (λ, u) be an eigenpair of (3.3). Define B : V Q suc tat B v, q Q = b(v, q) for all q Q. Since B : Z B (V ) is an isomorpism, p B (V ) Q is well-defined by p, q Q = λ u, B 1 q V a(u, B 1 q) q B (V ). Ten, for any v Z, p satisfies b(v, p) = B v, p Q = λ u, v V a(u, v). Furter, by definition b(v, p) = 0 for any v Z. Hence, by te assumption tat (λ, u) is an eigenpair of (3.3), (λ, (u, p)) satisfies (3.4). Te converse statement is obvious. Finally, letting u = 0 in (3.4), we see tat (λ, (0, p) satisfies (3.4) if and only if p N, but for any λ R. Note tat, as a consequence of te last observation in Lemma 3.2, if N is nontrivial, te generalized eigenvalue problem (3.4) is computationally not well-posed since any scalar λ is an eigenvalue. In te subsequent section, we sall numerically investigate te stability of families of finite element discretizations V Q suc tat div V Q for te mixed Laplacian, using te eigenvalue problem (3.2) in terms of standard bases for te spaces V and Q. Te eigenvalue problem (3.4) does not enter, but would, were we to investigate discretizations were div V Q A caracterization of te mixed Laplacian Brezzi inf- constant. We now turn from te general setting to consider te H(div) L 2 formulation of te mixed Laplacian (1.1). In Lemma 3.3 below, we sow tat te Brezzi inf- constant can be identified wit te smallest eigenvalue of te negative Laplacian. Consequently, if a discretization family {V Q } guarantees eigenvalue convergence for te mixed Laplace eigenvalue problem and is suc tat div V Q, te Brezzi inf- constant of te discretization will converge to te continuous Brezzi inf- constant. Lemma 3.3. Let V H(div, Ω) and Q L 2 (Ω) be suc tat tat div V Q. Consider te Brezzi inf- eigenvalue problem (3.2) applied to (1.1): (3.5) u, v div + div v, p + div u, q = λ p, q (v, q) V Q. Consider also te mixed Laplace eigenvalue problem: (3.6) û, v + div v, ˆp + div û, q = ˆλ ˆp, q (v, q) V Q. Ten, (λ, (u, p)) is an eigenpair of (3.5) if and only if (ˆλ, (û, ˆp)) is an eigenpair of (3.6) were ˆλ = λ(1 λ) 1, û = u and ˆp = (1 λ)p. Moreover, in tis case, 0 λ < 1 and ˆλ 0. Also, p 0, ˆp 0. Proof. Assume tat (λ, (u, p)) is an eigenpair of (3.5). First, note tat λ 1. Letting λ = 1, v = u and q = div u in (3.5), implies tat u = 0. Furter, v = 0 and q = p gives tat p = 0. Hence, λ = 1 is only associated wit te zero solution, wic by definition, cannot form an eigenpair. Next, note tat p 0, since oterwise implies tat u = 0, wic again is impossible. By te assumption div V Q, (3.7) div u, div v = λ p, div v v V. Taking v = u in (3.7) and letting v = u and q = (λ 1)p in (3.5) sow tat u 2 = λ(1 λ) p 2. So, 0 λ < 1. Te combination of (3.7) and (3.5) gives u, v + div u, q + (1 λ) div v, p = λ p, q. Finally, letting û = u, ˆp = (1 λ)p and ˆλ = λ(1 λ) 1, gives tat (ˆλ, (û, ˆp)) solves (3.6). Te converse olds by similar arguments.

8 8 DOUGLAS N. ARNOLD AND MARIE E. ROGNES Te equivalence demonstrated in te lemma above affords a simple caracterization of te Brezzi inf- constant for te mixed Laplacian. Te eigenvalue problem (3.5) wit V = H(div, Ω) and Q = L 2 (Ω) is te eigenvalue problem associated wit te continuous Brezzi inf- constant β div, cf. Lemma 3.1. Hence, β div is te square-root of te smallest eigenvalue of (3.5). On te oter and, te eigenvalue problem (3.6) is a mixed weak formulation of te standard eigenvalue problem for te negative Laplacian wit Diriclet boundary conditions, given in strong form below: (3.8) ˆp = ˆλ ˆp in Ω, ˆp = 0 on Ω. Tus, if µ is te smallest eigenvalue of (3.8), β div = µ(1 + µ) 1. Remark. An alternative eigenvalue problem arises from noting tat, under te assumption div V Q, (3.6) wit v = 0 implies tat ˆλ ˆp = div û. Hence, if (ˆλ, (û, ˆp)) solves (3.6), ten eiter ˆλ = 0, û = 0 and ˆp div V, or ˆλ > 0, û 0, and (ˆλ, û) will be an eigenpair of te problem: Tis eigenproblem was studied in [5]. div u, div v = ˆλ u, v v V. Now, consider a stable discretization family V Q of (1.1) suc tat div V Q, wit Brezzi inf- constants β div. Let µ denote te smallest eigenvalue approximation of (3.6) by V Q. As a consequence of te preceding considerations, if µ µ, ten β div β div. In oter words, if te discretization family is stable, satisfies div V Q, and gives eigenvalue convergence, ten te Brezzi inf- constant will converge to te continuous Brezzi inf- constant. Note owever, tat te discrete stability conditions are not sufficient to ensure te convergence of approximations to te eigenvalue problem (3.6) [1, 4]. Mixed finite element discretizations of (1.1) based on te Raviart-Tomas [18] and Brezzi-Douglas-Marini [8] families of H(div) conforming elements are known to give eigenvalue convergence, and ence β div β div. Some cases were β div seems to be uniformly bounded in, but β div β div are exemplified in te subsequent section. Finally, note tat if Ω is te unit square: Ω = (0, 1) 2, te smallest eigenvalue of (3.8) is 2π 2 and so 2π (3.9) β div 2 = 1 + 2π Lower order Lagrange elements for te mixed Laplacian From ere on, we restrict our attention to finite element discretizations of te mixed Laplacian (1.1) on a polygonal domain Ω R 2. Te primary aim is to examine te stability, or reduced stability, and convergence properties of Lagrange elements, tat is, continuous piecewise polynomials for te vector variable and discontinuous piecewise polynomials for te scalar variable: (4.1) V Q = P c r(t ; V) P r 1 (T ), for r = 1, 2,.... Altoug te Brezzi conditions are in general not satisfied for tese discretizations, stability or reduced stability may be identified on families of structured triangulations. Te pair (4.1) is clearly suc tat div V Q. Terefore, te stability of te discretization relies on a uniform bound for te Brezzi inf- constant only. Furter, a uniform lower bound on te Brezzi inf- constant for te Stokes equations induces te corresponding bound for te mixed Laplacian. Hence, te results on te reduced stability of tis element pair for te

9 STABILITY OF LAGRANGE ELEMENTS FOR THE MIXED LAPLACIAN 9 Stokes equations can be directly applied to te mixed Laplacian. In te following, new numerical evidence is presented and compared to te known results. Te stability of te Pr(V) c P r 1 family of elements, for bot te Stokes equations and te mixed Laplacian, depends on te polynomial degree r and te structure of te triangulation T. For triangulations tat ave interior singular vertices, te space of spurious modes N, defined by (2.10) applied to (1.1), will be nontrivial. Here, an interior vertex is labeled singular if te edges meeting at tat vertex fall on two straigt lines. Let x denote an interior singular vertex and let ω x be te star of x. For any r 1, tere exists a p N suc tat p is ported in ω x [16, 15]. Consequently, letting σ denote te number of interior singular vertices of a triangulation, dim N σ. Scott and Vogelius sowed, for r 4 tat if tere are no interior singular vertices, ten dim N = 0 and so β div β 1 > 0. Moreover, tey proved tat for a family of meses witout interior singular vertices, β 1 remains bounded above zero as long as te meses do not tend to singularity as 0. For te precise statement and more details, see [19] or [6, Section 12.6]. As we sall see below, for r < 4, te space of spurious modes may be non-trivial even wen tere are no singular vertices. Furter, for te Stokes equations, more restrictive conditions tan te above must be placed on te triangulations in order to obtain a uniform bound for te Stokes Brezzi inf- constant [17]. Te stability properties of tese lower order discretizations for te mixed Laplacian is te main question of interest in te following. Remark. We sall not consider te pairing Pr(V) c P s except for s = r 1. Tis coice is easily motivated. First, a dimension count sows tat te pairing of continuous piecewise polynomial vector fields wit discontinuous piecewise polynomials of te same or iger degree must ave a non-trivial space of spurious modes. Second, altoug te Brezzi inf- constant is uniformly bounded for te pairs Pr(V) c P r 2, r = 2, 3,..., te Brezzi coercivity constant for te mixed Laplacian is not uniformly bounded, and tus stability fails Stability. In te spirit of [17, Section 5], we aim to numerically investigate te stability of P c r(v) P r 1 for r = 1, 2, 3 on certain families of structured triangulations of te unit square. Te triangulation patterns considered are illustrated and labeled in Figure 1. For n even, an n n triangulation of eac family is constructed by first partitioning te domain into n n squares, and subsequently dividing eac block of 2 2 squares into triangles by te respective patterns. For instance, an n n diagonal triangulation is formed by dividing te unit square into n n subsquares, and dividing eac subsquare into triangles by te positive diagonal. Trougout, we identify = 1/n and assume tat n > 2. Observe tat te diagonal, flipped, and zigzag triangulations contain no interior singular vertices, wile te crisscross and te Union Jack triangulation contain n 2 and n(n 2)/2 interior singular vertices, respectively. Tis is summarized in te first row of Table 1. Recall tat dim N σ for r 1 and equality olds for r 4. Qin proved tat equality olds for 1 r 3 in te case of te diagonal and te crisscross meses and numerically observed equality for te flipped mes for r = 2 [17]. Our own experiments sow tat equality olds for te zigzag and Union Jack meses for 1 r 3. Equality also olds for te flipped mes wen r = 2, 3, but not for r = 1. Tese results are summarized in te second row of Table 1. We continue by studying te beaviour of te Brezzi inf- constants on te above triangulations. Te cases r = 2, 3 are considered first, but we will return to te case r = 1 below. For te Stokes equations, it is known tat te diagonal and crisscross triangulation families exibit very different beaviour for r = 2, 3 [17]. Namely, altoug tere are non-trivial spurious modes on te crisscross triangulation family, te reduced Brezzi inf- constant is uniformly bounded. In contrast,

10 10 DOUGLAS N. ARNOLD AND MARIE E. ROGNES (a) Diagonal (b) Flipped (c) Zigzag (d) Crisscross (e) Union Jack Figure 1. Structured 2 2 triangulations of te unit square. Diagonal Flipped Zigzag Crisscross Union Jack σ 0 { 0 0 n 2 n(n 2)/2 (n/2 1) 2 r = 1 dim N 0 0 n 2 n(n 2)/2 0 r = 2, 3 Table 1. Te number of interior singular vertices σ and te dimension of te space of spurious modes, dim N, for labeled families of n n triangulations of te unit square, cf. Figure 1. For te flipped, zigzag and Union Jack meses, dim N is conjectural. for te diagonal family, te Brezzi inf- constant decays as approximately O(). As te discretization is reduced stable for te Stokes equations on crisscross triangulations, it is also reduced stable for te mixed Laplacian. A natural question becomes weter te lack of stability on diagonal triangulations for te Stokes equations is also present for te mixed Laplacian. In view of Lemma 3.1, we sall make an attempt at answering tis question troug a set of numerical experiments. For a given r and a given T, te smallest, and smallest non-zero, eigenvalue of (3.5) for V = Pr(T c, V), Q = P r 1 (T ) give te Brezzi inf- and reduced Brezzi inf- constant. Tese eigenvalues for te triangulation families considered, computed using LAPACK, SLEPc [11] and DOLFIN [12], are given for r = 1, 2, 3 in Tables 2, 3 and 4. For te purpose of identifying spurious modes, eigenvalues below a tresold of 10 4 ave been tabulated to zero 1. For te diagonal meses, te numerical experiments indicate tat in contrast to Stokes, te mixed Laplacian Brezzi inf- constants are bounded from below for bot r = 2, 3. For te flipped and zigzag meses, experiments give similar 1 Had a smaller tresold been cosen, some of te zero eigenvalues associated to interior singular vertices would ave been missed for r = 2 on te Union Jack mes of size n = 6.

11 STABILITY OF LAGRANGE ELEMENTS FOR THE MIXED LAPLACIAN 11 β div β div (dim N ) n Diagonal Zigzag Flipped Union Jack (1) (4) (4) (12) (9) (24) (16) (40) (25) (60) (36) (84) (49) (112) Table 2. Te mixed Laplacian (reduced) Brezzi inf- constant for P c 1(T ; V) P 0 (T ) on labeled structured families of triangulations T. Te dimension of te space of spurious modes in parentesis if non-trivial. β div β div (dim N ) n Diagonal Zigzag Flipped Union Jack (4) (12) (24) (40) (60) Table 3. Te mixed Laplacian (reduced) Brezzi inf- constant for P c 2(T ; V) P 1 (T ) on labeled structured families of triangulations T. Te dimension of te space of spurious modes in parentesis if non-trivial. β div β div (dim N ) n Diagonal Zigzag Flipped Union Jack (4) (12) (24) Table 4. Te mixed Laplacian (reduced) Brezzi inf- constant for P c 3(T ; V) P 2 (T ) on labeled structured families of triangulations T. Te dimension of te space of spurious modes in parentesis if non-trivial.

12 12 DOUGLAS N. ARNOLD AND MARIE E. ROGNES results. Neiter exibits any spurious modes. Moreover, wile te Brezzi inf- constant decays approximately as O() for te Stokes equations [17], it appears to be uniformly bounded for te mixed Laplacian. For te Union Jack family, te same experiment gives n(n 2)/2 spurious modes, but te reduced Brezzi inf- constant again seems to be uniformly bounded. In summary for r = 2, 3, te Pr(V) c P r 1 elements appear to be at least reduced stable for all te families considered. Wit te discussion in Section 3.2 in mind, we also note tat te Brezzi inf- constant converges to te exact value, given by (3.9), for some, but not all, of tese meses. For r = 2, te Brezzi inf- constant seems to converge to te exact value on te diagonal meses, but not on te flipped or te zigzag meses. Te situation is te opposite for r = 3. Tere, te Brezzi inf- constant seems to converge to te exact value on te zigzag and flipped meses, but not for te diagonal meses. Te situation is different and more diverse in te lowest-order case: r = 1. Boffi et al. proved tat P1(V) c P 0 is in fact reduced stable for te mixed Laplacian on crisscross meses [4]. It is not reduced stable for Stokes [17]. However, te element pair does not seem to be stable on diagonal meses. Te values in te first column of Table 2 indicate tat te Brezzi inf- constant decays approximately as O(). Te same is te case for te zigzag meses. For te Union Jack meses, te situation is similar to te crisscross case. Tat is, te number of singular modes matc te number of interior singular vertices and te reduced Brezzi inf- constant appears to be bounded from below. Finally, te flipped meses display a surprising beaviour. Tere seem to be (n/2 1) 2 spurious modes, even toug tere are no singular vertices. Tis is te only case were we ave observed dim N > σ. However, te reduced Brezzi inf- constant appears to be uniformly bounded Convergence. In te previous subsections, we ave investigated te stability of te P c r(v) P r 1 elements. Now, we proceed to examine te convergence properties of tese elements on te diagonal meses. Conjecturing tat P c r(v) P r 1 is stable on tis mes family for r 2, in accordance wit te numerical evidence presented above, te standard teory gives te error estimate (4.2) u u div + p p 0 C r ( u r+1 + p r ). For r 4, te L 2 error estimate for te velocity can be improved [6, Teorem ], tus yielding: (4.3) u u 0 C r+1 u r+1. In order to verify (4.2) and to see weter (4.3) appears to be attained for r = 2, 3, we consider a standard smoot exact solution to te Laplacian wit pure Diriclet boundary conditions: (4.4) p(x, y) = sin(2πx) sin(2πy), u = grad p, g = div u. Te errors of te P c r(v) P r 1 approximations for 1 r 4 on diagonal meses can be examined in Figure 2. To compute te errors, bot te source function g and te exact solutions u, p ave been represented by sixt order piecewise polynomial interpolants Π g and Π u, Π p, wereupon te errors ave been calculated exactly (up to numerical precision). For r = 1, we observed te discretization to be unstable on diagonal meses. As expected in tis case, neiter te pressure nor te velocity approximation seems to converge in te L 2 norm. Tis indicates tat te estimates for te approximation error, based on te standard estimates and te decaying Brezzi inf- constant, cannot be improved. On te oter and, for r > 1, we observed te metod to be stable. For r = 2, 3, 4, te orders of convergence in te H(div) norm of te velocity

13 STABILITY OF LAGRANGE ELEMENTS FOR THE MIXED LAPLACIAN 13 (a) Normalized pressure errors in 0. (b) Normalized velocity errors in div (c) Normalized velocity errors in 0 Figure 2. Te errors of P c r(v) P r 1 approximations for r = 1, 2, 3, 4 on diagonal meses versus mes number n. Te errors ave been normalized, tat is, multiplied by te inverse of te error at te smallest mes n = 4.

14 14 DOUGLAS N. ARNOLD AND MARIE E. ROGNES and te L 2 norm of te pressure approximations are indeed optimal, as predicted by (4.2). Te situation seems different for te convergence of te velocity approximation in te L 2 norm. For r 4, a convergence rate of order r + 1 is predicted by (4.3). Tis is also observed for r = 4 in Figure 2(c). On te oter and, for r = 2 and 3 te rate of convergence appears to be of order r, and tus one order suboptimal, indicating tat te estimate (4.3) does not old for r < 4. References [1] D. N. Arnold, R. S. Falk, and R. Winter. Finite element exterior calculus, omological tecniques, and applications. Acta Numer., 15:1 155, [2] I. Babuška. Te finite element metod wit Lagrangian multipliers. Numer. Mat., 20: , 1972/73. [3] D. Boffi, F. Brezzi, and M. Fortin. Finite elements for te Stokes problem. In Mixed Finite Elements, Compatibility Conditions, and Applications, volume 1939 of Lecture Notes in Matematics, pages Springer-Verlag, Berlin, Lectures given at te C.I.M.E. Summer Scool eld in Cetraro, June 26 July 1, 2006, Edited by D. Boffi and L. Gastaldi. [4] D. Boffi, F. Brezzi, and L. Gastaldi. On te problem of spurious eigenvalues in te approximation of linear elliptic problems in mixed form. Mat. Comp., 69(229): , [5] D. Boffi, R. G. Duran, and L. Gastaldi. A remark on spurious eigenvalues in a square. Appl. Mat. Lett., 12(3): , [6] S. C. Brenner and L. R. Scott. Te matematical teory of finite element metods, volume 15 of Texts in Applied Matematics. Springer, New York, tird edition, [7] F. Brezzi. On te existence, uniqueness and approximation of saddle-point problems arising from Lagrangian multipliers. Rev. Française Automat. Informat. Recerce Opérationnelle Sér. Rouge, 8(R-2): , [8] F. Brezzi, J. Douglas, Jr., and L. D. Marini. Two families of mixed finite elements for second order elliptic problems. Numer. Mat., 47(2): , [9] F. Brezzi and M. Fortin. Mixed and ybrid finite element metods, volume 15 of Springer Series in Computational Matematics. Springer-Verlag, New York, [10] D. Capelle and K.-J. Bate. Te inf- test. Comput. & Structures, 47(4-5): , [11] V. Hernandez, J. E. Roman, and V. Vidal. SLEPc: a scalable and flexible toolkit for te solution of eigenvalue problems. ACM Trans. Mat. Software, 31(3): , [12] A. Logg, G. N. Wells, et al. DOLFIN. URL: ttp// [13] D. S. Malkus. Eigenproblems associated wit te discrete LBB condition for incompressible finite elements. Internat. J. Engrg. Sci., 19(10): , [14] K. A. Mardal, X.-C. Tai, and R. Winter. A robust finite element metod for Darcy-Stokes flow. SIAM J. Numer. Anal., 40(5): (electronic), [15] J. Morgan and L. R. Scott. Te dimension of te space of C 1 piecewise polynomials. Researc Report UH/MD 78, University of Houston, Matematics Department, [16] J. Morgan and R. Scott. A nodal basis for C 1 piecewise polynomials of degree n 5. Mat. Comput., 29: , [17] J. Qin. On te Convergence of Some Low Order Mixed Finite Elements for Incompressible Fluids. PD tesis, Penn State, [18] P.-A. Raviart and J. M. Tomas. A mixed finite element metod for 2nd order elliptic problems. In Matematical aspects of finite element metods (Proc. Conf., Consiglio Naz. delle Ricerce (C.N.R.), Rome, 1975), pages Lecture Notes in Mat., Vol Springer, Berlin, [19] L. R. Scott and M. Vogelius. Norm estimates for a maximal rigt inverse of te divergence operator in spaces of piecewise polynomials. RAIRO Modél. Mat. Anal. Numér., 19(1): , Scool of Matematics, 512 Vincent Hall, 206 Curc Street S.E., University of Minnesota, Minneapolis, MN address: arnold@umn.edu Center for Biomedical Computing, Simula Researc Laboratory, P.O.Box 134, 1325 Lysaker, Norway. address: meg@simula.no