Transactions on Modelling and Simulation vol 12, 1996 WIT Press, ISSN X
|
|
- Maximilian Terry
- 5 years ago
- Views:
Transcription
1 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 of the boundary element method requires consideration of three types of integral - simple non-singular integrals, singular integrals of two kinds and those integrals which are 'nearly singular' (Brebbia and Dominguez 1989). If linear elements are used all integrals may be performed exactly (Jaswon and Symm 1977), but the additional complexity of quadratic or cubic element programs appears to give a significant gain in efficiency, and such elements are therefore often utilised; These elements may be curved and can represent complex shapes and changes in boundary variables more effectively. Unfortunately, more complex elements require more complex integration routines than the routines for straight elements. In dealing with BEM codes simplicity is clearly desirable, particularly when attempting to show others the value of boundary methods. Non-singular integrals are almost invariably dealt with by a simple Gauss- Legendre quadrature and problems with near-singular integrals may be overcome by using the same simple integration scheme given a sufficient number of integration points and some restriction on the minimum distance between a given element and the nearest off-element node. The strongly singular integrals which appear in the usual BEM direct method may be evaluated using the row sum technique which may not be the most efficient available but has the virtue of being easy to explain and implement. This
2 234 Boundary Elements paper implements a recently developed one part Gauss integration scheme for singular logarithmic integrals which is both simpler to implement and more accurate than the usual two-part Gauss approach with separate Gauss- Legendre and logarithmically weighted Gauss schemes. We aim to present a description for the quadratic element which is self-contained and suitable as an introduction to this aspect of BEM. Quadratic elements Consider an isoparametric quadratic element defined on [0,1] as in figure #, # Figure 1: A straight quadratic element The variation of x, y, potential u arid potential derivative q du/dn are all expressed in the same form in terms of their nodal values, for example: 4-92^2 + %Ms (1) where the Lagrangian interpolation functions MI, A/2,A/3 are M, = (2f -!)( -!) Mi = 4f(l - f) -1) In one dimension, a typical boundary element integral of q over an element Si on the boundary has the form (omitting a l/(2ir) multiplier) / = / ln(r) q ds (2) J S, where ds is the element of arc length around the boundary S and r = i(x - xo)2 + (y - 3/0)2 is the distance between the source point (xo,yo) and the boundary point (,r,y) on the current element.
3 Boundary Elements 235 If g, r and ds are now all expressed in terms of the parameter using an isoparametric quadratic element then the integral (2) becomes (3) where J(^) is the Jacobian transforming ck to Jo If the point (xo,yo) coincides with the first node of the quadratic element, then r 0 at this point and the integral including MI above is weakly singular. The other two shape functions MI and M^ are zero at f = 0 and those integrals are not therefore singular for this position of the source node - although they will be singular for the source node at node 2 and node 3 respectively. We discuss here only the case for a singularity on node 1, since integrals singular on the other two nodes are simple transformations/combinations of the node 1 routine. Numerical integration Very simple elements allow the use of exact integration but for more complex elements, numerical integration is used. In general, an n-point Gauss- Legendre formula with abscissae x, and weights W{ (%)(& = ]T %,,/(%,) (4) i=i will give a result for the integral which is exact for all polynomials of degree 2n 1 (neglecting any rounding errors). The high efficiency of such integration schemes mean that they are almost invariably used for numerical integration in FEM and BEM. One weighted Gauss scheme which is particularly useful in the BEM is for integrals which include a log( ) term on the interval [0,1]. With abscissae Xi and weights W{ this scheme is written In order to incorporate this weighted Gauss scheme into the integral (3) we must isolate the part of? ( ). One way of doing this (Smith and Mason
4 236 Boundary Elements 1982) is as follows: When x is expressed as an isoparametric quadratic on [0,1] we can write where z = A^ 4- Bf + zi; % = Cf 4- Df + 2/1 C = 2%/i - 4?/2 + 2?/3, D = -3m + 4?/2-7/3 and if the element is singular on node 1 then r% = ^[(A" +C')(" + (2AB 4- The ln(r) integral singular on node 1 can thus be written as Jo This kind of separation by direct factorisation into the two components avoids division by small quantities or taking square roots and would thus seem a more robust approach than is sometimes used. The first integral on the right-hand side of (6) is evaluated using the weighted scheme (5) and the second non-singular integral using the simple quadrature (4). This two-part gauss quadrature with n points for each part is, of course, exact for polynomials of degree (2n 1) and therefore for integrals on straight quadratic elements when n > 1. An alternative scheme for the singular logarithmic integral is based on subtracting out the singularity after expanding the integrand using Taylor series (Aliabadi et al 1985) Whilst this approach is effective, it is substantially more difficult to explain or implement than a two-part Gauss quadrature. Crow (1993) developed a gauss quadrature scheme for a function of the form: which is exact where a and b are polynomials of degree at most N. Choosing a 0 and 6=0 separately leads to the equations: TV N 0 0 = / Jo i
5 Boundary Elements 237 The choice of basis polynomials for this problem is arbitrary since the unknown abscissae are not the roots of orthogonal polynomials but Legendre polynomials P*(x] on [0,1] produce lower condition numbers than powers of x in the following scheme: E o giving 2N+2 equations for the 2N+2 unknowns: Given an initial estimate for, the first equation above is solved for w and the second equation for in an iterative fashion. The computation of quadrature points in this way seems to be particularly ill-conditioned and data is only available for N up to N = 7. The quadrature points for the case N 6 are given below for convenience: The above quadrature scheme X w X w Table 1: Gauss quadrature points covers the requirements for the BEM ln(r) singularity without the need to explicitly form the two components of the integrand in the singular case. Results The two-part log procedure and the single gauss quadrature are now compared for several elements with different profiles. Since exact solutions are not generally available, accurate comparison values were obtained from the adaptive numerical integration provided by the symbolic algebra package. The other results were obtained using single precision arithmetic (c± 7dp) on a DEC MICRO VAX and include a straightforward 6-point Gauss-Legendre
6 238 Boundary Elements rule which takes no account of the singularity. Both the singular rules also use 6 points, the two-part rule with 3 log -f 3 Legendre points. The integrals from all 3 nodes are treated in the same way although strictly only the first node integral is singular. nodes 1,2,3 x-coord 6-pt Gauss Two part log Crow 6-pt Table 2: log integrals on a straight element Table 2 gives integrals for a linear geometry quadratic element which can be evaluated exactly since the Jacobian is a constant. The singular numerical schemes should also be 'exact' to within numerical errors since the unaccounted for polynomial in the integrand is just the quadratic shape function of degree n = 2 and these schemes should therefore be exact with (2n 1) > 2. This is borne out by table 2 which shows all 3 methods apart from the simple Gauss rule agreeing to machine precision. Observe that although the integrals from node 2 and node 3 are non-singular, the singular rules appear to give significantly better estimates than simple Gauss- Legendre. Once the element deviates from a pure linear element, exact nodes 1,2,3 x-coord 6-pt Gauss Two part log Crow 6-pt Table 3: log integrals - small curvature integration is no longer generally possible and tables 3 and 4 show that the Crow rule performs significantly better than the other schemes; accurate to 5 or 6 places of decimals compared to 3 or 4 otherwise. As the element curvature increases, integrals become less accurate and in table 5 accuracy may be only 3 decimal places but the Crow scheme is still better than other rules considered. Although it is usual in two-part schemes to use the same number of points in each part, there is no particular reason for this, other than convenience. Indeed, the second non-singular part of the
7 Boundary Elements 239 nodes 1,2.3 x-coord 6-pt Gauss Two part log Crow 6-pt Table 4: log integrals - moderate curvature nodes 1,2,3 x-coord 6-pt Gauss Two part log Crow 6-pt Table 5: log integrals - high curvature integrand must represent a more complex expression because it includes an explicit log term which is represented implicitly by the weight in the singular part. Table 6 confirms this as 4 Gauss-Legendre (g) and 3 logarithmically weighted points (1) are much better than the converse. A Crow style scheme should be possible which takes into account the different numbers of points required. nodes 1,2,3 x-coord Two part log 3g+31 Two part log 3g+41 Two part log 4g-f Table 6: log integrals - unequal points
8 240 Boundary Elements Conclusions The special gauss quadrature used here for the evaluation of singular logarithmic integrals seems to be the simplest to teach and to program, as well as being more accurate than the conventional two-part gauss quadrature. However, abscissae and weights for larger numbers of integration points are needed for higher accuracy on quadratic elements or for integration on higher order elements. The logarithmic integrals on the other nodes of an element which are not strictly speaking singular, are nevertheless evaluated more accurately by the singular schemes used here than by a straightforward Gauss-Legendre scheme. More points will generally be needed for the nonsingular part of the integration scheme than for the singular part. Acknowledgements I would like to thank G J Fernandes for computing some example data. References Aliabadi M H, Hall W S and Phemister T G, (1985) Taylor expansions for singular kernels in the boundary element method. Int J Num Meth Engg, 21, Brebbia, C.A. and Dominguez, J. Boundary Elements: An Introductory Course, Computational Mechanics Publications, Southampton, Crow, J.A. (1993) Quadrature of integrands with a logarithmic singularity. Math. Comp. 60, 201, \ J as won, M.A. and Symm, G.T. (1977) Integral Equation Methods in Potential Theory and Elastostatics. Academic Press, London. Smith R.N.L. and Mason J C (1982) A Boundary element method for curved crack problems in two dimensions, Proc. 4th Int Conf. on Boundary Elements (Ed. Brebbia C A), Springer-Verlag Berlin, pp
The use of exact values at quadrature points in the boundary element method
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 E-mail:
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 informationJUST THE MATHS UNIT NUMBER 1.6. ALGEBRA 6 (Formulae and algebraic equations) A.J.Hobson
JUST THE MATHS UNIT NUMBER 1.6 ALGEBRA 6 (Formulae and algebraic equations) by A.J.Hobson 1.6.1 Transposition of formulae 1.6. of linear equations 1.6.3 of quadratic equations 1.6. Exercises 1.6.5 Answers
More informationINTERPOLATION. and y i = cos x i, i = 0, 1, 2 This gives us the three points. Now find a quadratic polynomial. p(x) = a 0 + a 1 x + a 2 x 2.
INTERPOLATION Interpolation is a process of finding a formula (often a polynomial) whose graph will pass through a given set of points (x, y). As an example, consider defining and x 0 = 0, x 1 = π/4, x
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 informationTransactions on Modelling and Simulation vol 18, 1997 WIT Press, ISSN X
A direct formulation for three-dimensional elastoplastic boundary elements A. P. Cisilino*,**, M. H. Aliabadi* and J. L. Otegui** *Wessex Institute of Technology, As hurst Lodge, As hurst, Southampton
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 informationTransactions on Modelling and Simulation vol 8, 1994 WIT Press, ISSN X
Analysis of singular and near-singular integrals by the continuation approach D. Rosen & D.E. Cormack Department of Chemical Engineering and Applied Chemistry, University of Toronto, Toronto M5S 1A4, Canada
More informationLEAST SQUARES APPROXIMATION
LEAST SQUARES APPROXIMATION One more approach to approximating a function f (x) on an interval a x b is to seek an approximation p(x) with a small average error over the interval of approximation. A convenient
More informationDepartment of Structural, Faculty of Civil Engineering, Architecture and Urban Design, State University of Campinas, Brazil
Blucher Mechanical Engineering Proceedings May 2014, vol. 1, num. 1 www.proceedings.blucher.com.br/evento/10wccm A SIMPLIFIED FORMULATION FOR STRESS AND TRACTION BOUNDARY IN- TEGRAL EQUATIONS USING THE
More informationThe iterated sinh transformation
The iterated sinh transformation Author Elliott, David, Johnston, Peter Published 2008 Journal Title International Journal for Numerical Methods in Engineering DOI https://doi.org/10.1002/nme.2244 Copyright
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 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 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 informationNumerical Integration (Quadrature) Another application for our interpolation tools!
Numerical Integration (Quadrature) Another application for our interpolation tools! Integration: Area under a curve Curve = data or function Integrating data Finite number of data points spacing specified
More informationPart IB - Easter Term 2003 Numerical Analysis I
Part IB - Easter Term 2003 Numerical Analysis I 1. Course description Here is an approximative content of the course 1. LU factorization Introduction. Gaussian elimination. LU factorization. Pivoting.
More informationelastoplastic contact problems D. Martin and M.H. Aliabadi Wessex Institute of Technology, Ashurst Lodge, Ashurst, Southampton, SO40 7AA, UK
Non-conforming BEM elastoplastic contact problems D. Martin and M.H. Aliabadi discretisation in Wessex Institute of Technology, Ashurst Lodge, Ashurst, Southampton, SO40 7AA, UK Abstract In this paper,
More informationChapter 10 INTEGRATION METHODS Introduction Unit Coordinate Integration
Chapter 0 INTEGRATION METHODS 0. Introduction The finite element analysis techniques are always based on an integral formulation. At the very minimum it will always be necessary to integrate at least an
More informationƒ f(x)dx ~ X) ^i,nf(%i,n) -1 *=1 are the zeros of P«(#) and where the num
ZEROS OF THE HERMITE POLYNOMIALS AND WEIGHTS FOR GAUSS' MECHANICAL QUADRATURE FORMULA ROBERT E. GREENWOOD AND J. J. MILLER In the numerical integration of a function ƒ (x) it is very desirable to choose
More informationExample 1 Which of these functions are polynomials in x? In the case(s) where f is a polynomial,
1. Polynomials A polynomial in x is a function of the form p(x) = a 0 + a 1 x + a 2 x 2 +... a n x n (a n 0, n a non-negative integer) where a 0, a 1, a 2,..., a n are constants. We say that this polynomial
More informationInterpolation and Approximation
Interpolation and Approximation The Basic Problem: Approximate a continuous function f(x), by a polynomial p(x), over [a, b]. f(x) may only be known in tabular form. f(x) may be expensive to compute. Definition:
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 informationAlgebraic Statistics progress report
Algebraic Statistics progress report Joe Neeman December 11, 2008 1 A model for biochemical reaction networks We consider a model introduced by Craciun, Pantea and Rempala [2] for identifying biochemical
More informationScientific Computing
2301678 Scientific Computing Chapter 2 Interpolation and Approximation Paisan Nakmahachalasint Paisan.N@chula.ac.th Chapter 2 Interpolation and Approximation p. 1/66 Contents 1. Polynomial interpolation
More informationJUST THE MATHS UNIT NUMBER 1.5. ALGEBRA 5 (Manipulation of algebraic expressions) A.J.Hobson
JUST THE MATHS UNIT NUMBER 1.5 ALGEBRA 5 (Manipulation of algebraic expressions) by A.J.Hobson 1.5.1 Simplification of expressions 1.5.2 Factorisation 1.5.3 Completing the square in a quadratic expression
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 informationIntroduction Linear system Nonlinear equation Interpolation
Interpolation Interpolation is the process of estimating an intermediate value from a set of discrete or tabulated values. Suppose we have the following tabulated values: y y 0 y 1 y 2?? y 3 y 4 y 5 x
More informationNumerical Integration for Multivariable. October Abstract. We consider the numerical integration of functions with point singularities over
Numerical Integration for Multivariable Functions with Point Singularities Yaun Yang and Kendall E. Atkinson y October 199 Abstract We consider the numerical integration of functions with point singularities
More informationNUMERICAL COMPUTATION IN SCIENCE AND ENGINEERING
NUMERICAL COMPUTATION IN SCIENCE AND ENGINEERING C. Pozrikidis University of California, San Diego New York Oxford OXFORD UNIVERSITY PRESS 1998 CONTENTS Preface ix Pseudocode Language Commands xi 1 Numerical
More informationNumerical Methods of Approximation
Contents 31 Numerical Methods of Approximation 31.1 Polynomial Approximations 2 31.2 Numerical Integration 28 31.3 Numerical Differentiation 58 31.4 Nonlinear Equations 67 Learning outcomes In this Workbook
More informationCOMPUTATION OF BESSEL AND AIRY FUNCTIONS AND OF RELATED GAUSSIAN QUADRATURE FORMULAE
BIT 6-85//41-11 $16., Vol. 4, No. 1, pp. 11 118 c Swets & Zeitlinger COMPUTATION OF BESSEL AND AIRY FUNCTIONS AND OF RELATED GAUSSIAN QUADRATURE FORMULAE WALTER GAUTSCHI Department of Computer Sciences,
More informationAN ALTERNATIVE TECHNIQUE FOR TANGENTIAL STRESS CALCULATION IN DISCONTINUOUS BOUNDARY ELEMENTS
th Pan-American Congress of Applied Mechanics January 04-08, 00, Foz do Iguaçu, PR, Brazil AN ALTERNATIVE TECHNIQUE FOR TANGENTIAL STRESS CALCULATION IN DISCONTINUOUS BOUNDARY ELEMENTS Otávio Augusto Alves
More informationAn Arithmetic-Geometric Mean Inequality
EXPOSITIONES MATHEMATICAE Mathematical Note Expo. Math. 19 (2001):273-279 Urban & Fischer Verlag www. u rbanfischer.de/jou rnals/expomath An Arithmetic-Geometric Mean Inequality Paul Bracken Centre de
More informationJUST THE MATHS UNIT NUMBER 1.9. ALGEBRA 9 (The theory of partial fractions) A.J.Hobson
JUST THE MATHS UNIT NUMBER 1. ALGEBRA (The theory of partial fractions) by A.J.Hobson 1..1 Introduction 1..2 Standard types of partial fraction problem 1.. Exercises 1..4 Answers to exercises UNIT 1. -
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 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 informationA Gauss Lobatto quadrature method for solving optimal control problems
ANZIAM J. 47 (EMAC2005) pp.c101 C115, 2006 C101 A Gauss Lobatto quadrature method for solving optimal control problems P. Williams (Received 29 August 2005; revised 13 July 2006) Abstract This paper proposes
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 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 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 informationChapter 1 Error Analysis
Chapter 1 Error Analysis Several sources of errors are important for numerical data processing: Experimental uncertainty: Input data from an experiment have a limited precision. Instead of the vector of
More informationFirst- and second-order sensitivity analysis schemes by collocation-type Trefftz method
Computer Assisted Mechanics and Engineering Sciences, 4: 477 489, 1997 Copyright c 1997 by Polska Akademia Nauk First- and second-order sensitivity analysis schemes by collocation-type Trefftz method Eisuke
More informationHIGHER ORDER RIGIDITY - WHAT IS THE PROPER DEFINITION?
HIGHER ORDER RIGIDITY - WHAT IS THE PROPER DEFINITION? ROBERT CONNELLY AND HERMAN SERVATIUS Abstract. We show that there is a bar and joint framework G(p) which has a configuration p in the plane such
More information22.3. Repeated Eigenvalues and Symmetric Matrices. Introduction. Prerequisites. Learning Outcomes
Repeated Eigenvalues and Symmetric Matrices. Introduction In this Section we further develop the theory of eigenvalues and eigenvectors in two distinct directions. Firstly we look at matrices where one
More informationMidterm Review. Igor Yanovsky (Math 151A TA)
Midterm Review Igor Yanovsky (Math 5A TA) Root-Finding Methods Rootfinding methods are designed to find a zero of a function f, that is, to find a value of x such that f(x) =0 Bisection Method To apply
More informationThe Normal Equations. For A R m n with m > n, A T A is singular if and only if A is rank-deficient. 1 Proof:
Applied Math 205 Homework 1 now posted. Due 5 PM on September 26. Last time: piecewise polynomial interpolation, least-squares fitting Today: least-squares, nonlinear least-squares The Normal Equations
More informationHIGHER MATHS REVISION CHECKLIST (Grades 9 4)
HIGHER MATHS REVISION CHECKLIST 2017+ (s 9 4) Geometry and Measures Circle theorems 8 Vector arguments and proof 8 Area of a triangle 7 Cosine Rule 7 Pythagoras and trig 2D and 3D 7 Sine Rule 7 Combined
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 informationTransactions on Modelling and Simulation vol 9, 1995 WIT Press, ISSN X
A path-independent integral for the calculation of stress intensity factors in three-dimensional bodies C.M. Bainbridge," M.H. Aliabadi," D.P. Rooke* "Wessex Institute of Technology, Ashurst Lodge, Ashurst,
More information(Mathematical Operations with Arrays) Applied Linear Algebra in Geoscience Using MATLAB
Applied Linear Algebra in Geoscience Using MATLAB (Mathematical Operations with Arrays) Contents Getting Started Matrices Creating Arrays Linear equations Mathematical Operations with Arrays Using Script
More informationIntroduction to Quantitative Techniques for MSc Programmes SCHOOL OF ECONOMICS, MATHEMATICS AND STATISTICS MALET STREET LONDON WC1E 7HX
Introduction to Quantitative Techniques for MSc Programmes SCHOOL OF ECONOMICS, MATHEMATICS AND STATISTICS MALET STREET LONDON WC1E 7HX September 2007 MSc Sep Intro QT 1 Who are these course for? The September
More informationFINITE-DIMENSIONAL LINEAR ALGEBRA
DISCRETE MATHEMATICS AND ITS APPLICATIONS Series Editor KENNETH H ROSEN FINITE-DIMENSIONAL LINEAR ALGEBRA Mark S Gockenbach Michigan Technological University Houghton, USA CRC Press Taylor & Francis Croup
More informationDifferentiation and Integration
Differentiation and Integration (Lectures on Numerical Analysis for Economists II) Jesús Fernández-Villaverde 1 and Pablo Guerrón 2 February 12, 2018 1 University of Pennsylvania 2 Boston College Motivation
More informationNon-polynomial Least-squares fitting
Applied Math 205 Last time: piecewise polynomial interpolation, least-squares fitting Today: underdetermined least squares, nonlinear least squares Homework 1 (and subsequent homeworks) have several parts
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 informationGaussian interval quadrature rule for exponential weights
Gaussian interval quadrature rule for exponential weights Aleksandar S. Cvetković, a, Gradimir V. Milovanović b a Department of Mathematics, Faculty of Mechanical Engineering, University of Belgrade, Kraljice
More informationSecondary Math 3 Honors - Polynomial and Polynomial Functions Test Review
Name: Class: Date: Secondary Math 3 Honors - Polynomial and Polynomial Functions Test Review 1 Write 3x 2 ( 2x 2 5x 3 ) in standard form State whether the function is even, odd, or neither Show your work
More informationCh. 03 Numerical Quadrature. Andrea Mignone Physics Department, University of Torino AA
Ch. 03 Numerical Quadrature Andrea Mignone Physics Department, University of Torino AA 2017-2018 Numerical Quadrature In numerical analysis quadrature refers to the computation of definite integrals. y
More informationStudies on Barlow points, Gauss points and Superconvergent points in 1D with Lagrangian and Hermitian finite element basis
Volume 7, N., pp. 75 0, 008 Copyright 008 SBMAC ISSN 0101-805 www.scielo.br/cam Studies on Barlow points, Gauss points and Superconvergent points in 1D with Lagrangian and Hermitian finite element basis
More informationA Method for Reducing Ill-Conditioning of Polynomial Root Finding Using a Change of Basis
Portland State University PDXScholar University Honors Theses University Honors College 2014 A Method for Reducing Ill-Conditioning of Polynomial Root Finding Using a Change of Basis Edison Tsai Portland
More informationMATH 304 Linear Algebra Lecture 20: Review for Test 1.
MATH 304 Linear Algebra Lecture 20: Review for Test 1. Topics for Test 1 Part I: Elementary linear algebra (Leon 1.1 1.4, 2.1 2.2) Systems of linear equations: elementary operations, Gaussian elimination,
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 informationOn the use of multipole methods for domain integration in the BEM
On the use of multipole methods for domain integration in the BEM A.A. Mammoii, M.S. Ingber & M.J. Brown Department of Mechanical Engineering, University of New Mexico, USA Abstract Often, a single application
More informationFOUNDATION MATHS REVISION CHECKLIST (Grades 5 1)
FOUNDATION MATHS REVISION CHECKLIST 2017+ (s 5 1) Geometry and Measures Arc lengths and sectors 5 Derive triangle results 5 Enlargements and negative SF 5 Loci 5 Pythagoras 5 Similarity and Congruence
More informationLEAST SQUARES DATA FITTING
LEAST SQUARES DATA FITTING Experiments generally have error or uncertainty in measuring their outcome. Error can be human error, but it is more usually due to inherent limitations in the equipment being
More informationTransactions on Modelling and Simulation vol 12, 1996 WIT Press, ISSN X
Plate-soil elastodynamic coupling using analysis S.F.A. Baretto, H.B. Coda, W.S. Venturini Sao Carlos School of Engineering, University ofsao Paulo, Sao Carlos - SP, Brazil BEM Abstract The aim of this
More informationEngg. Math. II (Unit-IV) Numerical Analysis
Dr. Satish Shukla of 33 Engg. Math. II (Unit-IV) Numerical Analysis Syllabus. Interpolation and Curve Fitting: Introduction to Interpolation; Calculus of Finite Differences; Finite Difference and Divided
More informationCS 450 Numerical Analysis. Chapter 5: Nonlinear Equations
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 informationComputer Arithmetic. MATH 375 Numerical Analysis. J. Robert Buchanan. Fall Department of Mathematics. J. Robert Buchanan Computer Arithmetic
Computer Arithmetic MATH 375 Numerical Analysis J. Robert Buchanan Department of Mathematics Fall 2013 Machine Numbers When performing arithmetic on a computer (laptop, desktop, mainframe, cell phone,
More informationNewton-Raphson Type Methods
Int. J. Open Problems Compt. Math., Vol. 5, No. 2, June 2012 ISSN 1998-6262; Copyright c ICSRS Publication, 2012 www.i-csrs.org Newton-Raphson Type Methods Mircea I. Cîrnu Department of Mathematics, Faculty
More informationLecture Notes to Accompany. Scientific Computing An Introductory Survey. by Michael T. Heath. Chapter 5. Nonlinear Equations
Lecture Notes to Accompany Scientific Computing An Introductory Survey Second Edition by Michael T Heath Chapter 5 Nonlinear Equations Copyright c 2001 Reproduction permitted only for noncommercial, educational
More informationNumerical Integra/on
Numerical Integra/on Applica/ons The Trapezoidal Rule is a technique to approximate the definite integral where For 1 st order: f(a) f(b) a b Error Es/mate of Trapezoidal Rule Truncation error: From Newton-Gregory
More informationLECTURE NOTES ELEMENTARY NUMERICAL METHODS. Eusebius Doedel
LECTURE NOTES on ELEMENTARY NUMERICAL METHODS Eusebius Doedel TABLE OF CONTENTS Vector and Matrix Norms 1 Banach Lemma 20 The Numerical Solution of Linear Systems 25 Gauss Elimination 25 Operation Count
More informationLinear Algebra. Min Yan
Linear Algebra Min Yan January 2, 2018 2 Contents 1 Vector Space 7 1.1 Definition................................. 7 1.1.1 Axioms of Vector Space..................... 7 1.1.2 Consequence of Axiom......................
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 informationTwo hours. To be provided by Examinations Office: Mathematical Formula Tables. THE UNIVERSITY OF MANCHESTER. 29 May :45 11:45
Two hours MATH20602 To be provided by Examinations Office: Mathematical Formula Tables. THE UNIVERSITY OF MANCHESTER NUMERICAL ANALYSIS 1 29 May 2015 9:45 11:45 Answer THREE of the FOUR questions. If more
More informationOutline. Scientific Computing: An Introductory Survey. Nonlinear Equations. Nonlinear Equations. Examples: Nonlinear Equations
Methods for Systems of Methods for Systems of Outline Scientific Computing: An Introductory Survey Chapter 5 1 Prof. Michael T. Heath Department of Computer Science University of Illinois at Urbana-Champaign
More informationChordal cubic spline interpolation is fourth order accurate
Chordal cubic spline interpolation is fourth order accurate Michael S. Floater Abstract: It is well known that complete cubic spline interpolation of functions with four continuous derivatives is fourth
More informationAn Introduction to Differential Algebra
An Introduction to Differential Algebra Alexander Wittig1, P. Di Lizia, R. Armellin, et al. 1 ESA Advanced Concepts Team (TEC-SF) SRL, Milan Dinamica Outline 1 Overview Five Views of Differential Algebra
More informationModule 6: Implicit Runge-Kutta Methods Lecture 17: Derivation of Implicit Runge-Kutta Methods(Contd.) The Lecture Contains:
The Lecture Contains: We continue with the details about the derivation of the two stage implicit Runge- Kutta methods. A brief description of semi-explicit Runge-Kutta methods is also given. Finally,
More informationwhich arises when we compute the orthogonal projection of a vector y in a subspace with an orthogonal basis. Hence assume that P y = A ij = x j, x i
MODULE 6 Topics: Gram-Schmidt orthogonalization process We begin by observing that if the vectors {x j } N are mutually orthogonal in an inner product space V then they are necessarily linearly independent.
More informationEigenvalues and Eigenvectors
Contents Eigenvalues and Eigenvectors. Basic Concepts. Applications of Eigenvalues and Eigenvectors 8.3 Repeated Eigenvalues and Symmetric Matrices 3.4 Numerical Determination of Eigenvalues and Eigenvectors
More informationSolving Equations. Solving Equations - decimal coefficients and constants. 2) Solve for x: 3(3x 6) = 3(x -2) 1) Solve for x: 5 x 2 28 x
Level C Review Packet This packet briefly reviews the topics covered on the Level A Math Skills Assessment. If you need additional study resources and/or assistance with any of the topics below, please
More informationCOMP 558 lecture 18 Nov. 15, 2010
Least squares We have seen several least squares problems thus far, and we will see more in the upcoming lectures. For this reason it is good to have a more general picture of these problems and how to
More informationSemi-implicit Krylov Deferred Correction Methods for Ordinary Differential Equations
Semi-implicit Krylov Deferred Correction Methods for Ordinary Differential Equations Sunyoung Bu University of North Carolina Department of Mathematics CB # 325, Chapel Hill USA agatha@email.unc.edu Jingfang
More informationTECHNIQUES IN FACTORISATION
TECHNIQUES IN FACTORISATION The process where brackets are inserted into an equation is referred to as factorisation. Factorisation is the opposite process to epansion. METHOD: Epansion ( + )( 5) 15 Factorisation
More informationJim Lambers MAT 610 Summer Session Lecture 2 Notes
Jim Lambers MAT 610 Summer Session 2009-10 Lecture 2 Notes These notes correspond to Sections 2.2-2.4 in the text. Vector Norms Given vectors x and y of length one, which are simply scalars x and y, the
More informationSIMULATION OF PLANE STRAIN FIBER COMPOSITE PLATES IN BENDING THROUGH A BEM/ACA/HM FORMULATION
8 th GRACM International Congress on Computational Mechanics Volos, 12 July 15 July 2015 SIMULATION OF PLANE STRAIN FIBER COMPOSITE PLATES IN BENDING THROUGH A BEM/ACA/HM FORMULATION Theodore V. Gortsas
More information= W z1 + W z2 and W z1 z 2
Math 44 Fall 06 homework page Math 44 Fall 06 Darij Grinberg: homework set 8 due: Wed, 4 Dec 06 [Thanks to Hannah Brand for parts of the solutions] Exercise Recall that we defined the multiplication of
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 informationA matrix over a field F is a rectangular array of elements from F. The symbol
Chapter MATRICES Matrix arithmetic A matrix over a field F is a rectangular array of elements from F The symbol M m n (F ) denotes the collection of all m n matrices over F Matrices will usually be denoted
More informationZero Energy Modes in One Dimension: An Introduction to Hourglass Modes
Zero Energy Modes in One Dimension: An Introduction to Hourglass Modes David J. Benson March 9, 2003 Reduced integration does a lot of good things for an element: it reduces the computational cost, it
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 informationVibration of Thin Beams by PIM and RPIM methods. *B. Kanber¹, and O. M. Tufik 1
APCOM & ISCM -4 th December, 23, Singapore Vibration of Thin Beams by PIM and RPIM methods *B. Kanber¹, and O. M. Tufik Mechanical Engineering Department, University of Gaziantep, Turkey. *Corresponding
More informationIntroduction to Techniques for Counting
Introduction to Techniques for Counting A generating function is a device somewhat similar to a bag. Instead of carrying many little objects detachedly, which could be embarrassing, we put them all in
More informationApplied Numerical Analysis
Applied Numerical Analysis Using MATLAB Second Edition Laurene V. Fausett Texas A&M University-Commerce PEARSON Prentice Hall Upper Saddle River, NJ 07458 Contents Preface xi 1 Foundations 1 1.1 Introductory
More informationAbsolute Value Programming
O. L. Mangasarian Absolute Value Programming Abstract. We investigate equations, inequalities and mathematical programs involving absolute values of variables such as the equation Ax + B x = b, where A
More informationExpressing a Rational Fraction as the sum of its Partial Fractions
PARTIAL FRACTIONS Dear Reader An algebraic fraction or a rational fraction can be, often, expressed as the algebraic sum of relatively simpler fractions called partial fractions. The application of partial
More informationOn the positivity of linear weights in WENO approximations. Abstract
On the positivity of linear weights in WENO approximations Yuanyuan Liu, Chi-Wang Shu and Mengping Zhang 3 Abstract High order accurate weighted essentially non-oscillatory (WENO) schemes have been used
More informationBrockington College Mathematics Personal Learning Checklist
Brockington College Mathematics Personal Learning Checklist To help you use this personal learning checklist, the target levels for each topic have given to help you decide what to focus on for your tier
More information