Polynomial-degree-robust liftings for potentials and fluxes
|
|
- Lora Stokes
- 5 years ago
- Views:
Transcription
1 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
2 Outline 1. Motivation: a posteriori error estimates 2. Main results on H 1 and H(div) polynomial liftings 3. Main ingredients of proofs 4. Numerical results (with V. Dolejší) 5. Global H(div) polynomial liftings (with I. Smears)
3 A posteriori error estimates Model problem in Ω R d with data f L 2 (Ω) u H 1 0 (Ω) s.t. ( u, v) Ω = (f, v) Ω, v H 1 0 (Ω) H 1 0 -conforming FEM solution u h on a simplicial mesh T h Fully computable error upper bound (u u h ) 2 Ω K T h η 2 K Local error indicators η K represent local error lower bounds up to data oscillation η K C eff (u u h ) ωk + osc(f, ω K ) Pioneered by [Babuška, Rheinbolt 78]; recent textbook [Verfürth 13] classical technique to compute the ηk s is residual-based drawback: upper bound features an undetermined constant
4 Equilibrated flux reconstruction Exact flux σ := u H(div, Ω) s.t. σ = f (equilibrium) An equilibrated flux reconstruction is s.t. σ h H(div, Ω) ( σ h, 1) K = (f, 1) K, K T h Setting η F,K := u h + σ h K and η osc,k := h K π f σ h K, (u u h ) 2 Ω K T h (η F,K + η osc,k }{{} =η K Literature: hypercircle method [Prager, Synge 47]; computational mechanics [Ladevèze et al. 75]; textbooks [Ainsworth, Oden 00; Repin 08] ) 2
5 A simple proof Residual ρ(u h ) H 1 (Ω) s.t. ρ(u h ), ϕ Ω := (f, ϕ) Ω ( u h, ϕ) Ω, ϕ H 1 0 (Ω) (u u h ) Ω = ρ(u h ) H 1 (Ω) = Introduce H(div, Ω)-flux and use Green s formula sup ρ(u h ), ϕ Ω ϕ H0 1(Ω); ϕ Ω=1 ρ(u h ), ϕ Ω = (f σ h, ϕ) Ω ( u h + σ h, ϕ) Ω Cauchy Schwarz and Poincaré Steklov inequalities (equilibration) (f σ h, ϕ) Ω = K T h (f σ h, ϕ ϕ K ) K K T h h K π f σ h K ϕ K ( u h + σ h, ϕ) Ω K T h u h + σ h K ϕ K Poincaré (1894) [eigenvalue pb], Steklov (1897) [d = 1], Payne, Weinberger (60) [d = 2], Bebendorf (03) [d 3]
6 a5 a4 a2 a3 How to build σ h? Global flux equilibration (... expensive) σ h := arg min u h + v h Ω v h V h, v h =Π ( Vh )f for some global Raviart Thomas spaces V h H(div, Ω) Cheap local flux equilibration by working on FE stars Ta T h : patch of elements sharing vertex a V h ; domain ω a a Vh a1 patch ωa ψa(a) =1,ψa(a ) =0 Literature on a posteriori analysis on FE stars: [Babuška, Miller 87; Morin, Nochetto, Siebert 03; Cohen, DeVore, Nochetto 12; Veeser, Verfürth 12]
7 Local flux equilibration on FE stars Hat basis functions {ψ a } a Vh with local PU a V K ψ a K = 1 We focus for simplicity on FE stars around interior vertices flux equilibration for boundary vertices depends on the BC s for model problem, see [Dolejší, AE, MV, 16] Raviart Thomas space V a h H 0(div, ω a ) (Neumann BC s on ω a ) σ a h := arg min ψ a u h + v h ωa v h Vh a, v h=gh a σ h := a V h σ a h with data gh a := Π ( Vh a) (f ψ a u h ψ a ) a V K gh a = Π ( Vh a) f (PU), (gh a, 1) ω a = 0 (Galerkin s orthogonality on hat basis functions) Literature [Destuynder, Métivet 99; Braess, Schöberl 08; MV 08; AE, MV 10]
8 p-robust local efficiency Let p be the polynomial degree used to compute u h Local flux equilibration using RT spaces of order p Key result: Local lower error bound (ω K = a VK ω a ) η F,K = u h + σ h K C eff (u u h ) ωk + osc(f, ω K ) C eff depends on patch geometry (mesh regularity) Ceff is p-robust for d = 2 [Braess, Pillwein, Schöberl 09] Ceff is p-robust for d = 3 [AE, MV 16]... in contrast to residual-based estimators [Melenk, Wohlmuth 01]
9 A two-step proof Two-step proof using -dimensional, local problems replaces classical bubble-function argument by Verfürth see [AE, MV 15] We have C eff = C st C PS Cst results from p-robust stability properties of RT spaces [BPS 09] CPS results from Poincaré Steklov inequalities on FE stars [Carstensen, Funken 00], values estimated in [Veeser, Verfürth 12]
10 Step 1 Set τ a h := ψ a u h and V a h H 0(div, ω a ) = V a Local PU yields u h + σ h K a V K τ a h + σa h ω a Recall that local flux equilibration delivers σ a h := arg min τ v h Vh a, v h a + v h ωa h=gh a Local -dimensional constrained min. pb. with polynomial data σ a := arg min τh a + v ωa v V a, v=gh a Thm. The following (nontrivial) result holds: τh a + σh a ωa C st τh a + σ a ωa Note the (trivial) converse bound τ a h + σa ωa τ a h + σa h ωa
11 Step 2 Evaluate τh a + σa ωa using equivalent problem in primal form r a H (ω 1 a) = {v H 1 (ω a) (v, 1) ωa = 0} s.t. ( r a, v) ωa = (τh a, v) ωa + (gh a, v) ωa v H (ω 1 a) Equivalence of primal/dual energies τ a h + σ a ωa = min τ a v V a, v=g h a h + v ωa = r a ωa Since ( r a, v) ωa = ( (u u h ), (ψ av)) ωa + osc(f, ω a) (ψav) ωa (1 + C PS,ωa h ωa ψ a L (ω a)) v ωa we can conclude that τh a + σ a ωa C PS (u u h ) ωa + osc(f, ω a ) with C PS = max a Vh (1 + C PS,ωa h ωa ψ a L (ω a))
12 Nonconforming case (u h H 1 0 (Ω)) Error measure w.r.t. broken gradient T of discrete solution Additional nonconformity estimator in error upper bound H 1 0 -potential reconstruction s h H 1 0 (Ω) Setting η NC,K := T (u h s h ), T (u u h ) 2 ( ) 2 ( ) 2 ηf,k + η osc,k + ηnc,k K T h K T h Typically, s h is built by prescribing its nodal values as averages see [Achdou, Bernardi, Coquel 03; Karakashian, Pascal, 03] for Crouzeix Raviart and IPDG residual-based estimates... not p-robust [Burman, AE 07; Houston, Schötzau, Wihler 07]
13 Potential reconstruction Local FEM solves of order (p + 1) on vertex-based patches s a h P p+1 (T a) H0 1 (ω a) set sh = a V h sh a p-robust local efficiency proved in [AE, MV 15] in 2D assuming [[uh ]], 1 F = 0, η NC,K C eff T (u u h ) ωa, C eff = C stc bps a V K two-step proof as above (no oscillation here) 1st step: mixed RT solve of order p for rotated gradient of s a h 2nd step: broken PS inequalities, see also [Carstensen, Merdon 13] jump seminorm added to error and estimator if [[uh ]], 1 F 0 use discrete gradient for DG (instead of broken gdt); broken PS and jump seminorm can be avoided for symmetric IPDG [AE, MV 15]
14 Main results [AE, MV 16] 3D, p-robust flux and potential reconstructions Thm. 1 p-robust, H 1 -stable polynomial lifting for potentials Thm. 2 p-robust, H(div)-stable polynomial lifting for fluxes Our proofs are constructive (as 2D proofs) e.g., build σ a h Vh a s.t. σ h a = gh a, τh a + σ h a ωa C st τh a + σ a ωa then τ a h + σh a ωa τh a + σ h a ωa C st τh a + σ a ωa Main challenges how to enumerate tetrahedra in 3D star (triangles in 2D star are enumerated by circling around the vertex) need to work with potentials (and not with rotated gradients) We focus for simplicity on FE stars around interior vertices
15 Some notation T a T h : FE star (elements sharing vertex a V h ); domain ω a F a = F s a F b a : faces of the elements in the star T a (2D) skeletal faces F s a (2D) boundary faces F b a Broken H 1 - and H(div)-spaces (broken gradient T ) H 1 (T a) := {v L 2 (ω a) v K H 1 (K) K T a} H(div, T a) := {v L 2 (ω a) v K H(div, K) K T a} Broken polynomial subspaces P p (T a) := {v L 2 (ω a) v K P p (K) K T a} RT p (T a) := {v L 2 (ω a) v K RT p K) K T a}
16 Some notation (cont d) In 3D, a FE star ω a is homeomorphic to a ball in R 3 We can look at the star boundary ω a the traces of the tetrahedra in Ta form a triangulation of ω a every triangle is a boundary face F F b a every edge is the trace of a skeletal face F F s a every point is the trace of a skeletal edge e E a Orientation every skeletal face is oriented so as to define a jump across it every skeletal edge is oriented so as to circle around it incidence coefficients ιf,e = ±1, for all F F e and e E a
17 H 1 -stable polynomial lifting Let p 1 Let r F P p (F ) for all F F a Assume r F 0 on boundary faces F Fa b, and assume on skeletal faces F Fa s the compatibility conditions r F F ωa = 0, ι F,e r F e = 0 e E a F F e Then, min v h P p (T a) v h =r F F F b a [v h ]=r F F F s a T v h ωa C st min v H 1 (T a) v=r F F F b a [v ]=r F F F s a T v ωa with p-robust constant C st only depending on mesh regularity
18 H(div)-stable polynomial lifting Let p 0 Let r K P p (K) for all K T a and let r F P p (F ) for all F F a Assume the compatibility condition (r K, 1) K (r F, 1) F = 0 K T a F F a Then, min v h RT p (T a) v h n F =r F F F b a [v h n F ]=r F F F s a T v h K =r K K T a v h ωa C st min v H(div,T a) v n F =r F F F b a [v n F ]=r F F F s a T v K =r K K T a v ωa with p-robust constant C st only depending on mesh regularity
19 Shifted reformulation: potential Let ξ a h Pp (T a ) be any function from the minimization set ξ a h = r F F F b a, [[ξ a h]] = r F F F s a An equivalent reformulation of Thm. 1 is min T (ξ a v h P p (T a) H0 1(ωa) h v h ) ωa C st min T (ξ a v H0 1(ωa) h v) ωa Application to a posteriori error analysis: ξ a h = ψ au h and r F = 0 F F b a, r F = ψ a [[u h ]] F F s a The compatibility conditions F F e ι F,e r F e = 0, e E a, follow from algebraic properties of jump operator
20 Shifted reformulation: flux Let τ a h RTp (T a ) be any function s.t. τ a h n F = r F F F b a, [[τ a h ]] n F = r F F F s a An equivalent reformulation of Thm. 2 is min τ v h RT p h a +v h ωa C st min (T a) H 0(div,ω a) v h K =r K τ a h K K T a v H 0(div,ω a) v K =r K τ a h K Application to a posteriori error analysis: τ a h = ψ a T u h and r F = 0 F Fa b, r F = ψ a [[ T u h ]] n F F Fa s r K = ψ a (f + T u h ) K T a (f pcw. polynomial) K T a τ a h +v ωa The compatibility condition K T a (r K, 1) K F F a (r F, 1) F = 0 is nothing but Galerkin s orthogonality on the hat basis function ψ a
21 Main ingredients of proofs On a fixed tetrahedron K T a, we can Lift the prescribed divergence of the flux using [Costabel, McIntosh 10] (valid in any space dimension) Lift prescribed polynomials for the flux normal component using [Demkowicz, Gopalakrishnan, Schöberl 12] (proved for d = 3) (a compatibility condition is required if the prescription is on all faces) Lift prescribed compatible polynomials for the potential trace using [DGS 09] (proved for d = 3) We are left with the p-robust lifting of the prescribed jumps... but this requires a careful enumeration of the tetrahedra in the star
22 2D enumeration Circle around the interior vertex a K1: do nothing Kn: fix jump on face touching K n 1, n {2...4} K5: fix last two jumps (possible owing to compatibility condition)
23 3D enumeration: shellability Let T a be a star of tetrahedra around the interior vertex a Consider the triangulation of the sphere S (2) R 3 with same connectivity Enumerate surface triangles so that j i T (2) j is connected for all i The notion of shellability of polytopes shows that this is possible [Ziegler, Lectures on Polytopes, Chap. 8, Springer, 2006] Enumerate patch tetrahedra following surface triangle enumeration and, for each K i T a, fix jump on the skeletal faces of K i touching any tetrahedron K j for j < i
24 Numerical results (with V. Dolejší) Smooth analytical solution in Ω = (0, 1) 2 u(x1, x 2) = sin(2πx 1) sin(2πx 2) Uniformly refined, non-nested unstructured triangulations Discretization by symmetric IPDG method asymptotic exactness observed for pol. degrees p {1...6} similar results for incomplete version of IPDG slightly larger effectivity indices for nonsymmetric version and even p
25 Errors, estimators, effectivity indices h/h 0 p e j(e) e + j(e) η F ηosc η NC η η + j(u h ) I eff I j eff E E E E E E E E / E E E E E E E E / E E E E E E E E / E E E E E E E E E E E E E E E E / E E E E E E E E / E E E E E E E E / E E E E E E E E E E E E E E E E / E E E E E E E E / E E E E E E E E / E E E E E E E E E E E E E E E E / E E E E E E E E / E E E E E E E E / E E E E E E E E E E E E E E E E / E E E E E E E E / E E E E E E E E / E E E E E E E E E E E E E E E E / E E E E E E E E / E E E E E E E E
26 hp-adaptivity Nested simplicial meshes allowing for hanging nodes Extension of reconstruction procedures to hanging nodes only matching refinement of individual patches is needed Bulk chasing criterion based on local p-robust estimators hp-refinement decision criteria inspired from [Mitchell, McClain 14] Algebraic convergence w.r.t. dof s observed on several 2D benchmark problems from [Mitchell 13]
27 Numerical examples Re-entrant corner singularity (u H 1+t (Ω), t < 2 3 ) error in energy norm 1e-01 1e-02 1e-03 1e-04 ideal param, γ =0.30 param, γ =0.60 prior hp decay h-adapt 1e Ε Ε Ε Ε Ε Ε 02 DoF 1/3 Multiple difficulties (re-entrant corner, point sing, circular wave) error in energy norm param, γ =0.30 param, γ =0.60 prior hp decay DoF 1/3 1.0Ε Ε Ε+00
28 Global H(div)-liftings (with I. Smears) So far, we devised local liftings of polynomial data on FE stars We now consider global H(div)-liftings of polynomial data One important application is to devise liftings on patches of mesh cells that are (much) larger than a star This allows us to remove some theoretical barriers on number of hanging nodes for elliptic problems level of mesh coarsening between time-steps in parabolic problems T n 1 \ T n a a T n T a,n ωa ω a
29 Main result Lipschitz domain Ω R d ; boundary partition Γ = Γ D Γ N Let T h be a simplicial mesh of Ω, possibly locally refined Thm. Let p 1, f P p 1 (T h ), ξ RT p 1 (T h ) (broken spaces) (with (f, 1) Ω = 0 if Γ N = Γ). Then, min v h RT p (T h ) H(div,Ω) v h =f in Ω v h n=0 on Γ N ξ + v h Ω C st min v H(div,Ω) v=f in Ω v n=0 on Γ N ξ + v Ω with p-robust constant C st only depending on mesh regularity (Note that the discrete minimizer is one polynomial order higher than the data) See [Ainsworth, Guzman, Sayas 16] for zero interior source terms, nonzero boundary traces, and fixed polynomial degree
30 Main idea in proof Primal problem with H 1 data u H (Ω) 1 s.t. ( u, v) Ω = (f, v) Ω (ξ, v) Ω, v H (Ω) 1 H (Ω) 1 = {v H 1 (Ω) (v, 1) Ω = 0} if Γ N = Γ H 1(Ω) = {v H1 (Ω) v ΓD = 0} otherwise Equivalence of primal/dual energies (f, v) Ω (ξ, v) Ω min ξ + v Ω = max = u Ω v H(div,Ω) v H 1 (Ω) v Ω v=f in Ω v n=0 on Γ N Need to construct σ h RT p (T h ) H(div, Ω) s.t. ξ + σ h Ω C st u Ω
31 A posteriori error estimate with H 1 data Consider general data f L 2 (Ω) and ξ L 2 (Ω) u H (Ω) 1 s.t. ( u, v) Ω = (f, v) Ω (ξ, v) Ω, v H (Ω) 1 H -conforming 1 approximation u h P p (T h ) H (Ω), 1 p 1 Equilibrated flux reconstruction σ h RT p (T h ) H(div, Ω), p p (u u h ) 2 Ω ( uh + ξ + σ h K + osc(f, K) ) 2 K T h u h + ξ + σ h K C eff (u u h ) ωk + osc(f, ξ, ω K ) σ h constructed locally from local shifted flux equilibration in V a h = RTp (T a ) H(div, ω a ) (+ Neumann BC s) σ a h := arg min τ v h Vh a, v h a + v h ωa h=gh a with data τ a h := ψ a(ξ + u h ), g a h := Π ( V a h ) (f ψ a (ξ + u h ) ψ a )
32 Conclusion of proof Consider polynomial data f P p 1 (T h ), ξ RT p 1 (T h ), p 1 Consider H 1 -conforming FEM approximation with 1 = p p The data oscillation term osc(f, ξ, ω K ) vanishes so that u h + ξ + σ h Ω C eff (u u h ) Ω Combined with the basic bound u h Ω u Ω, we conclude that ξ + σ h Ω u h + ξ + σ h Ω + u h Ω (2C eff + 1) u Ω
33 Sharper stability result for special data Let ψ be continuous, pcw. affine on T h ψ can be viewed as a large-scale hat basis function Let Γ N be the subset of Γ where ψ vanishes Let p 1, f P p 1 (T h ), ξ RT p 1 (T h ) (with (f, 1) Ω = 0 if Γ N = Γ) Thm. The following holds: min v h RT p (T h ) H(div,Ω) v h =ψ f ψ ξ in Ω v h n=0 on Γ N ξ + v h Ω C st C(ψ ) min v H(div,Ω) v=f in Ω no BC ξ + v Ω with p-robust constant C st only depending on mesh regularity and C(ψ ) = ψ L (Ω) + h Ω ψ L (Ω) lifting and data have matching polynomial order the function ψ is factored out from RHS
34 Conclusions Equilibrated-flux estimates offer several benefits guaranteed (fully computable) upper bounds p-robust local efficiency adaptive inexact Newton solvers [AE, MV 13] Unified analysis for p-robust H 1 - and H(div)-polynomial liftings New local efficiency proofs for arbitrary-level of hanging nodes (elliptic PDEs) and no coarsening restriction (parabolic PDEs) Thank you for your attention 1. AE, MV, SINUM (2015), 53, V. Dolejší, AE, MV, SISC (2016), 38 A3220 A AE, MV (2016), arxiv AE, I. Smears, MV, Calcolo (2017); DOI /s
Polynomial-degree-robust a posteriori estimates in a unified setting for conforming, nonconforming, discontinuous Galerkin, and mixed discretizations
Polynomial-degree-robust a posteriori estimates in a unified setting for conforming, nonconforming, discontinuous Galerkin, and mixed discretizations Alexandre Ern, Martin Vohralík To cite this version:
More informationA 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 informationEquilibrated flux a posteriori error estimates in L 2 (H 1 )-norms for high-order discretizations of parabolic problems
Equilibrated flux a posteriori error estimates in L 2 (H 1 )-norms for high-order discretizations of parabolic problems Alexandre Ern, Iain Smears, Martin Vohralík To cite this version: Alexandre Ern,
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 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 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 informationA NOTE ON CONSTANT-FREE A POSTERIORI ERROR ESTIMATES
A NOTE ON CONSTANT-FREE A POSTERIORI ERROR ESTIMATES R. VERFÜRTH Abstract. In this note we look at constant-free a posteriori error estimates from a different perspective. We show that they can be interpreted
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 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 informationA 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 informationEquilibrated Residual Error Estimates are p-robust
Equilibrated Residual Error Estimates are p-robust Dietrich Braess, Veronika Pillwein and Joachim Schöberl April 16, 2008 Abstract Equilibrated residual error estimators applied to high order finite elements
More informationAdaptive 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 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 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 informationDiscontinuous Galerkin Methods: Theory, Computation and Applications
Discontinuous Galerkin Methods: Theory, Computation and Applications Paola. Antonietti MOX, Dipartimento di Matematica Politecnico di Milano.MO. X MODELLISTICA E CALCOLO SCIENTIICO. MODELING AND SCIENTIIC
More informationSymmetry-Free, p-robust Equilibrated Error Indication for the hp-version of the FEM in Nearly Incompressible Linear Elasticity
Computational Methods in Applied Mathematics Vol. 13 (2013), No. 3, pp. 291 304 c 2013 Institute of Mathematics, NAS of Belarus Doi: 10.1515/cmam-2013-0007 Symmetry-Free, p-robust Equilibrated Error Indication
More informationarxiv: 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 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 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 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 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 informationAdaptive 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 informationETNA 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 informationA 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 informationAMS subject classifications. Primary, 65N15, 65N30, 76D07; Secondary, 35B45, 35J50
A SIMPLE FINITE ELEMENT METHOD FOR THE STOKES EQUATIONS LIN MU AND XIU YE Abstract. The goal of this paper is to introduce a simple finite element method to solve the Stokes equations. This method is in
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 informationA POSTERIORI ERROR ESTIMATES INCLUDING ALGEBRAIC ERROR: COMPUTABLE UPPER BOUNDS AND STOPPING CRITERIA FOR ITERATIVE SOLVERS
A POSTERIORI ERROR ESTIMATES INCLUDING ALGEBRAIC ERROR: COMPUTABLE UPPER BOUNDS AND STOPPING CRITERIA FOR ITERATIVE SOLVERS PAVEL JIRÁNEK, ZDENĚK STRAKOŠ, AND MARTIN VOHRALÍK Abstract. We consider the
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 informationMaximum-norm a posteriori estimates for discontinuous Galerkin methods
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
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 informationGuaranteed and robust a posteriori bounds for Laplace eigenvalues and eigenvectors: conforming approximations
Guaranteed and robust a posteriori bounds for Laplace eigenvalues and eigenvectors: conforming approximations Eric Cancès, Geneviève Dusson, Yvon Maday, Benjamin Stamm, Martin Vohralík To cite this version:
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 informationAxioms of Adaptivity (AoA) in Lecture 1 (sufficient for optimal convergence rates)
Axioms of Adaptivity (AoA) in Lecture 1 (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 informationAxioms of Adaptivity
Axioms of Adaptivity Carsten Carstensen Humboldt-Universität zu Berlin Open-Access Reference: C-Feischl-Page-Praetorius: Axioms of Adaptivity. Computer & Mathematics with Applications 67 (2014) 1195 1253
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 informationA 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 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 informationRESIDUAL BASED ERROR ESTIMATES FOR THE SPACE-TIME DISCONTINUOUS GALERKIN METHOD APPLIED TO NONLINEAR HYPERBOLIC EQUATIONS
Proceedings of ALGORITMY 2016 pp. 113 124 RESIDUAL BASED ERROR ESTIMATES FOR THE SPACE-TIME DISCONTINUOUS GALERKIN METHOD APPLIED TO NONLINEAR HYPERBOLIC EQUATIONS VÍT DOLEJŠÍ AND FILIP ROSKOVEC Abstract.
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 informationA posteriori error estimates in FEEC for the de Rham complex
A posteriori error estimates in FEEC for the de Rham complex Alan Demlow Texas A&M University joint work with Anil Hirani University of Illinois Urbana-Champaign Partially supported by NSF DMS-1016094
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 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 informationBasic Principles of Weak Galerkin Finite Element Methods for PDEs
Basic Principles of Weak Galerkin Finite Element Methods for PDEs Junping Wang Computational Mathematics Division of Mathematical Sciences National Science Foundation Arlington, VA 22230 Polytopal Element
More informationEdge finite element approximation of Maxwell s equations with low regularity solutions
Edge finite element approximation of Maxwell s equations with low regularity solutions and Jean-Luc Guermond Université Paris-Est, CERMICS,, Paris Vienna, 18 July 2018 Outline Maxwell s equations Edge
More informationCell Conservative Flux Recovery and A Posteriori Error Estimate of Vertex-Centered Finite Volume Methods
Advances in Applied Mathematics and Mechanics Adv. Appl. Math. Mech., Vol. 5, No. 5, pp. 705-727 DOI: 10.4208/aamm.12-m1279 October 2013 Cell Conservative Flux Recovery and A Posteriori Error Estimate
More informationA posteriori error estimates for a Maxwell type problem
Russ. J. Numer. Anal. Math. Modelling, Vol. 24, No. 5, pp. 395 408 (2009) DOI 0.55/ RJNAMM.2009.025 c de Gruyter 2009 A posteriori error estimates for a Maxwell type problem I. ANJAM, O. MALI, A. MUZALEVSKY,
More informationGoal. 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 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 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 informationTowards pressure-robust mixed methods for the incompressible Navier Stokes equations
Mohrenstrasse 39 10117 Berlin Germany Tel. +49 30 20372 0 www.wias-berlin.de 2016-05-05, Pittsburgh Weierstrass Institute for Applied Analysis and Stochastics Towards pressure-robust mixed methods for
More informationAn Iterative Substructuring Method for Mortar Nonconforming Discretization of a Fourth-Order Elliptic Problem in two dimensions
An Iterative Substructuring Method for Mortar Nonconforming Discretization of a Fourth-Order Elliptic Problem in two dimensions Leszek Marcinkowski Department of Mathematics, Warsaw University, Banacha
More informationComplementarity 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 informationNumerical 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 informationMIXED FINITE ELEMENT METHODS FOR PROBLEMS WITH ROBIN BOUNDARY CONDITIONS
MIXED FINITE ELEMENT METHODS FOR PROBLEMS WITH ROBIN BOUNDARY CONDITIONS JUHO KÖNNÖ, DOMINIK SCHÖTZAU, AND ROLF STENBERG Abstract. We derive new a-priori and a-posteriori error estimates for mixed nite
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 informationFind (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 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 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 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 Adaptive Mixed Finite Element Method using the Lagrange Multiplier Technique
An Adaptive Mixed Finite Element Method using the Lagrange Multiplier Technique by Michael Gagnon A Project Report Submitted to the Faculty of the WORCESTER POLYTECHNIC INSTITUTE In partial fulfillment
More informationA WEAK GALERKIN MIXED FINITE ELEMENT METHOD FOR BIHARMONIC EQUATIONS
A WEAK GALERKIN MIXED FINITE ELEMENT METHOD FOR BIHARMONIC EQUATIONS LIN MU, JUNPING WANG, YANQIU WANG, AND XIU YE Abstract. This article introduces and analyzes a weak Galerkin mixed finite element method
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 informationTrefftz-discontinuous Galerkin methods for time-harmonic wave problems
Trefftz-discontinuous Galerkin methods for time-harmonic wave problems Ilaria Perugia Dipartimento di Matematica - Università di Pavia (Italy) http://www-dimat.unipv.it/perugia Joint work with Ralf Hiptmair,
More informationMultigrid Methods for Saddle Point Problems
Multigrid Methods for Saddle Point Problems Susanne C. Brenner Department of Mathematics and Center for Computation & Technology Louisiana State University Advances in Mathematics of Finite Elements (In
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 informationAnalysis 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 informationDiscontinuous Petrov-Galerkin Methods
Discontinuous Petrov-Galerkin Methods Friederike Hellwig 1st CENTRAL School on Analysis and Numerics for Partial Differential Equations, November 12, 2015 Motivation discontinuous Petrov-Galerkin (dpg)
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 informationRecovery-Based a Posteriori Error Estimators for Interface Problems: Mixed and Nonconforming Elements
Recovery-Based a Posteriori Error Estimators for Interface Problems: Mixed and Nonconforming Elements Zhiqiang Cai Shun Zhang Department of Mathematics Purdue University Finite Element Circus, Fall 2008,
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 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 informationHybridized DG methods
Hybridized DG methods University of Florida (Banff International Research Station, November 2007.) Collaborators: Bernardo Cockburn University of Minnesota Raytcho Lazarov Texas A&M University Thanks:
More informationQuasi-Optimal Cardinality of AFEM Driven by Nonresidual Estimators
IMA Journal of Numerical Analysis (010) Page 1 of 8 doi:10.1093/imanum/drnxxx Quasi-Optimal Cardinality of AFEM Driven by Nonresidual Estimators J. MANUEL CASCÓN Departamento de Economía e Historia Económica,
More informationNon-Conforming Finite Element Methods for Nonmatching Grids in Three Dimensions
Non-Conforming Finite Element Methods for Nonmatching Grids in Three Dimensions Wayne McGee and Padmanabhan Seshaiyer Texas Tech University, Mathematics and Statistics (padhu@math.ttu.edu) Summary. In
More informationAn optimal adaptive finite element method. Rob Stevenson Korteweg-de Vries Institute for Mathematics Faculty of Science University of Amsterdam
An optimal adaptive finite element method Rob Stevenson Korteweg-de Vries Institute for Mathematics Faculty of Science University of Amsterdam Contents model problem + (A)FEM newest vertex bisection convergence
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 informationarxiv: v1 [math.na] 11 Jul 2011
Multigrid Preconditioner for Nonconforming Discretization of Elliptic Problems with Jump Coefficients arxiv:07.260v [math.na] Jul 20 Blanca Ayuso De Dios, Michael Holst 2, Yunrong Zhu 2, and Ludmil Zikatanov
More informationarxiv: v2 [math.na] 2 Nov 2018
DISCREIZION-CCURE SOPPING CRIERION FOR IERIVE SOLVERS ON FINIE ELEMEN PPROXIMION ZHIQING CI, SHUHO CO, ND ROB FLGOU arxiv:161108967v2 [mathn] 2 Nov 2018 bstract his paper introduces a localized duality-type
More informationNumerische Mathematik
Numer. Math. (2013) :425 459 DOI 10.1007/s00211-012-0494-4 Numerische Mathematik Effective postprocessing for equilibration a posteriori error estimators C. Carstensen C. Merdon Received: 22 June 2011
More informationChapter 2. The constitutive relation error method. for linear problems
2. The constitutive relation error method for linear problems Chapter 2 The constitutive relation error method for linear problems 2.1 Introduction The reference problem considered here is the linear problem
More informationOn angle conditions in the finite element method. Institute of Mathematics, Academy of Sciences Prague, Czech Republic
On angle conditions in the finite element method Michal Křížek Institute of Mathematics, Academy of Sciences Prague, Czech Republic Joint work with Jan Brandts (University of Amsterdam), Antti Hannukainen
More informationA POSTERIORI ERROR ESTIMATES INCLUDING ALGEBRAIC ERROR AND STOPPING CRITERIA FOR ITERATIVE SOLVERS
A POSTERIORI ERROR ESTIMATES INCLUDING ALGEBRAIC ERROR AND STOPPING CRITERIA FOR ITERATIVE SOLVERS PAVEL JIRÁNEK, ZDENĚK STRAKOŠ, AND MARTIN VOHRALÍK Abstract. For the finite volume discretization of a
More informationBUBBLE STABILIZED DISCONTINUOUS GALERKIN METHOD FOR PARABOLIC AND ELLIPTIC PROBLEMS
BUBBLE STABILIZED DISCONTINUOUS GALERKIN METHOD FOR PARABOLIC AND ELLIPTIC PROBLEMS ERIK BURMAN AND BENJAMIN STAMM Abstract. In this paper we give an analysis of a bubble stabilized discontinuous Galerkin
More informationMultilevel Preconditioning of Graph-Laplacians: Polynomial Approximation of the Pivot Blocks Inverses
Multilevel Preconditioning of Graph-Laplacians: Polynomial Approximation of the Pivot Blocks Inverses P. Boyanova 1, I. Georgiev 34, S. Margenov, L. Zikatanov 5 1 Uppsala University, Box 337, 751 05 Uppsala,
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 informationFINITE ELEMENT APPROXIMATION OF STOKES-LIKE SYSTEMS WITH IMPLICIT CONSTITUTIVE RELATION
Proceedings of ALGORITMY pp. 9 3 FINITE ELEMENT APPROXIMATION OF STOKES-LIKE SYSTEMS WITH IMPLICIT CONSTITUTIVE RELATION JAN STEBEL Abstract. The paper deals with the numerical simulations of steady flows
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 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 informationA posteriori error estimates for space-time domain decomposition method for two-phase flow problem
A posteriori error estimates for space-time domain decomposition method for two-phase flow problem Sarah Ali Hassan, Elyes Ahmed, Caroline Japhet, Michel Kern, Martin Vohralík INRIA Paris & ENPC (project-team
More informationc 2017 Society for Industrial and Applied Mathematics
SIAM J. NUMER. ANAL. Vol. 55, No. 6, pp. 2811 2834 c 2017 Society for Industrial and Applied Matematics GUARANTEED, LOCALLY SPACE-TIME EFFICIENT, AND POLYNOMIAL-DEGREE ROBUST A POSTERIORI ERROR ESTIMATES
More informationSpace-time Finite Element Methods for Parabolic Evolution Problems
Space-time Finite Element Methods for Parabolic Evolution Problems with Variable Coefficients Ulrich Langer, Martin Neumüller, Andreas Schafelner Johannes Kepler University, Linz Doctoral Program Computational
More informationWEAK GALERKIN FINITE ELEMENT METHODS ON POLYTOPAL MESHES
INERNAIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume 12, Number 1, Pages 31 53 c 2015 Institute for Scientific Computing and Information WEAK GALERKIN FINIE ELEMEN MEHODS ON POLYOPAL MESHES LIN
More informationA posteriori error analysis of an augmented mixed finite element method for Darcy flow
A posteriori error analysis of an augmented mixed finite element method for Darcy flow Tomás P. Barrios, J. Manuel Cascón and María González Abstract We develop an a posteriori error analysis of residual
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 informationhp-adaptation driven by polynomial-degree-robust a posteriori error estimates for elliptic problems
p-adaptation driven by polynomial-degree-robust a posteriori error estimates for elliptic problems Vít Dolejší, Alexandre Ern, Martin Voralík To cite tis version: Vít Dolejší, Alexandre Ern, Martin Voralík.
More informationAdaptive tree approximation with finite elements
Adaptive tree approximation with finite elements Andreas Veeser Università degli Studi di Milano (Italy) July 2015 / Cimpa School / Mumbai Outline 1 Basic notions in constructive approximation 2 Tree approximation
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 informationWeak Galerkin Finite Element Methods and Applications
Weak Galerkin Finite Element Methods and Applications Lin Mu mul1@ornl.gov Computational and Applied Mathematics Computationa Science and Mathematics Division Oak Ridge National Laboratory Georgia Institute
More information