Boundary Element Method (BEM)

Size: px
Start display at page:

Download "Boundary Element Method (BEM)"

Transcription

1 Boudary Elemet Method BEM Zora Ilievski Wedesday 8 th Jue 006 HG 6.96 TU/e

2 Talk Overview The idea of BEM ad its advatages The D potetial problem Numerical implemetatio

3 Idea of BEM 3

4 Idea of BEM 4

5 Advatages of BEM Reductio of problem dimesio by Less data preparatio time. Easier to chage the applied mesh. Useful for problems that require re-meshig. 5

6 Advatages of BEM High Accuracy Stresses are accurate as there are o approximatios imposed o the solutio i iterior domai poits. Suitable for modelig problems of rapidly chagig stresses. 6

7 Advatages of BEM 3 Less computer time ad storage For the same level of accuracy as other methods BEM uses less umber of odes ad elemets. 7

8 Advatages of BEM 4 Filter out uwated iformatio. Iteral poits of the domai are optioal. Focus o particular iteral regio. Further reduces computer time. 8

9 Advatages of BEM. Reductio of problem dimesio by.. High Accuracy. 3. Less computer time ad storage. 4. Filter out uwated iformatio ad so focus o sectio of the domai you are iterested i. BEM is a attractive optio. 9

10 The D potetial problem Where ca BEM be applied? Two importat fuctios. Descriptio of the domai. Mappig of higher to lower dimesios. Satisfactio of the Laplace equatios ad how to deal with a sigularity. The boudary itegral equatio BIE 0

11 The D potetial problem Where ca BEM be applied? Where ay potetial problem is govered by a differetial equatio that satisfies the Laplace equatio. or ay other behavior that has a related fudametal solutio e.g. The followig ca be aalyzed with the Laplace equatio: fluid flow torsio of bars diffusio ad steady state heat coductio.

12 The D potetial problem The Laplace equatio for D + x y 0. x y Laplacia operator otetial fuctio Cartesia coordiate axis

13 The D potetial problem Two importat fuctios. The fuctio describig the property uder aalysis. e.g. heat. Ukow λ The fudametal solutio of the Laplace equatio. These are well kow 3

14 4 The D potetial problem l p r p π λ l p r p π λ Fudametal solutio of the D Laplace equatio for a cocetrated source poit at p is p p y Y x X p r + Where Descriptio of the domai

15 The D potetial problem Mappig of higher to lower dimesios Boudary of ay domai is of a dimesio less tha of the domai. I BEM the problem is moved from withi the domai to its boudary. This meas you must i this case map Area to Lie. The well kow Grees Secod Idetity is used to do this. λ λ λ λ da λ λ d A have cotiuous st ad d derivatives. ukow potetial at ay poit. kow fudametal solutio at ay poit. uit outward ormal. derivative i the directio of ormal. 5

16 The D potetial problem Satisfyig the Laplace equatio The ukow will satisfy The kow fudametal solutio satisfies λ 0 λ except the poit p where it is sigular. 0 everywhere i the solutio domai. everywhere λ p l π r p X x + Y y r p p p 6

17 The D potetial problem How to deal with the sigularity Surroud p with a small circle of radius ε the examie solutio as ε à 0 New area is A Aε New boudary is + ε λ λ da λ d A Aε + ε Withi area A Aε 0 & λ 0 λ The left had side of the equatio is ow 0 ad the right is ow 7

18 8 The D potetial problem How to deal with the sigularity + d d ε λ λ λ λ 0 The secod term must be evaluated ad to do this let α d εd r r r π λ λ. Ad use the fact that

19 9 The D potetial problem How to deal with the sigularity π π α ε ε ε π λ λ π ε. l 0 d d 0 / C Evaluated with p i the domai o the boudary Smooth surface ad outside the boudary. π θ C For coarse surfaces

20 0 The D potetial problem The boudary itegral equatio d K d K C K K λ λ Where K ad K are the kow fudametal solutios ad are equal to π θ C

21 The D potetial problem BEM ca be applied where ay potetial problem is govered by a differetial equatio that satisfies the Laplace equatio. I this case the D form. A potetial problem ca be mapped from higher to lower dimesio usig Gree s secod idetity. Show how to deal with the case of the sigularity poit. Derived the boudary itegral equatio BIE C K d K d

22 Numerical Implemetatio Dirichlet Neuma ad mixed case. Discretisatio Reductio to a form AxB

23 Numerical Implemetatio Dirichlet Neuma ad mixed case. C K d K d The ukows of the above are values o the boudary ad are Dirichlet roblem is give every poit o the boudary. Neuma roblem is give every poit o the boudary. Mixed case Either are give at poit 3

24 4 Numerical Implemetatio Discretisatio N i N i i d K d K Ukows

25 Numerical Implemetatio Discretisatio Let K i K i d K i K d i Ukows N N i Ki K i 5

26 Numerical Implemetatio i whe i N K + i δi i N K Ax Bz 6

27 Numerical Implemetatio Neuma roblem Ax c Matrix A ad vector C are kow Dirichlet roblem c Bz Matrix B ad vector C are kow Mixed case Ax Bz Ukows ad kows ca be separated i to same form as above 7

28 8 Numerical Implemetatio As each poit p i the domai is expressed i terms of the boudary values oce all boudary values are kow ANY potetial value withi the domai ca ow be foud. d K d K C

29 Boudary Elemet Method BEM THE END Book: The Boudary Elemet Method i Egieerig A.A.BECKER 9

PAPER : IIT-JAM 2010

PAPER : IIT-JAM 2010 MATHEMATICS-MA (CODE A) Q.-Q.5: Oly oe optio is correct for each questio. Each questio carries (+6) marks for correct aswer ad ( ) marks for icorrect aswer.. Which of the followig coditios does NOT esure

