Chapter 6 A posteriori error estimates for finite element approximations 6.1 Introduction

Size: px
Start display at page:

Download "Chapter 6 A posteriori error estimates for finite element approximations 6.1 Introduction"

Transcription

1 Chapter 6 A posteriori error estimates for finite element approximations 6.1 Introduction The a posteriori error estimation of finite element approximations of elliptic boundary value problems has reached some state of maturity, as it is documented by a variety of monographs on this subject (cf., e.g., [1, 2, 3, 4, 5]). There are different concepts such as residual type a posteriori error estimators, hierarchical type a posteriori error estimators, error estimators based on local averaging, error estimators based on the goal oriented weighted dual approach (cf., in particular, [3]). In this chapter, we will focus on residual type a posteriori error estimators and follow the exposition in [5]. We shall deal with the following model problem: Let Ω be a bounded simply-connected polygonal domain in uclidean space lr 2 with boundary Γ = Γ D Γ N, Γ D Γ N = and consider the elliptic boundary value problem (6.1) Lu := div (a grad u) = f in Ω, u = 0 on Γ D, n a grad u = g on Γ N, where f L 2 (Ω), g H 1/2 (Γ N ) and a = (a ij ) 2 i,j=1 is supposed to be a matrix-valued function with entries a ij L (Ω), that is symmetric, i.e., a ij (x) = a ji (x) f.a.a. x Ω, 1 i, j 2, and uniformly positive definite in the sense that for almost all x Ω 2 a ij (x) ξ i ξ j α ξ 2, ξ R 2, α > 0. i,j=1 The vector n denotes the exterior unit normal vector on Γ N. We further set ᾱ := max 1 i,j 2 a ij. Setting H 1 0,Γ D (Ω) := { v H 1 (Ω) v ΓD = 0 }, the weak formulation of (6.1) is as follows: Find u H 1 0,Γ D (Ω) such that (6.2) a(u, v) = l(v), v H0,Γ 1 D (Ω), 1

2 2 Ronald H.W. Hoppe where a(v, w) := a grad v grad w dx, v, w H 1 0,Γ D (Ω), l(v) := Ω Ω f v dx + Γ N g v dσ, v H 1 0,Γ D (Ω). Given a geometrically conforming simplicial triangulation T h of Ω, we denote by S 1,ΓD (Ω; T h ) := { v h H 1 0,Γ D (Ω) v h P 1 (), T h } the trial space of continuous, piecewise linear finite elements with respect to T h. Note that P k (), k 0, denotes the linear space of polynomials of degree k on. In the sequel we will refer to N h (D) and h (D), D Ω as the sets of vertices and edges of T h on D. The conforming P1 approximation of (6.1) reads as follows: Find u h S 1,ΓD (Ω; T h ) such that (6.3) a(u h, v h ) = l(v h ), v h S 1,ΓD (Ω; T h ). Now, assuming that ũ h S 1,ΓD (Ω; T h ) is some iterative approximation of u h S 1,ΓD (Ω; T h ), we are interested in the total error e := u ũ h, which is the sum of the discretization error e d := u u h and the iteration error e it := u h ũ h. It is easy to see that the total error e is in V := H 1 0,Γ D (Ω) and satisfies the error equation (6.4) a(e, v) = r(v), v V, where r( ) stands for the residual with respect to the computed approximation ũ h (6.5) r(v) := f v dx + g vdσ a(ũ h, v), v V. Ω Γ N We are interested in a cheaply computable a posteriori error estimator η of the total error e consisting of elementwise error contributions η, T h and edgewise error contributions η, h, in the sense that (6.6) η 2 = T h η 2 + h η 2, which does provide a lower and an upper bound for e according to (6.7) γ η e 1,Ω Γ η

3 Finite lement Methods 3 with constants 0 < γ Γ depending only on the ellipticity constants and on the shape regularity of the triangulation T h. We may use the local error terms η and η as a criterion for local refinements of the elements T h. Among several refinement strategies, the so-called mean-value strategy is as follows: Compute the mean values η := 1 (6.8) η, n T h η := 1 n h η, where n := card T h and n := card h. Mark an element T h and an edge h for refinement, if (6.9) η σ η, η σ η, where 0 < σ 1 is some appropriate safety factor, e.g., σ = 0.9. Definition 6.1 fficient and reliable a posteriori error estimators An a posteriori error estimator η satisfying (6.10) e 1,Ω Γ η is called reliable, since it ensures a sufficient refinement in the sense that the H 1 -norm of the total error e will be bounded by a quantity of the same order of magnitude as a user-prescribed accuracy, if this accuracy is tested by η. On the other hand, an a posteriori error estimator η for which (6.11) γ η e 1,Ω is said to be efficient, since it underestimates the H 1 -norm of the total error e and thus prevents too much refinement. 6.2 Residual based a posteriori error estimators The residual based a posteriori error estimator can be derived by viewing the residual as an element of the dual space V and evaluating it with respect to the dual norm Upper bound for the total error An important tool in the construction of an upper bound for the total error is Clément s interpolation operator which is defined as follows:

4 4 Ronald H.W. Hoppe Definition 6.2 Clément s interpolation operator For p N h (Ω) N h (Γ N ) we denote by ϕ p the basis function in S 1,ΓD (Ω; T h ) with supporting point p and we refer to D p as the set (6.12) D p := { T h p N h (T ) }.0 D p p Fig. 6.1: Clément s interpolation operator (definition) We refer to π p as the L 2 -projection onto P 1 (D p ), i.e., (π p (v), w) 0,Dp = (v, w) 0,Dp, w P 1 (D p ), where (, ) 0,Dp stands for the L 2 -inner product on L 2 (D p ) L 2 (D p ). Then, Clément s interpolation operator P C is defined as follows (6.13) P C : L 2 (Ω) S 1,ΓD (Ω, T h ), P C v := π P (v) ϕ P. p N h (Ω) N h (Γ N ) In order to establish local approximation properties of Clément s interpolation operator, for T h and h (Ω) h (Γ N ) we introduce the sets (6.14) (6.15) D (1) := { T h N h ( ) N h () }, D (1) := { T h N h () N h ( ) }.

5 Finite lement Methods 5 Fig. 6.2: Clément s interpolation operator (properties) Using the affine equivalence of the elements and the Bramble- Hilbert Lemma one can show: Theorem 6.1 Approximation properties of Clément s interpolation operator For T h and h (Ω) h (Γ N ) let D (1) and D(1) be given by (6.14) and let P C be Clément s interpolation operator as given by (6.13). Then, there exist constants C ν > 0, 1 ν 5, depending only on the shape regularity of T h such that for all v S 1,ΓD (Ω; T h ): (6.16) (6.17) (6.18) (6.19) (6.20) P C v 0, C 1 v (1) 0,D, P C v 0, C 2 v (1) 0,D, grad P C v 0, C 3 grad v (1) 0,D v P C v 0, C 4 h v (1) 1,D, v P C v 0, C 5 h 1/2 v 1,D (1),, where h := diam and h :=. Further, there exist constants C 6, C 7 > 0 depending only on the shape regularity of T h such that (6.21) (6.22)( h (Ω) h (Γ N ) ( v 2 ) 1/2 C µ,d (1) 6 v µ,ω, 0 µ 1, T h v 2 ) 1/2 C µ,d (1) 7 v µ,ω, 0 µ 1. Proof. We refer to [5]. We have now provided all prerequisites to establish an upper bound for the total error e measured in the H 1 -norm. For functions v h W 0 (Ω; T h ) we further refer to [v h ] J as the jump across the common

