An optimal adaptive finite element method. Rob Stevenson Korteweg-de Vries Institute for Mathematics Faculty of Science University of Amsterdam
|
|
- Stuart French
- 5 years ago
- Views:
Transcription
1 An optimal adaptive finite element method Rob Stevenson Korteweg-de Vries Institute for Mathematics Faculty of Science University of Amsterdam
2 Contents model problem + (A)FEM newest vertex bisection convergence of AFEM quasi-optimal computational complexity 1/16
3 Optimal adaptive finite element methods Model problem: Poisson, 2D, newest vertex bisection. [Generalizations: A with A symm. pos. def. piecewise constant, any space dimension n, red-refinements.] Given f H 1 (Ω), find u H 1 0(Ω), a(u, v) := := a(, ) 1 2. Ω u v = Ω fv =: f(v) (v H 1 0(Ω)). τ conforming part. of Ω into triangles, V τ = H 1 0(Ω) T τ P d 1(T ). Find u τ V τ, a(u τ, v τ ) = f(v τ ) (v τ V τ ). AFEM: GALSOLVE, compute a post. error est, MARK, REFINE. 2/16
4 Initial, conforming partition τ 0 Newest vertex bisection τ := REFINE[τ, M] : Bisect all T M τ a few times. Complete. Th 1 (Binev, Dahmen, DeVore 04). With suitable assignment of newest vertices in initial mesh, #τ K #τ 0 K 1 i=0 #M i (unif. in K). (can be generalized to any space dimension [St. 08]) 3/16
5 Perspective: Approximation classes A s := {u H 1 0(Ω) : u A s := sup N N s inf {τ P:#τ #τ 0 N} u u τ H 1 < }, [i.e. u u τ H 1 ε generally requires #τ #τ 0 ε 1/s u 1/s A s. ] When P unif. refs, then u A (d 1)/n = u H d (Ω), i.e., α u L 2 (Ω), α d With P all partitions created by newest vertex bisection: u A (d 1)/n = u B d q(l p (Ω)), any p > ( d 1 n ) 1 i.e., α u L p (Ω) α d ([Binev, Dahmen, DeVore, Petrushev 02]) Regul. th: For suff sm f, n = 2, u in such a Besov space for any d ([Dahlke, DeVore 97]). 4/16
6 A posteriori error estimator η T (f, u τ ) := diam(t ) 2 f + u τ 2 L 2 (T ) + diam(t ) u τ n T 2 L 2 ( T ), Th 2 (Babu ska, Rheinboldt 78; Verfürth 96). u u τ C 1 E(τ, f, u τ ) := C 1 [ η T (f, u τ )] 1 2. T τ Proof. u u τ = sup v H 1 0 (Ω) a(u u τ,v) v. a(u u τ, v) = a(u u τ, v v τ ) = = T τ T (f + u τ )(v v τ ) Ω T f(v v τ ) a(u τ, v v τ ) u τ n(v v τ ). Th 3 (St. 06, Cascon,Kreuzer,Nochetto,Siebert 07). σ τ, R τ σ := {T τ : T σ}, u σ u τ C 1 [ η T (f, u τ )] 1 T R τ σ Note #R τ σ #σ #τ. 5/16 2.
7 A posteriori error estimator We call σ τ a full refinement with respect to T τ, when T and all its neighbours in τ, as well as all faces of T contain a vertex of σ in their interiors. T τσ Th 4 (Morin, Nochetto, Siebert 00). Let f T τ P d 2(T ), σ τ full ref w.r.t. T τ. Then η T (f, u τ ) T {T } {neighbours} u σ u τ 2 H 1 ( T ) 6/16
8 Corol 5. σ τ full ref w.r.t. T M τ. Then c 2 [ η T (f, u τ )] 2 1 uσ u τ (not true for any f L 2 (Ω)) T M In part, c 2 E(τ, f, u τ ) u u τ. AFEM converges (Dörfler 96, Morin, Nochetto, Siebert 00) MARK[τ, f, u τ ] M: Let θ (0, 1]. Select smallest M τ s.t. [ η T (τ, f, u τ )] 1 2 θ E(τ, f, uτ ). T M REFINE[τ, M] σ: Construct smallest (conforming) σ τ that is a full refinement w.r.t. all T M. u u τ 2 = u u σ 2 + u σ u τ 2 u σ u τ c 2 θe(τ, f, u τ ) c 2θ C 1 u u τ u u σ (1 c2 2θ 2 ) 2 1 u uτ C 2 1 u τ u u σ 7/16 V σ
9 AFEM converges with optimal rate [St 06] ([Binev,Dahmen,DeVore 04] using coarsening) AFEM[f, ε] [τ m, u τm ] % let f T τ 0 P d 2 (T ) compute Gal sol u τ0 V τ0 ; k := 0 while C 1 E(τ k, f, u τk ) > ε do M k := MARK[τ k, f, u τk ] τ k+1 := REFINE[τ k, M k ] compute Gal sol u τk+1 V τk+1 k := k + 1 end do m := k Th 6. If θ < c 2 C 1, then u u τm ε, #τ m #τ 0 m 1 k=0 #M k #M k inf { #ρ #τ 0 : u u ρ [1 C2 1 θ2 ] 1 2 u u c 2 τk } 2 [ [1 C2 1 θ2 ] 1 2 u u c 2 τk ] 1/s 1/s u A s when u As 2. #τ m #τ 0 m 1 k=0 u u τk 1/s u 1/s A s u u τ m 1 1/s u 1/s A s ε 1/s u 1/s A s. 8/16
10 AFEM converges with opt rate Pr. Th. 6. Let σ τ k s.t. u u σ [1 C2 1 θ2 ] 1 2 u u c 2 τk. Th.3 shows 2 C 2 1 T R τk σ η T (f, u τk ) u σ u τk 2 = u u τk 2 u u σ 2 i.e., [ C2 1θ 2 c 2 u u τk 2 C1θ 2 2 E(τ k, f, u τk ) 2, 2 η T (f, u τk )] 1 2 θe(τk, f, u τk ). T R τk σ M := MARK[τ k, f, u τk ] is smallest set with this prop, so #M #R τk σ #σ #τ k. For arb. ρ with u u ρ [1 C2 1 θ2 ] 1 2 u u c 2 τk, take σ := τ k ρ. Then 2 #M #σ #τ k #ρ #τ 0. 9/16
11 General right hand sides and inexact solves RHS[τ, f, δ] [σ, f σ ] % In: τ a partition, f H 1 (Ω) and δ > 0. % Out: f σ T τ P d 2(T ), where σ = τ, or, if necessary, σ τ, % such that f f σ H 1 (Ω) δ. If u A s, cost of RHS will not dominate if c f s.t. #σ #τ c 1/s f δ 1/s, and cost #σ. Such a pair (f, RHS) is called s-optimal. For suff. smooth f, s-optimality with s = d n (> d 1 n ) can be realized. GALSOLVE[τ, f τ, u (0) τ, δ] w τ % In: τ, f τ (S τ ), and u (0) τ S τ. % Out: w τ S τ with u τ w τ H 1 (Ω) δ. % The call should require max{1, log(δ 1 u τ u (0) τ H 1 (Ω))}#τ ops. 10/16
12 AFEM has opt compl AFEM[f, ε] [τ m, w τm ] % ω, β > 0 suff small constants. w τ0 := 0; k := 0; δ 0 f H 1 (Ω) do do δ k := δ k /2 [τ k, f τk ] := RHS[τ k, f, δ k /2] w τk := GALSOLVE[τ k, f τk, w τk, δ k /2] if η k := (2 + C 1 c 1 2 )δ k/2 + C 1 E(τ k, f τk, w τk ) ε then stop endif until δ k ωe(τ k, f τk, w τk ). M k := MARK[τ k, f τk, w τk ] τ k+1 := REFINE[τ k, M k ] w τk+1 := w τk, δ k+1 := 2βη k, k := k + 1 enddo m := k Th 7. u w τm H 1 (Ω) ε. If u A s, and (f, RHS) is s-optimal, then both #τ m and work ε 1/s ( u 1/s A s + c1/s f ). 11/16
13 Realizations of RHS ex. f L 2 (Ω). Q σ being L 2 -orth.proj. onto T σ P d 2(T ). f Q σ f H 1 (Ω) T σ vol(t ) 2 n f Q T f 2 L 2 (T ) =: osc(f, σ) Given θ (0, 1), run MARK and REFINE algorithm on osc until δ. Is quasi-optimal in the sense that whenever f Ās := {f L 2 (Ω) : sup N N s inf osc(f, τ) < }, {τ:#τ #τ 0 N} then (f, RHS) is s-optimal (assuming T fp d 2 exact in O(1) operations) So greedy works (thanks to factor vol(t ) 2 n). 12/16
14 Realizations of RHS ex. General d, n > 1. ( n ) 1 < q 2, f L q (Ω) ( H 1 (Ω)). Generalization of previous case with osc(f, σ) := T σ vol(t ) 1 2 q + 2 n inf p P d 2 (T ) f p 2 L q (T ). ex. n = 2, d = 2, f(v) = K v, K a smooth curve. RHS[τ, f, δ] [σ, f σ]: Refine those T that intersect K until their diameters δ 2. (f σ ) T := length(k T ) vol(t ). s-optimal for s = 1 2 (=d 1 n ). 13/16
15 References [1] P. Binev, W. Dahmen, and R. DeVore. Adaptive finite element methods with convergence rates. Numer. Math., 97(2): , [2] P. Binev, W. Dahmen, R. DeVore, and P. Petruchev. Approximation classes for adaptive methods. Serdica Math. J., 28: , [3] S. Dahlke and R. DeVore. Besov regularity for elliptic boundary value problems. Comm. Partial Differential Equations, 22(1 & 2):1 16, [4] W. Dörfler. A convergent adaptive algorithm for Poisson s equation. SIAM J. Numer. Anal., 33: , [5] P. Morin, R. Nochetto, and K. Siebert. Data oscillation and convergence of adaptive FEM. SIAM J. Numer. Anal., 38(2): , /16
16 [6] R.S. Optimality of a standard adaptive finite element method. Found. Comput. Math., 7(2): , [7] R.S. The completion of locally refined simplicial partitions created by bisection. Math. Comp., 77: , [8] R. Verfürth. A Review of A Posteriori Error Estimation and Adaptive Mesh-Refinement Techniques. Wiley-Teubner, Chichester, /16
17 Extensions [Cascon,Kreuzer,Nochetto,Siebert 07]: One bisection of marked cells suffices (no interior node). Marking for reducing osc can be omitted, assuming f L 2 (Ω), exact integration, and with exact solving of Gal systems. L. Chen, M.J. Holst, and J. Xu. Convergence and optimality of adaptive mixed finite element methods. Accepted by Mathematics of Computation, Y. Kondratyuk, R. S. An optimal adaptive algorithm for the Stokes problem. To appear in SIAM J. Numer. Anal. 16/16
Convergence and optimality of an adaptive FEM for controlling L 2 errors
Convergence and optimality of an adaptive FEM for controlling L 2 errors Alan Demlow (University of Kentucky) joint work with Rob Stevenson (University of Amsterdam) Partially supported by NSF DMS-0713770.
More 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 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 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 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 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 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 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 informationAxioms of adaptivity. ASC Report No. 38/2013. C. Carstensen, M. Feischl, M. Page, D. Praetorius
ASC Report No. 38/2013 Axioms of adaptivity C. Carstensen, M. Feischl, M. Page, D. Praetorius Institute for Analysis and Scientific Computing Vienna University of Technology TU Wien www.asc.tuwien.ac.at
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 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 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 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 informationADAPTIVE MULTILEVEL CORRECTION METHOD FOR FINITE ELEMENT APPROXIMATIONS OF ELLIPTIC OPTIMAL CONTROL PROBLEMS
ADAPTIVE MULTILEVEL CORRECTION METHOD FOR FINITE ELEMENT APPROXIMATIONS OF ELLIPTIC OPTIMAL CONTROL PROBLEMS WEI GONG, HEHU XIE, AND NINGNING YAN Abstract: In this paper we propose an adaptive multilevel
More informationComputation of operators in wavelet coordinates
Computation of operators in wavelet coordinates Tsogtgerel Gantumur and Rob Stevenson Department of Mathematics Utrecht University Tsogtgerel Gantumur - Computation of operators in wavelet coordinates
More informationAdaptive Finite Element Algorithms
Adaptive Finite Element Algorithms of optimal complexity Adaptieve Eindige Elementen Algoritmen van optimale complexiteit (met een samenvatting in het Nederlands) PROEFSCHRIF ter verkrijging van de graad
More informationConstruction of wavelets. Rob Stevenson Korteweg-de Vries Institute for Mathematics University of Amsterdam
Construction of wavelets Rob Stevenson Korteweg-de Vries Institute for Mathematics University of Amsterdam Contents Stability of biorthogonal wavelets. Examples on IR, (0, 1), and (0, 1) n. General domains
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 informationAdaptive Boundary Element Methods Part 1: Newest Vertex Bisection. Dirk Praetorius
CENTRAL School on Analysis and Numerics for PDEs November 09-12, 2015 Part 1: Newest Vertex Bisection Dirk Praetorius TU Wien Institute for Analysis and Scientific Computing Outline 1 Regular Triangulations
More informationEnhancing eigenvalue approximation by gradient recovery on adaptive meshes
IMA Journal of Numerical Analysis Advance Access published October 29, 2008 IMA Journal of Numerical Analysis Page 1 of 15 doi:10.1093/imanum/drn050 Enhancing eigenvalue approximation by gradient recovery
More informationSpace-Time Adaptive Wavelet Methods for Parabolic Evolution Problems
Space-Time Adaptive Wavelet Methods for Parabolic Evolution Problems Rob Stevenson Korteweg-de Vries Institute for Mathematics University of Amsterdam joint work with Christoph Schwab (ETH, Zürich), Nabi
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 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 informationConvergence Rates for AFEM: PDE on Parametric Surfaces
Department of Mathematics and Institute for Physical Science and Technology University of Maryland Joint work with A. Bonito and R. DeVore (Texas a&m) J.M. Cascón (Salamanca, Spain) K. Mekchay (Chulalongkorn,
More informationOPTIMAL MULTILEVEL METHODS FOR GRADED BISECTION GRIDS
OPTIMAL MULTILEVEL METHODS FOR GRADED BISECTION GRIDS LONG CHEN, RICARDO H. NOCHETTO, AND JINCHAO XU Abstract. We design and analyze optimal additive and multiplicative multilevel methods for solving H
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 informationSchur Complements on Hilbert Spaces and Saddle Point Systems
Schur Complements on Hilbert Spaces and Saddle Point Systems Constantin Bacuta Mathematical Sciences, University of Delaware, 5 Ewing Hall 976 Abstract For any continuous bilinear form defined on a pair
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 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 informationKey words. a posteriori error estimators, quasi-orthogonality, adaptive mesh refinement, error and oscillation reduction estimates, optimal meshes.
CONVERGENCE OF ADAPTIVE FINITE ELEMENT METHOD FOR GENERAL ECOND ORDER LINEAR ELLIPTIC PDE KHAMRON MEKCHAY AND RICARDO H. NOCHETTO Abstract. We prove convergence of adaptive finite element methods (AFEM)
More 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 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 informationAdaptive Boundary Element Methods Part 2: ABEM. Dirk Praetorius
CENTRAL School on Analysis and Numerics for PDEs November 09-12, 2015 Adaptive Boundary Element Methods Part 2: ABEM Dirk Praetorius TU Wien Institute for Analysis and Scientific Computing Outline 1 Boundary
More informationAdaptive Finite Element Approximations for Density Functional Models
HPPC Methods and Algorithms for Materials Defects, IMS Singapore, February 12, 2015 Adaptive Finite Element Approximations for Density Functional Models Aihui Zhou LSEC, Institute of Computational Mathematics
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 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 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 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 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 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 informationelement stiffness matrix, 21
Bibliography [1] R. Adams, Sobolev Spaces, Academic Press, 1975. [2] C. Amrouche, C. Bernardi, M. Dauge, and V. Girault, Vector potentials in three-diemnsional nonsmooth domains, Math. Meth. Appl. Sci.,
More informationIsogeometric mortaring
Isogeometric mortaring E. Brivadis, A. Buffa, B. Wohlmuth, L. Wunderlich IMATI E. Magenes - Pavia Technical University of Munich A. Buffa (IMATI-CNR Italy) IGA mortaring 1 / 29 1 Introduction Splines Approximation
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 informationPolynomial Preserving Recovery for Quadratic Elements on Anisotropic Meshes
Polynomial Preserving Recovery for Quadratic Elements on Anisotropic Meshes Can Huang, 1 Zhimin Zhang 1, 1 Department of Mathematics, Wayne State University, Detroit, Michigan 480 College of Mathematics
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 informationA Posteriori Error Estimation for Highly Indefinite Helmholtz Problems
Computational Methods in Applied Mathematics Vol. 13 013), No. 3, pp. 333 347 c 013 Institute of Mathematics, NAS of Belarus Doi: 10.1515/cmam-013-0008 A Posteriori Error Estimation for Highly Indefinite
More 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 u + b u + cu = f in Ω, (1.1)
A WEIGHTED H(div) LEAST-SQUARES METHOD FOR SECOND-ORDER ELLIPTIC PROBLEMS Z. CAI AND C. R. WESTPHAL Abstract. This paper presents analysis of a weighted-norm least squares finite element method for elliptic
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 informationTwo-Scale Composite Finite Element Method for Dirichlet Problems on Complicated Domains
Numerische Mathematik manuscript No. will be inserted by the editor) Two-Scale Composite Finite Element Method for Dirichlet Problems on Complicated Domains M. Rech, S. Sauter, A. Smolianski Institut für
More informationAdaptive Uzawa algorithm for the Stokes equation
ASC Report No. 34/2018 Adaptive Uzawa algorithm for the Stokes equation G. Di Fratta. T. Fu hrer, G. Gantner, and D. Praetorius Institute for Analysis and Scientific Computing Vienna University of Technology
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 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 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 information1. Introduction. In this paper, we shall design and analyze finite element approximations. q i δ i in R d,d=2, 3,
SIAM J. NUMER. ANAL. Vol. 45, No. 6, pp. 2298 2320 c 2007 Society for Industrial and Applied Mathematics THE FINITE ELEMENT APPROXIMATION OF THE NONLINEAR POISSON BOLTZMANN EQUATION LONG CHEN, MICHAEL
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 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 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 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 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 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 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 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 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 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 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 for Weak Galerkin Finite Element Methods for Second Order Elliptic Problems
DOI 10.1007/s10915-013-9771-3 A Posteriori Error Estimates for Weak Galerkin Finite Element Methods for Second Order Elliptic Problems Long Chen Junping Wang Xiu Ye Received: 29 January 2013 / Revised:
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 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 informationSchemes. Philipp Keding AT15, San Antonio, TX. Quarklet Frames in Adaptive Numerical. Schemes. Philipp Keding. Philipps-University Marburg
AT15, San Antonio, TX 5-23-2016 Joint work with Stephan Dahlke and Thorsten Raasch Let H be a Hilbert space with its dual H. We are interested in the solution u H of the operator equation Lu = f, with
More informationRELIABLE A POSTERIORI ERROR CONTROL FOR NONCONFORMING FINITE ELEMENT APPROXIMATION OF STOKES FLOW
MATHEMATICS OF COMPUTATION Volume 74, Number 252, Pages 1599 1619 S 0025-5718(05)01743-6 Article electronically published on January 3, 2005 RELIABLE A POSTERIORI ERROR CONTROL FOR NONCONFORMING FINITE
More informationPARTITION OF UNITY FOR THE STOKES PROBLEM ON NONMATCHING GRIDS
PARTITION OF UNITY FOR THE STOES PROBLEM ON NONMATCHING GRIDS CONSTANTIN BACUTA AND JINCHAO XU Abstract. We consider the Stokes Problem on a plane polygonal domain Ω R 2. We propose a finite element method
More informationA posteriori error estimates, stopping criteria, and adaptivity for multiphase compositional Darcy flows in porous media
A posteriori error estimates, stopping criteria, and adaptivity for multiphase compositional Darcy flows in porous media D. A. Di Pietro, M. Vohraĺık, and S. Yousef Université Montpellier 2 Marseille,
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 informationNumerische Mathematik
Numer. Math. (2004) 97: 193 217 Digital Object Identifier (DOI) 10.1007/s00211-003-0493-6 Numerische Mathematik Fast computation in adaptive tree approximation Peter Binev, Ronald DeVore Department of
More informationASYMPTOTICALLY EXACT A POSTERIORI ESTIMATORS FOR THE POINTWISE GRADIENT ERROR ON EACH ELEMENT IN IRREGULAR MESHES. PART II: THE PIECEWISE LINEAR CASE
MATEMATICS OF COMPUTATION Volume 73, Number 246, Pages 517 523 S 0025-5718(0301570-9 Article electronically published on June 17, 2003 ASYMPTOTICALLY EXACT A POSTERIORI ESTIMATORS FOR TE POINTWISE GRADIENT
More 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 Xuejun Xu, Huangxin Chen, Ronald H.W. Hoppe Optimality of Local Multilevel Methods on Adaptively Refined Meshes for Elliptic Boundary Value
More informationABSTRACT CONVERGENCE OF ADAPTIVE FINITE ELEMENT METHODS. Khamron Mekchay, Doctor of Philosophy, 2005
ABSTRACT Title of dissertation: CONVERGENCE OF ADAPTIVE FINITE ELEMENT METHODS Khamron Mekchay, Doctor of Philosophy, 2005 Dissertation directed by: Professor Ricardo H. Nochetto Department of Mathematics
More 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 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 Finite Element Method Using Singular Functions for Poisson Equations: Mixed Boundary Conditions
A Finite Element Method Using Singular Functions for Poisson Equations: Mixed Boundary Conditions Zhiqiang Cai Seokchan Kim Sangdong Kim Sooryun Kong Abstract In [7], we proposed a new finite element method
More informationarxiv: v2 [math.na] 17 Jun 2010
Numerische Mathematik manuscript No. (will be inserted by the editor) Local Multilevel Preconditioners for Elliptic Equations with Jump Coefficients on Bisection Grids Long Chen 1, Michael Holst 2, Jinchao
More informationAdaptive Approximation of Curves
Adaptive Approximation of Curves Peter Binev, Wolfgang Dahmen, Ronald DeVore, Nira Dyn July 15, 2004 Abstract We propose adaptive multiscale refinement algorithms for approximating and encoding curves
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 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 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 informationA gradient recovery method based on an oblique projection and boundary modification
ANZIAM J. 58 (CTAC2016) pp.c34 C45, 2017 C34 A gradient recovery method based on an oblique projection and boundary modification M. Ilyas 1 B. P. Lamichhane 2 M. H. Meylan 3 (Received 24 January 2017;
More informationc 2008 Society for Industrial and Applied Mathematics
SIAM J. NUMER. ANAL. Vol. 46, No. 3, pp. 640 65 c 2008 Society for Industrial and Applied Mathematics A WEIGHTED H(div) LEAST-SQUARES METHOD FOR SECOND-ORDER ELLIPTIC PROBLEMS Z. CAI AND C. R. WESTPHAL
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 Ronald H.W. Hoppe, Irwin Yousept Adaptive Edge Element Approximation of H(curl-Elliptic Optimal Control Problems with Control Constraints Preprint
More informationConvergence of an adaptive hp finite element strategy in higher space-dimensions
Convergence of an adaptive hp finite element strategy in higher space-dimensions M. Bürg W. Dörfler Preprint Nr. 1/6 INSTITUT FÜR WISSENSCHAFTLICHES RECHNEN UND MATHEMATISCHE MODELLBILDUNG Anschriften
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 informationDISCRETE MAXIMUM PRINCIPLES in THE FINITE ELEMENT SIMULATIONS
DISCRETE MAXIMUM PRINCIPLES in THE FINITE ELEMENT SIMULATIONS Sergey Korotov BCAM Basque Center for Applied Mathematics http://www.bcamath.org 1 The presentation is based on my collaboration with several
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 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 informationA mixed finite element approximation of the Stokes equations with the boundary condition of type (D+N)
wwwijmercom Vol2, Issue1, Jan-Feb 2012 pp-464-472 ISSN: 2249-6645 A mixed finite element approximation of the Stokes equations with the boundary condition of type (D+N) Jaouad El-Mekkaoui 1, Abdeslam Elakkad
More informationOn a residual-based a posteriori error estimator for the total error
On a residual-based a posteriori error estimator for the total error J. Papež Z. Strakoš December 28, 2016 Abstract A posteriori error analysis in numerical PDEs aims at providing sufficiently accurate
More information1 Definition of the Riemann integral
MAT337H1, Introduction to Real Analysis: notes on Riemann integration 1 Definition of the Riemann integral Definition 1.1. Let [a, b] R be a closed interval. A partition P of [a, b] is a finite set of
More informationDaniel Kessler 1, Ricardo H. Nochetto 1, 2 and Alfred Schmidt 3
ESAIM: MAN MAN, Vol. 38, N o 1, 004, pp. 19 14 DOI: 10.1051/man:004006 ESAIM: Mathematical Modelling and Numerical Analysis A POSTERIORI ERROR CONTROL FOR THE ALLEN CAHN PROBLEM: CIRCUMVENTING GRONWALL
More informationA POSTERIORI ERROR ESTIMATES BY RECOVERED GRADIENTS IN PARABOLIC FINITE ELEMENT EQUATIONS
A POSTERIORI ERROR ESTIMATES BY RECOVERED GRADIENTS IN PARABOLIC FINITE ELEMENT EQUATIONS D. LEYKEKHMAN AND L. B. WAHLBIN Abstract. This paper considers a posteriori error estimates by averaged gradients
More informationQuasi-optimal convergence rates for adaptive boundary element methods with data approximation, Part II: Hyper-singular integral equation
ASC Report No. 30/2013 Quasi-optimal convergence rates for adaptive boundary element methods with data approximation, Part II: Hyper-singular integral equation M. Feischl, T. Führer, M. Karkulik, J.M.
More information