More information

Finite Difference Derivations for Spreadsheet Modeling John C. Walton Modified: November 15, 2007 jcw

Finite Difference Derivations for Spreadsheet Modeling John C. Walton Modified: November 15, 2007 jcw Fiite Differece Derivatios for Spreadsheet Modelig Joh C. Walto Modified: November 15, 2007 jcw Figure 1. Suset with 11 swas o Little Platte Lake, Michiga. Page 1 Modificatio Date: November 15, 2007 Review

More information

Castiel, Supernatural, Season 6, Episode 18

Castiel, Supernatural, Season 6, Episode 18 13 Differetial Equatios the aswer to your questio ca best be epressed as a series of partial differetial equatios... Castiel, Superatural, Seaso 6, Episode 18 A differetial equatio is a mathematical equatio

More information

A numerical Technique Finite Volume Method for Solving Diffusion 2D Problem

A numerical Technique Finite Volume Method for Solving Diffusion 2D Problem The Iteratioal Joural Of Egieerig d Sciece (IJES) Volume 4 Issue 10 Pages PP -35-41 2015 ISSN (e): 2319 1813 ISSN (p): 2319 1805 umerical Techique Fiite Volume Method for Solvig Diffusio 2D Problem 1 Mohammed

More information

Lecture 4 Conformal Mapping and Green s Theorem. 1. Let s try to solve the following problem by separation of variables

Lecture 4 Conformal Mapping and Green s Theorem. 1. Let s try to solve the following problem by separation of variables Lecture 4 Coformal Mappig ad Gree s Theorem Today s topics. Solvig electrostatic problems cotiued. Why separatio of variables does t always work 3. Coformal mappig 4. Gree s theorem The failure of separatio

More information

Math 257: Finite difference methods

Math 257: Finite difference methods Math 257: Fiite differece methods 1 Fiite Differeces Remember the defiitio of a derivative f f(x + ) f(x) (x) = lim 0 Also recall Taylor s formula: (1) f(x + ) = f(x) + f (x) + 2 f (x) + 3 f (3) (x) +...

More information

Analytic Continuation

Analytic Continuation Aalytic Cotiuatio The stadard example of this is give by Example Let h (z) = 1 + z + z 2 + z 3 +... kow to coverge oly for z < 1. I fact h (z) = 1/ (1 z) for such z. Yet H (z) = 1/ (1 z) is defied for

More information

Math 220B Final Exam Solutions March 18, 2002