6 6 Ronald H.W. Hoppe edge h (Ω) of two adjacent elements 1, 2 T h [v h ] J := v h 1 v h 2. Theorem 6.2 Upper bound for the total error There exist constants Γ R, Γ osc > 0 and Γ it > 0 depending only on the elipticity constants and the shape regularity of T h such that (6.23) e 1,Ω Γ R η R + Γ osc osc + η it e it 1,Ω, where the element and edge residuals are given by η R := 3 ν=1 η (ν) R, η (1) R := ( T h h 2 T π h f Lũ h 2 0,) 1/2, η R := ( h (Γ N ) η (3) R := ( h (Ω) h π h g n a grad ũ h 2 0,) 1/2, h [n a grad ũ h ] J 2 0,) 1/2, and osc stands for the data oscillations osc := ( T h osc 2 + h (Γ N ) osc 2 ) 1/2, osc := h T f π h f 0,, osc := h g π h g 0,. Proof. Setting v = e in (6.4), we obtain (6.24) α e 2 1,Ω a(e, e) = r(e) = r(p C e) + r(e P C e). Taking advantage of (6.3), for the first term on the right-hand side of (6.24) we get r(p C e) = f P C e dx + g P C e dσ a(ũ h, P C e) = = Ω Γ N a (u h ũ h, P C e). T h

7 Finite lement Methods 7 Using (6.19), the Schwarz inequality, and observing (6.22), it follows that (6.25) r(p C e) α C 3 u h ũ h 1, e (1) 1,D T h α C 3 ( u h ũ h 2 1,) 1/2 ( e 2 ) 1/2 1,D (1) T h T h α C 3 C 6 e it 1,Ω e 1,Ω. On the other hand, for the second term on the right-hand side of (6.24), Green s formula yields r(e P C e) = f (e P C e) dx + g (e P C e) dσ + Ω Γ N + div a grad ũ h (e P }{{} C e) dx T h = Lũ h n a grad ũ h (e P C e) dσ = T h = T h (π h f Lũ h ) (e P C e) dx + + h (Ω) + h (Γ N ) [n a grad ũ h ] J (e P C e) dσ + (π h g n a grad ũ h ) (e P C e) dσ + + (f π h f) (e P C e) dx + T h + h (Γ N ) (g π h g) (e P C e) dσ. In view of (6.17),(6.18) and (6.22),(6.23), it follows that (6.26) r(e P C e) C 1 C 6 ( h 2 π h f Lũ h 2 0,) 1/2 e 1,Ω + T h + C 2 C 7 ( h [n a grad ũ h ] J 2 0,) 1/2 e 1,Ω + h (Ω)

8 8 Ronald H.W. Hoppe + C 2 C 7 ( h (Γ N ) h π h g n a grad ũ h 2 0,) 1/2 e 1,Ω + + C 1 C 6 ( T h h 2 f π h f 2 0,) 1/2 e 1,Ω + + C 2 C 7 ( h (Γ N ) h g π h g 2 0,) 1/2 e 1,Ω. Using (6.25),(6.26) in (6.24), the assertion follows with Γ R = Γ osc := α 1 max(c 1 C 6, C 2 C 7 ) and Γ it := α 1 α C 3 C 6. For the construction of a lower bound we will now show that the local contributions η (ν) R, := η(ν) R, T h, 1 ν 3 of the residual based error estimator η R do locally provide lower bounds for the total error e. For this purpose we need appropriate localized polynomial functions defined on the elements of the triangulation and the edges h (Ω) h (Γ N ), respectively. Such functions are given by the triangle-bubble functions ψ and the edge-bubble functions ψ. In particular, denoting by λ i, 1 i 3, the barycentric coordinates of T h, then the triangle-bubble function ψ is defined by means of (6.27) ψ := 27 λ 1 λ 2 λ 3. Note that supp ψ = int, i.e., ψ = 0, T h. On the other hand, for h (Ω) h (Γ N ) and T h such that and p i N h (T ), 1 i 2, we introduce the edge-bubble functions ψ according to (6.28) ψ := 4 λ 1 λ 2. Note that ψ = 0 for h (),. The bubble functions ψ and ψ have the following important properties that can be easily verified taking advantage of the affine equivalence of the elements:

9 Finite lement Methods 9 Lemma 6.1 Basic properties of the bubble functions. Part I There exist constants C ν > 0, 8 ν 12, depending only on the shape regularity of the triangulations T h such that (6.29) p h 2 0, C 8 p 2 h ψ dx, p h P 1 (), (6.30) p h 2 0, C 9 p 2 h ψ dσ, p h P 1 (), (6.31) (6.32) (6.33) p h ψ 1, C 10 h 1 p h 0,, p h P 1 (), p h ψ 0, C 11 p h 0,, p h P 1 (), p h ψ 0, C 12 p h 0,, p h P 1 (). For functions p h P 1 (), h (Ω) h (Γ N ) we further need an extension p h L2 () where T h such that. For this purpose we fix some,, and for x denote by x that point on such that (x x ). For p h P 1 () we then set (6.34) p h := p h (x ). Fig. 6.3: Level lines of the extension p h Further, for h (Ω) h (Γ N ) we define D as the union of elements T h containing as a common edge (6.35) D := { T h h () }.

10 10 Ronald H.W. Hoppe Fig. 6.4: The set D We have the following properties of the extensions: Lemma 6.2 Basic properties of the bubble functions. Part II There exist constants C ν > 0, 13 ν 14, depending only on the shape regularity of the triangulations T h such that (6.36) (6.37) p h ψ 1,D p h ψ 0,D C 13 h 1/2 p h 0,e, p h P 1 (), C 14 h 1/2 p h 0,, p h P 1 (). Further, there exists a constant C 15 > 0 independent of h, h such that for all v V and µ = 0, 1: (6.38) ( h 1 µ ) 1/2 C 15 ( h 1 µ v 2 µ,) 1/2. T h h (Ω) h (Γ N ) v 2 µ,d We are now able to prove that up to higher order terms the estimator η R also does provide a lower bound for the total error e: Theorem 6.3 Lower bound for the total error There exist constants γ R, γ > 0, depending only on ᾱ and on the shape regularity of T h such that (6.39) γ R η R γ osc e 1,Ω, where η R and osc are given as in the previous theorem. The theorem can be proved by a series of results which establish upper bounds for the local contributions η (ν) R,T, ν 3 of the estimator η R. Lemma 6.3 Upper bounds for the local contributions (i) Let T h. Then there holds: (6.40) h π h f L ũ h 0, ᾱ C 8 C 10 e 1, + C 8 C 11 h f π h f 0,.

