DISCRETE MAXIMUM PRINCIPLES in THE FINITE ELEMENT SIMULATIONS
|
|
- Sharon Singleton
- 5 years ago
- Views:
Transcription
1 DISCRETE MAXIMUM PRINCIPLES in THE FINITE ELEMENT SIMULATIONS Sergey Korotov BCAM Basque Center for Applied Mathematics 1
2 The presentation is based on my collaboration with several colleagues - I will mention their names (and papers) in the appropriate places of the talk. 2
3 Model Elliptic Problem Find a function u such that div(a u)+cu = f in Ω, (1) u = g on Ω. (2) Ω R d is a bounded polytope with Lipschitz boundary Ω. The diffusive tensor A is assumed to be a symmetric and uniformly positive definite matrix. The reactive coefficient c is assumed to be nonnegative in Ω. 3
4 Continuous Maximum Principle The classical (smooth) solutions of many elliptic/parabolic problems are known to satisfy various (continuous) maximum principles (or CMPs in short). For our model elliptic problem standard CMP reads as follows: f 0 = max x Ω u(x) max{0, maxg(s)}. (3) s Ω CMPs are not only pure mathematical properties of PDE models, they also reflect well physical behaviour of phenomena modelled. 4
5 Preliminaries Some reasonable discrete analogue of CMP (which may depend, in general, on the nature of numerical technique used) is often called the discrete maximum principle (or DMP in short). Conditions on parameters of computational schemes (e.g. on mesh shape and size, time-step values, etc) a priori providing a validity of DMPs by computed approximations are highly desired In what follows we shall only consider elliptic problems 5
6 The first DMP and certain conditions providing its validity were (probably?) formulated by R. Varga for the finite difference method (FDM) in On a discrete maximum principle, J. SIAM Numer. Anal. 3 (1966), However, only a special case f = 0 was considered there and CMP was of the following form: max x Ω u(x) max s Ω g(s). 6
7 Later, in 1970 Ph. Ciarlet presented a more sophisticated form of DMP suitable for both, finite element (FE) and finite difference (FD) types of discretization. He also proposed a practically convenient set of (sufficient) conditions on matrix blocks involved, which provides a validity of his DMP. Since that time Ciarlet conditions became popular in numerical community for proving DMPs for various problems of elliptic type. We shall present the conditions proposed by Ciarlet, sketch appearing practical problems with them when FEM is used as a main numerical technique and discuss some generalizations. 7
8 Weak Formulation Standard (linear) FE schemes are based on weak formulations. For our test problem it reads as follows: Find u g +H 1 0(Ω) such that a(u,v) = F(v) v H 1 0(Ω), where a(u,v) = Ω A u vdx+ Ω cuvdx and F(v) = Ω fvdx. Here, the matrix A is assumed to be in [L (Ω)] d d, c L (Ω), g H 1 (Ω), and f L 2 (Ω). The existence and uniqueness of the weak solution u is provided by the standard Lax-Milgram lemma. 8
9 FE Discretization: Mesh & Basis Functions Let T h be a FE mesh of Ω with interior nodes B 1,...,B N lying in Ω and boundary nodes B N+1,...,B N+N lying on Ω. Let V h be a finite-dimensional subspace of H 1 (Ω), associated with T h and its nodes, being spanned by basis functions φ 1,φ 2,...,φ N+N s.t. φ i 0 in Ω, i = 1,...,N +N, and N+N i=1 φ i 1 in Ω. We also assume that φ 1,φ 2,...,φ N vanish on the boundary Ω. Thus, they span a finite-dimensional subspace V 0 h of H 1 0(Ω). 9
10 FE Discretization: Matrix Equation Let g h = N+N i=n+1 g iφ i V h be a suitable approximation of the function g, for example its nodal interpolant. FE approximation is defined as a function u h g h +V 0 h such that a(u h,v h ) = F(v h ) v h V 0 h. In fact, u h = N+N i=1 y i φ i, where ȳ = [y 1,...,y N+N ] is the solution of of the following square system of equations Āȳ = F, (4) with Ā = A A 0 I, F = F F. (5) 10
11 Ā = A A 0 I, F = F, and Ā 1 = F A 1 A 1 A 0 I. Blocks A and A are matrices of size N N and N N, respectively, I stands for the unit matrix, and 0 for the zero matrix. Entries of Ā are a ij = a(φ j,φ i ), i = 1,...,N, j = 1,...,N +N. The block-vector F consists of entries f i = F(φ i ), i = 1,...,N, and the block-vector F has entries f i = g i, i = N +1,...,N +N. For the later reference we include the formula for Ā 1. Ā is nonsingular if and only if A is nonsingular. 11
12 Sufficient Algebraic Conditions of Ph. Ciarlet We should distinguish two basically different types of DMPs. Algebraic DMP: A natural algebraic analogue of CMP is F 0 = max y i max{0, max y j }. i=1,...,n+n j=n+1,...,n+n Functional DMP: A natural functional imitation of CMP is f 0 = max Ω u h max{0,max Ω u h}. It is easy to see that the above types of DMPs are actually equivalent e.g. for linear, multilinear, or prismatic FEs. (However, these DMPs are not equivalent, in general, for higher-order FEs.) 12
13 Ā = A A 0 I, F = F F, and Ā 1 = A 1 A 1 A 0 I. In the paper Discrete maximum principle for finite-difference operators, Aequationes Math. 4 (1970), , Ciarlet proved Theorem 1. The algebraic DMP is satisfied if and only if (A) Ā is monotone (i.e., Ā nonsingular and Ā 1 0) (B) ξ +A 1 A ξ 0, where ξ and ξ are vectors of all ones of sizes N and N, respectively. Since conditions (A) and (B) are difficult to verify in practice, Ciarlet proposed in the same work a more convenient set of sufficient conditions. 13
14 Ā = A A 0 I, F = F F, and Ā 1 = A 1 A 1 A 0 I. Theorem 2. The algebraic DMP is valid provided the matrix Ā satisfies (a) a ii > 0, i = 1,...,N, (b) a ij 0, i j, i = 1,...,N, j = 1,...,N +N, a ij 0, i = 1,...,N, (c) N+N j=1 (d) A is irreducibly diagonally dominant. Ciarlet essentially proposed the above conditions in order to utilize the following result of R. Varga from his book Matrix Iterative Analysis, 1962: 14
15 Lemma 3. If A R N N is an irreducibly diagonally dominant matrix with strictly positive diagonal and nonpositive off-diagonal entries then A 1 > 0. Ā = A A 0 I, F = F, and Ā 1 = F A 1 A 1 A 0 I. Remark 1. However, in the proof we only need monotonicity of A, i.e. A 1 0, not more. Remark 2. In case g = 0, the system reduces to Ay = F, and CMP is f 0 = max x Ω To immitate this, we also need A 1 0, not more. u(x) 0. (6) 15
16 Analysis of Ciarlet Conditions Condition (a) is always satisfied for elliptic problems. The row sums in (c) are nonnegative automatically for our test problem, because the basis functions are nonnegative and form a partition of unity: N+N j=1 a ij = a ( N+N j=1 φ j,φ i ) = a(1,φ i ) = Ω cφ i 0, i = 1,...,N. (7) On the other hand, the irreducibility of A required in (d) is not always obvious. 16
17 Problems with Irreducibility Figure 1: [Hannukainen, Korotov, Vejchodský, 2009]. Meshes leading to reducible matrix A for the Poisson problem with Dirichlet boundary conditions. Angles with dots are right. 17
18 Associated Geometric Conditions For some types of FEs, Ā can be computed explicitly, therefore condition (b): a ij 0, i j, can often be guaranteed a priori by imposing suitable geometrical requirements on the shape and size of FE meshes. For example, let A = I, then there are some known geometrical conditions providing (b). 18
19 Simplicial FE Meshes For simplicity we assume that c is constant and d = 2 If nodes B i and B j are not connected by an edge then a ij = 0 B i T 1 α 2 T 2 α 1 B j Otherwise, a ij always consists of two contributions computed over two (adjacent) triangles from suppφ i suppφ j : a ij = 1 ) (ctgα 1 +ctgα 2 + c ) (meast 1 +meast [Fujii, 1973], [Ruas Santos, 1982], [Ciarlet-book, 1978]. 19
20 a ij = 1 2 (ctgα 1 +ctgα 2 ) + c ) (meast 1 +meast 2 12 In case c = 0, we get a ij 0 if ctgα 1 +ctgα 2 0 It holds if e.g. α k π 2. This observation leads to the concept of nonobtuse triangulations In case c > 0, we get a ij 0 if c ) (meast 1 +meast 2 1 ) (ctgα 1 +ctgα It is possible if e.g. α k π 2 ε (ε > 0) and triangulations are sufficiently fine. This leads to the concept of acute triangulations 20
21 Similarly for d 3: B i B j K r B i S ir B j S jr α K r ij φ j φ i dx = meas d 1S ir meas d 1S jr K r d 2 meas d K r d! φ j φ i dx = K r (d+2)! meas dk r cosα K r ij [Křížek, Lin Qun, 1995], [Xu, Zikatanov, 1999], [Brandts, Korotov, Křížek, 2007] 21
22 Block FE Meshes [Karátson, Korotov, Křížek, 2007], [Korotov, Vejchodský, 2010]: the results strongly depend on dimension Let K be a block of a d-dimensional block mesh with edges of lengths b 1,b 2,...,b d, and let B i and B j be its vertices connected by b 1, then a K ij = b 1b 2...b d 3 d 1 ( d k=2 1 2b 2 k 1 b 2 1 ), i j d = 2: all rectangular elements have to be nonnarrow for (b), i.e. 2/2 b1 /b 2 2 d = 3: trilinear elements have to be cubes d 4: condition (b) is not valid even on hyper-cubes 22
23 Prismatic FE Meshes [Hannukainen, Korotov, Vejchodský, 2009]: for 3D meshes consiting of right triangular prisms the altitudes of prisms should be limited from both sides by certain quantities dependent on the area and angles of the triangular base (and the magnitude of coefficient c) to provide (b). F D E C d b γ a α β A c B 23
24 No Irreducibility-Proof is Needed A real square matrix A is an M-matrix if all its off-diagonal entries are nonpositive and if it is nonsingular and A 1 0. A real square matrix A is a Stieltjes matrix if all its off-diagonal entries are nonpositive and if it is symmetric and positive definite. The following lemma from [Varga-book, 1962] enables to eliminate (a), (c), and (d) from the set of Ciarlet conditions in our case: Lemma 4. If A is a Stieltjes matrix then it is also an M-matrix. Theorem 5. If the finite element matrix Ā, associated to our model problem, satisfies condition (b), a ij 0, i j, then the algebraic DMP is valid. For the proof see [Hannukainen, Korotov, Vejchodský, 2009]. 24
25 Geometrical Conditions are Essential [Brandts, Korotov, Křížek, Šolc, 2009]: one single obtuse triangle in the FE mesh can destroy DMP (e.g. for the Poisson equation with zero Dirichlet b.c.) 25
26 Problems with Mesh Reconstructions To improve the accuracy of FE computations we need to refine meshes globally and locally. How to preserve the desired geometric properties then? red refinement For global mesh refinements 2d red technique always preserves properties of acuteness and nonobtuseness 26
27 For local mesh refinements green post-refining (to preserve conformity) is often needed, but it generally leads to obtuse angles. Only in special cases we can preserve nonobtuseness for local mesh refinements using both, red and green, refinements: red refinement green refinement red + green producing nonobtuse local refinements 27
28 Problems with 3D FE meshes Existence of a face-to-face partition of a cube into acute tetrahedra was surprisingly (!) an open problem till D red (global) refinement technique does not generally preserve nonobtuseness or acuteness. green refinement red green refinement Problems with local refinements (especially for acuteness)... 28
29 More Problems Local refinements of block and prismatic FE meshes... Full diffusive tensor... Higher-dimensional problems (blocks?, acute simplices?)... 29
30 THANK YOU FOR YOUR ATTENTION! 30
31 On Sharpness of Theoretical Conditions Stieltjes matrices form only a certain subclass of monotone matrices. In what follows we test how sharp the conditions based on the Stieltjes matrix concept are, i.e., we try to test how wide is the class of meshes which lead to A being monotone but not Stieltjes. For this purpose we solve the 2D Poisson problem with homogenous Dirichlet boundary conditions on various domains using various triangulations. We present three representative tests... 31
32 In each test we construct a simple triangulation which is characterized by two parameters (angles) 0 < α < π and 0 < γ < π: γ α α α+γ α α α α γ γ γ γ γ Figure 2: The basic meshes for the three tests. The meshes used for the actual computations are 10-fold refined basic meshes. 32
33 Figure 3: Domain 1: triangulations with non-obtuse maximal angle (i.e. matrix A is a Stieltjes matrix). Domain 2: triangulations with obtuse maximal angle providing the DMP (i.e. matrix A is monotone but not Stieltjes). Domain 3: triangulations with obtuse maximal angle, DMP is not valid (matrix A is not monotone). 33
34 Conclusions If we compare the sizes of domains 1 with domains 2 we may conclude that the standard sufficient condition (b) is not very sharp. There is a wide space for its generalization. However, any generalization have to utilize more general criteria for the monotonicity of a matrix. These criteria are often quite complicated and their application for proving DMPs is not so straightforward. 34
On angle conditions in the finite element method. Institute of Mathematics, Academy of Sciences Prague, Czech Republic
On angle conditions in the finite element method Michal Křížek Institute of Mathematics, Academy of Sciences Prague, Czech Republic Joint work with Jan Brandts (University of Amsterdam), Antti Hannukainen
More 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 informationOn the equivalence of regularity criteria for triangular and tetrahedral finite element partitions
Computers and Mathematics with Applications 55 (2008) 2227 2233 www.elsevier.com/locate/camwa On the equivalence of regularity criteria for triangular and tetrahedral finite element partitions Jan Brandts
More informationA very short introduction to the Finite Element Method
A very short introduction to the Finite Element Method Till Mathis Wagner Technical University of Munich JASS 2004, St Petersburg May 4, 2004 1 Introduction This is a short introduction to the finite element
More informationThe maximum angle condition is not necessary for convergence of the finite element method
1 2 3 4 The maximum angle condition is not necessary for convergence of the finite element method Antti Hannukainen 1, Sergey Korotov 2, Michal Křížek 3 October 19, 2010 5 6 7 8 9 10 11 12 13 14 15 16
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 informationPDE & FEM TERMINOLOGY. BASIC PRINCIPLES OF FEM.
PDE & FEM TERMINOLOGY. BASIC PRINCIPLES OF FEM. Sergey Korotov Basque Center for Applied Mathematics / IKERBASQUE http://www.bcamath.org & http://www.ikerbasque.net 1 Introduction The analytical solution
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 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 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 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 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 informationPARTITION OF UNITY FOR THE STOKES PROBLEM ON NONMATCHING GRIDS
PARTITION OF UNITY FOR THE STOES PROBLEM ON NONMATCHING GRIDS CONSTANTIN BACUTA AND JINCHAO XU Abstract. We consider the Stokes Problem on a plane polygonal domain Ω R 2. We propose a finite element method
More informationA 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 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 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 Robust Preconditioner for the Hessian System in Elliptic Optimal Control Problems
A Robust Preconditioner for the Hessian System in Elliptic Optimal Control Problems Etereldes Gonçalves 1, Tarek P. Mathew 1, Markus Sarkis 1,2, and Christian E. Schaerer 1 1 Instituto de Matemática Pura
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 informationFinite difference method for elliptic problems: I
Finite difference method for elliptic problems: I Praveen. C praveen@math.tifrbng.res.in Tata Institute of Fundamental Research Center for Applicable Mathematics Bangalore 560065 http://math.tifrbng.res.in/~praveen
More informationChapter 3 Conforming Finite Element Methods 3.1 Foundations Ritz-Galerkin Method
Chapter 3 Conforming Finite Element Methods 3.1 Foundations 3.1.1 Ritz-Galerkin Method Let V be a Hilbert space, a(, ) : V V lr a bounded, V-elliptic bilinear form and l : V lr a bounded linear functional.
More informationFEniCS Course. Lecture 0: Introduction to FEM. Contributors Anders Logg, Kent-Andre Mardal
FEniCS Course Lecture 0: Introduction to FEM Contributors Anders Logg, Kent-Andre Mardal 1 / 46 What is FEM? The finite element method is a framework and a recipe for discretization of mathematical problems
More information1. Let a(x) > 0, and assume that u and u h are the solutions of the Dirichlet problem:
Mathematics Chalmers & GU TMA37/MMG800: Partial Differential Equations, 011 08 4; kl 8.30-13.30. Telephone: Ida Säfström: 0703-088304 Calculators, formula notes and other subject related material are not
More informationFrom Completing the Squares and Orthogonal Projection to Finite Element Methods
From Completing the Squares and Orthogonal Projection to Finite Element Methods Mo MU Background In scientific computing, it is important to start with an appropriate model in order to design effective
More informationFinite Elements. Colin Cotter. January 15, Colin Cotter FEM
Finite Elements January 15, 2018 Why Can solve PDEs on complicated domains. Have flexibility to increase order of accuracy and match the numerics to the physics. has an elegant mathematical formulation
More informationLocal Mesh Refinement with the PCD Method
Advances in Dynamical Systems and Applications ISSN 0973-5321, Volume 8, Number 1, pp. 125 136 (2013) http://campus.mst.edu/adsa Local Mesh Refinement with the PCD Method Ahmed Tahiri Université Med Premier
More 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 informationHelsinki University of Technology Institute of Mathematics Research Reports
Helsinki University of Technology Institute of Mathematics Research Reports Espoo 2008 A550 DISCRETE MAXIMUM PRINCIPLES FOR FE SOLUTIONS OF NONSTATIONARY DIFFUSION-REACTION PROBLEMS WITH MIXED BOUNDARY
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 informationNumerical methods for PDEs FEM convergence, error estimates, piecewise polynomials
Platzhalter für Bild, Bild auf Titelfolie hinter das Logo einsetzen Numerical methods for PDEs FEM convergence, error estimates, piecewise polynomials Dr. Noemi Friedman Contents of the course Fundamentals
More informationOverlapping Schwarz preconditioners for Fekete spectral elements
Overlapping Schwarz preconditioners for Fekete spectral elements R. Pasquetti 1, L. F. Pavarino 2, F. Rapetti 1, and E. Zampieri 2 1 Laboratoire J.-A. Dieudonné, CNRS & Université de Nice et Sophia-Antipolis,
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 informationFinite Elements. Colin Cotter. February 22, Colin Cotter FEM
Finite Elements February 22, 2019 In the previous sections, we introduced the concept of finite element spaces, which contain certain functions defined on a domain. Finite element spaces are examples of
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 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 informationDISCRETE MAXIMUM PRINCIPLE AND ADEQUATE DISCRETIZATIONS OF LINEAR PARABOLIC PROBLEMS
DISCRETE MAXIMUM PRINCIPLE AND ADEQUATE DISCRETIZATIONS OF LINEAR PARABOLIC PROBLEMS ISTVÁN FARAGÓ AND RÓBERT HORVÁTH Abstract. In this paper, we analyze the connections between the different qualitative
More information[2] (a) Develop and describe the piecewise linear Galerkin finite element approximation of,
269 C, Vese Practice problems [1] Write the differential equation u + u = f(x, y), (x, y) Ω u = 1 (x, y) Ω 1 n + u = x (x, y) Ω 2, Ω = {(x, y) x 2 + y 2 < 1}, Ω 1 = {(x, y) x 2 + y 2 = 1, x 0}, Ω 2 = {(x,
More informationETNA Kent State University
Electronic Transactions on Numerical Analysis. Volume 37, pp. 166-172, 2010. Copyright 2010,. ISSN 1068-9613. ETNA A GRADIENT RECOVERY OPERATOR BASED ON AN OBLIQUE PROJECTION BISHNU P. LAMICHHANE Abstract.
More 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 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 informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Numerical Methods for Partial Differential Equations Finite Difference Methods
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Variational Problems of the Dirichlet BVP of the Poisson Equation 1 For the homogeneous
More informationComparative Analysis of Mesh Generators and MIC(0) Preconditioning of FEM Elasticity Systems
Comparative Analysis of Mesh Generators and MIC(0) Preconditioning of FEM Elasticity Systems Nikola Kosturski and Svetozar Margenov Institute for Parallel Processing, Bulgarian Academy of Sciences Abstract.
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 informationMath 660-Lecture 15: Finite element spaces (I)
Math 660-Lecture 15: Finite element spaces (I) (Chapter 3, 4.2, 4.3) Before we introduce the concrete spaces, let s first of all introduce the following important lemma. Theorem 1. Let V h consists of
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 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 informationOn the design of higher-order FEM satisfying the discrete maximum principle
On the design of higher-order FEM satisfying the discrete maximum principle Dmitri Kuzmin Institute of Applied Mathematics (LS III), University of Dortmund Vogelpothsweg 87, D-44227, Dortmund, Germany
More informationOn Surface Meshes Induced by Level Set Functions
On Surface Meshes Induced by Level Set Functions Maxim A. Olshanskii, Arnold Reusken, and Xianmin Xu Bericht Nr. 347 Oktober 01 Key words: surface finite elements, level set function, surface triangulation,
More informationMultigrid absolute value preconditioning
Multigrid absolute value preconditioning Eugene Vecharynski 1 Andrew Knyazev 2 (speaker) 1 Department of Computer Science and Engineering University of Minnesota 2 Department of Mathematical and Statistical
More informationNUMERICAL SOLUTION OF CONVECTION DIFFUSION EQUATIONS USING UPWINDING TECHNIQUES SATISFYING THE DISCRETE MAXIMUM PRINCIPLE
Proceedings of the Czech Japanese Seminar in Applied Mathematics 2005 Kuju Training Center, Oita, Japan, September 15-18, 2005 pp. 69 76 NUMERICAL SOLUTION OF CONVECTION DIFFUSION EQUATIONS USING UPWINDING
More informationNew class of finite element methods: weak Galerkin methods
New class of finite element methods: weak Galerkin methods Xiu Ye University of Arkansas at Little Rock Second order elliptic equation Consider second order elliptic problem: a u = f, in Ω (1) u = 0, on
More informationConstruction of `Wachspress type' rational basis functions over rectangles
Proc. Indian Acad. Sci. (Math. Sci.), Vol. 110, No. 1, February 2000, pp. 69±77. # Printed in India Construction of `Wachspress type' rational basis functions over rectangles P L POWAR and S S RANA Department
More informationCONVERGENCE THEORY. G. ALLAIRE CMAP, Ecole Polytechnique. 1. Maximum principle. 2. Oscillating test function. 3. Two-scale convergence
1 CONVERGENCE THEOR G. ALLAIRE CMAP, Ecole Polytechnique 1. Maximum principle 2. Oscillating test function 3. Two-scale convergence 4. Application to homogenization 5. General theory H-convergence) 6.
More informationFinite Elements for Nonlinear Problems
Finite Elements for Nonlinear Problems Computer Lab 2 In this computer lab we apply finite element method to nonlinear model problems and study two of the most common techniques for solving the resulting
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 informationSolving PDEs with Multigrid Methods p.1
Solving PDEs with Multigrid Methods Scott MacLachlan maclachl@colorado.edu Department of Applied Mathematics, University of Colorado at Boulder Solving PDEs with Multigrid Methods p.1 Support and Collaboration
More 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 informationThe Dirichlet boundary problems for second order parabolic operators satisfying a Carleson condition
The Dirichlet boundary problems for second order parabolic operators satisfying a Carleson condition Sukjung Hwang CMAC, Yonsei University Collaboration with M. Dindos and M. Mitrea The 1st Meeting of
More informationConvergence analysis of finite element approximations of the Joule heating problem in three spatial dimensions
BIT Numer Math (2010) 50: 781 795 DOI 10.1007/s10543-010-0287-z Convergence analysis of finite element approximations of the Joule heating problem in three spatial dimensions Michael J. Holst Mats G. Larson
More informationNumerical Solution I
Numerical Solution I Stationary Flow R. Kornhuber (FU Berlin) Summerschool Modelling of mass and energy transport in porous media with practical applications October 8-12, 2018 Schedule Classical Solutions
More informationA simple FEM solver and its data parallelism
A simple FEM solver and its data parallelism Gundolf Haase Institute for Mathematics and Scientific Computing University of Graz, Austria Chile, Jan. 2015 Partial differential equation Considered Problem
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 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 informationTime-dependent variational forms
Time-dependent variational forms Hans Petter Langtangen 1,2 1 Center for Biomedical Computing, Simula Research Laboratory 2 Department of Informatics, University of Oslo Oct 30, 2015 PRELIMINARY VERSION
More informationSolving Symmetric Indefinite Systems with Symmetric Positive Definite Preconditioners
Solving Symmetric Indefinite Systems with Symmetric Positive Definite Preconditioners Eugene Vecharynski 1 Andrew Knyazev 2 1 Department of Computer Science and Engineering University of Minnesota 2 Department
More informationWeak Galerkin Finite Element Scheme and Its Applications
Weak Galerkin Finite Element Scheme and Its Applications Ran Zhang Department of Mathematics Jilin University, China IMS, Singapore February 6, 2015 Talk Outline Motivation WG FEMs: Weak Operators + Stabilizer
More informationMultigrid Method ZHONG-CI SHI. Institute of Computational Mathematics Chinese Academy of Sciences, Beijing, China. Joint work: Xuejun Xu
Multigrid Method ZHONG-CI SHI Institute of Computational Mathematics Chinese Academy of Sciences, Beijing, China Joint work: Xuejun Xu Outline Introduction Standard Cascadic Multigrid method(cmg) Economical
More informationIsogeometric Analysis:
Isogeometric Analysis: some approximation estimates for NURBS L. Beirao da Veiga, A. Buffa, Judith Rivas, G. Sangalli Euskadi-Kyushu 2011 Workshop on Applied Mathematics BCAM, March t0th, 2011 Outline
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 informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Nonconformity and the Consistency Error First Strang Lemma Abstract Error Estimate
More informationWELL POSEDNESS OF PROBLEMS I
Finite Element Method 85 WELL POSEDNESS OF PROBLEMS I Consider the following generic problem Lu = f, where L : X Y, u X, f Y and X, Y are two Banach spaces We say that the above problem is well-posed (according
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 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 informationHelsinki University of Technology Institute of Mathematics Research Reports
Helsinki University of Technology Institute of Mathematics Research Reports Espoo 2008 A551 DISCRETE MAXIMUM PRINCIPLES FOR THE FEM SOLUTION OF SOME NONLINEAR PARABOLIC PROBLEMS István Faragó János Karátson
More informationFinite Element Method (FEM)
Finite Element Method (FEM) The finite element method (FEM) is the oldest numerical technique applied to engineering problems. FEM itself is not rigorous, but when combined with integral equation techniques
More informationA Finite Element Method for an Ill-Posed Problem. Martin-Luther-Universitat, Fachbereich Mathematik/Informatik,Postfach 8, D Halle, Abstract
A Finite Element Method for an Ill-Posed Problem W. Lucht Martin-Luther-Universitat, Fachbereich Mathematik/Informatik,Postfach 8, D-699 Halle, Germany Abstract For an ill-posed problem which has its origin
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 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 informationINF-SUP CONDITION FOR OPERATOR EQUATIONS
INF-SUP CONDITION FOR OPERATOR EQUATIONS LONG CHEN We study the well-posedness of the operator equation (1) T u = f. where T is a linear and bounded operator between two linear vector spaces. We give equivalent
More informationResearch Article A Two-Grid Method for Finite Element Solutions of Nonlinear Parabolic Equations
Abstract and Applied Analysis Volume 212, Article ID 391918, 11 pages doi:1.1155/212/391918 Research Article A Two-Grid Method for Finite Element Solutions of Nonlinear Parabolic Equations Chuanjun Chen
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 informationSparse Linear Systems. Iterative Methods for Sparse Linear Systems. Motivation for Studying Sparse Linear Systems. Partial Differential Equations
Sparse Linear Systems Iterative Methods for Sparse Linear Systems Matrix Computations and Applications, Lecture C11 Fredrik Bengzon, Robert Söderlund We consider the problem of solving the linear system
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 informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Discretization of Boundary Conditions Discretization of Boundary Conditions On
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 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 informationLinear Algebra Section 2.6 : LU Decomposition Section 2.7 : Permutations and transposes Wednesday, February 13th Math 301 Week #4
Linear Algebra Section. : LU Decomposition Section. : Permutations and transposes Wednesday, February 1th Math 01 Week # 1 The LU Decomposition We learned last time that we can factor a invertible matrix
More information(f P Ω hf) vdx = 0 for all v V h, if and only if Ω (f P hf) φ i dx = 0, for i = 0, 1,..., n, where {φ i } V h is set of hat functions.
Answer ALL FOUR Questions Question 1 Let I = [, 2] and f(x) = (x 1) 3 for x I. (a) Let V h be the continous space of linear piecewise function P 1 (I) on I. Write a MATLAB code to compute the L 2 -projection
More informationYongdeok Kim and Seki Kim
J. Korean Math. Soc. 39 (00), No. 3, pp. 363 376 STABLE LOW ORDER NONCONFORMING QUADRILATERAL FINITE ELEMENTS FOR THE STOKES PROBLEM Yongdeok Kim and Seki Kim Abstract. Stability result is obtained for
More 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 informationTopological Data Analysis - Spring 2018
Topological Data Analysis - Spring 2018 Simplicial Homology Slightly rearranged, but mostly copy-pasted from Harer s and Edelsbrunner s Computational Topology, Verovsek s class notes. Gunnar Carlsson s
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 informationOptimal Interface Conditions for an Arbitrary Decomposition into Subdomains
Optimal Interface Conditions for an Arbitrary Decomposition into Subdomains Martin J. Gander and Felix Kwok Section de mathématiques, Université de Genève, Geneva CH-1211, Switzerland, Martin.Gander@unige.ch;
More informationFinite Elements. Colin Cotter. January 18, Colin Cotter FEM
Finite Elements January 18, 2019 The finite element Given a triangulation T of a domain Ω, finite element spaces are defined according to 1. the form the functions take (usually polynomial) when restricted
More informationSimplicial finite elements in higher dimensions
Helsinki University of Technology, Institute of Mathematics, Research Reports Teknillisen korkeakoulun matematiikan laitoksen tutkimusraporttisarja Espoo 2007 A517 Simplicial finite elements in higher
More informationAn additive average Schwarz method for the plate bending problem
J. Numer. Math., Vol. 10, No. 2, pp. 109 125 (2002) c VSP 2002 Prepared using jnm.sty [Version: 02.02.2002 v1.2] An additive average Schwarz method for the plate bending problem X. Feng and T. Rahman Abstract
More informationFriedrich symmetric systems
viii CHAPTER 8 Friedrich symmetric systems In this chapter, we describe a theory due to Friedrich [13] for positive symmetric systems, which gives the existence and uniqueness of weak solutions of boundary
More informationAn inverse problem for a system of steady-state reaction-diffusion equations on a porous domain using a collage-based approach
Journal of Physics: Conference Series PAPER OPEN ACCESS An inverse problem for a system of steady-state reaction-diffusion equations on a porous domain using a collage-based approach To cite this article:
More informationMaximum norm estimates for energy-corrected finite element method
Maximum norm estimates for energy-corrected finite element method Piotr Swierczynski 1 and Barbara Wohlmuth 1 Technical University of Munich, Institute for Numerical Mathematics, piotr.swierczynski@ma.tum.de,
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 information