Math 220B Final Exam Solutions March 18, 2002 Math 0B Fial Exam Solutios March 18, 00 1. (1 poits) (a) (6 poits) Fid the Gree s fuctio for the tilted half-plae {(x 1, x ) R : x 1 + x > 0}. For x (x 1, x ), y (y 1, y ), express your Gree s fuctio G(x,

More information

1 1 2 = show that: over variables x and y. [2 marks] Write down necessary conditions involving first and second-order partial derivatives for ( x0, y

1 1 2 = show that: over variables x and y. [2 marks] Write down necessary conditions involving first and second-order partial derivatives for ( x0, y Questio (a) A square matrix A= A is called positive defiite if the quadratic form waw > 0 for every o-zero vector w [Note: Here (.) deotes the traspose of a matrix or a vector]. Let 0 A = 0 = show that:

More information

Salmon: Lectures on partial differential equations. 3. First-order linear equations as the limiting case of second-order equations

Salmon: Lectures on partial differential equations. 3. First-order linear equations as the limiting case of second-order equations 3. First-order liear equatios as the limitig case of secod-order equatios We cosider the advectio-diffusio equatio (1) v = 2 o a bouded domai, with boudary coditios of prescribed. The coefficiets ( ) (2)

More information

The picture in figure 1.1 helps us to see that the area represents the distance traveled. Figure 1: Area represents distance travelled

The picture in figure 1.1 helps us to see that the area represents the distance traveled. Figure 1: Area represents distance travelled 1 Lecture : Area Area ad distace traveled Approximatig area by rectagles Summatio The area uder a parabola 1.1 Area ad distace Suppose we have the followig iformatio about the velocity of a particle, how

More information

Assignment 1 : Real Numbers, Sequences. for n 1. Show that (x n ) converges. Further, by observing that x n+2 + x n+1

Assignment 1 : Real Numbers, Sequences. for n 1. Show that (x n ) converges. Further, by observing that x n+2 + x n+1 Assigmet : Real Numbers, Sequeces. Let A be a o-empty subset of R ad α R. Show that α = supa if ad oly if α is ot a upper boud of A but α + is a upper boud of A for every N. 2. Let y (, ) ad x (, ). Evaluate

More information

Appendix A. Nabla and Friends

Appendix A. Nabla and Friends Appedix A Nabla ad Frieds A1 Notatio for Derivatives The partial derivative u u(x + he i ) u(x) (x) = lim x i h h of a scalar fuctio u : R R is writte i short otatio as Similarly we have for the higher

More information

We will conclude the chapter with the study a few methods and techniques which are useful

We will conclude the chapter with the study a few methods and techniques which are useful Chapter : Coordiate geometry: I this chapter we will lear about the mai priciples of graphig i a dimesioal (D) Cartesia system of coordiates. We will focus o drawig lies ad the characteristics of the graphs

More information

Rouché s Theorem 7/13/2011

Rouché s Theorem 7/13/2011 7/13/211 Rouché s Theorem Prelimiary Iformatio Complex Aalysis I complex aalysis, variables have a real part (x) ad a imagiary part (iy), where i = 1 z = x + iy Fuctios of complex variables also have a

More information

Chapter 10: Power Series

Chapter 10: Power Series Chapter : Power Series 57 Chapter Overview: Power Series The reaso series are part of a Calculus course is that there are fuctios which caot be itegrated. All power series, though, ca be itegrated because

More information

MATH 10550, EXAM 3 SOLUTIONS

MATH 10550, EXAM 3 SOLUTIONS MATH 155, EXAM 3 SOLUTIONS 1. I fidig a approximate solutio to the equatio x 3 +x 4 = usig Newto s method with iitial approximatio x 1 = 1, what is x? Solutio. Recall that x +1 = x f(x ) f (x ). Hece,

More information

PRELIM PROBLEM SOLUTIONS

PRELIM PROBLEM SOLUTIONS PRELIM PROBLEM SOLUTIONS THE GRAD STUDENTS + KEN Cotets. Complex Aalysis Practice Problems 2. 2. Real Aalysis Practice Problems 2. 4 3. Algebra Practice Problems 2. 8. Complex Aalysis Practice Problems

More information

A COMPLEX VARIABLE BOUNDARY ELEMENT METHOD FOR AN ELLIPTIC PARTIAL DIFFERENTIAL EQUATION WITH VARIABLE COEFFICIENTS

A COMPLEX VARIABLE BOUNDARY ELEMENT METHOD FOR AN ELLIPTIC PARTIAL DIFFERENTIAL EQUATION WITH VARIABLE COEFFICIENTS A COMPLEX VARIABLE BOUNDARY ELEMENT METHOD FOR AN ELLIPTIC PARTIAL DIFFERENTIAL EQUATION WITH VARIABLE COEFFICIENTS Y. S. PARK a adw.t.ang b Computatioal ad Mathematical Scieces Group Faculty of Iformatio

More information

Now we are looking to find a volume of solid S that lies below a surface z = f(x,y) and R= ab, cd,,[a,b] is the interval over

Now we are looking to find a volume of solid S that lies below a surface z = f(x,y) and R= ab, cd,,[a,b] is the interval over Multiple Itegratio Double Itegrals, Volume, ad Iterated Itegrals I sigle variable calculus we looked to fid the area uder a curve f(x) bouded by the x- axis over some iterval usig summatios the that led

More information

Numerical Conformal Mapping via a Fredholm Integral Equation using Fourier Method ABSTRACT INTRODUCTION

Numerical Conformal Mapping via a Fredholm Integral Equation using Fourier Method ABSTRACT INTRODUCTION alaysia Joural of athematical Scieces 3(1): 83-93 (9) umerical Coformal appig via a Fredholm Itegral Equatio usig Fourier ethod 1 Ali Hassa ohamed urid ad Teh Yua Yig 1, Departmet of athematics, Faculty

More information

A Lattice Green Function Introduction. Abstract

A Lattice Green Function Introduction. Abstract August 5, 25 A Lattice Gree Fuctio Itroductio Stefa Hollos Exstrom Laboratories LLC, 662 Nelso Park Dr, Logmot, Colorado 853, USA Abstract We preset a itroductio to lattice Gree fuctios. Electroic address:

More information

Lecture 7: Polar representation of complex numbers

Lecture 7: Polar representation of complex numbers Lecture 7: Polar represetatio of comple umbers See FLAP Module M3.1 Sectio.7 ad M3. Sectios 1 ad. 7.1 The Argad diagram I two dimesioal Cartesia coordiates (,), we are used to plottig the fuctio ( ) with

More information

NBHM QUESTION 2007 Section 1 : Algebra Q1. Let G be a group of order n. Which of the following conditions imply that G is abelian?

NBHM QUESTION 2007 Section 1 : Algebra Q1. Let G be a group of order n. Which of the following conditions imply that G is abelian? NBHM QUESTION 7 NBHM QUESTION 7 NBHM QUESTION 7 Sectio : Algebra Q Let G be a group of order Which of the followig coditios imply that G is abelia? 5 36 Q Which of the followig subgroups are ecesarily

More information

Math 5C Discussion Problems 2

Math 5C Discussion Problems 2 Math iscussio Problems Path Idepedece. Let be the striaght-lie path i R from the origi to (3, ). efie f(x, y) = xye xy. (a) Evaluate f dr. (b) Evaluate ((, 0) + f) dr. (c) Evaluate ((y, 0) + f) dr.. Let

More information

Streamfunction-Vorticity Formulation

Streamfunction-Vorticity Formulation Streamfuctio-Vorticity Formulatio A. Salih Departmet of Aerospace Egieerig Idia Istitute of Space Sciece ad Techology, Thiruvaathapuram March 2013 The streamfuctio-vorticity formulatio was amog the first

More information

Chimica Inorganica 3

Chimica Inorganica 3 himica Iorgaica Irreducible Represetatios ad haracter Tables Rather tha usig geometrical operatios, it is ofte much more coveiet to employ a ew set of group elemets which are matrices ad to make the rule

More information

MTH Assignment 1 : Real Numbers, Sequences

MTH Assignment 1 : Real Numbers, Sequences MTH -26 Assigmet : Real Numbers, Sequeces. Fid the supremum of the set { m m+ : N, m Z}. 2. Let A be a o-empty subset of R ad α R. Show that α = supa if ad oly if α is ot a upper boud of A but α + is a

More information

Math 5C Discussion Problems 2 Selected Solutions

Math 5C Discussion Problems 2 Selected Solutions Math 5 iscussio Problems 2 elected olutios Path Idepedece. Let be the striaght-lie path i 2 from the origi to (3, ). efie f(x, y) = xye xy. (a) Evaluate f dr. olutio. (b) Evaluate olutio. (c) Evaluate

More information

Math 113 Exam 3 Practice

Math 113 Exam 3 Practice Math Exam Practice Exam 4 will cover.-., 0. ad 0.. Note that eve though. was tested i exam, questios from that sectios may also be o this exam. For practice problems o., refer to the last review. This

More information

1 Adiabatic and diabatic representations

1 Adiabatic and diabatic representations 1 Adiabatic ad diabatic represetatios 1.1 Bor-Oppeheimer approximatio The time-idepedet Schrödiger equatio for both electroic ad uclear degrees of freedom is Ĥ Ψ(r, R) = E Ψ(r, R), (1) where the full molecular

More information

FINALTERM EXAMINATION Fall 9 Calculus & Aalytical Geometry-I Questio No: ( Mars: ) - Please choose oe Let f ( x) is a fuctio such that as x approaches a real umber a, either from left or right-had-side,

More information

Poincaré Problem for Nonlinear Elliptic Equations of Second Order in Unbounded Domains

Poincaré Problem for Nonlinear Elliptic Equations of Second Order in Unbounded Domains Advaces i Pure Mathematics 23 3 72-77 http://dxdoiorg/4236/apm233a24 Published Olie Jauary 23 (http://wwwscirporg/oural/apm) Poicaré Problem for Noliear Elliptic Equatios of Secod Order i Ubouded Domais

More information

Numerical Methods in Geophysics: Implicit Methods

Numerical Methods in Geophysics: Implicit Methods Numerical Methods i Geophysics: What is a implicit scheme? Explicit vs. implicit scheme for Newtoia oolig rak-nicholso Scheme (mixed explicit-implicit Explicit vs. implicit for the diffusio equatio Relaxatio

More information

Problem Cosider the curve give parametrically as x = si t ad y = + cos t for» t» ß: (a) Describe the path this traverses: Where does it start (whe t =

Problem Cosider the curve give parametrically as x = si t ad y = + cos t for» t» ß: (a) Describe the path this traverses: Where does it start (whe t = Mathematics Summer Wilso Fial Exam August 8, ANSWERS Problem 1 (a) Fid the solutio to y +x y = e x x that satisfies y() = 5 : This is already i the form we used for a first order liear differetial equatio,

More information

Brief Review of Functions of Several Variables

Brief Review of Functions of Several Variables Brief Review of Fuctios of Several Variables Differetiatio Differetiatio Recall, a fuctio f : R R is differetiable at x R if ( ) ( ) lim f x f x 0 exists df ( x) Whe this limit exists we call it or f(

More information

THE NUMERICAL SOLUTION OF THE NEWTONIAN FLUIDS FLOW DUE TO A STRETCHING CYLINDER BY SOR ITERATIVE PROCEDURE ABSTRACT

THE NUMERICAL SOLUTION OF THE NEWTONIAN FLUIDS FLOW DUE TO A STRETCHING CYLINDER BY SOR ITERATIVE PROCEDURE ABSTRACT Europea Joural of Egieerig ad Techology Vol. 3 No., 5 ISSN 56-586 THE NUMERICAL SOLUTION OF THE NEWTONIAN FLUIDS FLOW DUE TO A STRETCHING CYLINDER BY SOR ITERATIVE PROCEDURE Atif Nazir, Tahir Mahmood ad

More information

The Stokes Theorem. (Sect. 16.7) The curl of a vector field in space

The Stokes Theorem. (Sect. 16.7) The curl of a vector field in space The tokes Theorem. (ect. 6.7) The curl of a vector field i space. The curl of coservative fields. tokes Theorem i space. Idea of the proof of tokes Theorem. The curl of a vector field i space Defiitio

More information

Fluid Physics 8.292J/12.330J % (1)

Fluid Physics 8.292J/12.330J % (1) Fluid Physics 89J/133J Problem Set 5 Solutios 1 Cosider the flow of a Euler fluid i the x directio give by for y > d U = U y 1 d for y d U + y 1 d for y < This flow does ot vary i x or i z Determie the

More information

SECTION 2 Electrostatics

SECTION 2 Electrostatics SECTION Electrostatics This sectio, based o Chapter of Griffiths, covers effects of electric fields ad forces i static (timeidepedet) situatios. The topics are: Electric field Gauss s Law Electric potetial

More information

9.3 Power Series: Taylor & Maclaurin Series

9.3 Power Series: Taylor & Maclaurin Series 9.3 Power Series: Taylor & Maclauri Series If is a variable, the a ifiite series of the form 0 is called a power series (cetered at 0 ). a a a a a 0 1 0 is a power series cetered at a c a a c a c a c 0

More information

Estimation for Complete Data

Estimation for Complete Data Estimatio for Complete Data complete data: there is o loss of iformatio durig study. complete idividual complete data= grouped data A complete idividual data is the oe i which the complete iformatio of

More information

HOMEWORK #10 SOLUTIONS

HOMEWORK #10 SOLUTIONS Math 33 - Aalysis I Sprig 29 HOMEWORK # SOLUTIONS () Prove that the fuctio f(x) = x 3 is (Riema) itegrable o [, ] ad show that x 3 dx = 4. (Without usig formulae for itegratio that you leart i previous

More information

Teaching Mathematics Concepts via Computer Algebra Systems

Teaching Mathematics Concepts via Computer Algebra Systems Iteratioal Joural of Mathematics ad Statistics Ivetio (IJMSI) E-ISSN: 4767 P-ISSN: - 4759 Volume 4 Issue 7 September. 6 PP-- Teachig Mathematics Cocepts via Computer Algebra Systems Osama Ajami Rashaw,

More information

( ) (( ) ) ANSWERS TO EXERCISES IN APPENDIX B. Section B.1 VECTORS AND SETS. Exercise B.1-1: Convex sets. are convex, , hence. and. (a) Let.

( ) (( ) ) ANSWERS TO EXERCISES IN APPENDIX B. Section B.1 VECTORS AND SETS. Exercise B.1-1: Convex sets. are convex, , hence. and. (a) Let. Joh Riley 8 Jue 03 ANSWERS TO EXERCISES IN APPENDIX B Sectio B VECTORS AND SETS Exercise B-: Covex sets (a) Let 0 x, x X, X, hece 0 x, x X ad 0 x, x X Sice X ad X are covex, x X ad x X The x X X, which

More information

Complex Analysis Spring 2001 Homework I Solution

Complex Analysis Spring 2001 Homework I Solution Complex Aalysis Sprig 2001 Homework I Solutio 1. Coway, Chapter 1, sectio 3, problem 3. Describe the set of poits satisfyig the equatio z a z + a = 2c, where c > 0 ad a R. To begi, we see from the triagle

More information

1 Cabin. Professor: What is. Student: ln Cabin oh Log Cabin! Professor: No. Log Cabin + C = A Houseboat!

1 Cabin. Professor: What is. Student: ln Cabin oh Log Cabin! Professor: No. Log Cabin + C = A Houseboat! MATH 4 Sprig 0 Exam # Tuesday March st Sectios: Sectios 6.-6.6; 6.8; 7.-7.4 Name: Score: = 00 Istructios:. You will have a total of hour ad 50 miutes to complete this exam.. A No-Graphig Calculator may

More information

6.3 Testing Series With Positive Terms

6.3 Testing Series With Positive Terms 6.3. TESTING SERIES WITH POSITIVE TERMS 307 6.3 Testig Series With Positive Terms 6.3. Review of what is kow up to ow I theory, testig a series a i for covergece amouts to fidig the i= sequece of partial

More information

a b c d e f g h Supplementary Information

a b c d e f g h Supplementary Information Supplemetary Iformatio a b c d e f g h Supplemetary Figure S STM images show that Dark patters are frequetly preset ad ted to accumulate. (a) mv, pa, m ; (b) mv, pa, m ; (c) mv, pa, m ; (d) mv, pa, m ;

More information

Fundamental Concepts: Surfaces and Curves

Fundamental Concepts: Surfaces and Curves UNDAMENTAL CONCEPTS: SURACES AND CURVES CHAPTER udametal Cocepts: Surfaces ad Curves. INTRODUCTION This chapter describes two geometrical objects, vi., surfaces ad curves because the pla a ver importat

More information

j=1 dz Res(f, z j ) = 1 d k 1 dz k 1 (z c)k f(z) Res(f, c) = lim z c (k 1)! Res g, c = f(c) g (c)

j=1 dz Res(f, z j ) = 1 d k 1 dz k 1 (z c)k f(z) Res(f, c) = lim z c (k 1)! Res g, c = f(c) g (c) Problem. Compute the itegrals C r d for Z, where C r = ad r >. Recall that C r has the couter-clockwise orietatio. Solutio: We will use the idue Theorem to solve this oe. We could istead use other (perhaps

More information

Lecture 7: Density Estimation: k-nearest Neighbor and Basis Approach

Lecture 7: Density Estimation: k-nearest Neighbor and Basis Approach STAT 425: Itroductio to Noparametric Statistics Witer 28 Lecture 7: Desity Estimatio: k-nearest Neighbor ad Basis Approach Istructor: Ye-Chi Che Referece: Sectio 8.4 of All of Noparametric Statistics.

More information

An efficient time integration method for extra-large eddy simulations

An efficient time integration method for extra-large eddy simulations A efficiet time itegratio method for extra-large eddy simulatios M.A. Scheibeler Departmet of Mathematics Master s Thesis A efficiet time itegratio method for extra-large eddy simulatios M.A. Scheibeler

More information

: Transforms and Partial Differential Equations

: Transforms and Partial Differential Equations Trasforms ad Partial Differetial Equatios 018 SUBJECT NAME : Trasforms ad Partial Differetial Equatios SUBJECT CODE : MA 6351 MATERIAL NAME : Part A questios REGULATION : R013 WEBSITE : wwwharigaeshcom

More information

ECE 669 Parallel Computer Architecture

ECE 669 Parallel Computer Architecture ECE 669 Parallel Computer Architecture Lecture 4 Parallel Applicatios Outlie Motivatig Problems (applicatio case studies) Classifyig problems Parallelizig applicatios Examiig tradeoffs Uderstadig commuicatio

More information

MIDTERM 3 CALCULUS 2. Monday, December 3, :15 PM to 6:45 PM. Name PRACTICE EXAM SOLUTIONS

MIDTERM 3 CALCULUS 2. Monday, December 3, :15 PM to 6:45 PM. Name PRACTICE EXAM SOLUTIONS MIDTERM 3 CALCULUS MATH 300 FALL 08 Moday, December 3, 08 5:5 PM to 6:45 PM Name PRACTICE EXAM S Please aswer all of the questios, ad show your work. You must explai your aswers to get credit. You will

More information

Apply change-of-basis formula to rewrite x as a linear combination of eigenvectors v j.

Apply change-of-basis formula to rewrite x as a linear combination of eigenvectors v j. Eigevalue-Eigevector Istructor: Nam Su Wag eigemcd Ay vector i real Euclidea space of dimesio ca be uiquely epressed as a liear combiatio of liearly idepedet vectors (ie, basis) g j, j,,, α g α g α g α

More information

Lecture 2: Finite Difference Methods in Heat Transfer

Lecture 2: Finite Difference Methods in Heat Transfer Lecture 2: Fiite Differece Methods i Heat Trasfer V.Vuorie Aalto Uiversity School of Egieerig Heat ad Mass Trasfer Course, Autum 2016 November 1 st 2017, Otaiemi ville.vuorie@aalto.fi Overview Part 1 (

More information

ME NUMERICAL METHODS Fall 2007

ME NUMERICAL METHODS Fall 2007 ME - 310 NUMERICAL METHODS Fall 2007 Group 02 Istructor: Prof. Dr. Eres Söylemez (Rm C205, email:eres@metu.edu.tr ) Class Hours ad Room: Moday 13:40-15:30 Rm: B101 Wedesday 12:40-13:30 Rm: B103 Course

More information

Definitions and Theorems. where x are the decision variables. c, b, and a are constant coefficients.

Definitions and Theorems. where x are the decision variables. c, b, and a are constant coefficients. Defiitios ad Theorems Remember the scalar form of the liear programmig problem, Miimize, Subject to, f(x) = c i x i a 1i x i = b 1 a mi x i = b m x i 0 i = 1,2,, where x are the decisio variables. c, b,

More information

Regn. No. North Delhi : 33-35, Mall Road, G.T.B. Nagar (Opp. Metro Gate No. 3), Delhi-09, Ph: ,

Regn. No. North Delhi : 33-35, Mall Road, G.T.B. Nagar (Opp. Metro Gate No. 3), Delhi-09, Ph: , . Sectio-A cotais 30 Multiple Choice Questios (MCQ). Each questio has 4 choices (a), (b), (c) ad (d), for its aswer, out of which ONLY ONE is correct. From Q. to Q.0 carries Marks ad Q. to Q.30 carries

More information

MATH 2411 Spring 2011 Practice Exam #1 Tuesday, March 1 st Sections: Sections ; 6.8; Instructions:

MATH 2411 Spring 2011 Practice Exam #1 Tuesday, March 1 st Sections: Sections ; 6.8; Instructions: MATH 411 Sprig 011 Practice Exam #1 Tuesday, March 1 st Sectios: Sectios 6.1-6.6; 6.8; 7.1-7.4 Name: Score: = 100 Istructios: 1. You will have a total of 1 hour ad 50 miutes to complete this exam.. A No-Graphig

More information

Math 113 Exam 3 Practice

Math 113 Exam 3 Practice Math Exam Practice Exam will cover.-.9. This sheet has three sectios. The first sectio will remid you about techiques ad formulas that you should kow. The secod gives a umber of practice questios for you

More information

Z - Transform. It offers the techniques for digital filter design and frequency analysis of digital signals.

Z - Transform. It offers the techniques for digital filter design and frequency analysis of digital signals. Z - Trasform The -trasform is a very importat tool i describig ad aalyig digital systems. It offers the techiques for digital filter desig ad frequecy aalysis of digital sigals. Defiitio of -trasform:

More information

We are mainly going to be concerned with power series in x, such as. (x)} converges - that is, lims N n

We are mainly going to be concerned with power series in x, such as. (x)} converges - that is, lims N n Review of Power Series, Power Series Solutios A power series i x - a is a ifiite series of the form c (x a) =c +c (x a)+(x a) +... We also call this a power series cetered at a. Ex. (x+) is cetered at

More information

Section 1.4. Power Series

Section 1.4. Power Series Sectio.4. Power Series De itio. The fuctio de ed by f (x) (x a) () c 0 + c (x a) + c 2 (x a) 2 + c (x a) + ::: is called a power series cetered at x a with coe ciet sequece f g :The domai of this fuctio

More information

Analysis of composites with multiple rigid-line reinforcements by the BEM

Analysis of composites with multiple rigid-line reinforcements by the BEM Aalysis of composites with multiple rigid-lie reiforcemets by the BEM Piotr Fedeliski* Departmet of Stregth of Materials ad Computatioal Mechaics, Silesia Uiversity of Techology ul. Koarskiego 18A, 44-100

More information

3. Z Transform. Recall that the Fourier transform (FT) of a DT signal xn [ ] is ( ) [ ] = In order for the FT to exist in the finite magnitude sense,

3. Z Transform. Recall that the Fourier transform (FT) of a DT signal xn [ ] is ( ) [ ] = In order for the FT to exist in the finite magnitude sense, 3. Z Trasform Referece: Etire Chapter 3 of text. Recall that the Fourier trasform (FT) of a DT sigal x [ ] is ω ( ) [ ] X e = j jω k = xe I order for the FT to exist i the fiite magitude sese, S = x [

More information

MATHEMATICS. 61. The differential equation representing the family of curves where c is a positive parameter, is of

MATHEMATICS. 61. The differential equation representing the family of curves where c is a positive parameter, is of MATHEMATICS 6 The differetial equatio represetig the family of curves where c is a positive parameter, is of Order Order Degree (d) Degree (a,c) Give curve is y c ( c) Differetiate wrt, y c c y Hece differetial

More information

Lecture 8: Solving the Heat, Laplace and Wave equations using finite difference methods

Lecture 8: Solving the Heat, Laplace and Wave equations using finite difference methods Itroductory lecture otes o Partial Differetial Equatios - c Athoy Peirce. Not to be copied, used, or revised without explicit writte permissio from the copyright ower. 1 Lecture 8: Solvig the Heat, Laplace

More information

ES.182A Topic 40 Notes Jeremy Orloff

ES.182A Topic 40 Notes Jeremy Orloff ES.182A opic 4 Notes Jeremy Orloff 4 Flux: ormal form of Gree s theorem Gree s theorem i flux form is formally equivalet to our previous versio where the lie itegral was iterpreted as work. Here we will

More information

Stochastic Simulation

Stochastic Simulation Stochastic Simulatio 1 Itroductio Readig Assigmet: Read Chapter 1 of text. We shall itroduce may of the key issues to be discussed i this course via a couple of model problems. Model Problem 1 (Jackso

More information

G-2 Applied Computational Aerodynamics. G-1 Program THINFOIL

G-2 Applied Computational Aerodynamics. G-1 Program THINFOIL G-2 Applied Computatioal Aerodyamics G-1 Program THINFOIL THINFOIL solves Laplace s Equatio by fiite differeces usig a variety of iteratio methods. It was writte by Valery Razgoyaev. The iteratio optios

More information

Math 105: Review for Final Exam, Part II - SOLUTIONS

Math 105: Review for Final Exam, Part II - SOLUTIONS Math 5: Review for Fial Exam, Part II - SOLUTIONS. Cosider the fuctio f(x) = x 3 lx o the iterval [/e, e ]. (a) Fid the x- ad y-coordiates of ay ad all local extrema ad classify each as a local maximum

More information

TEACHER CERTIFICATION STUDY GUIDE

TEACHER CERTIFICATION STUDY GUIDE COMPETENCY 1. ALGEBRA SKILL 1.1 1.1a. ALGEBRAIC STRUCTURES Kow why the real ad complex umbers are each a field, ad that particular rigs are ot fields (e.g., itegers, polyomial rigs, matrix rigs) Algebra

More information

Lecture 2: Monte Carlo Simulation

Lecture 2: Monte Carlo Simulation STAT/Q SCI 43: Itroductio to Resamplig ethods Sprig 27 Istructor: Ye-Chi Che Lecture 2: ote Carlo Simulatio 2 ote Carlo Itegratio Assume we wat to evaluate the followig itegratio: e x3 dx What ca we do?

More information

PARTIAL DIFFERENTIAL EQUATIONS SEPARATION OF VARIABLES

PARTIAL DIFFERENTIAL EQUATIONS SEPARATION OF VARIABLES Diola Bagayoko (0 PARTAL DFFERENTAL EQUATONS SEPARATON OF ARABLES. troductio As discussed i previous lectures, partial differetial equatios arise whe the depedet variale, i.e., the fuctio, varies with

More information

MATH 1A FINAL (7:00 PM VERSION) SOLUTION. (Last edited December 25, 2013 at 9:14pm.)

MATH 1A FINAL (7:00 PM VERSION) SOLUTION. (Last edited December 25, 2013 at 9:14pm.) MATH A FINAL (7: PM VERSION) SOLUTION (Last edited December 5, 3 at 9:4pm.) Problem. (i) Give the precise defiitio of the defiite itegral usig Riema sums. (ii) Write a epressio for the defiite itegral

More information

THE SOLUTION OF NONLINEAR EQUATIONS f( x ) = 0.

THE SOLUTION OF NONLINEAR EQUATIONS f( x ) = 0. THE SOLUTION OF NONLINEAR EQUATIONS f( ) = 0. Noliear Equatio Solvers Bracketig. Graphical. Aalytical Ope Methods Bisectio False Positio (Regula-Falsi) Fied poit iteratio Newto Raphso Secat The root of

More information

CO-LOCATED DIFFUSE APPROXIMATION METHOD FOR TWO DIMENSIONAL INCOMPRESSIBLE CHANNEL FLOWS

CO-LOCATED DIFFUSE APPROXIMATION METHOD FOR TWO DIMENSIONAL INCOMPRESSIBLE CHANNEL FLOWS CO-LOCATED DIFFUSE APPROXIMATION METHOD FOR TWO DIMENSIONAL INCOMPRESSIBLE CHANNEL FLOWS C.PRAX ad H.SADAT Laboratoire d'etudes Thermiques,URA CNRS 403 40, Aveue du Recteur Pieau 86022 Poitiers Cedex,

More information

Areas and Distances. We can easily find areas of certain geometric figures using well-known formulas:

Areas and Distances. We can easily find areas of certain geometric figures using well-known formulas: Areas ad Distaces We ca easily fid areas of certai geometric figures usig well-kow formulas: However, it is t easy to fid the area of a regio with curved sides: METHOD: To evaluate the area of the regio

More information

Statistical Pattern Recognition

Statistical Pattern Recognition Statistical Patter Recogitio Classificatio: No-Parametric Modelig Hamid R. Rabiee Jafar Muhammadi Sprig 2014 http://ce.sharif.edu/courses/92-93/2/ce725-2/ Ageda Parametric Modelig No-Parametric Modelig

More information

METHOD OF FUNDAMENTAL SOLUTIONS FOR HELMHOLTZ EIGENVALUE PROBLEMS IN ELLIPTICAL DOMAINS

METHOD OF FUNDAMENTAL SOLUTIONS FOR HELMHOLTZ EIGENVALUE PROBLEMS IN ELLIPTICAL DOMAINS Please cite this article as: Staisław Kula, Method of fudametal solutios for Helmholtz eigevalue problems i elliptical domais, Scietific Research of the Istitute of Mathematics ad Computer Sciece, 009,

More information

1.3 Convergence Theorems of Fourier Series. k k k k. N N k 1. With this in mind, we state (without proof) the convergence of Fourier series.

1.3 Convergence Theorems of Fourier Series. k k k k. N N k 1. With this in mind, we state (without proof) the convergence of Fourier series. .3 Covergece Theorems of Fourier Series I this sectio, we preset the covergece of Fourier series. A ifiite sum is, by defiitio, a limit of partial sums, that is, a cos( kx) b si( kx) lim a cos( kx) b si(

More information

1 Approximating Integrals using Taylor Polynomials

1 Approximating Integrals using Taylor Polynomials Seughee Ye Ma 8: Week 7 Nov Week 7 Summary This week, we will lear how we ca approximate itegrals usig Taylor series ad umerical methods. Topics Page Approximatig Itegrals usig Taylor Polyomials. Defiitios................................................

More information

For example suppose we divide the interval [0,2] into 5 equal subintervals of length

For example suppose we divide the interval [0,2] into 5 equal subintervals of length Math 1206 Calculus Sec 1: Estimatig with Fiite Sums Abbreviatios: wrt with respect to! for all! there exists! therefore Def defiitio Th m Theorem sol solutio! perpedicular iff or! if ad oly if pt poit

More information

SOLVED EXAMPLES

SOLVED EXAMPLES Prelimiaries Chapter PELIMINAIES Cocept of Divisibility: A o-zero iteger t is said to be a divisor of a iteger s if there is a iteger u such that s tu I this case we write t s (i) 6 as ca be writte as

More information

EDEXCEL STUDENT CONFERENCE 2006 A2 MATHEMATICS STUDENT NOTES

EDEXCEL STUDENT CONFERENCE 2006 A2 MATHEMATICS STUDENT NOTES EDEXCEL STUDENT CONFERENCE 006 A MATHEMATICS STUDENT NOTES South: Thursday 3rd March 006, Lodo EXAMINATION HINTS Before the eamiatio Obtai a copy of the formulae book ad use it! Write a list of ad LEARN

More information

10-701/ Machine Learning Mid-term Exam Solution

10-701/ Machine Learning Mid-term Exam Solution 0-70/5-78 Machie Learig Mid-term Exam Solutio Your Name: Your Adrew ID: True or False (Give oe setece explaatio) (20%). (F) For a cotiuous radom variable x ad its probability distributio fuctio p(x), it

More information

f(x) dx as we do. 2x dx x also diverges. Solution: We compute 2x dx lim

f(x) dx as we do. 2x dx x also diverges. Solution: We compute 2x dx lim Math 3, Sectio 2. (25 poits) Why we defie f(x) dx as we do. (a) Show that the improper itegral diverges. Hece the improper itegral x 2 + x 2 + b also diverges. Solutio: We compute x 2 + = lim b x 2 + =

More information

b i u x i U a i j u x i u x j

b i u x i U a i j u x i u x j M ath 5 2 7 Fall 2 0 0 9 L ecture 1 9 N ov. 1 6, 2 0 0 9 ) S ecod- Order Elliptic Equatios: Weak S olutios 1. Defiitios. I this ad the followig two lectures we will study the boudary value problem Here

More information

f(w) w z =R z a 0 a n a nz n Liouville s theorem, we see that Q is constant, which implies that P is constant, which is a contradiction.

f(w) w z =R z a 0 a n a nz n Liouville s theorem, we see that Q is constant, which implies that P is constant, which is a contradiction. Theorem 3.6.4. [Liouville s Theorem] Every bouded etire fuctio is costat. Proof. Let f be a etire fuctio. Suppose that there is M R such that M for ay z C. The for ay z C ad R > 0 f (z) f(w) 2πi (w z)

More information

Exact Solutions for a Class of Nonlinear Singular Two-Point Boundary Value Problems: The Decomposition Method

Exact Solutions for a Class of Nonlinear Singular Two-Point Boundary Value Problems: The Decomposition Method Exact Solutios for a Class of Noliear Sigular Two-Poit Boudary Value Problems: The Decompositio Method Abd Elhalim Ebaid Departmet of Mathematics, Faculty of Sciece, Tabuk Uiversity, P O Box 741, Tabuki

More information

Optimization of the Brownie Pan

Optimization of the Brownie Pan Optimizatio of the Browie Pa Briaa Oshiro bsoshiro@uw.edu Patrick Larso palarso@uw.edu February 4, 2013 Sujay Cauligi sujayc@uw.edu Abstract I this paper we address the effect that pa shape has o uiformity

More information

Vector Quantization: a Limiting Case of EM

Vector Quantization: a Limiting Case of EM . Itroductio & defiitios Assume that you are give a data set X = { x j }, j { 2,,, }, of d -dimesioal vectors. The vector quatizatio (VQ) problem requires that we fid a set of prototype vectors Z = { z

More information

L 5 & 6: RelHydro/Basel. f(x)= ( ) f( ) ( ) ( ) ( ) n! 1! 2! 3! If the TE of f(x)= sin(x) around x 0 is: sin(x) = x - 3! 5!

L 5 & 6: RelHydro/Basel. f(x)= ( ) f( ) ( ) ( ) ( ) n! 1! 2! 3! If the TE of f(x)= sin(x) around x 0 is: sin(x) = x - 3! 5! aylor epasio: Let ƒ() be a ifiitely differetiable real fuctio. At ay poit i the eighbourhood of =0, the fuctio ca be represeted as a power series of the followig form: X 0 f(a) f() ƒ() f()= ( ) f( ) (

More information

Machine Learning for Data Science (CS 4786)

Machine Learning for Data Science (CS 4786) Machie Learig for Data Sciece CS 4786) Lecture & 3: Pricipal Compoet Aalysis The text i black outlies high level ideas. The text i blue provides simple mathematical details to derive or get to the algorithm

More information

Failure Theories Des Mach Elem Mech. Eng. Department Chulalongkorn University

Failure Theories Des Mach Elem Mech. Eng. Department Chulalongkorn University Failure Theories Review stress trasformatio Failure theories for ductile materials Maimum-Shear-Stress Theor Distortio-Eerg Theor Coulomb-Mohr Theor Failure theories for brittle materials Maimum-Normal-Stress

More information

6 Integers Modulo n. integer k can be written as k = qn + r, with q,r, 0 r b. So any integer.

6 Integers Modulo n. integer k can be written as k = qn + r, with q,r, 0 r b. So any integer. 6 Itegers Modulo I Example 2.3(e), we have defied the cogruece of two itegers a,b with respect to a modulus. Let us recall that a b (mod ) meas a b. We have proved that cogruece is a equivalece relatio

More information

Ellipsoid Method for Linear Programming made simple

Ellipsoid Method for Linear Programming made simple Ellipsoid Method for Liear Programmig made simple Sajeev Saxea Dept. of Computer Sciece ad Egieerig, Idia Istitute of Techology, Kapur, INDIA-08 06 December 3, 07 Abstract I this paper, ellipsoid method

More information