WRT in 2D: Poisson Example
|
|
- Maximillian Lawrence
- 5 years ago
- Views:
Transcription
1 WRT in 2D: Poisson Example Consider 2 u f on [, L x [, L y with u. WRT: For all v X N, find u X N a(v, u) such that v u dv v f dv. Follows from strong form plus integration by parts: ( ) 2 u v + 2 u dx dy }{{ 2 }}{{ 2 } I x I y v f dv. First integral: I x : Ly [ Lx Ly [ Lx v 2 u dx dy 2 v u dx v(x, y) u xl x,y x,y } {{ } dy Ly Lx v u dx dy. Second integral, by symmetry: I y Ly Lx v u dx dy.
2 Evaluating Integrals Using Basis Functions For all v, u X N, v(x, y) u(x, y) N N v ij φ i (x) ψ j (y) i j N N u pq φ p (x) ψ q (y). p q Here, we take φ and ψ to be the bases on [, L x [, L y, φ i (x) : l i (r(x)), r(x) : + 2 L x x, ψ j (y) : l j (s(y)), s(y) : + 2 L y y, and the l i s and l j s are the (same) Lagrange interpolating polynomials based on the Gauss-Lobatto-Legendre (GLL) points r j [, : ˆ. Derivative with respect to x: Derivative with respect to y: N v i N u N φ i(x) ψ j (y) v ij j p q N v i N u N φ p(x) ψ q (y) u pq. N φ i (x) ψ j(y) v ij j p q N φ p (x) ψ q(y) u pq.
3 First integral, I x Ly Lx pq pq pq ( ) ( v ij φ i(x) ψ j (y)v ij pq ij Ly v ij ij ( Ly v ij ij v ij By,jq Ā x,ip u pq ij v T ( By Āx) ū. Second integral, by symmetry, I y v T ( Ā y B x ) ū. [ Lx ψ j (y) ψ q (y) ) ( Lx ψ j (y) ψ q (y) dy u pq φ p(x) ψ q (y)u pq ) φ i(x) φ p(x) dx dy u pq ) φ i(x) φ p(x) dx u pq, dx dy
4 Following as we did in D, we end up with the linear system Au Bf A B y A x + A y B x. Let, as in D, Ā x : 2 L x ˆDT ˆB ˆD, [Āx ip : Lx φ i(x) φ p(x) dx N φ i(x k ) k [ N 2 k L x 2 L x N k 2 L x N k ( ) Lx 2 ρ k φ p(x k ) dl i dr [ 2 L x ˆDT ˆB ˆD rk (Lx ˆD ki ρ k ˆDkp ˆD ki ˆBkk ˆDkp. ip 2 ρ k ) [ 2 L x (GLL quadrature) dl p dr rk This gives us the bar matrix, Ā x, which includes the endpoints.
5 To get the interior points only, we (re-)introduce the restriction matrix, R, and prolongation matrix, P : R T. Recall P T I T lr (N+) (N ). Q: What is the role of P in a matrix-vector product? Q: What is the role of R in a matrix-vector product?
6 Matlab Example wrttest.m N5; % Polynomial Order [Ah,Bh,Ch,Dh,z,w semhat(n); R speye(n+); RR(2:(end-),:); % Dirichlet Restriction Matrix A R*Ah*R ; f+*z; u A \ (R*(Bh*f)); u[;u;; ue.5*(-z.*z); err max(abs(u-ue)) After executing this code we can look at R: R full(r) R Ah 5.667e e e e-.97e e e e e e e-.97e e e+ 6.4e+ -3.4e e e e e- -3.4e+ 6.4e e e-.97e e e e e e e-2.97e e e e e+ >> A A 9.267e e e e e+ 6.4e+ -3.4e e e- -3.4e+ 6.4e e e e e e+ >> Bh Bh e e e e e e-2 >> R*Bh*R ans e e e e- We see that R does just what we need. It picks out the active rows (i.e., equations) in the ODE/PDE.
7 R T u extends u by zero, so that we can apply quadrature and/or differentiate u(x) on x, x,..., x N, as needed. >> u[ u >> R *u ans 2 3 4
8 Note: If we have Neumann conditions in the y-direction (say), then y-indices range from to N. For example, LHS of the ij-th equation will be I N q N p [B y A x + A y B x u A u. (B y,jq A x,ip + A y,jq B x,ip ) u pq { i,..., N j,..., N The mass matrices, B x, B y, and stiffness matrices, A x, A y, will have the correct size (and hence, correct index range in the first expression), if we set B x R x BR T x A x R x ĀR T x B y R y BR T y A y R y ĀR T y, with R T x T I T lr (N+) (N ), R T y I lr (N+) (N+). We typically refer to Ā and B as the Neumann operators. They are the full (N + ) (N + ) matrices corresponding to the Neumann problem. In terms of the generic ( hat ) matrices, they are (for x, say) Ā 2 L x Â, B L x 2 ˆB. Example: Solve 2 u on [, 2 [, u(x, y) at x, 2. u at y,.
9 Variable Coefficients in 2D: Recall D example: d dx q(x)du dx { u() f(x) + BCs u () WRT: v d dx q(x)du dx dx dv dx q(x)du dx dx v f(x) dx v f(x) dx (Dv) T QBDu v T Bf v T D T QBDu v T Bf Au Bf A : D T QBD. With (highly accurate) GLL quadrature, it is sufficient to evaluate q(x) at the GLL nodal points, x k, k,..., N The matrices Q and B are thus both diagonal: Q kk q(x k ), B kk ρ k, k,..., N.
10 2D example: [q(x, y) u f(x, y) + BCs v [q u dx dy v f(x, y) dx dy v [ q u + q u dx dy v f(x, y) dx dy First term on left, itegrated by parts [ Ly Lx I x Let Γ i denote the ith faces of, i: left face i2: right face i3: lower face i4: upper face v q u u dx v q Let ˆn i be the associated outward-pointing unit-normals on these faces: i: ˆn (ˆn x,, ˆn y, ) (, ) i2: ˆn 2 (ˆn x,2, ˆn y,2 ) (, ) i3: ˆn 3 (ˆn x,3, ˆn y,3 ) (, ) i4: ˆn 4 (ˆn x,4, ˆn y,4 ) (, ) We can rewrite I x as: I x Ly Lx v q u 2 dx dy i x x Γ i v q u ˆn x,i dy. dy Similarly, I y is I y Ly Lx v q u 4 dx dy i3 Γ i v q u ˆn y,i dy. Combining I x and I y, I I x + I y Ly Lx v q u + v q u 4 dx dy v (q u) dv v q u ˆn ds. i Γ i v q ( u ˆn x,i + v q u ) ˆn y,i ds.
11 If v q u ˆn, then u satisfies v (q u) dv v f dv for all v X N. Note: For any u X N we have the gradient ) u u u ( Dx D y u I 2 L x ˆD ˆD I 2 L y u. Example: x L x 2 (z + ), y L y (z + ) 2 [X, Y ndgrid(x, y) U sin πx U x D x U, U y D y U D x Î 2 ˆD, Dy 2 ˆD Î L x L y Plot U x, U y.
12 Returning to our weighted residual form, insertion of the bases and using quadrature yields I Ly Lx v q u + v q u dx dy (D x v) T Q B D x u + (D y v) T Q B D y u v T D T x Q B D x u + v T D T y Q B D y u. Note that if the range of indices on v and u includes all nodes, including the boundaries, then we would have the form of our Neumann operator, Ā D T x Q B D x + D T y Q B D y. To construct the A-matrix that governs the interior degrees of freedom, apply the restriction operator R to Ā: A R Ā RT. In the case of a rectangular domain, R restricts index sets in both the x and y directions and it can be written in tensor-product form: R R y R x. As always, the R x matrix is on the right, because it is acting on the fastest changing index in the u T [u u... u NN T lexicographical ordering of the unknowns.
13 Gradient-Based Form It is sometimes convenient to write Ā in a higher-level form. We can also write Ā [ Dx D y T [ QB QB [ Dx D y D T Q D. Here, we use the bold font to indicate that the matrices are working on or producing vector fields, i.e., inputs and/or outputs that have two components in 2D, three components in 3D, etc. Obviously, the D matrix is the discrete gradient operator introduced above.
14 Deformed Geometries If is deformed, ˆ x(r, s) ij x ij l i (r) l j (s) N N x ij l i (r) l j (s). j i Derivative is computed using the chain rule [w pq u u rps q + u (r x ) pq (u r ) pq + (s x ) pq (u s ) pq w [r x (Î ˆD)u + [s x ( ˆD Î)u. The full gradient (vector field) is given by, w u [ [rx [s x [r y [s y }{{} R x [ Dr D s u R x D r u. Here, we use bold font to indicate a vector field or a matrix that operates on and/or produces a vector field.
15 Metrics What is? rpsq pq ij pq ij dl i x ij dr l j (s q ) D r x (Î ˆD)x p x ij l i (r p ) dl j ds q D s x ( ˆD Î)x. Chain Rule: Set u x: u u + u At each point, (r p, s q ): [ [ [. Thus we find the inverse metrics, i j, by inverting the 2 2 matrix of computable derivatives i j, at each gridpoint (r p, s q ).
16 Integration in Similarly, we must transform f da ˆ f J dr ds. Here, J is the Jacobian: It is a scalar field corresponding to the amount of area in that is associated with a unit area in ˆ. An example is shown in the figure below. The blue area is given by the cross product, da dr ds J dr ds [ dr ds dr ds. ˆ That is, J pq is the determinant of the 2 2 metric tensor at each quadrature point. Mass Matrix: We note that the mass matrix in the deformed geometry is diagonal, evaluated at the GLL points: B diag(ρ p ρ q J pq J ( ˆB ˆB), if we view J as a diagonal matrix, following standard practice with GLL quadrature.
17 Bilinear Form in Deformed Coordinates As usual, we evaluate the bilinear form, a(v, u) v u da [Dv T (B) [Du [R x D r v T (B) [R x D r u v T ( Dr D s ) T [ [rx [s x [r y [s y T ( J ( ˆB ˆB) J ( ˆB ˆB) ) [ [rx [s x [r y [s y ( Dr D s ) u. v T D T r GDu. In d space dimensions, ( d ) i j G ij : J B, i,..., d, j,..., d, d 2 or 3. k k k Example, d 2: [ Grr G rs [ G G 2 G G sr G ss G 2 G 22. Note: G ij G ji. Each G.. is a diagonal matrix with (N + ) 2 (or, (N + ) d ) entries.
Simple 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 informationMA8502 Numerical solution of partial differential equations. The Poisson problem: Mixed Dirichlet/Neumann boundary conditions along curved boundaries
MA85 Numerical solution of partial differential equations The Poisson problem: Mied Dirichlet/Neumann boundar conditions along curved boundaries Fall c Einar M. Rønquist Department of Mathematical Sciences
More informationChapter 1: Systems of linear equations and matrices. Section 1.1: Introduction to systems of linear equations
Chapter 1: Systems of linear equations and matrices Section 1.1: Introduction to systems of linear equations Definition: A linear equation in n variables can be expressed in the form a 1 x 1 + a 2 x 2
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 informationFinal Exam May 4, 2016
1 Math 425 / AMCS 525 Dr. DeTurck Final Exam May 4, 2016 You may use your book and notes on this exam. Show your work in the exam book. Work only the problems that correspond to the section that you prepared.
More informationA High-Order Galerkin Solver for the Poisson Problem on the Surface of the Cubed Sphere
A High-Order Galerkin Solver for the Poisson Problem on the Surface of the Cubed Sphere Michael Levy University of Colorado at Boulder Department of Applied Mathematics August 10, 2007 Outline 1 Background
More informationMath 5630: Iterative Methods for Systems of Equations Hung Phan, UMass Lowell March 22, 2018
1 Linear Systems Math 5630: Iterative Methods for Systems of Equations Hung Phan, UMass Lowell March, 018 Consider the system 4x y + z = 7 4x 8y + z = 1 x + y + 5z = 15. We then obtain x = 1 4 (7 + y z)
More informationQualifying Examination
Summer 24 Day. Monday, September 5, 24 You have three hours to complete this exam. Work all problems. Start each problem on a All problems are 2 points. Please email any electronic files in support of
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 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 informationWeak form of Boundary Value Problems. Simulation Methods in Acoustics
Weak form of Boundary Value Problems Simulation Methods in Acoustics Note on finite dimensional description of functions Approximation: N f (x) ˆf (x) = q j φ j (x) j=1 Residual function: r(x) = f (x)
More informationInverses. Stephen Boyd. EE103 Stanford University. October 28, 2017
Inverses Stephen Boyd EE103 Stanford University October 28, 2017 Outline Left and right inverses Inverse Solving linear equations Examples Pseudo-inverse Left and right inverses 2 Left inverses a number
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 informationRecall: Dot product on R 2 : u v = (u 1, u 2 ) (v 1, v 2 ) = u 1 v 1 + u 2 v 2, u u = u u 2 2 = u 2. Geometric Meaning:
Recall: Dot product on R 2 : u v = (u 1, u 2 ) (v 1, v 2 ) = u 1 v 1 + u 2 v 2, u u = u 2 1 + u 2 2 = u 2. Geometric Meaning: u v = u v cos θ. u θ v 1 Reason: The opposite side is given by u v. u v 2 =
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 information8 A pseudo-spectral solution to the Stokes Problem
8 A pseudo-spectral solution to the Stokes Problem 8.1 The Method 8.1.1 Generalities We are interested in setting up a pseudo-spectral method for the following Stokes Problem u σu p = f in Ω u = 0 in Ω,
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 informationAlgorithms for Scientific Computing
Algorithms for Scientific Computing Finite Element Methods Michael Bader Technical University of Munich Summer 2016 Part I Looking Back: Discrete Models for Heat Transfer and the Poisson Equation Modelling
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 informationHIGHER-ORDER THEORIES
HIGHER-ORDER THEORIES THIRD-ORDER SHEAR DEFORMATION PLATE THEORY LAYERWISE LAMINATE THEORY J.N. Reddy 1 Third-Order Shear Deformation Plate Theory Assumed Displacement Field µ u(x y z t) u 0 (x y t) +
More informationCSL361 Problem set 4: Basic linear algebra
CSL361 Problem set 4: Basic linear algebra February 21, 2017 [Note:] If the numerical matrix computations turn out to be tedious, you may use the function rref in Matlab. 1 Row-reduced echelon matrices
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 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 information3. Numerical integration
3. Numerical integration... 3. One-dimensional quadratures... 3. Two- and three-dimensional quadratures... 3.3 Exact Integrals for Straight Sided Triangles... 5 3.4 Reduced and Selected Integration...
More informationLecture 21: Isoparametric Formulation of Plane Elements.
6.6. Rectangular Plane Stress/Strain Element. he CS element uses a triangular shape. he 3 nodes of the CS element allow us to employ linear multivariate approximations to the u and v displacements. he
More informationMultivariable Calculus
2 Multivariable Calculus 2.1 Limits and Continuity Problem 2.1.1 (Fa94) Let the function f : R n R n satisfy the following two conditions: (i) f (K ) is compact whenever K is a compact subset of R n. (ii)
More informationClassical Field Theory: Electrostatics-Magnetostatics
Classical Field Theory: Electrostatics-Magnetostatics April 27, 2010 1 1 J.D.Jackson, Classical Electrodynamics, 2nd Edition, Section 1-5 Electrostatics The behavior of an electrostatic field can be described
More informationTransformations from R m to R n.
Transformations from R m to R n 1 Differentiablity First of all because of an unfortunate combination of traditions (the fact that we read from left to right and the way we define matrix multiplication
More information12. Stresses and Strains
12. Stresses and Strains Finite Element Method Differential Equation Weak Formulation Approximating Functions Weighted Residuals FEM - Formulation Classification of Problems Scalar Vector 1-D T(x) u(x)
More informationMathematical Tripos Part IA Lent Term Example Sheet 1. Calculate its tangent vector dr/du at each point and hence find its total length.
Mathematical Tripos Part IA Lent Term 205 ector Calculus Prof B C Allanach Example Sheet Sketch the curve in the plane given parametrically by r(u) = ( x(u), y(u) ) = ( a cos 3 u, a sin 3 u ) with 0 u
More informationA Fast Fourier transform based direct solver for the Helmholtz problem. AANMPDE 2017 Palaiochora, Crete, Oct 2, 2017
A Fast Fourier transform based direct solver for the Helmholtz problem Jari Toivanen Monika Wolfmayr AANMPDE 2017 Palaiochora, Crete, Oct 2, 2017 Content 1 Model problem 2 Discretization 3 Fast solver
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 information2.20 Fall 2018 Math Review
2.20 Fall 2018 Math Review September 10, 2018 These notes are to help you through the math used in this class. This is just a refresher, so if you never learned one of these topics you should look more
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 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 informationOR MSc Maths Revision Course
OR MSc Maths Revision Course Tom Byrne School of Mathematics University of Edinburgh t.m.byrne@sms.ed.ac.uk 15 September 2017 General Information Today JCMB Lecture Theatre A, 09:30-12:30 Mathematics revision
More informationVector Calculus. A primer
Vector Calculus A primer Functions of Several Variables A single function of several variables: f: R $ R, f x (, x ),, x $ = y. Partial derivative vector, or gradient, is a vector: f = y,, y x ( x $ Multi-Valued
More informationMath 340 Final Exam December 16, 2006
Math 34 Final Exam December 6, 6. () Suppose A 3 4. a) Find the row-reduced echelon form of A. 3 4 so the row reduced echelon form is b) What is rank(a)? 3 4 4 The rank is two since there are two pivots.
More informationCIV-E1060 Engineering Computation and Simulation Examination, December 12, 2017 / Niiranen
CIV-E16 Engineering Computation and Simulation Examination, December 12, 217 / Niiranen This examination consists of 3 problems rated by the standard scale 1...6. Problem 1 Let us consider a long and tall
More informationONE - DIMENSIONAL INTEGRATION
Chapter 4 ONE - DIMENSIONA INTEGRATION 4. Introduction Since the finite element method is based on integral relations it is logical to expect that one should strive to carry out the integrations as efficiently
More informationPartial Differential Equations
Partial Differential Equations Xu Chen Assistant Professor United Technologies Engineering Build, Rm. 382 Department of Mechanical Engineering University of Connecticut xchen@engr.uconn.edu Contents 1
More informationChapter 4. The First Fundamental Form (Induced Metric)
Chapter 4. The First Fundamental Form (Induced Metric) We begin with some definitions from linear algebra. Def. Let V be a vector space (over IR). A bilinear form on V is a map of the form B : V V IR which
More informationSolution to Homework 1
Solution to Homework Sec 2 (a) Yes It is condition (VS 3) (b) No If x, y are both zero vectors Then by condition (VS 3) x = x + y = y (c) No Let e be the zero vector We have e = 2e (d) No It will be false
More informationProperties of Transformations
6. - 6.4 Properties of Transformations P. Danziger Transformations from R n R m. General Transformations A general transformation maps vectors in R n to vectors in R m. We write T : R n R m to indicate
More informationG H. Extended Unit Tests A L L. Higher Still Advanced Higher Mathematics. (more demanding tests covering all levels) Contents. 3 Extended Unit Tests
M A T H E M A T I C S H I G H E R Higher Still Advanced Higher Mathematics S T I L L Extended Unit Tests A (more demanding tests covering all levels) Contents Extended Unit Tests Detailed marking schemes
More informationChapter 3. Differentiable Mappings. 1. Differentiable Mappings
Chapter 3 Differentiable Mappings 1 Differentiable Mappings Let V and W be two linear spaces over IR A mapping L from V to W is called a linear mapping if L(u + v) = Lu + Lv for all u, v V and L(λv) =
More informationFinite-Elements Method 2
Finite-Elements Method 2 January 29, 2014 2 From Applied Numerical Analysis Gerald-Wheatley (2004), Chapter 9. Finite-Elements Method 3 Introduction Finite-element methods (FEM) are based on some mathematical
More informationEngg. Math. I. Unit-I. Differential Calculus
Dr. Satish Shukla 1 of 50 Engg. Math. I Unit-I Differential Calculus Syllabus: Limits of functions, continuous functions, uniform continuity, monotone and inverse functions. Differentiable functions, Rolle
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 informationSummary of time independent electrodynamics
hapter 10 Summary of time independent electrodynamics 10.1 Electrostatics Physical law oulomb s law charges as origin of electric field Superposition principle ector of the electric field E(x) in vacuum
More informationPART IV Spectral Methods
PART IV Spectral Methods Additional References: R. Peyret, Spectral methods for incompressible viscous flow, Springer (2002), B. Mercier, An introduction to the numerical analysis of spectral methods,
More informationPDEs, part 1: Introduction and elliptic PDEs
PDEs, part 1: Introduction and elliptic PDEs Anna-Karin Tornberg Mathematical Models, Analysis and Simulation Fall semester, 2013 Partial di erential equations The solution depends on several variables,
More information2. Polynomial interpolation
2. Polynomial interpolation Contents 2. POLYNOMIAL INTERPOLATION... 1 2.1 TYPES OF INTERPOLATION... 1 2.2 LAGRANGE ONE-DIMENSIONAL INTERPOLATION... 2 2.3 NATURAL COORDINATES*... 15 2.4 HERMITE ONE-DIMENSIONAL
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 informationPreliminary Examination, Numerical Analysis, August 2016
Preliminary Examination, Numerical Analysis, August 2016 Instructions: This exam is closed books and notes. The time allowed is three hours and you need to work on any three out of questions 1-4 and any
More informationFinite Element Method-Part II Isoparametric FE Formulation and some numerical examples Lecture 29 Smart and Micro Systems
Finite Element Method-Part II Isoparametric FE Formulation and some numerical examples Lecture 29 Smart and Micro Systems Introduction Till now we dealt only with finite elements having straight edges.
More informationGetting started: CFD notation
PDE of p-th order Getting started: CFD notation f ( u,x, t, u x 1,..., u x n, u, 2 u x 1 x 2,..., p u p ) = 0 scalar unknowns u = u(x, t), x R n, t R, n = 1,2,3 vector unknowns v = v(x, t), v R m, m =
More informationImplementation of the dg method
Implementation of the dg method Matteo Cicuttin CERMICS - École Nationale des Ponts et Chaussées Marne-la-Vallée March 6, 2017 Outline Model problem, fix notation Representing polynomials Computing integrals
More information256 Summary. D n f(x j ) = f j+n f j n 2n x. j n=1. α m n = 2( 1) n (m!) 2 (m n)!(m + n)!. PPW = 2π k x 2 N + 1. i=0?d i,j. N/2} N + 1-dim.
56 Summary High order FD Finite-order finite differences: Points per Wavelength: Number of passes: D n f(x j ) = f j+n f j n n x df xj = m α m dx n D n f j j n= α m n = ( ) n (m!) (m n)!(m + n)!. PPW =
More informationMATH H53 : Final exam
MATH H53 : Final exam 11 May, 18 Name: You have 18 minutes to answer the questions. Use of calculators or any electronic items is not permitted. Answer the questions in the space provided. If you run out
More informationIntroduction. J.M. Burgers Center Graduate Course CFD I January Least-Squares Spectral Element Methods
Introduction In this workshop we will introduce you to the least-squares spectral element method. As you can see from the lecture notes, this method is a combination of the weak formulation derived from
More informationSummary for Vector Calculus and Complex Calculus (Math 321) By Lei Li
Summary for Vector alculus and omplex alculus (Math 321) By Lei Li 1 Vector alculus 1.1 Parametrization urves, surfaces, or volumes can be parametrized. Below, I ll talk about 3D case. Suppose we use e
More informationTMA4220 Numerical Solution of Partial Differential Equations Using Element Methods Høst 2014
Norwegian University of Science and Technology Department of Mathematical Sciences TMA422 Numerical Solution of Partial Differential Equations Using Element Methods Høst 214 Exercise set 2 1 Consider the
More informationMathematical Methods wk 2: Linear Operators
John Magorrian, magog@thphysoxacuk These are work-in-progress notes for the second-year course on mathematical methods The most up-to-date version is available from http://www-thphysphysicsoxacuk/people/johnmagorrian/mm
More informationLEAST-SQUARES FINITE ELEMENT MODELS
LEAST-SQUARES FINITE ELEMENT MODELS General idea of the least-squares formulation applied to an abstract boundary-value problem Works of our group Application to Poisson s equation Application to flows
More informationApplied Mathematics Masters Examination Fall 2016, August 18, 1 4 pm.
Applied Mathematics Masters Examination Fall 16, August 18, 1 4 pm. Each of the fifteen numbered questions is worth points. All questions will be graded, but your score for the examination will be the
More informationLemma 8: Suppose the N by N matrix A has the following block upper triangular form:
17 4 Determinants and the Inverse of a Square Matrix In this section, we are going to use our knowledge of determinants and their properties to derive an explicit formula for the inverse of a square matrix
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 information18.06 Problem Set 2 Solution
18.06 Problem Set 2 Solution Total: 100 points Section 2.5. Problem 24: Use Gauss-Jordan elimination on [U I] to find the upper triangular U 1 : 1 a b 1 0 UU 1 = I 0 1 c x 1 x 2 x 3 = 0 1 0. 0 0 1 0 0
More informationMatrices A brief introduction
Matrices A brief introduction Basilio Bona DAUIN Politecnico di Torino Semester 1, 2014-15 B. Bona (DAUIN) Matrices Semester 1, 2014-15 1 / 44 Definitions Definition A matrix is a set of N real or complex
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 informationCS137 Introduction to Scientific Computing Winter Quarter 2004 Solutions to Homework #3
CS137 Introduction to Scientific Computing Winter Quarter 2004 Solutions to Homework #3 Felix Kwok February 27, 2004 Written Problems 1. (Heath E3.10) Let B be an n n matrix, and assume that B is both
More informationLinear Algebra. Week 7
Linear Algebra. Week 7 Dr. Marco A Roque Sol 10 / 09 / 2018 If {v 1, v 2,, v n } is a basis for a vector space V, then any vector v V, has a unique representation v = x 1 v 1 + x 2 v 2 + + x n v n where
More information12. Thermal and scalar field analysis
2. THERMAL AND SCALAR FIELD ANALYSIS... 2. INTRODUCTION... 2.2 GENERAL FIELD PROBLEM... 2 2.3 COMMON FLUX COMPONENTS... 5 2. GALERKIN INTEGRAL FORM... 6 2.5 GALERKIN INTEGRAL FORM*... 7 2.6 ORTHOTROPIC
More informationTopic 2-2: Derivatives of Vector Functions. Textbook: Section 13.2, 13.4
Topic 2-2: Derivatives of Vector Functions Textbook: Section 13.2, 13.4 Warm-Up: Parametrization of Circles Each of the following vector functions describe the position of an object traveling around 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 information. D Matrix Calculus D 1
D Matrix Calculus D 1 Appendix D: MATRIX CALCULUS D 2 In this Appendix we collect some useful formulas of matrix calculus that often appear in finite element derivations D1 THE DERIVATIVES OF VECTOR FUNCTIONS
More informationLinear Equations and Matrix
1/60 Chia-Ping Chen Professor Department of Computer Science and Engineering National Sun Yat-sen University Linear Algebra Gaussian Elimination 2/60 Alpha Go Linear algebra begins with a system of linear
More informationMath 23b Practice Final Summer 2011
Math 2b Practice Final Summer 211 1. (1 points) Sketch or describe the region of integration for 1 x y and interchange the order to dy dx dz. f(x, y, z) dz dy dx Solution. 1 1 x z z f(x, y, z) dy dx dz
More informationIntroduction to Boundary Value Problems
Chapter 5 Introduction to Boundary Value Problems When we studied IVPs we saw that we were given the initial value of a function and a differential equation which governed its behavior for subsequent times.
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 information1 Linear transformations; the basics
Linear Algebra Fall 2013 Linear Transformations 1 Linear transformations; the basics Definition 1 Let V, W be vector spaces over the same field F. A linear transformation (also known as linear map, or
More informationMath 550, Exam 1. 2/10/11.
Math 55, Exam. //. Read problems carefully. Show all work. No notes, calculator, or text. The exam is approximately 5 percent of the total grade. There are points total. Partial credit may be given. Write
More informationElectromagnetism HW 1 math review
Electromagnetism HW math review Problems -5 due Mon 7th Sep, 6- due Mon 4th Sep Exercise. The Levi-Civita symbol, ɛ ijk, also known as the completely antisymmetric rank-3 tensor, has the following properties:
More informationNumerical Methods for Engineers, Second edition: Chapter 1 Errata
Numerical Methods for Engineers, Second edition: Chapter 1 Errata 1. p.2 first line, remove the Free Software Foundation at 2. p.2 sixth line of the first proper paragraph, fe95.res should be replaced
More informationChapter 2. General concepts. 2.1 The Navier-Stokes equations
Chapter 2 General concepts 2.1 The Navier-Stokes equations The Navier-Stokes equations model the fluid mechanics. This set of differential equations describes the motion of a fluid. In the present work
More informationThe Spectral-Element Method: Introduction
The Spectral-Element Method: Introduction Heiner Igel Department of Earth and Environmental Sciences Ludwig-Maximilians-University Munich Computational Seismology 1 / 59 Outline 1 Introduction 2 Lagrange
More informationMath 265H: Calculus III Practice Midterm II: Fall 2014
Name: Section #: Math 65H: alculus III Practice Midterm II: Fall 14 Instructions: This exam has 7 problems. The number of points awarded for each question is indicated in the problem. Answer each question
More information1. (a) (5 points) Find the unit tangent and unit normal vectors T and N to the curve. r (t) = 3 cos t, 0, 3 sin t, r ( 3π
1. a) 5 points) Find the unit tangent and unit normal vectors T and N to the curve at the point P 3, 3π, r t) 3 cos t, 4t, 3 sin t 3 ). b) 5 points) Find curvature of the curve at the point P. olution:
More informationCourse Summary Math 211
Course Summary Math 211 table of contents I. Functions of several variables. II. R n. III. Derivatives. IV. Taylor s Theorem. V. Differential Geometry. VI. Applications. 1. Best affine approximations.
More informationAlgebra Workshops 10 and 11
Algebra Workshops 1 and 11 Suggestion: For Workshop 1 please do questions 2,3 and 14. For the other questions, it s best to wait till the material is covered in lectures. Bilinear and Quadratic Forms on
More informationFast Numerical Methods for Stochastic Computations
Fast AreviewbyDongbinXiu May 16 th,2013 Outline Motivation 1 Motivation 2 3 4 5 Example: Burgers Equation Let us consider the Burger s equation: u t + uu x = νu xx, x [ 1, 1] u( 1) =1 u(1) = 1 Example:
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 informationME FINITE ELEMENT ANALYSIS FORMULAS
ME 2353 - FINITE ELEMENT ANALYSIS FORMULAS UNIT I FINITE ELEMENT FORMULATION OF BOUNDARY VALUE PROBLEMS 01. Global Equation for Force Vector, {F} = [K] {u} {F} = Global Force Vector [K] = Global Stiffness
More informationHIGHER-ORDER THEORIES
HIGHER-ORDER THEORIES Third-order Shear Deformation Plate Theory Displacement and strain fields Equations of motion Navier s solution for bending Layerwise Laminate Theory Interlaminar stress and strain
More informationMatrix Differentiation
Matrix Differentiation CS5240 Theoretical Foundations in Multimedia Leow Wee Kheng Department of Computer Science School of Computing National University of Singapore Leow Wee Kheng (NUS) Matrix Differentiation
More informationMathematical Concepts & Notation
Mathematical Concepts & Notation Appendix A: Notation x, δx: a small change in x t : the partial derivative with respect to t holding the other variables fixed d : the time derivative of a quantity that
More informationx + ye z2 + ze y2, y + xe z2 + ze x2, z and where T is the
1.(8pts) Find F ds where F = x + ye z + ze y, y + xe z + ze x, z and where T is the T surface in the pictures. (The two pictures are two views of the same surface.) The boundary of T is the unit circle
More informationL. Vandenberghe EE133A (Spring 2017) 3. Matrices. notation and terminology. matrix operations. linear and affine functions.
L Vandenberghe EE133A (Spring 2017) 3 Matrices notation and terminology matrix operations linear and affine functions complexity 3-1 Matrix a rectangular array of numbers, for example A = 0 1 23 01 13
More informationReview of Strain Energy Methods and Introduction to Stiffness Matrix Methods of Structural Analysis
uke University epartment of Civil and Environmental Engineering CEE 42L. Matrix Structural Analysis Henri P. Gavin Fall, 22 Review of Strain Energy Methods and Introduction to Stiffness Matrix Methods
More information