arxiv: v1 [math.na] 1 May 2013
|
|
- Clinton Sherman
- 6 years ago
- Views:
Transcription
1 arxiv: v [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 on projecting the discrete gradient onto the standard finite element space is considered We use an oblique projection, where the test and trial spaces are different, and the bases of these two spaces form a biorthogonal system Biorthogonality allows efficient computation of the recovery operator We analyse the approximation properties of the gradient recovery operator Key words Gradient reconstruction, oblique projection, biorthogonal system AMS subject classification 65N30, 65N5, 65N50 Introduction The gradient reconstruction is a popular technique to develop a reliable a posteriori error estimators for approximating the solution of partial differential equations using adaptive finite element methods [, 3,, 9, 3, ] Recently we have presented a gradient reconstruction operator based on an oblique projection [5] The oblique projection operator is constructed by using a biorthogonal system In fact, for the linear finite element in simplicial meshes, this approach reproduces the so-called gradient reconstruction scheme by the weighted averaging [4, 9, ] We proved that the approximation property of the recovered gradient for any finite element space is similar to the one obtained by using the orthogonal projection with respect to L -norm [9, 3] In this article, we aim at analysing the approximation property of the recovered gradient in one School of Mathematical and Physical Sciences, University of Newcastle, Callaghan, NSW 308, bishnulamichhane@newcastleeduau, adammcneilly@uoneduau
2 dimension using the oblique projection This construction is quite useful in extending the weighted average gradient recovery of linear finite elements [9] to quadrilaterals and hexahedras Let = (α, β) with α, β R and α < β Let = {α = < x < < x n = β} be a partition of the interval We define the interior of the grid, denoted int ( ), as int ( ) = {x i : i n } We also define the set of intervals in the partition as {I i } n i=0, where I i = [x i, x i+ ) Two sets A n and B n of indices are also defined as A n = {i N : 0 i n} and B n = A n \ {0, n}, respectively A piecewise linear interpolant of a continuous function u is written as I h u V h with I h u(x) = n u(x i )φ i (x), i=0 where φ i is the standard hat function associated with the point x i, 0 i n We define a discrete space, V h = span {φ 0,, φ n } H () The linear interpolant of u H () is the continuous function defined by I h u = n i=0 u(x i)φ i However, if we compute the derivative of this interpolant I h u, the resulting function will not be continuous To make the derivative continuous we project the derivative of the interpolant, u h x = n i=0 u(x i) φi x, onto the discrete space V h There are two different types of projection One is an orthogonal projection and the other is an oblique projection The orthogonal projection operator, P h, that projects u h x onto V u h is to find a g h = P h h x V h that satisfies: u h g h φ j dx = x φ j dx () Since g h V h, we can represent it as an (n + )-dimensional vector: g = g 0 g n with g h = n g i φ i Now the requirement given in equation () is equivalent to a linear system: M g = f, where M is a mass matrix, and f = f 0 f n with f j = i=0 u h x φ j dx
3 Here the mass matrix M is tridiagonal We can reduce computation time greatly if we have a diagonal mass matrix This can be done if we use a suitable oblique projection instead of an orthogonal projection We consider the projection Q h : L () V h, du which is defined as the problem of finding g h = Q h dx V h such that g h λ h dx = dx λ h dx, λ h M h, where M h is another piecewise polynomial space, not orthogonal to V h, with dim (M h ) = dim (V h ), see [8] In fact, the projection operator Q h is well-defined due to the following stability condition There is a constant β > 0 independent of the mesh-size h such that [5, 6] φ h L () β sup µ h M h \{0} µ hφ h dx, φ h V h µ h L () In order to achieve that the mass matrix M is diagonal we need to define a new set of basis functions for M h, {λ 0,, λ n }, that are biorthogonal to the standard hat basis function we used previously This biorthogonality relation is defined as: λ i φ j dx = c j δ ij, c j 0, i, j n, () where δ ij is the Kronecker delta function: {, if i = j δ ij = 0, otherwise, and c j is a positive scaling factor The basis functions for M h are simply given by { (x x)+(x ) λ 0 (x) = x, x x 0, otherwise { (x xn )+(x x n) x λ n (x) = n x n, x n x x n, 0, otherwise and for i n (x x i )+(x x i) x i x i, x i x x i (x x λ i (x) = i+)+(x x i) x i x i+, x i x x i+ 0, otherwise 3
4 Figure : The hat basis function (left) and biorthogonal basis function (right) with stepsize h = By using an oblique projection Q h the mass matrix will be diagonal We let the diagonal mass matrix be D, so that our system is D g = f The values g i are our estimates of the gradient of u at the point x i So, we estimate the gradient by finding g = D f, where g i = dx λ idx φ i λ i dx We want to calculate the error in this approximation, and find out when g i approximates u (x i ) exactly for each x i As in [0, ] we want to see if g i approximates u (x i ) exactly when u is a quadratic polynomial Superconvergence Theorem Let u C 0 () Then we have g i = u(x i+) u(x i ) x i+ x i, i B n, and Proof: g 0 = u(x ) u( ) x, g n = u(x n) u(x n ) x n x n We note that (x x i )+(x x i) x i x i, x i x x i (x x λ i (x) = i+)+(x x i) x i x i+ x i x x i+, for all i B n 0, otherwise Now, we calculate g i for i B n : g i = dx λ i dx φ i λ i dx, 4
5 where φ i λ i dx = + xi x i xi+ x i ( ) ( x xi (x xi ) + (x x i ) x i x i x i x i ( ) ( x xi+ (x xi+ ) + (x x i ) x i x i+ x i x i+ ) dx ) dx = (x i x i+ ), and dx λ i dx = n j=0 = u(x i ) dφ j u(x j ) dx λ i dx xi x i dφ i dx λ i dx + u(x i ) ( xi x i xi+ dφ i+ + u(x i+ ) x i dx λ i dx (since φ j and λ i overlap only when j {i, i, i + } ) Therefore, g i = u(x i+) u(x i ) x i+ x i ) dφ xi+ i dx λ dφ i i dx + x i dx λ i dx Now we look at the end-points We note that g 0 = Computing as before we get dx λ 0 dx φ 0 λ 0 dx, and g n = dx λ n dx φ n λ n dx g 0 = u(x ) u( ) x, g n = u(x n) u(x n ) x n x n We have the following super-convergence in L -norm This is proved as in [7, 9] Theorem Let h i = x i+ x i for 0 i n, h = ma i n h i, and 0 = (x, x n ) If the point distribution satisfies h i+ h i = O(h ) for 0 i n Then we have the estimate du dx Q h dx h u W 3, ( 0), u W 3 () L ( 0) For the tensor product meshes in two or three dimensions satisfying the above mesh condition this theorem has an easy extension 5
6 Application to quadratic functions Corollary 3 Let u P (R) Then g i reproduces u ( x i ) exactly for all x i, where: +x, i = 0 x x i = i +x i+ x, i n, and g i = 0 dx λ i dx x n +x n x, i = n 0 φ i λ i dx Figure : A non-uniform grid with 8 nodes (vertical lines) The points x i (dots) are also shown Proof: We use the result of the previous theorem to get ( ax i + bx i ax i+ bx i+) g i = (x i x i+ ) = a ( x i ) x i+ + b (xi x i+ ) x i x i+ = a (x i x i+ ) (x i + x i+ ) + b (x i x i+ ) x i x i+ = a (x i + x i+ ) + b On the other hand, ( ) u xi + x i+ ( x i ) = a + b = a (x i + x i+ ) + b = g i 6
7 So, g i reproduces u ( x i ) exactly for i B n Now for i = 0 and i = n, we have ( ) ax g 0 = + bx ax 0 b (x = a ( ) x x 0 + b (x ) ) x = a (x ) (x + ) + b (x ) = a (x + ) + b, x x ( 0 ) ax and g n = n + bx n ax n bx n (x = a ( ) x n x n + b (xn x n ) n x n ) x n x n = a (x n x n ) (x n + x n ) + b (x n x n ) = a (x n + x n ) + b x n x n Since u ( ) = a ( + x ) + b, and u ( x n ) = a (x n + x n ) + b, we have g 0 and g n reproduce u ( ) and u ( x n ), respectively, exactly Remark 4 (Uniform Grid) Let = {α = < x < < x n = β} be a uniform grid on the interval so that x i x i = h, i A n \ {0}, where h is some constant, called the stepsize We note that if our grid is uniform, then x i = x i x i int ( ) So, our gradient recovery operator will reproduce the exact gradient of any quadratic function on the interior of a uniform grid We cannot recover the gradients at the endpoints exactly, however, since and x n x n Corollary 5 Let u P (R) with u(x) = ax +bx+c, and let the grid be uniform with stepsize h and = 0 Then g i u (x i ) = ah for i = 0, n (ie for the endpoints of the grid) Proof: We will start with the case where i = 0 (ie the left endpoint) We know from Theorem 3 that g 0 = a (x ) + b Since our grid is uniform with stepsize h, this simplifies to g 0 = ah + b u ( ) = a + b = b, since = 0 Therefore g 0 u ( ) = (ah + b) b = ah The case for i = n (ie the right endpoint) is proven similarly For a non-uniform grid, we cannot simplify our approximations using the stepsize h, since the spacing between each adjacent node is not always equal We did not make any assumption about the uniformity of the grid in Theorem 3 Thus g i u (x i ), i A n, is not zero for a non-uniform grid This is estimated in the following corollary 7
8 Corollary 6 Let u (x) = ax + bx + c Then, g i u (x i ) = a (x i + x i+ x i ) x i int ( ), g 0 u ( ) = a (x ), and g n u (x n ) = a (x n x n ) Remark 7 For i B n let h i = x i+ x i Then we have g i u (x i ) = a (x i + x i+ x i ) = a(h i+ h i ) We still get superapproximation of the gradient recovery when h i+ h i = O(h ) when i B n Application to cubic functions Corollary 8 Let u P 3 (R) with u(x) = ax 3 + bx + cx + d Then, g i u ( x i ) = a 4 (x i x i+ ) for all i B n, and g 0 u ( ) = a 4 ( x ) g n u ( x n ) = a 4 (x n x n ) where x i is defined as in Theorem 3 Similarly, for all i B n we have g i u (x i ) = a ( x i + x i x i+ + x i+ 3x ) i + b (xi + x i+ x i ), and g 0 u ( ) = a ( x + x x 0) + b (x ), g n u (x n ) = a ( x n + x n x n x n) + b (xn x n ) Proof: The proof of this theorem is similar to Theorem 3 3 Conclusion We have presented an analysis of approximation property of the reconstructed gradient using an oblique projection The reconstruction of the gradient is numerically efficient due to the use of a biorthogonal system It is useful to investigate the extension to higher order finite elements 8
9 References [] M Ainsworth and J T Oden A Posteriori Error Estimation in Finite Element Analysis Wiley Interscience, New York, 000 [] J Chen and D Wang Three-dimensional finite element superconvergent gradient recovery on par6 patterns Numerical Marthematics: Theory, Methods and Applications, 3:78 94, 00 [3] L Chen Superconvergence of tetrahedral linear finite elements International Journal of Numerical Analysis and Modeling,, 3:73 8, 006 [4] J Goodsell Pointwise superconvergence of the gradient for the linear tetrahedral element Numerical Methods for Partial Differential Equations, 0:65 666, 994 [5] BP Lamichhane A gradient recovery operator based on an oblique projection Electronic Transactions on Numerical Analysis, 37:66 7, 00 [6] BP Lamichhane A stabilized mixed finite element method for the biharmonic equation based on biorthogonal systems Journal of Computational and Applied Mathematics, 35: , 0 [7] B Li and Z Zhang Analysis of a class of superconvergence patch recovery techniques for linear and bilinear finite elements Numerical Methods for Partial Differential Equations, 5:5 67, 999 [8] BI Wohlmuth Discretization Methods and Iterative Solvers Based on Domain Decomposition, volume 7 of LNCS Springer Heidelberg, 00 [9] J Xu and Z Zhang Analysis of recovery type a posteriori error estimators for mildly structured grids Mathematics of Computation, 73:39 5, 004 [0] Z Zhang Ultraconvergence of the patch recovery technique Mathematics of Computation, 65:43 437, 996 [] Z Zhang Ultraconvergence of the patch recovery technique ii Mathematics of Computation, 69:4 58, 000 [] OC Zienkiewicz and JZ Zhu The superconvergent patch recovery and a posteriori error estimates part : The recovery technique International Journal for Numerical Methods in Engineering, 33:33 364, 99 [3] OC Zienkiewicz and JZ Zhu The superconvergent patch recovery and a posteriori error estimates part : Error estimates and adaptivity International Journal for Numerical Methods in Engineering, 33:365 38, 99 9
ETNA Kent State University
Electronic Transactions on Numerical Analysis. Volume 37, pp. 166-172, 2010. Copyright 2010,. ISSN 1068-9613. ETNA A GRADIENT RECOVERY OPERATOR BASED ON AN OBLIQUE PROJECTION BISHNU P. LAMICHHANE Abstract.
More informationA 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 Gradient Recovery Operator Based on an Oblique Projection
A Gradient Recovery Operator Based on an Bishnu P. Lamichhane, bishnu.lamichhane@newcastle.edu.au School of Mathematical and Physical Sciences, University of Newcastle, Australia CARMA Retreat July 18,
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 informationPolynomial Preserving Recovery for Quadratic Elements on Anisotropic Meshes
Polynomial Preserving Recovery for Quadratic Elements on Anisotropic Meshes Can Huang, 1 Zhimin Zhang 1, 1 Department of Mathematics, Wayne State University, Detroit, Michigan 480 College of Mathematics
More informationTHE PATCH RECOVERY FOR FINITE ELEMENT APPROXIMATION OF ELASTICITY PROBLEMS UNDER QUADRILATERAL MESHES. Zhong-Ci Shi and Xuejun Xu.
DISCRETE AND CONTINUOUS Website: http://aimsciences.org DYNAMICAL SYSTEMS SERIES B Volume 9, Number, January 2008 pp. 63 82 THE PATCH RECOVERY FOR FINITE ELEMENT APPROXIMATION OF ELASTICITY PROBLEMS UNDER
More informationStudies on Barlow points, Gauss points and Superconvergent points in 1D with Lagrangian and Hermitian finite element basis
Volume 7, N., pp. 75 0, 008 Copyright 008 SBMAC ISSN 0101-805 www.scielo.br/cam Studies on Barlow points, Gauss points and Superconvergent points in 1D with Lagrangian and Hermitian finite element basis
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 informationInterpolation. 1. Judd, K. Numerical Methods in Economics, Cambridge: MIT Press. Chapter
Key References: Interpolation 1. Judd, K. Numerical Methods in Economics, Cambridge: MIT Press. Chapter 6. 2. Press, W. et. al. Numerical Recipes in C, Cambridge: Cambridge University Press. Chapter 3
More informationChapter 4: Interpolation and Approximation. October 28, 2005
Chapter 4: Interpolation and Approximation October 28, 2005 Outline 1 2.4 Linear Interpolation 2 4.1 Lagrange Interpolation 3 4.2 Newton Interpolation and Divided Differences 4 4.3 Interpolation Error
More informationNodal O(h 4 )-superconvergence of piecewise trilinear FE approximations
Preprint, Institute of Mathematics, AS CR, Prague. 2007-12-12 INSTITTE of MATHEMATICS Academy of Sciences Czech Republic Nodal O(h 4 )-superconvergence of piecewise trilinear FE approximations Antti Hannukainen
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 informationon! 0, 1 and 2 In the Zienkiewicz-Zhu SPR p 1 and p 2 are obtained by solving the locally discrete least-squares p
Analysis of a Class of Superconvergence Patch Recovery Techniques for Linear and Bilinear Finite Elements Bo Li Zhimin Zhang y Abstract Mathematical proofs are presented for the derivative superconvergence
More informationA C 0 Linear Finite Element Method for Biharmonic Problems
DOI 10.1007/s10915-017-0501-0 A C 0 Linear Finite Element Method for Biharmonic Problems Hailong Guo 1,2 Zhimin Zhang 1,3 Qingsong Zou 4 Received: 8 December 2016 / Revised: 29 June 2017 / Accepted: 11
More informationRecovery Techniques in Finite Element Methods
Chapter 8 Recovery Techniques in Finite Element Methods Zhimin Zhang Department of Mathematics Wayne State University Detroit, MI 48202 USA zhang@math.wayne.edu Contents 8. Introduction and preliminary.............................
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 brief introduction to finite element methods
CHAPTER A brief introduction to finite element methods 1. Two-point boundary value problem and the variational formulation 1.1. The model problem. Consider the two-point boundary value problem: Given a
More informationOutline. 1 Boundary Value Problems. 2 Numerical Methods for BVPs. Boundary Value Problems Numerical Methods for BVPs
Boundary Value Problems Numerical Methods for BVPs Outline Boundary Value Problems 2 Numerical Methods for BVPs Michael T. Heath Scientific Computing 2 / 45 Boundary Value Problems Numerical Methods for
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 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 informationPOLYNOMIAL PRESERVING GRADIENT RECOVERY AND A POSTERIORI ESTIMATE FOR BILINEAR ELEMENT ON IRREGULAR QUADRILATERALS
INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume, Number, Pages 24 c 2004 Institute for Scientific Computing and Information POLYNOMIAL PRESERVING GRADIENT RECOVERY AND A POSTERIORI ESTIMATE
More informationA posteriori error estimation in the FEM
A posteriori error estimation in the FEM Plan 1. Introduction 2. Goal-oriented error estimates 3. Residual error estimates 3.1 Explicit 3.2 Subdomain error estimate 3.3 Self-equilibrated residuals 3.4
More informationChapter Two: Numerical Methods for Elliptic PDEs. 1 Finite Difference Methods for Elliptic PDEs
Chapter Two: Numerical Methods for Elliptic PDEs Finite Difference Methods for Elliptic PDEs.. Finite difference scheme. We consider a simple example u := subject to Dirichlet boundary conditions ( ) u
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 informationc 2004 Society for Industrial and Applied Mathematics
SIAM J. NUMER. ANAL. Vol. 42, No. 4, pp. 1780 1800 c 2004 Society for Industrial and Applied Mathematics A POSTERIORI ERROR ESTIMATES BASED ON THE POLYNOMIAL PRESERVING RECOVERY AHMED NAGA AND ZHIMIN ZHANG
More informationValidation of the a posteriori error estimator based on polynomial preserving recovery for linear elements
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING Int. J. Numer. Meth. Engng 2004; 61:1860 1893 ublished online 18 October 2004 in Wiley InterScience (www.interscience.wiley.com). DOI: 10.1002/nme.1134
More informationNumerical Methods for Partial Differential Equations
Numerical Methods for Partial Differential Equations Eric de Sturler University of Illinois at Urbana-Champaign Read section 8. to see where equations of type (au x ) x = f show up and their (exact) solution
More informationScientific Computing I
Scientific Computing I Module 8: An Introduction to Finite Element Methods Tobias Neckel Winter 2013/2014 Module 8: An Introduction to Finite Element Methods, Winter 2013/2014 1 Part I: Introduction to
More informationLecture 10 Polynomial interpolation
Lecture 10 Polynomial interpolation Weinan E 1,2 and Tiejun Li 2 1 Department of Mathematics, Princeton University, weinan@princeton.edu 2 School of Mathematical Sciences, Peking University, tieli@pku.edu.cn
More informationLehrstuhl Informatik V. Lehrstuhl Informatik V. 1. solve weak form of PDE to reduce regularity properties. Lehrstuhl Informatik V
Part I: Introduction to Finite Element Methods Scientific Computing I Module 8: An Introduction to Finite Element Methods Tobias Necel Winter 4/5 The Model Problem FEM Main Ingredients Wea Forms and Wea
More informationarxiv: v3 [math.na] 27 Aug 2018
THE HESSIAN DISCRETISATION METHOD FOR FOURTH ORDER LINEAR ELLIPTIC EQUATIONS arxiv:1803.06985v3 [math.na] 27 Aug 2018 Abstract. In this paper, we propose a unified framework, the Hessian discretisation
More informationInput: A set (x i -yy i ) data. Output: Function value at arbitrary point x. What for x = 1.2?
Applied Numerical Analysis Interpolation Lecturer: Emad Fatemizadeh Interpolation Input: A set (x i -yy i ) data. Output: Function value at arbitrary point x. 0 1 4 1-3 3 9 What for x = 1.? Interpolation
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 informationScientific Computing: Numerical Integration
Scientific Computing: Numerical Integration Aleksandar Donev Courant Institute, NYU 1 donev@courant.nyu.edu 1 Course MATH-GA.2043 or CSCI-GA.2112, Fall 2015 Nov 5th, 2015 A. Donev (Courant Institute) Lecture
More informationSection 8.3 Partial Fraction Decomposition
Section 8.6 Lecture Notes Page 1 of 10 Section 8.3 Partial Fraction Decomposition Partial fraction decomposition involves decomposing a rational function, or reversing the process of combining two or more
More informationUniversität Stuttgart
Universität Stuttgart Second Order Lagrange Multiplier Spaces for Mortar Finite Elements in 3D Bishnu P. Lamichhane, Barbara I. Wohlmuth Berichte aus dem Institut für Angewandte Analysis und Numerische
More informationCubic Splines. Antony Jameson. Department of Aeronautics and Astronautics, Stanford University, Stanford, California, 94305
Cubic Splines Antony Jameson Department of Aeronautics and Astronautics, Stanford University, Stanford, California, 94305 1 References on splines 1. J. H. Ahlberg, E. N. Nilson, J. H. Walsh. Theory of
More information2 Two-Point Boundary Value Problems
2 Two-Point Boundary Value Problems Another fundamental equation, in addition to the heat eq. and the wave eq., is Poisson s equation: n j=1 2 u x 2 j The unknown is the function u = u(x 1, x 2,..., x
More informationThe Finite Element Method for the Wave Equation
The Finite Element Method for the Wave Equation 1 The Wave Equation We consider the scalar wave equation modelling acoustic wave propagation in a bounded domain 3, with boundary Γ : 1 2 u c(x) 2 u 0, in
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 informationASYMPTOTICALLY EXACT A POSTERIORI ESTIMATORS FOR THE POINTWISE GRADIENT ERROR ON EACH ELEMENT IN IRREGULAR MESHES. PART II: THE PIECEWISE LINEAR CASE
MATEMATICS OF COMPUTATION Volume 73, Number 246, Pages 517 523 S 0025-5718(0301570-9 Article electronically published on June 17, 2003 ASYMPTOTICALLY EXACT A POSTERIORI ESTIMATORS FOR TE POINTWISE GRADIENT
More 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 information1 Discretizing BVP with Finite Element Methods.
1 Discretizing BVP with Finite Element Methods In this section, we will discuss a process for solving boundary value problems numerically, the Finite Element Method (FEM) We note that such method is a
More informationPartial Fractions. Calculus 2 Lia Vas
Calculus Lia Vas Partial Fractions rational function is a quotient of two polynomial functions The method of partial fractions is a general method for evaluating integrals of rational function The idea
More informationDS-GA 1002 Lecture notes 0 Fall Linear Algebra. These notes provide a review of basic concepts in linear algebra.
DS-GA 1002 Lecture notes 0 Fall 2016 Linear Algebra These notes provide a review of basic concepts in linear algebra. 1 Vector spaces You are no doubt familiar with vectors in R 2 or R 3, i.e. [ ] 1.1
More informationFixed point iteration and root finding
Fixed point iteration and root finding The sign function is defined as x > 0 sign(x) = 0 x = 0 x < 0. It can be evaluated via an iteration which is useful for some problems. One such iteration is given
More informationA POSTERIORI ERROR ESTIMATES BY RECOVERED GRADIENTS IN PARABOLIC FINITE ELEMENT EQUATIONS
A POSTERIORI ERROR ESTIMATES BY RECOVERED GRADIENTS IN PARABOLIC FINITE ELEMENT EQUATIONS D. LEYKEKHMAN AND L. B. WAHLBIN Abstract. This paper considers a posteriori error estimates by averaged gradients
More informationNumerical Analysis of Differential Equations Numerical Solution of Elliptic Boundary Value
Numerical Analysis of Differential Equations 188 5 Numerical Solution of Elliptic Boundary Value Problems 5 Numerical Solution of Elliptic Boundary Value Problems TU Bergakademie Freiberg, SS 2012 Numerical
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 informationAND BARBARA I. WOHLMUTH
A QUASI-DUAL LAGRANGE MULTIPLIER SPACE FOR SERENDIPITY MORTAR FINITE ELEMENTS IN 3D BISHNU P. LAMICHHANE AND BARBARA I. WOHLMUTH Abstract. Domain decomposition techniques provide a flexible tool for the
More informationLecture 1 INF-MAT3350/ : Some Tridiagonal Matrix Problems
Lecture 1 INF-MAT3350/4350 2007: Some Tridiagonal Matrix Problems Tom Lyche University of Oslo Norway Lecture 1 INF-MAT3350/4350 2007: Some Tridiagonal Matrix Problems p.1/33 Plan for the day 1. Notation
More informationIterative Methods for Linear Systems
Iterative Methods for Linear Systems 1. Introduction: Direct solvers versus iterative solvers In many applications we have to solve a linear system Ax = b with A R n n and b R n given. If n is large the
More informationA note on accurate and efficient higher order Galerkin time stepping schemes for the nonstationary Stokes equations
A note on accurate and efficient higher order Galerkin time stepping schemes for the nonstationary Stokes equations S. Hussain, F. Schieweck, S. Turek Abstract In this note, we extend our recent work for
More informationFinite Difference Methods for Boundary Value Problems
Finite Difference Methods for Boundary Value Problems October 2, 2013 () Finite Differences October 2, 2013 1 / 52 Goals Learn steps to approximate BVPs using the Finite Difference Method Start with two-point
More informationAN ALTERNATIVE TECHNIQUE FOR TANGENTIAL STRESS CALCULATION IN DISCONTINUOUS BOUNDARY ELEMENTS
th Pan-American Congress of Applied Mechanics January 04-08, 00, Foz do Iguaçu, PR, Brazil AN ALTERNATIVE TECHNIQUE FOR TANGENTIAL STRESS CALCULATION IN DISCONTINUOUS BOUNDARY ELEMENTS Otávio Augusto Alves
More informationC e n t r u m v o o r W i s k u n d e e n I n f o r m a t i c a
C e n t r u m v o o r W i s k u n d e e n I n f o r m a t i c a Modelling, Analysis and Simulation Modelling, Analysis and Simulation Bilinear forms for the recovery-based discontinuous Galerkin method
More informationFinite Element Method for Ordinary Differential Equations
52 Chapter 4 Finite Element Method for Ordinary Differential Equations In this chapter we consider some simple examples of the finite element method for the approximate solution of ordinary differential
More informationBoundary Value Problems - Solving 3-D Finite-Difference problems Jacob White
Introduction to Simulation - Lecture 2 Boundary Value Problems - Solving 3-D Finite-Difference problems Jacob White Thanks to Deepak Ramaswamy, Michal Rewienski, and Karen Veroy Outline Reminder about
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 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 informationPreliminary Examination in Numerical Analysis
Department of Applied Mathematics Preliminary Examination in Numerical Analysis August 7, 06, 0 am pm. Submit solutions to four (and no more) of the following six problems. Show all your work, and justify
More informationCubic Splines MATH 375. J. Robert Buchanan. Fall Department of Mathematics. J. Robert Buchanan Cubic Splines
Cubic Splines MATH 375 J. Robert Buchanan Department of Mathematics Fall 2006 Introduction Given data {(x 0, f(x 0 )), (x 1, f(x 1 )),...,(x n, f(x n ))} which we wish to interpolate using a polynomial...
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 informationFinite Difference and Finite Element Methods
Finite Difference and Finite Element Methods Georgy Gimel farb COMPSCI 369 Computational Science 1 / 39 1 Finite Differences Difference Equations 3 Finite Difference Methods: Euler FDMs 4 Finite Element
More informationSupraconvergence of a Non-Uniform Discretisation for an Elliptic Third-Kind Boundary-Value Problem with Mixed Derivatives
Supraconvergence of a Non-Uniform Discretisation for an Elliptic Third-Kind Boundary-Value Problem with Mixed Derivatives Etienne Emmrich Technische Universität Berlin, Institut für Mathematik, Straße
More 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 informationA NEW COLLOCATION METHOD FOR SOLVING CERTAIN HADAMARD FINITE-PART INTEGRAL EQUATION
INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume 6, Number, Pages 4 54 c 9 Institute for Scientific Computing and Information A NEW COLLOCATION METHOD FOR SOLVING CERTAIN HADAMARD FINITE-PART
More informationBlock-Structured Adaptive Mesh Refinement
Block-Structured Adaptive Mesh Refinement Lecture 2 Incompressible Navier-Stokes Equations Fractional Step Scheme 1-D AMR for classical PDE s hyperbolic elliptic parabolic Accuracy considerations Bell
More informationPIECEWISE LINEAR FINITE ELEMENT METHODS ARE NOT LOCALIZED
PIECEWISE LINEAR FINITE ELEMENT METHODS ARE NOT LOCALIZED ALAN DEMLOW Abstract. Recent results of Schatz show that standard Galerkin finite element methods employing piecewise polynomial elements of degree
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 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 informationFinite Element Methods
Solving Operator Equations Via Minimization We start with several definitions. Definition. Let V be an inner product space. A linear operator L: D V V is said to be positive definite if v, Lv > for every
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 informationQuadratic forms. Here. Thus symmetric matrices are diagonalizable, and the diagonalization can be performed by means of an orthogonal matrix.
Quadratic forms 1. Symmetric matrices An n n matrix (a ij ) n ij=1 with entries on R is called symmetric if A T, that is, if a ij = a ji for all 1 i, j n. We denote by S n (R) the set of all n n symmetric
More informationSUPERCONVERGENCE ANALYSIS OF THE LINEAR FINITE ELEMENT METHOD AND A GRADIENT RECOVERY POSTPROCESSING ON ANISOTROPIC MESHES
MATHEMATICS OF COMPUTATION Volume 84, Number 29, January 205, Pages 89 7 S 0025-578(204)02846-9 Article electronically published on May 28, 204 SUPERCONVERGENCE ANALYSIS OF THE LINEAR FINITE ELEMENT METHOD
More informationBoundary Value Problems and Iterative Methods for Linear Systems
Boundary Value Problems and Iterative Methods for Linear Systems 1. Equilibrium Problems 1.1. Abstract setting We want to find a displacement u V. Here V is a complete vector space with a norm v V. In
More informationWe consider the problem of finding a polynomial that interpolates a given set of values:
Chapter 5 Interpolation 5. Polynomial Interpolation We consider the problem of finding a polynomial that interpolates a given set of values: x x 0 x... x n y y 0 y... y n where the x i are all distinct.
More informationA C 0 linear finite element method for two fourth-order eigenvalue problems
IMA Journal of Numerical Analysis (2017) 37, 2120 2138 doi: 10.1093/imanum/drw051 Advance Access publication on November 2, 2016 A C 0 linear finite element method for two fourth-order eigenvalue problems
More informationAMS 529: Finite Element Methods: Fundamentals, Applications, and New Trends
AMS 529: Finite Element Methods: Fundamentals, Applications, and New Trends Lecture 3: Finite Elements in 2-D Xiangmin Jiao SUNY Stony Brook Xiangmin Jiao Finite Element Methods 1 / 18 Outline 1 Boundary
More informationThe Plane Stress Problem
The Plane Stress Problem Martin Kronbichler Applied Scientific Computing (Tillämpad beräkningsvetenskap) February 2, 2010 Martin Kronbichler (TDB) The Plane Stress Problem February 2, 2010 1 / 24 Outline
More informationTexas at Austin 2 Nazarbayev University. January 15, 2019
Recovery of the Interface Velocity for the Incompressible Flow in Enhanced Velocity Mixed Finite Element Method arxiv:1901.04401v1 [math.na] 14 Jan 2019 Yerlan Amanbek 1,2, Gurpreet Singh 1, and Mary F.
More informationPhD dissertation defense
Isogeometric mortar methods with applications in contact mechanics PhD dissertation defense Brivadis Ericka Supervisor: Annalisa Buffa Doctor of Philosophy in Computational Mechanics and Advanced Materials,
More informationIn practice one often meets a situation where the function of interest, f(x), is only represented by a discrete set of tabulated points,
1 Interpolation 11 Introduction In practice one often meets a situation where the function of interest, f(x), is only represented by a discrete set of tabulated points, {x i, y i = f(x i ) i = 1 n, obtained,
More informationThe Closed Form Reproducing Polynomial Particle Shape Functions for Meshfree Particle Methods
The Closed Form Reproducing Polynomial Particle Shape Functions for Meshfree Particle Methods by Hae-Soo Oh Department of Mathematics, University of North Carolina at Charlotte, Charlotte, NC 28223 June
More informationSecond Order Elliptic PDE
Second Order Elliptic PDE T. Muthukumar tmk@iitk.ac.in December 16, 2014 Contents 1 A Quick Introduction to PDE 1 2 Classification of Second Order PDE 3 3 Linear Second Order Elliptic Operators 4 4 Periodic
More informationA truly meshless Galerkin method based on a moving least squares quadrature
A truly meshless Galerkin method based on a moving least squares quadrature Marc Duflot, Hung Nguyen-Dang Abstract A new body integration technique is presented and applied to the evaluation of the stiffness
More informationExact a posteriori error analysis of the least squares nite element method 1
Applied Mathematics and Computation 116 (2000) 297±305 www.elsevier.com/locate/amc Exact a posteriori error analysis of the least squares nite element method 1 Jinn-Liang Liu * Department of Applied Mathematics,
More informationHamburger Beiträge zur Angewandten Mathematik
Hamburger Beiträge zur Angewandten Mathematik Numerical analysis of a control and state constrained elliptic control problem with piecewise constant control approximations Klaus Deckelnick and Michael
More informationProjected Surface Finite Elements for Elliptic Equations
Available at http://pvamu.edu/aam Appl. Appl. Math. IN: 1932-9466 Vol. 8, Issue 1 (June 2013), pp. 16 33 Applications and Applied Mathematics: An International Journal (AAM) Projected urface Finite Elements
More informationAlgorithms for Scientific Computing
Algorithms for Scientific Computing Hierarchical Methods and Sparse Grids d-dimensional Hierarchical Basis Michael Bader Technical University of Munich Summer 208 Intermezzo/ Big Picture : Archimedes Quadrature
More informationA Stabilized Finite Element Method for the Darcy Problem on Surfaces
arxiv:1511.03747v1 [math.na] 12 Nov 2015 A Stabilized Finite Element Method for the Darcy Problem on Surfaces Peter Hansbo Mats G. Larson October 8, 2018 Abstract We consider a stabilized finite element
More informationA posteriori error estimates for Maxwell Equations
www.oeaw.ac.at A posteriori error estimates for Maxwell Equations J. Schöberl RICAM-Report 2005-10 www.ricam.oeaw.ac.at A POSTERIORI ERROR ESTIMATES FOR MAXWELL EQUATIONS JOACHIM SCHÖBERL Abstract. Maxwell
More informationSuperconvergent Two-Grid Methods for Elliptic Eigenvalue Problems
DOI 1.17/s1915-16-245-2 Superconvergent Two-Grid Methods for Elliptic Eigenvalue Problems Hailong Guo 2,3 Zhimin Zhang 1,2 Ren Zhao 2 Received: 4 April 215 / Revised: 12 May 216 / Accepted: 28 June 216
More informationLecture 9 Approximations of Laplace s Equation, Finite Element Method. Mathématiques appliquées (MATH0504-1) B. Dewals, C.
Lecture 9 Approximations of Laplace s Equation, Finite Element Method Mathématiques appliquées (MATH54-1) B. Dewals, C. Geuzaine V1.2 23/11/218 1 Learning objectives of this lecture Apply the finite difference
More informationNon-Conforming Finite Element Methods for Nonmatching Grids in Three Dimensions
Non-Conforming Finite Element Methods for Nonmatching Grids in Three Dimensions Wayne McGee and Padmanabhan Seshaiyer Texas Tech University, Mathematics and Statistics (padhu@math.ttu.edu) Summary. In
More informationLECTURE 16 GAUSS QUADRATURE In general for Newton-Cotes (equispaced interpolation points/ data points/ integration points/ nodes).
CE 025 - Lecture 6 LECTURE 6 GAUSS QUADRATURE In general for ewton-cotes (equispaced interpolation points/ data points/ integration points/ nodes). x E x S fx dx hw' o f o + w' f + + w' f + E 84 f 0 f
More informationTaylor Series and Numerical Approximations
Taylor Series and Numerical Approximations Hilary Weller h.weller@reading.ac.uk August 7, 05 An introduction to the concept of a Taylor series and how these are used in numerical analysis to find numerical
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 informationAn Efficient Algorithm Based on Quadratic Spline Collocation and Finite Difference Methods for Parabolic Partial Differential Equations.
An Efficient Algorithm Based on Quadratic Spline Collocation and Finite Difference Methods for Parabolic Partial Differential Equations by Tong Chen A thesis submitted in conformity with the requirements
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 information