Unified A Posteriori Error Control for all Nonstandard Finite Elements 1
|
|
- Philip Lynch
- 5 years ago
- Views:
Transcription
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 we know of Guidelines for Applicants M. Eigel (Humboldt) Unified ErrorThese Analysis guidelines should help you EFEF, as you Warwick think about 2010 applying 1 / to 22the
2 Model Example Flux or stress field p in equilibrium equation g + div p = 0 is approximated by piecewise constant p l and yields equilibrium residual Res(v) := (g v p l : D l v) dx for v V := H0 1 (Ω; R m ) Ω Endow V with norm v V := Dv L 2 (Ω) s.t. V H 1 (Ω) and Res(v) div(p p l ) V = Res V := sup v 0 v V Estimation by edge residuals η E := E (p l T+ ν T+ + p l T ν T ) on each interior edge E = T + T with outer unit normals ν T± ( ) 1/2 p l η E := η E 2 T+ E E ν T+ E ν T p l T M. Eigel (Humboldt) Unified Error Analysis EFEF, Warwick / 22
3 Estimation of Equilibrium Error Up to data oscillation osc(g, ), edge estimator η E is reliable & efficient Res V η E ± osc(g, elements edges) Remark: Edge residual estimator is equivalent to many others (cf. Ainsworth/Oden, Babuška/Strouboulis, Verfürth) Example: For triangulation of Ω R 2 and first-order conforming or nonconforming FEM p l := D l u l and V c l := P 1(T l ; R m ) V ker Res M. Eigel (Humboldt) Unified Error Analysis EFEF, Warwick / 22
4 Schemes & Applications CR Wilson Han NR (M) NR (A) CNR DSSY Laplace CR KS Han NR (M) NR (A) HMS CJY Stokes BS KS Zhang Ming LLS HMS Elasticity cf. C. Carstensen, Hu Jun, A. Orlando M. Eigel (Humboldt) Unified Error Analysis EFEF, Warwick / 22
5 Unified Analysis of... applications: Laplace, Stokes, Navier-Lamé, Maxwell equations... schemes: (all?!) conforming, nonconforming, tri/quad, mixed, mortar elements, dg, hanging nodes... M. Eigel (Humboldt) Unified Error Analysis EFEF, Warwick / 22
6 Goals of Unified Analysis Generalise analysis to cover many different discretisation schemes and applications in one framework Reduce repetition of similar mathematical arguments and focus on specific properties/difficulties No optimal constants but common point of departure and guiding principles M. Eigel (Humboldt) Unified Error Analysis EFEF, Warwick / 22
7 A Unified Approach Generic Approach For each Application: Verify error residual For each Scheme: Determine discrete space V l ker(residual) and design computable lower/upper bounds of residual Topics Mixed Setting for Unified A Posteriori Error Control Unified Equilibrium Estimator Analysis of Consistency Residual Applications: Poisson, Lamé, Stokes,... M. Eigel (Humboldt) Unified Error Analysis EFEF, Warwick / 22
8 Mixed Setting for Unified Analysis Abstract problem formulation: Given X l Y l X Y, A : X Y (X Y ) l := l X + l Y (X Y ), find (x, y) X Y s.t. (PM) A(x, y)(ξ, η) = l(ξ, η) for all (ξ, η) X Y Given a (X X ), Λ : Y X, b (X Y ) with b(, v) := a(λv, ) c (Y Y ), A(x, y) := a(x, ) b(, y) + b(x, ) + c(y, ) (PM) then reads a(x, ξ) b(ξ, y) = l X (ξ) b(x, η) + c(y, η) = l Y (η) for all ξ X for all η Y M. Eigel (Humboldt) Unified Error Analysis EFEF, Warwick / 22
9 Mixed Setting for Unified Analysis Abstract problem formulation: Given X l Y l X Y, A : X Y (X Y ) l := l X + l Y (X Y ), find (x, y) X Y s.t. (PM) A(x, y)(ξ, η) = l(ξ, η) for all (ξ, η) X Y Given a (X X ), Λ : Y X, b (X Y ) with b(, v) := a(λv, ) c (Y Y ), A(x, y) := a(x, ) b(, y) + b(x, ) + c(y, ) (PM) then reads a(x, ξ) b(ξ, y) = l X (ξ) b(x, η) + c(y, η) = l Y (η) for all ξ X for all η Y M. Eigel (Humboldt) Unified Error Analysis EFEF, Warwick / 22
10 Errors & Residuals in Unified Analysis Given approx. (x l, ỹ l ) X l Y to (x, y), define Res := Res X + Res Y (consistency) (equilibrium) Res X := l X a(x l, ) + b(, ỹ l ) = l X a(x l Λỹ l ) Res Y := l Y b(x l, ) c(ỹ l, ) Remarks: ỹ l Y close to y l, not necessarily discrete Res X involves piecewise gradient D l (y l ỹ l ) of y l ỹ l / Y Since A isomorphism, error residual, i.e. x x l X + y ỹ l Y Res X X + Res Y Y M. Eigel (Humboldt) Unified Error Analysis EFEF, Warwick / 22
11 Errors & Residuals in Unified Analysis Given approx. (x l, ỹ l ) X l Y to (x, y), define Res := Res X + Res Y (consistency) (equilibrium) Res X := l X a(x l, ) + b(, ỹ l ) = l X a(x l Λỹ l ) Res Y := l Y b(x l, ) c(ỹ l, ) Remarks: ỹ l Y close to y l, not necessarily discrete Res X involves piecewise gradient D l (y l ỹ l ) of y l ỹ l / Y Since A isomorphism, error residual, i.e. x x l X + y ỹ l Y Res X X + Res Y Y M. Eigel (Humboldt) Unified Error Analysis EFEF, Warwick / 22
12 [Preparation] Generic Equilibrium Error Analysis Two discrete spaces V c l V V and V nc l H 1 (T l ; R m ) with J V c l Π V nc l (H1) H 1 -stable Clément-type operator J : V Vl c into (conforming) subspace Vl c V with first-order approximation property (H2) Vl c and V nc l piecewise smooth w.r.t. shape-regular T l (H3) For p l L 2 (Ω; R m n ) Π : Vl c V l nc s.t. v l Vl c T T l D(Πv l ) L 2 (T ) Dv l L 2 (ω T ) v l dx = Πv l dx T T p l : D l v l dx = p l : D l (Πv l ) dx Ω Ω M. Eigel (Humboldt) Unified Error Analysis EFEF, Warwick / 22
13 [Result] Generic Equilibrium Error Analysis Thm.[CEHL10 + ]: Suppose R T L 2 (T l ), R E L 2 ( E l ) and Res : Vl nc + V R reads Res(v) := R T v dx + S R E v ds El Suppose (H1)-(H3) and Vl nc ker Res Then η l := h E R E 2 L 2 (E) E El is reliable and efficient in the sense Ω 1/2 = h 1/2 E R E L 2 ( S E l ) η l osc(r T, T l ) osc(r E, E l ) Res V η l + osc(r T, {ω z : z K l }) M. Eigel (Humboldt) Unified Error Analysis EFEF, Warwick / 22
14 [Result] Generic Equilibrium Error Analysis Thm.[CEHL10 + ]: Suppose R T L 2 (T l ), R E L 2 ( E l ) and Res : Vl nc + V R reads Res(v) := R T v dx + S R E v ds El Suppose (H1)-(H3) and Vl nc ker Res Then η l := h E R E 2 L 2 (E) E El is reliable and efficient in the sense Ω 1/2 = h 1/2 E R E L 2 ( S E l ) η l osc(r T, T l ) osc(r E, E l ) Res V η l + osc(r T, {ω z : z K l }) M. Eigel (Humboldt) Unified Error Analysis EFEF, Warwick / 22
15 Analysis of Consistency Error Thm.[CEHL10 + ]: Given X = L 2 (Ω), x l = D l y l X l = P 0 (T l ; R n ), Y = H 1 0 (Ω), y l Y l = CR 1 (T l ) and µ l := min η Y x l Dη L 2 (Ω) µ l min η l P 1 (T l ) Y x l D l η l L 2 (Ω) min η l P 1 (T l ) Y h 1 T (y l η l ) L 2 (Ω) h 1 T (y l Sy l ) L 2 (Ω) h 1/2 E [y l ] L 2 ( S E l ) h 1/2 E [D ly l τ E ] L 2 ( S E l ) M. Eigel (Humboldt) Unified Error Analysis EFEF, Warwick / 22
16 Analysis of Consistency Error Thm.[CEHL10 + ]: Given X = L 2 (Ω), x l = D l y l X l = P 0 (T l ; R n ), Y = H 1 0 (Ω), y l Y l = CR 1 (T l ) and µ l := min η Y x l Dη L 2 (Ω) µ l min η l P 1 (T l ) Y x l D l η l L 2 (Ω) min η l P 1 (T l ) Y h 1 T (y l η l ) L 2 (Ω) h 1 T (y l Sy l ) L 2 (Ω) h 1/2 E [y l ] L 2 ( S E l ) h 1/2 E [D ly l τ E ] L 2 ( S E l ) M. Eigel (Humboldt) Unified Error Analysis EFEF, Warwick / 22
17 [Application] Laplace Equation Poisson model problem: u = g in Ω and u = 0 on Ω leads to primal mixed formulation with a(p, q) := p q dx, Λu := D l u, b(q, u) = Ω Ω D l u q dx Given g L 2 (Ω), find (p, u) X Y := L 2 (Ω; R n ) H0 1 (Ω) s.t. A(p, u)(q, v) = (g, v) L 2 (Ω) for all (q, v) X Y Since A isomorphism p p l L 2 + u ũ l H 1 Res X X + Res Y Y M. Eigel (Humboldt) Unified Error Analysis EFEF, Warwick / 22
18 [Application (cont.)] Laplace Equation Treatment of different schemes in the unified framework: Conforming FEM Discrete solution u l = ũ l Y defines discrete flux p l := Λu l = Du l. Consistency error vanishes, equilibrium residual treated as suggested previously. Nonconforming FEM Discrete solution u l defines discrete flux p l := D l u l. Same analysis for equilibrium residual while consistency residual p l Dũ l L 2 (Ω) can be bounded as before. Mixed FEM Discrete solution is discrete flux p l and u l / Y is Lagrange multiplier. Equilibrium residual reduces to osc(g, elements), consistency residual minũl H 1 0 (Ω) p l Dũ l L 2 (Ω) can be bounded as in [Carstensen, Math. Comp. 1997]. M. Eigel (Humboldt) Unified Error Analysis EFEF, Warwick / 22
19 [Application] Stokes Equations Let a(u, v) := Ω uv dx, Λu := div u, c(u, v) := Ω 2µε(u) : ε(v) dx l X (p, q) := Ω pq dx and X := L2 0 (Ω; Rn ), Y := H 1 0 (Ω; Rn ) With ε(v) := symdv and deviatoric operator dev σ := σ 1 n (tr(σ))1 the linear operator A : X Y (X Y ) reads A(σ, u)(τ, v) := 1 2µ (dev σ, dev τ) L 2 (Ω) (σ, ε(v)) L 2 (Ω) (τ, ε(u)) L 2 (Ω) A isomorphism... M. Eigel (Humboldt) Unified Error Analysis EFEF, Warwick / 22
20 [Application (cont.)] Stokes Equations Stokes problem: Given g L 2 (Ω; R n ), find (σ, u) X Y s.t. A(σ, u)(τ, v) = (g, v) L 2 (Ω) for all (τ, v) X Y Conforming or nonconforming FEM yield u l and p l with σ = 2µε(u) p1 and σ l = 2µε l (u l ) p l 1 Error control: Given FE solution (σ l, u l ) to (σ, u), then σ σ l L 2 (Ω) min ũ l Y 2µ ε l(u l ũ l ) L 2 (Ω) + Res Y Y + div l u l L 2 (Ω) Examples: all known conforming and nonconforming FEM M. Eigel (Humboldt) Unified Error Analysis EFEF, Warwick / 22
21 [Application] Linear Elasticity Fourth-order elasticity tensor for λ, µ > 0, Cτ := λ tr(τ)1 + 2µτ for all τ R n n Let a(σ, τ) := Ω (C 1 σ) : τ dx, X := L 2 (Ω; R n n sym ), Y := H 1 0 (Ω; Rn ) Λu := Cε(u) and The linear operator A : X Y (X Y ) then reads A(σ, u)(τ, v) := (C 1 σ, τ) L 2 (Ω) (σ, ε(v)) L 2 (Ω) (τ, ε(u)) L 2 (Ω) A isomorphism, λ-independent operator norms of A and A 1 M. Eigel (Humboldt) Unified Error Analysis EFEF, Warwick / 22
22 [Application (cont.)] Linear Elasticity Navier-Lamé problem: Given g L 2 (Ω; R n ), find (σ, u) X Y s.t. A(σ, u)(τ, v) = (g, v) L 2 (Ω) for (τ, v) X Y For conforming or nonconforming FEM, σ = Cε(u) and σ l = Cε l (u l ) Error control: Given FE solution (σ l, u l ) to (σ, u), is robust for λ σ σ l L 2 (Ω) min ũ l Y ε l(u l ũ l ) L 2 (Ω) + Res Y Y Examples: Conforming FEM, Kouhia&Stenberg, PEERS M. Eigel (Humboldt) Unified Error Analysis EFEF, Warwick / 22
23 Conclusions for Unified Analysis Framework for unified a posteriori analysis covers large class of applications & discretisation schemes Similarities are exposed and encountered obstacles/challenges guide development of specific error estimators Approach doesn t strive for most accurate/efficient estimators but rather provides an initial line of attack for as many problems as possible M. Eigel (Humboldt) Unified Error Analysis EFEF, Warwick / 22
24 Wrap-Up of Unified A Posteriori Analysis C. Carstensen: A unifying theory of a posteriori finite element error control (Numer. Math. 2005) C. Carstensen, Jun Hu, A. Orlando: Framework for the a posteriori error analysis of nonconforming finite elements (SINUM 2007) C. Carstensen, Jun Hu: A Unifying Theory of A Posteriori Error Control for Nonconforming FEM (Numer. Math. 2007) C. Carstensen, ME, R.H.W. Hoppe, C. Löbhard: A Unifying Theory of A Posteriori Error Control ( ) M. Eigel (Humboldt) Unified Error Analysis EFEF, Warwick / 22
25 Extensions of Unified Analysis Unified Analysis has also been applied/extended to Maxwell Equations [C. Carstensen, R.H.W. Hoppe 2009] Mortar FEM and higher order FEM dg FEM [C. Carstensen, T. Gudi, M. Jensen 2009] Hanging nodes [C. Carstensen, Jun Hu 2008] M. Eigel (Humboldt) Unified Error Analysis EFEF, Warwick / 22
26 References C. Carstensen, ME, R.H.W. Hoppe, C. Löbhard: Unified A Posteriori Error Control, in preparation C. Carstensen: A unifying theory of a posteriori finite element error control, Numer. Math. 100 (2005) C. Carstensen, Jun Hu, A. Orlando: Framework for the a posteriori error analysis of nonconforming finite elements, SINUM 45, 1 (2007) C. Carstensen, Jun Hu: A Unifying Theory of A Posteriori Error Control for Nonconforming Finite Element Methods, Numer. Math. (2007), DOI /s z C. Carstensen, T. Gudi, M. Jensen: A unifying theory of a posteriori error control for discontinuous Galerkin FEM C. Carstensen, R.H.W. Hoppe: Unified framework for an a posteriori analysis of non-standard finite element approximations of H(curl)-elliptic problems M. Eigel (Humboldt) Unified Error Analysis EFEF, Warwick / 22
A UNIFYING THEORY OF A POSTERIORI FINITE ELEMENT ERROR CONTROL
A UNIFYING THEORY OF A POSTERIORI FINITE ELEMENT ERROR CONTROL C. CARSTENSEN Abstract. Residual-based a posteriori error estimates are derived within a unified setting for lowest-order conforming, nonconforming,
More 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 informationc 2007 Society for Industrial and Applied Mathematics
SIAM J. NUMR. ANAL. Vol. 45, No. 1, pp. 68 82 c 2007 Society for Industrial and Applied Mathematics FRAMWORK FOR TH A POSTRIORI RROR ANALYSIS OF NONCONFORMING FINIT LMNTS CARSTN CARSTNSN, JUN HU, AND ANTONIO
More informationFive Recent Trends in Computational PDEs
Five Recent Trends in Computational PDEs Carsten Carstensen Center Computational Science Adlershof and Department of Mathematics, Humboldt-Universität zu Berlin C. Carstensen (CCSA & HU Berlin) 5 Trends
More informationMedius analysis and comparison results for first-order finite element methods in linear elasticity
IMA Journal of Numerical Analysis Advance Access published November 7, 2014 IMA Journal of Numerical Analysis (2014) Page 1 of 31 doi:10.1093/imanum/dru048 Medius analysis and comparison results for first-order
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 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 informationSpline Element Method for Partial Differential Equations
for Partial Differential Equations Department of Mathematical Sciences Northern Illinois University 2009 Multivariate Splines Summer School, Summer 2009 Outline 1 Why multivariate splines for PDEs? Motivation
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 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 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 informationYongdeok Kim and Seki Kim
J. Korean Math. Soc. 39 (00), No. 3, pp. 363 376 STABLE LOW ORDER NONCONFORMING QUADRILATERAL FINITE ELEMENTS FOR THE STOKES PROBLEM Yongdeok Kim and Seki Kim Abstract. Stability result is obtained for
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 informationAxioms of Adaptivity (AoA) in Lecture 2 (sufficient for optimal convergence rates)
Axioms of Adaptivity (AoA) in Lecture 2 (sufficient for optimal convergence rates) Carsten Carstensen Humboldt-Universität zu Berlin 2018 International Graduate Summer School on Frontiers of Applied and
More informationTHE PATCH RECOVERY FOR FINITE ELEMENT APPROXIMATION OF ELASTICITY PROBLEMS UNDER QUADRILATERAL MESHES. Zhong-Ci Shi and Xuejun Xu.
DISCRETE AND CONTINUOUS Website: http://aimsciences.org DYNAMICAL SYSTEMS SERIES B Volume 9, Number, January 2008 pp. 63 82 THE PATCH RECOVERY FOR FINITE ELEMENT APPROXIMATION OF ELASTICITY PROBLEMS UNDER
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 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 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 informationarxiv: v1 [math.na] 19 Dec 2017
Arnold-Winther Mixed Finite Elements for Stokes Eigenvalue Problems Joscha Gedicke Arbaz Khan arxiv:72.0686v [math.na] 9 Dec 207 Abstract This paper is devoted to study the Arnold-Winther mixed finite
More informationNUMERICAL EXPERIMENTS FOR THE ARNOLD WINTHER MIXED FINITE ELEMENTS FOR THE STOKES PROBLEM
NUMERICAL EXPERIMENTS FOR THE ARNOLD WINTHER MIXED FINITE ELEMENTS FOR THE STOKES PROBLEM CARSTEN CARSTENSEN, JOSCHA GEDICKE, AND EUN-JAE PARK Abstract. The stress-velocity formulation of the stationary
More informationTraces and Duality Lemma
Traces and Duality Lemma Recall the duality lemma with H / ( ) := γ 0 (H ()) defined as the trace space of H () endowed with minimal extension norm; i.e., for w H / ( ) L ( ), w H / ( ) = min{ ŵ H () ŵ
More informationTheory of PDE Homework 2
Theory of PDE Homework 2 Adrienne Sands April 18, 2017 In the following exercises we assume the coefficients of the various PDE are smooth and satisfy the uniform ellipticity condition. R n is always an
More informationNumerische Mathematik
Numer. Math. (2017) 136:1097 1115 DOI 10.1007/s00211-017-0866-x Numerische Mathematik Convergence of natural adaptive least squares finite element methods Carsten Carstensen 1 Eun-Jae Park 2 Philipp Bringmann
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 informationModel Reduction via Discretization : Elasticity and Boltzmann
Model Reduction via Discretization : Elasticity and Boltzmann Joachim Schöberl Computational Mathematics in Engineering Institute for Analysis and Scientific Computing Vienna University of Technology IMA
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 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 informationApplied Numerical Mathematics
Applied Numerical Mathematics 59 009) 119 130 Contents lists available at ScienceDirect Applied Numerical Mathematics www.elsevier.com/locate/apnum Convergence of adaptive finite element methods in computational
More 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 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 informationChapter 1: The Finite Element Method
Chapter 1: The Finite Element Method Michael Hanke Read: Strang, p 428 436 A Model Problem Mathematical Models, Analysis and Simulation, Part Applications: u = fx), < x < 1 u) = u1) = D) axial deformation
More 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 informationMixed Finite Element Methods. Douglas N. Arnold, University of Minnesota The 41st Woudschoten Conference 5 October 2016
Mixed Finite Element Methods Douglas N. Arnold, University of Minnesota The 41st Woudschoten Conference 5 October 2016 Linear elasticity Given the load f : Ω R n, find the displacement u : Ω R n and the
More informationSome remarks on the history and future of averaging techniques in a posteriori finite element error analysis
ZAMM Z. Angew. Math. Mech. 84, No. 1, 3 21 (2004) / DOI 10.1002/amm.200410101 Some remarks on the history and future of averaging techniques in a posteriori finite element error analysis Plenary lecture
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 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 informationA Mixed Nonconforming Finite Element for Linear Elasticity
A Mixed Nonconforming Finite Element for Linear Elasticity Zhiqiang Cai, 1 Xiu Ye 2 1 Department of Mathematics, Purdue University, West Lafayette, Indiana 47907-1395 2 Department of Mathematics and Statistics,
More 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 informationProblems of Corner Singularities
Stuttgart 98 Problems of Corner Singularities Monique DAUGE Institut de Recherche MAthématique de Rennes Problems of Corner Singularities 1 Vertex and edge singularities For a polyhedral domain Ω and a
More informationNumerische Mathematik
Numer. Math. (203) 24:309 335 DOI 0.007/s002-02-053-5 Numerische Mathematik A posteriori error estimates for nonconforming finite element methods for fourth-order problems on rectangles Carsten Carstensen
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 informationA primer on Numerical methods for elasticity
A primer on Numerical methods for elasticity Douglas N. Arnold, University of Minnesota Complex materials: Mathematical models and numerical methods Oslo, June 10 12, 2015 One has to resort to the indignity
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 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 informationWeak Convergence Methods for Energy Minimization
Weak Convergence Methods for Energy Minimization Bo Li Department of Mathematics University of California, San Diego E-mail: bli@math.ucsd.edu June 3, 2007 Introduction This compact set of notes present
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 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 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 informationON SINGULAR PERTURBATION OF THE STOKES PROBLEM
NUMERICAL ANALYSIS AND MATHEMATICAL MODELLING BANACH CENTER PUBLICATIONS, VOLUME 9 INSTITUTE OF MATHEMATICS POLISH ACADEMY OF SCIENCES WARSZAWA 994 ON SINGULAR PERTURBATION OF THE STOKES PROBLEM G. M.
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 informationOverview. A Posteriori Error Estimates for the Biharmonic Equation. Variational Formulation and Discretization. The Biharmonic Equation
Overview A Posteriori rror stimates for the Biharmonic quation R Verfürth Fakultät für Mathematik Ruhr-Universität Bochum wwwruhr-uni-bochumde/num1 Milan / February 11th, 013 The Biharmonic quation Summary
More informationLECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES)
LECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES) RAYTCHO LAZAROV 1 Notations and Basic Functional Spaces Scalar function in R d, d 1 will be denoted by u,
More informationA note on discontinuous Galerkin divergence-free solutions of the Navier-Stokes equations
A note on discontinuous Galerkin divergence-free solutions of the Navier-Stokes equations Bernardo Cockburn Guido anschat Dominik Schötzau June 1, 2007 Journal of Scientific Computing, Vol. 31, 2007, pp.
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 informationAposteriorierrorestimatesinFEEC for the de Rham complex
AposteriorierrorestimatesinFEEC 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 and a
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 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 informationFinite Element Exterior Calculus and Applications Part V
Finite Element Exterior Calculus and Applications Part V Douglas N. Arnold, University of Minnesota Peking University/BICMR August 15 18, 2015 De Rham complex 0 H 1 (Ω) grad H(curl, Ω) curl H(div, Ω) div
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 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 informationApproximation of fluid-structure interaction problems with Lagrange multiplier
Approximation of fluid-structure interaction problems with Lagrange multiplier Daniele Boffi Dipartimento di Matematica F. Casorati, Università di Pavia http://www-dimat.unipv.it/boffi May 30, 2016 Outline
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 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 informationA STOKES INTERFACE PROBLEM: STABILITY, FINITE ELEMENT ANALYSIS AND A ROBUST SOLVER
European Congress on Computational Methods in Applied Sciences and Engineering ECCOMAS 2004 P. Neittaanmäki, T. Rossi, K. Majava, and O. Pironneau (eds.) O. Nevanlinna and R. Rannacher (assoc. eds.) Jyväskylä,
More informationu( x) = g( y) ds y ( 1 ) U solves u = 0 in U; u = 0 on U. ( 3)
M ath 5 2 7 Fall 2 0 0 9 L ecture 4 ( S ep. 6, 2 0 0 9 ) Properties and Estimates of Laplace s and Poisson s Equations In our last lecture we derived the formulas for the solutions of Poisson s equation
More informationThree remarks on anisotropic finite elements
Three remarks on anisotropic finite elements Thomas Apel Universität der Bundeswehr München Workshop Numerical Analysis for Singularly Perturbed Problems dedicated to the 60th birthday of Martin Stynes
More informationLecture Note III: Least-Squares Method
Lecture Note III: Least-Squares Method Zhiqiang Cai October 4, 004 In this chapter, we shall present least-squares methods for second-order scalar partial differential equations, elastic equations of solids,
More 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 informationRobust Monolithic - Multigrid FEM Solver for Three Fields Formulation Rising from non-newtonian Flow Problems
Robust Monolithic - Multigrid FEM Solver for Three Fields Formulation Rising from non-newtonian Flow Problems M. Aaqib Afaq Institute for Applied Mathematics and Numerics (LSIII) TU Dortmund 13 July 2017
More informationA UNIFYING THEORY OF A POSTERIORI ERROR CONTROL FOR NONCONFORMING FINITE ELEMENT METHODS
A UNIFYING THORY OF A POSTRIORI RROR CONTROL FOR NONCONFORMING FINIT LMNT MTHODS C. CARSTNSN AND JUN HU Abstract. Residual-based a posteriori error estimates were derived witin one unifying framework for
More informationSome Remarks on the Reissner Mindlin Plate Model
Some Remarks on the Reissner Mindlin Plate Model Alexandre L. Madureira LNCC Brazil Talk at the Workshop of Numerical Analysis of PDEs LNCC February 10, 2003 1 Outline The 3D Problem and its Modeling Full
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 informationMIXED FINITE ELEMENT METHODS FOR LINEAR VISCOELASTICITY USING WEAK SYMMETRY
MIXED FINITE ELEMENT METHODS FOR LINEAR VISCOELASTICITY USING WEAK SYMMETRY MARIE E. ROGNES AND RAGNAR WINTHER Abstract. Small deformations of a viscoelastic body are considered through the linear Maxwell
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 informationNONCONFORMING MIXED ELEMENTS FOR ELASTICITY
Mathematical Models and Methods in Applied Sciences Vol. 13, No. 3 (2003) 295 307 c World Scientific Publishing Company NONCONFORMING MIXED ELEMENTS FOR ELASTICITY DOUGLAS N. ARNOLD Institute for Mathematics
More informationNumerische Mathematik
Numer. Math. DOI 10.1007/s00211-015-0738-1 Numerische Mathematik A natural nonconforming FEM for the Bingham flow problem is quasi-optimal C. Carstensen 1,2 B. D. Reddy 3 M. Schedensack 4 Received: 30
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 informationHamburger Beiträge zur Angewandten Mathematik
Hamburger Beiträge zur Angewandten Mathematik Numerical analysis of a control and state constrained elliptic control problem with piecewise constant control approximations Klaus Deckelnick and Michael
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 informationLecture 9 Approximations of Laplace s Equation, Finite Element Method. Mathématiques appliquées (MATH0504-1) B. Dewals, C.
Lecture 9 Approximations of Laplace s Equation, Finite Element Method Mathématiques appliquées (MATH54-1) B. Dewals, C. Geuzaine V1.2 23/11/218 1 Learning objectives of this lecture Apply the finite difference
More informationA posteriori error estimates for Maxwell Equations
www.oeaw.ac.at A posteriori error estimates for Maxwell Equations J. Schöberl RICAM-Report 2005-10 www.ricam.oeaw.ac.at A POSTERIORI ERROR ESTIMATES FOR MAXWELL EQUATIONS JOACHIM SCHÖBERL Abstract. Maxwell
More informationROBUST RESIDUAL-BASED A POSTERIORI ARNOLD-WINTHER MIXED FINITE ELEMENT ANALYSIS IN ELASTICITY
ROBUS RSIDUAL-BASD A POSRIORI ARNOLD-WINHR MIXD FINI LMN ANALYSIS IN LASICIY CARSN CARSNSN AND JOSCHA GDICK Abstract. his paper presents a residual-based a posteriori error estimator for the Arnold-Winther
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 informationABHELSINKI UNIVERSITY OF TECHNOLOGY
ABHELSINKI UNIVERSITY OF TECHNOLOGY TECHNISCHE UNIVERSITÄT HELSINKI UNIVERSITE DE TECHNOLOGIE D HELSINKI A posteriori error analysis for the Morley plate element Jarkko Niiranen Department of Structural
More informationConservation Laws & Applications
Rocky Mountain Mathematics Consortium Summer School Conservation Laws & Applications Lecture V: Discontinuous Galerkin Methods James A. Rossmanith Department of Mathematics University of Wisconsin Madison
More informationA LAGRANGE MULTIPLIER METHOD FOR ELLIPTIC INTERFACE PROBLEMS USING NON-MATCHING MESHES
A LAGRANGE MULTIPLIER METHOD FOR ELLIPTIC INTERFACE PROBLEMS USING NON-MATCHING MESHES P. HANSBO Department of Applied Mechanics, Chalmers University of Technology, S-4 96 Göteborg, Sweden E-mail: hansbo@solid.chalmers.se
More informationDivergence-conforming multigrid methods for incompressible flow problems
Divergence-conforming multigrid methods for incompressible flow problems Guido Kanschat IWR, Universität Heidelberg Prague-Heidelberg-Workshop April 28th, 2015 G. Kanschat (IWR, Uni HD) Hdiv-DG Práha,
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 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 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 informationSupraconvergence of a Non-Uniform Discretisation for an Elliptic Third-Kind Boundary-Value Problem with Mixed Derivatives
Supraconvergence of a Non-Uniform Discretisation for an Elliptic Third-Kind Boundary-Value Problem with Mixed Derivatives Etienne Emmrich Technische Universität Berlin, Institut für Mathematik, Straße
More 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 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 informationALL FIRST-ORDER AVERAGING TECHNIQUES FOR A POSTERIORI FINITE ELEMENT ERROR CONTROL ON UNSTRUCTURED GRIDS ARE EFFICIENT AND RELIABLE
MATHEMATICS OF COMPUTATION Volume 73, Number 247, Pages 1153 1165 S 0025-5718(03)01580-1 Article electronically published on August 12, 2003 ALL FIRST-ORDER AVERAGING TECHNIQUES FOR A POSTERIORI FINITE
More informationFEniCS Course. Lecture 0: Introduction to FEM. Contributors Anders Logg, Kent-Andre Mardal
FEniCS Course Lecture 0: Introduction to FEM Contributors Anders Logg, Kent-Andre Mardal 1 / 46 What is FEM? The finite element method is a framework and a recipe for discretization of mathematical problems
More informationLaplace s Equation. Chapter Mean Value Formulas
Chapter 1 Laplace s Equation Let be an open set in R n. A function u C 2 () is called harmonic in if it satisfies Laplace s equation n (1.1) u := D ii u = 0 in. i=1 A function u C 2 () is called subharmonic
More informationProjection Methods for Rotating Flow
Projection Methods for Rotating Flow Daniel Arndt Gert Lube Georg-August-Universität Göttingen Institute for Numerical and Applied Mathematics IACM - ECCOMAS 2014 Computational Modeling of Turbulent and
More informationPhD dissertation defense
Isogeometric mortar methods with applications in contact mechanics PhD dissertation defense Brivadis Ericka Supervisor: Annalisa Buffa Doctor of Philosophy in Computational Mechanics and Advanced Materials,
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 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