Scientific Computing WS 2017/2018. Lecture 18. Jürgen Fuhrmann Lecture 18 Slide 1
|
|
- Logan Douglas
- 6 years ago
- Views:
Transcription
1 Scientific Computing WS 2017/2018 Lecture 18 Jürgen Fuhrmann Lecture 18 Slide 1
2 Lecture 18 Slide 2 Weak formulation of homogeneous Dirichlet problem Search u H0 1 (Ω) (here, tr u = 0) such that λ u v dx = fv dx v H0 1 (Ω) Then, Ω a(u, v) := Ω Ω λ u v dx is a self-adjoint bilinear form defined on the Hilbert space H0 1(Ω). It is bounded due to Cauchy-Schwarz: a(u, v) = λ u v dx u H 1 0 (Ω) v H 1 0 (Ω) Ω f (v) = Ω fv dx is a linear functional on H1 0 (Ω). For Hilbert spaces V the dual space V (the space of linear functionals) can be identified with the space itself. Lecture 14 Slide 23
3 Lecture 18 Slide 3 The Lax-Milgram lemma Let V be a Hilbert space. Let a : V V R be a self-adjoint bilinear form, and f a linear functional on V. Assume a is coercive, i.e. α > 0 : u V, a(u, u) α u 2 V. Then the problem: find u V such that a(u, v) = f (v) v V admits one and only one solution with an a priori estimate u V 1 α f V Lecture 14 Slide 24
4 Lecture 18 Slide 4 The Galerkin method II Let V be a Hilbert space. Let a : V V R be a self-adjoint bilinear form, and f a linear functional on V. Assume a is coercive with coercivity constant α, and continuity constant γ. Continuous problem: search u V such that a(u, v) = f (v) v V Let V h V be a finite dimensional subspace of V Discrete problem Galerkin approximation: Search u h V h such that a(u h, v h ) = f (v h ) v h V h By Lax-Milgram, this problem has a unique solution as well. Lecture 15 Slide 30
5 Lecture 18 Slide 5 Céa s lemma What is the connection between u and u h? Let v h V h be arbitrary. Then α u u h 2 a(u u h, u u h ) (Coercivity) = a(u u h, u v h ) + a(u u h, v h u h ) = a(u u h, u v h ) (Galerkin Orthogonality) γ u u h u v h (Boundedness) As a result u u h γ α inf v h V h u v h Up to a constant, the error of the Galerkin approximation is the error of the best approximation of the solution in the subspace V h. Lecture 15 Slide 31
6 Lecture 18 Slide 6 From the Galerkin method to the matrix equation Let φ 1... φ n be a set of basis functions of V h. Then, we have the representation u h = n j=1 u jφ j In order to search u h V h such that it is actually sufficient to require a(u h, v h ) = f (v h ) v h V h a(u h, φ i = f (φ i ) (i = 1... n) ( n ) a u j φ j, φ i = f (φ i ) (i = 1... n) j=1 n a(φ j, φ i )u j = f (φ i ) (i = 1... n) j=1 AU = F with A = (a ij ), a ij = a(φ i, φ j ), F = (f i ), f i = F (φ i ), U = (u i ). Matrix dimension is n n. Matrix sparsity? Lecture 15 Slide 32
7 Matrix properties a(, ) symmetric bilinear form A symmetric Coercivity A positive definite Lecture 18 Slide 7
8 Lecture 18 Slide 8 Barycentric coordinates Let K be a simplex. Functions λ i (i = 0... d): λ i : R d R x λ i (x) = 1 (x a i) n i (a j a i ) n i where a j is any vertex of K situated in F i. For x K, one has 1 (x a i) n i = (a j a i ) n i (x a i ) n i (a j a i ) n i (a j a i ) n i = (a j x) n i = dist(x, F i) (a j a i ) n i dist(a i, F i ) = dist(x, F i) F i /d dist(a i, F i ) F i /d = dist(x, F i) F i K i.e. λ i (x) is the ratio of the volume of the simplex K i (x) made up of x and the vertices of F i to the volume of K. Lecture 16 Slide 22
9 Lecture 18 Slide 9 Polynomial space P k Space of polynomials in x 1... x d of total degree k with real coefficients α i1...i d : Dimension: P k = p(x) = 0 i 1...i d k i 1+ +i d k α i1...i d x i x id d dim P k = ( ) d + k = k k + 1, d = (k + 1)(k + 2), d = 2 (k + 1)(k + 2)(k + 3), d = dim P 1 = d + 1 3, d = 1 dim P 2 = 6, d = 2 10, d = 3 Lecture 16 Slide 24
10 Lecture 18 Slide 10 P k simplex finite elements K: simplex spanned by a 0... a d in R d P = P k, such that s = dim P k For 0 i 0... i d k, i i d = k, let the set of nodes be defined by the points a i1...i d ;k with barycentric coordinates ( i0 k... id k ). Define Σ by σ i1...i d ;k(p) = p(a i1...i d ;k). Lecture 16 Slide 25
11 Lecture 18 Slide 11 Conformal triangulations Let T h be a subdivision of the polygonal domain Ω R d into non-intersecting compact simplices K m, m = 1... n e : Ω = n e m=1 Each simplex can be seen as the image of a affine transformation of a reference (e.g. unit) simplex K: K m K m = T m ( K) We assume that it is conformal, i.e. if K m, K n have a d 1 dimensional intersection F = K m K n, then there is a face F of K and renumberings of the vertices of K n, K m such that F = T m ( F ) = T n ( F ) and T m F = T n F Lecture 16 Slide 32
12 Lecture 18 Slide 12 Conformal triangulations II d = 1 : Each intersection F = K m K n is either empty or a common vertex d = 2 : Each intersection F = K m K n is either empty or a common vertex or a common edge d = 3 : Each intersection F = K m K n is either empty or a common vertex or a common edge or a common face Delaunay triangulations are conformal Lecture 16 Slide 33
13 H 1 -Conformal approximation using Lagrangian finite elemenents Approximation spacev h with zero jumps at element faces: V h = {v h C 0 (Ω) : m, v h Km P k } Zero jumps at interfaces: n, m, K m K n 0 (v h Km ) Km K n = (v h Kn ) Km K n V h H 1 (Ω). Lecture 18 Slide 13
14 Lecture 18 Slide 14 Zero jump at interfaces with Lagrangian finite elements Assume geometrically conformal mesh Assume all faces of K have the same number of nodes s For any face F = K 1 K 2 there are renumberings of the nodes of K 1 and K 2 such that for i = 1... s, a K1,i = a K2,i Then, v h K1 and v h K2 match at the interface K 1 K 2 if and only if they match at the common nodes v h K1 (a K1,i) = v h K2 (a K2,i) (i = 1... s ) Lecture 17 Slide 23
15 Lecture 18 Slide 15 Global degrees of freedom Let {a 1... a N } = {a K,1... a K,s } K T h Degree of freedom map j : T h {1... s} {1... N} (K, m) j(k, m) the global degree of freedom number Global shape functions φ 1,..., φ N W h defined by φ i K (a K,m ) = { δ mn if n {1... s} : j(k, n) = i 0 otherwise Global degrees of freedom γ 1,..., γ N : V h R defined by γ i (v h ) = v h (a i ) Lecture 17 Slide 24
16 Lecture 18 Slide 16 Lagrange finite element basis {φ 1,..., φ N } is a basis of V h, and γ 1... γ N is a basis of L(V h, R). Proof: {φ 1,..., φ N } are linearly independent: if N j=1 α jφ j = 0 then evaluation at a 1... a N yields that α 1... α N = 0. Let v h V h. It is single valued in a 1... a N. Let w h = N j=1 v h(a j )φ j. Then for all K T h, v h K and w h K coincide in the local nodes a K,1... a K,2, and by unisolvence, v h K = w h K. Lecture 17 Slide 25
17 Lecture 18 Slide 17 P 1 global shape functions Lecture 17 Slide 27
18 Lecture 18 Slide 18 P 2 global shape functions Node based Edge based Lecture 17 Slide 28
19 Lecture 18 Slide 19 Affine transformation estimates I K: reference element Let K T h. Affine mapping: T K : K K x J K x + b K with J K R d,d, b K R d, J K nonsingular Diameter of K: h K = max x1,x 2 K x 1 x 2 ρ K diameter of largest ball that can be inscribed into K σ K = hk ρ K : local shape regularity Lecture 17 Slide 31
20 Lecture 18 Slide 20 Shape regularity Now we discuss a family of meshes T h for h 0. We want to estimate global interpolation errors and see how they possibly diminuish For given T h, assume that h = max K Th h j A family of meshes is called shape regular if In 1D, σ K = 1 In 2D, σ K 2 sin θ K h, K T h, σ K = h K ρ K σ 0 where θ K is the smallest angle Lecture 17 Slide 37
21 Lecture 18 Slide 21 Error estimates for homogeneous Dirichlet problem Search u H0 1 (Ω) such that Ω λ u v dx = Ω fv dx v H 1 0 (Ω) Then, lim h 0 u u h 1,Ω = 0. If u H 2 (Ω) (e.g. on convex domains) then u u h 1,Ω ch u 2,Ω u u h 0,Ω ch 2 u 2,Ω Under certain conditions (convex domain, smooth coefficients) one has u u h 0,Ω ch u 1,Ω ( Aubin-Nitsche-Lemma ) Lecture 17 Slide 40
22 Lecture 18 Slide 22 H 2 -Regularity u H 2 (Ω) may be not fulfilled e.g. if Ω has re-entrant corners if on a smooth part of the domain, the boundary condition type changes if problem coefficients (λ) are discontinuos Situations differ as well between two and three space dimensions Delicate theory, ongoing research in functional analysis Consequence for simuations Deterioration of convergence rate Remedy: local refinement of the discretization mesh using a priori information using a posteriori error estimators + automatic refinement of discretizatiom mesh Lecture 17 Slide 41
23 More complicated integrals Assume non-constant right hand side f, space dependent heat conduction coefficient κ. Right hand side integrals f i = P 1 stiffness matrix elements a ij = K K f (x)λ i (x) dx κ(x) λ i λ j dx P k stiffness matrix elements created from higher order ansatz functions Lecture 18 Slide 23
24 Quadrature rules Quadrature rule: K g(x) dx K l q l=1 ω l g(ξ l ) ξ l : nodes, Gauss points ω l : weights The largest number k such that the quadrature is exact for polynomials of order k is called order k q of the quadrature rule, i.e. k k q, p P k Error estimate: φ C kq+1 (K), 1 K K K p(x) dx = K φ(x) dx l q l=1 ch kq+1 K l q l=1 ω l g(ξ l ) ω l p(ξ l ) sup x K, α =k q+1 α φ(x) Lecture 18 Slide 24
25 Some common quadrature rules Nodes are characterized by the barycentric coordinates d k q l q Nodes Weights ( 1, 1 ) , (1, 0), (0, 1) 3 2 ( ( 1, ), ( 1 + 3, 1 3 ), ( 1 3, ) 8, 5, ( 1, 1, 1 ) (1, 0, 0), (0, 1, 0), (0, 0, 1), 1, ( 1, 1, 0), ( 1, 0, 1 ), (0, 1, 1 ) 1, 1, ( 1, 1, 1 ), ( 1, 1, 3 ), ( 1, 3, 1 ), ( 3, 1, 1 ), 9, 25, 25, ( 1, 1, 1, 1 ) (1, 0, 0, 0), (0, 1, 0, 0), (0, 0, 1, 0), (0, 0, 0, 1) 2 4 ( 5 5, 5 5, 5 5, ) , 1, 1, , 1, 1, Lecture 18 Slide 25
26 Matching of approximation order and quadrature order Variational crime : instead of we solve a(u h, v h ) = f (v h ) v h V h a h (u h, v h ) = f h (v h ) v h V h where a h, f h are derived from their exact counterparts by quadrature For P 1 finite elements, zero order quadrature for volume integrals and first order quadrature for surface intergals is sufficient to keep the convergence order estimates stated before The rule of thumb for the volume quadrature is that the highest order terms must be evaluated exactly if the coefficients of the PDE are constant. Lecture 18 Slide 26
27 Practical realization of integrals Integral over barycentric coordinate function λ i (x) dx = 1 3 K K Right hand side integrals. Assume f (x) is given as a piecewise linear function with given values in the nodes of the triangulation f i = f (x)λ i (x) dx 1 3 K f (a i) K Integral over space dependent heat conduction coefficient: Assume κ(x) is given as a piecewise linear function with given values in the nodes of the triangulation a ij = κ(x) λ i λ j dx = 1 3 (κ(a 0)+κ(a 1 )+κ(a 2 )) λ i λ j dx K K Lecture 18 Slide 27
28 Stiffness matrix for Laplace operator for P1 FEM Element-wise calculation: a ij = a(φ i, φ j ) = φ i φ j dx = Standard assembly loop: for i, j = 1... N do set a ij = 0 end for K T h do for m,n=0... d do s mn = λ m λ n dx end end K Ω Ω a jdof (K,m),j dof (K,n) = a jdof (K,m),j dof (K,n) + s mn K T h φ i K φ j K dx Local stiffness matrix: S K = (s K;m,n ) = λ m λ n dx K Lecture 18 Slide 28
29 Local stiffness matrix calculation for P1 FEM a 0... a d : vertices of the simplex K, a K. Barycentric coordinates: λ j (a) = Kj (a) K For indexing modulo d+1 we can write K = 1 d! det ( a j+1 a j,... a j+d a j ) K j (a) = 1 d! det ( a j+1 a,... a j+d a ) From this information, we can calculate explicitely λ j (x) (which are constant vectors due to linearity) and the corresponding entries of the local stiffness s ij = λ i λ j dx K Lecture 18 Slide 29
30 Local stiffness matrix calculation for P1 FEM in 2D a 0 = (x 0, y 0 )... a d = (x 2, y 2 ): vertices of the simplex K, a = (x, y) K. Barycentric coordinates: λ j (x, y) = Kj (x,y) K For indexing modulo d+1 we can write K = 1 ( ) 2 det xj+1 x j x j+2 x j y j+1 y j y j+2 y j K j (x, y) = 1 ( ) 2 det xj+1 x x j+2 x y j+1 y y j+2 y Therefore, we have K j (x, y) = 1 2 ((x j+1 x)(y j+2 y) (x j+2 x)(y j+1 y)) x K j (x, y) = 1 2 ((y j+1 y) (y j+2 y)) = 1 2 (y j+1 y j+2 ) y K j (x, y) = 1 2 ((x j+2 x) (x j+1 x)) = 1 2 (x j+2 x j+1 ) Lecture 18 Slide 30
31 Local stiffness matrix calculation for P1 FEM in 2D II s ij = K λ i λ j dx = ( ) x1 x So, let V = 0 x 2 x 0 y 1 y 0 y 2 y 0 Then K ( ) ( ) yj+1 y yi+1 y 4 K 2 i+2, x i+2 x j+2 i+1 x j+2 x j+1 x 1 x 2 = V 00 V 01 and y 1 y 2 = V 10 V 11 ( ) ( ) y1 y 2 K λ 0 = 2 V10 V = 11 x 2 x 1 V 01 V 00 ( ) ( ) y2 y 2 K λ 1 = 0 V11 = x 0 x 2 V 01 ( ) ( ) y0 y 2 K λ 2 = 1 V10 = x 1 x 0 V 00 Lecture 18 Slide 31
32 Degree of freedom map representation for P1 finite elements List of global nodes a 0... a N : two dimensional array of coordinate values with N rows and d columns Local-global degree of freedom map: two-dimensional array C of index values with N el rows and d + 1 columns such that C(i, m) = j dof (K i, m). The mesh generator triangle generates this information directly Lecture 18 Slide 32
33 Practical realization of boundary conditions Robin boundary value problem κ u = f κ u + α(u g) = 0 in Ω on Ω Weak formulation: search u H 1 (Ω) such that Ω κ u v dx + αuv ds = fv dx + αgv ds v H 1 (Ω) Ω Ω Ω In 2D, for P 1 FEM, boundary integrals can be calculated by trapezoidal rule without sacrificing approximation order Use Dirichlet penalty method to handle Dirichlet boundary conditions Lecture 18 Slide 33
34 Next lecture Next lecture: Jan. 9, Happy holidays! Lecture 18 Slide 34
Scientific 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 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 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 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 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 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 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 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 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 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 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 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 informationPREPRINT 2010:25. Fictitious domain finite element methods using cut elements: II. A stabilized Nitsche method ERIK BURMAN PETER HANSBO
PREPRINT 2010:25 Fictitious domain finite element methods using cut elements: II. A stabilized Nitsche method ERIK BURMAN PETER HANSBO Department of Mathematical Sciences Division of Mathematics CHALMERS
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 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 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 informationSolution of Non-Homogeneous Dirichlet Problems with FEM
Master Thesis Solution of Non-Homogeneous Dirichlet Problems with FEM Francesco Züger Institut für Mathematik Written under the supervision of Prof. Dr. Stefan Sauter and Dr. Alexander Veit August 27,
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 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 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 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 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 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 informationChapter 1 Foundations of Elliptic Boundary Value Problems 1.1 Euler equations of variational problems
Chapter 1 Foundations of Elliptic Boundary Value Problems 1.1 Euler equations of variational problems Elliptic boundary value problems often occur as the Euler equations of variational problems the latter
More informationVariational Formulations
Chapter 2 Variational Formulations In this chapter we will derive a variational (or weak) formulation of the elliptic boundary value problem (1.4). We will discuss all fundamental theoretical results that
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 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 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 informationBuilding a finite element program in MATLAB Linear elements in 1d and 2d
Building a finite element program in MATLAB Linear elements in 1d and 2d D. Peschka TU Berlin Supplemental material for the course Numerische Mathematik 2 für Ingenieure at the Technical University Berlin,
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 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 informationApplied/Numerical Analysis Qualifying Exam
Applied/Numerical Analysis Qualifying Exam August 9, 212 Cover Sheet Applied Analysis Part Policy on misprints: The qualifying exam committee tries to proofread exams as carefully as possible. Nevertheless,
More informationYour first day at work MATH 806 (Fall 2015)
Your first day at work MATH 806 (Fall 2015) 1. Let X be a set (with no particular algebraic structure). A function d : X X R is called a metric on X (and then X is called a metric space) when d satisfies
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 informationNumerical Analysis of Higher Order Discontinuous Galerkin Finite Element Methods
Numerical Analysis of Higher Order Discontinuous Galerkin Finite Element Methods Contents Ralf Hartmann Institute of Aerodynamics and Flow Technology DLR (German Aerospace Center) Lilienthalplatz 7, 3808
More informationApproximation in Banach Spaces by Galerkin Methods
2 Approximation in Banach Spaces by Galerkin Methods In this chapter, we consider an abstract linear problem which serves as a generic model for engineering applications. Our first goal is to specify the
More informationLecture Notes: African Institute of Mathematics Senegal, January Topic Title: A short introduction to numerical methods for elliptic PDEs
Lecture Notes: African Institute of Mathematics Senegal, January 26 opic itle: A short introduction to numerical methods for elliptic PDEs Authors and Lecturers: Gerard Awanou (University of Illinois-Chicago)
More 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 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 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 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 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 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 informationSubdiffusion in a nonconvex polygon
Subdiffusion in a nonconvex polygon Kim Ngan Le and William McLean The University of New South Wales Bishnu Lamichhane University of Newcastle Monash Workshop on Numerical PDEs, February 2016 Outline Time-fractional
More informationTMA4220: Programming project - part 1
TMA4220: Programming project - part 1 TMA4220 - Numerical solution of partial differential equations using the finite element method The programming project will be split into two parts. This is the first
More informationThe Mortar Boundary Element Method
The Mortar Boundary Element Method A Thesis submitted for the degree of Doctor of Philosophy by Martin Healey School of Information Systems, Computing and Mathematics Brunel University March 2010 Abstract
More informationNUMERICAL PARTIAL DIFFERENTIAL EQUATIONS
NUMERICAL PARTIAL DIFFERENTIAL EQUATIONS Lecture Notes, Winter 2011/12 Christian Clason January 31, 2012 Institute for Mathematics and Scientific Computing Karl-Franzens-Universität Graz CONTENTS I BACKGROUND
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 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 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 informationDIFFERENTIAL EQUATIONS
NUMERICAL PARTIAL DIFFERENTIAL EQUATIONS Lecture Notes, Winter 2013/14 Christian Clason July 15, 2014 Institute for Mathematics and Scientific Computing Karl-Franzens-Universität Graz CONTENTS I BACKGROUND
More informationLecture Note III: Least-Squares Method
Lecture Note III: Least-Squares Method Zhiqiang Cai October 4, 004 In this chapter, we shall present least-squares methods for second-order scalar partial differential equations, elastic equations of solids,
More 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 informationFinite Element Methods for Fourth Order Variational Inequalities
Louisiana State University LSU Digital Commons LSU Doctoral Dissertations Graduate School 2013 Finite Element Methods for Fourth Order Variational Inequalities Yi Zhang Louisiana State University and Agricultural
More informationNumerical analysis and a-posteriori error control for a new nonconforming quadrilateral linear finite element
Numerical analysis and a-posteriori error control for a new nonconforming quadrilateral linear finite element M. Grajewski, J. Hron, S. Turek Matthias.Grajewski@math.uni-dortmund.de, jaroslav.hron@math.uni-dortmund.de,
More informationMultigrid Methods for Saddle Point Problems
Multigrid Methods for Saddle Point Problems Susanne C. Brenner Department of Mathematics and Center for Computation & Technology Louisiana State University Advances in Mathematics of Finite Elements (In
More informationEfficient hp-finite elements
Efficient hp-finite elements Ammon Washburn July 31, 2015 Abstract Ways to make an hp-finite element method efficient are presented. Standard FEMs and hp-fems with their various strengths and weaknesses
More informationDISCRETE MAXIMUM PRINCIPLES in THE FINITE ELEMENT SIMULATIONS
DISCRETE MAXIMUM PRINCIPLES in THE FINITE ELEMENT SIMULATIONS Sergey Korotov BCAM Basque Center for Applied Mathematics http://www.bcamath.org 1 The presentation is based on my collaboration with several
More informationA non-standard Finite Element Method based on boundary integral operators
A non-standard Finite Element Method based on boundary integral operators Clemens Hofreither Ulrich Langer Clemens Pechstein June 30, 2010 supported by Outline 1 Method description Motivation Variational
More informationPartial Differential Equations and the Finite Element Method
Partial Differential Equations and the Finite Element Method Pavel Solin The University of Texas at El Paso Academy of Sciences ofthe Czech Republic iwiley- INTERSCIENCE A JOHN WILEY & SONS, INC, PUBLICATION
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 informationChapter 6. Finite Element Method. Literature: (tiny selection from an enormous number of publications)
Chapter 6 Finite Element Method Literature: (tiny selection from an enormous number of publications) K.J. Bathe, Finite Element procedures, 2nd edition, Pearson 2014 (1043 pages, comprehensive). Available
More informationACM/CMS 107 Linear Analysis & Applications Fall 2017 Assignment 2: PDEs and Finite Element Methods Due: 7th November 2017
ACM/CMS 17 Linear Analysis & Applications Fall 217 Assignment 2: PDEs and Finite Element Methods Due: 7th November 217 For this assignment the following MATLAB code will be required: Introduction http://wwwmdunloporg/cms17/assignment2zip
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 informationSpace-time sparse discretization of linear parabolic equations
Space-time sparse discretization of linear parabolic equations Roman Andreev August 2, 200 Seminar for Applied Mathematics, ETH Zürich, Switzerland Support by SNF Grant No. PDFMP2-27034/ Part of PhD thesis
More information2 Elliptic Differential Equations
2 Elliptic Differential Equations 2.1 Classical solutions As far as existence and uniqueness results for classical solutions are concerned, we restrict ourselves to linear elliptic second order elliptic
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 informationVirtual Element Methods for general second order elliptic problems
Virtual Element Methods for general second order elliptic problems Alessandro Russo Department of Mathematics and Applications University of Milano-Bicocca, Milan, Italy and IMATI-CNR, Pavia, Italy Workshop
More 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 informationTechnische Universität Graz
Technische Universität Graz Error Estimates for Neumann Boundary Control Problems with Energy Regularization T. Apel, O. Steinbach, M. Winkler Berichte aus dem Institut für Numerische Mathematik Bericht
More informationAn Extended Finite Element Method for a Two-Phase Stokes problem
XFEM project An Extended Finite Element Method for a Two-Phase Stokes problem P. Lederer, C. Pfeiler, C. Wintersteiger Advisor: Dr. C. Lehrenfeld August 5, 2015 Contents 1 Problem description 2 1.1 Physics.........................................
More informationNumerical Solution of Partial Differential Equations
Numerical Solution of Partial Differential Equations Prof. Ralf Hiptmair, Prof. Christoph Schwab und Dr. H. Harbrecht V1.0: summer term 2004, V2.0: winter term 2005/2006 Draft version October 26, 2005
More informationAn introduction to the mathematical theory of finite elements
Master in Seismic Engineering E.T.S.I. Industriales (U.P.M.) Discretization Methods in Engineering An introduction to the mathematical theory of finite elements Ignacio Romero ignacio.romero@upm.es October
More informationGraded Mesh Refinement and Error Estimates of Higher Order for DGFE-solutions of Elliptic Boundary Value Problems in Polygons
Graded Mesh Refinement and Error Estimates of Higher Order for DGFE-solutions of Elliptic Boundary Value Problems in Polygons Anna-Margarete Sändig, Miloslav Feistauer University Stuttgart, IANS Journées
More informationPREPRINT 2010:23. A nonconforming rotated Q 1 approximation on tetrahedra PETER HANSBO
PREPRINT 2010:23 A nonconforming rotated Q 1 approximation on tetrahedra PETER HANSBO Department of Mathematical Sciences Division of Mathematics CHALMERS UNIVERSITY OF TECHNOLOGY UNIVERSITY OF GOTHENBURG
More informationSpline Element Method for Partial Differential Equations
for Partial Differential Equations Department of Mathematical Sciences Northern Illinois University 2009 Multivariate Splines Summer School, Summer 2009 Outline 1 Why multivariate splines for PDEs? Motivation
More informationDiscontinuous Galerkin Methods
Discontinuous Galerkin Methods Joachim Schöberl May 20, 206 Discontinuous Galerkin (DG) methods approximate the solution with piecewise functions (polynomials), which are discontinuous across element interfaces.
More informationA COLLOCATION METHOD FOR SOLVING THE EXTERIOR NEUMANN PROBLEM
TUDIA UNIV. BABEŞ BOLYAI, MATHEMATICA, Volume XLVIII, Number 3, eptember 2003 A COLLOCATION METHOD FOR OLVING THE EXTERIOR NEUMANN PROBLEM ANDA MICULA Dedicated to Professor Gheorghe Micula at his 60 th
More informationICES REPORT Direct Serendipity Finite Elements on Convex Quadrilaterals
ICES REPORT 17-8 October 017 Direct Serendipity Finite Elements on Convex Quadrilaterals by Todd Arbogast and Zhen Tao The Institute for Computational Engineering and Sciences The University of Texas at
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 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 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 informationCHAPTER 3 Further properties of splines and B-splines
CHAPTER 3 Further properties of splines and B-splines In Chapter 2 we established some of the most elementary properties of B-splines. In this chapter our focus is on the question What kind of functions
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 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 BIVARIATE SPLINE METHOD FOR SECOND ORDER ELLIPTIC EQUATIONS IN NON-DIVERGENCE FORM
A BIVARIAE SPLINE MEHOD FOR SECOND ORDER ELLIPIC EQUAIONS IN NON-DIVERGENCE FORM MING-JUN LAI AND CHUNMEI WANG Abstract. A bivariate spline method is developed to numerically solve second order elliptic
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 informationarxiv: v5 [math.na] 1 Sep 2018
High-Order Adaptive Extended Stencil FEM (AES-FEM) Part I: Convergence and Superconvergence Xiangmin Jiao Rebecca Conley Tristan J. Delaney arxiv:1603.09325v5 [math.na] 1 Sep 2018 Abstract Finite elements
More informationTD 1: Hilbert Spaces and Applications
Université Paris-Dauphine Functional Analysis and PDEs Master MMD-MA 2017/2018 Generalities TD 1: Hilbert Spaces and Applications Exercise 1 (Generalized Parallelogram law). Let (H,, ) be a Hilbert space.
More informationLocal radial basis function approximation on the sphere
Local radial basis function approximation on the sphere Kerstin Hesse and Q. T. Le Gia School of Mathematics and Statistics, The University of New South Wales, Sydney NSW 05, Australia April 30, 007 Abstract
More informationNumerical Integration in Meshfree Methods
Numerical Integration in Meshfree Methods Pravin Madhavan New College University of Oxford A thesis submitted for the degree of Master of Science in Mathematical Modelling and Scientific Computing Trinity
More informationarxiv: v2 [math.ag] 24 Jun 2015
TRIANGULATIONS OF MONOTONE FAMILIES I: TWO-DIMENSIONAL FAMILIES arxiv:1402.0460v2 [math.ag] 24 Jun 2015 SAUGATA BASU, ANDREI GABRIELOV, AND NICOLAI VOROBJOV Abstract. Let K R n be a compact definable set
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 informationFINITE ELEMENT METHODS
FINITE ELEMENT METHODS Lecture notes arxiv:1709.08618v1 [math.na] 25 Sep 2017 Christian Clason September 25, 2017 christian.clason@uni-due.de https://udue.de/clason CONTENTS I BACKGROUND 1 overview of
More informationConsistency Estimates for gfd Methods and Selection of Sets of Influence
Consistency Estimates for gfd Methods and Selection of Sets of Influence Oleg Davydov University of Giessen, Germany Localized Kernel-Based Meshless Methods for PDEs ICERM / Brown University 7 11 August
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 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 informationSYLLABUS. 1 Linear maps and matrices
Dr. K. Bellová Mathematics 2 (10-PHY-BIPMA2) SYLLABUS 1 Linear maps and matrices Operations with linear maps. Prop 1.1.1: 1) sum, scalar multiple, composition of linear maps are linear maps; 2) L(U, V
More informationA FINITE ELEMENT METHOD FOR ELLIPTIC EQUATIONS ON SURFACES
A FINITE ELEMENT METHOD FOR ELLIPTIC EQUATIONS ON SURFACES MAXIM A. OLSHANSKII, ARNOLD REUSKEN, AND JÖRG GRANDE Abstract. In this paper a new finite element approach for the discretization of elliptic
More informationTwo-Scale Composite Finite Element Method for Dirichlet Problems on Complicated Domains
Numerische Mathematik manuscript No. will be inserted by the editor) Two-Scale Composite Finite Element Method for Dirichlet Problems on Complicated Domains M. Rech, S. Sauter, A. Smolianski Institut für
More information