Recovery-Based A Posteriori Error Estimation
|
|
- Bernice Hensley
- 6 years ago
- Views:
Transcription
1 Recovery-Based A Posteriori Error Estimation Zhiqiang Cai Purdue University Department of Mathematics, Purdue University Slide 1, March 2, 2011
2 Outline Introduction Diffusion Problems Higher Order Elements Convection-Diffusion-Reaction Problems Acknowledgement This work is supported in part by the National Science Foundation and by Lawrence Livermore National Laboratory. Department of Mathematics, Purdue University Slide 2, March 2, 2011
3 Recovery Estimators of ZZ Type recovery-based estimators (Zienkiewicz & Zhu 87) ξ = G( u h ) u h where G : L 2 (Ω) d U h C 0 (Ω) d is gradient recovery operator saturation assumption: u G( u h ) β u u h = 1 β efficiency and reliability bounds reliability bound efficiency bound there exists a constant β [0, 1) s.t. ξ G( u h ) u 1 + β (Carstensen, Zhou, etc.) u u h C r η η K C e u u h ωk K T Department of Mathematics, Purdue University Slide 3, March 2, 2011
4 Recovery Estimators of ZZ Type recovery-based estimators ξ = G( u h ) u h where G : L 2 (Ω) d U h C 0 (Ω) d is gradient recovery operator recovery operators Zienkiewicz & Zhu estimator (87) G(u h )(z) = 1 u h dx ω z ω z z N L 2 -projection find G(u h ) U h such that (G(u h ), τ ) = ( u, τ ) τ U h other recovery techniques survey article by Z. Zhang 07 Department of Mathematics, Purdue University Slide 4, March 2, 2011
5 Recovery Error Estimators of ZZ Type recovery-based estimators + simple, universal, asymptotically exact inefficiency for nonsmooth problems, unreliable on coarse meshes, higher-order finite elements, complex systems, etc. Department of Mathematics, Purdue University Slide 5, March 2, 2011
6 Diffusion Problems elliptic interface problems (A(x) u) = f in Ω R d u = g on Γ D and n (A(x) u) = h on Γ N where A(x) = a(x) I d d and a(x) is positive piecewise constant w.r.t Ω = n i=1 Ω i : a(x) = a i > 0 in Ω i smoothness u H 1+α (Ω) where α > 0 could be very small. Department of Mathematics, Purdue University Slide 6, March 2, 2011
7 A Test Problem with Intersecting Interfaces the test problem Ω = ( 1, 1) 2, Γ D = Ω, f = 0 and a(x) = in (0, 1) 2 ( 1, 0) 2 1 in Ω \ ([0, 1] 2 [ 1, 0] 2 ) exact solution with µ(θ) being smooth u(r, θ) = r 0.1 µ(θ) H 1.1 ɛ (Ω) Department of Mathematics, Purdue University Slide 7, March 2, 2011
8 3443 nodes mesh generated by η ZZ for 50% relative error Department of Mathematics, Purdue University Slide 8, March 2, 2011
9 Why Does It Fail? true gradient and flux for interface problems u / C 0 (Ω) d and σ = A u C 0 (Ω) d recovery spaces G( u h ) C 0 (Ω) d and G( A u h ) C 0 (Ω) d the reason of the failure approximating discontinuous functions by continuous functions Department of Mathematics, Purdue University Slide 9, March 2, 2011
10 How to Fix It? true gradient and flux u H 1 (Ω) = u H(curl, Ω) σ = A u H(div, Ω) conforming elements u h H 1 (Ω) = u h H(curl, Ω) σ h = A u h / H(div, Ω) = G(σ h) RT 0 or BDM 1 mixed elements u h L 2 (Ω) σ h H(div, Ω) = G( A 1 σ h ) D 1 or N 1 Department of Mathematics, Purdue University Slide 10, March 2, 2011
11 Efficient Recovery Estimator for Linear Elements flux σ = A u H(div; Ω) (Cai-Zhang 09, SINUM) L 2 flux recovery find σ h V = RT 0 or BDM 1 H(div; Ω) s.t. (A 1 σ h, τ ) = ( u h, τ ) τ V efficient recovery estimator ξ 1,K = A 1/2 ( σ h + A ũ h ) 0,K K T ξ 1 = A 1/2 ( σ h + A ũ h ) 0,Ω reliability bound u h u h C (ξ 1 + H f ) efficiency bound ξ 1 C u ũ h + C ( ) h 2 1/2 K f α f K 2 0,K K Department of Mathematics, Purdue University Slide 11, March 2, 2011
12 3557 nodes mesh generated by ξ RT for 10% relative error Department of Mathematics, Purdue University Slide 12, March 2, 2011
13 Explicit L 2 Flux Recovery RT 0 (Cai & Zhang SINUM 09) Let τ = a(x) u h, then σ T = e E Ω E D σ e φ e (x), where σ e = γ e (τ + e n e ) + (1 γ e ) (τ e n e ) for e E Ω, τ e n e for e E D for some constant γ e [0, 1], e.g., a e γ e = a + e + a e, γ e = a e φ e 2 K + e a e φ e 2 K + e + a + e φ e 2 K e Department of Mathematics, Purdue University Slide 13, March 2, 2011
14 Recovery Estimators using the Hypercircle Method (Lade veze & Leguillon 83, Vejchodsky 04, Braess & Schöberl 07) Prager-Synge identity A 1/2 (u u h ) 2 + A 1/2 (u + A 1 τ ) 2 = A 1/2 (u h + A 1 τ ) 2 where τ H(div; Ω) satisfying τ = f recovery estimators by the hypercircle method construct σ h = η K = A 1/2 (u h + A 1 σ h ) K Department of Mathematics, Purdue University Slide 14, March 2, 2011
15 Recovery Estimators for Higher-Order Elements existing recovery estimators: Naga and Zhang (05), Bank, Xu, and Zheng (07) Department of Mathematics, Purdue University Slide 15, March 2, 2011
16 Recovery Estimators for Higher-Order Elements (Cai-Zhang 10, SINUM) sources of errors for conforming elements: the element residual h K f + div (A u h ) 0,K the edge jump residual h 1/2 e [n e (A u h )] 0,e H(div) flux recovery find σ 2 h V = RT k 1 or BDM k s.t. (A 1 σ 2 h, τ ) + ( σ 2 h, τ ) = ( u h, τ ) + (f, τ ) τ V H(div) recovery-based estimator ξ 2,K = A 1/2 (σ 2 h + A u h ) 0,K K T ξ 2 = A 1/2 (σ 2 h + A u h ) 0,Ω Department of Mathematics, Purdue University Slide 16, March 2, 2011
17 Recovery Estimators for Higher-Order Elements Set e = u ũ h, E = σ σ 2 h, and (E, e) 1,Ω = ( A 1/2 E 2 0,K + A1/2 e 2 0,K ) reliability bound (E, e) 1,Ω C ( ξ 2 + h(f div σ 2 ) h Xũ h ) 0,Ω ) = C (ξ 2 + h(a 1/2 E + A 1/2 e) 0,Ω efficiency bound ξ 2,K A 1/2 E 0,K + A 1/2 e 0,K 2 (E, e) 1,K Department of Mathematics, Purdue University Slide 17, March 2, 2011
18 Convection-Diffusion-Reaction Problems convection-diffusion-reaction problems (A(x) u) + b u + a 0 (x)u = f in Ω R d u = g on Γ D and n (A(x) u) = h on Γ N where A d d (x) is uniformly elliptic in Ω possible computational difficulties corner and interface singularities discontinuities in the form of shock-like fronts, and of interior and boundary layers oscillations of various scales Department of Mathematics, Purdue University Slide 18, March 2, 2011
19 Efficient and Reliable Recovery Estimators (Cai-Zhang 10, SINUM) L 2 -based recovery estimator η 1,K = ( ξk 2 + βk 2 σ 1 h + Xũ h f 2 1/2 0, K) K T ( ) 1/2 η 1 = ξ 2 + K T β 2 K σ 1 h + Xũ h f 2 0, K (a similar but different idea used by Fierro & Veeser 2006) H(div)-based recovery estimator η 2,K = ( ξ 2 K + σ 2 h + Xũ h f 2 0, K) 1/2 K T ( ) 1/2 η 2 = ξ 2 + K T σ 2 h + Xũ h f 2 0, K Department of Mathematics, Purdue University Slide 19, March 2, 2011
20 Efficient and Reliable Recovery Estimators true errors e = u ũ h, E 1 = σ σ 1 h, and E 2 = σ σ 2 h reliability bound for ũ h = u h (E 1, e) 1,Ω C r η 1 efficiency bound exactness η 1,K C e (E 1, e) 1,K + β K f f T 0,K ( ) 1/2 η 1 Ĉe e Ω + βk f 2 f T 2 0,K K T η 2,K = (E 2, e) K and η 2 = (E 2, e) Ω where (E 2, e) Ω = ( A 1/2 (E 2 + A e) 2 0,Ω + E 2 + Xe 2 0,Ω ) 1 2 Department of Mathematics, Purdue University Slide 20, March 2, 2011
21 A Test Problem with Highly Oscillatory Solution Poisson equation u = µ sin(2 m πx) in I = (0, 1), and u(0) = u(1) = 0. where m is a fixed integer and µ is an arbitrary constant exact solution u = µ 4 m π 2 sin(2m πx) which is highly oscillatory for large m existing recovery-based estimators mesh x k = k/2 n for k = 0,..., 2 n, when m > n, estimators of recovery type equal to zero starting with uniform the true error is proportional to µ, arbitrarily large Department of Mathematics, Purdue University Slide 21, March 2, 2011
22 Department of Mathematics, Purdue University Slide 22, March 2, 2011
23 A Test Problem with Intersecting Interfaces the test problem Ω = ( 1, 1) 2, Γ D = Ω, f = 0, a 0 = 0, A(x) = a(x)i and a(x) = in (0, 1) 2 ( 1, 0) 2 1 in Ω \ ([0, 1] 2 [ 1, 0] 2 ) exact solution with µ(θ) being smooth u(r, θ) = r 0.25 µ(θ) H 1.25 ɛ (Ω) Department of Mathematics, Purdue University Slide 23, March 2, 2011
24 number of V-cycle iterations and number of unknowns (edges) (different from observation of Cascon-Nochetto-Siebert 07) Department of Mathematics, Purdue University Slide 24, March 2, 2011
25 Figure 1: mesh generated by ξ with H(div) MG solver Figure 2: ξ vs error Department of Mathematics, Purdue University Slide 25, March 2, 2011
26 Figure 3: mesh generated by ξ with one step V-cycle MG solver Figure 4: ξ with one step V- cycle MG solver and error Department of Mathematics, Purdue University Slide 26, March 2, 2011
27 A Test Problem with Interior Layer reaction-dominant-diffusion problem u + κ 2 u = f in Ω = ( 1, 1) 2 with κ = 100 u = 0 on Ω exact solution u = tanh(κ(x 2 + y 2 1/4)) 1 and u r ( ) 1 2 = 2κ which has an interior layer along the circle centered at the origin with radius 1/2 Department of Mathematics, Purdue University Slide 27, March 2, 2011
28 Figure 5: mesh generated by ξ Figure 6: ξ vs error Department of Mathematics, Purdue University Slide 28, March 2, 2011
29 Figure 7: mesh generated by η Figure 8: η vs error Department of Mathematics, Purdue University Slide 29, March 2, 2011
30 Concluding Remarks recovery-based estimator + detecting interface singularities accurate in the energy norm unreliable on coarse meshes reliable estimator + local and global exactness on any given mesh detecting interface singularities, interior/boundary layers, etc. working well for highly oscillatory solutions requiring a fast solver Department of Mathematics, Purdue University Slide 30, March 2, 2011
A Posteriori Error Estimation Techniques for Finite Element Methods. Zhiqiang Cai Purdue University
A Posteriori Error Estimation Techniques for Finite Element Methods Zhiqiang Cai Purdue University Department of Mathematics, Purdue University Slide 1, March 16, 2017 Books Ainsworth & Oden, A posteriori
More 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 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 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 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 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 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 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 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 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 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 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 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 informationarxiv: v3 [math.na] 8 Sep 2015
A Recovery-Based A Posteriori Error Estimator for H(curl) Interface Problems arxiv:504.00898v3 [math.na] 8 Sep 205 Zhiqiang Cai Shuhao Cao Abstract This paper introduces a new recovery-based a posteriori
More informationNitsche-type Mortaring for Maxwell s Equations
Nitsche-type Mortaring for Maxwell s Equations Joachim Schöberl Karl Hollaus, Daniel Feldengut Computational Mathematics in Engineering at the Institute for Analysis and Scientific Computing Vienna University
More informationAdaptive C1 Macroelements for Fourth Order and Divergence-Free Problems
Adaptive C1 Macroelements for Fourth Order and Divergence-Free Problems Roy Stogner Computational Fluid Dynamics Lab Institute for Computational Engineering and Sciences University of Texas at Austin March
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 informationA POSTERIORI ERROR ESTIMATES OF TRIANGULAR MIXED FINITE ELEMENT METHODS FOR QUADRATIC CONVECTION DIFFUSION OPTIMAL CONTROL PROBLEMS
A POSTERIORI ERROR ESTIMATES OF TRIANGULAR MIXED FINITE ELEMENT METHODS FOR QUADRATIC CONVECTION DIFFUSION OPTIMAL CONTROL PROBLEMS Z. LU Communicated by Gabriela Marinoschi In this paper, we discuss a
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 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 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 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 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 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 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 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 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 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 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 informationTHE CONVECTION DIFFUSION EQUATION
3 THE CONVECTION DIFFUSION EQUATION We next consider the convection diffusion equation ɛ 2 u + w u = f, (3.) where ɛ>. This equation arises in numerous models of flows and other physical phenomena. The
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 informationDiscontinuous Galerkin Methods
Discontinuous Galerkin Methods Joachim Schöberl May 20, 206 Discontinuous Galerkin (DG) methods approximate the solution with piecewise functions (polynomials), which are discontinuous across element interfaces.
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 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 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 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 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 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 informationA Stable Spectral Difference Method for Triangles
A Stable Spectral Difference Method for Triangles Aravind Balan 1, Georg May 1, and Joachim Schöberl 2 1 AICES Graduate School, RWTH Aachen, Germany 2 Institute for Analysis and Scientific Computing, Vienna
More informationAnisotropic meshes for PDEs: a posteriori error analysis and mesh adaptivity
ICAM 2013, Heraklion, 17 September 2013 Anisotropic meshes for PDEs: a posteriori error analysis and mesh adaptivity Simona Perotto MOX-Modeling and Scientific Computing Department of Mathematics, Politecnico
More informationMultigrid Method ZHONG-CI SHI. Institute of Computational Mathematics Chinese Academy of Sciences, Beijing, China. Joint work: Xuejun Xu
Multigrid Method ZHONG-CI SHI Institute of Computational Mathematics Chinese Academy of Sciences, Beijing, China Joint work: Xuejun Xu Outline Introduction Standard Cascadic Multigrid method(cmg) Economical
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 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 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 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 informationGeometric Multigrid Methods
Geometric Multigrid Methods Susanne C. Brenner Department of Mathematics and Center for Computation & Technology Louisiana State University IMA Tutorial: Fast Solution Techniques November 28, 2010 Ideas
More 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 informationHybrid Discontinuous Galerkin methods for incompressible flow problems
Hybrid Discontinuous Galerkin methods incompressible flow problems Christoph Lehrenfeld, Joachim Schöberl MathCCES Computational Mathematics in Engineering Workshop Linz, May 31 - June 1, 2010 Contents
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 informationMesh Grading towards Singular Points Seminar : Elliptic Problems on Non-smooth Domain
Mesh Grading towards Singular Points Seminar : Elliptic Problems on Non-smooth Domain Stephen Edward Moore Johann Radon Institute for Computational and Applied Mathematics Austrian Academy of Sciences,
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 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 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 informationarxiv: v1 [math.na] 1 May 2013
arxiv:3050089v [mathna] May 03 Approximation Properties of a Gradient Recovery Operator Using a Biorthogonal System Bishnu P Lamichhane and Adam McNeilly May, 03 Abstract A gradient recovery operator based
More 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 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 DELTA-REGULARIZATION FINITE ELEMENT METHOD FOR A DOUBLE CURL PROBLEM WITH DIVERGENCE-FREE CONSTRAINT
A DELTA-REGULARIZATION FINITE ELEMENT METHOD FOR A DOUBLE CURL PROBLEM WITH DIVERGENCE-FREE CONSTRAINT HUOYUAN DUAN, SHA LI, ROGER C. E. TAN, AND WEIYING ZHENG Abstract. To deal with the divergence-free
More informationEfficient numerical solution of the Biot poroelasticity system in multilayered domains
Efficient numerical solution of the Biot poroelasticity system Anna Naumovich Oleg Iliev Francisco Gaspar Fraunhofer Institute for Industrial Mathematics Kaiserslautern, Germany, Spain Workshop on Model
More informationAn A Posteriori Error Estimate for Discontinuous Galerkin Methods
An A Posteriori Error Estimate for Discontinuous Galerkin Methods Mats G Larson mgl@math.chalmers.se Chalmers Finite Element Center Mats G Larson Chalmers Finite Element Center p.1 Outline We present an
More 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 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 informationA MULTIGRID ALGORITHM FOR. Richard E. Ewing and Jian Shen. Institute for Scientic Computation. Texas A&M University. College Station, Texas SUMMARY
A MULTIGRID ALGORITHM FOR THE CELL-CENTERED FINITE DIFFERENCE SCHEME Richard E. Ewing and Jian Shen Institute for Scientic Computation Texas A&M University College Station, Texas SUMMARY In this article,
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 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 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 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 informationSolving PDEs with Multigrid Methods p.1
Solving PDEs with Multigrid Methods Scott MacLachlan maclachl@colorado.edu Department of Applied Mathematics, University of Colorado at Boulder Solving PDEs with Multigrid Methods p.1 Support and Collaboration
More informationComparison of V-cycle Multigrid Method for Cell-centered Finite Difference on Triangular Meshes
Comparison of V-cycle Multigrid Method for Cell-centered Finite Difference on Triangular Meshes Do Y. Kwak, 1 JunS.Lee 1 Department of Mathematics, KAIST, Taejon 305-701, Korea Department of Mathematics,
More informationA finite difference Poisson solver for irregular geometries
ANZIAM J. 45 (E) ppc713 C728, 2004 C713 A finite difference Poisson solver for irregular geometries Z. Jomaa C. Macaskill (Received 8 August 2003, revised 21 January 2004) Abstract The motivation for this
More informationStandard Finite Elements and Weighted Regularization
Standard Finite Elements and Weighted Regularization A Rehabilitation Martin COSTABEL & Monique DAUGE Institut de Recherche MAthématique de Rennes http://www.maths.univ-rennes1.fr/~dauge Slides of the
More informationOn Nonlinear Dirichlet Neumann Algorithms for Jumping Nonlinearities
On Nonlinear Dirichlet Neumann Algorithms for Jumping Nonlinearities Heiko Berninger, Ralf Kornhuber, and Oliver Sander FU Berlin, FB Mathematik und Informatik (http://www.math.fu-berlin.de/rd/we-02/numerik/)
More informationConvergence and Error Bound Analysis for the Space-Time CESE Method
Convergence and Error Bound Analysis for the Space-Time CESE Method Daoqi Yang, 1 Shengtao Yu, Jennifer Zhao 3 1 Department of Mathematics Wayne State University Detroit, MI 480 Department of Mechanics
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 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 informationChapter 1. Introduction and Background. 1.1 Introduction
Chapter 1 Introduction and Background 1.1 Introduction Over the past several years the numerical approximation of partial differential equations (PDEs) has made important progress because of the rapid
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Numerical Methods for Partial Differential Equations Finite Difference Methods
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 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 informationGeneralized Finite Element Methods for Three Dimensional Structural Mechanics Problems. C. A. Duarte. I. Babuška and J. T. Oden
Generalized Finite Element Methods for Three Dimensional Structural Mechanics Problems C. A. Duarte COMCO, Inc., 7800 Shoal Creek Blvd. Suite 290E Austin, Texas, 78757, USA I. Babuška and J. T. Oden TICAM,
More informationFEniCS Course. Lecture 8: A posteriori error estimates and adaptivity. Contributors André Massing Marie Rognes
FEniCS Course Lecture 8: A posteriori error estimates and adaptivity Contributors André Massing Marie Rognes 1 / 24 A priori estimates If u H k+1 (Ω) and V h = P k (T h ) then u u h Ch k u Ω,k+1 u u h
More 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 informationLocal Mesh Refinement with the PCD Method
Advances in Dynamical Systems and Applications ISSN 0973-5321, Volume 8, Number 1, pp. 125 136 (2013) http://campus.mst.edu/adsa Local Mesh Refinement with the PCD Method Ahmed Tahiri Université Med Premier
More informationExponential integrators for semilinear parabolic problems
Exponential integrators for semilinear parabolic problems Marlis Hochbruck Heinrich-Heine University Düsseldorf Germany Innsbruck, October 2004 p. Outline Exponential integrators general class of methods
More informationIterative methods for positive definite linear systems with a complex shift
Iterative methods for positive definite linear systems with a complex shift William McLean, University of New South Wales Vidar Thomée, Chalmers University November 4, 2011 Outline 1. Numerical solution
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 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 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 informationA FAMILY OF MULTISCALE HYBRID-MIXED FINITE ELEMENT METHODS FOR THE DARCY EQUATION WITH ROUGH COEFFICIENTS. 1. Introduction
A FAMILY OF MULTISCALE HYBRID-MIXED FINITE ELEMENT METHODS FOR THE DARCY EQUATION WITH ROUGH COEFFICIENTS CHRISTOPHER HARDER, DIEGO PAREDES 2, AND FRÉDÉRIC VALENTIN Abstract. We aim at proposing novel
More informationMultiscale Analysis of Vibrations of Streamers
Multiscale Analysis of Vibrations of Streamers Leszek Demkowicz Joint work with S. Prudhomme, W. Rachowicz, W. Qiu and L. Chamoin Institute for Computational Engineering and Sciences (ICES) The University
More informationHp-Adaptivity on Anisotropic Meshes for Hybridized Discontinuous Galerkin Scheme
Hp-Adaptivity on Anisotropic Meshes for Hybridized Discontinuous Galerkin Scheme Aravind Balan, Michael Woopen and Georg May AICES Graduate School, RWTH Aachen University, Germany 22nd AIAA Computational
More informationUltraconvergence of ZZ Patch Recovery at Mesh Symmetry Points
Ultraconvergence of ZZ Patch Recovery at Mesh Symmetry Points Zhimin Zhang and Runchang Lin Department of Mathematics, Wayne State University Abstract. The ultraconvergence property of the Zienkiewicz-Zhu
More 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 informationFEM Convergence for PDEs with Point Sources in 2-D and 3-D
FEM Convergence for PDEs with Point Sources in -D and 3-D Kourosh M. Kalayeh 1, Jonathan S. Graf, and Matthias K. Gobbert 1 Department of Mechanical Engineering, University of Maryland, Baltimore County
More information1. Introduction. We consider the model problem: seeking an unknown function u satisfying
A DISCONTINUOUS LEAST-SQUARES FINITE ELEMENT METHOD FOR SECOND ORDER ELLIPTIC EQUATIONS XIU YE AND SHANGYOU ZHANG Abstract In tis paper, a discontinuous least-squares (DLS) finite element metod is introduced
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 informationRemarks on the analysis of finite element methods on a Shishkin mesh: are Scott-Zhang interpolants applicable?
Remarks on the analysis of finite element methods on a Shishkin mesh: are Scott-Zhang interpolants applicable? Thomas Apel, Hans-G. Roos 22.7.2008 Abstract In the first part of the paper we discuss minimal
More informationDiscontinuous Galerkin Method for interface problem of coupling different order elliptic equations
Discontinuous Galerkin Method for interface problem of coupling different order elliptic equations Igor Mozolevski, Endre Süli Federal University of Santa Catarina, Brazil Oxford University Computing Laboratory,
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University The Implicit Schemes for the Model Problem The Crank-Nicolson scheme and θ-scheme
More informationNumerical Solution I
Numerical Solution I Stationary Flow R. Kornhuber (FU Berlin) Summerschool Modelling of mass and energy transport in porous media with practical applications October 8-12, 2018 Schedule Classical Solutions
More informationFEM: Domain Decomposition and Homogenization for Maxwell s Equations of Large Scale Problems
FEM: and Homogenization for Maxwell s Equations of Large Scale Problems Karl Hollaus Vienna University of Technology, Austria Institute for Analysis and Scientific Computing February 13, 2012 Outline 1
More information