Numerical lower bounds on eigenvalues of elliptic operators
|
|
- Merryl Walsh
- 5 years ago
- Views:
Transcription
1 Numerical lower bounds on eigenvalues of elliptic operators Tomáš Vejchodský Institute of Mathematics Czech Academy of Sciences INSTITUTEofMATHEMATICS Czech Academy of Sciences EMS School in Applied Mathematics, Kácov, May 29 June 3, 2016
2 What do we do in numerics?
3 What do we do in numerics? We solve problems approximately.
4 What do we do in numerics? We solve problems approximately. What should we do?
5 What do we do in numerics? We solve problems approximately. What should we do? Solve problems within a given accuracy.
6 Example dumbbell u i = λ i u i in Ω u i = 0 on Ω π π π 4 [Trefethen, Betcke 2006]
7 Example dumbbell u i = λ i u i u i = 0 in Ω on Ω λ λ
8 Example dumbbell u i = λ i u i u i = 0 in Ω on Ω λ λ λ λ
9 Example dumbbell u i = λ i u i u i = 0 in Ω on Ω λ λ λ λ λ λ
10 Example dumbbell u i = λ i u i u i = 0 in Ω on Ω N tri = λ λ λ λ λ λ λ λ
11 Example dumbbell u i = λ i u i u i = 0 in Ω on Ω N tri = λ λ λ λ λ λ λ λ λ λ
12 Example dumbbell u i = λ i u i u i = 0 in Ω on Ω λ λ
13 Example dumbbell u i = λ i u i u i = 0 in Ω on Ω λ λ λ λ
14 Example dumbbell u i = λ i u i u i = 0 in Ω on Ω λ λ λ λ λ λ
15 Example dumbbell u i = λ i u i u i = 0 in Ω on Ω N tri = λ λ λ λ λ λ λ λ
16 Example dumbbell u i = λ i u i u i = 0 in Ω on Ω N tri = λ λ λ λ λ λ λ λ λ λ
17 Upper bounds on eigenvalues Laplace eigenvalue problem u i = λ i u i u i = 0 in Ω on Ω Weak formulation λ i R, u i H0 1 (Ω) : ( u i, v) = λ i (u i, v) v H0 1 (Ω) Finite element method V h = {v h H 1 0 (Ω) : v h K P 1 (K), K T h } Λ h,i R, u h,i V h : ( u h,i, v h ) = Λ h,i (u h,i, v h ) v h V h Lower bound?? λ i Λ h,i, i = 1, 2,..., m
18 Lower bounds on eigenvalues Old problem: Temple 1928, Kato 1949, Lehmann 1949, 1950, Harrell 1978,... Methods based on FEM: 1. Eigenvalue inclusions [Behnke, Mertins, Plum, Wieners 2000] based on [Behnke, Goerish 1994] and [Plum 1997] 2. Crouzeix Raviart elements [Carstensen, Gedicke 2013] 3. Complementarity based [Šebestová, Vejchodský 2016]
19 Method 1. Eigenvalue inclusions Input: Rough lower bounds: λ 2 λ 2,..., λ m+1 λ m+1, Algorithm: FEM eigenpairs: Λ h,i R, u h,i V h, i = 1, 2,..., m Mixed FEM problem: σ h,i W h, q h,i Q h, i = 1, 2,..., m W h = {σ h H(div, Ω) : σ h K RT k (K) K T h } Q h = {q h L 2 (Ω) : q h K P k (K) K T h } (σ h,i, w h ) + (q h,i, div w h ) = 0 w h W h, (div σ h,i, ϕ h ) = ( u h,i, ϕ h ) ϕ h Q h,
20 Method 1. Eigenvalue inclusions Input: Rough lower bounds: λ 2 λ 2,..., λ m+1 λ m+1, Algorithm: FEM eigenpairs: Λ h,i R, u h,i V h, i = 1, 2,..., m Mixed FEM problem: σ h,i W h, q h,i Q h, i = 1, 2,..., m For n = 1, 2,..., m do γ = u h,n + div σ h,n L 2 (Ω), ρ = λ n+1 + γ M ij = ( u h,i, u h,j ) + (γ ρ)(u h,i, u h,j ) N ij = ( u h,i, u h,j ) + (γ 2ρ)(u h,i, u h,j ) + ρ 2 (σ h,i, σ h,j ) +(ρ 2 /γ)(u h,i + div σ h,i, u h,j + div σ h,j ) µ 1 µ n : My i = µ i Ny i, i = 1, 2,..., n If N is s.p.d. and if µ i < 0 then l incl j,n = ρ γ ρ/ (1 µ n+1 j) λ j, j = 1, 2,..., n. end for l incl j = max{l incl j,n, n = j, j + 1,..., m} λ j, j = 1, 2,..., m
21 Method 2. Crouzeix Raviart elements Crouzeix Raviart finite elements V CR h = {v h P 1 (T h ) : v h continuous in midpoints of all γ E h } Find 0 uh,i CR Vh CR, λ CR h,i R : ( u CR h,i, v h ) = λ CR h,i (ucr h,i, v h ) Lower bound (no round-off errors) v h V CR h. l CR i = λ CR h,i 1 + κ 2 λ CR h,i h2 max λ i i = 1, 2,... where κ 2 = 1/8 + j 2 1, h max = max K Th diam K
22 Method 2. Crouzeix Raviart elements Crouzeix Raviart finite elements V CR h = {v h P 1 (T h ) : v h continuous in midpoints of all γ E h } Find 0 uh,i CR Vh CR, λ CR h,i R : ( u CR h,i, v h ) = λ CR h,i (ucr h,i, v h ) v h V CR h. Lower bound (inexact solver: Aũ CR i where l CR i = λ CR h,i 1 + κ 2 ( λ CR h,i r B 1 r B 1 κ 2 = 1/8 + j 2 1, h max = max K Th diam K r = Aũ CR i λ CR h,i BũCR i λ CR h,i BũCR i ) ) h 2 max Provided r B 1 < λ CR h,i λ CR h,i λ i i = 1, 2,... is closer to λ CR h,i than to any other discrete eigenvalue λ CR h,j, j i
23 Method 2. Crouzeix Raviart elements Upper bound T h is the red refinement of T h for i = 1, 2,..., m u h,i = I CM ũ CR h,i S, Q R m m with entries S j,k = ( u h,j, u h,k ) and Q j,k = (u h,j, u h,k ) Sy i = Λ i Qy i, i = 1, 2,..., m Λ 1 Λ 2 Λ m λ i Λ i for i = 1, 2,..., m
24 Method 3. Complementarity based FEM eigenpairs: Λ h,i R, u h,i V h, i = 1, 2,..., m Flux reconstruction: q h,i = z N h q z,i Local mixed FEM: q z,i W z, d z,i P 1 (T z) (q z,i, w h ) ωz (d z,i, div w h ) ωz = (ψ z u h,i, w h ) ωz w h W z (div q z,i, ϕ h ) ωz = (r z,i, ϕ h ) ωz ϕ h P1 (T z ) where ωz is the patch of elements around vertex z N h Tz is the set of elements in ω z Wz = {w h H(div, ω z ) : w h K RT 1 (K) K T z and w h n ωz = 0 on Γ ext } { ω z P {vh P 1 (T z ) : 1 (T z ) = ω z v h dx = 0} for z N h \ Ω P 1 (T z ) for z N h Ω rz,i = Λ h,i ψ z u h,i ψ z u h,i
25 Method 3. Complementarity based FEM eigenpairs: Λ h,i R, u h,i V h, i = 1, 2,..., m Flux reconstruction: q h,i = z N h q z,i Local mixed FEM: q z,i W z, d z,i P 1 (T z) (q z,i, w h ) ωz (d z,i, div w h ) ωz = (ψ z u h,i, w h ) ωz w h W z (div q z,i, ϕ h ) ωz = (r z,i, ϕ h ) ωz ϕ h P1 (T z ) Error estimator: η i = u h,i q h,i L 2 (Ω) ( Lower bound: l cmpl 1 = η 1 + Provided Λ h,i 2 ( λ 1 i l cmpl i = Λ h,i ( 1 + λ 1/2 + λ 1 ) 1 i+1 η Λ h,1) 2 /4 1 η i ) 1, i = 2, 3,...
26 Comparison 1. Inclusions 2. CR elements 3. Complementarity convergence ** *** * generality *** * ** a priori info * *** ** DOFs needed * ** *** algebraic err. ** *** *** adaptivity * ** ***
27 Thank you for your attention Tomáš Vejchodský Institute of Mathematics Czech Academy of Sciences INSTITUTEofMATHEMATICS Czech Academy of Sciences EMS School in Applied Mathematics, Kácov, May 29 June 3, 2016
28 Square λ 1 = 2, u 1 (x, y) = sin(x) sin(y) λ 2 = 5, u 2 (x, y) = sin(2x) sin(y) λ 3 = 5, u 3 (x, y) = sin(x) sin(2y)
29 Two squares λ 1 = 2 λ 2 = 5 λ 3 = 5
30 Two squares λ 1 = 2 λ 2 = 2 λ 3 = 5 λ 4 = 5 λ 5 = 5 λ 6 = 5
31 Dumbbell λ λ λ λ λ λ
INSTITUTE of MATHEMATICS. ACADEMY of SCIENCES of the CZECH REPUBLIC
INSTITUTE of MATHEMATICS Academy of Sciences Czech Republic INSTITUTE of MATHEMATICS ACADEMY of SCIENCES of the CZECH REPUBLIC Two-sided bounds for eigenvalues of differential operators with applications
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 informationGuaranteed and robust a posteriori bounds for Laplace eigenvalues and eigenvectors: conforming approximations
Guaranteed and robust a posteriori bounds for Laplace eigenvalues and eigenvectors: conforming approximations Eric Cancès, Geneviève Dusson, Yvon Maday, Benjamin Stamm, Martin Vohralík To cite this version:
More informationDiscrete Maximum Principle for a 1D Problem with Piecewise-Constant Coefficients Solved by hp-fem
The University of Texas at El Paso Department of Mathematical Sciences Research Reports Series El Paso, Texas Research Report No. 2006-10 Discrete Maximum Principle for a 1D Problem with Piecewise-Constant
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 informationWeak Galerkin Finite Element Methods and Applications
Weak Galerkin Finite Element Methods and Applications Lin Mu mul1@ornl.gov Computational and Applied Mathematics Computationa Science and Mathematics Division Oak Ridge National Laboratory Georgia Institute
More 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 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 informationOptimal multilevel preconditioning of strongly anisotropic problems.part II: non-conforming FEM. p. 1/36
Optimal multilevel preconditioning of strongly anisotropic problems. Part II: non-conforming FEM. Svetozar Margenov margenov@parallel.bas.bg Institute for Parallel Processing, Bulgarian Academy of Sciences,
More 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 informationRobust Domain Decomposition Preconditioners for Abstract Symmetric Positive Definite Bilinear Forms
www.oeaw.ac.at Robust Domain Decomposition Preconditioners for Abstract Symmetric Positive Definite Bilinear Forms Y. Efendiev, J. Galvis, R. Lazarov, J. Willems RICAM-Report 2011-05 www.ricam.oeaw.ac.at
More informationAn a posteriori error estimate and a Comparison Theorem for the nonconforming P 1 element
Calcolo manuscript No. (will be inserted by the editor) An a posteriori error estimate and a Comparison Theorem for the nonconforming P 1 element Dietrich Braess Faculty of Mathematics, Ruhr-University
More informationAn 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 informationTHREE MATLAB IMPLEMENTATIONS OF THE LOWEST-ORDER RAVIART-THOMAS MFEM WITH A POSTERIORI ERROR CONTROL
COMPUAIONAL MEHODS IN APPLIED MAHEMAICS, Vol.5(25), No.4, pp.333 36 c 25 Institute of Mathematics of the National Academy of Sciences of Belarus HREE MALAB IMPLEMENAIONS OF HE LOWES-ORDER RAVIAR-HOMAS
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 informationIntroduction to finite element exterior calculus
Introduction to finite element exterior calculus Ragnar Winther CMA, University of Oslo Norway Why finite element exterior calculus? Recall the de Rham complex on the form: R H 1 (Ω) grad H(curl, Ω) curl
More informationarxiv: v2 [math.na] 2 Nov 2018
DISCREIZION-CCURE SOPPING CRIERION FOR IERIVE SOLVERS ON FINIE ELEMEN PPROXIMION ZHIQING CI, SHUHO CO, ND ROB FLGOU arxiv:161108967v2 [mathn] 2 Nov 2018 bstract his paper introduces a localized duality-type
More 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 informationMODEL OF GROUNDWATER FLOW IN FRACTURED ENVIRONMENT
Proceedings of ALGORITMY 2002 Conference on Scientific Computing, pp. 138 145 MODEL OF GROUNDWATER FLOW IN FRACTURED ENVIRONMENT JIŘÍ MARYŠKA, OTTO SEVERÝN AND MARTIN VOHRALÍK Abstract. A stochastic discrete
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 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 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 informationAN ADAPTIVE HOMOTOPY APPROACH FOR NON-SELFADJOINT EIGENVALUE PROBLEMS
AN ADAPTIVE HOMOTOPY APPROACH FOR NON-SELFADJOINT EIGENVALUE PROBLEMS C. CARSTENSEN, J. GEDICKE, V. MEHRMANN, AND A. MIEDLAR Abstract. This paper presents three different adaptive algorithms for eigenvalue
More informationIsogEometric Tearing and Interconnecting
IsogEometric Tearing and Interconnecting Christoph Hofer and Ludwig Mitter Johannes Kepler University, Linz 26.01.2017 Doctoral Program Computational Mathematics Numerical Analysis and Symbolic Computation
More informationFULLY COMPUTABLE ERROR BOUNDS FOR EIGENVALUE PROBLEM
INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume 15, Number 1-2, Pages 260 276 c 2018 Institute for Scientific Computing and Information FULLY COMPUTABLE ERROR BOUNDS FOR EIGENVALUE PROBLEM
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 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 informationA Domain Decomposition Method for Quasilinear Elliptic PDEs Using Mortar Finite Elements
W I S S E N T E C H N I K L E I D E N S C H A F T A Domain Decomposition Method for Quasilinear Elliptic PDEs Using Mortar Finite Elements Matthias Gsell and Olaf Steinbach Institute of Computational Mathematics
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 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 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 informationarxiv: v2 [math.na] 16 Nov 2016
BOOTSTRAP MULTIGRID FOR THE SHIFTED LAPLACE-BELTRAMI EIGENVALUE PROBLEM JAMES BRANNICK AND SHUHAO CAO arxiv:1511.07042v2 [math.na] 16 Nov 2016 Abstract. This paper introduces bootstrap two-grid and multigrid
More informationFinal Ph.D. Progress Report. Integration of hp-adaptivity with a Two Grid Solver: Applications to Electromagnetics. David Pardo
Final Ph.D. Progress Report Integration of hp-adaptivity with a Two Grid Solver: Applications to Electromagnetics. David Pardo Dissertation Committee: I. Babuska, L. Demkowicz, C. Torres-Verdin, R. Van
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 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 information1. Introduction. The Stokes problem seeks unknown functions u and p satisfying
A DISCRETE DIVERGENCE FREE WEAK GALERKIN FINITE ELEMENT METHOD FOR THE STOKES EQUATIONS LIN MU, JUNPING WANG, AND XIU YE Abstract. A discrete divergence free weak Galerkin finite element method is developed
More informationSolving linear equations with Gaussian Elimination (I)
Term Projects Solving linear equations with Gaussian Elimination The QR Algorithm for Symmetric Eigenvalue Problem The QR Algorithm for The SVD Quasi-Newton Methods Solving linear equations with Gaussian
More informationEXPLICIT ERROR ESTIMATES FOR COURANT, CROUZEIX-RAVIART AND RAVIART-THOMAS FINITE ELEMENT METHODS
EXPLICI ERROR ESIMAES FOR COURAN, CROUZEIX-RAVIAR AND RAVIAR-HOMAS FINIE ELEMEN MEHODS CARSEN CARSENSEN, JOSCHA GEDICKE, AND DONSUB RIM Abstract. he elementary analysis of this paper presents explicit
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 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 informationAn Iterative Substructuring Method for Mortar Nonconforming Discretization of a Fourth-Order Elliptic Problem in two dimensions
An Iterative Substructuring Method for Mortar Nonconforming Discretization of a Fourth-Order Elliptic Problem in two dimensions Leszek Marcinkowski Department of Mathematics, Warsaw University, Banacha
More informationLecture Notes: African Institute of Mathematics Senegal, January Topic Title: A short introduction to numerical methods for elliptic PDEs
Lecture Notes: African Institute of Mathematics Senegal, January 26 opic itle: A short introduction to numerical methods for elliptic PDEs Authors and Lecturers: Gerard Awanou (University of Illinois-Chicago)
More informationINSTITUTE of MATHEMATICS. ACADEMY of SCIENCES of the CZECH REPUBLIC
INSTITUTEofMATHEMATICS Academy of Sciences Czech Republic INSTITUTE of MATHEMATICS ACADEMY of SCIENCES of the CZECH REPUBLIC Error estimates for a numerical method for the compressible Navier-Stokes system
More informationSolving the curl-div system using divergence-free or curl-free finite elements
Solving the curl-div system using divergence-free or curl-free finite elements Alberto Valli Dipartimento di Matematica, Università di Trento, Italy or: Why I say to my students that divergence-free finite
More informationAMS subject classifications. Primary, 65N15, 65N30, 76D07; Secondary, 35B45, 35J50
A SIMPLE FINITE ELEMENT METHOD FOR THE STOKES EQUATIONS LIN MU AND XIU YE Abstract. The goal of this paper is to introduce a simple finite element method to solve the Stokes equations. This method is in
More informationA Posteriori Existence in Numerical Computations
Oxford University Computing Laboratory OXMOS: New Frontiers in the Mathematics of Solids OXPDE: Oxford Centre for Nonlinear PDE January, 2007 Introduction Non-linear problems are handled via the implicit
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 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 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 informationApplied/Numerical Analysis Qualifying Exam
Applied/Numerical Analysis Qualifying Exam August 9, 212 Cover Sheet Applied Analysis Part Policy on misprints: The qualifying exam committee tries to proofread exams as carefully as possible. Nevertheless,
More informationMultilevel Preconditioning of Graph-Laplacians: Polynomial Approximation of the Pivot Blocks Inverses
Multilevel Preconditioning of Graph-Laplacians: Polynomial Approximation of the Pivot Blocks Inverses P. Boyanova 1, I. Georgiev 34, S. Margenov, L. Zikatanov 5 1 Uppsala University, Box 337, 751 05 Uppsala,
More informationDivergence-free or curl-free finite elements for solving the curl-div system
Divergence-free or curl-free finite elements for solving the curl-div system Alberto Valli Dipartimento di Matematica, Università di Trento, Italy Joint papers with: Ana Alonso Rodríguez Dipartimento di
More informationTHREE MATLAB IMPLEMENTATIONS OF THE LOWEST-ORDER RAVIART-THOMAS MFEM WITH A POSTERIORI ERROR CONTROL
HREE MALAB IMPLEMENAIONS OF HE LOWES-ORDER RAVIAR-HOMAS MFEM WIH A POSERIORI ERROR CONROL C. BAHRIAWAI AND C. CARSENSEN ABSRAC. he numerical approximation of the Laplace equation with inhomogeneous mixed
More informationEXPLICIT ERROR ESTIMATES FOR COURANT, CROUZEIX-RAVIART AND RAVIART-THOMAS FINITE ELEMENT METHODS
EXPLICI ERROR ESIMAES FOR COURAN, CROUZEIX-RAVIAR AND RAVIAR-HOMAS FINIE ELEMEN MEHODS CARSEN CARSENSEN, JOSCHA GEDICKE, AND DONSUB RIM Abstract. he elementary analysis of this paper presents explicit
More informationStructure of the Hessian
Structure of the Hessian Graham C. Goodwin September 2004 1. Introduction We saw in earlier lectures that a core ingredient in quadratic constrained optimisation problems is the Hessian matrix H. So far
More informationA posteriori error estimates for non conforming approximation of eigenvalue problems
A posteriori error estimates for non conforming approximation of eigenvalue problems E. Dari a, R. G. Durán b and C. Padra c, a Centro Atómico Bariloche, Comisión Nacional de Energía Atómica and CONICE,
More 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 informationSome Domain Decomposition Methods for Discontinuous Coefficients
Some Domain Decomposition Methods for Discontinuous Coefficients Marcus Sarkis WPI RICAM-Linz, October 31, 2011 Marcus Sarkis (WPI) Domain Decomposition Methods RICAM-2011 1 / 38 Outline Discretizations
More informationFinite Element Clifford Algebra: A New Toolkit for Evolution Problems
Finite Element Clifford Algebra: A New Toolkit for Evolution Problems Andrew Gillette joint work with Michael Holst Department of Mathematics University of California, San Diego http://ccom.ucsd.edu/ agillette/
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 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 informationA multipoint flux mixed finite element method on hexahedra
A multipoint flux mixed finite element method on hexahedra Ross Ingram Mary F. Wheeler Ivan Yotov Abstract We develop a mixed finite element method for elliptic problems on hexahedral grids that reduces
More informationOblique derivative problems for elliptic and parabolic equations, Lecture II
of the for elliptic and parabolic equations, Lecture II Iowa State University July 22, 2011 of the 1 2 of the of the As a preliminary step in our further we now look at a special situation for elliptic.
More informationEigenvalue Problems in Surface Acoustic Wave Filter Simulations
Eigenvalue Problems in Surface Acoustic Wave Filter Simulations S. Zaglmayr, J. Schöberl, U. Langer FWF Start Project Y-192 3D hp-finite Elements: Fast Solvers and Adaptivity JKU + RICAM, Linz in collaboration
More informationAlgebraic Multigrid Methods for the Oseen Problem
Algebraic Multigrid Methods for the Oseen Problem Markus Wabro Joint work with: Walter Zulehner, Linz www.numa.uni-linz.ac.at This work has been supported by the Austrian Science Foundation Fonds zur Förderung
More informationMIXED FINITE ELEMENTS FOR PLATES. Ricardo G. Durán Universidad de Buenos Aires
MIXED FINITE ELEMENTS FOR PLATES Ricardo G. Durán Universidad de Buenos Aires - Necessity of 2D models. - Reissner-Mindlin Equations. - Finite Element Approximations. - Locking. - Mixed interpolation or
More informationMath 1080: Numerical Linear Algebra Chapter 4, Iterative Methods
Math 1080: Numerical Linear Algebra Chapter 4, Iterative Methods M. M. Sussman sussmanm@math.pitt.edu Office Hours: MW 1:45PM-2:45PM, Thack 622 March 2015 1 / 70 Topics Introduction to Iterative Methods
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 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 informationFinite Element Methods for Optical Device Design
DFG Research Center Matheon mathematics for key technologies and Zuse Institute Berlin Finite Element Methods for Optical Device Design Frank Schmidt Sven Burger, Roland Klose, Achim Schädle, Lin Zschiedrich
More informationPolynomial-degree-robust a posteriori estimates in a unified setting for conforming, nonconforming, discontinuous Galerkin, and mixed discretizations
Polynomial-degree-robust a posteriori estimates in a unified setting for conforming, nonconforming, discontinuous Galerkin, and mixed discretizations Alexandre Ern, Martin Vohralík To cite this version:
More informationarxiv: v1 [math.na] 4 Jun 2014
MIXE GENERALIZE MULTISCALE FINITE ELEMENT METHOS AN APPLICATIONS ERIC T. CHUNG, YALCHIN EFENIEV, AN CHA SHING LEE arxiv:1406.0950v1 [math.na] 4 Jun 2014 Abstract. In this paper, we present a mixed Generalized
More informationIntroduction. Math 1080: Numerical Linear Algebra Chapter 4, Iterative Methods. Example: First Order Richardson. Strategy
Introduction Math 1080: Numerical Linear Algebra Chapter 4, Iterative Methods M. M. Sussman sussmanm@math.pitt.edu Office Hours: MW 1:45PM-2:45PM, Thack 622 Solve system Ax = b by repeatedly computing
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 informationBUBBLE STABILIZED DISCONTINUOUS GALERKIN METHOD FOR PARABOLIC AND ELLIPTIC PROBLEMS
BUBBLE STABILIZED DISCONTINUOUS GALERKIN METHOD FOR PARABOLIC AND ELLIPTIC PROBLEMS ERIK BURMAN AND BENJAMIN STAMM Abstract. In this paper we give an analysis of a bubble stabilized discontinuous Galerkin
More informationCOMPLEMENTARITY THE WAY TOWARDS GUARANTEED ERROR ESTIMATES
Preprint, Institute of Mathematics, AS CR, Prague. 2010-12-27 INSTITUTE of MATHEMATICS Academy of Sciences Czech Republic COMPLEMENTARITY THE WAY TOWARDS GUARANTEED ERROR ESTIMATES Tomáš Vejchodský Abstract
More informationIMPLEMENTATION OF FUNCTIONAL-TYPE A POSTERIORI ERROR ESTIMATES FOR LINEAR PROBLEMS IN SOLID MECHANICS. Maksim Frolov
IMPLEMENTATION OF FUNCTIONAL-TYPE A POSTERIORI ERROR ESTIMATES FOR LINEAR PROBLEMS IN SOLID MECHANICS Maksim Frolov Peter the Great St.Petersburg Polytechnic University AANMPDE11, FINLAND 06 10 August
More informationFinite element approximation of eigenvalue problems
Acta NumeTica (2010), pp. 1-120 doi: 10.1017/80962492910000012 Cambridge University Press, 2010 Finite element approximation of eigenvalue problems Daniele Boffi Dipartimento di Matematica, Universita
More informationChapter 5 A priori error estimates for nonconforming finite element approximations 5.1 Strang s first lemma
Chapter 5 A priori error estimates for nonconforming finite element approximations 51 Strang s first lemma We consider the variational equation (51 a(u, v = l(v, v V H 1 (Ω, and assume that the conditions
More informationA POSTERIORI ERROR ESTIMATORS FOR CONVECTION-DIFFUSION EIGENVALUE PROBLEMS
A POSERIORI ERROR ESIMAORS FOR CONVECION-DIFFUSION EIGENVALUE PROBLEMS JOSCHA GEDICKE AND CARSEN CARSENSEN Abstract. A posteriori error estimators for convection-diffusion eigenvalue model problems are
More informationConferences Fictitious domain methods for numerical realization of unilateral problems. J.H. (jointly with T. Kozubek and R. Kučera) p.
p. /2 Conferences 28 Fictitious domain methods for numerical realization of unilateral problems J.H. (jointly with T. Kozubek and R. Kučera) Lyon 23-24.6. 28 p. 2/2 Outline Fictitious domain methods -
More informationA computer-assisted Band-Gap Proof for 3D Photonic Crystals
A computer-assisted Band-Gap Proof for 3D Photonic Crystals Henning Behnke Vu Hoang 2 Michael Plum 2 ) TU Clausthal, Institute for Mathematics, 38678 Clausthal, Germany 2) Univ. of Karlsruhe (TH), Faculty
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 informationModular hp-fem System HERMES and Its Application to the Maxwell s Equations 1
Modular hp-fem System HERMES and Its Application to the Maxwell s Equations 1 Tomáš Vejchodský a Pavel Šolín b Martin Zítka b a Mathematical Institute, Academy of Sciences, Žitná 25, 11567 Praha 1, Czech
More informationA NONCONFORMING PENALTY METHOD FOR A TWO DIMENSIONAL CURL-CURL PROBLEM
A NONCONFORMING PENALTY METHOD FOR A TWO DIMENSIONAL CURL-CURL PROBLEM SUSANNE C. BRENNER, FENGYAN LI, AND LI-YENG SUNG Abstract. A nonconforming penalty method for a two-dimensional curl-curl problem
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 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 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 informationb i (x) u + c(x)u = f in Ω,
SIAM J. NUMER. ANAL. Vol. 39, No. 6, pp. 1938 1953 c 2002 Society for Industrial and Applied Mathematics SUBOPTIMAL AND OPTIMAL CONVERGENCE IN MIXED FINITE ELEMENT METHODS ALAN DEMLOW Abstract. An elliptic
More 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 informationAn Introduction to the Discontinuous Galerkin Method
An Introduction to the Discontinuous Galerkin Method Krzysztof J. Fidkowski Aerospace Computational Design Lab Massachusetts Institute of Technology March 16, 2005 Computational Prototyping Group Seminar
More informationBasic Aspects of Discretization
Basic Aspects of Discretization Solution Methods Singularity Methods Panel method and VLM Simple, very powerful, can be used on PC Nonlinear flow effects were excluded Direct numerical Methods (Field Methods)
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 information. D CR Nomenclature D 1
. D CR Nomenclature D 1 Appendix D: CR NOMENCLATURE D 2 The notation used by different investigators working in CR formulations has not coalesced, since the topic is in flux. This Appendix identifies the
More informationMixed Hybrid Finite Element Method: an introduction
Mixed Hybrid Finite Element Method: an introduction First lecture Annamaria Mazzia Dipartimento di Metodi e Modelli Matematici per le Scienze Applicate Università di Padova mazzia@dmsa.unipd.it Scuola
More informationVirtual Element Methods for general second order elliptic problems
Virtual Element Methods for general second order elliptic problems Alessandro Russo Department of Mathematics and Applications University of Milano-Bicocca, Milan, Italy and IMATI-CNR, Pavia, Italy Workshop
More informationIII. Iterative Monte Carlo Methods for Linear Equations
III. Iterative Monte Carlo Methods for Linear Equations In general, Monte Carlo numerical algorithms may be divided into two classes direct algorithms and iterative algorithms. The direct algorithms provide
More informationNumerical Programming I (for CSE)
Technische Universität München WT 1/13 Fakultät für Mathematik Prof. Dr. M. Mehl B. Gatzhammer January 1, 13 Numerical Programming I (for CSE) Tutorial 1: Iterative Methods 1) Relaxation Methods a) Let
More informationMatrix Iteration. Giacomo Boffi.
http://intranet.dica.polimi.it/people/boffi-giacomo Dipartimento di Ingegneria Civile Ambientale e Territoriale Politecnico di Milano April 12, 2016 Outline Second -Ritz Method Dynamic analysis of MDOF
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 information