The use of exact values at quadrature points in the boundary element method
|
|
- Lucas Powers
- 5 years ago
- Views:
Transcription
1 The use of exact values at quadrature points in the boundary element method N.D. Stringfellow, R.N.L. Smith Applied Mathematics & OR Group, Cranfield University, Shrivenham, Swindon, U.K., SN6 SLA N. D. String cranfield. ac. uk Abstract This paper describes the use of exact values in the solution of boundary element problems and compares the accuracy of the corresponding numerical solutions with those obtained using the traditional interpolation method. The new scheme improves accuracy significantly, perhaps by several orders of magnitude, at little computational cost. 1 Introduction The traditional approach to solving a boundary element problem over a domain il with boundary conditions $ defined over a boundary F has been to take the geometry of F, discretise this geometry and use a set of polynomial interpolation functions to represent each section of the discretised shape. The assembly of entries in the coefficient matrices then consists of using these interpolants to calculate the x and y coordinates (in 2 dimensions) of an approximation to the boundary F together with information related to the boundary conditions and the particular type of problem (potential, elasticity, etc.) which is to be solved. A similar interpolation scheme is used to approximate the values of the boundary conditions along the discretised region.
2 240 Boundary Elements Figure 1: Example domain for interior problem For many real problems, the geometry may be described exactly in terms of simple shapes such as straight lines, circular arcs and splines. The boundary conditions also frequently consist of a combination of simple functions such as trigonometric functions, polynomials and logarithmic or exponential functions. Using these functions it is possible to calculate the exact values of the x and y coordinates of the geometry and the exact values of the boundary conditions at each quadrature point along an element, and to use these exact values in the calculation of the coefficient matrices. The use of these values immediately remove much of the approximation which is usually carried out before the problem is solved numerically. We first consider the collocation formulation for the 2-dimensional potential boundary element method using polynomial interpolation functions to represent the geometry and boundary conditions. We show how removing interpolation functions from the formulation removes inaccuracies and present numerical results to demonstrate improved convergence. 2 Polynomial Collocation Method Here we consider the traditional collocation approach to the boundary element method for potential problems over a domain fi. For interior problems an example domain is shown in figure 1. The general form of the boundary integral equation for potential problems is:- cu+ / up*dr= I u*pdt (1) Jr Jr
3 Boundary Elements 241 Figure 2: Discretised boundary where:- c a constant dependent upon the boundary shape u the function values p the derivative values n* and p* fundamental solutions of the PDE F the boundary The boundary is chopped up into "elements" as shown in figure 2 to give the new equation c%+ ]T / up*df= ^2 I u*pdr (2) At this point we would normally assume that u and p may be represented by a polynomial of order n varying over the element F^ and that the geometrical shape of the element can be represented by a polynomial of order ra varying over the element, with ra n for isoparametric elements. Introducing a parameterisation of each element F^ using a parameter t and representing the boundary condition shape functions as </>(t) and the geometry shape functions as ip(t) we obtain a new equation with p* and u* being functions of i >i(t) and J(t) the Jacobian of :,he transformation. Note that equation 3 is not the same as equation 2 but is an approximation to it.
4 242 Boundary Elements For potential problems we have the kernel functions p* and u* given by p*, u* = log r (4) where r is the vector from the collocation point at which u is the potential, to the point on the boundary of I\, and n is the normal vector to the boundary on IV However both r and n are now only approximated in the terms J?* and u* because the geometry is approximated at these points using the geometry shape functions if)i(t) and the Jacobian J(t) evaluated using the i/ji(t) gives an approximation J(t)dt to dt. Our equation may now be written in which we still have to evaluate the integrals. In order to evaluate the integrals we use quadrature schemes - Gauss quadrature for ordinary non-singular integrals and schemes such as those devised by Smith [1] for the singular u* integrals - and the row-sum property ^y r*dr = o (6) for the p* integrals. After evaluating the integrals we are left with a system of equations Hu = Gt. (7) Collecting the known values on the right hand side and multiplying through we are left with a linear system of equations Ax = b (8) from which we may retrieve x using a linear solver such as Gaussian elimination. 3 Exact Geometry For most practical problems the geometrical composition of the boundary F is defined in terms of simple line sections such as straight lines, circular and elliptical sections. If the boundary is composed of straight lines any interpolant involving at least two points can produce any point x(to) along an element for a given to. A polynomial interpolant of a circular section x(to) will not lie on the circle except at the sampling points. This inaccuracy gives rise to numerical errors in the integrations. As an example, consider a quadratic element passing through the points ( 1,2 v^) (0,0) (1,2 - A/3) which can be used to interpolate a circle
5 Boundary Elements 243 Figure 3: Interpolant and circle centre (0,2) with radius 2 apparently quite well as shown in figure 3. Closer inspection of the region x G [-0.6, 0.4] as shown in figure 4 reveals that the interpolant is some distance from the "true" boundary. The actual values of the x and y coordinates used in the numerical integrations of the numerical scheme already described are interpolated values at the quadrature points using the polynomial interpolation functions ^ (t). For simple geometries (e.g. circular sections) we know the true x and y coordinates at the quadrature points, and by using these values we remove the interpolation error from the kernel or Green's functions and from the Jacobian. Equation 5 may now be represented as cu + E E (Y ^ re r*=i ^Te = E E (9) where we now use p* and u* instead of p* and u* since we are modelling the exact geometry up to numerical accuracy. Note also that the Jacobian J is not now a function of interpolation functions. 4 Exact Boundary Conditions The function values used in the numerical scheme are also interpolated values at the quadrature points using the shape functions (fri(t). However, for most problems of interest we have analytical expressions for the boundary conditions which can be written as some function of the coordinates of the boundary points. We may therefore use the exact values of the boundary conditions at the quadrature points in order to form the 6 vector in the linear system which we wish to solve. Note that the exact values for the geometry
6 244 Boundary Elements Figure 4: Close up of interpolant and circle and boundary condition functions may be given at the quadrature points used to evaluate the u* singularities as well as for the Gauss quadrature points. Our code allows the boundary conditions to be denned as symbolic expressions consisting of common functions such as sin(x) or 32 on each section of the boundary. 5 Results In order to analyse the accuracy of the numerical schemes we use two exterior Neumann potential problems. The use of exterior problems ensures that imposing pure Neumann boundary conditions does not result in an underdetermined system. The use of Neumann boundary conditions avoids the necessity to use the row-sum singularity in the right hand side vector b where boundary condition information along the element is reduced to a single boundary condition at the singular node. The particular problems to be considered are those analysed by Kirkup and Henwood [2] which provide varying boundary conditions or geometries that use functions other than polynomials, enabling an effective comparison of the methods. For the two potential problems for which we have analytical expressions for the function values, we use as the error estimator the mean residual error (MRE), defined by 1 " %^# (10) for TV collocation points, with HJ the computed value, and Uj the exact value at the jth collocation point.
7 Boundary Elements 245 Figure 5: Square Example Multiple runs were performed, using isoparametric and subparametric elements, with the geometry varying from linear to cubic in the interpolant, and the function values ranging from constant to the same order as the geometry shape. Thefirstexample is a potential problem with Neumann boundary condition *± = * (11) dm?r(d2+4) ^ / imposed, where d is the distance from the centre of the side as shown in figure 5. The second example is a potential problem with Neumann boundary condition du 4 - cos(0) (12) dn 27r(17-8cos(<9)) imposed, where 9 is as shown infigure6. In order to minimise quadrature errors the number of Gauss points used was increased until the solution appeared to be stable. Sixteen points were used for the solutions given here. Although results were produced for different configurations of boundary shape functions, geometry shape functions and numbers of nodes, the general trend can be identified by consideration of just one combination of geometry and boundary condition shape functions. In order to provide a good comparison with standard methods, quadratic isoparametric elements were chosen as the polynomial interpolant elements. Since the mean relative error diminishes rapidly logarithmic scales (base 10) were used when plotting. For thefirstexample shown infigure5 the quadratic polynomial is able to interpolate the linear geometry exactly so that the polynomial interpolant and exact geometry methods are the same. As can be seen infigure7 the
8 246 Boundary Elements Figure 6: Circle Example interpolant/exact geom i exact geom & b.c.s x e-05 1e-06 1e Number of Nodes Figure 7: Square convergence (log^q scale)
9 Boundary Elements interpolant Hexact geom x- exact geom & b.c.s - -* 1e-06 1e-08 1e-10 1e-12 1e-14 1e Number of Nodes Figure 8: Circle convergence scale) case using the exact boundary condition values at the quadrature points has a consistently lower mean relative error than the other two cases with the same accuracy being achieved with just over half the number of elements used in the polynomial interpolant case. For the second example shown in figure 6 neither the geometry or boundary condition function can be exactly represented using polynomial interpolants and therefore all three methods produce different results. Figure 8 shows that using the exact geometry in this case produces better results than using a polynomial interpolant as shown by Kirkup & Henwood and our results for the use of exact geometry agree closely with theirs. However, the best results are produced by using the exact values of geometry and boundary condition function at the quadrature points. Here the mean relative error (and other error measures used) soon approach machine accuracy, in this case a double precision implementation with 16 digits of accuracy, a dramatic improvement in accuracy over the conventional BEM implementation. 6 Conclusions Good convergence can normally be achieved using polynomial interpolants for the geometry of the boundary and the boundary condition function. Using exact values of x and y coordinates at the quadrature points produces increased accuracy when the boundary cannot be modelled exactly using a
10 24 8 Boundary Elements polynomial interpolant. Where the boundary condition is a known function which cannot be modelled accurately using a polynomial interpolant, the use of exact boundary condition values at the quadrature points produces smaller errors and perhaps dramatically smaller errors. The use of exact values of x and y coordinates together with exact values of the boundary condition function at quadrature points can lead to a reduction in the number of elements required for a given error tolerance. As the matrix assembly process is O(n^) and the linear solver is O(n?), using exact values at quadrature points can lead to a large reduction in solution time for little extra initial computation. There would seem to be a clear need for more research in this area and we intend to consider a wider range of examples if different boundary conditions and also to consider the effect of exact boundary conditions on vector problems. References [1] Smith R.N.L. Specialised Quadrature Schemes for Boundary Integral Methods. First UK Conference on Boundary Integral Methods, pages , [2] Kirkup S.M., Henwood D.J. An Empirical Error Analysis of the Boundary Element Method Applied to Laplace's Equation. Applied Mathematical Modelling, 18:62-38, 1994.
Transactions on Modelling and Simulation vol 12, 1996 WIT Press, ISSN X
Simplifying integration for logarithmic singularities R.N.L. Smith Department ofapplied Mathematics & OR, Cranfield University, RMCS, Shrivenham, Swindon, Wiltshire SN6 SLA, UK Introduction Any implementation
More informationn 1 f n 1 c 1 n+1 = c 1 n $ c 1 n 1. After taking logs, this becomes
Root finding: 1 a The points {x n+1, }, {x n, f n }, {x n 1, f n 1 } should be co-linear Say they lie on the line x + y = This gives the relations x n+1 + = x n +f n = x n 1 +f n 1 = Eliminating α and
More informationCubic Splines MATH 375. J. Robert Buchanan. Fall Department of Mathematics. J. Robert Buchanan Cubic Splines
Cubic Splines MATH 375 J. Robert Buchanan Department of Mathematics Fall 2006 Introduction Given data {(x 0, f(x 0 )), (x 1, f(x 1 )),...,(x n, f(x n ))} which we wish to interpolate using a polynomial...
More informationx x2 2 + x3 3 x4 3. Use the divided-difference method to find a polynomial of least degree that fits the values shown: (b)
Numerical Methods - PROBLEMS. The Taylor series, about the origin, for log( + x) is x x2 2 + x3 3 x4 4 + Find an upper bound on the magnitude of the truncation error on the interval x.5 when log( + x)
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 informationFinite Mathematics Chapter 2. where a, b, c, d, h, and k are real numbers and neither a and b nor c and d are both zero.
Finite Mathematics Chapter 2 Section 2.1 Systems of Linear Equations: An Introduction Systems of Equations Recall that a system of two linear equations in two variables may be written in the general form
More informationPARTIAL DIFFERENTIAL EQUATIONS
MATHEMATICAL METHODS PARTIAL DIFFERENTIAL EQUATIONS I YEAR B.Tech By Mr. Y. Prabhaker Reddy Asst. Professor of Mathematics Guru Nanak Engineering College Ibrahimpatnam, Hyderabad. SYLLABUS OF MATHEMATICAL
More informationChapter 5 HIGH ACCURACY CUBIC SPLINE APPROXIMATION FOR TWO DIMENSIONAL QUASI-LINEAR ELLIPTIC BOUNDARY VALUE PROBLEMS
Chapter 5 HIGH ACCURACY CUBIC SPLINE APPROXIMATION FOR TWO DIMENSIONAL QUASI-LINEAR ELLIPTIC BOUNDARY VALUE PROBLEMS 5.1 Introduction When a physical system depends on more than one variable a general
More informationCHAPTER 4 NUMERICAL INTEGRATION METHODS TO EVALUATE TRIPLE INTEGRALS USING GENERALIZED GAUSSIAN QUADRATURE
CHAPTER 4 NUMERICAL INTEGRATION METHODS TO EVALUATE TRIPLE INTEGRALS USING GENERALIZED GAUSSIAN QUADRATURE 4.1 Introduction Many applications in science and engineering require the solution of three dimensional
More informationJEPPIAAR ENGINEERING COLLEGE
JEPPIAAR ENGINEERING COLLEGE Jeppiaar Nagar, Rajiv Gandhi Salai 600 119 DEPARTMENT OFMECHANICAL ENGINEERING QUESTION BANK VI SEMESTER ME6603 FINITE ELEMENT ANALYSIS Regulation 013 SUBJECT YEAR /SEM: III
More informationAlgebraic. techniques1
techniques Algebraic An electrician, a bank worker, a plumber and so on all have tools of their trade. Without these tools, and a good working knowledge of how to use them, it would be impossible for them
More informationTABLE OF CONTENTS INTRODUCTION, APPROXIMATION & ERRORS 1. Chapter Introduction to numerical methods 1 Multiple-choice test 7 Problem set 9
TABLE OF CONTENTS INTRODUCTION, APPROXIMATION & ERRORS 1 Chapter 01.01 Introduction to numerical methods 1 Multiple-choice test 7 Problem set 9 Chapter 01.02 Measuring errors 11 True error 11 Relative
More informationExact and Approximate Numbers:
Eact and Approimate Numbers: The numbers that arise in technical applications are better described as eact numbers because there is not the sort of uncertainty in their values that was described above.
More informationNumerical Analysis of Electromagnetic Fields
Pei-bai Zhou Numerical Analysis of Electromagnetic Fields With 157 Figures Springer-Verlag Berlin Heidelberg New York London Paris Tokyo Hong Kong Barcelona Budapest Contents Part 1 Universal Concepts
More informationTHE BOUNDARY ELEMENT METHOD
THE BOUNDARY ELEMENT METHOD SOLID MECHANICS AND ITS APPLICATIONS Volume 27 Series Editor: G.M.L. GLADWELL Solid Mechanics Division, Faculty of Engineering University of Waterloo Waterloo, Ontario, Canada
More informationMMJ1153 COMPUTATIONAL METHOD IN SOLID MECHANICS PRELIMINARIES TO FEM
B Course Content: A INTRODUCTION AND OVERVIEW Numerical method and Computer-Aided Engineering; Phsical problems; Mathematical models; Finite element method;. B Elements and nodes, natural coordinates,
More informationDepartment of Aerospace Engineering AE602 Mathematics for Aerospace Engineers Assignment No. 4
Department of Aerospace Engineering AE6 Mathematics for Aerospace Engineers Assignment No.. Decide whether or not the following vectors are linearly independent, by solving c v + c v + c 3 v 3 + c v :
More informationJim Lambers MAT 419/519 Summer Session Lecture 13 Notes
Jim Lambers MAT 419/519 Summer Session 2011-12 Lecture 13 Notes These notes correspond to Section 4.1 in the text. Least Squares Fit One of the most fundamental problems in science and engineering is data
More informationc 2007 Society for Industrial and Applied Mathematics
SIAM J. SCI. COMPUT. Vol. 9, No. 3, pp. 07 6 c 007 Society for Industrial and Applied Mathematics NUMERICAL QUADRATURES FOR SINGULAR AND HYPERSINGULAR INTEGRALS IN BOUNDARY ELEMENT METHODS MICHAEL CARLEY
More informationMATHEMATICAL METHODS INTERPOLATION
MATHEMATICAL METHODS INTERPOLATION I YEAR BTech By Mr Y Prabhaker Reddy Asst Professor of Mathematics Guru Nanak Engineering College Ibrahimpatnam, Hyderabad SYLLABUS OF MATHEMATICAL METHODS (as per JNTU
More informationCOMPUTATIONAL ELASTICITY
COMPUTATIONAL ELASTICITY Theory of Elasticity and Finite and Boundary Element Methods Mohammed Ameen Alpha Science International Ltd. Harrow, U.K. Contents Preface Notation vii xi PART A: THEORETICAL ELASTICITY
More informationMatrices. A matrix is a method of writing a set of numbers using rows and columns. Cells in a matrix can be referenced in the form.
Matrices A matrix is a method of writing a set of numbers using rows and columns. 1 2 3 4 3 2 1 5 7 2 5 4 2 0 5 10 12 8 4 9 25 30 1 1 Reading Information from a Matrix Cells in a matrix can be referenced
More informationChapter 3. Formulation of FEM for Two-Dimensional Problems
Chapter Formulation of FEM for Two-Dimensional Problems.1 Two-Dimensional FEM Formulation Many details of 1D and 2D formulations are the same. To demonstrate how a 2D formulation works we ll use the following
More informationThe detection of subsurface inclusions using internal measurements and genetic algorithms
The detection of subsurface inclusions using internal measurements and genetic algorithms N. S. Meral, L. Elliott2 & D, B, Ingham2 Centre for Computational Fluid Dynamics, Energy and Resources Research
More informationDifferential equations
Differential equations Math 7 Spring Practice problems for April Exam Problem Use the method of elimination to find the x-component of the general solution of x y = 6x 9x + y = x 6y 9y Soln: The system
More informationCHAPTER II MATHEMATICAL BACKGROUND OF THE BOUNDARY ELEMENT METHOD
CHAPTER II MATHEMATICAL BACKGROUND OF THE BOUNDARY ELEMENT METHOD For the second chapter in the thesis, we start with surveying the mathematical background which is used directly in the Boundary Element
More informationThomas Algorithm for Tridiagonal Matrix
P a g e 1 Thomas Algorithm for Tridiagonal Matrix Special Matrices Some matrices have a particular structure that can be exploited to develop efficient solution schemes. Two of those such systems are banded
More information7.6 The Inverse of a Square Matrix
7.6 The Inverse of a Square Matrix Copyright Cengage Learning. All rights reserved. What You Should Learn Verify that two matrices are inverses of each other. Use Gauss-Jordan elimination to find inverses
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 informationLectures 9-10: Polynomial and piecewise polynomial interpolation
Lectures 9-1: Polynomial and piecewise polynomial interpolation Let f be a function, which is only known at the nodes x 1, x,, x n, ie, all we know about the function f are its values y j = f(x j ), j
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 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 informationData Fitting and Uncertainty
TiloStrutz Data Fitting and Uncertainty A practical introduction to weighted least squares and beyond With 124 figures, 23 tables and 71 test questions and examples VIEWEG+ TEUBNER IX Contents I Framework
More informationDiscrete Projection Methods for Integral Equations
SUB Gttttingen 7 208 427 244 98 A 5141 Discrete Projection Methods for Integral Equations M.A. Golberg & C.S. Chen TM Computational Mechanics Publications Southampton UK and Boston USA Contents Sources
More informationNumerical integration in the axisymmetric finite element formulation
Advances in Engineering Software 31 (2000) 137 141 Short Communication Numerical integration in the axisymmetric finite element formulation J.D. Clayton a, J.J. Rencis b, * a The George W. Woodruff School
More informationxvi xxiii xxvi Construction of the Real Line 2 Is Every Real Number Rational? 3 Problems Algebra of the Real Numbers 7
About the Author v Preface to the Instructor xvi WileyPLUS xxii Acknowledgments xxiii Preface to the Student xxvi 1 The Real Numbers 1 1.1 The Real Line 2 Construction of the Real Line 2 Is Every Real
More informationUniversity of Hertfordshire Department of Mathematics. Study on the Dual Reciprocity Boundary Element Method
University of Hertfordshire Department of Mathematics Study on the Dual Reciprocity Boundary Element Method Wattana Toutip Technical Report 3 July 999 Preface The boundary Element method (BEM) is now recognised
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 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 informationCandidates are expected to have available a calculator. Only division by (x + a) or (x a) will be required.
Revision Checklist Unit C2: Core Mathematics 2 Unit description Algebra and functions; coordinate geometry in the (x, y) plane; sequences and series; trigonometry; exponentials and logarithms; differentiation;
More information8.8. Applications of Taylor Polynomials. Infinite Sequences and Series 8
8.8 Applications of Taylor Polynomials Infinite Sequences and Series 8 Applications of Taylor Polynomials In this section we explore two types of applications of Taylor polynomials. First we look at how
More informationResearch Statement. James Bremer Department of Mathematics, University of California, Davis
Research Statement James Bremer Department of Mathematics, University of California, Davis Email: bremer@math.ucdavis.edu Webpage: https.math.ucdavis.edu/ bremer I work in the field of numerical analysis,
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 informationPROJECTION METHODS FOR DYNAMIC MODELS
PROJECTION METHODS FOR DYNAMIC MODELS Kenneth L. Judd Hoover Institution and NBER June 28, 2006 Functional Problems Many problems involve solving for some unknown function Dynamic programming Consumption
More information4 The matrix approach
4 The matrix approach 4. Introduction It is not a coincidence that the finite element method blossomed at the time when computers were being developed, because it calls for a vast amount of computation
More information7.5 Operations with Matrices. Copyright Cengage Learning. All rights reserved.
7.5 Operations with Matrices Copyright Cengage Learning. All rights reserved. What You Should Learn Decide whether two matrices are equal. Add and subtract matrices and multiply matrices by scalars. Multiply
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 informationDM559 Linear and Integer Programming. Lecture 2 Systems of Linear Equations. Marco Chiarandini
DM559 Linear and Integer Programming Lecture Marco Chiarandini Department of Mathematics & Computer Science University of Southern Denmark Outline 1. Outline 1. 3 A Motivating Example You are organizing
More informationRegion 16 Board of Education. Precalculus Curriculum
Region 16 Board of Education Precalculus Curriculum 2008 1 Course Description This course offers students an opportunity to explore a variety of concepts designed to prepare them to go on to study calculus.
More informationAccelerated Advanced Algebra. Chapter 1 Patterns and Recursion Homework List and Objectives
Chapter 1 Patterns and Recursion Use recursive formulas for generating arithmetic, geometric, and shifted geometric sequences and be able to identify each type from their equations and graphs Write and
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 informationNumerical Analysis: Solutions of System of. Linear Equation. Natasha S. Sharma, PhD
Mathematical Question we are interested in answering numerically How to solve the following linear system for x Ax = b? where A is an n n invertible matrix and b is vector of length n. Notation: x denote
More informationNumerical Methods. Equations and Partial Fractions. Jaesung Lee
Numerical Methods Equations and Partial Fractions Jaesung Lee Solving linear equations Solving linear equations Introduction Many problems in engineering reduce to the solution of an equation or a set
More informationQuaternion Cubic Spline
Quaternion Cubic Spline James McEnnan jmcennan@mailaps.org May 28, 23 1. INTRODUCTION A quaternion spline is an interpolation which matches quaternion values at specified times such that the quaternion
More informationAlgebra I. CORE TOPICS (Key Concepts & Real World Context) UNIT TITLE
CHAPTER 1 Using variables Exponents and Order of Operations Exploring real numbers Adding real numbers Subtracting real numbers Multiplying and dividing real numbers The Distributive Property Properties
More informationComputational Methods. Least Squares Approximation/Optimization
Computational Methods Least Squares Approximation/Optimization Manfred Huber 2011 1 Least Squares Least squares methods are aimed at finding approximate solutions when no precise solution exists Find the
More informationFitting a Natural Spline to Samples of the Form (t, f(t))
Fitting a Natural Spline to Samples of the Form (t, f(t)) David Eberly, Geometric Tools, Redmond WA 9852 https://wwwgeometrictoolscom/ This work is licensed under the Creative Commons Attribution 4 International
More information10.2 ITERATIVE METHODS FOR SOLVING LINEAR SYSTEMS. The Jacobi Method
54 CHAPTER 10 NUMERICAL METHODS 10. ITERATIVE METHODS FOR SOLVING LINEAR SYSTEMS As a numerical technique, Gaussian elimination is rather unusual because it is direct. That is, a solution is obtained after
More informationCLASS NOTES Computational Methods for Engineering Applications I Spring 2015
CLASS NOTES Computational Methods for Engineering Applications I Spring 2015 Petros Koumoutsakos Gerardo Tauriello (Last update: July 27, 2015) IMPORTANT DISCLAIMERS 1. REFERENCES: Much of the material
More informationAPPENDIX A. Background Mathematics. A.1 Linear Algebra. Vector algebra. Let x denote the n-dimensional column vector with components x 1 x 2.
APPENDIX A Background Mathematics A. Linear Algebra A.. Vector algebra Let x denote the n-dimensional column vector with components 0 x x 2 B C @. A x n Definition 6 (scalar product). The scalar product
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 informationMEAN VALUE THEOREMS FUNCTIONS OF SINGLE & SEVERAL VARIABLES
MATHEMATICS-I MEAN VALUE THEOREMS FUNCTIONS OF SINGLE & SEVERAL VARIABLES I YEAR B.TECH By Y. Prabhaker Reddy Asst. Professor of Mathematics Guru Nanak Engineering College Ibrahimpatnam, Hyderabad. Name
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 informationCAM Ph.D. Qualifying Exam in Numerical Analysis CONTENTS
CAM Ph.D. Qualifying Exam in Numerical Analysis CONTENTS Preliminaries Round-off errors and computer arithmetic, algorithms and convergence Solutions of Equations in One Variable Bisection method, fixed-point
More informationSolving the complete-electrode direct model of ERT using the boundary element method and the method of fundamental solutions
Solving the complete-electrode direct model of ERT using the boundary element method and the method of fundamental solutions T. E. Dyhoum 1,2, D. Lesnic 1, and R. G. Aykroyd 2 1 Department of Applied Mathematics,
More information6.2 Modeling of Systems and Components
Chapter 6 Modelling of Equipment, Processes, and Systems 61 Introduction Modeling is one of the most important elements of thermal system design Most systems are analyzed by considering equations which
More informationInfinite series, improper integrals, and Taylor series
Chapter 2 Infinite series, improper integrals, and Taylor series 2. Introduction to series In studying calculus, we have explored a variety of functions. Among the most basic are polynomials, i.e. functions
More informationY. M. TEMIS AND V. V. KARABAN 1
J. KSIAM Vol.5, No.2, 39-51, 2001 BOUNDARY ELEMENT TECHNIQUE IN TORSION PROBLEMS OF BEAMS WITH MULTIPLY CONNECTED CROSS-SECTIONS Y. M. TEMIS AND V. V. KARABAN 1 Abstract. This paper shows how boundary
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 informationThe matrix will only be consistent if the last entry of row three is 0, meaning 2b 3 + b 2 b 1 = 0.
) Find all solutions of the linear system. Express the answer in vector form. x + 2x + x + x 5 = 2 2x 2 + 2x + 2x + x 5 = 8 x + 2x + x + 9x 5 = 2 2 Solution: Reduce the augmented matrix [ 2 2 2 8 ] to
More informationIntegral Section List for OCR Specification
Integral Section List for OCR Specification 161 Sections are listed below, covering all of AS/A level Mathematics and Further Mathematics, categorised by Module and Topic Module Topic Section Name OCR
More informationIntegrated Math 3 Math 3 Course Description:
Course Description: Integrated Math 3 Math 3 Course Description: Integrated strands include algebra, functions, geometry, trigonometry, statistics, probability and discrete math. Scope and sequence includes
More informationAn Introduction to Numerical Methods for Differential Equations. Janet Peterson
An Introduction to Numerical Methods for Differential Equations Janet Peterson Fall 2015 2 Chapter 1 Introduction Differential equations arise in many disciplines such as engineering, mathematics, sciences
More informationSec 3 Express E-Math & A-Math Syllabus (For Class 301 to 305)
Sec 3 Express E-Math & A-Math Syllabus (For Class 301 to 305) Chapter 1 (EM) Quadratic Equations and Chapter 4 (EM) Coordinate Geometry Chapter 6 (EM) Further Trigonometry Chapter 2 (EM) Linear Inequalities
More informationME 1401 FINITE ELEMENT ANALYSIS UNIT I PART -A. 2. Why polynomial type of interpolation functions is mostly used in FEM?
SHRI ANGALAMMAN COLLEGE OF ENGINEERING AND TECHNOLOGY (An ISO 9001:2008 Certified Institution) SIRUGANOOR, TIRUCHIRAPPALLI 621 105 Department of Mechanical Engineering ME 1401 FINITE ELEMENT ANALYSIS 1.
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 informationName of the Student: Unit I (Solution of Equations and Eigenvalue Problems)
Engineering Mathematics 8 SUBJECT NAME : Numerical Methods SUBJECT CODE : MA6459 MATERIAL NAME : University Questions REGULATION : R3 UPDATED ON : November 7 (Upto N/D 7 Q.P) (Scan the above Q.R code for
More informationReview. Numerical Methods Lecture 22. Prof. Jinbo Bi CSE, UConn
Review Taylor Series and Error Analysis Roots of Equations Linear Algebraic Equations Optimization Numerical Differentiation and Integration Ordinary Differential Equations Partial Differential Equations
More informationSTEM-Prep Pathway SLOs
STEM-Prep Pathway SLOs Background: The STEM-Prep subgroup of the MMPT adopts a variation of the student learning outcomes for STEM from the courses Reasoning with Functions I and Reasoning with Functions
More informationPreface. 2 Linear Equations and Eigenvalue Problem 22
Contents Preface xv 1 Errors in Computation 1 1.1 Introduction 1 1.2 Floating Point Representation of Number 1 1.3 Binary Numbers 2 1.3.1 Binary number representation in computer 3 1.4 Significant Digits
More information11.2 Basic First-order System Methods
112 Basic First-order System Methods 797 112 Basic First-order System Methods Solving 2 2 Systems It is shown here that any constant linear system u a b = A u, A = c d can be solved by one of the following
More informationTaylor series. Chapter Introduction From geometric series to Taylor polynomials
Chapter 2 Taylor series 2. Introduction The topic of this chapter is find approximations of functions in terms of power series, also called Taylor series. Such series can be described informally as infinite
More informationChapter 1. Condition numbers and local errors in the boundary element method
Chapter 1 Condition numbers and local errors in the boundary element method W. Dijkstra, G. Kakuba and R. M. M. Mattheij Eindhoven University of Technology, Department of Mathematics and Computing Science,
More informationNumerical Methods for Engineers and Scientists
Numerical Methods for Engineers and Scientists Second Edition Revised and Expanded Joe D. Hoffman Department of Mechanical Engineering Purdue University West Lafayette, Indiana m MARCEL D E К К E R MARCEL
More informationFinal year project. Methods for solving differential equations
Final year project Methods for solving differential equations Author: Tej Shah Supervisor: Dr. Milan Mihajlovic Date: 5 th May 2010 School of Computer Science: BSc. in Computer Science and Mathematics
More informationLearning Module 1 - Basic Algebra Review (Appendix A)
Learning Module 1 - Basic Algebra Review (Appendix A) Element 1 Real Numbers and Operations on Polynomials (A.1, A.2) Use the properties of real numbers and work with subsets of the real numbers Determine
More informationChapter P: Preliminaries
Chapter P: Preliminaries Winter 2016 Department of Mathematics Hong Kong Baptist University 1 / 59 Preliminaries The preliminary chapter reviews the most important things that you should know before beginning
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 informationThe Learning Objectives of the Compulsory Part Notes:
17 The Learning Objectives of the Compulsory Part Notes: 1. Learning units are grouped under three strands ( Number and Algebra, Measures, Shape and Space and Data Handling ) and a Further Learning Unit.
More informationFast multipole boundary element method for the analysis of plates with many holes
Arch. Mech., 59, 4 5, pp. 385 401, Warszawa 2007 Fast multipole boundary element method for the analysis of plates with many holes J. PTASZNY, P. FEDELIŃSKI Department of Strength of Materials and Computational
More informationLinear Systems. Math A Bianca Santoro. September 23, 2016
Linear Systems Math A4600 - Bianca Santoro September 3, 06 Goal: Understand how to solve Ax = b. Toy Model: Let s study the following system There are two nice ways of thinking about this system: x + y
More informationA THEORETICAL INTRODUCTION TO NUMERICAL ANALYSIS
A THEORETICAL INTRODUCTION TO NUMERICAL ANALYSIS Victor S. Ryaben'kii Semyon V. Tsynkov Chapman &. Hall/CRC Taylor & Francis Group Boca Raton London New York Chapman & Hall/CRC is an imprint of the Taylor
More information(Linear equations) Applied Linear Algebra in Geoscience Using MATLAB
Applied Linear Algebra in Geoscience Using MATLAB (Linear equations) Contents Getting Started Creating Arrays Mathematical Operations with Arrays Using Script Files and Managing Data Two-Dimensional Plots
More informationUNC Charlotte Super Competition Level 3 Test March 4, 2019 Test with Solutions for Sponsors
. Find the minimum value of the function f (x) x 2 + (A) 6 (B) 3 6 (C) 4 Solution. We have f (x) x 2 + + x 2 + (D) 3 4, which is equivalent to x 0. x 2 + (E) x 2 +, x R. x 2 + 2 (x 2 + ) 2. How many solutions
More informationlecture 5: Finite Difference Methods for Differential Equations
lecture : Finite Difference Methods for Differential Equations 1 Application: Boundary Value Problems Example 1 (Dirichlet boundary conditions) Suppose we want to solve the differential equation u (x)
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 informationPre-Calculus and Trigonometry Capacity Matrix
Review Polynomials A1.1.4 A1.2.5 Add, subtract, multiply and simplify polynomials and rational expressions Solve polynomial equations and equations involving rational expressions Review Chapter 1 and their
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 informationCLASS NOTES Models, Algorithms and Data: Introduction to computing 2018
CLASS NOTES Models, Algorithms and Data: Introduction to computing 2018 Petros Koumoutsakos, Jens Honore Walther (Last update: June 11, 2018) IMPORTANT DISCLAIMERS 1. REFERENCES: Much of the material (ideas,
More informationIterative Methods for Ax=b
1 FUNDAMENTALS 1 Iterative Methods for Ax=b 1 Fundamentals consider the solution of the set of simultaneous equations Ax = b where A is a square matrix, n n and b is a right hand vector. We write the iterative
More information