Boundary Element Method for the Helmholtz Equation
|
|
- Hilary Wood
- 5 years ago
- Views:
Transcription
1 Journal of Mathematical Extension Vol. 9, No. 4, (205), ISSN: URL: Boundary Element Method for the Helmholtz Equation A. Mesforush Shahrood University Z. Mehraban Shahrood University Abstract. The boundary element methos is applied to the Helholtltz equation. To this end we discretized the Helmholtz equation over the boundary and conclude a system of equations. By applying the boundary condition we get a new system which gives the solution of the equation. AMS Subject Classification: 65L80; 35Axx Keywords and Phrases: Laplace equation, Helmholtz equation, boundary integral equation, discretization, boundary element method. Introduction Numerical approximation for solving a differential equation has a wide range of application in engineering and mathematics. Some of these methods are: finite difference methods, finite element methods and boundary element method. In these methods we usually discritize the equation and apply the governing equations on the suitable mesh. The most important thing in the numerical approximation is determination of the size of mesh such that the approximation solution has the suitable accuracy. Received: January 205; Accepted: May 205 Corresponding author 65
2 66 A. MESFORUSH AND Z. MEHRABAN The boundary element method is a powerful tool for numerical approximation for a differential equation. For the first time, in 872, Betti presented this method in elasticity theory []. In 903 Fredholm extended this method [2], but in 977 the boundary element method appears in some publications of Banerjee, Butterfield, Brebbia and Dominguez [3, 4]. In this paper we want to find an approximate solution for the Helmholtz equation by the boundary element method. This paper organized in three sections: in the first section, we introduce some preliminaries and also describe the structure of the boundary element method. In the second section, we applied boundary element method for Laplace equation which is a special case of helmholtz equation. Finally, in the third section we approximat the boundary element solution of Helmholtz equation. In boundary element method we use three important things. - The Gauss divergence theorem: If R d be a bounded domain and w = (w,, w d ) C () d be a vector-valued function, then for x R d we have [8] w dv = w n ds, where n = (n,, n d ) is the outward unit normal to boundary Γ. 2- The second Green identity: Let R d be a bounded domain and v, w C 2 (), then [8] (w 2 u u 2 w)dv = (w n u w n )ds. 3- The Dirac delta function: The Dirac delta function defined as [] {, X = X 0, δ(x X 0 ) = 0, X X 0, and has two following properties: δ(x X 0 ) dx =, Γ Γ f(x)δ(x X 0 ) dx = f(x 0 ).
3 BOUNDARY ELEMENT METHOD FOR The structure of the boundary element method can be decribed as follows. Consider the following boundary element value problem Au := Au xx + 2Bu xy + Cu yy + Du = f, in, u = ū, on Γ, q = n = q, on Γ 2, where the coefficients A, B, C, D are constant and Γ, Γ 2 are boundaries of domain, q is the derivative of potential function in direction of n and ū, q are given values of the flux and potential function on the boundary. In boundary element method we multiply the first equation of () by weight function w C 2 () and integrate over to get the following integral form (Au)wdv = fwdv. (2) by choosing the fundamental solution Aw = δ(x X 0 ) as the weight function and using the second Green identity, we can convert (2) to a continuous integral equation over the boundary. For applying the boundary element method we discretize the problem over the boundary of. To this end we consider = N j= j and conclude the system of equations HU = GQ. If we apply the boundary conditions, we get a new system of equations as AX = B, which gives the flux and potential u on the boundary. In this paper we want to apply the bounday element method for the Helmholtz equation. In () if we consider A =, B = 0, C =, D = k 2 and by considering 2 u = u xx + u yy, we get the Helmholtz equation as where k is the wave number. Au := 2 u + k 2 u = 0, in R 2, u = ū, on Γ, q = n = q, on Γ 2, () (3)
4 68 A. MESFORUSH AND Z. MEHRABAN 2. Boundary Element Method for the Laplace Problem In this section we apply the boundary element method to the Laplace equation which is the special case of Helmholtz equation. In equation (3) if we set k = 0, the following laplace equation is resulted: 2 u = 0, in R 2, u = ū, on Γ, q = (4) n = q, on Γ 2. Consider weight function w such that it satisfies 2 w + δ(x X 0 ) = 0, (5) where X 0 = (ξ, η) is a singular point and δ is the Dirac delta function. Considering the polar form 2 w = w (r ), where r = X r r r X 0 = (x ξ) 2 + (y η) 2. In this case for r > 0, δ(x X 0 ) = 0, one obtain w (r r r r ) = 0, where A and B are arbitrary constants. Since we look for a particular solution, we may set B = 0. Which gives after integrating twice w = A ln r + B. By considering a neighborhood of X 0 and using the integral property for dirac delta function, we obtain A =. Hence, the fundamental 2π solution for the Laplac equation becomes [] w = ( ) 2π ln. r Also note that w n = w r r n = ( r n) = 2πr 2πr 2 (r xn x + r y n y ),
5 BOUNDARY ELEMENT METHOD FOR where r x = x ξ and r y = y η. Let X 0 and consider a neighborhood of X 0 with radius ε, namely ε and consider the new domain ε with boundary ε. In this domain by letting q = n and q = w the Green formula can be n written as lim ε (u 2 w w 2 u)dv = lim (uq wq)ds + lim (uq wq)ds, ε by using (5) and the definition of the Dirac delta and X 0 / ε, one obtain lim u 2 wdv = lim uδ(x X 0 )dv = 0, ε ε using the Laplac equation (4), result in lim ε (u 2 w w 2 u)dv = 0. On the other hand, [9] lim (uq wq)ds = lim (u u(x 0 ) + u(x 0 )) q ds ε ε lim ln ε lim ln ε lim ln ε 2π ε n ds 2π ε n ds 2π ε n ds 2π 2π = lim u(x = lim u(x 0 ) 0 2π dθ + lim ε ln ε 0 ) 2π 0 2π dθ + lim ε ln ε (X 0 ) 2π ln ( ) dθ 2π n 0 = u(x 0 ). = u(x 0 ). Therefore, the boundary integral equation in is: u(x 0 ) = (wq uq )ds, (X 0 ) n which means that for given u, q on, the values of u can be determinded for the interior point of the domain.
6 70 A. MESFORUSH AND Z. MEHRABAN Let X 0 and consider a neighborhood of X 0, namely i ε and consider the domain i ε with boundary i ε ( ε ). In this domain the Green formula can be written as, lim i ε (u 2 w w 2 u)dv = lim (uq ( wq)ds ε) + lim (uq wq)ds. i ε Since the governing equation is Laplace equation and the X 0 is out of the domain i ε,we have: lim (u 2 w w 2 u)dv = 0. ε On the other hand lim (uq wq)ds = lim (u u(x 0 ) + u(x 0 )) q ds i ε ε lim ln ε lim ln ε 2π ε n ds 2π ε n ds = lim u(x 0 ) i 2π dθ + lim ε ln ε = lim u(x 0 ) 2π i 2π dθ + lim ε ln ε (X 0 ) 2π dθ 2πε n ε 0 = C(X = C(X 0 )u(x 0 ), 0 )u(x 0 ), (X 0 ) n where C(X 0 ) is a geometry-dependent parameter. It varies between zero and unity and equals 2. If the point is on a smooth part of the boundary at edges, it is related to the angle of the joining surfaces and equals internal angle C(X 0 ) =. Also 2π lim (uq wq)ds = lim (uq ε wq)ds. So the boundary integral equation in is: C(α)u(X 0 ) + uq ds = wqds. 2π 0
7 BOUNDARY ELEMENT METHOD FOR... 7 To apply the boundary element method, the boundary is divided into small elements as = N j= j where N represents the number of elements that form the boundary. The mesh is generated only on the boundary and it is more flexible than in the case of finite element mesh. The integrals are then expressed for each element and summed up in order to evaluate the integral on the whole boundary. Considering constant element discretization and the source points situated on a smooth boundary, the boundary integral equation can be written as 2 u(x 0) + uq ds = wqds. N j= j N j= j where by definition Ĥ ij, i j, H ij = 2 + Ĥij, i = j, and we get w Ĥ ij = j n ds j, G ij = w ds j, j N H ij u j = j= N G ij q j. j= Here we know some of u j s and some of q j s, i.e.,there are N known quantities and N unknowns. Applying the boundary conditions to every midpoint of the elements, we obtain the equations in the matrix-vector form as HU = GQ. Rearranging this system with known quantities on one side and unknowns on the other, we get a linear system of equations AX = B where X is unknown values. Example 2.. Consider the Laplace equation in the rectangle plate with
8 72 A. MESFORUSH AND Z. MEHRABAN mixed boundary conditions as [7] 2 T = 0, 0 < x < 2, 0 < y <, (x, 0) = 0, T (2, y) = 0, n (x, ) = 0, T (0, y) =. n Using four constant element discretization with the collocation points at the center of the elements, H ij and G ij for i j, one obtain: H ij = q ds j = Xj+ j 2π X j r 2 (r xn x + r y n y )dx = l 4 e w k 4π rk 2 (r xk n xk + r yk n yk ), k= and G ij = wds j = j 2π X j = l 4 ( ) e w k ln. 4π r k k= Xj+ ln ( ) dx r also for i = j, we have H ij = 2, and G ij = π le 2 0 ln ( ) dx = l ( e ln r 2π ( 2 l e ) ) +, where l e = X j+ X j is the length of the element on which the integral is calculated, w k represent the weight functions from the Gauss integration scheme and r k is the distance from the i th source point to each integration points on the j th element. So H = ,
9 BOUNDARY ELEMENT METHOD FOR G = From these equations unknow values in system AX = B will be calculated as T = , T 3 = , q 2 = , q 4 = Consider these values of T in interior points, result in as where So P = (0.75, 0.5), P 2 = (, 0.5), P 3 = (.5, 0.5), T (i) = N G ik q k k= N H ik T k, k= w i G ik = w i ds k, H ik = Γ k Γ k n ds k, T P = , T P2 = , T P3 = 0.38, Therefore the exact values are T e P = 0.625, T e P 2 = 0.5, T e P 3 = Compared with the exact values, the solution shows quite a good approximation for the internal points not only for the boundary points. Figure, 2 Figure. The effect of increasing the number of elements on approximate solution of the Laplace equation
10 74 A. MESFORUSH AND Z. MEHRABAN Figure 2. The effect of increasing the number of elements on approximate solution of the Laplace equation (zoom) 3. Boundary Element Method for the Helmholtz Equation In this section we apply the boundary element method to the Helmholtz equation (3) for k = 0. Consider weight function w such that it satisfies 2 w + k 2 w + δ(x X 0 ) = 0, (6) where X 0 = (ξ, η) is a singular point. The solutions for the homogeneous case are the Hankel functions of the first and second kinds of order zero. So we write w(r) = AH (2) 0 (kr), where r = (x ξ) 2 + (y η) 2 and H (2) 0 is the Hankel function of the second kind of order zero and it is representable in terms of the 0th-order Bessel functions of the first and second kinds J 0, Y 0 as H (2) 0 (z) = J 0(z) iy 0 (z) = H () 0 (z). By using the small-argument approximation to the Hankel function as [2] 0 (kr) i 2 π log(γkr ), r 0, 2 H (2)
11 BOUNDARY ELEMENT METHOD FOR where γ =.78 and integrating (6) over a very small circle of radius ε centered at the origin, one obtain A [ + k 2 ][ i 2 π log(γkr )]dv =, 2 using the divergence theorem, the first term is converted to a line integral as i 2 2π ( log( γkr ) π 2 ) rdφ = 4i, the second term goes to zero, so 0 A = i 4, Therefore, we get the following fundamental solution for helmholtz equation as [2] Also we note that [0] q = w n = w r r n = 4i kh(2) w = 4i H(2) 0 (kr). (kr)( r n) = n xr x + n y r y 4ir kh (2) (kr), where r x = x ξ and r y = y η. To obtain the boundary integral equation, we integrate (3) over the domain with a weighting function w to get, ( 2 u + k 2 u)wdv = 0. The second Green identity, implies that w 2 udv = u 2 wdv + so (uq wq)ds, ( 2 w + k 2 w)udv + (uq wq)ds = 0. (7) By using the properties of Dirac delta gives ( 2 w + k 2 w)udv = u(x 0 ).
12 76 A. MESFORUSH AND Z. MEHRABAN So the boundary integral equation for all the interior points is u(x 0 ) + uq ds = wqds. The boundary integral equation in general case is [3] C(X 0 )u(x 0 ) + uq ds = wqds, where C(X 0 ) = α is a geometry-dependent parameter where α is internal angle. Assume that the boundary is discretized and represented by 2π N linear constant elements as = N j= j and the source points situated on a smooth boundary (α = π), so the boundary integral equation can be written as 2 u(x 0) + u w N j= n ds = w j N j= n ds. j Therefore where and N H ij u j = j= N G ij q j, j= Ĥ ij, i j, H ij = 2 + Ĥij, i = j, w Ĥ ij = j n ds j, G ij = w ds j. j By applying these boundary conditions to every midpoint of elements, we obtain the equations in the matrix-vector form HU = GQ. Rearranging this system, a linear system of equations AX = B, where X is unknow values will be resulted. Example 3.. Consider the following Helmholtz equation with k = and the domain 0 < x <, 0 < y <, the boundary conditions in [4] 2π 4
13 BOUNDARY ELEMENT METHOD FOR u + k 2 u = 0, 0 < x <, 0 < y <, u(0, y) = 0, 0 < y <, 2 u(, y) = 2 cos(πy ), 0 < y <, 4 n = 0, 0 < x <, y=0 n = π 2 y= 8 sin(πx ), 0 < x <. 4 The exact solution of this particular boundary value problem is u(x, y) = sin( πx 4 ) cos(πy 4 ). In Table, a comparison between the exact solutions and the approximate solutions for differente values of elements of equation (6) of this element is given. Also in figure 3 and 4, these results are depicted. Table. The effect of increasing the number of elements on approximate solution helmholtz equation for interior points Figure 3. The effect of increasing the number of elements on approximate solution helmholtz equation
14
15 BOUNDARY ELEMENT METHOD FOR [9] Y. Khalighi and D. J. Bodony, Improved near-wall accuracy for solutions of the Helmholtz equation using the boundary element method, Center for Turbulence Research Annual Research Briefs, (2006), [0] K. Mohan Kadalbajoo and K. Arumugam, Boundary element method for elliptic differential equations with a small parameter, Applied mathematics and computation, (995), [] C. Walton Gibson, The Method of Moments in Electromagnetics, Chapman and Hall/CRC, [2] L. C. Wrobel, The Boundary Element Method: Applications in Thermo- Fluids and Acoustics, UK, [3] K. Yukio, S. Yonghao, and M. Zaheed, Regular boundary integral formulation for the analysis of open dielectric/optical waveguides, Memoirs of the Faculty of Engineering, 30 (995), [4] A. Ywhye-Teong, A beginner s course in boundary element methods, Universal Publishers, (2007), Ali Mesforush Department of Mathematics Assistant Professor of Mathematics School of Mathematics Shahrood University of Technology Shahrood, Iran ali.mesforush@shahroodut.ac.ir Zohreh Mehraban Department of Mathematics M.Sc Student of Mathematics School of Mathematics Shahrood University of Technology Shahrood, Iran zohre.mehraban@ymail.com
CHAPTER 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 informationImproved near-wall accuracy for solutions of the Helmholtz equation using the boundary element method
Center for Turbulence Research Annual Research Briefs 2006 313 Improved near-wall accuracy for solutions of the Helmholtz equation using the boundary element method By Y. Khalighi AND D. J. Bodony 1. Motivation
More informationu xx + u yy = 0. (5.1)
Chapter 5 Laplace Equation The following equation is called Laplace equation in two independent variables x, y: The non-homogeneous problem u xx + u yy =. (5.1) u xx + u yy = F, (5.) where F is a function
More informationSpotlight on Laplace s Equation
16 Spotlight on Laplace s Equation Reference: Sections 1.1,1.2, and 1.5. Laplace s equation is the undriven, linear, second-order PDE 2 u = (1) We defined diffusivity on page 587. where 2 is the Laplacian
More informationSOLUTIONS TO THE FINAL EXAM. December 14, 2010, 9:00am-12:00 (3 hours)
SOLUTIONS TO THE 18.02 FINAL EXAM BJORN POONEN December 14, 2010, 9:00am-12:00 (3 hours) 1) For each of (a)-(e) below: If the statement is true, write TRUE. If the statement is false, write FALSE. (Please
More informationCLASSIFICATION AND PRINCIPLE OF SUPERPOSITION FOR SECOND ORDER LINEAR PDE
CLASSIFICATION AND PRINCIPLE OF SUPERPOSITION FOR SECOND ORDER LINEAR PDE 1. Linear Partial Differential Equations A partial differential equation (PDE) is an equation, for an unknown function u, that
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 informationMath 234 Final Exam (with answers) Spring 2017
Math 234 Final Exam (with answers) pring 217 1. onsider the points A = (1, 2, 3), B = (1, 2, 2), and = (2, 1, 4). (a) [6 points] Find the area of the triangle formed by A, B, and. olution: One way to solve
More informationFundamental Solutions and Green s functions. Simulation Methods in Acoustics
Fundamental Solutions and Green s functions Simulation Methods in Acoustics Definitions Fundamental solution The solution F (x, x 0 ) of the linear PDE L {F (x, x 0 )} = δ(x x 0 ) x R d Is called the fundamental
More informationMath 106 Fall 2014 Exam 2.1 October 31, ln(x) x 3 dx = 1. 2 x 2 ln(x) + = 1 2 x 2 ln(x) + 1. = 1 2 x 2 ln(x) 1 4 x 2 + C
Math 6 Fall 4 Exam. October 3, 4. The following questions have to do with the integral (a) Evaluate dx. Use integration by parts (x 3 dx = ) ( dx = ) x3 x dx = x x () dx = x + x x dx = x + x 3 dx dx =
More informationLecture notes: Introduction to Partial Differential Equations
Lecture notes: Introduction to Partial Differential Equations Sergei V. Shabanov Department of Mathematics, University of Florida, Gainesville, FL 32611 USA CHAPTER 1 Classification of Partial Differential
More informationMath 113 Winter 2005 Key
Name Student Number Section Number Instructor Math Winter 005 Key Departmental Final Exam Instructions: The time limit is hours. Problem consists of short answer questions. Problems through are multiple
More informationMath 4263 Homework Set 1
Homework Set 1 1. Solve the following PDE/BVP 2. Solve the following PDE/BVP 2u t + 3u x = 0 u (x, 0) = sin (x) u x + e x u y = 0 u (0, y) = y 2 3. (a) Find the curves γ : t (x (t), y (t)) such that that
More informationIn this chapter we study elliptical PDEs. That is, PDEs of the form. 2 u = lots,
Chapter 8 Elliptic PDEs In this chapter we study elliptical PDEs. That is, PDEs of the form 2 u = lots, where lots means lower-order terms (u x, u y,..., u, f). Here are some ways to think about the physical
More informationAPPM 2350 FINAL EXAM FALL 2017
APPM 25 FINAL EXAM FALL 27. ( points) Determine the absolute maximum and minimum values of the function f(x, y) = 2 6x 4y + 4x 2 + y. Be sure to clearly give both the locations and values of the absolute
More informationPDEs, part 1: Introduction and elliptic PDEs
PDEs, part 1: Introduction and elliptic PDEs Anna-Karin Tornberg Mathematical Models, Analysis and Simulation Fall semester, 2013 Partial di erential equations The solution depends on several variables,
More information3 2 6 Solve the initial value problem u ( t) 3. a- If A has eigenvalues λ =, λ = 1 and corresponding eigenvectors 1
Math Problem a- If A has eigenvalues λ =, λ = 1 and corresponding eigenvectors 1 3 6 Solve the initial value problem u ( t) = Au( t) with u (0) =. 3 1 u 1 =, u 1 3 = b- True or false and why 1. if A is
More informationIntroduction to Differentials
Introduction to Differentials David G Radcliffe 13 March 2007 1 Increments Let y be a function of x, say y = f(x). The symbol x denotes a change or increment in the value of x. Note that a change in the
More informationINTRODUCTION TO PDEs
INTRODUCTION TO PDEs In this course we are interested in the numerical approximation of PDEs using finite difference methods (FDM). We will use some simple prototype boundary value problems (BVP) and initial
More informationChapter 3 Second Order Linear Equations
Partial Differential Equations (Math 3303) A Ë@ Õæ Aë áöß @. X. @ 2015-2014 ú GA JË@ É Ë@ Chapter 3 Second Order Linear Equations Second-order partial differential equations for an known function u(x,
More informationMultiple Integrals and Vector Calculus (Oxford Physics) Synopsis and Problem Sets; Hilary 2015
Multiple Integrals and Vector Calculus (Oxford Physics) Ramin Golestanian Synopsis and Problem Sets; Hilary 215 The outline of the material, which will be covered in 14 lectures, is as follows: 1. Introduction
More informationVariational Methods & Optimal Control
Variational Methods & Optimal Control lecture 12 Matthew Roughan Discipline of Applied Mathematics School of Mathematical Sciences University of Adelaide April 14, 216
More informationMath 350 Solutions for Final Exam Page 1. Problem 1. (10 points) (a) Compute the line integral. F ds C. z dx + y dy + x dz C
Math 35 Solutions for Final Exam Page Problem. ( points) (a) ompute the line integral F ds for the path c(t) = (t 2, t 3, t) with t and the vector field F (x, y, z) = xi + zj + xk. (b) ompute the line
More informationLaplace s Equation. Chapter Mean Value Formulas
Chapter 1 Laplace s Equation Let be an open set in R n. A function u C 2 () is called harmonic in if it satisfies Laplace s equation n (1.1) u := D ii u = 0 in. i=1 A function u C 2 () is called subharmonic
More informationMATH 425, FINAL EXAM SOLUTIONS
MATH 425, FINAL EXAM SOLUTIONS Each exercise is worth 50 points. Exercise. a The operator L is defined on smooth functions of (x, y by: Is the operator L linear? Prove your answer. L (u := arctan(xy u
More informationDifferential equations, comprehensive exam topics and sample questions
Differential equations, comprehensive exam topics and sample questions Topics covered ODE s: Chapters -5, 7, from Elementary Differential Equations by Edwards and Penney, 6th edition.. Exact solutions
More informationApplied Numerical Mathematics. High-order numerical schemes based on difference potentials for 2D elliptic problems with material interfaces
Applied Numerical Mathematics 111 (2017) 64 91 Contents lists available at ScienceDirect Applied Numerical Mathematics www.elsevier.com/locate/apnum High-order numerical schemes based on difference potentials
More information3 Green s functions in 2 and 3D
William J. Parnell: MT34032. Section 3: Green s functions in 2 and 3 57 3 Green s functions in 2 and 3 Unlike the one dimensional case where Green s functions can be found explicitly for a number of different
More information[2] (a) Develop and describe the piecewise linear Galerkin finite element approximation of,
269 C, Vese Practice problems [1] Write the differential equation u + u = f(x, y), (x, y) Ω u = 1 (x, y) Ω 1 n + u = x (x, y) Ω 2, Ω = {(x, y) x 2 + y 2 < 1}, Ω 1 = {(x, y) x 2 + y 2 = 1, x 0}, Ω 2 = {(x,
More information16.5 Surface Integrals of Vector Fields
16.5 Surface Integrals of Vector Fields Lukas Geyer Montana State University M73, Fall 011 Lukas Geyer (MSU) 16.5 Surface Integrals of Vector Fields M73, Fall 011 1 / 19 Parametrized Surfaces Definition
More informationDO NOT BEGIN THIS TEST UNTIL INSTRUCTED TO START
Math 265 Student name: KEY Final Exam Fall 23 Instructor & Section: This test is closed book and closed notes. A (graphing) calculator is allowed for this test but cannot also be a communication device
More informationPractice Exam 1 Solutions
Practice Exam 1 Solutions 1a. Let S be the region bounded by y = x 3, y = 1, and x. Find the area of S. What is the volume of the solid obtained by rotating S about the line y = 1? Area A = Volume 1 1
More informationSimple Examples on Rectangular Domains
84 Chapter 5 Simple Examples on Rectangular Domains In this chapter we consider simple elliptic boundary value problems in rectangular domains in R 2 or R 3 ; our prototype example is the Poisson equation
More informationSample Final Questions: Solutions Math 21B, Winter y ( y 1)(1 + y)) = A y + B
Sample Final Questions: Solutions Math 2B, Winter 23. Evaluate the following integrals: tan a) y y dy; b) x dx; c) 3 x 2 + x dx. a) We use partial fractions: y y 3 = y y ) + y)) = A y + B y + C y +. Putting
More informationTheory of PDE Homework 2
Theory of PDE Homework 2 Adrienne Sands April 18, 2017 In the following exercises we assume the coefficients of the various PDE are smooth and satisfy the uniform ellipticity condition. R n is always an
More informationPractice Problems for Exam 3 (Solutions) 1. Let F(x, y) = xyi+(y 3x)j, and let C be the curve r(t) = ti+(3t t 2 )j for 0 t 2. Compute F dr.
1. Let F(x, y) xyi+(y 3x)j, and let be the curve r(t) ti+(3t t 2 )j for t 2. ompute F dr. Solution. F dr b a 2 2 F(r(t)) r (t) dt t(3t t 2 ), 3t t 2 3t 1, 3 2t dt t 3 dt 1 2 4 t4 4. 2. Evaluate the line
More informationSeries Solutions of ODEs. Special Functions
C05.tex 6/4/0 3: 5 Page 65 Chap. 5 Series Solutions of ODEs. Special Functions We continue our studies of ODEs with Legendre s, Bessel s, and the hypergeometric equations. These ODEs have variable coefficients
More informationJim Lambers MAT 280 Summer Semester Practice Final Exam Solution. dy + xz dz = x(t)y(t) dt. t 3 (4t 3 ) + e t2 (2t) + t 7 (3t 2 ) dt
Jim Lambers MAT 28 ummer emester 212-1 Practice Final Exam olution 1. Evaluate the line integral xy dx + e y dy + xz dz, where is given by r(t) t 4, t 2, t, t 1. olution From r (t) 4t, 2t, t 2, we obtain
More informationStokes Theorem. MATH 311, Calculus III. J. Robert Buchanan. Summer Department of Mathematics. J. Robert Buchanan Stokes Theorem
tokes Theorem MATH 311, alculus III J. Robert Buchanan Department of Mathematics ummer 2011 Background (1 of 2) Recall: Green s Theorem, M(x, y) dx + N(x, y) dy = R ( N x M ) da y where is a piecewise
More informationCalculus II - Fall 2013
Calculus II - Fall Midterm Exam II, November, In the following problems you are required to show all your work and provide the necessary explanations everywhere to get full credit.. Find the area between
More informationLecture Notes on PDEs
Lecture Notes on PDEs Alberto Bressan February 26, 2012 1 Elliptic equations Let IR n be a bounded open set Given measurable functions a ij, b i, c : IR, consider the linear, second order differential
More informationThe Helmholtz theorem at last!
Problem. The Helmholtz theorem at last! Recall in class the Helmholtz theorem that says that if if E =0 then E can be written as E = φ () if B =0 then B can be written as B = A (2) (a) Let n be a unit
More informationMath 342 Partial Differential Equations «Viktor Grigoryan
Math 342 Partial ifferential Equations «Viktor Grigoryan 3 Green s first identity Having studied Laplace s equation in regions with simple geometry, we now start developing some tools, which will lead
More informationAPPENDIX 2.1 LINE AND SURFACE INTEGRALS
2 APPENDIX 2. LINE AND URFACE INTEGRAL Consider a path connecting points (a) and (b) as shown in Fig. A.2.. Assume that a vector field A(r) exists in the space in which the path is situated. Then the line
More informationUNIVERSITY of LIMERICK OLLSCOIL LUIMNIGH
UNIVERSITY of LIMERICK OLLSCOIL LUIMNIGH Faculty of Science and Engineering Department of Mathematics and Statistics END OF SEMESTER ASSESSMENT PAPER MODULE CODE: MA4006 SEMESTER: Spring 2011 MODULE TITLE:
More informationGREEN S IDENTITIES AND GREEN S FUNCTIONS
GREEN S IENTITIES AN GREEN S FUNCTIONS Green s first identity First, recall the following theorem. Theorem: (ivergence Theorem) Let be a bounded solid region with a piecewise C 1 boundary surface. Let
More informationGreen s Functions and Distributions
CHAPTER 9 Green s Functions and Distributions 9.1. Boundary Value Problems We would like to study, and solve if possible, boundary value problems such as the following: (1.1) u = f in U u = g on U, where
More informationMath Exam III - Spring
Math 3 - Exam III - Spring 8 This exam contains 5 multiple choice questions and hand graded questions. The multiple choice questions are worth 5 points each and the hand graded questions are worth a total
More informationFirst order Partial Differential equations
First order Partial Differential equations 0.1 Introduction Definition 0.1.1 A Partial Deferential equation is called linear if the dependent variable and all its derivatives have degree one and not multiple
More informationENGI Multiple Integration Page 8-01
ENGI 345 8. Multiple Integration Page 8-01 8. Multiple Integration This chapter provides only a very brief introduction to the major topic of multiple integration. Uses of multiple integration include
More informationMethod of Green s Functions
Method of Green s Functions 18.33 Linear Partial ifferential Equations Matthew J. Hancock Fall 24 Ref: Haberman, h 9, 11 We introduce another powerful method of solving PEs. First, we need to consider
More informationPrelim 2 Math Please show your reasoning and all your work. This is a 90 minute exam. Calculators are not needed or permitted. Good luck!
April 4, Prelim Math Please show your reasoning and all your work. This is a 9 minute exam. Calculators are not needed or permitted. Good luck! Trigonometric Formulas sin x sin x cos x cos (u + v) cos
More informationQualifying Examination
Summer 24 Day. Monday, September 5, 24 You have three hours to complete this exam. Work all problems. Start each problem on a All problems are 2 points. Please email any electronic files in support of
More informationPRACTICE PROBLEMS. Please let me know if you find any mistakes in the text so that i can fix them. 1. Mixed partial derivatives.
PRACTICE PROBLEMS Please let me know if you find any mistakes in the text so that i can fix them. 1.1. Let Show that f is C 1 and yet How is that possible? 1. Mixed partial derivatives f(x, y) = {xy x
More informationMath 31CH - Spring Final Exam
Math 3H - Spring 24 - Final Exam Problem. The parabolic cylinder y = x 2 (aligned along the z-axis) is cut by the planes y =, z = and z = y. Find the volume of the solid thus obtained. Solution:We calculate
More informationFrom the Circle to the Square: Symmetry and Degeneracy in Quantum Mechanics. Dahyeon Lee Department of Physics and Astronomy, Oberlin College
From the Circle to the Square: Symmetry and Degeneracy in Quantum Mechanics Dahyeon Lee Department of Physics and Astronomy, Oberlin College March 31, 2017 Executive Summary The relationship between degeneracy
More informationMath Homework 2
Math 73 Homework Due: September 8, 6 Suppose that f is holomorphic in a region Ω, ie an open connected set Prove that in any of the following cases (a) R(f) is constant; (b) I(f) is constant; (c) f is
More informationElectromagnetism HW 1 math review
Electromagnetism HW math review Problems -5 due Mon 7th Sep, 6- due Mon 4th Sep Exercise. The Levi-Civita symbol, ɛ ijk, also known as the completely antisymmetric rank-3 tensor, has the following properties:
More informationUNIVERSITY of LIMERICK OLLSCOIL LUIMNIGH
UNIVERSITY of LIMERICK OLLSCOIL LUIMNIGH College of Informatics and Electronics END OF SEMESTER ASSESSMENT PAPER MODULE CODE: MS425 SEMESTER: Autumn 25/6 MODULE TITLE: Applied Analysis DURATION OF EXAMINATION:
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 informationMATH 228: Calculus III (FALL 2016) Sample Problems for FINAL EXAM SOLUTIONS
MATH 228: Calculus III (FALL 216) Sample Problems for FINAL EXAM SOLUTIONS MATH 228 Page 2 Problem 1. (2pts) Evaluate the line integral C xy dx + (x + y) dy along the parabola y x2 from ( 1, 1) to (2,
More informationPartial Differential Equations
M3M3 Partial Differential Equations Solutions to problem sheet 3/4 1* (i) Show that the second order linear differential operators L and M, defined in some domain Ω R n, and given by Mφ = Lφ = j=1 j=1
More informationMath RE - Calculus II Antiderivatives and the Indefinite Integral Page 1 of 5
Math 201-203-RE - Calculus II Antiderivatives and the Indefinite Integral Page 1 of 5 What is the Antiderivative? In a derivative problem, a function f(x) is given and you find the derivative f (x) using
More informationPartial Differential Equations
Partial Differential Equations Xu Chen Assistant Professor United Technologies Engineering Build, Rm. 382 Department of Mechanical Engineering University of Connecticut xchen@engr.uconn.edu Contents 1
More informationFysikens matematiska metoder week: Semester: Spring 2015 (1FA121) Lokal: Bergsbrunnagatan 15, Sal 2 Tid: 08:00-13:00 General remark: Material al
Fysikens matematiska metoder week: 13-23 Semester: Spring 2015 (1FA121) Lokal: Bergsbrunnagatan 15, Sal 2 Tid: 08:00-13:00 General remark: Material allowed to be taken with you to the exam: Physics Handbook,
More informationCircular Membranes. Farlow, Lesson 30. November 21, Department of Mathematical Sciences Florida Atlantic University Boca Raton, FL 33431
The Problem Polar coordinates Solving the Problem by Separation of Variables Circular Membranes Farlow, Lesson 30 Department of Mathematical Sciences Florida Atlantic University Boca Raton, FL 33431 November
More informationMathematics (Course B) Lent Term 2005 Examples Sheet 2
N12d Natural Sciences, Part IA Dr M. G. Worster Mathematics (Course B) Lent Term 2005 Examples Sheet 2 Please communicate any errors in this sheet to Dr Worster at M.G.Worster@damtp.cam.ac.uk. Note that
More informationThe double layer potential
The double layer potential In this project, our goal is to explain how the Dirichlet problem for a linear elliptic partial differential equation can be converted into an integral equation by representing
More informationx 2 y = 1 2. Problem 2. Compute the Taylor series (at the base point 0) for the function 1 (1 x) 3.
MATH 8.0 - FINAL EXAM - SOME REVIEW PROBLEMS WITH SOLUTIONS 8.0 Calculus, Fall 207 Professor: Jared Speck Problem. Consider the following curve in the plane: x 2 y = 2. Let a be a number. The portion of
More informationPower Series. Part 1. J. Gonzalez-Zugasti, University of Massachusetts - Lowell
Power Series Part 1 1 Power Series Suppose x is a variable and c k & a are constants. A power series about x = 0 is c k x k A power series about x = a is c k x a k a = center of the power series c k =
More informationThe purpose of this lecture is to present a few applications of conformal mappings in problems which arise in physics and engineering.
Lecture 16 Applications of Conformal Mapping MATH-GA 451.001 Complex Variables The purpose of this lecture is to present a few applications of conformal mappings in problems which arise in physics and
More informationWave Field Synthesis
Wave Field ynthesis Ramani Duraiswami April 26, 2006 Abstract ome background material for the technique called wave-field synthesis is provided. THen we move on to discuss some papers from the literature.
More informationProblem Set 1. This week. Please read all of Chapter 1 in the Strauss text.
Math 425, Spring 2015 Jerry L. Kazdan Problem Set 1 Due: Thurs. Jan. 22 in class. [Late papers will be accepted until 1:00 PM Friday.] This is rust remover. It is essentially Homework Set 0 with a few
More informationElementary Differential Equations, Section 2 Prof. Loftin: Practice Test Problems for Test Find the radius of convergence of the power series
Elementary Differential Equations, Section 2 Prof. Loftin: Practice Test Problems for Test 2 SOLUTIONS 1. Find the radius of convergence of the power series Show your work. x + x2 2 + x3 3 + x4 4 + + xn
More informationS12.1 SOLUTIONS TO PROBLEMS 12 (ODD NUMBERS)
OLUTION TO PROBLEM 2 (ODD NUMBER) 2. The electric field is E = φ = 2xi + 2y j and at (2, ) E = 4i + 2j. Thus E = 2 5 and its direction is 2i + j. At ( 3, 2), φ = 6i + 4 j. Thus the direction of most rapid
More information1 Review of di erential calculus
Review of di erential calculus This chapter presents the main elements of di erential calculus needed in probability theory. Often, students taking a course on probability theory have problems with concepts
More information2.20 Fall 2018 Math Review
2.20 Fall 2018 Math Review September 10, 2018 These notes are to help you through the math used in this class. This is just a refresher, so if you never learned one of these topics you should look more
More informationWe denote the space of distributions on Ω by D ( Ω) 2.
Sep. 1 0, 008 Distributions Distributions are generalized functions. Some familiarity with the theory of distributions helps understanding of various function spaces which play important roles in the study
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Numerical Methods for Partial Differential Equations Finite Difference Methods
More informationLECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES)
LECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES) RAYTCHO LAZAROV 1 Notations and Basic Functional Spaces Scalar function in R d, d 1 will be denoted by u,
More informationFunctions of Several Variables
Functions of Several Variables Partial Derivatives Philippe B Laval KSU March 21, 2012 Philippe B Laval (KSU) Functions of Several Variables March 21, 2012 1 / 19 Introduction In this section we extend
More informationMath 263 Final. (b) The cross product is. i j k c. =< c 1, 1, 1 >
Math 63 Final Problem 1: [ points, 5 points to each part] Given the points P : (1, 1, 1), Q : (1,, ), R : (,, c 1), where c is a parameter, find (a) the vector equation of the line through P and Q. (b)
More information(a) x cos 3x dx We apply integration by parts. Take u = x, so that dv = cos 3x dx, v = 1 sin 3x, du = dx. Thus
Math 128 Midterm Examination 2 October 21, 28 Name 6 problems, 112 (oops) points. Instructions: Show all work partial credit will be given, and Answers without work are worth credit without points. You
More informationy = x 3 and y = 2x 2 x. 2x 2 x = x 3 x 3 2x 2 + x = 0 x(x 2 2x + 1) = 0 x(x 1) 2 = 0 x = 0 and x = (x 3 (2x 2 x)) dx
Millersville University Name Answer Key Mathematics Department MATH 2, Calculus II, Final Examination May 4, 2, 8:AM-:AM Please answer the following questions. Your answers will be evaluated on their correctness,
More informationMATH COURSE NOTES - CLASS MEETING # Introduction to PDEs, Spring 2018 Professor: Jared Speck
MATH 8.52 COURSE NOTES - CLASS MEETING # 6 8.52 Introduction to PDEs, Spring 208 Professor: Jared Speck Class Meeting # 6: Laplace s and Poisson s Equations We will now study the Laplace and Poisson equations
More informationG G. G. x = u cos v, y = f(u), z = u sin v. H. x = u + v, y = v, z = u v. 1 + g 2 x + g 2 y du dv
1. Matching. Fill in the appropriate letter. 1. ds for a surface z = g(x, y) A. r u r v du dv 2. ds for a surface r(u, v) B. r u r v du dv 3. ds for any surface C. G x G z, G y G z, 1 4. Unit normal N
More informationlim = F F = F x x + F y y + F z
Physics 361 Summary of Results from Lecture Physics 361 Derivatives of Scalar and Vector Fields The gradient of a scalar field f( r) is given by g = f. coordinates f g = ê x x + ê f y y + ê f z z Expressed
More informationMATH COURSE NOTES - CLASS MEETING # Introduction to PDEs, Fall 2011 Professor: Jared Speck
MATH 8.52 COURSE NOTES - CLASS MEETING # 6 8.52 Introduction to PDEs, Fall 20 Professor: Jared Speck Class Meeting # 6: Laplace s and Poisson s Equations We will now study the Laplace and Poisson equations
More informationMath 233. Practice Problems Chapter 15. i j k
Math 233. Practice Problems hapter 15 1. ompute the curl and divergence of the vector field F given by F (4 cos(x 2 ) 2y)i + (4 sin(y 2 ) + 6x)j + (6x 2 y 6x + 4e 3z )k olution: The curl of F is computed
More informationFinal Exam. Monday March 19, 3:30-5:30pm MAT 21D, Temple, Winter 2018
Name: Student ID#: Section: Final Exam Monday March 19, 3:30-5:30pm MAT 21D, Temple, Winter 2018 Show your work on every problem. orrect answers with no supporting work will not receive full credit. Be
More informationFinal Exam 2011 Winter Term 2 Solutions
. (a Find the radius of convergence of the series: ( k k+ x k. Solution: Using the Ratio Test, we get: L = lim a k+ a k = lim ( k+ k+ x k+ ( k k+ x k = lim x = x. Note that the series converges for L
More informationMcGill University April 16, Advanced Calculus for Engineers
McGill University April 16, 2014 Faculty of cience Final examination Advanced Calculus for Engineers Math 264 April 16, 2014 Time: 6PM-9PM Examiner: Prof. R. Choksi Associate Examiner: Prof. A. Hundemer
More informationMath Practice Exam 2 - solutions
C Roettger, Fall 205 Math 66 - Practice Exam 2 - solutions State clearly what your result is. Show your work (in particular, integrand and limits of integrals, all substitutions, names of tests used, with
More informationNote: Final Exam is at 10:45 on Tuesday, 5/3/11 (This is the Final Exam time reserved for our labs). From Practice Test I
MA Practice Final Answers in Red 4/8/ and 4/9/ Name Note: Final Exam is at :45 on Tuesday, 5// (This is the Final Exam time reserved for our labs). From Practice Test I Consider the integral 5 x dx. Sketch
More informationx+ y = 50 Dividing both sides by 2 : ( ) dx By (7.2), x = 25m gives maximum area. Substituting this value into (*):
Solutions 7(b 1 Complete solutions to Exercise 7(b 1. Since the perimeter 100 we have x+ y 100 [ ] ( x+ y 50 ividing both sides by : y 50 x * The area A xy, substituting y 50 x gives: A x( 50 x A 50x x
More informationEXPOSITORY NOTES ON DISTRIBUTION THEORY, FALL 2018
EXPOSITORY NOTES ON DISTRIBUTION THEORY, FALL 2018 While these notes are under construction, I expect there will be many typos. The main reference for this is volume 1 of Hörmander, The analysis of liner
More informationMath 11 Fall 2018 Practice Final Exam
Math 11 Fall 218 Practice Final Exam Disclaimer: This practice exam should give you an idea of the sort of questions we may ask on the actual exam. Since the practice exam (like the real exam) is not long
More informationMATH 131P: PRACTICE FINAL SOLUTIONS DECEMBER 12, 2012
MATH 3P: PRACTICE FINAL SOLUTIONS DECEMBER, This is a closed ook, closed notes, no calculators/computers exam. There are 6 prolems. Write your solutions to Prolems -3 in lue ook #, and your solutions to
More informationFinal exam (practice) UCLA: Math 31B, Spring 2017
Instructor: Noah White Date: Final exam (practice) UCLA: Math 3B, Spring 207 This exam has 8 questions, for a total of 80 points. Please print your working and answers neatly. Write your solutions in the
More informationQ ( q(m, t 0 ) n) S t.
THE HEAT EQUATION The main equations that we will be dealing with are the heat equation, the wave equation, and the potential equation. We use simple physical principles to show how these equations are
More information