11 (ii) Let h (Ω). Then there holds: Finite lement Methods 11 (6.41) h 1/2 [n a grad ũ h ] J 0, ᾱc 9 C 13 e 1,D + +C 9 C 14 h f π h f 0,D + C 9 C 14 h π h f Lũ h 0,D. (iii) Let h (Γ N ). Then there holds: (6.42) h 1/2 π h g n grad ũ h 0, ᾱc 9 C 13 e 1,D +C 9 C 14 h f π h f 0,D + C 9 C 12 h 1/2 g π h g 0, + + C 9 C 14 h π h f Lũ h 0,D. Proof of (i) For the proof of (6.40) we set p h := π h f. Observing ψ = 0, by Green s formula (6.43) a (ũ h, p h ψ ) = + n a grad ũ h p h ψ }{{} = 0 div (a grad ũ h ) p h ψ dx + dσ. Then, using (6.30),(6.32) and (6.33) and taking advantage of (6.4),(6.43), it follows that π h f L ũ h 2 0, C 8 (π h f L ũ h ) π h ψ dx = = C 8 ( f π h ψ dx a (ũ h, π h ψ ) + + ) (π h f f) π h ψ dx = = C 8 (a (e, π h ψ ) + ) (π h f f) π h ψ dx C 8 C 10 α h 1 e 1, p h 0, + C 8 C 11 π h f f 0, p h 0,, from which (6.40) can be easily deduced.

12 12 Ronald H.W. Hoppe Proof of (ii) We set p h := [n a grad ũ h ] J. Again, in view of ψ = 0,,, Green s formula gives (6.44) D n D a grad ũ h p h ψ dσ = = a D (ũ h, p h ψ ) + D div a grad ũ h }{{} = Lũ h p h ψ dx. If we use (6.31),(6.37) and (6.38) and observe (6.4) and (6.44), it follows that [n D C 9 = D a grad ũ h ] J 2 0, [n D [n D a grad ũ h ] J p h ψ dσ = a grad ũ h ] J p h ψ dσ = = C 9 (a D (ũ h, p h ψ ) + (f π h f) p h ψ dx + D f p h ψ dx + ) (π h f Lũ h ) p h ψ dx = D D = C 9 a D (e, p h ψ ) + ( + C 9 (f π h f) p h ψ dx + ) (π h f Lũ h ) p h ψ dx D D C 9 C 13 α h 1/2 e 1,D p h 0, + + C 9 C 14 h 1/2 f π hf 0,D p h 0, + + C 9 C 14 h 1/2 π h Lũ h 0,D p h 0,, from which we readily deduce (6.41).

13 Finite lement Methods 13 Proof of (iii) We set p h := π hg n a grad ũ h. Observing ψ = 0,, by Green s formula we obtain (6.45) n a grad ũ h p h ψ dσ = = D n D = a D (ũ h, p h ψ ) + a grad ũ h p h ψ dσ = D div a grad ũ h }{{} = Lũ h p h ψ dx. Now, using (6.31),(6.34),(6.37) and (6.38) and taking advantage of (6.4),(6.45), we get π h g n a grad ũ h 2 0, C 9 (π h g n a grad ũ h ) p h ψ dσ = = C 9 ( + D f p h ψ dx + (π h f f) p h ψ dx + g p h ψ dσ a D (ũ h, p h ψ ) + (π h g g) f p h ψ dσ D ) (π h f Lũ h ) p h ψ dx = D = C 9 a D (e, p h + ψ ) + C 9 ( (π h g g) f p h ψ dσ D (π h f f) p h ψ dx + ) (π h f Lũ h ) p h ψ dx D C 9 C 13 α h 1/2 e 1,D p h 0, + + C 9 C 14 h 1/2 π hf f 0,D p h 0, + + C 9 C 12 π h g g 0, p h 0, + + C 9 C 14 h 1/2 π hf Lũ h 0,D p h 0,,

14 14 Ronald H.W. Hoppe from which (6.41) follows easily. References [1] Ainsworth, M. and Oden, J.T. (2000). A Posteriori rror stimation in Finite lement Analysis. Wiley, Chichester. [2] Babuska, I. and Strouboulis, T. (2001). The Finite lement Method and its Reliability. Clarendon Press, Oxford. [3] Bangerth, W. and Rannacher, R. (2003). Adaptive Finite lement Methods for Differential quations. Lectures in Mathematics. TH-Zürich. Birkhäuser, Basel. [4] riksson,. and step, D. and Hansbo, P. and Johnson, C. (1995). Computational Differential quations. Cambridge University Press, Cambridge. [5] Verfürth, R. (1996). A Review of A Posteriori stimation and Adaptive Mesh- Refinement Techniques. Wiley-Teubner, New York, Stuttgart.

Basic Concepts of Adaptive Finite Element Methods for Elliptic Boundary Value Problems

Basic Concepts of Adaptive Finite Element Methods for Elliptic Boundary Value Problems Basic Concepts of Adaptive Finite lement Methods for lliptic Boundary Value Problems Ronald H.W. Hoppe 1,2 1 Department of Mathematics, University of Houston 2 Institute of Mathematics, University of Augsburg

More information

Numerical Solutions to Partial Differential Equations

Numerical Solutions to Partial Differential Equations Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University The Residual and Error of Finite Element Solutions Mixed BVP of Poisson Equation

More information

Chapter 3 Conforming Finite Element Methods 3.1 Foundations Ritz-Galerkin Method

Chapter 3 Conforming Finite Element Methods 3.1 Foundations Ritz-Galerkin Method Chapter 3 Conforming Finite Element Methods 3.1 Foundations 3.1.1 Ritz-Galerkin Method Let V be a Hilbert space, a(, ) : V V lr a bounded, V-elliptic bilinear form and l : V lr a bounded linear functional.

More information

A Posteriori Estimates for Cost Functionals of Optimal Control Problems

A Posteriori Estimates for Cost Functionals of Optimal Control Problems A Posteriori Estimates for Cost Functionals of Optimal Control Problems Alexandra Gaevskaya, Ronald H.W. Hoppe,2 and Sergey Repin 3 Institute of Mathematics, Universität Augsburg, D-8659 Augsburg, Germany

More information

Energy norm a-posteriori error estimation for divergence-free discontinuous Galerkin approximations of the Navier-Stokes equations

Energy norm a-posteriori error estimation for divergence-free discontinuous Galerkin approximations of the Navier-Stokes equations INTRNATIONAL JOURNAL FOR NUMRICAL MTHODS IN FLUIDS Int. J. Numer. Meth. Fluids 19007; 1:1 [Version: 00/09/18 v1.01] nergy norm a-posteriori error estimation for divergence-free discontinuous Galerkin approximations

More information

Institut für Mathematik

Institut für Mathematik U n i v e r s i t ä t A u g s b u r g Institut für Mathematik Dietrich Braess, Carsten Carstensen, and Ronald H.W. Hoppe Convergence Analysis of a Conforming Adaptive Finite Element Method for an Obstacle

More information

Chapter 5 A priori error estimates for nonconforming finite element approximations 5.1 Strang s first lemma

Chapter 5 A priori error estimates for nonconforming finite element approximations 5.1 Strang s first lemma Chapter 5 A priori error estimates for nonconforming finite element approximations 51 Strang s first lemma We consider the variational equation (51 a(u, v = l(v, v V H 1 (Ω, and assume that the conditions

More information

A posteriori error estimation for elliptic problems

A posteriori error estimation for elliptic problems A posteriori error estimation for elliptic problems Praveen. C praveen@math.tifrbng.res.in Tata Institute of Fundamental Research Center for Applicable Mathematics Bangalore 560065 http://math.tifrbng.res.in

More information

ETNA Kent State University

ETNA Kent State University Electronic Transactions on Numerical Analysis. Volume 37, pp. 166-172, 2010. Copyright 2010,. ISSN 1068-9613. ETNA A GRADIENT RECOVERY OPERATOR BASED ON AN OBLIQUE PROJECTION BISHNU P. LAMICHHANE Abstract.

More information

A posteriori error estimates applied to flow in a channel with corners

A posteriori error estimates applied to flow in a channel with corners Mathematics and Computers in Simulation 61 (2003) 375 383 A posteriori error estimates applied to flow in a channel with corners Pavel Burda a,, Jaroslav Novotný b, Bedřich Sousedík a a Department of Mathematics,

More information

Downloaded 07/25/13 to Redistribution subject to SIAM license or copyright; see

Downloaded 07/25/13 to Redistribution subject to SIAM license or copyright; see SIAM J. NUMR. ANAL. Vol. 39, No. 6, pp. 2034 2044 c 2002 Society for Industrial and Applied Mathematics RSIDUAL-BASD A POSRIORI RROR SIMA FOR A NONCONFORMING RISSNR MINDLIN PLA FINI LMN CARSN CARSNSN Abstract.

More information

ENERGY NORM A POSTERIORI ERROR ESTIMATES FOR MIXED FINITE ELEMENT METHODS

ENERGY NORM A POSTERIORI ERROR ESTIMATES FOR MIXED FINITE ELEMENT METHODS ENERGY NORM A POSTERIORI ERROR ESTIMATES FOR MIXED FINITE ELEMENT METHODS CARLO LOVADINA AND ROLF STENBERG Abstract The paper deals with the a-posteriori error analysis of mixed finite element methods

More information

A posteriori error estimation in the FEM

A posteriori error estimation in the FEM A posteriori error estimation in the FEM Plan 1. Introduction 2. Goal-oriented error estimates 3. Residual error estimates 3.1 Explicit 3.2 Subdomain error estimate 3.3 Self-equilibrated residuals 3.4

More information

Finite Elements. Colin Cotter. February 22, Colin Cotter FEM

Finite Elements. Colin Cotter. February 22, Colin Cotter FEM Finite Elements February 22, 2019 In the previous sections, we introduced the concept of finite element spaces, which contain certain functions defined on a domain. Finite element spaces are examples of

More information

Numerical Solutions to Partial Differential Equations

Numerical Solutions to Partial Differential Equations Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Nonconformity and the Consistency Error First Strang Lemma Abstract Error Estimate

More information

R T (u H )v + (2.1) J S (u H )v v V, T (2.2) (2.3) H S J S (u H ) 2 L 2 (S). S T

R T (u H )v + (2.1) J S (u H )v v V, T (2.2) (2.3) H S J S (u H ) 2 L 2 (S). S T 2 R.H. NOCHETTO 2. Lecture 2. Adaptivity I: Design and Convergence of AFEM tarting with a conforming mesh T H, the adaptive procedure AFEM consists of loops of the form OLVE ETIMATE MARK REFINE to produce

More information

Adaptive Finite Element Methods Lecture 1: A Posteriori Error Estimation

Adaptive Finite Element Methods Lecture 1: A Posteriori Error Estimation Adaptive Finite Element Methods Lecture 1: A Posteriori Error Estimation Department of Mathematics and Institute for Physical Science and Technology University of Maryland, USA www.math.umd.edu/ rhn 7th

More information

A GAGLIARDO NIRENBERG INEQUALITY, WITH APPLICATION TO DUALITY-BASED A POSTERIORI ESTIMATION IN THE L 1 NORM

A GAGLIARDO NIRENBERG INEQUALITY, WITH APPLICATION TO DUALITY-BASED A POSTERIORI ESTIMATION IN THE L 1 NORM 27 Kragujevac J. Math. 30 (2007 27 43. A GAGLIARDO NIRENBERG INEQUALITY, WITH APPLICATION TO DUALITY-BASED A POSTERIORI ESTIMATION IN THE L 1 NORM Endre Süli Oxford University Computing Laboratory, Wolfson

More information

Basic Concepts of Adaptive Finite Element Methods for Elliptic Boundary Value Problems

Basic Concepts of Adaptive Finite Element Methods for Elliptic Boundary Value Problems Department of Mathematics, University of Houston Basic Concepts of Adaptive Finite Element Methods for Elliptic Boundary Value Problems Ronald H.W. Hoppe 1,2 1 Department of Mathematics, University of

More information

On an Approximation Result for Piecewise Polynomial Functions. O. Karakashian

On an Approximation Result for Piecewise Polynomial Functions. O. Karakashian BULLETIN OF THE GREE MATHEMATICAL SOCIETY Volume 57, 010 (1 7) On an Approximation Result for Piecewise Polynomial Functions O. arakashian Abstract We provide a new approach for proving approximation results

More information

ELLIPTIC RECONSTRUCTION AND A POSTERIORI ERROR ESTIMATES FOR PARABOLIC PROBLEMS

ELLIPTIC RECONSTRUCTION AND A POSTERIORI ERROR ESTIMATES FOR PARABOLIC PROBLEMS ELLIPTIC RECONSTRUCTION AND A POSTERIORI ERROR ESTIMATES FOR PARABOLIC PROBLEMS CHARALAMBOS MAKRIDAKIS AND RICARDO H. NOCHETTO Abstract. It is known that the energy technique for a posteriori error analysis

More information

LECTURE 1: SOURCES OF ERRORS MATHEMATICAL TOOLS A PRIORI ERROR ESTIMATES. Sergey Korotov,

LECTURE 1: SOURCES OF ERRORS MATHEMATICAL TOOLS A PRIORI ERROR ESTIMATES. Sergey Korotov, LECTURE 1: SOURCES OF ERRORS MATHEMATICAL TOOLS A PRIORI ERROR ESTIMATES Sergey Korotov, Institute of Mathematics Helsinki University of Technology, Finland Academy of Finland 1 Main Problem in Mathematical

More information

arxiv: v2 [math.na] 23 Apr 2016

arxiv: v2 [math.na] 23 Apr 2016 Improved ZZ A Posteriori Error Estimators for Diffusion Problems: Conforming Linear Elements arxiv:508.009v2 [math.na] 23 Apr 206 Zhiqiang Cai Cuiyu He Shun Zhang May 2, 208 Abstract. In [8], we introduced

More information

An a posteriori error estimate and a Comparison Theorem for the nonconforming P 1 element

An a posteriori error estimate and a Comparison Theorem for the nonconforming P 1 element Calcolo manuscript No. (will be inserted by the editor) An a posteriori error estimate and a Comparison Theorem for the nonconforming P 1 element Dietrich Braess Faculty of Mathematics, Ruhr-University

More information

A posteriori error estimation of approximate boundary fluxes

A posteriori error estimation of approximate boundary fluxes COMMUNICATIONS IN NUMERICAL METHODS IN ENGINEERING Commun. Numer. Meth. Engng 2008; 24:421 434 Published online 24 May 2007 in Wiley InterScience (www.interscience.wiley.com)..1014 A posteriori error estimation

More information

A POSTERIORI ERROR ESTIMATION FOR NON-CONFORMING QUADRILATERAL FINITE ELEMENTS

A POSTERIORI ERROR ESTIMATION FOR NON-CONFORMING QUADRILATERAL FINITE ELEMENTS INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume 2, Number, Pages 8 c 2005 Institute for Scientific Computing and Information A POSTERIORI ERROR ESTIMATION FOR NON-CONFORMING QUADRILATERAL

More information

Chapter 1 Mathematical Foundations

Chapter 1 Mathematical Foundations Computational Electromagnetics; Chapter 1 1 Chapter 1 Mathematical Foundations 1.1 Maxwell s Equations Electromagnetic phenomena can be described by the electric field E, the electric induction D, the

More information

INTRODUCTION TO FINITE ELEMENT METHODS

INTRODUCTION TO FINITE ELEMENT METHODS INTRODUCTION TO FINITE ELEMENT METHODS LONG CHEN Finite element methods are based on the variational formulation of partial differential equations which only need to compute the gradient of a function.

More information

Find (u,p;λ), with u 0 and λ R, such that u + p = λu in Ω, (2.1) div u = 0 in Ω, u = 0 on Γ.

Find (u,p;λ), with u 0 and λ R, such that u + p = λu in Ω, (2.1) div u = 0 in Ω, u = 0 on Γ. A POSTERIORI ESTIMATES FOR THE STOKES EIGENVALUE PROBLEM CARLO LOVADINA, MIKKO LYLY, AND ROLF STENBERG Abstract. We consider the Stokes eigenvalue problem. For the eigenvalues we derive both upper and

More information

ENERGY NORM A POSTERIORI ERROR ESTIMATES FOR MIXED FINITE ELEMENT METHODS

ENERGY NORM A POSTERIORI ERROR ESTIMATES FOR MIXED FINITE ELEMENT METHODS MATHEMATICS OF COMPUTATION Volume 75, Number 256, October 2006, Pages 1659 1674 S 0025-57180601872-2 Article electronically published on June 26, 2006 ENERGY NORM A POSTERIORI ERROR ESTIMATES FOR MIXED

More information

c 2007 Society for Industrial and Applied Mathematics

c 2007 Society for Industrial and Applied Mathematics SIAM J. NUMR. ANAL. Vol. 45, No. 1, pp. 68 82 c 2007 Society for Industrial and Applied Mathematics FRAMWORK FOR TH A POSTRIORI RROR ANALYSIS OF NONCONFORMING FINIT LMNTS CARSTN CARSTNSN, JUN HU, AND ANTONIO

More information

Adaptive Finite Element Methods Lecture Notes Winter Term 2017/18. R. Verfürth. Fakultät für Mathematik, Ruhr-Universität Bochum

Adaptive Finite Element Methods Lecture Notes Winter Term 2017/18. R. Verfürth. Fakultät für Mathematik, Ruhr-Universität Bochum Adaptive Finite Element Methods Lecture Notes Winter Term 2017/18 R. Verfürth Fakultät für Mathematik, Ruhr-Universität Bochum Contents Chapter I. Introduction 7 I.1. Motivation 7 I.2. Sobolev and finite

More information

An interpolation operator for H 1 functions on general quadrilateral and hexahedral meshes with hanging nodes

An interpolation operator for H 1 functions on general quadrilateral and hexahedral meshes with hanging nodes An interpolation operator for H 1 functions on general quadrilateral and hexahedral meshes with hanging nodes Vincent Heuveline Friedhelm Schieweck Abstract We propose a Scott-Zhang type interpolation

More information

Chapter 1 Foundations of Elliptic Boundary Value Problems 1.1 Euler equations of variational problems

Chapter 1 Foundations of Elliptic Boundary Value Problems 1.1 Euler equations of variational problems Chapter 1 Foundations of Elliptic Boundary Value Problems 1.1 Euler equations of variational problems Elliptic boundary value problems often occur as the Euler equations of variational problems the latter

More information

A POSTERIORI ERROR ESTIMATION OF FINITE ELEMENT APPROXIMATIONS OF POINTWISE STATE CONSTRAINED DISTRIBUTED CONTROL PROBLEMS

A POSTERIORI ERROR ESTIMATION OF FINITE ELEMENT APPROXIMATIONS OF POINTWISE STATE CONSTRAINED DISTRIBUTED CONTROL PROBLEMS A POSTERIORI ERROR ESTIMATION OF FINITE ELEMENT APPROXIMATIONS OF POINTWISE STATE CONSTRAINED DISTRIBUTED CONTROL PROBLEMS RONALD H.W. HOPPE AND MICHAEL KIEWEG Abstract. We provide an a posteriori error

More information

Axioms of Adaptivity (AoA) in Lecture 3 (sufficient for optimal convergence rates)

Axioms of Adaptivity (AoA) in Lecture 3 (sufficient for optimal convergence rates) Axioms of Adaptivity (AoA) in Lecture 3 (sufficient for optimal convergence rates) Carsten Carstensen Humboldt-Universität zu Berlin 2018 International Graduate Summer School on Frontiers of Applied and

More information

Enhancing eigenvalue approximation by gradient recovery on adaptive meshes

Enhancing eigenvalue approximation by gradient recovery on adaptive meshes IMA Journal of Numerical Analysis Advance Access published October 29, 2008 IMA Journal of Numerical Analysis Page 1 of 15 doi:10.1093/imanum/drn050 Enhancing eigenvalue approximation by gradient recovery

More information

A Posteriori Error Estimation Techniques for Finite Element Methods. Zhiqiang Cai Purdue University

A Posteriori Error Estimation Techniques for Finite Element Methods. Zhiqiang Cai Purdue University A Posteriori Error Estimation Techniques for Finite Element Methods Zhiqiang Cai Purdue University Department of Mathematics, Purdue University Slide 1, March 16, 2017 Books Ainsworth & Oden, A posteriori

More information

PREPRINT 2010:25. Fictitious domain finite element methods using cut elements: II. A stabilized Nitsche method ERIK BURMAN PETER HANSBO

PREPRINT 2010:25. Fictitious domain finite element methods using cut elements: II. A stabilized Nitsche method ERIK BURMAN PETER HANSBO PREPRINT 2010:25 Fictitious domain finite element methods using cut elements: II. A stabilized Nitsche method ERIK BURMAN PETER HANSBO Department of Mathematical Sciences Division of Mathematics CHALMERS

More information

1. Introduction. We consider the following convection-diffusion-reaction problem: Find u = u(x) such that

1. Introduction. We consider the following convection-diffusion-reaction problem: Find u = u(x) such that SIAM J. UMR. AAL. Vol. 42, o. 6, pp. 2496 2521 c 2005 Society for Industrial and Applied Mathematics XPLICIT AD AVRAGIG A POSTRIORI RROR STIMATS FOR ADAPTIV FIIT VOLUM MTHODS C. CARSTS, R. LAZAROV, AD

More information

A robust a-posteriori error estimate for hp-adaptive DG methods for convection-diffusion equations

A robust a-posteriori error estimate for hp-adaptive DG methods for convection-diffusion equations IMA Journal of Numerical Analysis Page 1 of 3 doi:10.1093/imanum/drnxxx A robust a-posteriori error estimate for hp-adaptive DG methods for convection-diffusion equations LIANG ZHU AND DOMINIK SCHÖTZAU

More information

Adaptive Finite Element Methods for Elliptic Optimal Control Problems

Adaptive Finite Element Methods for Elliptic Optimal Control Problems Adaptive Finite Element Methods for Elliptic Optimal Control Problems Ronald H.W. Hoppe 1,2 1 2 Institute of Mathematics, University of Augsburg Summer School Optimal Control of PDEs Cortona, July 12-17,

More information

A mixed finite element approximation of the Stokes equations with the boundary condition of type (D+N)

A mixed finite element approximation of the Stokes equations with the boundary condition of type (D+N) wwwijmercom Vol2, Issue1, Jan-Feb 2012 pp-464-472 ISSN: 2249-6645 A mixed finite element approximation of the Stokes equations with the boundary condition of type (D+N) Jaouad El-Mekkaoui 1, Abdeslam Elakkad

More information

ADAPTIVE HYBRIDIZED INTERIOR PENALTY DISCONTINUOUS GALERKIN METHODS FOR H(CURL)-ELLIPTIC PROBLEMS

ADAPTIVE HYBRIDIZED INTERIOR PENALTY DISCONTINUOUS GALERKIN METHODS FOR H(CURL)-ELLIPTIC PROBLEMS ADAPTIVE HYBRIDIZED INTERIOR PENALTY DISCONTINUOUS GALERKIN METHODS OR HCURL-ELLIPTIC PROBLEMS C. CARSTENSEN, R. H. W. HOPPE, N. SHARMA, AND T. WARBURTON Abstract. We develop and analyze an adaptive hybridized

More information

Graded Mesh Refinement and Error Estimates of Higher Order for DGFE-solutions of Elliptic Boundary Value Problems in Polygons

Graded Mesh Refinement and Error Estimates of Higher Order for DGFE-solutions of Elliptic Boundary Value Problems in Polygons Graded Mesh Refinement and Error Estimates of Higher Order for DGFE-solutions of Elliptic Boundary Value Problems in Polygons Anna-Margarete Sändig, Miloslav Feistauer University Stuttgart, IANS Journées

More information

ASYMPTOTICALLY EXACT A POSTERIORI ESTIMATORS FOR THE POINTWISE GRADIENT ERROR ON EACH ELEMENT IN IRREGULAR MESHES. PART II: THE PIECEWISE LINEAR CASE

ASYMPTOTICALLY EXACT A POSTERIORI ESTIMATORS FOR THE POINTWISE GRADIENT ERROR ON EACH ELEMENT IN IRREGULAR MESHES. PART II: THE PIECEWISE LINEAR CASE MATEMATICS OF COMPUTATION Volume 73, Number 246, Pages 517 523 S 0025-5718(0301570-9 Article electronically published on June 17, 2003 ASYMPTOTICALLY EXACT A POSTERIORI ESTIMATORS FOR TE POINTWISE GRADIENT

More information

A POSTERIORI FINITE ELEMENT ERROR ESTIMATION FOR SECOND-ORDER HYPERBOLIC PROBLEMS

A POSTERIORI FINITE ELEMENT ERROR ESTIMATION FOR SECOND-ORDER HYPERBOLIC PROBLEMS A POSTERIORI FINITE ELEMENT ERROR ESTIMATION FOR SECOND-ORDER HYPERBOLIC PROBLEMS Slimane Adjerid Department of Mathematics and Interdisciplinary Center for Applied Mathematics Virginia Polytechnic Institute

More information

Lecture Note III: Least-Squares Method

Lecture Note III: Least-Squares Method Lecture Note III: Least-Squares Method Zhiqiang Cai October 4, 004 In this chapter, we shall present least-squares methods for second-order scalar partial differential equations, elastic equations of solids,

More information

Geometric Multigrid Methods

Geometric Multigrid Methods Geometric Multigrid Methods Susanne C. Brenner Department of Mathematics and Center for Computation & Technology Louisiana State University IMA Tutorial: Fast Solution Techniques November 28, 2010 Ideas

More information

AN EQUILIBRATED A POSTERIORI ERROR ESTIMATOR FOR THE INTERIOR PENALTY DISCONTINUOUS GALERKIN METHOD

AN EQUILIBRATED A POSTERIORI ERROR ESTIMATOR FOR THE INTERIOR PENALTY DISCONTINUOUS GALERKIN METHOD AN EQUILIBRATED A POSTERIORI ERROR ESTIMATOR FOR THE INTERIOR PENALTY DISCONTINUOUS GALERIN METHOD D. BRAESS, T. FRAUNHOLZ, AND R. H. W. HOPPE Abstract. Interior Penalty Discontinuous Galerkin (IPDG) methods

More information

Chapter 1: The Finite Element Method

Chapter 1: The Finite Element Method Chapter 1: The Finite Element Method Michael Hanke Read: Strang, p 428 436 A Model Problem Mathematical Models, Analysis and Simulation, Part Applications: u = fx), < x < 1 u) = u1) = D) axial deformation

More information

Analysis of an Adaptive Finite Element Method for Recovering the Robin Coefficient

Analysis of an Adaptive Finite Element Method for Recovering the Robin Coefficient Analysis of an Adaptive Finite Element Method for Recovering the Robin Coefficient Yifeng Xu 1 Jun Zou 2 Abstract Based on a new a posteriori error estimator, an adaptive finite element method is proposed

More information

A Posteriori Error Estimation for Highly Indefinite Helmholtz Problems

A Posteriori Error Estimation for Highly Indefinite Helmholtz Problems Computational Methods in Applied Mathematics Vol. 13 013), No. 3, pp. 333 347 c 013 Institute of Mathematics, NAS of Belarus Doi: 10.1515/cmam-013-0008 A Posteriori Error Estimation for Highly Indefinite

More information

PREPRINT 2010:23. A nonconforming rotated Q 1 approximation on tetrahedra PETER HANSBO

PREPRINT 2010:23. A nonconforming rotated Q 1 approximation on tetrahedra PETER HANSBO PREPRINT 2010:23 A nonconforming rotated Q 1 approximation on tetrahedra PETER HANSBO Department of Mathematical Sciences Division of Mathematics CHALMERS UNIVERSITY OF TECHNOLOGY UNIVERSITY OF GOTHENBURG

More information

arxiv: v1 [math.na] 17 Nov 2017 Received: date / Accepted: date

arxiv: v1 [math.na] 17 Nov 2017 Received: date / Accepted: date Noname manuscript No. (will be inserted by the editor) Maximum norm a posteriori error estimates for an optimal control problem Alejandro Allendes Enrique Otárola Richard Rankin Abner J. algado arxiv:1711.06707v1

More information

A POSTERIORI ERROR ESTIMATES FOR FOURTH-ORDER ELLIPTIC PROBLEMS

A POSTERIORI ERROR ESTIMATES FOR FOURTH-ORDER ELLIPTIC PROBLEMS A POSTERIORI ERROR ESTIMATES FOR FOURTH-ORDER ELLIPTIC PROBLEMS Slimane Adjerid Department of Mathematics and Interdisciplinary Center for Applied Mathematics Virginia Polytechnic Institute and State University

More information

Scientific Computing WS 2017/2018. Lecture 18. Jürgen Fuhrmann Lecture 18 Slide 1

Scientific Computing WS 2017/2018. Lecture 18. Jürgen Fuhrmann Lecture 18 Slide 1 Scientific Computing WS 2017/2018 Lecture 18 Jürgen Fuhrmann juergen.fuhrmann@wias-berlin.de Lecture 18 Slide 1 Lecture 18 Slide 2 Weak formulation of homogeneous Dirichlet problem Search u H0 1 (Ω) (here,

More information

Complementarity based a posteriori error estimates and their properties

Complementarity based a posteriori error estimates and their properties Complementarity based a posteriori error estimates and their properties Tomáš Vejchodský November 3, 2009 Institute of Mathematics, Czech Academy of Sciences Žitná 25, CZ-115 67 Praha 1, Czech Republic

More information

Overview. A Posteriori Error Estimates for the Biharmonic Equation. Variational Formulation and Discretization. The Biharmonic Equation

Overview. A Posteriori Error Estimates for the Biharmonic Equation. Variational Formulation and Discretization. The Biharmonic Equation Overview A Posteriori rror stimates for the Biharmonic quation R Verfürth Fakultät für Mathematik Ruhr-Universität Bochum wwwruhr-uni-bochumde/num1 Milan / February 11th, 013 The Biharmonic quation Summary

More information

A Least-Squares Finite Element Approximation for the Compressible Stokes Equations

A Least-Squares Finite Element Approximation for the Compressible Stokes Equations A Least-Squares Finite Element Approximation for the Compressible Stokes Equations Zhiqiang Cai, 1 Xiu Ye 1 Department of Mathematics, Purdue University, 1395 Mathematical Science Building, West Lafayette,

More information

ALL FIRST-ORDER AVERAGING TECHNIQUES FOR A POSTERIORI FINITE ELEMENT ERROR CONTROL ON UNSTRUCTURED GRIDS ARE EFFICIENT AND RELIABLE

ALL FIRST-ORDER AVERAGING TECHNIQUES FOR A POSTERIORI FINITE ELEMENT ERROR CONTROL ON UNSTRUCTURED GRIDS ARE EFFICIENT AND RELIABLE MATHEMATICS OF COMPUTATION Volume 73, Number 247, Pages 1153 1165 S 0025-5718(03)01580-1 Article electronically published on August 12, 2003 ALL FIRST-ORDER AVERAGING TECHNIQUES FOR A POSTERIORI FINITE

More information

MULTIGRID PRECONDITIONING IN H(div) ON NON-CONVEX POLYGONS* Dedicated to Professor Jim Douglas, Jr. on the occasion of his seventieth birthday.

MULTIGRID PRECONDITIONING IN H(div) ON NON-CONVEX POLYGONS* Dedicated to Professor Jim Douglas, Jr. on the occasion of his seventieth birthday. MULTIGRID PRECONDITIONING IN H(div) ON NON-CONVEX POLYGONS* DOUGLAS N ARNOLD, RICHARD S FALK, and RAGNAR WINTHER Dedicated to Professor Jim Douglas, Jr on the occasion of his seventieth birthday Abstract

More information

Applied Numerical Mathematics

Applied Numerical Mathematics Applied Numerical Mathematics 59 009) 119 130 Contents lists available at ScienceDirect Applied Numerical Mathematics www.elsevier.com/locate/apnum Convergence of adaptive finite element methods in computational

More information

Adaptive methods for control problems with finite-dimensional control space

Adaptive methods for control problems with finite-dimensional control space Adaptive methods for control problems with finite-dimensional control space Saheed Akindeinde and Daniel Wachsmuth Johann Radon Institute for Computational and Applied Mathematics (RICAM) Austrian Academy

More information

Supraconvergence of a Non-Uniform Discretisation for an Elliptic Third-Kind Boundary-Value Problem with Mixed Derivatives

Supraconvergence of a Non-Uniform Discretisation for an Elliptic Third-Kind Boundary-Value Problem with Mixed Derivatives Supraconvergence of a Non-Uniform Discretisation for an Elliptic Third-Kind Boundary-Value Problem with Mixed Derivatives Etienne Emmrich Technische Universität Berlin, Institut für Mathematik, Straße

More information

A posteriori estimators for obstacle problems by the hypercircle method

A posteriori estimators for obstacle problems by the hypercircle method A posteriori estimators for obstacle problems by the hypercircle method Dietrich Braess 1 Ronald H.W. Hoppe 2,3 Joachim Schöberl 4 January 9, 2008 Abstract A posteriori error estimates for the obstacle

More information

Error estimates for the Raviart-Thomas interpolation under the maximum angle condition

Error estimates for the Raviart-Thomas interpolation under the maximum angle condition Error estimates for the Raviart-Thomas interpolation under the maximum angle condition Ricardo G. Durán and Ariel L. Lombardi Abstract. The classical error analysis for the Raviart-Thomas interpolation

More information

Numerical Methods for Large-Scale Nonlinear Systems

Numerical Methods for Large-Scale Nonlinear Systems Numerical Methods for Large-Scale Nonlinear Systems Handouts by Ronald H.W. Hoppe following the monograph P. Deuflhard Newton Methods for Nonlinear Problems Springer, Berlin-Heidelberg-New York, 2004 Num.

More information

A UNIFYING THEORY OF A POSTERIORI FINITE ELEMENT ERROR CONTROL

A UNIFYING THEORY OF A POSTERIORI FINITE ELEMENT ERROR CONTROL A UNIFYING THEORY OF A POSTERIORI FINITE ELEMENT ERROR CONTROL C. CARSTENSEN Abstract. Residual-based a posteriori error estimates are derived within a unified setting for lowest-order conforming, nonconforming,

More information

Some New Elements for the Reissner Mindlin Plate Model

Some New Elements for the Reissner Mindlin Plate Model Boundary Value Problems for Partial Differential Equations and Applications, J.-L. Lions and C. Baiocchi, eds., Masson, 1993, pp. 287 292. Some New Elements for the Reissner Mindlin Plate Model Douglas

More information

Convergence and optimality of an adaptive FEM for controlling L 2 errors

Convergence and optimality of an adaptive FEM for controlling L 2 errors Convergence and optimality of an adaptive FEM for controlling L 2 errors Alan Demlow (University of Kentucky) joint work with Rob Stevenson (University of Amsterdam) Partially supported by NSF DMS-0713770.

More information

A posteriori error estimates for Maxwell Equations

A posteriori error estimates for Maxwell Equations www.oeaw.ac.at A posteriori error estimates for Maxwell Equations J. Schöberl RICAM-Report 2005-10 www.ricam.oeaw.ac.at A POSTERIORI ERROR ESTIMATES FOR MAXWELL EQUATIONS JOACHIM SCHÖBERL Abstract. Maxwell

More information

LOCAL MULTILEVEL METHODS FOR ADAPTIVE NONCONFORMING FINITE ELEMENT METHODS

LOCAL MULTILEVEL METHODS FOR ADAPTIVE NONCONFORMING FINITE ELEMENT METHODS LOCAL MULTILEVEL METHODS FOR ADAPTIVE NONCONFORMING FINITE ELEMENT METHODS XUEUN XU, HUANGXIN CHEN, AND RONALD H.W. HOPPE Abstract. In this paper, we propose a local multilevel product algorithm and its

More information

arxiv: v2 [math.na] 28 Feb 2018

arxiv: v2 [math.na] 28 Feb 2018 A POSTRIORI STIMATS FOR CONFORMING KIRCHHOFF PLAT LMNTS TOM GUSTAFSSON, ROLF STNBRG, AND JUHA VIDMAN arxiv:1707.08396v [math.na] 8 Feb 018 Abstract. We derive a residual a posteriori estimator for the

More information

Multigrid Methods for Elliptic Obstacle Problems on 2D Bisection Grids

Multigrid Methods for Elliptic Obstacle Problems on 2D Bisection Grids Multigrid Methods for Elliptic Obstacle Problems on 2D Bisection Grids Long Chen 1, Ricardo H. Nochetto 2, and Chen-Song Zhang 3 1 Department of Mathematics, University of California at Irvine. chenlong@math.uci.edu

More information

Institut für Mathematik

Institut für Mathematik U n i v e r s i t ä t A u g s b u r g Institut für Mathematik Xuejun Xu, Huangxin Chen, Ronald H.W. Hoppe Local Multilevel Methods for Adaptive Nonconforming Finite Element Methods Preprint Nr. 21/2009

More information

QUASI-OPTIMAL CONVERGENCE RATE OF AN ADAPTIVE DISCONTINUOUS GALERKIN METHOD

QUASI-OPTIMAL CONVERGENCE RATE OF AN ADAPTIVE DISCONTINUOUS GALERKIN METHOD QUASI-OPIMAL CONVERGENCE RAE OF AN ADAPIVE DISCONINUOUS GALERKIN MEHOD ANDREA BONIO AND RICARDO H. NOCHEO Abstract. We analyze an adaptive discontinuous finite element method (ADFEM) for symmetric second

More information

An A Posteriori Error Estimate for Discontinuous Galerkin Methods

An A Posteriori Error Estimate for Discontinuous Galerkin Methods An A Posteriori Error Estimate for Discontinuous Galerkin Methods Mats G Larson mgl@math.chalmers.se Chalmers Finite Element Center Mats G Larson Chalmers Finite Element Center p.1 Outline We present an

More information

Goal. Robust A Posteriori Error Estimates for Stabilized Finite Element Discretizations of Non-Stationary Convection-Diffusion Problems.

Goal. Robust A Posteriori Error Estimates for Stabilized Finite Element Discretizations of Non-Stationary Convection-Diffusion Problems. Robust A Posteriori Error Estimates for Stabilized Finite Element s of Non-Stationary Convection-Diffusion Problems L. Tobiska and R. Verfürth Universität Magdeburg Ruhr-Universität Bochum www.ruhr-uni-bochum.de/num

More information

A Mixed Nonconforming Finite Element for Linear Elasticity

A Mixed Nonconforming Finite Element for Linear Elasticity A Mixed Nonconforming Finite Element for Linear Elasticity Zhiqiang Cai, 1 Xiu Ye 2 1 Department of Mathematics, Purdue University, West Lafayette, Indiana 47907-1395 2 Department of Mathematics and Statistics,

More information

ABSTRACT CONVERGENCE OF ADAPTIVE FINITE ELEMENT METHODS. Khamron Mekchay, Doctor of Philosophy, 2005

ABSTRACT CONVERGENCE OF ADAPTIVE FINITE ELEMENT METHODS. Khamron Mekchay, Doctor of Philosophy, 2005 ABSTRACT Title of dissertation: CONVERGENCE OF ADAPTIVE FINITE ELEMENT METHODS Khamron Mekchay, Doctor of Philosophy, 2005 Dissertation directed by: Professor Ricardo H. Nochetto Department of Mathematics

More information

Numerical methods for PDEs FEM convergence, error estimates, piecewise polynomials

Numerical methods for PDEs FEM convergence, error estimates, piecewise polynomials Platzhalter für Bild, Bild auf Titelfolie hinter das Logo einsetzen Numerical methods for PDEs FEM convergence, error estimates, piecewise polynomials Dr. Noemi Friedman Contents of the course Fundamentals

More information

High order, finite volume method, flux conservation, finite element method

High order, finite volume method, flux conservation, finite element method FLUX-CONSERVING FINITE ELEMENT METHODS SHANGYOU ZHANG, ZHIMIN ZHANG, AND QINGSONG ZOU Abstract. We analyze the flux conservation property of the finite element method. It is shown that the finite element

More information

A Finite Element Method Using Singular Functions for Poisson Equations: Mixed Boundary Conditions

A Finite Element Method Using Singular Functions for Poisson Equations: Mixed Boundary Conditions A Finite Element Method Using Singular Functions for Poisson Equations: Mixed Boundary Conditions Zhiqiang Cai Seokchan Kim Sangdong Kim Sooryun Kong Abstract In [7], we proposed a new finite element method

More information

FEniCS Course. Lecture 8: A posteriori error estimates and adaptivity. Contributors André Massing Marie Rognes

FEniCS Course. Lecture 8: A posteriori error estimates and adaptivity. Contributors André Massing Marie Rognes FEniCS Course Lecture 8: A posteriori error estimates and adaptivity Contributors André Massing Marie Rognes 1 / 24 A priori estimates If u H k+1 (Ω) and V h = P k (T h ) then u u h Ch k u Ω,k+1 u u h

More information

Recovery-Based A Posteriori Error Estimation

Recovery-Based A Posteriori Error Estimation Recovery-Based A Posteriori Error Estimation Zhiqiang Cai Purdue University Department of Mathematics, Purdue University Slide 1, March 2, 2011 Outline Introduction Diffusion Problems Higher Order Elements

More information

PARTITION OF UNITY FOR THE STOKES PROBLEM ON NONMATCHING GRIDS

PARTITION OF UNITY FOR THE STOKES PROBLEM ON NONMATCHING GRIDS PARTITION OF UNITY FOR THE STOES PROBLEM ON NONMATCHING GRIDS CONSTANTIN BACUTA AND JINCHAO XU Abstract. We consider the Stokes Problem on a plane polygonal domain Ω R 2. We propose a finite element method

More information

Constrained H 1 -interpolation on quadrilateral and hexahedral meshes with hanging nodes

Constrained H 1 -interpolation on quadrilateral and hexahedral meshes with hanging nodes Constrained H 1 -interpolation on quadrilateral and hexahedral meshes with hanging nodes Vincent Heuveline Friedhelm Schieweck Abstract We propose a Scott-Zhang type finite element interpolation operator

More information

A posteriori error estimator based on gradient recovery by averaging for discontinuous Galerkin methods

A posteriori error estimator based on gradient recovery by averaging for discontinuous Galerkin methods A posteriori error estimator based on gradient recovery by averaging for discontinuous Galerkin methods Emmanuel Creusé, Serge Nicaise October 7, 29 Abstract We consider some (anisotropic and piecewise

More information

CONVERGENCE ANALYSIS OF AN ADAPTIVE INTERIOR PENALTY DISCONTINUOUS GALERKIN METHOD FOR THE HELMHOLTZ EQUATION

CONVERGENCE ANALYSIS OF AN ADAPTIVE INTERIOR PENALTY DISCONTINUOUS GALERKIN METHOD FOR THE HELMHOLTZ EQUATION CONVERGENCE ANALYSIS OF AN ADAPTIVE INTERIOR PENALTY DISCONTINUOUS GALERKIN METHOD FOR THE HELMHOLTZ EQUATION A Dissertation Presented to the Faculty of the Department of Mathematics University of Houston

More information

L 2 and pointwise a posteriori error estimates for FEM for elliptic PDEs on surfaces

L 2 and pointwise a posteriori error estimates for FEM for elliptic PDEs on surfaces IMA Journal of Numerical Analysis 214) Page 1 of 27 doi:1.193/imanum/drn L 2 and pointwise a posteriori error estimates for FEM for elliptic PDEs on surfaces Fernando Camacho and Alan Demlow, Department

More information

Institut für Mathematik

Institut für Mathematik U n i v e r s i t ä t A u g s b u r g Institut für Mathematik Ronald H.W. Hoppe, Irwin Yousept Adaptive Edge Element Approximation of H(curl-Elliptic Optimal Control Problems with Control Constraints Preprint

More information

arxiv: v1 [math.na] 19 Dec 2017

arxiv: v1 [math.na] 19 Dec 2017 Arnold-Winther Mixed Finite Elements for Stokes Eigenvalue Problems Joscha Gedicke Arbaz Khan arxiv:72.0686v [math.na] 9 Dec 207 Abstract This paper is devoted to study the Arnold-Winther mixed finite

More information

Key words. a posteriori error estimators, quasi-orthogonality, adaptive mesh refinement, error and oscillation reduction estimates, optimal meshes.

Key words. a posteriori error estimators, quasi-orthogonality, adaptive mesh refinement, error and oscillation reduction estimates, optimal meshes. CONVERGENCE OF ADAPTIVE FINITE ELEMENT METHOD FOR GENERAL ECOND ORDER LINEAR ELLIPTIC PDE KHAMRON MEKCHAY AND RICARDO H. NOCHETTO Abstract. We prove convergence of adaptive finite element methods (AFEM)

More information

arxiv: v3 [math.na] 8 Sep 2015

arxiv: v3 [math.na] 8 Sep 2015 A Recovery-Based A Posteriori Error Estimator for H(curl) Interface Problems arxiv:504.00898v3 [math.na] 8 Sep 205 Zhiqiang Cai Shuhao Cao Abstract This paper introduces a new recovery-based a posteriori

More information

1. Let a(x) > 0, and assume that u and u h are the solutions of the Dirichlet problem:

1. Let a(x) > 0, and assume that u and u h are the solutions of the Dirichlet problem: Mathematics Chalmers & GU TMA37/MMG800: Partial Differential Equations, 011 08 4; kl 8.30-13.30. Telephone: Ida Säfström: 0703-088304 Calculators, formula notes and other subject related material are not

More information

Finite Elements. Colin Cotter. January 15, Colin Cotter FEM

Finite Elements. Colin Cotter. January 15, Colin Cotter FEM Finite Elements January 15, 2018 Why Can solve PDEs on complicated domains. Have flexibility to increase order of accuracy and match the numerics to the physics. has an elegant mathematical formulation

More information

Numerische Mathematik

Numerische Mathematik Numer. Math. (2012) 122:61 99 DOI 10.1007/s00211-012-0456-x Numerische Mathematik C 0 elements for generalized indefinite Maxwell equations Huoyuan Duan Ping Lin Roger C. E. Tan Received: 31 July 2010

More information

Ultraconvergence of ZZ Patch Recovery at Mesh Symmetry Points

Ultraconvergence of ZZ Patch Recovery at Mesh Symmetry Points Ultraconvergence of ZZ Patch Recovery at Mesh Symmetry Points Zhimin Zhang and Runchang Lin Department of Mathematics, Wayne State University Abstract. The ultraconvergence property of the Zienkiewicz-Zhu

More information

Scientific Computing WS 2018/2019. Lecture 15. Jürgen Fuhrmann Lecture 15 Slide 1

Scientific Computing WS 2018/2019. Lecture 15. Jürgen Fuhrmann Lecture 15 Slide 1 Scientific Computing WS 2018/2019 Lecture 15 Jürgen Fuhrmann juergen.fuhrmann@wias-berlin.de Lecture 15 Slide 1 Lecture 15 Slide 2 Problems with strong formulation Writing the PDE with divergence and gradient

More information