Maximum-norm a posteriori estimates for discontinuous Galerkin methods
|
|
- Mitchell Fields
- 6 years ago
- Views:
Transcription
1 Maximum-norm a posteriori estimates for discontinuous Galerkin methods Emmanuil Georgoulis Department of Mathematics, University of Leicester, UK Based on joint work with Alan Demlow (Kentucky, USA) DG Workshop, Ηράκλειο, Κρήτη 1
2 Model Problem Find u H 1 (Ω) such that u = f in Ω u = g on Ω with: Ω R d, d = 2, 3, is a polyhedral domain (not necessarily Lipschitz!). f L (Ω). g C 0 ( Ω) with extension in Wq 1 (Ω) for some q > d. 2
3 Overview 1 Motivation 2 Review of max norm a posteriori estimates, strategies and difficulties 3 Reliable and efficient maximum norm a posteriori bounds for dg 4 Application of these ideas to standard (conforming) FEM 5 Numerical experiments 6 Conclusions 3
4 Motivation Energy norm a posteriori bounds for dg well understood. Becker, Hansbo & Larson ( 03), Karakashian & Pascal ( 03, 07), Houston, Schötzau & Wihler ( 07), Ainsworth ( 07), Houston, Süli & Wihler ( 08), Ern & Stephansen ( 08), Cochez-Dhondt & Nicaise ( 08), Carstensen, Gudi & Jensen ( 08), Vohralik ( 08), Creusé & Nicaise ( 10), Ainsworth & Rankin ( 10), Bonito & Nochetto ( 10), Ern, Stephansen & Vohralik ( 10)... 4
5 Motivation Energy norm a posteriori bounds for dg well understood. Becker, Hansbo & Larson ( 03), Karakashian & Pascal ( 03, 07), Houston, Schötzau & Wihler ( 07), Ainsworth ( 07), Houston, Süli & Wihler ( 08), Ern & Stephansen ( 08), Cochez-Dhondt & Nicaise ( 08), Carstensen, Gudi & Jensen ( 08), Vohralik ( 08), Creusé & Nicaise ( 10), Ainsworth & Rankin ( 10), Bonito & Nochetto ( 10), Ern, Stephansen & Vohralik ( 10)... Controlling max norm error is of natural interest in some applications, such as obstacle problems, contact problems, etc.,... French, Larson & Nochetto ( 01), Nochetto, Siebert & Veeser ( 03, 05), Nochetto, Schmidt, Siebert & Veeser ( 06)... 4
6 Motivation Energy norm a posteriori bounds for dg well understood. Becker, Hansbo & Larson ( 03), Karakashian & Pascal ( 03, 07), Houston, Schötzau & Wihler ( 07), Ainsworth ( 07), Houston, Süli & Wihler ( 08), Ern & Stephansen ( 08), Cochez-Dhondt & Nicaise ( 08), Carstensen, Gudi & Jensen ( 08), Vohralik ( 08), Creusé & Nicaise ( 10), Ainsworth & Rankin ( 10), Bonito & Nochetto ( 10), Ern, Stephansen & Vohralik ( 10)... Controlling max norm error is of natural interest in some applications, such as obstacle problems, contact problems, etc.,... French, Larson & Nochetto ( 01), Nochetto, Siebert & Veeser ( 03, 05), Nochetto, Schmidt, Siebert & Veeser ( 06)... A posteriori max norm error estimates can be included in conditional estimates even for integral norm errors for non-linear evolution problems. Bartels & Müller ( 11) 4
7 Literature review: previous a posteriori L results for FEM Notation: u h denotes FE solution, Γ int interior skeleton, and [ u h ] the normal jump across element interfaces. 5
8 Literature review: previous a posteriori L results for FEM Notation: u h denotes FE solution, Γ int interior skeleton, and [ u h ] the normal jump across element interfaces. Known approaches: 1 Eriksson ( 94) 2D results (more later...) 5
9 Literature review: previous a posteriori L results for FEM Notation: u h denotes FE solution, Γ int interior skeleton, and [ u h ] the normal jump across element interfaces. Known approaches: 1 Eriksson ( 94) 2D results (more later...) 2 Nochetto ( 95), Dari, Duran & Padra ( 00): If g = 0 and the non-degeneracy condition h max h γ min holds, then u u h L (Ω) (ln 1/h min ) α d ( h 2 (f + h u h ) L (Ω) + h u h L (Γ int)) where α 2 = 2 and α 3 = 4/3. 5
10 Literature review: previous a posteriori L results for FEM Notation: u h denotes FE solution, Γ int interior skeleton, and [ u h ] the normal jump across element interfaces. Known approaches: 1 Eriksson ( 94) 2D results (more later...) 2 Nochetto ( 95), Dari, Duran & Padra ( 00): If g = 0 and the non-degeneracy condition h max h γ min holds, then u u h L (Ω) (ln 1/h min ) α d ( h 2 (f + h u h ) L (Ω) + h u h L (Γ int)) where α 2 = 2 and α 3 = 4/3. 3 Nochetto, Schmidt, Siebert & Veeser ( 06): If Ω is a Lipschitz domain, u u h L (Ω) g u h L ( Ω) + C(ln 1/h min ) α d ( h 2 (f + h u h ) L (Ω) + h u h L (Γ int)) 5
11 Literature review: previous a posteriori L proofs for FEM The basic approach: Respresent pointwise errors using a Green s function G, then prove regularity bounds for G in W 1 1, W 2 1,... 6
12 Literature review: previous a posteriori L proofs for FEM The basic approach: Respresent pointwise errors using a Green s function G, then prove regularity bounds for G in W 1 1, W 2 1,... Eriksson ( 94) (for d = 2, g = 0): Use the actual Green s function G. u(x) u h (x) =: e(x) = G(x, y) e(y)dy. Assumption: Singularities can be additively split into fundamental solution singularities and corner singularities. Ω 6
13 Literature review: previous a posteriori L proofs for FEM The basic approach: Respresent pointwise errors using a Green s function G, then prove regularity bounds for G in W 1 1, W 2 1,... Eriksson ( 94) (for d = 2, g = 0): Use the actual Green s function G. u(x) u h (x) =: e(x) = G(x, y) e(y)dy. Assumption: Singularities can be additively split into fundamental solution singularities and corner singularities. Problems: Correct in 2D? (cf. Maz ya & Rossmann ( 10), etc...) Very difficult to extend to 3D. Ω 6
14 Literature review: previous a posteriori L proofs for FEM cont d Nochetto ( 95), Dari, Duran & Padra ( 00), Nochetto, Schmidt, Siebert & Veeser ( 06): Use a regularized Green s function G with regularization parameter ρ: e(x) G e dy + Cρα e C 0,α (B ρ(x)) Ω 7
15 Literature review: previous a posteriori L proofs for FEM cont d Nochetto ( 95), Dari, Duran & Padra ( 00), Nochetto, Schmidt, Siebert & Veeser ( 06): Use a regularized Green s function G with regularization parameter ρ: e(x) G e dy + Cρα e C 0,α (B ρ(x)) Problem: e C 0,α (B ρ(x)) is an a priori term! Ω 7
16 Literature review: previous a posteriori L proofs for FEM cont d Nochetto ( 95), Dari, Duran & Padra ( 00), Nochetto, Schmidt, Siebert & Veeser ( 06): Use a regularized Green s function G with regularization parameter ρ: e(x) G e dy + Cρα e C 0,α (B ρ(x)) Problem: e C 0,α (B ρ(x)) is an a priori term! Remedies: Ω 1 Nochetto ( 95), Dari, Duran & Padra ( 00) use mesh nondegeneracy. 2 Nochetto, Schmidt, Siebert & Veeser ( 06) use Hölder regularity results to bound the a priori term: e C 0,α (Ω) e W 1 q (Ω) something computable. 7
17 Literature review: previous a posteriori L proofs for FEM cont d Nochetto ( 95), Dari, Duran & Padra ( 00), Nochetto, Schmidt, Siebert & Veeser ( 06): Use a regularized Green s function G with regularization parameter ρ: e(x) G e dy + Cρα e C 0,α (B ρ(x)) Problem: e C 0,α (B ρ(x)) is an a priori term! Remedies: Ω 1 Nochetto ( 95), Dari, Duran & Padra ( 00) use mesh nondegeneracy. 2 Nochetto, Schmidt, Siebert & Veeser ( 06) use Hölder regularity results to bound the a priori term: The latter raises two new problems: e C 0,α (Ω) e W 1 q (Ω) something computable. The analysis excludes non-lipschitz domains. 7
18 Literature review: previous a posteriori L proofs for FEM cont d Nochetto ( 95), Dari, Duran & Padra ( 00), Nochetto, Schmidt, Siebert & Veeser ( 06): Use a regularized Green s function G with regularization parameter ρ: e(x) G e dy + Cρα e C 0,α (B ρ(x)) Problem: e C 0,α (B ρ(x)) is an a priori term! Remedies: Ω 1 Nochetto ( 95), Dari, Duran & Padra ( 00) use mesh nondegeneracy. 2 Nochetto, Schmidt, Siebert & Veeser ( 06) use Hölder regularity results to bound the a priori term: The latter raises two new problems: e C 0,α (Ω) e W 1 q (Ω) something computable. The analysis excludes non-lipschitz domains. Not straightforward for dg, since u u h is not globally Hölder continuous. 7
19 Discontinuous finite element space Consider subdivision T of Ω into triangular or quadrilateral elements, which are images of the pullbacks F κ : ˆκ κ. Define piecewise polynomial space S = {v L 2 (Ω) : v F κ P r (ˆκ)} where P r (ˆκ) = 0 i+j r 0 i,j r α ij x i y j, α ij x i y j, if ˆκ triangle if ˆκ quad Ω Γ int Ω T is a shape-regular mesh derived by systematic refinement from an initial conforming mesh. Hanging nodes are OK, with the number per face bounded uniformly. Karakashian & Pascal ( 03), Bonito & Nochetto ( 10), Ainsworth & Rankin ( 10) 8
20 Notation T subdivision of Ω. Consider adjacent elements κ i and κ j [u] := u i n i + u j n j scalar jump [r] := r i n i + r j n j vector jump {u} := 1 2 (u i + u j ) scalar average {r} := 1 2 (r i + r j ) vector average κ i r i u i n j κ j r j u j n i 9
21 Interior penalty dg method Interior penalty dg method: Find u h S s. t. B(u h, v) = l(v h ) v h S, 10
22 Interior penalty dg method Interior penalty dg method: Find u h S s. t. B(u h, v) = l(v h ) v h S, where, for Γ := Ω Γ int, B(u h, v h ) := h u h h v h dx Ω Γ ( ) { u h } [v h ] + { v h } [u h ] σ[u h ] [v h ] ds 10
23 Interior penalty dg method Interior penalty dg method: Find u h S s. t. B(u h, v) = l(v h ) v h S, where, for Γ := Ω Γ int, B(u h, v h ) := h u h h v h dx Ω l(v h ) := Ω Γ ( ) { u h } [v h ] + { v h } [u h ] σ[u h ] [v h ] ds fv dx ( v h n σv h )gds Ω For h κ := diam(κ), and for e Γ, h e := {h}, we define, for C σ > 0, σ = C σ h 1 10
24 Main Result Theorem Assume that Ω R d, d = 2, 3, is an arbitrary polyhedral domain. Then ( u u h L (Ω) (ln 1/h) α d h 2 (f + h u h ) L (Ω) + h[ u h ] L (Γ int) + C σ r 2( [u h ] L (Γ int) + g u h L ( Ω)) ), where α 2 = 2 and α 3 = 1, with h := min Ω h. 11
25 Main Result Theorem Assume that Ω R d, d = 2, 3, is an arbitrary polyhedral domain. Then ( u u h L (Ω) (ln 1/h) α d h 2 (f + h u h ) L (Ω) + h[ u h ] L (Γ int) + C σ r 2( [u h ] L (Γ int) + g u h L ( Ω)) ), where α 2 = 2 and α 3 = 1, with h := min Ω h. Discussion: Results improve upon conforming FEM results in some respects. 11
26 Main Result Theorem Assume that Ω R d, d = 2, 3, is an arbitrary polyhedral domain. Then ( u u h L (Ω) (ln 1/h) α d h 2 (f + h u h ) L (Ω) + h[ u h ] L (Γ int) + C σ r 2( [u h ] L (Γ int) + g u h L ( Ω)) ), where α 2 = 2 and α 3 = 1, with h := min Ω h. Discussion: Results improve upon conforming FEM results in some respects. Improvements result from handling the Green s function differently. 11
27 Idea of proof Idea: Return to Eriksson s approach and use the standard Green s function, but with different use of regularity properties. 12
28 Sketch of proof Regularity estimates for Green s functions Lemma There exists a unique function G : Ω Ω R such that for each y Ω, each 1 q < d d 1, each 1 p < 4 3, and any r > 0, G(, y) W 1,q 0 (Ω) H 1 (Ω \ B r (y)) W 2,p (Ω \ B r (y)), G(, y) v dx = 0, v H0 1 (Ω), v = 0 in B r (y). Ω 13
29 Sketch of proof Regularity estimates for Green s functions Lemma There exists a unique function G : Ω Ω R such that for each y Ω, each 1 q < d d 1, each 1 p < 4 3, and any r > 0, G(, y) W 1,q 0 (Ω) H 1 (Ω \ B r (y)) W 2,p (Ω \ B r (y)), G(, y) v dx = 0, v H0 1 (Ω), v = 0 in B r (y). Ω In addition, for any v W 1,q 0 (Ω), q > d, we have v(y) = v G(, y) dx. Ω 13
30 Sketch of proof Regularity estimates for Green s functions Lemma There exists a unique function G : Ω Ω R such that for each y Ω, each 1 q < d d 1, each 1 p < 4 3, and any r > 0, G(, y) W 1,q 0 (Ω) H 1 (Ω \ B r (y)) W 2,p (Ω \ B r (y)), G(, y) v dx = 0, v H0 1 (Ω), v = 0 in B r (y). Ω In addition, for any v W 1,q 0 (Ω), q > d, we have v(y) = v G(, y) dx. Ω Also, for any fixed 1 q < d d 1, G(, y) W 1,q (Ω) C, where the constant C depend on Ω but is independent of y. Finally, for (x, y) Ω Ω, x y, we have { C x y 1, d = 3, G(x, y) C ln C x y, d = 2. 13
31 Sketch of proof Regularity estimates for Green s functions Lemma Assume that f L p (Ω) for some 1 < p < 4 3 and that g = 0. Then there is a unique solution u W 2,p (Ω) to the model problem and u W 2,p (Ω) C p f Lp(Ω). In addition, u W 1,q 0 (Ω) for some q > d. Grisvard ( 85), Dauge ( 88, 92), Maz ya & Rossmann ( 10) 14
32 Sketch of proof Regularity estimates for Green s functions Lemma Assume that f L p (Ω) for some 1 < p < 4 3 and that g = 0. Then there is a unique solution u W 2,p (Ω) to the model problem and u W 2,p (Ω) C p f Lp(Ω). In addition, u W 1,q 0 (Ω) for some q > d. Grisvard ( 85), Dauge ( 88, 92), Maz ya & Rossmann ( 10) Hence, point values of u are well defined. 14
33 Sketch of proof Regularity estimates for Green s functions Lemma Assume that f L p (Ω) for some 1 < p < 4 3 and that g = 0. Then there is a unique solution u W 2,p (Ω) to the model problem and u W 2,p (Ω) C p f Lp(Ω). In addition, u W 1,q 0 (Ω) for some q > d. Grisvard ( 85), Dauge ( 88, 92), Maz ya & Rossmann ( 10) Hence, point values of u are well defined. Also, we consider the boundary correction for the non-homogeneous Dirichlet conditions. 14
34 Sketch of proof cont d Problem: Green s function is not regular enough regarding dg face integrals... 15
35 Sketch of proof cont d Problem: Green s function is not regular enough regarding dg face integrals... Solution: Consider inconsistent extension B : (H 1 (Ω) + S) (H 1 (Ω) + S) R: B(w, v) := h w h v dx Ω ( ) { πv} [w] + { πw} [v] σ[w] [v] ds, Γ 15
36 Sketch of proof cont d Problem: Green s function is not regular enough regarding dg face integrals... Solution: Consider inconsistent extension B : (H 1 (Ω) + S) (H 1 (Ω) + S) R: B(w, v) := h w h v dx Ω ( ) { πv} [w] + { πw} [v] σ[w] [v] ds, Γ where π : L 1 (Ω) S is a discontinuous version of the Scott-Zhang interpolation operator onto S. 15
37 Sketch of proof cont d We have, respectively, (u u h )(x 0 ) = u(x 0 ) u h (x 0 ) 16
38 Sketch of proof cont d We have, respectively, (u u h )(x 0 ) = u(x 0 ) u h (x 0 ) = u G dx g n G ds u h (x 0 ) Ω Ω 16
39 Sketch of proof cont d We have, respectively, (u u h )(x 0 ) = u(x 0 ) u h (x 0 ) = u G dx g n G ds u h (x 0 ) Ω Ω = l(g) B(u h, G) + B(u h, G) u h (x 0 ) 16
40 Sketch of proof cont d We have, respectively, (u u h )(x 0 ) = u(x 0 ) u h (x 0 ) = u G dx g n G ds u h (x 0 ) Ω Ω = l(g) B(u h, G) + B(u h, G) u h (x 0 ) = l(g πg) B(u h, G πg) + B(u h, G) u h (x 0 ) 16
41 Sketch of proof cont d We have, respectively, (u u h )(x 0 ) = u(x 0 ) u h (x 0 ) = u G dx Ω Ω g n G ds u h (x 0 ) = l(g) B(u h, G) + B(u h, G) u h (x 0 ) = l(g πg) B(u h, G πg) + B(u h, G) u h (x 0 ) [ ] = f (G πg) dx B(u h, G πg) + σg(g πg) ds Ω Ω [ ] + g n (G πg) ds + B(u h, G) u h (x 0 ) Ω 16
42 Sketch of proof cont d We have, respectively, (u u h )(x 0 ) = u(x 0 ) u h (x 0 ) = u G dx Ω Ω g n G ds u h (x 0 ) = l(g) B(u h, G) + B(u h, G) u h (x 0 ) = l(g πg) B(u h, G πg) + B(u h, G) u h (x 0 ) [ ] = f (G πg) dx B(u h, G πg) + σg(g πg) ds Ω Ω [ ] + g n (G πg) ds + B(u h, G) u h (x 0 ) Ω =: [I ] + [II ]. 16
43 Sketch of proof cont d We have, respectively, (u u h )(x 0 ) = u(x 0 ) u h (x 0 ) = u G dx Ω Ω g n G ds u h (x 0 ) = l(g) B(u h, G) + B(u h, G) u h (x 0 ) = l(g πg) B(u h, G πg) + B(u h, G) u h (x 0 ) [ ] = f (G πg) dx B(u h, G πg) + σg(g πg) ds Ω Ω [ ] + g n (G πg) ds + B(u h, G) u h (x 0 ) Ω =: [I ] + [II ]. Remark: Term II measures in part the degree to which u h is not a conforming approximation to u. 16
44 Sketch of proof cont d Term I is more or less standard, giving, ( I h 2 (f + h u h ) L (Ω) + h[ u h ] L (Γ int) + C σ r 2 [u h ] L (Γ int) ) + C σ r 2 (u h g) L ( Ω) ( ) h 2 (G πg) L1(Ω) + h 1 (G πg) L1(Ω). 17
45 Sketch of proof cont d Term I is more or less standard, giving, ( I h 2 (f + h u h ) L (Ω) + h[ u h ] L (Γ int) + C σ r 2 [u h ] L (Γ int) ) + C σ r 2 (u h g) L ( Ω) ( ) h 2 (G πg) L1(Ω) + h 1 (G πg) L1(Ω). The second term in the product can now be bounded using the approximation properties of π. 17
46 Sketch of proof cont d For term II, we separate cases: 1 elements containing or near x 0 ; 2 elements away from the element contaning x 0. 18
47 Sketch of proof cont d For term II, we separate cases: 1 elements containing or near x 0 ; 2 elements away from the element contaning x 0. and we use the splitting into conforming and non-conforming parts à la Karakashian & Pascal ( 03) Lemma Let T be the minimal conforming refinement of T. For every v h S, there exists a function χ W 1, (Ω) which is a piecewise polynomial on T such that v h χ L (Ω) [v h ] L (Γ int) + g v h L ( Ω). 18
48 Sketch of proof cont d For term II, we separate cases: 1 elements containing or near x 0 ; 2 elements away from the element contaning x 0. and we use the splitting into conforming and non-conforming parts à la Karakashian & Pascal ( 03) Lemma Let T be the minimal conforming refinement of T. For every v h S, there exists a function χ W 1, (Ω) which is a piecewise polynomial on T such that v h χ L (Ω) [v h ] L (Γ int) + g v h L ( Ω). With a great deal of massaging and using: the stability of π in various norms; diadic decomposition about x 0, 18
49 Sketch of proof cont d For term II, we separate cases: 1 elements containing or near x 0 ; 2 elements away from the element contaning x 0. and we use the splitting into conforming and non-conforming parts à la Karakashian & Pascal ( 03) Lemma Let T be the minimal conforming refinement of T. For every v h S, there exists a function χ W 1, (Ω) which is a piecewise polynomial on T such that v h χ L (Ω) [v h ] L (Γ int) + g v h L ( Ω). With a great deal of massaging and using: the stability of π in various norms; diadic decomposition about x 0, a bound also for II can be obtained, to conclude the proof. 18
50 Efficiency bounds Theorem Under the same assumptions, the following bounds hold: h 2 (πf + u h ) L (κ) u u h L (κ) + h 2 (f πf ) L (κ), for all κ T, and h[ u h ] L (e) u u h L (κ 1 κ 2) + h 2 (f πf ) L (κ 1 κ 2), for all faces κ 1 κ 2 =: e Γ int, for all neighbouring κ 1, κ 2 T. 19
51 A side benefit: Improvement for FEM New approach yields modest improvement for max-norm a posteriori estimates for conforming FEM. 20
52 A side benefit: Improvement for FEM New approach yields modest improvement for max-norm a posteriori estimates for conforming FEM. Corollary Assume that u h is the CG solution and that Ω R d, d = 2, 3, is an arbitrary polyhedral domain. Then u u h L (Ω) g u h L ( Ω) + C(ln 1/h min ) α d ( h 2 (f + h u h ) L (Ω) + h u h L (Γ int)), where α 2 = 2 and α 3 = 1. 20
53 A side benefit: Improvement for FEM New approach yields modest improvement for max-norm a posteriori estimates for conforming FEM. Corollary Assume that u h is the CG solution and that Ω R d, d = 2, 3, is an arbitrary polyhedral domain. Then u u h L (Ω) g u h L ( Ω) + C(ln 1/h min ) α d ( h 2 (f + h u h ) L (Ω) + h u h L (Γ int)), where α 2 = 2 and α 3 = 1. Discussion: non-lipschitz domains are now included; 20
54 A side benefit: Improvement for FEM New approach yields modest improvement for max-norm a posteriori estimates for conforming FEM. Corollary Assume that u h is the CG solution and that Ω R d, d = 2, 3, is an arbitrary polyhedral domain. Then u u h L (Ω) g u h L ( Ω) + C(ln 1/h min ) α d ( h 2 (f + h u h ) L (Ω) + h u h L (Γ int)), where α 2 = 2 and α 3 = 1. Discussion: non-lipschitz domains are now included; α 3 = 1 now (as opposed to α 3 = 4/3 before). 20
55 Numerical Experiments "sol-20.gnuplot" 2D quad elements in deal.ii on an L-shaped domain: Figure: Adaptive mesh using cubic elements and 20 adaptive steps 21
56 Numerical Experiments 150 r=1 r=2 r=3 Effectivity Indices for r=1, 2, 3 Effectivity Index dof 22
57 Numerical Experiments L error dof Note: No log factors in effectivity indices or observed in error reduction. 23
58 Numerical Experiments Ω is the non-lipschitz two-brick domain of Maz ya (or Wendland?): e 1 e 2 O We take f = 1, g = 0, with unknown exact solution. 24
59 Numerical Experiments: Solution properties From theory of elliptic problems on polyhedral domains (Dauge, Maz ya,...), we know: u r 2/3 e i near edges e i (i = 1, 2), where r ei =distance to e i. Known from theory. 25
60 Numerical Experiments: Solution properties From theory of elliptic problems on polyhedral domains (Dauge, Maz ya,...), we know: u r 2/3 e i near edges e i (i = 1, 2), where r ei =distance to e i. Known from theory. u ρ λ near O, where ρ=distance to O and λ Found by numerically solving a Laplace-Beltrami eigenvalue problem on a spherical cap. 25
61 Numerical Experiments: Solution properties From theory of elliptic problems on polyhedral domains (Dauge, Maz ya,...), we know: u r 2/3 e i near edges e i (i = 1, 2), where r ei =distance to e i. Known from theory. u ρ λ near O, where ρ=distance to O and λ Found by numerically solving a Laplace-Beltrami eigenvalue problem on a spherical cap. Effects on convergence rate: The non-lipschitz corner O has little effect on AFEM convergence rate. 25
62 Numerical Experiments: Solution properties From theory of elliptic problems on polyhedral domains (Dauge, Maz ya,...), we know: u r 2/3 e i near edges e i (i = 1, 2), where r ei =distance to e i. Known from theory. u ρ λ near O, where ρ=distance to O and λ Found by numerically solving a Laplace-Beltrami eigenvalue problem on a spherical cap. Effects on convergence rate: The non-lipschitz corner O has little effect on AFEM convergence rate. Anisotropic singularities at e 1, e 2 limit convergence in L to O(DOF 2/3 ) for all r. 25
63 Numerical Experiments Figure: Eigensolution on the spherical cap for O with adaptively refined mesh. 26
64 Numerical Experiments 10 0 L error 10 1 P2 estimator P3 estimator slope= 2/ dof 27
65 Numerical Experiments Figure: Adaptive mesh using cubic elements and 4 adaptive steps 28
66 Numerical Experiments Figure: Adaptive mesh using cubic elements and 5 adaptive steps 29
67 Numerical Experiments Figure: Adaptive mesh using cubic elements and 6 adaptive steps 30
68 Numerical Experiments Figure: Adaptive mesh using cubic elements and 7 adaptive steps 31
69 Numerical Experiments Figure: Adaptive mesh using cubic elements and 8 adaptive steps 32
70 Outlook Conclusions: Reliable and efficiient max norm a posteriori estimates for (IP) dg; 33
71 Outlook Conclusions: Reliable and efficiient max norm a posteriori estimates for (IP) dg; modest improvements to the corresponding conforming FEM result; 33
72 Outlook Conclusions: Reliable and efficiient max norm a posteriori estimates for (IP) dg; modest improvements to the corresponding conforming FEM result; new theoretical side-developments and method of proof. 33
73 Outlook Conclusions: Reliable and efficiient max norm a posteriori estimates for (IP) dg; modest improvements to the corresponding conforming FEM result; new theoretical side-developments and method of proof. For the future: 3D computations with known test functions; look at non-linear problems 33
u = f in Ω, u = g on Ω, (1.1)
POINTWISE A POSTERIORI ERROR CONTROL FOR DISCONTINUOUS GALERKIN METHODS FOR ELLIPTIC PROBLEMS ALAN DEMLOW AND EMMANUIL H. GEORGOULIS Abstract. An a posteriori error bound for the maximum (pointwise) error
More informationA 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 informationConvergence 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 informationLocal pointwise a posteriori gradient error bounds for the Stokes equations. Stig Larsson. Heraklion, September 19, 2011 Joint work with A.
Local pointwise a posteriori gradient error bounds for the Stokes equations Stig Larsson Department of Mathematical Sciences Chalmers University of Technology and University of Gothenburg Heraklion, September
More informationAn 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 informationAdaptive 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 informationA Posteriori Existence in Adaptive Computations
Report no. 06/11 A Posteriori Existence in Adaptive Computations Christoph Ortner This short note demonstrates that it is not necessary to assume the existence of exact solutions in an a posteriori error
More informationThe Discontinuous Galerkin Finite Element Method
The Discontinuous Galerkin Finite Element Method Michael A. Saum msaum@math.utk.edu Department of Mathematics University of Tennessee, Knoxville The Discontinuous Galerkin Finite Element Method p.1/41
More informationEnergy 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 informationUNIFIED A POSTERIORI ERROR ESTIMATOR FOR FINITE ELEMENT METHODS FOR THE STOKES EQUATIONS
INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume 10, Number 3, Pages 551 570 c 013 Institute for Scientific Computing and Information UNIFIED A POSTERIORI ERROR ESTIMATOR FOR FINITE ELEMENT
More informationAN 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 informationAdaptive approximation of eigenproblems: multiple eigenvalues and clusters
Adaptive approximation of eigenproblems: multiple eigenvalues and clusters Francesca Gardini Dipartimento di Matematica F. Casorati, Università di Pavia http://www-dimat.unipv.it/gardini Banff, July 1-6,
More informationScientific 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 informationInstitut 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 informationL 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 informationOn a Discontinuous Galerkin Method for Surface PDEs
On a Discontinuous Galerkin Method for Surface PDEs Pravin Madhavan (joint work with Andreas Dedner and Bjo rn Stinner) Mathematics and Statistics Centre for Doctoral Training University of Warwick Applied
More informationRecovery-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 informationLocal discontinuous Galerkin methods for elliptic problems
COMMUNICATIONS IN NUMERICAL METHODS IN ENGINEERING Commun. Numer. Meth. Engng 2002; 18:69 75 [Version: 2000/03/22 v1.0] Local discontinuous Galerkin methods for elliptic problems P. Castillo 1 B. Cockburn
More informationA 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 informationA Posteriori Error Estimates For Discontinuous Galerkin Methods Using Non-polynomial Basis Functions
Lin Lin A Posteriori DG using Non-Polynomial Basis 1 A Posteriori Error Estimates For Discontinuous Galerkin Methods Using Non-polynomial Basis Functions Lin Lin Department of Mathematics, UC Berkeley;
More informationAMS 529: Finite Element Methods: Fundamentals, Applications, and New Trends
AMS 529: Finite Element Methods: Fundamentals, Applications, and New Trends Lecture 25: Introduction to Discontinuous Galerkin Methods Xiangmin Jiao SUNY Stony Brook Xiangmin Jiao Finite Element Methods
More informationAxioms 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 informationA non-standard Finite Element Method based on boundary integral operators
A non-standard Finite Element Method based on boundary integral operators Clemens Hofreither Ulrich Langer Clemens Pechstein June 30, 2010 supported by Outline 1 Method description Motivation Variational
More informationBasic 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 informationarxiv: v1 [math.na] 29 Feb 2016
EFFECTIVE IMPLEMENTATION OF THE WEAK GALERKIN FINITE ELEMENT METHODS FOR THE BIHARMONIC EQUATION LIN MU, JUNPING WANG, AND XIU YE Abstract. arxiv:1602.08817v1 [math.na] 29 Feb 2016 The weak Galerkin (WG)
More informationQUASI-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 informationScientific 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 informationA-posteriori error analysis of hp-version discontinuous Galerkin finite element methods for second-order quasilinear elliptic problems
Report no. NA-0/1 A-posteriori error analysis of hp-version discontinuous Galerkin finite element methods for second-order quasilinear elliptic problems Paul Houston 1, Endre Süli, and Thomas P. Wihler
More informationThe CG1-DG2 method for conservation laws
for conservation laws Melanie Bittl 1, Dmitri Kuzmin 1, Roland Becker 2 MoST 2014, Germany 1 Dortmund University of Technology, Germany, 2 University of Pau, France CG1-DG2 Method - Motivation hp-adaptivity
More informationOptimal Error Estimates for the hp Version Interior Penalty Discontinuous Galerkin Finite Element Method
Report no. 03/06 Optimal Error Estimates for the hp Version Interior Penalty Discontinuous Galerkin Finite Element Method Emmanuil H. Georgoulis and Endre Süli Oxford University Computing Laboratory Numerical
More informationA space-time Trefftz method for the second order wave equation
A space-time Trefftz method for the second order wave equation Lehel Banjai The Maxwell Institute for Mathematical Sciences Heriot-Watt University, Edinburgh & Department of Mathematics, University of
More informationFinite 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 informationA Multigrid Method for Two Dimensional Maxwell Interface Problems
A Multigrid Method for Two Dimensional Maxwell Interface Problems Susanne C. Brenner Department of Mathematics and Center for Computation & Technology Louisiana State University USA JSA 2013 Outline A
More informationA Framework for Analyzing and Constructing Hierarchical-Type A Posteriori Error Estimators
A Framework for Analyzing and Constructing Hierarchical-Type A Posteriori Error Estimators Jeff Ovall University of Kentucky Mathematics www.math.uky.edu/ jovall jovall@ms.uky.edu Kentucky Applied and
More informationBesov regularity for the solution of the Stokes system in polyhedral cones
the Stokes system Department of Mathematics and Computer Science Philipps-University Marburg Summer School New Trends and Directions in Harmonic Analysis, Fractional Operator Theory, and Image Analysis
More informationPREPRINT 2010:27. Local Pointwise A Posteriori Gradient Error Bounds for the Stokes Equation ALAN DEMLOW STIG LARSSON
PREPRIN 2010:27 Local Pointwise A Posteriori Gradient Error Bounds for the Stokes Equation ALAN DEMLOW SIG LARSSON Department of Mathematical Sciences Division of Mathematics CHALMERS UNIVERSIY OF ECHNOLOGY
More informationProjected Surface Finite Elements for Elliptic Equations
Available at http://pvamu.edu/aam Appl. Appl. Math. IN: 1932-9466 Vol. 8, Issue 1 (June 2013), pp. 16 33 Applications and Applied Mathematics: An International Journal (AAM) Projected urface Finite Elements
More informationA space-time Trefftz method for the second order wave equation
A space-time Trefftz method for the second order wave equation Lehel Banjai The Maxwell Institute for Mathematical Sciences Heriot-Watt University, Edinburgh Rome, 10th Apr 2017 Joint work with: Emmanuil
More informationA posteriori error estimates for non conforming approximation of eigenvalue problems
A posteriori error estimates for non conforming approximation of eigenvalue problems E. Dari a, R. G. Durán b and C. Padra c, a Centro Atómico Bariloche, Comisión Nacional de Energía Atómica and CONICE,
More informationwhere N h stands for all nodes of T h, including the boundary nodes. We denote by N h the set of interior nodes. We have
48 R.H. NOCHETTO 5. Lecture 5. Pointwise Error Control and Applications In this lecture we extend the theory of Lecture 2 to the maximum norm following [38, 39]. We do this in sections 5.1 to 5.3. The
More informationFinite element approximation of eigenvalue problems
Acta NumeTica (2010), pp. 1-120 doi: 10.1017/80962492910000012 Cambridge University Press, 2010 Finite element approximation of eigenvalue problems Daniele Boffi Dipartimento di Matematica, Universita
More informationLOCAL A POSTERIORI ESTIMATES FOR POINTWISE GRADIENT ERRORS IN FINITE ELEMENT METHODS FOR ELLIPTIC PROBLEMS
LOCAL A POSTERIORI ESTIMATES FOR POINTWISE GRADIENT ERRORS IN FINITE ELEMENT METHODS FOR ELLIPTIC PROBLEMS ALAN DEMLOW Abstract. We prove local a posteriori error estimates for pointwise gradient errors
More informationNumerical 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 informationANALYSIS OF THE FEM AND DGM FOR AN ELLIPTIC PROBLEM WITH A NONLINEAR NEWTON BOUNDARY CONDITION
Proceedings of EQUADIFF 2017 pp. 127 136 ANALYSIS OF THE FEM AND DGM FOR AN ELLIPTIC PROBLEM WITH A NONLINEAR NEWTON BOUNDARY CONDITION MILOSLAV FEISTAUER, ONDŘEJ BARTOŠ, FILIP ROSKOVEC, AND ANNA-MARGARETE
More informationOn 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 informationA POSTERIORI ENERGY-NORM ERROR ESTIMATES FOR ADVECTION-DIFFUSION EQUATIONS APPROXIMATED BY WEIGHTED INTERIOR PENALTY METHODS
A posteriori error estimates for weighted interior penalty methods A POSERIORI ENERGY-NORM ERROR ESIMAES FOR ADVECION-DIFFUSION EQUAIONS APPROXIMAED BY WEIGHED INERIOR PENALY MEHODS Alexandre Ern (Université
More informationRemarks on the analysis of finite element methods on a Shishkin mesh: are Scott-Zhang interpolants applicable?
Remarks on the analysis of finite element methods on a Shishkin mesh: are Scott-Zhang interpolants applicable? Thomas Apel, Hans-G. Roos 22.7.2008 Abstract In the first part of the paper we discuss minimal
More informationarxiv: 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 informationA 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 informationPolynomial-degree-robust liftings for potentials and fluxes
Polynomial-degree-robust liftings for potentials and fluxes and Martin Vohraĺık Université Paris-Est, CERMICS (ENPC) and SERENA Project-Team, INRIA, France RAFEM, Hong-Kong, March 2017 Outline 1. Motivation:
More informationENERGY 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 informationhp-version Discontinuous Galerkin Finite Element Methods for Semilinear Parabolic Problems
Report no. 3/11 hp-version Discontinuous Galerkin Finite Element Methods for Semilinear Parabolic Problems Andris Lasis Endre Süli We consider the hp version interior penalty discontinuous Galerkin finite
More informationb i (x) u + c(x)u = f in Ω,
SIAM J. NUMER. ANAL. Vol. 39, No. 6, pp. 1938 1953 c 2002 Society for Industrial and Applied Mathematics SUBOPTIMAL AND OPTIMAL CONVERGENCE IN MIXED FINITE ELEMENT METHODS ALAN DEMLOW Abstract. An elliptic
More informationarxiv: v1 [math.na] 19 Nov 2018
QUASINORMS IN SEMILINEAR ELLIPTIC PROBLEMS JAMES JACKAMAN AND TRISTAN PRYER arxiv:1811.07816v1 [math.na] 19 Nov 2018 (1.1) Abstract. In this note we examine the a priori and a posteriori analysis of discontinuous
More informationA posteriori energy-norm error estimates for advection-diffusion equations approximated by weighted interior penalty methods
A posteriori energy-norm error estimates for advection-diffusion equations approximated by weighted interior penalty methods Alexandre Ern, Annette Stephansen o cite this version: Alexandre Ern, Annette
More informationDifferent Approaches to a Posteriori Error Analysis of the Discontinuous Galerkin Method
WDS'10 Proceedings of Contributed Papers, Part I, 151 156, 2010. ISBN 978-80-7378-139-2 MATFYZPRESS Different Approaces to a Posteriori Error Analysis of te Discontinuous Galerkin Metod I. Šebestová Carles
More informationAn 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 informationAn hp-adaptive Mixed Discontinuous Galerkin FEM for Nearly Incompressible Linear Elasticity
An hp-adaptive Mixed Discontinuous Galerkin FEM for Nearly Incompressible Linear Elasticity Paul Houston 1 Department of Mathematics, University of Leicester, Leicester LE1 7RH, UK (email: Paul.Houston@mcs.le.ac.uk)
More informationSECOND ORDER TIME DISCONTINUOUS GALERKIN METHOD FOR NONLINEAR CONVECTION-DIFFUSION PROBLEMS
Proceedings of ALGORITMY 2009 pp. 1 10 SECOND ORDER TIME DISCONTINUOUS GALERKIN METHOD FOR NONLINEAR CONVECTION-DIFFUSION PROBLEMS MILOSLAV VLASÁK Abstract. We deal with a numerical solution of a scalar
More informationRate optimal adaptive FEM with inexact solver for strongly monotone operators
ASC Report No. 19/2016 Rate optimal adaptive FEM with inexact solver for strongly monotone operators G. Gantner, A. Haberl, D. Praetorius, B. Stiftner Institute for Analysis and Scientific Computing Vienna
More informationKey words. Laplace-Beltrami operator, surface finite element methods, a posteriori error estimates, adaptive finite element methods
A POSTERIORI ERROR ESTIMATES FOR THE LAPLACE-BELTRAMI OPERATOR ON PARAMETRIC C 2 SURFACES ANDREA BONITO AND ALAN DEMLOW Key words. Laplace-Beltrami operator, surface finite element methods, a posteriori
More informationNew class of finite element methods: weak Galerkin methods
New class of finite element methods: weak Galerkin methods Xiu Ye University of Arkansas at Little Rock Second order elliptic equation Consider second order elliptic problem: a u = f, in Ω (1) u = 0, on
More informationINTRODUCTION 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 informationThe Mortar Boundary Element Method
The Mortar Boundary Element Method A Thesis submitted for the degree of Doctor of Philosophy by Martin Healey School of Information Systems, Computing and Mathematics Brunel University March 2010 Abstract
More informationChapter 6 A posteriori error estimates for finite element approximations 6.1 Introduction
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
More informationhp-dgfem FOR SECOND-ORDER MIXED ELLIPTIC PROBLEMS IN POLYHEDRA
MATHEMATICS OF COMPUTATION Volume 00, Number 0, Pages 000 000 S 0025-5718(XX)0000-0 hp-dgfem FOR SECOND-ORDER MIXED ELLIPTIC PROBLEMS IN POLYHEDRA DOMINIK SCHÖTZAU, CHRISTOPH SCHWAB, AND THOMAS P. WIHLER
More informationMaximum-norm stability of the finite element Ritz projection with mixed boundary conditions
Noname manuscript No. (will be inserted by the editor) Maximum-norm stability of the finite element Ritz projection with mixed boundary conditions Dmitriy Leykekhman Buyang Li Received: date / Accepted:
More informationAn a posteriori error indicator for discontinuous Galerkin discretizations of H(curl) elliptic partial differential equations
IMA Journal of Numerical Analysis 005) Page of 7 doi: 0.093/imanum/ An a posteriori error indicator for discontinuous Galerkin discretizations of Hcurl) elliptic partial differential equations PAUL HOUSTON
More informationENERGY 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 informationGraded 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 informationPREPRINT 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 informationAn a posteriori error estimator for the weak Galerkin least-squares finite-element method
An a posteriori error estimator for the weak Galerkin least-squares finite-element method James H. Adler a, Xiaozhe Hu a, Lin Mu b, Xiu Ye c a Department of Mathematics, Tufts University, Medford, MA 02155
More informationInstitut de Recherche MAthématique de Rennes
LMS Durham Symposium: Computational methods for wave propagation in direct scattering. - July, Durham, UK The hp version of the Weighted Regularization Method for Maxwell Equations Martin COSTABEL & Monique
More informationMaximum norm estimates for energy-corrected finite element method
Maximum norm estimates for energy-corrected finite element method Piotr Swierczynski 1 and Barbara Wohlmuth 1 Technical University of Munich, Institute for Numerical Mathematics, piotr.swierczynski@ma.tum.de,
More information1. Introduction. The Stokes problem seeks unknown functions u and p satisfying
A DISCRETE DIVERGENCE FREE WEAK GALERKIN FINITE ELEMENT METHOD FOR THE STOKES EQUATIONS LIN MU, JUNPING WANG, AND XIU YE Abstract. A discrete divergence free weak Galerkin finite element method is developed
More informationarxiv: v1 [math.na] 1 May 2013
arxiv:3050089v [mathna] May 03 Approximation Properties of a Gradient Recovery Operator Using a Biorthogonal System Bishnu P Lamichhane and Adam McNeilly May, 03 Abstract A gradient recovery operator based
More informationADAPTIVE 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 informationAn 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 informationCorner Singularities
Corner Singularities Martin Costabel IRMAR, Université de Rennes 1 Analysis and Numerics of Acoustic and Electromagnetic Problems Linz, 10 15 October 2016 Martin Costabel (Rennes) Corner Singularities
More informationPREPRINT 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 informationHybridized Discontinuous Galerkin Methods
Hybridized Discontinuous Galerkin Methods Theory and Christian Waluga Joint work with Herbert Egger (Uni Graz) 1st DUNE User Meeting, Stuttgart Christian Waluga (AICES) HDG Methods October 6-8, 2010 1
More informationGreen s Functions and Distributions
CHAPTER 9 Green s Functions and Distributions 9.1. Boundary Value Problems We would like to study, and solve if possible, boundary value problems such as the following: (1.1) u = f in U u = g on U, where
More informationSpace-time sparse discretization of linear parabolic equations
Space-time sparse discretization of linear parabolic equations Roman Andreev August 2, 200 Seminar for Applied Mathematics, ETH Zürich, Switzerland Support by SNF Grant No. PDFMP2-27034/ Part of PhD thesis
More informationMultigrid 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 informationScientific Computing I
Scientific Computing I Module 8: An Introduction to Finite Element Methods Tobias Neckel Winter 2013/2014 Module 8: An Introduction to Finite Element Methods, Winter 2013/2014 1 Part I: Introduction to
More informationMultigrid Methods for Maxwell s Equations
Multigrid Methods for Maxwell s Equations Jintao Cui Institute for Mathematics and Its Applications University of Minnesota Outline Nonconforming Finite Element Methods for a Two Dimensional Curl-Curl
More informationWeighted Regularization of Maxwell Equations Computations in Curvilinear Polygons
Weighted Regularization of Maxwell Equations Computations in Curvilinear Polygons Martin Costabel, Monique Dauge, Daniel Martin and Gregory Vial IRMAR, Université de Rennes, Campus de Beaulieu, Rennes,
More informationOptimal multilevel preconditioning of strongly anisotropic problems.part II: non-conforming FEM. p. 1/36
Optimal multilevel preconditioning of strongly anisotropic problems. Part II: non-conforming FEM. Svetozar Margenov margenov@parallel.bas.bg Institute for Parallel Processing, Bulgarian Academy of Sciences,
More informationR 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 informationNumerical Analysis of Higher Order Discontinuous Galerkin Finite Element Methods
Numerical Analysis of Higher Order Discontinuous Galerkin Finite Element Methods Contents Ralf Hartmann Institute of Aerodynamics and Flow Technology DLR (German Aerospace Center) Lilienthalplatz 7, 3808
More informationRobust exponential convergence of hp-fem for singularly perturbed systems of reaction-diffusion equations
Robust exponential convergence of hp-fem for singularly perturbed systems of reaction-diffusion equations Christos Xenophontos Department of Mathematics and Statistics University of Cyprus joint work with
More informationAnalysis of Hybrid Discontinuous Galerkin Methods for Incompressible Flow Problems
Analysis of Hybrid Discontinuous Galerkin Methods for Incompressible Flow Problems Christian Waluga 1 advised by Prof. Herbert Egger 2 Prof. Wolfgang Dahmen 3 1 Aachen Institute for Advanced Study in Computational
More informationBUBBLE STABILIZED DISCONTINUOUS GALERKIN METHOD FOR STOKES PROBLEM
BUBBLE STABILIZED DISCONTINUOUS GALERKIN METHOD FOR STOKES PROBLEM ERIK BURMAN AND BENJAMIN STAMM Abstract. We propose a low order discontinuous Galerkin method for incompressible flows. Stability of the
More information[2] (a) Develop and describe the piecewise linear Galerkin finite element approximation of,
269 C, Vese Practice problems [1] Write the differential equation u + u = f(x, y), (x, y) Ω u = 1 (x, y) Ω 1 n + u = x (x, y) Ω 2, Ω = {(x, y) x 2 + y 2 < 1}, Ω 1 = {(x, y) x 2 + y 2 = 1, x 0}, Ω 2 = {(x,
More informationA Posteriori Error Control for Discontinuous Galerkin Methods for Parabolic Problems
A Posteriori Error Control for Discontinuous Galerkin Methods for Parabolic Problems EMMANUIL H. GEORGOULIS AND OMAR LAKKIS Abstract. We derive energy-norm a posteriori error bounds for an Euler timestepping
More informationUnified A Posteriori Error Control for all Nonstandard Finite Elements 1
Unified A Posteriori Error Control for all Nonstandard Finite Elements 1 Martin Eigel C. Carstensen, C. Löbhard, R.H.W. Hoppe Humboldt-Universität zu Berlin 19.05.2010 1 we know of Guidelines for Applicants
More informationChapter 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 informationELLIPTIC 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 informationLECTURE 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 informationOn Surface Meshes Induced by Level Set Functions
On Surface Meshes Induced by Level Set Functions Maxim A. Olshanskii, Arnold Reusken, and Xianmin Xu Bericht Nr. 347 Oktober 01 Key words: surface finite elements, level set function, surface triangulation,
More information