arxiv: v1 [math.na] 17 Jul 2014

Transcription

1 Div First-Order System LL* FOSLL* for Second-Order Elliptic Partial Differential Equations Ziqiang Cai Rob Falgout Sun Zang arxiv: v1 [mat.na] 17 Jul 2014 February 13, 2018 Abstract. Te first-order system LL* FOSLL* approac for general second-order elliptic partial differential equations was proposed and analyzed in [10], in order to retain te full efficiency of te L 2 norm first-order system least-squares FOSLS approac wile exibiting te generality of te inverse-norm FOSLS approac. Te FOSLL* approac in [10] was applied to te div-curl system wit added slack variables, and ence it is quite complicated. In tis paper, we apply te FOSLL* approac to te div system and establis its well-posedness. For te corresponding finite element approximation, we obtain a quasi-optimal a priori error bound under te same regularity assumption as te standard Galerkin metod, but witout te restriction to sufficiently small mes size. Unlike te FOSLS approac, te FOSLL* approac does not ave a free a posteriori error estimator, we ten propose an explicit residual error estimator and establis its reliability and efficiency bounds. Key words. LL* metod, least-squares metod, a priori error estimate, a posteriori error estimate, elliptic equations. AMSMOS subject classifications. 65M60, 65M15 1 Introduction Tere are substantial interests in te use of least squares principles for te approximate solution of partial differential equations wit applications in bot solid and fluid mecanics. Many least-squares metods for te scalar elliptic partial differential equations ave been proposed and analyzed, [3, 15]. Teir numerical properties depend on coices suc as te first-order system and te least-squares norm. Loosely speaking, tere are tree types of least-squares metods: te inverse approac, te div approac, and te div-curl approac. Te inverse approac employs an inverse norm tat is furter replaced by eiter te weigted mes-dependent norm see [2] or te discrete H 1 norm see [7] for computational feasibility. Bot te div and te div-curl approaces use te L 2 norm and Tis work performed under te auspices of te U.S. Department of Energy by Lawrence Livermore National Laboratory under Contract DE-AC52-07NA27344 LLNL-JRNL Tis work was supported in part by te National Science Foundation under grant DMS Department of Matematics, Purdue University, 150 N. University Street, West Lafayette, IN , caiz@purdue.edu. Center for Applied Scientific Computing, Lawrence Livermore National Laboratory, Livermore, CA , falgout2@llnl.gov. Department of Matematics, City University of Hong Kong, Hong Kong, sun.zang@cityu.edu.k. 1

2 Div FOSLL* for 2nd-oder elliptic PDEs 2 te corresponding omogeneous least-squares functionals are equivalent to te Hdiv H 1 and te Hdiv Hcurl H 1 norms, respectively. Te div approac based on te fluxpressure formulation as been studied by many researcers see, e.g., [4, 8, 14, 16]. Te div-curl approac [9] as also been well studied. In order to retain te full efficiency of te L 2 norm first-order system least-squares FOSLS approac wile exibiting te generality of te inverse-norm FOSLS approac, te first-order system LL* FOSLL* approac for general second-order elliptic partial differential equations was proposed and analyzed in [10]. Te FOSLL* approac was applied to te div-curl system, wose adjoint system is an underdetermined system and ence is not suitable for FOSLL*. Tis difficulty was overcome by carefully adding slack variables to te div-curl system. But te resulting approac is quite complicated. Our purpose ere is to study te FOSLL* approac applying to te div system. Witout adding any slack variables to te div system, te resulting approac is muc simpler tan tat in [10]. By sowing tat te bilinear form of te FOSLL* approac is coercive and bounded and tat te linear form is bounded wit respect to a weigted Hdiv H 1 norm, we establis te well-posedness of te FOSLL* approac. Under te same regularity assumption as te standard Galerkin metod, but witout te restriction to sufficiently small mes size, we obtain a quasi-optimal a priori error bound for te corresponding finite element approximation. Note tat tis assumption is weaker tan tat for te div FOSLS [11]. Unlike te FOSLS approac, te FOSLL* approac does not ave a free a posteriori error estimator, tus we study an explicit residual error estimator and establis its reliability and efficiency bounds. Te paper is organized as follows. In Section 2 we introduce matematical equations for te second-order scalar elliptic partial differential equations and its div first-order system, and we ten derive te FOSLL* variational formulation and establis its wellposedness. In Section 3, te FOSLL* finite element approximation is described. A priori and a posteriori error estimations are obtained in Sections 4 and 5 respectively. In Section 6, we present numerical results. 1.1 Notation We use te standard notations and definitions for te Sobolev spaces H s Ω d and H s Ω d for s 0. Te standard associated inner products are denoted by, s,ω and, s, Ω, and teir respective norms are denoted by s,ω and s, Ω. We suppress te superscript d because te dependence on dimension will be clear by context. We also omit te subscript Ω from te inner product and norm designation wen tere is no risk of confusion. For s = 0, H s Ω d coincides wit L 2 Ω d. In tis case, te inner product and norm will be denoted by and,, respectively. Set H 1 DΩ := {q H 1 Ω : q = 0 on Γ D }. Wen Γ = Ω, denote HD 1 Ω by H1 0 Ω. Finally, set wic is a Hilbert space under te norm Hdiv; Ω = {v L 2 Ω d : v L 2 Ω}, v Hdiv; Ω = v 2 + v 2 1/2,

