Method of Finite Elements I
|
|
- Leo Holt
- 5 years ago
- Views:
Transcription
1 Method of Finite Elements I PhD Candidate - Charilaos Mylonas HIL H33.1 and Boundary Conditions, 26 March, 2018 Institute of Structural Engineering Method of Finite Elements I 1
2 Outline 1 2 Penalty method Lagrange multiplier method Institute of Structural Engineering Method of Finite Elements I 2
3 Outline 1 2 Penalty method Lagrange multiplier method Institute of Structural Engineering Method of Finite Elements I 3
4 Integration Basic numerical integration Riemann sum Trapezoidal rule Comments: Both simple but inefficient! useful techniques in other contexts we use many more points than needed in order to achieve good accuracy Institute of Structural Engineering Method of Finite Elements I 4
5 Numerical example: Euler-Bernoulli beam We are going to apply the Galerkin method to the Euler Bernoulli beam, but with the basis N 1(x) = x 2 N 2(x) = x 3 N 3(x) = x 4 Instead of the Hermite basis functions, and discretize only displacements! Strong form: d 2 ( EI(x) d2 w ) = f(x) dx 2 dx 2 multiply with test function u(x) and integrate by parts twice a 0 EI(x) d2 w(x) d 2 u(x) dx = dx 2 dx 2 f(x)u(x)dx 0 a boundary terms [ ] 1 0 zero by basis selection Institute of Structural Engineering Method of Finite Elements I 5
6 Numerical integration - Galerkin method code The full commented code of the following demo will be made available to you in case you would like to play around. This code is a comparisson between analytical and numerical evaluation of integrals needed in the implementation of the Galerkin procedure. In the following slides the most important parts of the code are highlighted. Namely: The definition of the load the assembly of the system with analytical integration (symbolic integration) the assembly of the system with numerical integration (Gauss-Legendre) Institute of Structural Engineering Method of Finite Elements I 6
7 Numerical results 1 syms x % Coordinate a x i s 2 syms q ( x ) 3 % The l o a d i n g o f t h e beam : 4 q ( x ) = x. ˆ x. ˆ 4 5 NBASIS = 4 ; 6 f o r i = 1 : NBASIS 7 Phi ( i ) = x. ˆ ( i +1) ; 8 end 9 % Compute t h e i n t e g r a l s needed f o r t h e s t i f f n e s s m a t r i x : 10 f o r i =1:NBASIS 11 f o r j = 1 : NBASIS 12 K( i, j )= i n t ( EI d i f f ( Phi ( i ), 2 ) d i f f ( Phi ( j ), 2 ), 0, L ) ; 14 end 15 % A n a l y t i c a l computation o f t h e l o a d v e c t o r : 16 f ( i ) = i n t ( q ( x ) Phi ( i ), 0, L ) ; 17 end Institute of Structural Engineering Method of Finite Elements I 7
8 Numerical integration for LHS 1 q u a d r o r d e r = QUAD ORD; 2 3 %Loading the Gauss Legendre q ua dr at ur e nodes 4 % and w e i g h t s f o r a f i l e : 5 n o d e s w e i g h t s = W{ q u a d r o r d e r 1}; 6 nodes = ( n o d e s w e i g h t s. nodes +1) L / 2 ; 7 weights = W{ quadr order 1}. weights ; 8 9 f o r i =1: N b a s i s 10 f o r j =1: N b a s i s 11 V = e v a l ( s u b s ( EI d i f f ( Phi ( i ), 2 ) d i f f ( Phi ( j ), 2 ), nodes ) ) ; 12 K approx ( i, j ) = V weights ; 13 % Note t h a t we don t i n t e g r a t e a n a l y t i c a l l y! 14 % we can get EXACT r e s u l t s without doing so 15 end f a p p r o x ( i ) = e v a l ( s u b s ( Phi ( i ) q ( x ), nodes ) ) w e i g h t s ; 18 end Institute of Structural Engineering Method of Finite Elements I 8
9 Numerical Results Institute of Structural Engineering Method of Finite Elements I 9
10 Conclussions Notice the errors in the approximation (1.7e 13 for quadrature of order 4!). This is close to machine precission (the eps command in matlab). 1D order N means we evaluate the function inside the integral N times, and we sum the results according to some weight (see next slide). Instead of integrating analytically a function, we can get EXACT results by carefully selecting points of evaluation (instead of regular intervals as in trapezoidal and Riemann sums). The accuracy of the approximation depends on the Gaussian quadrature order Institute of Structural Engineering Method of Finite Elements I 10
11 Optimal points and weights integration for polynomials: - Optimal weights given the nodes For a polynomial f(x) = a 0 + a 1x + + a nx n, the integral is approximated as N f(x)dx w if(x i) 1 i=1 Assume x i are given. Also, for simplicity, f(x) = a 0 + a 1x + a 2x 2 + a 3x 3. Then, a 0 dx + a 1 xdx + a 2 x 2 dx + a 3 x 3 dx = w i(a 0 + a 1x i + a 2x 2 i + a3x3 i ) i=1 We equate all the terms, for every a i term. we arrive at the following linear system: w 1 1 x 1 x 2 x 3 x 4 dx w 2 x 2 1 x 2 2 x 2 x 2 1 = xdx 4 w 1 3 x 3 1 x 3 2 x 3 3 x w 4 x2 dx }{{} 1 x3 dx Vandermonde matrix Institute of Structural Engineering Method of Finite Elements I 11
12 Non-Gaussian Numerical computation of weights N-point quadrature, extact for polynomials up to degree N-1 Weights can be negative, in practice they have values of very different scales and solution of the system for high order quadrature is numerically unstable! The technique will approximate well integrals of functions that are approximated well with the corresponding polynomials. Lower order quadrature works, but as we will see it is far from the best we can do! (half as good). Institute of Structural Engineering Method of Finite Elements I 12
13 Legendre polynomials In the following, we will see how the Legendre polynomials are important for quadrature in 1D. We are to use x i such that this integral is exact for every polynomial f(x) with degree n < 2N 1. For a given N, solution x i are roots of the Legendre orthogonal polynomial of degree N Institute of Structural Engineering Method of Finite Elements I 13
14 Orthogonal polynomials Legendre orthogonal polynomials L nsatisfy (by construction) a : 2 L k, L m = L k (x)l m(x)dx = 1 2m + 1 δ mk From now on this integral is going to be refered to as L 2 inner product. Formulas for the first few polynomials are b, L 0(x) = 1 L 1(x) = x L 2(x) = 1 2 (3x2 1) L 3(x) = 1 8 (5x3 3x) L 4(x) = 1 8 (35x4 30x 2 + 3) a There are more general definitions of orthogonality and Gaussian quadrature where f, g w = Ω f(x)g(x)w(x)dx. The resulting polynomials are useful in other important contexts (stochastic finite elements!) b They are computed with a process known as Gram-Schmidt orthogonalization Institute of Structural Engineering Method of Finite Elements I 14
15 Legendre orthogonal polynomials and quadrature The roots of Legendre polynomials of degree N, used as integration nodes, yield the exact integral for f polynomial of maximum degree 2N 1. Outline of proof: 1 consider a polynomial f 2N 1 (x) = a 0 + a 1x + + a 2N 1 x 2N 1 2 Apply polynomial divission with Legendre polynomial of degree N: f 2N 1 (x) = g(x) }{{} degree N 1 L N (x) + r(x) }{{} degree N 1 3 {x (N) 1,, x (N) N } are the N roots of the Legendre polynomial of degree N. For brevity, we ommit the superscript x (N) i = x i. Chose the quadrature nodes to be these exact points: g(x)l n(x)dx N r(x)dx = 0 i=1 w ig(x i)l n(x i) + N w ir(x i) i=1 4 Since L n(x i) = 0, the first term of the RHS vanishes. Moreover, since the degree of g(x) is lower than N, the first integral term on the left is also zero! That is, through orthogonality and the abillity to express polynomials of order N 1 with a polynomial basis that contains all the monomials (such as the Legendre basis order N)! Institute of Structural Engineering Method of Finite Elements I 15
16 Considering again that we are evaluating the quadrature rule on the roots of Legendre polynomial of degree N, the term L N becomes zero Finally we are left with f 2N 1 (x i) = g(x i) 0 L N (x i) + r(x i) = r(x i) N f 2N 1 (x)dx = r(x) = w if 2N 1 (x i) 1 1 i=1 And we are left with the task of determining the weights as in the non-gaussian case. Therefore, surprizingly, a N-point gauss quadrature rule is exact for polynomials up to degree 2N 1. In practice, the computational determination of weights w i and nodes x i of 1D Gaussian quadrature rules happens through the so-called Golub-Welsch algorithm. Institute of Structural Engineering Method of Finite Elements I 16
17 Quadrilaterals Tensor-product Take tensor product of the 1D rules. Example: in x: in y rule for quadrilateral: w (ξ) i = {w (ξ) 1, w(ξ) 2, w(ξ) 3, w(ξ) 4 } ξ i = {ξ 1, ξ 2, ξ 3, ξ 4} w (η) i = {w (η) 1, w(η) 2, w(η) 3, x(η) 4 } η i = {η 1, η 2, η 3, η 4} 4 4 f(ξ, η)dxdy = w (η) i w (ξ) j f(ξ i, η j) 1 1 i=1 j=1 Institute of Structural Engineering Method of Finite Elements I 17
18 Outline Penalty method Lagrange multiplier method 1 2 Penalty method Lagrange multiplier method Institute of Structural Engineering Method of Finite Elements I 18
19 Homogeneous constraints Penalty method Lagrange multiplier method Linear system from FE discretization Consider a linear system from a FE discretization: K K 1N K N1... K NN [K]{u} = {f} u 1 =. u N f 1.. f N Homogeneous constraints the constraint u n = 0, in a matrix form consistent with our FE linear system reads, { } N v iu i = 0 i=0 u 1. u ṇ. u N = 0 Institute of Structural Engineering Method of Finite Elements I 19
20 Penalty method Lagrange multiplier method Imposing homogeneous constr. with the Penalty method Imposing the constraint The constraint vector should be in a form consistent with the matrix (a bilinear form that corresponds to a matrix). In order to impose the constraint, we simply multiply with a large number p >> tr([k]) and add that to our linear system. [V] = V ij = v iv j [V] = {v} T {v} [V] = This forces the degrees of freedom in the system to satisfy the constraint numerically (not explicitly). Institute of Structural Engineering Method of Finite Elements I 20
21 Adding to the linear system Penalty method Lagrange multiplier method Adding the constraint Simply set, [K] [K] pen with [K] pen = [K] + [V] p with p >> tr([k]) and solve {u} = [K] 1 pen {f} Inhomogeneous constraints and MFC: What if, we need to set a displacement to a given value (inhomogeneous constraint), or a displacement the same as another displacement (multi-dof constraints)? Institute of Structural Engineering Method of Finite Elements I 21
22 Multi- degree of freedom constraints Penalty method Lagrange multiplier method Equation for a multi-dof constraint Consider we have a constraint of the form 2u 5 + 3u 7 2u 10 = 0. Again, written as an equation, { } u 1. u 5 u 6 u 7 u 8 u 9 u 10. u N = 0 Institute of Structural Engineering Method of Finite Elements I 22
23 Penalty, multi-dof constraints Penalty method Lagrange multiplier method Application of the constraints Not much change when being systematic about the application! [V] ={v} T {v} [V] = The rest are the same! (multiply by a large number and add to [K] before solution). Institute of Structural Engineering Method of Finite Elements I 23
24 Penalty method Lagrange multiplier method Non-homogeneous constraints with the Penalty method and multiple constraints Non-zero RHS Say, we want to constrain u 3 2u 5 = 3 and at the same time u 2 u 3 = 0. In matrix form, the constraints read: { } } {{ } {v 1 } { } } {{ } {v 2 } u 1 u 2 u 3 u 4 u 5. u 1 u 2 u 3 u 4 u 5. = 3 }{{} a 1 = 0 }{{} a 2 Institute of Structural Engineering Method of Finite Elements I 24
25 Application of BC and solution Penalty method Lagrange multiplier method Application of BC on system In order to impose the constraint on the system, we proceed as such: Multiplying both equations with the corresponding {v i } T and summing up the constraints for the left-hand side we get, [V ] = 2 {v i } T {v i } i=1 and for the right-hand-side: 2 {A} = {v i } T a i i=1 [K pen] = [K] + [V ] p {f pen} = {f} + {A} p What changes, is that we also have a right-hand-side due to the non-homogeneous constraints! Solution Simply compute {u} = [K pen] 1 {f pen} Institute of Structural Engineering Method of Finite Elements I 25
26 Lagrange multipliers Penalty method Lagrange multiplier method constraints for LHS Another method, just as general and simple to implement, is the Lagrange multiplier method: given the constraint matrix [v] from the previous slide, set [ ] [K] [v] T [K Lagr ] = [v] 0 K K 1N v 11 v [K Lagr ] = K N1... K NN v 1N v 2N v v 1N 0 0 v v 2N 0 0 Constraints at RHS Instead of adding, we are again appending. {f Lagr } = { {f} {P } Where, {P } is a vector containing { } the values of the non-homogeneous constraints. From the 3 previous example, {P } =. 0 } Institute of Structural Engineering Method of Finite Elements I 26
27 Comparisson Penalty method Lagrange multiplier method The clear advantage of Lagr. Multiplier method over penalty is accuracy. Penalty has the advantage that does not introduce degrees of freedom to the linear system. However, this rarely becomes a real problem with modern hardware. Some linear system solvers rely on the matrix being positive definite for fast convergence - Lagrange multipliers may not work with some linear system solvers Lagrange multiplier method provides directly the constraint forces (as the values corresponding to the rows and columns of the constraints). When constraints are linearly dependent (a constraint vector can be expressed as linear combination of other constraints in the same system) iterative solvers (i.e. GMRES) have to be used. Linear independence of constraints is hard to judge in some cases (however, there are linear algebra techniques one can perform to remedy that). The multi-freedom homogeneous case for the Master-slave elimination method, is described partially in the project description. This method is also exact (as Lagr. Multipliers). Institute of Structural Engineering Method of Finite Elements I 27
LECTURE 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 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 informationMulti-Point Constraints
Multi-Point Constraints Multi-Point Constraints Multi-Point Constraints Single point constraint examples Multi-Point constraint examples linear, homogeneous linear, non-homogeneous linear, homogeneous
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 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 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 informationNumerical Analysis Preliminary Exam 10.00am 1.00pm, January 19, 2018
Numerical Analysis Preliminary Exam 0.00am.00pm, January 9, 208 Instructions. You have three hours to complete this exam. Submit solutions to four (and no more) of the following six problems. Please start
More informationMA2501 Numerical Methods Spring 2015
Norwegian University of Science and Technology Department of Mathematics MA5 Numerical Methods Spring 5 Solutions to exercise set 9 Find approximate values of the following integrals using the adaptive
More informationSection 6.6 Gaussian Quadrature
Section 6.6 Gaussian Quadrature Key Terms: Method of undetermined coefficients Nonlinear systems Gaussian quadrature Error Legendre polynomials Inner product Adapted from http://pathfinder.scar.utoronto.ca/~dyer/csca57/book_p/node44.html
More informationIntroduction to Finite Element Method
Introduction to Finite Element Method Dr. Rakesh K Kapania Aerospace and Ocean Engineering Department Virginia Polytechnic Institute and State University, Blacksburg, VA AOE 524, Vehicle Structures Summer,
More informationCOURSE Numerical integration of functions (continuation) 3.3. The Romberg s iterative generation method
COURSE 7 3. Numerical integration of functions (continuation) 3.3. The Romberg s iterative generation method The presence of derivatives in the remainder difficulties in applicability to practical problems
More informationNumerical Analysis Preliminary Exam 10 am to 1 pm, August 20, 2018
Numerical Analysis Preliminary Exam 1 am to 1 pm, August 2, 218 Instructions. You have three hours to complete this exam. Submit solutions to four (and no more) of the following six problems. Please start
More informationCS 450 Numerical Analysis. Chapter 8: Numerical Integration and Differentiation
Lecture slides based on the textbook Scientific Computing: An Introductory Survey by Michael T. Heath, copyright c 2018 by the Society for Industrial and Applied Mathematics. http://www.siam.org/books/cl80
More informationAPPENDIX B GRAM-SCHMIDT PROCEDURE OF ORTHOGONALIZATION. Let V be a finite dimensional inner product space spanned by basis vector functions
301 APPENDIX B GRAM-SCHMIDT PROCEDURE OF ORTHOGONALIZATION Let V be a finite dimensional inner product space spanned by basis vector functions {w 1, w 2,, w n }. According to the Gram-Schmidt Process an
More informationChapter 3 Variational Formulation & the Galerkin Method
Institute of Structural Engineering Page 1 Chapter 3 Variational Formulation & the Galerkin Method Institute of Structural Engineering Page 2 Today s Lecture Contents: Introduction Differential formulation
More informationMCE503: Modeling and Simulation of Mechatronic Systems. Modeling of Continuous Beams: Finite-Mode Approach
MCE503: Modeling and Simulation of Mechatronic Systems MCE503 p./23 Modeling of Continuous Beams: Finite-Mode Approach Reading: KMR Chapter 0 Cleveland State University Mechanical Engineering Hanz Richter,
More informationIn numerical analysis quadrature refers to the computation of definite integrals.
Numerical Quadrature In numerical analysis quadrature refers to the computation of definite integrals. f(x) a x i x i+1 x i+2 b x A traditional way to perform numerical integration is to take a piece of
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 informationx n+1 = x n f(x n) f (x n ), n 0.
1. Nonlinear Equations Given scalar equation, f(x) = 0, (a) Describe I) Newtons Method, II) Secant Method for approximating the solution. (b) State sufficient conditions for Newton and Secant to converge.
More information12.0 Properties of orthogonal polynomials
12.0 Properties of orthogonal polynomials In this section we study orthogonal polynomials to use them for the construction of quadrature formulas investigate projections on polynomial spaces and their
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 4: Interpolation and Approximation. October 28, 2005
Chapter 4: Interpolation and Approximation October 28, 2005 Outline 1 2.4 Linear Interpolation 2 4.1 Lagrange Interpolation 3 4.2 Newton Interpolation and Divided Differences 4 4.3 Interpolation Error
More 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 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 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 informationClass notes: Approximation
Class notes: Approximation Introduction Vector spaces, linear independence, subspace The goal of Numerical Analysis is to compute approximations We want to approximate eg numbers in R or C vectors in R
More informationFixed point iteration and root finding
Fixed point iteration and root finding The sign function is defined as x > 0 sign(x) = 0 x = 0 x < 0. It can be evaluated via an iteration which is useful for some problems. One such iteration is given
More 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 informationInterpolation. 1. Judd, K. Numerical Methods in Economics, Cambridge: MIT Press. Chapter
Key References: Interpolation 1. Judd, K. Numerical Methods in Economics, Cambridge: MIT Press. Chapter 6. 2. Press, W. et. al. Numerical Recipes in C, Cambridge: Cambridge University Press. Chapter 3
More 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 informationSolutions: Problem Set 3 Math 201B, Winter 2007
Solutions: Problem Set 3 Math 201B, Winter 2007 Problem 1. Prove that an infinite-dimensional Hilbert space is a separable metric space if and only if it has a countable orthonormal basis. Solution. If
More informationYou may not start to read the questions printed on the subsequent pages until instructed to do so by the Invigilator.
MATHEMATICAL TRIPOS Part IB Thursday 7 June 2007 9 to 12 PAPER 3 Before you begin read these instructions carefully. Each question in Section II carries twice the number of marks of each question in Section
More informationSOLVING LINEAR SYSTEMS
SOLVING LINEAR SYSTEMS We want to solve the linear system a, x + + a,n x n = b a n, x + + a n,n x n = b n This will be done by the method used in beginning algebra, by successively eliminating unknowns
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 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 informationAPPM/MATH Problem Set 6 Solutions
APPM/MATH 460 Problem Set 6 Solutions This assignment is due by 4pm on Wednesday, November 6th You may either turn it in to me in class or in the box outside my office door (ECOT ) Minimal credit will
More informationChapter 3 Numerical Methods
Chapter 3 Numerical Methods Part 2 3.2 Systems of Equations 3.3 Nonlinear and Constrained Optimization 1 Outline 3.2 Systems of Equations 3.3 Nonlinear and Constrained Optimization Summary 2 Outline 3.2
More informationWRT in 2D: Poisson Example
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
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 informationLecture 3: QR-Factorization
Lecture 3: QR-Factorization This lecture introduces the Gram Schmidt orthonormalization process and the associated QR-factorization of matrices It also outlines some applications of this factorization
More informationApplied Linear Algebra in Geoscience Using MATLAB
Applied Linear Algebra in Geoscience Using MATLAB Contents Getting Started Creating Arrays Mathematical Operations with Arrays Using Script Files and Managing Data Two-Dimensional Plots Programming in
More informationMA22S3 Summary Sheet: Ordinary Differential Equations
MA22S3 Summary Sheet: Ordinary Differential Equations December 14, 2017 Kreyszig s textbook is a suitable guide for this part of the module. Contents 1 Terminology 1 2 First order separable 2 2.1 Separable
More informationNumerical Methods - Numerical Linear Algebra
Numerical Methods - Numerical Linear Algebra Y. K. Goh Universiti Tunku Abdul Rahman 2013 Y. K. Goh (UTAR) Numerical Methods - Numerical Linear Algebra I 2013 1 / 62 Outline 1 Motivation 2 Solving Linear
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 informationRitz Method. Introductory Course on Multiphysics Modelling
Ritz Method Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Institute of Fundamental Technological Research of the Polish Academy of Sciences Warsaw Poland
More informationPhysics 331 Introduction to Numerical Techniques in Physics
Physics 331 Introduction to Numerical Techniques in Physics Instructor: Joaquín Drut Lecture 12 Last time: Polynomial interpolation: basics; Lagrange interpolation. Today: Quick review. Formal properties.
More informationOutline. 1 Interpolation. 2 Polynomial Interpolation. 3 Piecewise Polynomial Interpolation
Outline Interpolation 1 Interpolation 2 3 Michael T. Heath Scientific Computing 2 / 56 Interpolation Motivation Choosing Interpolant Existence and Uniqueness Basic interpolation problem: for given data
More informationLinear DifferentiaL Equation
Linear DifferentiaL Equation Massoud Malek The set F of all complex-valued functions is known to be a vector space of infinite dimension. Solutions to any linear differential equations, form a subspace
More informationAn Analysis of Five Numerical Methods for Approximating Certain Hypergeometric Functions in Domains within their Radii of Convergence
An Analysis of Five Numerical Methods for Approximating Certain Hypergeometric Functions in Domains within their Radii of Convergence John Pearson MSc Special Topic Abstract Numerical approximations of
More informationNumerical Analysis Solution of Algebraic Equation (non-linear equation) 1- Trial and Error. 2- Fixed point
Numerical Analysis Solution of Algebraic Equation (non-linear equation) 1- Trial and Error In this method we assume initial value of x, and substitute in the equation. Then modify x and continue till we
More informationCalculus and linear algebra for biomedical engineering Week 3: Matrices, linear systems of equations, and the Gauss algorithm
Calculus and linear algebra for biomedical engineering Week 3: Matrices, linear systems of equations, and the Gauss algorithm Hartmut Führ fuehr@matha.rwth-aachen.de Lehrstuhl A für Mathematik, RWTH Aachen
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 informationHigher-order ordinary differential equations
Higher-order ordinary differential equations 1 A linear ODE of general order n has the form a n (x) dn y dx n +a n 1(x) dn 1 y dx n 1 + +a 1(x) dy dx +a 0(x)y = f(x). If f(x) = 0 then the equation is called
More informationIntroduction to Finite Element Modelling in Geosciences
Introduction to Finite Element Modelling in Geosciences Dave A. May (dave.may@erdw.ethz.ch) Marcel Frehner (marcel.frehner@erdw.ethz.ch) Mike Afanasiev (ETH Zürich) Patrick Sanan (USI Lugano) 651-4144-L
More informationNumerical Modelling in Geosciences
Numerical Modelling in Geosciences Dave A. May (ETHZ) dave.may@erdw.ethz.ch Laetitia Le Pourhiet (UPMC) Jonas Ruh (UPMC) Liang Zheng (ETHZ) Zooming In between Plates (FP7 64713) ZIP WP5.1 Short-course
More informationExamination paper for TMA4215 Numerical Mathematics
Department of Mathematical Sciences Examination paper for TMA425 Numerical Mathematics Academic contact during examination: Trond Kvamsdal Phone: 93058702 Examination date: 6th of December 207 Examination
More informationyou expect to encounter difficulties when trying to solve A x = b? 4. A composite quadrature rule has error associated with it in the following form
Qualifying exam for numerical analysis (Spring 2017) Show your work for full credit. If you are unable to solve some part, attempt the subsequent parts. 1. Consider the following finite difference: f (0)
More informationHomework and Computer Problems for Math*2130 (W17).
Homework and Computer Problems for Math*2130 (W17). MARCUS R. GARVIE 1 December 21, 2016 1 Department of Mathematics & Statistics, University of Guelph NOTES: These questions are a bare minimum. You should
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 information4. Numerical Quadrature. Where analytical abilities end... continued
4. Numerical Quadrature Where analytical abilities end... continued Where analytical abilities end... continued, November 30, 22 1 4.3. Extrapolation Increasing the Order Using Linear Combinations Once
More informationOrthogonal Polynomials and Gaussian Quadrature
Orthogonal Polynomials and Gaussian Quadrature 1. Orthogonal polynomials Given a bounded, nonnegative, nondecreasing function w(x) on an interval, I of the real line, we consider the Hilbert space L 2
More informationHierarchical Parallel Solution of Stochastic Systems
Hierarchical Parallel Solution of Stochastic Systems Second M.I.T. Conference on Computational Fluid and Solid Mechanics Contents: Simple Model of Stochastic Flow Stochastic Galerkin Scheme Resulting Equations
More informationMATH 235: Inner Product Spaces, Assignment 7
MATH 235: Inner Product Spaces, Assignment 7 Hand in questions 3,4,5,6,9, by 9:3 am on Wednesday March 26, 28. Contents Orthogonal Basis for Inner Product Space 2 2 Inner-Product Function Space 2 3 Weighted
More informationHow might we evaluate this? Suppose that, by some good luck, we knew that. x 2 5. x 2 dx 5
8.4 1 8.4 Partial Fractions Consider the following integral. 13 2x (1) x 2 x 2 dx How might we evaluate this? Suppose that, by some good luck, we knew that 13 2x (2) x 2 x 2 = 3 x 2 5 x + 1 We could then
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 informationLecture Note 3: Interpolation and Polynomial Approximation. Xiaoqun Zhang Shanghai Jiao Tong University
Lecture Note 3: Interpolation and Polynomial Approximation Xiaoqun Zhang Shanghai Jiao Tong University Last updated: October 10, 2015 2 Contents 1.1 Introduction................................ 3 1.1.1
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 informationApplied Numerical Analysis Quiz #2
Applied Numerical Analysis Quiz #2 Modules 3 and 4 Name: Student number: DO NOT OPEN UNTIL ASKED Instructions: Make sure you have a machine-readable answer form. Write your name and student number in the
More informationSpectral Representation of Random Processes
Spectral Representation of Random Processes Example: Represent u(t,x,q) by! u K (t, x, Q) = u k (t, x) k(q) where k(q) are orthogonal polynomials. Single Random Variable:! Let k (Q) be orthogonal with
More informationAProofoftheStabilityoftheSpectral Difference Method For All Orders of Accuracy
AProofoftheStabilityoftheSpectral Difference Method For All Orders of Accuracy Antony Jameson 1 1 Thomas V. Jones Professor of Engineering Department of Aeronautics and Astronautics Stanford University
More information5. Hand in the entire exam booklet and your computer score sheet.
WINTER 2016 MATH*2130 Final Exam Last name: (PRINT) First name: Student #: Instructor: M. R. Garvie 19 April, 2016 INSTRUCTIONS: 1. This is a closed book examination, but a calculator is allowed. The test
More informationSeparation of Variables in Linear PDE: One-Dimensional Problems
Separation of Variables in Linear PDE: One-Dimensional Problems Now we apply the theory of Hilbert spaces to linear differential equations with partial derivatives (PDE). We start with a particular example,
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 informationNumerical integration and differentiation. Unit IV. Numerical Integration and Differentiation. Plan of attack. Numerical integration.
Unit IV Numerical Integration and Differentiation Numerical integration and differentiation quadrature classical formulas for equally spaced nodes improper integrals Gaussian quadrature and orthogonal
More informationScientific Computing: An Introductory Survey
Scientific Computing: An Introductory Survey Chapter 7 Interpolation Prof. Michael T. Heath Department of Computer Science University of Illinois at Urbana-Champaign Copyright c 2002. Reproduction permitted
More information10. Applications of 1-D Hermite elements
10. Applications of 1-D Hermite elements... 1 10.1 Introduction... 1 10.2 General case fourth-order beam equation... 3 10.3 Integral form... 5 10.4 Element Arrays... 7 10.5 C1 Element models... 8 10.6
More informationElements of linear algebra
Elements of linear algebra Elements of linear algebra A vector space S is a set (numbers, vectors, functions) which has addition and scalar multiplication defined, so that the linear combination c 1 v
More informationOutline. 1 Numerical Integration. 2 Numerical Differentiation. 3 Richardson Extrapolation
Outline Numerical Integration Numerical Differentiation Numerical Integration Numerical Differentiation 3 Michael T. Heath Scientific Computing / 6 Main Ideas Quadrature based on polynomial interpolation:
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 informationQR Decomposition. When solving an overdetermined system by projection (or a least squares solution) often the following method is used:
(In practice not Gram-Schmidt, but another process Householder Transformations are used.) QR Decomposition When solving an overdetermined system by projection (or a least squares solution) often the following
More informationMiddle East Technical University Department of Mechanical Engineering ME 413 Introduction to Finite Element Analysis. Chapter 2 Introduction to FEM
Middle East Technical University Department of Mechanical Engineering ME 43 Introction to Finite Element Analysis Chapter 2 Introction to FEM These notes are prepared by Dr. Cüneyt Sert http://www.me.metu.e.tr/people/cuneyt
More informationDepartment of Applied Mathematics and Theoretical Physics. AMA 204 Numerical analysis. Exam Winter 2004
Department of Applied Mathematics and Theoretical Physics AMA 204 Numerical analysis Exam Winter 2004 The best six answers will be credited All questions carry equal marks Answer all parts of each question
More informationStabilization and Acceleration of Algebraic Multigrid Method
Stabilization and Acceleration of Algebraic Multigrid Method Recursive Projection Algorithm A. Jemcov J.P. Maruszewski Fluent Inc. October 24, 2006 Outline 1 Need for Algorithm Stabilization and Acceleration
More informationSolving the steady state diffusion equation with uncertainty Final Presentation
Solving the steady state diffusion equation with uncertainty Final Presentation Virginia Forstall vhfors@gmail.com Advisor: Howard Elman elman@cs.umd.edu Department of Computer Science May 6, 2012 Problem
More informationLU Factorization. Marco Chiarandini. DM559 Linear and Integer Programming. Department of Mathematics & Computer Science University of Southern Denmark
DM559 Linear and Integer Programming LU Factorization Marco Chiarandini Department of Mathematics & Computer Science University of Southern Denmark [Based on slides by Lieven Vandenberghe, UCLA] Outline
More informationFinite Element Modeling and Analysis. CE 595: Course Part 2 Amit H. Varma
Finite Element Modeling and Analysis CE 595: Course Part 2 Amit H. Varma Discussion of planar elements Constant Strain Triangle (CST) - easiest and simplest finite element Displacement field in terms of
More informationCS 257: Numerical Methods
CS 57: Numerical Methods Final Exam Study Guide Version 1.00 Created by Charles Feng http://www.fenguin.net CS 57: Numerical Methods Final Exam Study Guide 1 Contents 1 Introductory Matter 3 1.1 Calculus
More informationFinite-dimensional spaces. C n is the space of n-tuples x = (x 1,..., x n ) of complex numbers. It is a Hilbert space with the inner product
Chapter 4 Hilbert Spaces 4.1 Inner Product Spaces Inner Product Space. A complex vector space E is called an inner product space (or a pre-hilbert space, or a unitary space) if there is a mapping (, )
More informationIntroduction to Finite Element Method. Dr. Aamer Haque
Introduction to Finite Element Method 4 th Order Beam Equation Dr. Aamer Haque http://math.iit.edu/~ahaque6 ahaque7@iit.edu Illinois Institute of Technology July 1, 009 Outline Euler-Bernoulli Beams Assumptions
More informationExercises * on Linear Algebra
Exercises * on Linear Algebra Laurenz Wiskott Institut für Neuroinformatik Ruhr-Universität Bochum, Germany, EU 4 February 7 Contents Vector spaces 4. Definition...............................................
More informationComputational Linear Algebra
Computational Linear Algebra PD Dr. rer. nat. habil. Ralf Peter Mundani Computation in Engineering / BGU Scientific Computing in Computer Science / INF Winter Term 2017/18 Part 2: Direct Methods PD Dr.
More informationRomberg Integration and Gaussian Quadrature
Romberg Integration and Gaussian Quadrature P. Sam Johnson October 17, 014 P. Sam Johnson (NITK) Romberg Integration and Gaussian Quadrature October 17, 014 1 / 19 Overview We discuss two methods for integration.
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 informationPreliminary Examination in Numerical Analysis
Department of Applied Mathematics Preliminary Examination in Numerical Analysis August 7, 06, 0 am pm. Submit solutions to four (and no more) of the following six problems. Show all your work, and justify
More 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 informationChapter 5 Structural Elements: The truss & beam elements
Institute of Structural Engineering Page 1 Chapter 5 Structural Elements: The truss & beam elements Institute of Structural Engineering Page 2 Chapter Goals Learn how to formulate the Finite Element Equations
More information4 Finite Element Method for Trusses
4 Finite Element Method for Trusses To solve the system of linear equations that arises in IPM, it is necessary to assemble the geometric matrix B a. For the sake of simplicity, the applied force vector
More informationConcepts. 3.1 Numerical Analysis. Chapter Numerical Analysis Scheme
Chapter 3 Concepts The objective of this work is to create a framework to implement multi-disciplinary finite element applications. Before starting, it is necessary to explain some basic concepts of the
More informationQuintic beam closed form matrices (revised 2/21, 2/23/12) General elastic beam with an elastic foundation
General elastic beam with an elastic foundation Figure 1 shows a beam-column on an elastic foundation. The beam is connected to a continuous series of foundation springs. The other end of the foundation
More informationThis ODE arises in many physical systems that we shall investigate. + ( + 1)u = 0. (λ + s)x λ + s + ( + 1) a λ. (s + 1)(s + 2) a 0
Legendre equation This ODE arises in many physical systems that we shall investigate We choose We then have Substitution gives ( x 2 ) d 2 u du 2x 2 dx dx + ( + )u u x s a λ x λ a du dx λ a λ (λ + s)x
More information