Adaptive Cell-centered Finite Volume Method
|
|
- Alexia Hart
- 5 years ago
- Views:
Transcription
1 ASC Report No. 0/2008 Adaptive Cell-centered Finite Volume Method Christoph rath, Stefan Funken, Dirk Praetorius Institute for Analysis and Scientific Computing Vienna University of echnology U Wien ISBN
2 Most recent ASC Reports 2/2007 Othmar Koch, Roswitha März, Dirk Praetorius, wa B. Weinmüller Collocation for solving DAs with Singularities /2007 Anton Baranov, Harald Woracek Finite-Dimensional de Branges Subspaces Generated by Majorants 0/2007 Anton Baranov, Harald Woracek Admissible Majorants for de Branges Spaces of ntire Functions 29/2007 Vyacheslav Pivovarchik, Harald Woracek Shifted Hermite-Biehler Functions and their Applications 28/2007 Michael Kaltenbäck, Harald Woracek Canonical Differential quations of Hilbert-Schmidt ype 27/2007 Henrik Winkler, Harald Woracek On Semibounded Canonical Systems 26/2007 Irena Rachunková, Gernot Pulverer, wa B. Weinmüller A Unified Approach to Singular Problems Arising in the Membrane heory 2/2007 Roberta Bosi, Jean Dolbeault, Maria J. steban stimates for the Optimal Constants in Multipolar Hardy Inequalities for Schrödinger and Dirac Operators 2/2007 Xavier Antoine, Anton Arnold, Christophe Besse, Matthias hrhardt, Achim Schädle A Review of ransparent and Artificial Boundary Conditions echniques for Linear and Nonlinear Schrödinger quations 2/2007 ino ibner, Jens Markus Melenk p-fm Quadrature rror Analysis on etrahedra Institute for Analysis and Scientific Computing Vienna University of echnology Wiedner Hauptstraße Wien, Austria -Mail: admin@asc.tuwien.ac.at WWW: FAX: ISBN c Alle Rechte vorbehalten. Nachdruck nur mit Genehmigung des Autors. ASC U WIN
3 ADAPIV CLL-CNRD FINI VOLUM MHOD CHRISOPH RAH, SFAN FUNKN, AND DIRK PRAORIUS Abstract. In our talk, we propose an adaptive mesh-refining strategy for the cell-centered FVM based on some a posteriori error control for the quantity (u Iu h ) L 2. Here, u h P 0 ( ) denotes the FVM approximation of u and I is a certain interpolation operator. As model example serves the Laplace equation with mixed boundary conditions, where our contributions extend a result of [NIC0]. Moreover, this approach allows the coupling of finite volume schemes with the boundary element method, which is a new and fruitful combination of the FVM with ideas from [CAR99a, CAR99b].. Introduction and lliptic Model Problem Let Ω R 2 be a bounded and connected domain with Lipschitz boundary Γ := Ω, which is divided into a closed Dirichlet boundary Γ D Γ with positive surface measure and a Neumann boundary Γ N := Γ\Γ D. We consider the elliptic boundary value problem (.) u = f in Ω, u = u D on Γ D and u/ n = g on Γ N. Here, f L 2 (Ω), u D H (Γ D ), and g L 2 (Γ N ) are given data, and L 2 ( ) and H ( ) denote the standard Lebesgue- and Sobolev-spaces equipped with the usual norms L 2 ( ) and H ( ). We aim to approximate the unique weak solution u H (Ω) by a postprocessed finite volume scheme. hroughout, denotes a certain triangulation of Ω, where N and are the corresponding sets of nodes and edges, respectively. he set of nodes (edges) on the Dirichlet resp. Neumann boundary are N D resp. N N ( D resp. N ). For brevity, we assume that the elements are non-degenerate rectangles and refer to [RA07] for triangular elements. For, h := diam() denotes the uclidean diameter, and for an edge, we denote by h its length. We say that the triangulation is almost regular, if (i) the mixed boundary conditions are resolved, i.e. each edge with Γ satisfies either Γ D or Γ N. (ii) the intersection 2 of two elements, 2 with 2 is either empty or a node or an edge. (iii) the edge contains an interior (i.e. hanging) node, there are two edges, 2 with = 2, cf. Figure 2. for examples. A node a N \(N D N N ) is a hanging node provided that there are elements, 2 such that a 2 is a node of but not of 2. Let N H be the set of all hanging nodes. Date: January, Key words and phrases. finite volume method, cell-centered method, diamond path, a posteriori error estimate, adaptive algorithm.
4 Figure 2.. he weights ψ for an almost regular mesh of squares. 2. Cell-Centered Finite Volume Method We integrate the differential equation (.) over a control volume and use the Gauss divergence theorem to obtain, with the set of edges of, u f dx = ds = Φ, (u) ds for all. n Here, Φ, (u) = u/ n, ds is the diffusive flux and n, is the outer normal vector of on. Let u h P 0 ( ) be a piecewise constant approximation of u, namely u := u h u(x ), where x denotes the center of an element. For the cell-centered finite volume method, u h is computed by replacing the diffusion flux Φ, (u) by a discrete diffusion flux F, (u h ). First, Φ, (u) = g ds is known for a Neumann edge N. One therefore defines F, (u h ) := Φ, (u h ) = g ds for N. Second, for a non-elementary edge with = 2 and, 2, where denotes the set of elementary edges (i.e. neither nor 2 contain an interior node), there holds Φ, (u) = Φ, (u) + Φ,2 (u), which leads to the definition F, (u h ) := F, (u h ) + F,2 (u h ) for = 2 with, 2. Finally, it remains to define F, (u h ) for elementary and Dirichlet edges D. his is done by the diamond path method in the following way: For each node a N, we define ψ (a) u, for all a N \(N D N N ), eω a (2.) u a = u D (a), for all a N D, u a + g a, for all a N N, with certain weights { ψ (a), a N } for each element of the patch ωa := { a }. Here, N denotes the set of nodes of. For details on the computation of the weights, the reader is referred to [COU00]. Figure 2. gives the precise values for almost regular triangulations of squares. Details on the computation of u a and g a are found in [RA07]. For an elementary edge, let x a and x b be the starting and end point of 2
5 x b x d = (x x ) n h t, n, (x x ) t x x a Figure 2.2. Notation for the diamond path. D and, the unique elements with =, cf. Figure 2.2. hen, ( ) u u u b u a F, (u h ) := h α d h (2.2) with α = (x x ) t, (x x ) n,, d = (x x ) n,. Here, the additional unknowns u b and u a are located at the nodes x b and x a and are computed by (2.). Finally, the tangential vector t, is chosen orthogonal to n, in mathematically positive sense. For a Dirichlet edge D, we compute F, (u h ) by (2.2), where x is now replaced by the midpoint x m of and u becomes u D (x m ). Altogether, the discrete problem reads: Find u h P 0 ( ) such that F, (u h ) = f dx, for all. Note, that the conservativity of the continuous flux Φ, (u) = Φ,(u) also holds for the discrete flux F, (u h ) = F,(u h ). Remark. We stress, that even for an admissible triangular mesh in the sense of [YM00, Definition 9.], local mesh-refinement is nontrivial, since admissibility necessarily implies that all angles are strictly less than π/2. For rectangular meshes, local mesh-refinement cannot avoid hanging nodes and thus contradicts the admissibility condition.. A Posteriori rror stimate For a non-degenerate rectangular element = conv{a, a 2, a, a } R 2 with edges j = conv{a j, a j+ } and a := a, we define P = P 2 span{x xy 2, y yx 2 } and Σ = (S,...,S 8 ), where S j (p) = p(a j ), S j+ (p) = j p n,j ds for j =,..., (p P ), cf. [NIC0, Section.2]. hen, the Morley-type element (, P, Σ ) is a nonconforming finite element. he corresponding Morley interpolant Iu h is now defined -elementwise by (.) (.). he definition of which is an extension of the definition in [NIC0, Section
6 ] to the case of hanging nodes and mixed boundary conditions. For each free node a N ( N \(N D N N N H ) ), we enforce (.) (Iu h ) (a) = ψ a (a)u h a. a eω a For each boundary node, the value Iu h (a) is prescribed by { u D (a) for a N N D, (.2) (Iu h ) (a) = u a + g a for a N N N. For each hanging node a N N H, there holds (.) (Iu h ) (a) = (Iu h ) a (a), where a is the unique element with a int() for some (non-elementary) edge a. Finally, for each edge holds (Iu h ) (.) ds = F, (u h ). n, We stress, that the Morley interpolant Iu h is uniquely defined by (.) (.). Moreover, the definition ensures the following orthogonality properties, which are essential for the analysis of the proposed error estimator: First, the residual R := f + (Iu h ) is L 2 -orthogonal to P 0 ( ), i.e. ( f + (Iuh ) ) dx = 0 for all. In particular, R = f f, where f denotes the -piecewise integral mean, i.e. f := f dx. Second, for boundary edges hold (u Iu h ) (u Iu h ) ds = 0 ( D ) resp. ds = 0 ( N ). t, n, Finally, for interior edges hold [ (Iuh ) ] ds = 0 ( 0 ) n, resp. [ (Iuh ) ] ds = 0 ( ). t, Here, 0 denotes the set of all interior edges which are not part of a non-elementary edge. o introduce our error estimator, we define the normal jump over an edge by [ (Iu h )/ n, ] = (Iu h )/ n, + (Iu h )/ n,, where, and =. Note that n, = n, so that the sum in the definition is in fact a difference. For each non-elementary edge = 2 with, 2, we define the normal jump [ (Iu h )/ n, ] (x) := [ (Iu h )/ n, ] i (x) for all x i, i =, 2.
7 0 0 error and error estimator I (unif.) η (unif.) h (unif.) I (adap.) η (adap.) h (adap.) / / number of elements Figure.. rrors I and h as well as a posteriori error estimator η in Problem (.) for uniform and adaptive mesh-refinement. he tangential jump [ (Iu h )/ t, ] is defined analogously. For each element, we now define the refinement indicator η 2 := h2 f f 2 L 2 () + h [ (Iu h )] 2 L 2 () (.) + N h (u Iu h) 2 n, \( D N ) L 2 () + D h (u Iu h) 2 t, L 2 () With techniques known from the finite element method [AIN00], one proves reliability ( C rel /2. (.6) (u Iu h ) L 2 (Ω) η := η) 2 As for non-conforming FM the proof employs a Helmholtz decomposition. he Galerkin orthogonality is replaced by the orthogonality properties of Iu h stated before. he converse inequality holds even -elementwise, namely for all (.7) C loc η ( (u Iu h ) 2 L 2 (ω ) + h2 f f 2 L 2 (ω )) /2, where the involved patch reads ω := { }. his estimate follows by usual techniques and the use of appropriate bubble functions. In particular, we obtain efficiency of η up to oscillation terms, (.8) C eff η (u Iu h ) L 2 (Ω) + h(f f ) L 2 (Ω), where h L (Ω), h := h denotes the local mesh-width. he constants C rel, C loc, C eff > 0 only depend on the shape of the elements in. We refer to [RA07] for the detailed proofs of the estimates (.6) (.8)..
8 # = 69 # = 67 # = 629 Figure.2. Adaptively generated meshes in problem (.).. Numerical xperiments We consider the Laplace problem (.) on the L-shaped domain Ω = (, ) 2 \ ( [0, ] [, 0] ). he given exact solution is the harmonic function u(x, y) = I ( (x + iy) 2/) which reads in polar coordinates u(x, y) = r 2/ sin(2ϕ/) with (x, y) = r (cosϕ, sin ϕ). Note that u has a generic singularity at the reentrant corner (0, 0). For the numerical computation, we prescribe the exact Neumann and Dirichlet data, where Γ D = Γ\Γ N and Γ N := {0} (, 0) (0, ) {0}. Note that Γ N includes the reentrant corner, where the normal derivative u/ n is singular. We compute the energy error h := ( u u h 2 L 2 (Ω) + /2, u u h,h) 2 with u P 0 ( ) the -piecewise integral mean of u, and where the discrete H -seminorm is defined by ( 2 ) /2 v h,h = (v v )/d h d D for any v h P 0 ( ). Provided that u H 2 (Ω), the diamond path method leads to h = O(N /2 ) with respect to the number N = of elements [COU00]. Moreover, we compute the Morley error I = (u Iu h ) L 2 (Ω) and the corresponding error estimator η. All computations are done in Matlab in the spirit of [ALB99, RA08]. We employ an adaptive mesh-refining algorithm, which is steered by the local refinement indicators η from (.): An element j is marked for refinement provided η j θ max{η,...,η N }, where we use θ = 0 for uniform and θ = 0. for adaptive mesh-refinement, respectively. he numerical outcome is plotted in Figure. over the number N of elements. For uniform mesh-refinement, the energy error h converges with a suboptimal order which appears to be slightly better than O(N / ). he proposed adaptive strategy regains the optimal order of convergence O(N /2 ). As can be expected from the finite element method, the Morley error I decreases like O(N / ) for uniform mesh-refinement. he adaptive algorithm leads to an improved order of convergence O(N /2 ). As predicted by theory, the error estimator η is observed to be reliable and efficient.. Outlook he finite volume method is a well adapted method for the discretization of various convection dominated partial differential equations. It is very popular in the engineering community (fluid mechanics) because of its conservative properties. his method is contrary to 6
9 the boundary element method (BM), which can be applied to the most important linear and homogeneous partial differential equations with constant coefficients also in unbounded domains. he coupling of FVM and BM combines the advantages of both methods, e.g. in problems where stationary diffusive heat and convection in different media are coupled. While convection would be modelled by the finite volume method, diffusive heat (in a possibly unbounded domain) is solved using the boundary element method. First numerical examples show the efficiency of the symmetric coupling of FVM and BM, where we used the results given here for the FVM and results from [CAR99a, CAR99b] to develop a stable discretization. he resulting system of linear equations is now a block-system instead of a system for the FM NC -BM coupling. he coupling of FVM and BM involves two further continuous ansatz functions on the interface to link the discontinuous displacement field to necessarily continuous boundary ansatz functions on the boundary. Quasi-optimal a priori error estimates and sharp a posteriori error estimates are almost established which justify adaptive mesh-refining algorithms. Numerical experiments show the adaptive coupling as an efficient tool for the numerical treatment of transmission problems. References [AIN00] M. Ainsworth, J.. Oden: A posteriori error estimation in finite element analysis, John Wiley & Sons, [ALB99] J. Alberty, C. Carstensen, S.A. Funken: Remarks around 0 lines of Matlab: Short finite element implementation, Num. Alg., 20:7 7, 999. [CAR99a] C. Carstensen, S. A. Funken: Coupling of nonconforming finite elements and boundary elements I: A priori estimates, Computing, 62:229 2, 999. [CAR99b] C. Carstensen, S. A. Funken: Coupling of nonconforming finite elements and boundary elements II: A posteriori estimates and adaptive mesh-refinement, Computing, 62:2 29, 999. [CIA78] P.G. Ciarlet: he finite element method for elliptic problems, North-Holland, Amsterdam, 978. [COU00] Y. Coudiére, P. Villedieu: Convergence rate of a finite volume scheme for the linear convection-diffusion equation on locally refined meshes, SAIM: M2AN, (6):2 9, [RA08] C. rath, S. Funken. D. Praetorius: Matlab implementation of the finite volume method, part II: he cell centered FVM, in preparation, [RA07] C. rath, D. Praetorius: A posteriori error estimate and adaptive mesh-refinement for the cell-centered finite volume method for elliptic boundary value problems, submitted [YM00] R. ymard,. Gallouët, R. Herbin: Finite volume methods, Handbook of Numerical Analysis, Volume 7. lsevier Science B.V., first edition, 999. [NIC0] S. Nicaise: A posteriori error estimations of some cell-centered finite volume methods, SIAM J. Numer. Anal., (0):8 0, 200. University of Ulm, Institute for Numerical Mathematics, Helmholtzstraße 8, D Ulm, Germany -mail address: {Christoph.rath,Stefan.Funken}@uni-ulm.de Vienna University of echnology, Institute for Analysis and Scientific Computing, Wiedner Hauptstraße 8-0, A-00 Vienna, Austria -mail address: Dirk.Praetorius@tuwien.ac.at 7
Adaptive Cell-centered Finite Volume Method
ASC Report No. 0/2008 Adaptive Cell-centered Finite Volume Method Christoph rath, Stefan Funken, Dirk Praetorius Institute for Analysis and Scientific Computing Vienna University of Technology TU Wien
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 informationDownloaded 07/25/13 to Redistribution subject to SIAM license or copyright; see
SIAM J. NUMR. ANAL. Vol. 39, No. 6, pp. 2034 2044 c 2002 Society for Industrial and Applied Mathematics RSIDUAL-BASD A POSRIORI RROR SIMA FOR A NONCONFORMING RISSNR MINDLIN PLA FINI LMN CARSN CARSNSN Abstract.
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 informationOn Semibounded Canonical Systems
ASC Report No. 27/2007 On Semibounded Canonical Systems Henrik Winkler, Harald Woracek Institute for Analysis and Scientific Computing Vienna University of Technology TU Wien www.asc.tuwien.ac.at ISBN
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 2018/2019. Lecture 15. Jürgen Fuhrmann Lecture 15 Slide 1
Scientific Computing WS 2018/2019 Lecture 15 Jürgen Fuhrmann juergen.fuhrmann@wias-berlin.de Lecture 15 Slide 1 Lecture 15 Slide 2 Problems with strong formulation Writing the PDE with divergence and gradient
More 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 informationTechnische Universität Graz
Technische Universität Graz A non-symmetric coupling of the finite volume method and the boundary element method C. Erath, G. Of, F. J. Sayas Berichte aus dem Institut für Numerische Mathematik Bericht
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 informationc 2007 Society for Industrial and Applied Mathematics
SIAM J. NUMR. ANAL. Vol. 45, No. 1, pp. 68 82 c 2007 Society for Industrial and Applied Mathematics FRAMWORK FOR TH A POSTRIORI RROR ANALYSIS OF NONCONFORMING FINIT LMNTS CARSTN CARSTNSN, JUN HU, AND ANTONIO
More informationWhispering gallery modes in oblate spheroidal cavities: calculations with a variable stepsize
ASC Report No. 37/14 Whispering gallery modes in oblate spheroidal cavities: calculations with a variable stepsize P. Amodio, T. Levitina, G. Settanni, E. Weinmüller Institute for Analysis and Scientic
More informationSuperconvergence and H(div) Projection for Discontinuous Galerkin Methods
INTRNATIONAL JOURNAL FOR NUMRICAL MTHODS IN FLUIDS Int. J. Numer. Meth. Fluids 2000; 00:1 6 [Version: 2000/07/27 v1.0] Superconvergence and H(div) Projection for Discontinuous Galerkin Methods Peter Bastian
More informationBasic Concepts of Adaptive Finite Element Methods for Elliptic Boundary Value Problems
Basic Concepts of Adaptive Finite lement Methods for lliptic Boundary Value Problems Ronald H.W. Hoppe 1,2 1 Department of Mathematics, University of Houston 2 Institute of Mathematics, University of Augsburg
More informationA mixed 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 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 informationTechnische Universität Graz
Technische Universität Graz Stability of the Laplace single layer boundary integral operator in Sobolev spaces O. Steinbach Berichte aus dem Institut für Numerische Mathematik Bericht 2016/2 Technische
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 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 informationENERGY NORM A POSTERIORI ERROR ESTIMATES FOR MIXED FINITE ELEMENT METHODS
ENERGY NORM A POSTERIORI ERROR ESTIMATES FOR MIXED FINITE ELEMENT METHODS CARLO LOVADINA AND ROLF STENBERG Abstract The paper deals with the a-posteriori error analysis of mixed finite element methods
More informationA 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 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 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 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 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 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 informationEfficient Implementation of Adaptive P1-FEM in MATLAB
ASC Report No. 19/2008 Efficient Implementation of Adaptive P1-FEM in MATLAB Stefan Funken, Dirk Praetorius, Philipp Wissgott Institute for Analysis and Scientific Computing Vienna University of Technology
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 informationWavenumber-explicit hp-bem for High Frequency Scattering
ASC Report No. 2/2010 Wavenumber-explicit hp-bem for High Frequency Scattering Maike Löhndorf, Jens Markus Melenk Institute for Analysis and Scientific Computing Vienna University of Technology TU Wien
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 information1. Introduction. This paper is devoted to the numerical treatment of the hypersingular
SIAM J. SCI. COMPUT. Vol. 29, No. 2, pp. 782 80 c 2007 Society for Industrial and Applied Mathematics AVERAGING TECNIQUES FOR TE A POSTERIORI BEM ERROR CONTROL FOR A YPERSINGULAR INTEGRAL EQUATION IN TWO
More informationNonconforming FEM for the obstacle problem
IMA Journal of Numerical Analysis (2017) 37, 64 93 doi: 10.1093/imanum/drw005 Advance Access publication on May 10, 2016 Nonconforming FM for the obstacle problem C. Carstensen and K. Köhler Institut für
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 informationHigh order, finite volume method, flux conservation, finite element method
FLUX-CONSERVING FINITE ELEMENT METHODS SHANGYOU ZHANG, ZHIMIN ZHANG, AND QINGSONG ZOU Abstract. We analyze the flux conservation property of the finite element method. It is shown that the finite element
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 informationMapping Properties of Combined Field Helmholtz Boundary Integral Operators
ASC Report No. 1/2010 Mapping Properties of Combined Field Helmholtz Boundary Integral Operators Jens Markus Melenk Institute for Analysis and Scientific Computing Vienna University of Technology TU Wien
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 informationc 2006 Society for Industrial and Applied Mathematics
SIAM J. SCI. COMPUT. Vol. 27, No. 4, pp. 226 260 c 2006 Society for Industrial and Applied Mathematics AVERAGING TECNIQUES FOR TE EFFECTIVE NUMERICAL SOLUTION OF SYMM S INTEGRAL EQUATION OF TE FIRST KIND
More informationDifferent Approaches to a Posteriori Error Analysis of the Discontinuous Galerkin Method
WDS'10 Proceedings of Contributed Papers, Part I, 151 156, 2010. ISBN 978-80-7378-139-2 MATFYZPRESS Different Approaces to a Posteriori Error Analysis of te Discontinuous Galerkin Metod I. Šebestová Carles
More informationALL FIRST-ORDER AVERAGING TECHNIQUES FOR A POSTERIORI FINITE ELEMENT ERROR CONTROL ON UNSTRUCTURED GRIDS ARE EFFICIENT AND RELIABLE
MATHEMATICS OF COMPUTATION Volume 73, Number 247, Pages 1153 1165 S 0025-5718(03)01580-1 Article electronically published on August 12, 2003 ALL FIRST-ORDER AVERAGING TECHNIQUES FOR A POSTERIORI FINITE
More 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 informationTechnische Universität Graz
Technische Universität Graz A note on the stable coupling of finite and boundary elements O. Steinbach Berichte aus dem Institut für Numerische Mathematik Bericht 2009/4 Technische Universität Graz A
More informationSECOND ORDER TIME DISCONTINUOUS GALERKIN METHOD FOR NONLINEAR CONVECTION-DIFFUSION PROBLEMS
Proceedings of ALGORITMY 2009 pp. 1 10 SECOND ORDER TIME DISCONTINUOUS GALERKIN METHOD FOR NONLINEAR CONVECTION-DIFFUSION PROBLEMS MILOSLAV VLASÁK Abstract. We deal with a numerical solution of a scalar
More informationGeneralized DG-Methods for Highly Indefinite Helmholtz Problems based on the Ultra-Weak Variational Formulation
ASC Report No. 06/2012 Generalized DG-Methods for Highly Indefinite Helmholtz Problems based on the Ultra-Weak Variational Formulation J.M. Melenk, A. Parsania, and S. Sauter Institute for Analysis and
More informationThe Pole Condition: A Padé Approximation of the Dirichlet to Neumann Operator
The Pole Condition: A Padé Approximation of the Dirichlet to Neumann Operator Martin J. Gander and Achim Schädle Mathematics Section, University of Geneva, CH-, Geneva, Switzerland, Martin.gander@unige.ch
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 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 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 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 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 informationAutomatic Grid Control in Adaptive BVP Solvers
ASC Report o. /28 Automatic Grid Control in Adaptive BVP Solvers Gernot Pulverer, Gustaf Söderlind, Ewa Weinmüller Institute for Analysis and Scientific Computing Vienna University of Technology TU Wien
More informationIt is known that Morley element is not C 0 element and it is divergent for Poisson equation (see [6]). When Morley element is applied to solve problem
Modied Morley Element Method for a ourth Order Elliptic Singular Perturbation Problem Λ Wang Ming LMAM, School of Mathematical Science, Peking University Jinchao u School of Mathematical Science, Peking
More informationAveraging technique for FE ± a posteriori error control in elasticity. Part II: k-independent estimates
Comput. Methods Appl. Mech. Engrg. 190 2001) 4663±4675 www.elsevier.com/locate/cma Averaging technique for FE ± a posteriori error control in elasticity. Part II: k-independent estimates Carsten Carstensen,
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 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 informationWEAK GALERKIN FINITE ELEMENT METHODS ON POLYTOPAL MESHES
INERNAIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume 12, Number 1, Pages 31 53 c 2015 Institute for Scientific Computing and Information WEAK GALERKIN FINIE ELEMEN MEHODS ON POLYOPAL MESHES LIN
More informationA Balancing Algorithm for Mortar Methods
A Balancing Algorithm for Mortar Methods Dan Stefanica Baruch College, City University of New York, NY 11, USA Dan Stefanica@baruch.cuny.edu Summary. The balancing methods are hybrid nonoverlapping Schwarz
More informationarxiv: v1 [math.na] 19 Nov 2018
QUASINORMS IN SEMILINEAR ELLIPTIC PROBLEMS JAMES JACKAMAN AND TRISTAN PRYER arxiv:1811.07816v1 [math.na] 19 Nov 2018 (1.1) Abstract. In this note we examine the a priori and a posteriori analysis of discontinuous
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 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 informationABHELSINKI UNIVERSITY OF TECHNOLOGY
ABHELSINKI UNIVERSITY OF TECHNOLOGY TECHNISCHE UNIVERSITÄT HELSINKI UNIVERSITE DE TECHNOLOGIE D HELSINKI A posteriori error analysis for the Morley plate element Jarkko Niiranen Department of Structural
More 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 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 informationA local-structure-preserving local discontinuous Galerkin method for the Laplace equation
A local-structure-preserving local discontinuous Galerkin method for the Laplace equation Fengyan Li 1 and Chi-Wang Shu 2 Abstract In this paper, we present a local-structure-preserving local discontinuous
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 informationMedius analysis and comparison results for first-order finite element methods in linear elasticity
IMA Journal of Numerical Analysis Advance Access published November 7, 2014 IMA Journal of Numerical Analysis (2014) Page 1 of 31 doi:10.1093/imanum/dru048 Medius analysis and comparison results for first-order
More informationA posteriori error estimation of approximate boundary fluxes
COMMUNICATIONS IN NUMERICAL METHODS IN ENGINEERING Commun. Numer. Meth. Engng 2008; 24:421 434 Published online 24 May 2007 in Wiley InterScience (www.interscience.wiley.com)..1014 A posteriori error estimation
More informationFinite volume method for nonlinear transmission problems
Finite volume method for nonlinear transmission problems Franck Boyer, Florence Hubert To cite this version: Franck Boyer, Florence Hubert. Finite volume method for nonlinear transmission problems. Proceedings
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 informationarxiv: v1 [math.na] 29 Feb 2016
EFFECTIVE IMPLEMENTATION OF THE WEAK GALERKIN FINITE ELEMENT METHODS FOR THE BIHARMONIC EQUATION LIN MU, JUNPING WANG, AND XIU YE Abstract. arxiv:1602.08817v1 [math.na] 29 Feb 2016 The weak Galerkin (WG)
More 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 informationSome remarks on the history and future of averaging techniques in a posteriori finite element error analysis
ZAMM Z. Angew. Math. Mech. 84, No. 1, 3 21 (2004) / DOI 10.1002/amm.200410101 Some remarks on the history and future of averaging techniques in a posteriori finite element error analysis Plenary lecture
More informationIntroduction. Finite and Spectral Element Methods Using MATLAB. Second Edition. C. Pozrikidis. University of Massachusetts Amherst, USA
Introduction to Finite and Spectral Element Methods Using MATLAB Second Edition C. Pozrikidis University of Massachusetts Amherst, USA (g) CRC Press Taylor & Francis Group Boca Raton London New York CRC
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 informationDerivation and Analysis of Piecewise Constant Conservative Approximation for Anisotropic Diffusion Problems
Derivation Analysis of Piecewise Constant Conservative Approximation for Anisotropic Diffusion Problems A. Agouzal, Naïma Debit o cite this version: A. Agouzal, Naïma Debit. Derivation Analysis of Piecewise
More informationDiscontinuous Galerkin Methods: An Overview and Some Applications. Daya Reddy UNIVERSITY OF CAPE TOWN
Discontinuous Galerkin Methods: An Overview and Some Applications Daya Reddy UNIVRSITY OF CAP TOWN Joint work with Beverley Grieshaber and Andrew McBride SANUM, Stellenbosch, 22 24 March 2016 Structure
More informationWeierstraß-Institut. für Angewandte Analysis und Stochastik. im Forschungsverbund Berlin e.v. Preprint ISSN
Weierstraß-Institut für Angewandte Analysis und Stochastik im Forschungsverbund Berlin e.v. Preprint ISSN 946 8633 Computational comparison between the Taylor Hood and the conforming Crouzeix Raviart element
More informationWeierstraß Institut. für Angewandte Analysis und Stochastik. im Forschungsverbund Berlin e.v. Preprint ISSN
Weierstraß Institut für Angewandte Analysis und Stochastik im Forschungsverbund Berlin e.v. Preprint ISSN 0946 8633 Discrete Sobolev-Poincaré inequalities for Voronoi finite volume approximations Annegret
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 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 informationAnalysis of an Adaptive Finite Element Method for Recovering the Robin Coefficient
Analysis of an Adaptive Finite Element Method for Recovering the Robin Coefficient Yifeng Xu 1 Jun Zou 2 Abstract Based on a new a posteriori error estimator, an adaptive finite element method is proposed
More 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 informationSimple Examples on Rectangular Domains
84 Chapter 5 Simple Examples on Rectangular Domains In this chapter we consider simple elliptic boundary value problems in rectangular domains in R 2 or R 3 ; our prototype example is the Poisson equation
More informationA LAGRANGE MULTIPLIER METHOD FOR ELLIPTIC INTERFACE PROBLEMS USING NON-MATCHING MESHES
A LAGRANGE MULTIPLIER METHOD FOR ELLIPTIC INTERFACE PROBLEMS USING NON-MATCHING MESHES P. HANSBO Department of Applied Mechanics, Chalmers University of Technology, S-4 96 Göteborg, Sweden E-mail: hansbo@solid.chalmers.se
More informationPreconditioned space-time boundary element methods for the heat equation
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 Preconditioned space-time boundary element methods for the heat equation S. Dohr and O. Steinbach Institut für Numerische Mathematik Space-Time Methods
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 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 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 informationLECTURE 1: SOURCES OF ERRORS MATHEMATICAL TOOLS A PRIORI ERROR ESTIMATES. Sergey Korotov,
LECTURE 1: SOURCES OF ERRORS MATHEMATICAL TOOLS A PRIORI ERROR ESTIMATES Sergey Korotov, Institute of Mathematics Helsinki University of Technology, Finland Academy of Finland 1 Main Problem in Mathematical
More informationWeierstraß-Institut. für Angewandte Analysis und Stochastik. Leibniz-Institut im Forschungsverbund Berlin e. V. Preprint ISSN
Weierstraß-Institut für Angewandte Analysis und Stochastik Leibniz-Institut im Forschungsverbund Berlin e. V. Preprint ISSN 2198-5855 On the divergence constraint in mixed finite element methods for incompressible
More informationTechnische Universität Graz
Technische Universität Graz Robust boundary element domain decomposition solvers in acoustics O. Steinbach, M. Windisch Berichte aus dem Institut für Numerische Mathematik Bericht 2009/9 Technische Universität
More informationOn a Discontinuous Galerkin Method for Surface PDEs
On a Discontinuous Galerkin Method for Surface PDEs Pravin Madhavan (joint work with Andreas Dedner and Bjo rn Stinner) Mathematics and Statistics Centre for Doctoral Training University of Warwick Applied
More informationA P4 BUBBLE ENRICHED P3 DIVERGENCE-FREE FINITE ELEMENT ON TRIANGULAR GRIDS
A P4 BUBBLE ENRICHED P3 DIVERGENCE-FREE FINITE ELEMENT ON TRIANGULAR GRIDS SHANGYOU ZHANG DEDICATED TO PROFESSOR PETER MONK ON THE OCCASION OF HIS 6TH BIRTHDAY Abstract. On triangular grids, the continuous
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 informationThe mortar element method for quasilinear elliptic boundary value problems
The mortar element method for quasilinear elliptic boundary value problems Leszek Marcinkowski 1 Abstract We consider a discretization of quasilinear elliptic boundary value problems by the mortar version
More informationChapter 1: The Finite Element Method
Chapter 1: The Finite Element Method Michael Hanke Read: Strang, p 428 436 A Model Problem Mathematical Models, Analysis and Simulation, Part Applications: u = fx), < x < 1 u) = u1) = D) axial deformation
More informationUnique continuation for the Helmholtz equation using a stabilized finite element method
Unique continuation for the Helmholtz equation using a stabilized finite element method Lauri Oksanen University College London Based on a joint work with Erik Burman and Mihai Nechita Motivation: recovering
More informationNumerische Mathematik
Numer. Math. (203) 24:309 335 DOI 0.007/s002-02-053-5 Numerische Mathematik A posteriori error estimates for nonconforming finite element methods for fourth-order problems on rectangles Carsten Carstensen
More informationANALYSIS OF THE FEM AND DGM FOR AN ELLIPTIC PROBLEM WITH A NONLINEAR NEWTON BOUNDARY CONDITION
Proceedings of EQUADIFF 2017 pp. 127 136 ANALYSIS OF THE FEM AND DGM FOR AN ELLIPTIC PROBLEM WITH A NONLINEAR NEWTON BOUNDARY CONDITION MILOSLAV FEISTAUER, ONDŘEJ BARTOŠ, FILIP ROSKOVEC, AND ANNA-MARGARETE
More informationA Balancing Algorithm for Mortar Methods
A Balancing Algorithm for Mortar Methods Dan Stefanica Baruch College, City University of New York, NY, USA. Dan_Stefanica@baruch.cuny.edu Summary. The balancing methods are hybrid nonoverlapping Schwarz
More information