3 Div FOSLL* for 2nd-oder elliptic PDEs 3 and define te subspace H N div; Ω = {v Hdiv; Ω : n v = 0 on Γ N }. 2 First-Order System LL* Formulation Let Ω be a bounded, open, connected subset of R d d = 2 or 3 wit a Lipscitz continuous boundary Ω. Denote by n = n 1,..., n d te outward unit vector normal to te boundary. We partition te boundary of te domain Ω into two open subsets Γ D and Γ N suc tat Ω = Γ D Γ N and Γ D Γ N =. For simplicity, we assume tat Γ D is not empty i.e., mes Γ D 0 and is connected. 2.1 Second-Order Elliptic Problem Consider te following second-order elliptic boundary value problem: wit boundary conditions A u + b u + a u = f in Ω 2.1 u = g D on Γ D and A u n = g N on Γ N, 2.2 were te symbols and stand for te divergence and gradient operators, respectively; A is a given d d tensor-valued function; b L Ω d and a L Ω are given vectorand scalar-valued functions, respectively; and f is a given scalar function. Assume tat A is uniformly symmetric positive definite: tere exist positive constants 0 < Λ 0 Λ 1 suc tat Λ 0 ξ T ξ ξ T Aξ Λ 1 ξ T ξ for all ξ R d and almost all x Ω. Te corresponding variational form of system is to find u H 1 Ω suc tat u ΓD = g D and tat au, v = f, v g D v ds v HDΩ, Γ N were te bilinear form is defined by au, v = A u, v + b u + a u, v. Te dual problem of 2.3 is to find z H 1 Ω suc tat z ΓD = g D and tat âφ, z = f, v g N v ds Γ N φ HDΩ, were te bilinear form is defined by âφ, z = φ, A + b z + φ, a z.

4 Div FOSLL* for 2nd-oder elliptic PDEs 4 Assume tat bot problems 2.3 and 2.4 ave unique solutions and, for simplicity of te presentation, satisfy te full H 2 regularity estimates: u 2 C f and z 2 C f. 2.5 Here and tereafter, we use C wit or witout subscripts in tis paper to denote a generic positive constant, possibly different at different occurrences, tat is independent of te mes size but may depend on te domain Ω. 2.2 First-Order System Introducing te flux vector variable σ = A u, te scalar elliptic problem in may be rewritten as te following first-order partial differential system: { A 1 σ + u = 0 in Ω, 2.6 σ b A 1 σ + a u = f in Ω wit boundary conditions u = g D on Γ D and n σ = g N on Γ N. 2.7 Let A 1 L = b A 1 ten 2.6 may be rewritten as a σ 0, U =, and F =, u f L U = F Div FOSLL* Variational Formulation Multiplying test function V = τ, v t H N div; Ω HD 1 Ω, integrating over te domain Ω, and using integration by parts, we ave f, v = F, V = L U, V = A 1 σ + u, τ + σ b A 1 σ + a u, v = σ, A 1 τ σ, v + σ n v ds σ, A 1 b v Ω u, τ + τ n u ds + u, a v Ω = σ, A 1 τ + A 1 b v + u, a v τ + g N v ds + g D τ n ds Γ N Γ D = U, L V + gτ, v,

5 Div FOSLL* for 2nd-oder elliptic PDEs 5 were te formal adjoint of L and te boundary functional are defined by A 1 + A 1 b L = and gτ, v = g N v ds + g D τ n ds, a Γ N Γ D respectively. Witout loss of generality, we assume tat a 0 in tis paper and let A 0 A =. 0 a 1 Let W = η, w t satisfy U = A L W = η A + b w a 1 η + w, 2.9 ten we ave A L W, L V = f, v gτ, v fτ, v. Now, our div FOSLL* variational formulation is to find η, w H N div; Ω H 1 D Ω suc tat bη, w; τ, v = fτ, v, τ, v H N div; Ω H 1 DΩ, 2.10 were te bilinear form b, is defined by bη, w; τ, v = A L W, L V = η A + b w, A 1 τ + A 1 b v + a 1 η w, τ a v. Note tat bot non-omogenous Diriclet and Neumann boundary conditions are imposed weekly. Remark 2.1. For any η, w, τ, v H N div; Ω HD 1 Ω, integration by parts gives w, τ + w, τ = v, η + v, η = 0. Hence, te bilinear form b, as of te form bη, w; τ, v = A 1 η, τ + a 1 η, τ + A w, v + a w, v b w, A 1 τ A 1 η, b v + b w, v + w, b v + A 1 b w, b v. In te case tat b = 0 and a 0, i.e., te diffusion-reaction problem, te div FOSLL* problem in 2.10 is decoupled. More specifically, w HD 1 Ω is te solution of A w, v + a w, v = f, v g N v ds, v HDΩ, 1 Γ N and η H N div; Ω satisfies A 1 η, τ + a 1 η, τ = g D τ n ds. Γ D

6 Div FOSLL* for 2nd-oder elliptic PDEs 6 Note tat te problem for w is similar to te standard variational formulation for te diffusion-reaction problem, but te non-omogeneous Diriclet boundary condition is weakly imposed ere. Remark 2.2. In te case tat a = 0, let A = diag A, 1, ten W = η, w t satisfy η A + b w U = A L W = η Te corresponding bilinear form is modified as follows bη, w; τ, v = η A + b w, A 1 τ + A 1 b v + η, τ. 2.4 Well-Posedness Denote by v 1 = a 1/2 v 2 + A 1/2 v 2 1/2 and τ Hdiv = A 1/2 τ 2 + a 1/2 τ 2 1/2 te weigted H 1 Ω and Hdiv; Ω norms, respectively. Let τ, v = v τ 2 Hdiv 1/2. Te following teorem establises te coercivity and continuity of te bilinear form. Teorem 2.3. Te bilinear form b, is coercive and continuous in H N div; Ω H 1 D Ω, i.e., tere exist positive constant α and C, depending on bounds of te coefficients A, b, and a, suc tat α τ, v 2 bτ, v; τ, v 2.12 and tat for all η, u, τ, v H N div; Ω H 1 D Ω. bη, w; τ, v C η, w τ, v 2.13 A similar result to tat of Teorem 2.3 was proved in [8]. readers, we provide a compreensive proof ere. For te convenience of Proof is a direct consequence of te Caucy-Scwarz and te triangle inequalities and te bounds of te coefficients A, b, and a of te underlying problem. To sow te validity of 2.12, we first establis tat τ, v 2 C bτ, v; τ, v + v for all τ, v H N div; Ω HD 1 Ω. To tis end, integrating by parts gives τ, v = a 1/2 τ + a 1/2 v, a 1/2 v a v, v.

7 Div FOSLL* for 2nd-oder elliptic PDEs 7 It ten follows from te Caucy-Scwarz inequality tat = A 1/2 v 2 + a 1/2 v 2 A 1/2 v A 1/2 τ b v, A 1/2 v + a 1/2 v a 1/2 τ, a 1/2 v b v, v A 1/2 v A 1/2 τ b v + C v A 1/2 v + a 1/2 τ a 1/2 v a 1/2 v, wic implies A 1/2 v 2 + a 1/2 v 2 C bτ, v; τ, v + v By te triangle inequality and 2.15, we ave tat A 1/2 τ A 1/2 τ b v A 1/2 v + A 1/2 v + C v and tat C bτ, v; τ, v 1/2 + v. a 1/2 τ a 1/2 τ a 1/2 v + a 1/2 v C bτ, v; τ, v + v 2. Combining te above tree inequalities yields Wit 2.14, we now sow te validity of 2.12 by te standard compactness argument. To tis end, assume tat 2.12 is not true. Tis implies tat tere exists a sequence {τ n, v n } H N div; Ω HD 1 Ω suc tat τ n 2 Hdiv + v n 2 1 = 1 and bτ n, v n ; τ n, v n 1 n Since H 1 D Ω is compactly contained in L2 Ω, tere exists a subsequence {v nk } H 1 D Ω wic converges in L 2 Ω. For any k, l and τ nk, v nk, τ nl, v nl H N div; Ω H 1 D Ω, it follows from 2.14 and te triangle inequality tat τ nk τ nl 2 Hdiv + v nk v nl 2 1, Ω C bτ nk τ nl, v nk v nl ; τ nk τ nl, v nk v nl + v nk v nl 2 C bτ nk, v nk ; τ nk, v nk + bτ nl, v nl ; τ nl, v nl + v nk v nl 2 0, as k, l. Tis implies tat τ nk, v nk is a Caucy sequence in te complete space H N div; Ω HD 1 Ω. Hence, tere exists τ, v H Ndiv; Ω HD 1 Ω suc tat lim τ nk τ Hdiv + v nk v 1 = 0. k Next, we sow tat v = 0 and τ = 0, 2.17

8 Div FOSLL* for 2nd-oder elliptic PDEs 8 wic contradict wit 2.16 tat 0 = τ 2 Hdiv + v 2 1 = lim τ nk 2 Hdiv + v nk 2 1 = 1. k To tis end, for any φ HD 1 Ω, integration by parts, te Caucy-Scwarz inequality, and 2.16 give âφ, v nk = φ, A + b v nk + φ, a v nk = φ, A + b v nk τ nk + φ, a v nk τ nk 1 1/2 bτ nk, v nk ; τ nk, v nk 1/2 φ 1 φ 1. n k Since lim k v nk = v in H 1 Ω, we ten ave 1 1/2 âφ, v = lim âφ, v n k lim φ 1 = 0. k k n k Because 2.4 as a unique solution, we ave tat v = 0. Now, τ = 0 follows from 2.14: τ 2 Hdiv = lim τ nk 2 Hdiv C lim bτ nk, v nk ; τ nk, v nk + v nk 2 = 0. k k Tis completes te proof of 2.17 and, ence, te teorem. Teorem 2.4. Te variational formulation in 2.10 as a uniques solution η, w H N div; Ω HD 1 Ω satisfying te following a priori estimate η, w C f 1,Ω + g D 1/2,ΓD + g N 1/2,ΓN Proof. For all τ, v H N div; Ω HD 1 Ω, it follows from te definition of te dual norms and te trace teorem tat fτ, v f 1,Ω v 1,Ω + g D 1/2,ΓD τ n 1/2,ΓD + g N 1/2,ΓN v 1/2,ΓD C f 1,Ω + g D 1/2,ΓD + g N 1/2,ΓN τ, v. Now, by te Lax-Milgram lemma, te well possedness of 2.10 and te a priori estimate in 3.4 follow directly from Teorem 2.3. Remark 2.5. In te case tat a = 0, te norms are modified as follows v 1 = A 1/2 v 2 1/2 and τ Hdiv = A 1/2 τ 2 + τ 2 1/2. Wit tese norms, all te results obtained in tis paper for a 0 old.

9 Div FOSLL* for 2nd-oder elliptic PDEs 9 3 Div FOSLL* Finite Element Approximation Teorem 2.3 guarantees tat confirming finite element spaces of H N div; Ω HD 1 Ω for te vector and scalar variables, η and w, may be cosen independently. However, te only finite element spaces aving optimal approximations in terms of bot te regularity and te approximation property are te continuous piecewise polynomials for te scalar variable and te Raviart-Tomas or Brezzi-Douglas-Marini elements for te vector variable. Te BDM element as sligt more degrees of freedom tan tat of te RT element. Moreover, te system of algebraic equations resulting from tese elements can be solved efficiently by fast multigrid metods. For te above reasons, only tese elements are analyzed in tis paper. But it is easy to see tat our analysis does apply to any oter conforming finite element spaces wit no essential modifications. For simplicity of presentation, we consider only triangular and tetraedra elements for te respective two and tree dimensions. Assuming tat te domain Ω is polygonal, let T be a regular triangulation of Ω see [13] wit triangular/tetraedra elements of size O. Let P k K be te space of polynomials of degree k on triangle K and denote te local Raviart-Tomas space of order k on K: RT k K = P k K d + x P k K wit x = x 1,..., x d. Ten te standard Hdiv; Ω conforming Raviart-Tomas space of index k [17] and te standard conforming continuous piecewise polynomials of degree k + 1 are defined, respectively, by Σ k = {τ H N div; Ω : τ K RT k K K T }, and V k+1 = {v H 1 DΩ : v K P k+1 K K T } It is well-known see [13] tat V k+1 as te following approximation property: let k 0 be an integer and let l [0, k + 1] inf v V k+1 u v 1 C l u l+1, 3.1 for u H l+1 Ω HD 1 Ω. It is also well-known see [17] tat Σk approximation property: let k 0 be an integer and let l [1, k + 1] as te following inf σ τ Hdiv; Ω C l σ l + σ l 3.2 τ Σ k for σ H l Ω d H N div; Ω wit σ H l Ω m. Since σ and σ are one order less smoot tan u, we will coose k to be te smallest integer greater tan or equal to l 1. Te finite element discretization of te FOSLL* variational problem is: find η, w Σ k V k+1 suc tat bη, w ; τ, v = fτ, v, τ, v Σ k V k Since Σ k V k+1 is a subspace of H N div; Ω HD 1 Ω, te div FOSLL* problem in 3.3 is well-posed and te solution continuously depends on te data.

10 Div FOSLL* for 2nd-oder elliptic PDEs 10 Teorem 3.1. Te variational formulation in 3.3 as a uniques solution η, w Σ k V k+1 satisfying te following a priori estimate η, w C f 1,Ω + g D 1/2,ΓD + g N 1/2,ΓN. 3.4 Now, te finite element approximation to σ, u is defined as follows: σ = η A w b w and u = a 1 η + w. 3.5 Remark 3.2. Wen te coefficients A, b, and a are not polynomials, tey can be replaced by teir approximations of appropriate polynomials locally, if piecewise polynomial approximation to σ, u is desirable. Remark 3.3. Te FOSLL* approximation to te solution u is not continuous. To obtain a continuous approximation, one can simply project u onto appropriate continuous finite element space. Remark 3.4. In te case tat a = 0, te finite element approximation to σ, u is given by σ = η A w b w and u = η. 4 A Priori Error Estimate Difference between equations in 2.10 and 3.3 gives te error equation: bη η, w w ; τ, v = 0 τ, v Σ k V k Te following error estimation in te energy norm is a simple consequence of Teorem 3.1, te error equation in 4.6, te Caucy-Scwarz inequality, and te approximation properties in 3.1 and 3.2. Teorem 4.1. Assume tat te solution σ, u of is in H l Ω d H l+1 Ω and tat te solution η, w of 2.10 satisfies η l + η l + w l+1 C A 1/2 σ l + a 1/2 u l. 4.7 Let k be te smallest integer greater tan or equal to l 1. Ten te FOSLL* approximation σ, u defined in 3.5 as te following error estimate A 1/2 σ σ + a 1/2 u u C l σ l + u l C l u l Proof. Let η, w be te solution of 3.3. It follows from Teorem 2.1, te error equation in 4.6, and te approximation properties in 3.1 and 3.2 tat η η, w w C inf η τ Hdiv; Ω + inf w v 1 τ Σ k v V k+1 C l η l + η l + w l+1, wic, togeter wit 4.7, implies 4.8. Tis completes te proof of te teorem.

11 Div FOSLL* for 2nd-oder elliptic PDEs 11 Remark 4.2. Wen l = 0, Assumption 4.7 is te coercivity bound in Remark 4.3. In te case tat a = 0, 4.8 becomes A 1/2 σ σ + u u C l σ l + u l+1 C l u l A Posteriori Error Estimate Unlike te FOSLS approac, te FOSLL* approac does not ave a free a posteriori error estimator, tus in tis section we study an explicit residual error estimator and establis its reliability and efficiency bounds. 5.1 Local Indicator and Global Estimator Since te bilinear form b, is coercive and continuous in H N div; Ω HD 1 Ω see Teroem 2.3, te explicit residual a posteriori error estimator to be derived in tis paper is a combination of tose for te Hdiv and te elliptic problems see [12, 1, 18, 19]. To tis end, we first introduce some notations. Denote by E K te set of edges/faces of element K T and te set of edges/faces of te triangulation T by E := E I E D E N, were E I is te set of interior element edges, and E D and E N are te sets of boundary edges belonging to te respective Γ D and Γ N. For eac e E, denote by e te lengt/diameter of te edge/face e and by n e a unit vector normal to e. Let Ke and K e + be te two elements saring te common edge/face e suc tat te unit outward normal vector of Ke coincides wit n e. Wen e E D E N, n e is te unit outward vector normal to Ω and denote by Ke te element aving te edge/face e. For a function v defined on Ke K e +, denote its traces on F by v e and v + e, respectively. Te jump over te edge/face e is denoted by { v e v + e e E I, [v ] e := v e e E D E N. Wen tere is no ambiguity, te subscript or superscript e in te designation of te jump will be dropped. For a function v, we will use te following notations on te weigted L 2 norms: 1/2 v T = v 2 K were v K = K v 0,K K T, K T 1/2 and v E = v 2 e were v e = e v 0,e e E. e E Let η, w be te solution of 3.3 and let σ, u be te finite element approximation to σ, u defined in 3.5. On eac element K T, denote te following element residuals by r 1 K = f σ + b A 1 σ a u, r 2 K = A 1 σ + u, and r 3 K = A 1 σ.

12 Div FOSLL* for 2nd-oder elliptic PDEs 12 Denote te following edge jumps by J 1 e = [σ n], J 2 e = [u ], J 3 e = [A 1 σ t], on e E I, J 1 e = 0, J 2 e = u g D, J 3 e = g D t + A 1 σ t, on e E D, J 1 e = σ n g N, J 2 e = 0, J 3 e = 0, on e E N. Let r 2 K, r 1 K and r 3 K, and J i e i = 1, 2, 3 be te L 2 -projections of te respective r 2 K, r 1 K and r 3 K, and J i e i = 1, 2, 3 onto P k K 2, P k K, and P k e, respectively. Now, te local error indicator on eac element K T is defined by η 2 = r K 1 2 K + r 2 2 K + r 3 2 K + 1 1/2 J1 2 e + 1/2 J3 2 e 2 + e E D E K 1/2 J3 2 e + and te global error estimator is defined by e E I E K e E N E K 1/2 J1 2 e, 5.10 η 2 = K T η 2 K = r 1 2 T + r 2 2 T + r 3 2 T + 1/2 J1 2 E + 1/2 J3 2 E Te terms r 1 and r 2 are te residuals of te equations in 2.6. Te term r 3 measures te violation of te fact tat te exact quantity A 1 σ = u is in te kernel of operator. Te terms J 1, J 2, and J 3 are due to te facts tat te numerical flux σ, te numerical solution u, and te numerical gradient A 1 σ are not in Hdiv; Ω, H 1 Ω, and Hcurl; Ω, respectively. 5.2 Reliability and Efficiency Bounds For simplicity, we analyze only two dimensions ere since tere is no essential difficulties for tree dimensions. For a vector field τ = τ 1, τ 2 t and a scalar-value function v, define te respective curl operator and its formal adjoint by τ := τ 2 x 1 τ 1 x 2 and v := v x 1, v x 1 t. Denote by Π : H N div; Ω L s Ω 2 Σ 0 wit s > 2 te RT 0 interpolation operator; i.e., for all τ H 1 Ω 2, one as [5] τ Π τ n, v e = 0, v P 0 e and e E Let S be te standard continuous piecewise linear finite element space on te triangulation T. For B = D or N, denote by I B : HB 1 Ω H1 B Ω S te Clement interpolation operator wic satisfies te following local approximation property [6]: 1 v I B v T C v v H 1 BΩ, were B = D or N. It is easy to ceck tat I N v Σ 0.

13 Div FOSLL* for 2nd-oder elliptic PDEs 13 Let η, w be te solution of 2.10, and let E = η η H N div; Ω and e = w w H 1 DΩ By Lemmas 5.1 in [12], te E as te following quasi-helmoltz decomposition: E = φ + ψ in Ω, 5.14 were φ HN 1 Ω2 and ψ HN 1 Ω. Moreover, tere exists a constant C > 0 suc tat Let and let φ C E Hdiv and ψ C E Hdiv φ := Π φ, ψ := I N ψ, and e := I D e ẽ = e e, φ = φ φ, and ψ = ψ ψ. By te approximation properties of te interpolation operators and 5.15, we ave 1 ẽ T + 1/2 ẽ E C e C e 1, 1 φ T + 1/2 φ E C φ C E Hdiv, and 1 ψ T + 1/2 ψ E C ψ C E Hdiv. Lemma 5.1. We ave te following error representation: be, e; E, e = { r1, ẽ K r 2, φ K r 3, ψ } { K + J1, ẽ e + J 2, φ n e + J 3, ψ } e. K T e E 5.16 Proof. Since E := φ + ψ Σ 0 Σk and e V 1 V k+1, te error equation in 4.6 gives be, e; E, e = be, e; E E, e e = be, e; φ + ψ, ẽ By te fact tat ψ = 0, te definitions of te bilinear form b,, and te FOSLL* finite element approximation σ, u, we ave bη, w ; φ + ψ, ẽ = σ, A 1 φ + ψ ẽ A 1 b ẽ u, φ a ẽ = A 1 σ, φ + ψ u, φ σ, ẽ b t A 1 σ a u, ẽ. It follows from integration by parts and te boundary conditions tat A 1 σ, ψ = A 1 σ, ψ K [A 1 σ t], ψ e A 1 σ t, ψ e, e E I e E D K T u, φ = K T u, φ K + e EI σ, ẽ = K T σ, ẽ K + e EI [u ], φ n e + [σ n], ẽ e + e E D u, φ n e, e E N σ n, ẽ e,

14 Div FOSLL* for 2nd-oder elliptic PDEs 14 wic, togeter wit te definitions of te residuals and te jumps, lead to bη, w ; φ + ψ, ẽ = K T { r3, ψ K + r 2, φ K r 1, ẽ K } + K T f, ẽ K e E I J 3, ψ e A 1 σ t, ψ e 2, e E D e E J φ n e J 1, ẽ e g D, φ n e g N, ẽ e. e E e E D e E N By 2.10, integration by parts, and te boundary condition of ψ HN 1 Ω, we ave bη, w; φ + ψ, ẽ = f, ẽ gn, ẽ ΓN g D, φ + ψ n ΓD = f, ẽ g N, ẽ ΓN g D, φ n ΓD + g D t, ψ ΓD. Now, 5.16 is a direct consequence of 5.17 and te difference of te above two equalities. Tis completes te proof of te lemma. Define te local and global oscillations as follows: osc 2 k K = r 1 r 1 2 K + r 2 r 2 2 K + r 3 r 3 2 K + 1/2 J 1 J 1 2 K + 1/2 J 2 J 2 2 K + 1/2 J 3 J 3 2 K 5.18 and osc 2 k T = osc 2 k K, 5.19 K T respectively. Teorem 5.2. Reliability Bound Te global estimator η defined in 5.11 is reliable; i.e., tere exists a positive constant C suc tat A 1/2 σ σ + a 1/2 u u C E, e C η + osc k T Proof. Te first inequality in 5.20 is a direct consequence of te definition of σ, u and te triangle inequality. To sow te validity of te second inequality in 5.20, by te coercivity in 2.12, it suffices to sow tat be, e; E, e C η + osc k T E, e To tis end, first notice tat by te property in 5.12 and te definition of J 2, we ave J 2, φ φ n e = J 2 J 2, φ φ n e.

15 Div FOSLL* for 2nd-oder elliptic PDEs 15 Now, it follows from Lemma 5.1, te Caucy-Scwarz inequality, and te approximation properties of ẽ, φ, and ψ, and te triangle inequality tat be, e; E, e = r1, ẽ K + r 2, φ K + r 3, ψ K + J1, ẽ e + J 2 J 2, φ n e + J 3, ψ e K T K T r1 K 1 ẽ K + r 2 K 1 φ K + r 3 K 1 ψ K + e E 1/2 J 1 e 1/2 ẽ e + 1/2 J 2 J 2 e 1/2 φ n e + 1/2 J 3 e 1/2 ψ e e E C ri 2 T + 1/2 J i 2 E + r2 2 T + 1/2 J 2 J 2 2 E i=1,3 1/2 E, e C η + osc k T E, e, wic proves 5.21 and, ence, te teorem. Teorem 5.3. Local Efficiency Bound For all K T, te local error indicator η K defined in 5.10 is efficient; i.e., tere exists a positive constant C suc tat C η K σ σ ωk + u u ωk + osc k ω K, 5.22 were ω K is te union of elements in T saring an edge wit K. Te proof of te local efficiency bound in Teorem 5.3 is standard; i.e., it is proved by using local edge and element bubble functions, φ e and φ K see [18] for teir definitions and properties. For simplicity, we only sketc te proof below. Proof. For any τ, v H N div; Ω HD 1 Ω, by te quasi-helmoltz decomposition, we ave τ = φ + ψ H N div; Ω, were φ HN 1 Ω and ψ H1 N Ω. Te same argument as te proof of Lemma 5.1 gives be, e; τ, v = K T {r 1, v K + r 2, φ K + r 3, ψ K }+ e E { J 1, v e + J 2, φ n e + J 3, ψ e }, wic, togeter wit te definitions of σ and u, yields be, e; τ, v = σ σ, A 1 τ + A 1 b v u u, τ a v. Hence, = r 1, v K + r 2, φ K + r 3, ψ K + J 1, v e + J 2, φ n e + J 3, ψ e K T e E σ σ, A 1 φ + ψ + A 1 b v u u, φ a v. 5.23

16 Div FOSLL* for 2nd-oder elliptic PDEs 16 In 5.23, by coosing 1 φ = 0, ψ = 0, and v = r 1 φ K ; 2 φ = r 2 φ K, ψ = 0, and v = 0; and 3 φ = 0, ψ = r 3 φ K, and v = 0 and by te standard argument, we can ten establis upper bounds for te element residuals, r 1 K, r 2 K, and r 3 K, respectively. In a similar fasion, to bound te edge jumps J 1 e and J 3 e above, we coose 1 φ = 0, ψ = 0, and v = J 1 φ e and 2 φ = 0, ψ = J 3 φ e, and v = 0 in 5.23, respectively. 6 Numerical Results In tis section, numerical results for a second order elliptic partial differential equation are presented. We begin wit discretizations of a test problem on a sequence of uniform meses to verify te a priori error estimation. Te test problem is defined on Ω = 0, 1 2 wit coefficients A = I, b = 3, 2 t, and a = 2. Te exact solution of tis problem is u = sinπx sinπy wit omogeneous boundary condition on Ω. Finite element spaces Σ 0 and V 1 are used to approximate τ and w, respectively. Table 6.1 sows tat te convergence rates of te errors for te original variables σ and u in L 2 -norms are optimal. Table 1: Errors σ σ, u u, σ σ + u u, and convergence rates σ σ rate u u rate σ σ + u u rate 1/ E E E-1 1/ E E E / E E E / E E E / E E E Te next example is to test te a posteriori error estimator. Te test problem is te Laplace equation u = 0 defined on an L-saped domain Ω := 1, 1 2 \[0, 1 1, 0] wit a reentrant corner at te origin. Te Diriclet boundary condition u Ω = g D is cosen suc tat te exact solution is ur, θ = r 2/3 sin2θ in polar coordinates. Starting wit te coarsest triangulation T 0 obtained from alving 12 uniform squares, a sequence of meses is generated by using te standard adaptive mesing algoritm tat adopts te bulk marking strategy wit bulk marking parameter 0.5. Marked triangles are refined by bisection. Mes generated by η is sown in Figure 6.1. Te refinement is mainly around te reentrant corner. Comparison of te true error and te η is sown in Figure 6.2. Moreover, te slope of te logdof- logerror for te η and te true error is 1/2, wic indicates te optimal decay of te error wit respect to te number of unknowns.

17 Div FOSLL* for 2nd-oder elliptic PDEs Convergence of error error estimator N 0.5 error number of nodes Figure 1: mes generated by error estimator Figure 2: error σ σ + u u and estimator References [1] M. Ainswort and J. T. Oden, A Posteriori Error Estimation in Finite Element Analysis, Jon Wiley & Sons, Inc., [2] A. Aziz, R. Kellogg, and A. Stepens, Least-squares metods for elliptic systems, Mat. Comp., , [3] P. B. Bocev and M. D. Gunzburger, Least Squares Finite Element Metods, Springer, Berlin, [4] P. B. Bocev and M. D. Gunzburger, On least-squares finite element metods for te Poisson equation and teir connection to te Diriclet and Kelvin principles, SIAM J. Numer. Anal., 43:1 2005, [5] D. Boffi, F. Brezzi, and M. Fortin, Mixed Finite Element Metods and Applications, Springer Series in Computational Matematics, 44, Springer, [6] D. Braess, Finite Elements: Teory, Fast Solvers and Applications in Solid Mecanics, 3rd edition, Cambridge University Press, Cambridge, UK, [7] J. H. Bramble, R. D. Lazarov, and J. E. Pasciak, A least-squares approac based on a discrete minus one inner product for first order system, Mat. Comp., , [8] Z. Cai, R. Lazarov, T.A. Manteuffel and S.F. McCormick, First-order system least squares for second-order partial differential equations: part I., SIAM J. Numer. Anal., , , 2.4 [9] Z. Cai, T.A. Manteuffel and S.F. McCormick, First-order system least squares for second-order partial differential equations: Part II, SIAM J. Numer. Anal., 34:2 1997,

18 Div FOSLL* for 2nd-oder elliptic PDEs 18 [10] Z. Cai, T.A. Manteuffel, S.F. McCormick, and J. Ruge, First-order system LL* FOSLL* for scalar partial differential equations, SIAM J. Numer. Anal., , document, 1 [11] Z. Cai and Ku, Optimal error estimates for te div least-squares metod wit data f in L2 and application to nonlinear problems, SIAM J. Numer. Anal., 47:6 2010, [12] J. M. Cascon, R. H. Nocetto, and K. G. Siebert, Design and convergence of afem in Hdiv, Mat. Models and Metods in Applied Sciences, , , 5.2 [13] P. G. Ciarlet, Te Finite Element Metod for Elliptic Problems, Nort-Holland, Amsterdam, [14] D. C. Jespersen, A least-square decomposition metod for solving elliptic systems, Mat. Comp., , [15] B. N. Jiang, Te Least-Squares Finite Element Metod: Teory and Applications in Computational Fluid Dynamics and Electromagnectics, Spring, Berlin, [16] A. I. Pelivanov, G. F. Carey, and R. D. Lazarov, Least squares mixed finite elements for second order elliptic problems, SIAM J. Numer. Anal., , [17] P. A. Raviart and J. M. Tomas, A mixed finite element metod for 2-nd order elliptic problems, Matmatical Aspects of Finite Element Metods, Lecture Notes in Matematics, #606, I. Galligani and E. Magenes, eds., Springer-Verlag, New York, 1977, , 3 [18] R. Verfürt, A Review of A-Posteriori Error Estimation and Adaptive Mes Refinement Tecniques, Jon Wiley and Teubner Series, Advances in Numerical Matematics, , 5.2 [19] R. Verfürt, A Posteriori Error Estimation Tecniques for Finite Element Metods, Oxford Numerical Matematics and Scientific Computation, Oxford University Press,