Using Green s functions with inhomogeneous BCs
|
|
- Rodney McDowell
- 5 years ago
- Views:
Transcription
1 Using Green s functions with inhomogeneous BCs
2 Using Green s functions with inhomogeneous BCs Surprise: Although Green s functions satisfy homogeneous boundary conditions, they can be used for problems with inhomogeneous BCs!
3 Using Green s functions with inhomogeneous BCs Surprise: Although Green s functions satisfy homogeneous boundary conditions, they can be used for problems with inhomogeneous BCs! For self adjoint L and u, v with homogeneous boundary conditions it follows that (Lv)u dx (Lu)v dx = 0.
4 Using Green s functions with inhomogeneous BCs Surprise: Although Green s functions satisfy homogeneous boundary conditions, they can be used for problems with inhomogeneous BCs! For self adjoint L and u, v with homogeneous boundary conditions it follows that (Lv)u dx (Lu)v dx = 0. But if u, v don t satisfy homogeneous boundary conditions, get (Lv)u dx (Lu)v dx = boundary terms involving u and v. This is called the Green s formula, which depends on L and.
5 Using Green s formula for inhomogeneous boundary conditions Want to solve Lu(x) = f (x) + inhomogeneous boundary conditions.
6 Using Green s formula for inhomogeneous boundary conditions Want to solve Lu(x) = f (x) + inhomogeneous boundary conditions. Set v(x) = G(x, x 0 ) in Green s formula (LG)u dx (Lu)G dx = boundary terms.
7 Using Green s formula for inhomogeneous boundary conditions Want to solve Lu(x) = f (x) + inhomogeneous boundary conditions. Set v(x) = G(x, x 0 ) in Green s formula (LG)u dx (Lu)G dx = boundary terms. Since LG = δ(x x 0 ) and Lu = f, we have δ(x x 0 )u(x)dx f (x)g(x, x 0 )dx = boundary terms
8 Using Green s formula for inhomogeneous boundary conditions Want to solve Lu(x) = f (x) + inhomogeneous boundary conditions. Set v(x) = G(x, x 0 ) in Green s formula (LG)u dx (Lu)G dx = boundary terms. Since LG = δ(x x 0 ) and Lu = f, we have δ(x x 0 )u(x)dx f (x)g(x, x 0 )dx = boundary terms Collapsing the integral involving the δ function, u(x 0 ) = G(x, x 0 )f (x)dx + boundary terms
9 Green s formula for Laplacian Want to solve u = f in, u = h on,
10 Green s formula for Laplacian Want to solve u = f in, u = h on, For dimensions 2, the Green s formula is just Green s identity u v v u dx = u v ˆn v u ˆn dx.
11 Green s formula for Laplacian Want to solve u = f in, u = h on, For dimensions 2, the Green s formula is just Green s identity u v v u dx = u v ˆn v u ˆn dx. Let G solve G = δ(x x 0 ) and G = 0 on boundary.
12 Green s formula for Laplacian Want to solve u = f in, u = h on, For dimensions 2, the Green s formula is just Green s identity u v v u dx = u v ˆn v u ˆn dx. Let G solve G = δ(x x 0 ) and G = 0 on boundary. Substituting v(x) = G(x, x 0 ) into Green s formula, u(x)δ(x x 0 ) G(x, x 0 )f (x) dx = u(x) xg(x, x 0 ) ˆn(x) G(x, x 0 ) u(x) ˆn(x) dx
13 Green s formula for Laplacian Want to solve u = f in, u = h on, For dimensions 2, the Green s formula is just Green s identity u v v u dx = u v ˆn v u ˆn dx. Let G solve G = δ(x x 0 ) and G = 0 on boundary. Substituting v(x) = G(x, x 0 ) into Green s formula, u(x)δ(x x 0 ) G(x, x 0 )f (x) dx = u(x) xg(x, x 0 ) ˆn(x) G(x, x 0 ) u(x) ˆn(x) dx Simplifies to u(x 0 ) = G(x, x 0 )f (x)dx + h(x) x G(x, x 0 ) ˆn(x) dx,
14 Example: the Poisson integral formula revisited In the case that is a disk of radius a, Green s function is G(r, θ; r 0, θ 0 ) = 1 ( a 2 4π ln r 2 0 r 2 + r 2 0 2rr 0 cos(θ θ 0 ) r 2 + a 4 /r 2 0 2ra2 /r 0 cos(θ θ 0 ) ).
15 Example: the Poisson integral formula revisited In the case that is a disk of radius a, Green s function is G(r, θ; r 0, θ 0 ) = 1 ( a 2 4π ln The boundary value problem has a solution u(x 0 ) = r 2 0 u = 0, r 2 + r 2 0 2rr 0 cos(θ θ 0 ) r 2 + a 4 /r 2 0 2ra2 /r 0 cos(θ θ 0 ) u(a, θ) = h(θ) h(x) xg(x, x 0 ) ˆn(x) dx. ).
16 Example: the Poisson integral formula revisited In the case that is a disk of radius a, Green s function is G(r, θ; r 0, θ 0 ) = 1 ( a 2 4π ln The boundary value problem has a solution u(x 0 ) = Need normal derivative of G r 2 0 u = 0, xg(x, x 0 ) ˆn(x) = G r (r, θ; r 0, θ 0 ) = 1 4π which at r = a is r 2 + r 2 0 2rr 0 cos(θ θ 0 ) r 2 + a 4 /r 2 0 2ra2 /r 0 cos(θ θ 0 ) u(a, θ) = h(θ) h(x) xg(x, x 0 ) ˆn(x) dx. ). ( 2r0 2r cos(θ θ 0 ) r 2 + r 2 0 2rr 0 cos(θ θ 0 ) 2r 0 r 2 2ra 2 cos(θ θ 0 ) r 2 r a4 2rr 0 a 2 cos(θ θ 0 ) ( a 2π 1 (r/a) 2 r 2 + a 2 2ar cos(θ θ 0 ) ). ),
17 Example: the Poisson integral formula revisited In the case that is a disk of radius a, Green s function is G(r, θ; r 0, θ 0 ) = 1 ( a 2 4π ln The boundary value problem has a solution u(x 0 ) = Need normal derivative of G r 2 0 u = 0, xg(x, x 0 ) ˆn(x) = G r (r, θ; r 0, θ 0 ) = 1 4π r 2 + r 2 0 2rr 0 cos(θ θ 0 ) r 2 + a 4 /r 2 0 2ra2 /r 0 cos(θ θ 0 ) u(a, θ) = h(θ) h(x) xg(x, x 0 ) ˆn(x) dx. ). ( 2r0 2r cos(θ θ 0 ) r 2 + r 2 0 2rr 0 cos(θ θ 0 ) 2r 0 r 2 2ra 2 cos(θ θ 0 ) r 2 r a4 2rr 0 a 2 cos(θ θ 0 ) which at r = a is ( ) a 1 (r/a) 2. 2π r 2 + a 2 2ar cos(θ θ 0 ) Parameterize boundary integral using θ and dx = a dθ, u(r 0, θ 0 ) = 1 2π 2π 0 (a 2 r 2 0 )h(θ) a 2 + r 2 0 2ar 0 cos(θ θ 0 ) dθ. ),
18 Neumann boundary conditions Want to solve u = 0, lim u(x, y, z) = 0, z u z (x, y, 0) = h(x, y), in upper half space {(x, y, z) z > 0}.
19 Neumann boundary conditions Want to solve u = 0, lim u(x, y, z) = 0, z u z (x, y, 0) = h(x, y), in upper half space {(x, y, z) z > 0}. Green s formula u v v u dx = u v ˆn v u ˆn dx. has both Dirichlet and Neumann boundary terms in u, but only know u(x, y, 0) ˆn = u z (x, y, 0).
20 Neumann boundary conditions Want to solve u = 0, lim u(x, y, z) = 0, z u z (x, y, 0) = h(x, y), in upper half space {(x, y, z) z > 0}. Green s formula u v v u dx = u v ˆn v u ˆn dx. has both Dirichlet and Neumann boundary terms in u, but only know u(x, y, 0) ˆn = u z (x, y, 0). To make v ˆn = G ˆn vanish on boundary, need Green s function to respect boundary condition principle": The Green s function must have the same type of boundary conditions as the problem to be solved, and they must be homogeneous.
21 Neumann boundary conditions,cont. Method of images prescribes even reflection so G z = 0 when z = 0: G(x, y, z; x 0, y 0, z 0 ) = G 3 (x, y, z; x 0, y 0, z 0 ) + G 3 (x, y, z; x 0, y 0, z 0 ) ( ) = 1 1 4π (x x0 ) 2 + (y y 0 ) 2 + (z z 0 ) (x x0 ) 2 + (y y 0 ) 2 + (z + z 0 ) 2
22 Neumann boundary conditions,cont. Method of images prescribes even reflection so G z = 0 when z = 0: G(x, y, z; x 0, y 0, z 0 ) = G 3 (x, y, z; x 0, y 0, z 0 ) + G 3 (x, y, z; x 0, y 0, z 0 ) ( ) = 1 1 4π (x x0 ) 2 + (y y 0 ) 2 + (z z 0 ) (x x0 ) 2 + (y y 0 ) 2 + (z + z 0 ) 2 Substituting v(x) = G(x, x 0 ) into Green s formula and collapsing the δ-function integral, u(x 0 ) = G(x; x 0 ) u(x) ˆn(x) dx,
23 Neumann boundary conditions,cont. Method of images prescribes even reflection so G z = 0 when z = 0: G(x, y, z; x 0, y 0, z 0 ) = G 3 (x, y, z; x 0, y 0, z 0 ) + G 3 (x, y, z; x 0, y 0, z 0 ) ( ) = 1 1 4π (x x0 ) 2 + (y y 0 ) 2 + (z z 0 ) (x x0 ) 2 + (y y 0 ) 2 + (z + z 0 ) 2 Substituting v(x) = G(x, x 0 ) into Green s formula and collapsing the δ-function integral, u(x 0 ) = G(x; x 0 ) u(x) ˆn(x) dx, Boundary is both xy-plane and the effective boundary at infinity, but integrand vanishes on the latter.
24 Neumann boundary conditions,cont. Method of images prescribes even reflection so G z = 0 when z = 0: G(x, y, z; x 0, y 0, z 0 ) = G 3 (x, y, z; x 0, y 0, z 0 ) + G 3 (x, y, z; x 0, y 0, z 0 ) ( ) = 1 1 4π (x x0 ) 2 + (y y 0 ) 2 + (z z 0 ) (x x0 ) 2 + (y y 0 ) 2 + (z + z 0 ) 2 Substituting v(x) = G(x, x 0 ) into Green s formula and collapsing the δ-function integral, u(x 0 ) = G(x; x 0 ) u(x) ˆn(x) dx, Boundary is both xy-plane and the effective boundary at infinity, but integrand vanishes on the latter. Since ˆn is directed outward, u(x) ˆn(x) = u z (x, y, 0), and u(x 0, y 0, z 0 ) = 1 2π h(x, y) (x x 0 ) 2 + (y y 0 ) 2 + z0 2 dxdy.
25 Symmetry (reciprocity) of the Green s function If L is self-adjoint, might expect that its inverse to also be: G(x, x 0 ) = G(x 0, x), Reciprocity"
26 Symmetry (reciprocity) of the Green s function If L is self-adjoint, might expect that its inverse to also be: G(x, x 0 ) = G(x 0, x), Reciprocity" Proof: Insert v(x) = G(x, x 1 ), u(x) = G(x, x 2 ), Lv = δ(x x 1 ) and Lu = δ(x x 2 ) into Green s formula: δ(x x 1 )G(x, x 2 )dx G(x, x 1 )δ(x x 2 )dx = 0 which simplifies to G(x 1, x 2 ) G(x 2, x 1 ) = 0.
27 Symmetry (reciprocity) of the Green s function If L is self-adjoint, might expect that its inverse to also be: G(x, x 0 ) = G(x 0, x), Reciprocity" Proof: Insert v(x) = G(x, x 1 ), u(x) = G(x, x 2 ), Lv = δ(x x 1 ) and Lu = δ(x x 2 ) into Green s formula: δ(x x 1 )G(x, x 2 )dx G(x, x 1 )δ(x x 2 )dx = 0 which simplifies to G(x 1, x 2 ) G(x 2, x 1 ) = 0. Also: can interchange arguments of partial derivatives, e.g. xg(x, x 0 ) = lim h 0 (G(x + h, x 0 ) G(x, x 0 ))/h = lim h 0 (G(x 0, x + h) G(x 0, x))/h = x0 G(x 0, x).
28 Symmetry (reciprocity) of the Green s function If L is self-adjoint, might expect that its inverse to also be: G(x, x 0 ) = G(x 0, x), Reciprocity" Proof: Insert v(x) = G(x, x 1 ), u(x) = G(x, x 2 ), Lv = δ(x x 1 ) and Lu = δ(x x 2 ) into Green s formula: δ(x x 1 )G(x, x 2 )dx G(x, x 1 )δ(x x 2 )dx = 0 which simplifies to G(x 1, x 2 ) G(x 2, x 1 ) = 0. Also: can interchange arguments of partial derivatives, e.g. xg(x, x 0 ) = lim h 0 (G(x + h, x 0 ) G(x, x 0 ))/h = lim h 0 (G(x 0, x + h) G(x 0, x))/h = x0 G(x 0, x). For example, representation formula u(x 0 ) = G(x, x 0 )f (x)dx + h(x) xg(x, x 0 ) ˆn(x) dx, can be rewritten by exchanging the notation for x and x 0 and using reciprocity, u(x) = G(x, x 0 )f (x 0 )dx 0 + h(x 0 ) x0 G(x, x 0 ) ˆn(x 0 ) dx 0.
LAPLACE EQUATION. = 2 is call the Laplacian or del-square operator. In two dimensions, the expression of in rectangular and polar coordinates are
LAPLACE EQUATION If a diffusion or wave problem is stationary (time independent), the pde reduces to the Laplace equation u = u =, an archetype of second order elliptic pde. = 2 is call the Laplacian or
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 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 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 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 informationModule 7: The Laplace Equation
Module 7: The Laplace Equation In this module, we shall study one of the most important partial differential equations in physics known as the Laplace equation 2 u = 0 in Ω R n, (1) where 2 u := n i=1
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 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 information3.6 Applications of Poisson s equation
William J. Parnell: MT3432. Section 3: Green s functions in 2 and 3D 69 3.6 Applications of Poisson s equation The beaut of Green s functions is that the are useful! Let us consider an example here where
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 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 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 information18 Green s function for the Poisson equation
8 Green s function for the Poisson equation Now we have some experience working with Green s functions in dimension, therefore, we are ready to see how Green s functions can be obtained in dimensions 2
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 informationCalculus II Practice Test Problems for Chapter 7 Page 1 of 6
Calculus II Practice Test Problems for Chapter 7 Page of 6 This is a set of practice test problems for Chapter 7. This is in no way an inclusive set of problems there can be other types of problems on
More informationElements of linear algebra
Elements of linear algebra Elements of linear algebra A vector space S is a set (numbers, vectors, functions) which has addition and scalar multiplication defined, so that the linear combination c 1 v
More information13 PDEs on spatially bounded domains: initial boundary value problems (IBVPs)
13 PDEs on spatially bounded domains: initial boundary value problems (IBVPs) A prototypical problem we will discuss in detail is the 1D diffusion equation u t = Du xx < x < l, t > finite-length rod u(x,
More informationIntroduction and some preliminaries
1 Partial differential equations Introduction and some preliminaries A partial differential equation (PDE) is a relationship among partial derivatives of a function (or functions) of more than one variable.
More informationAssignment 2 - Complex Analysis
Assignment 2 - Complex Analysis MATH 440/508 M.P. Lamoureux Sketch of solutions. Γ z dz = Γ (x iy)(dx + idy) = (xdx + ydy) + i Γ Γ ( ydx + xdy) = (/2)(x 2 + y 2 ) endpoints + i [( / y) y ( / x)x]dxdy interiorγ
More informationThe integral test and estimates of sums
The integral test Suppose f is a continuous, positive, decreasing function on [, ) and let a n = f (n). Then the series n= a n is convergent if and only if the improper integral f (x)dx is convergent.
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 informationSuggested Solution to Assignment 6
MATH 4 (6-7) partial diferential equations Suggested Solution to Assignment 6 Exercise 6.. Note that in the spherical coordinates (r, θ, φ), Thus, Let u = v/r, we get 3 = r + r r + r sin θ θ sin θ θ +
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 informationJun-ichi Takada
Wave Theory II Numerical Simulation of Waves 2) Representation of Wave Function by Using Green Function I) Green s Theorem and Free Space Green Function Jun-ichi Takada takada@ide.titech.ac.jp)) This lecture
More informationSECTION 3.6 THE CONSTRUCTION OF THE GREEN'S FUNCTION FOR THE TWO DIMENSIONAL OPERATOR WITH DIRICHLET BOUNDARY CONDITIONS
SECTION 3.6 THE CONSTRUCTION OF THE GREEN'S FUNCTION FOR THE TWO DIMENSIONAL OPERATOR WITH DIRICHLET BOUNDARY CONDITIONS We now will construct a Green's function for the following problem: Find u such
More information21 Laplace s Equation and Harmonic Functions
2 Laplace s Equation and Harmonic Functions 2. Introductory Remarks on the Laplacian operator Given a domain Ω R d, then 2 u = div(grad u) = in Ω () is Laplace s equation defined in Ω. If d = 2, in cartesian
More informationTwo special equations: Bessel s and Legendre s equations. p Fourier-Bessel and Fourier-Legendre series. p
LECTURE 1 Table of Contents Two special equations: Bessel s and Legendre s equations. p. 259-268. Fourier-Bessel and Fourier-Legendre series. p. 453-460. Boundary value problems in other coordinate system.
More informationPreliminary Exam 2016 Solutions to Morning Exam
Preliminary Exam 16 Solutions to Morning Exam Part I. Solve four of the following five problems. Problem 1. Find the volume of the ice cream cone defined by the inequalities x + y + z 1 and x + y z /3
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 informationAero III/IV Conformal Mapping
Aero III/IV Conformal Mapping View complex function as a mapping Unlike a real function, a complex function w = f(z) cannot be represented by a curve. Instead it is useful to view it as a mapping. Write
More informationLet s estimate the volume under this surface over the rectangle R = [0, 4] [0, 2] in the xy-plane.
Math 54 - Vector Calculus Notes 3. - 3. Double Integrals Consider f(x, y) = 8 x y. Let s estimate the volume under this surface over the rectangle R = [, 4] [, ] in the xy-plane. Here is a particular estimate:
More informationLaplace s equation in polar coordinates. Boundary value problem for disk: u = u rr + u r r. r 2
Laplace s equation in polar coordinates Boundary value problem for disk: u = u rr + u r r + u θθ = 0, u(a, θ) = h(θ). r 2 Laplace s equation in polar coordinates Boundary value problem for disk: u = u
More informationIntegral Equations in Electromagnetics
Integral Equations in Electromagnetics Massachusetts Institute of Technology 6.635 lecture notes Most integral equations do not have a closed form solution. However, they can often be discretized and solved
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 informationMathematical Methods
Mathematical Methods Lecture 11. 1/12-28 1 Adjoint and self-adjoint operators We again define the adjoint operator L or L through (Lf, g) = (f, L g). We should think in the distribution-sense, so that
More informationExistence Theory: Green s Functions
Chapter 5 Existence Theory: Green s Functions In this chapter we describe a method for constructing a Green s Function The method outlined is formal (not rigorous) When we find a solution to a PDE by constructing
More informationM412 Assignment 5 Solutions
M41 Assignment 5 Solutions 1. [1 pts] Haberman.5.1 (a). Solution. Setting u(x, y) =X(x)Y(y), we find that X and Y satisfy X + λx = Y λy =, X () = X () = Y() =. In this case, we find from the X equation
More informationGeneral Technical Remarks on PDE s and Boundary Conditions Kurt Bryan MA 436
General Technical Remarks on PDE s and Boundary Conditions Kurt Bryan MA 436 1 Introduction You may have noticed that when we analyzed the heat equation on a bar of length 1 and I talked about the equation
More informationMATH 220: INNER PRODUCT SPACES, SYMMETRIC OPERATORS, ORTHOGONALITY
MATH 22: INNER PRODUCT SPACES, SYMMETRIC OPERATORS, ORTHOGONALITY When discussing separation of variables, we noted that at the last step we need to express the inhomogeneous initial or boundary data as
More informationLecture 13: Electromagnetic Theory Professor D. K. Ghosh, Physics Department, I.I.T., Bombay. Poisson s and Laplace s Equations
Poisson s and Laplace s Equations Lecture 13: Electromagnetic Theory Professor D. K. Ghosh, Physics Department, I.I.T., Bombay We will spend some time in looking at the mathematical foundations of electrostatics.
More informationBoundary value problems
1 Introduction Boundary value problems Lecture 5 We have found that the electric potential is a solution of the partial differential equation; 2 V = ρ/ǫ 0 The above is Poisson s equation where ρ is the
More informationcarries the circle w 1 onto the circle z R and sends w = 0 to z = a. The function u(s(w)) is harmonic in the unit circle w 1 and we obtain
4. Poisson formula In fact we can write down a formula for the values of u in the interior using only the values on the boundary, in the case when E is a closed disk. First note that (3.5) determines the
More informationModeling using conservation laws. Let u(x, t) = density (heat, momentum, probability,...) so that. u dx = amount in region R Ω. R
Modeling using conservation laws Let u(x, t) = density (heat, momentum, probability,...) so that u dx = amount in region R Ω. R Modeling using conservation laws Let u(x, t) = density (heat, momentum, probability,...)
More informationEXAM MATHEMATICAL METHODS OF PHYSICS. TRACK ANALYSIS (Chapters I-V). Thursday, June 7th,
EXAM MATHEMATICAL METHODS OF PHYSICS TRACK ANALYSIS (Chapters I-V) Thursday, June 7th, 1-13 Students who are entitled to a lighter version of the exam may skip problems 1, 8-11 and 16 Consider the differential
More informationBoundary Value Problems in Cylindrical Coordinates
Boundary Value Problems in Cylindrical Coordinates 29 Outline Differential Operators in Various Coordinate Systems Laplace Equation in Cylindrical Coordinates Systems Bessel Functions Wave Equation the
More informationare harmonic functions so by superposition
J. Rauch Applied Complex Analysis The Dirichlet Problem Abstract. We solve, by simple formula, the Dirichlet Problem in a half space with step function boundary data. Uniqueness is proved by complex variable
More informationCalculus II Practice Test 1 Problems: , 6.5, Page 1 of 10
Calculus II Practice Test Problems: 6.-6.3, 6.5, 7.-7.3 Page of This is in no way an inclusive set of problems there can be other types of problems on the actual test. To prepare for the test: review homework,
More informationMath 46, Applied Math (Spring 2009): Final
Math 46, Applied Math (Spring 2009): Final 3 hours, 80 points total, 9 questions worth varying numbers of points 1. [8 points] Find an approximate solution to the following initial-value problem which
More informationLegendre s Equation. PHYS Southern Illinois University. October 13, 2016
PHYS 500 - Southern Illinois University October 13, 2016 PHYS 500 - Southern Illinois University Legendre s Equation October 13, 2016 1 / 10 The Laplacian in Spherical Coordinates The Laplacian is given
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 informationc2 2 x2. (1) t = c2 2 u, (2) 2 = 2 x x 2, (3)
ecture 13 The wave equation - final comments Sections 4.2-4.6 of text by Haberman u(x,t), In the previous lecture, we studied the so-called wave equation in one-dimension, i.e., for a function It was derived
More informationPartial Differential Equations for Engineering Math 312, Fall 2012
Partial Differential Equations for Engineering Math 312, Fall 2012 Jens Lorenz July 17, 2012 Contents Department of Mathematics and Statistics, UNM, Albuquerque, NM 87131 1 Second Order ODEs with Constant
More informationChapter 12. Partial di erential equations Di erential operators in R n. The gradient and Jacobian. Divergence and rotation
Chapter 12 Partial di erential equations 12.1 Di erential operators in R n The gradient and Jacobian We recall the definition of the gradient of a scalar function f : R n! R, as @f grad f = rf =,..., @f
More information1 Distributions (due January 22, 2009)
Distributions (due January 22, 29). The distribution derivative of the locally integrable function ln( x ) is the principal value distribution /x. We know that, φ = lim φ(x) dx. x ɛ x Show that x, φ =
More informationELLIPTIC PARTIAL DIFFERENTIAL EQUATIONS EXERCISES I (HARMONIC FUNCTIONS)
ELLIPTIC PARTIAL DIFFERENTIAL EQUATIONS EXERCISES I (HARMONIC FUNCTIONS) MATANIA BEN-ARTZI. BOOKS [CH] R. Courant and D. Hilbert, Methods of Mathematical Physics, Vol. Interscience Publ. 962. II, [E] L.
More informationMy signature below certifies that I have complied with the University of Pennsylvania s Code of Academic Integrity in completing this exam.
My signature below certifies that I have complied with the University of Pennsylvania s Code of Academic Integrity in completing this exam. Signature Printed Name Math 241 Exam 1 Jerry Kazdan Feb. 17,
More informationLinear Operators, Eigenvalues, and Green s Operator
Chapter 10 Linear Operators, Eigenvalues, and Green s Operator We begin with a reminder of facts which should be known from previous courses. 10.1 Inner Product Space A vector space is a collection of
More informationFinite Difference Methods for Boundary Value Problems
Finite Difference Methods for Boundary Value Problems October 2, 2013 () Finite Differences October 2, 2013 1 / 52 Goals Learn steps to approximate BVPs using the Finite Difference Method Start with two-point
More informationNPTEL web course on Complex Analysis. A. Swaminathan I.I.T. Roorkee, India. and. V.K. Katiyar I.I.T. Roorkee, India
NPTEL web course on Complex Analysis A. Swaminathan I.I.T. Roorkee, India and V.K. Katiyar I.I.T. Roorkee, India A.Swaminathan and V.K.Katiyar (NPTEL) Complex Analysis 1 / 20 Complex Analysis Module: 10:
More informationAnnouncements. Topics: Homework:
Announcements Topics: - sections 7.5 (additional techniques of integration), 7.6 (applications of integration), * Read these sections and study solved examples in your textbook! Homework: - review lecture
More informationTD M1 EDP 2018 no 2 Elliptic equations: regularity, maximum principle
TD M EDP 08 no Elliptic equations: regularity, maximum principle Estimates in the sup-norm I Let be an open bounded subset of R d of class C. Let A = (a ij ) be a symmetric matrix of functions of class
More informationMATH 31BH Homework 5 Solutions
MATH 3BH Homework 5 Solutions February 4, 204 Problem.8.2 (a) Let x t f y = x 2 + y 2 + 2z 2 and g(t) = t 2. z t 3 Then by the chain rule a a a D(g f) b = Dg f b Df b c c c = [Dg(a 2 + b 2 + 2c 2 )] [
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 information(You may need to make a sin / cos-type trigonometric substitution.) Solution.
MTHE 7 Problem Set Solutions. As a reminder, a torus with radii a and b is the surface of revolution of the circle (x b) + z = a in the xz-plane about the z-axis (a and b are positive real numbers, with
More information51. General Surface Integrals
51. General urface Integrals The area of a surface in defined parametrically by r(u, v) = x(u, v), y(u, v), z(u, v) over a region of integration in the input-variable plane is given by d = r u r v da.
More informationSuggested Solution to Assignment 7
MATH 422 (25-6) partial diferential equations Suggested Solution to Assignment 7 Exercise 7.. Suppose there exists one non-constant harmonic function u in, which attains its maximum M at x. Then by the
More informationSeparation of Variables
Separation of Variables A typical starting point to study differential equations is to guess solutions of a certain form. Since we will deal with linear PDEs, the superposition principle will allow us
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 informationThe continuity method
The continuity method The method of continuity is used in conjunction with a priori estimates to prove the existence of suitably regular solutions to elliptic partial differential equations. One crucial
More informationPartial Differential Equations
Part II Partial Differential Equations Year 2015 2014 2013 2012 2011 2010 2009 2008 2007 2006 2005 2015 Paper 4, Section II 29E Partial Differential Equations 72 (a) Show that the Cauchy problem for u(x,
More informationIntegral Representation Formula, Boundary Integral Operators and Calderón projection
Integral Representation Formula, Boundary Integral Operators and Calderón projection Seminar BEM on Wave Scattering Franziska Weber ETH Zürich October 22, 2010 Outline Integral Representation Formula Newton
More informationElliptic Equations. Laplace equation on bounded domains Circular Domains
Elliptic Equtions Lplce eqution on bounded domins Sections 7.7.2, 7.7.3, An Introduction to Prtil Differentil Equtions, Pinchover nd Rubinstein 1.2 Circulr Domins We study the two-dimensionl Lplce eqution
More informationAn Introduction to Partial Differential Equations in the Undergraduate Curriculum
An Introduction to Partial Differential Equations in the Undergraduate Curriculum John C. Polking LECTURE 12 Heat Transfer in the Ball 12.1. Outline of Lecture The problem The problem for radial temperatures
More informationMath 251 December 14, 2005 Final Exam. 1 18pt 2 16pt 3 12pt 4 14pt 5 12pt 6 14pt 7 14pt 8 16pt 9 20pt 10 14pt Total 150pt
Math 251 December 14, 2005 Final Exam Name Section There are 10 questions on this exam. Many of them have multiple parts. The point value of each question is indicated either at the beginning of each question
More informationElliptic Equations. Chapter Definitions. Contents. 4.2 Properties of Laplace s and Poisson s Equations
5 4. Properties of Laplace s and Poisson s Equations Chapter 4 Elliptic Equations Contents. Neumann conditions the normal derivative, / = n u is prescribed on the boundar second BP. In this case we have
More informationPartial Differential Equations (TATA27)
Partial Differential Equations (TATA7) David Rule Spring 9 Contents Preliminaries. Notation.............................................. Differential equations...................................... 3
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 informationMATH 220: MIDTERM OCTOBER 29, 2015
MATH 22: MIDTERM OCTOBER 29, 25 This is a closed book, closed notes, no electronic devices exam. There are 5 problems. Solve Problems -3 and one of Problems 4 and 5. Write your solutions to problems and
More information36. Double Integration over Non-Rectangular Regions of Type II
36. Double Integration over Non-Rectangular Regions of Type II When establishing the bounds of a double integral, visualize an arrow initially in the positive x direction or the positive y direction. A
More informationIntroduction to the Calculus of Variations
236861 Numerical Geometry of Images Tutorial 1 Introduction to the Calculus of Variations Alex Bronstein c 2005 1 Calculus Calculus of variations 1. Function Functional f : R n R Example: f(x, y) =x 2
More informationIntegration by Substitution
November 22, 2013 Introduction 7x 2 cos(3x 3 )dx =? 2xe x2 +5 dx =? Chain rule The chain rule: d dx (f (g(x))) = f (g(x)) g (x). Use the chain rule to find f (x) and then write the corresponding anti-differentiation
More informationBoundary value problems for partial differential equations
Boundary value problems for partial differential equations Henrik Schlichtkrull March 11, 213 1 Boundary value problem 2 1 Introduction This note contains a brief introduction to linear partial differential
More informationConnection to Laplacian in spherical coordinates (Chapter 13)
Connection to Laplacian in spherical coordinates (Chapter 13) We might often encounter the Laplace equation and spherical coordinates might be the most convenient 2 u(r, θ, φ) = 0 We already saw in Chapter
More informationMATH 261 FINAL EXAM PRACTICE PROBLEMS
MATH 261 FINAL EXAM PRACTICE PROBLEMS These practice problems are pulled from the final exams in previous semesters. The 2-hour final exam typically has 8-9 problems on it, with 4-5 coming from the post-exam
More informationIntroduction to Partial Differential Equations
Introduction to Partial Differential Equations by Peter J Olver Page ix Corrections to First Printing (2014) Last updated: May 2, 2017 Line 7: Change In order the make progress to In order to make progress
More informationMATH-UA 263 Partial Differential Equations Recitation Summary
MATH-UA 263 Partial Differential Equations Recitation Summary Yuanxun (Bill) Bao Office Hour: Wednesday 2-4pm, WWH 1003 Email: yxb201@nyu.edu 1 February 2, 2018 Topics: verifying solution to a PDE, dispersion
More informationClass Meeting # 1: Introduction to PDEs
MATH 18.152 COURSE NOTES - CLASS MEETING # 1 18.152 Introduction to PDEs, Spring 2017 Professor: Jared Speck Class Meeting # 1: Introduction to PDEs 1. What is a PDE? We will be studying functions u =
More informationTopic: Introduction to Green s functions
Lecture notes on Variational and Approximate Methods in Applied Mathematics - A Peirce UBC 1 Topic: Introduction to Green s functions (Compiled 2 September 212) In this lecture we provide a brief introduction
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 informationFinal Exam May 4, 2016
1 Math 425 / AMCS 525 Dr. DeTurck Final Exam May 4, 2016 You may use your book and notes on this exam. Show your work in the exam book. Work only the problems that correspond to the section that you prepared.
More informationSturm-Liouville operators have form (given p(x) > 0, q(x)) + q(x), (notation means Lf = (pf ) + qf ) dx
Sturm-Liouville operators Sturm-Liouville operators have form (given p(x) > 0, q(x)) L = d dx ( p(x) d ) + q(x), (notation means Lf = (pf ) + qf ) dx Sturm-Liouville operators Sturm-Liouville operators
More informationMath 265 (Butler) Practice Midterm III B (Solutions)
Math 265 (Butler) Practice Midterm III B (Solutions). Set up (but do not evaluate) an integral for the surface area of the surface f(x, y) x 2 y y over the region x, y 4. We have that the surface are is
More informationOBJECTIVES Use the area under a graph to find total cost. Use rectangles to approximate the area under a graph.
4.1 The Area under a Graph OBJECTIVES Use the area under a graph to find total cost. Use rectangles to approximate the area under a graph. 4.1 The Area Under a Graph Riemann Sums (continued): In the following
More informationYou can learn more about the services offered by the teaching center by visiting
MAC 232 Exam 3 Review Spring 209 This review, produced by the Broward Teaching Center, contains a collection of questions which are representative of the type you may encounter on the exam. Other resources
More information6.3 Fundamental solutions in d
6.3. Fundamental solutions in d 5 6.3 Fundamental solutions in d Since Lapalce equation is invariant under translations, and rotations (see Exercise 6.4), we look for solutions to Laplace equation having
More information1 Integration in many variables.
MA2 athaye Notes on Integration. Integration in many variables.. Basic efinition. The integration in one variable was developed along these lines:. I f(x) dx, where I is any interval on the real line was
More informationVarious lecture notes for
Various lecture notes for 18306. R. R. Rosales (MIT, Math. Dept., 2-337) April 14, 2014 Abstract Notes for MIT s 18.306 Advanced PDE with applications. Contents 1 First lecture. 3 Tomography. Use filtering
More informationThe Maximum and Minimum Principle
MODULE 5: HEAT EQUATION 5 Lecture 2 The Maximum and Minimum Principle In this lecture, we shall prove the maximum and minimum properties of the heat equation. These properties can be used to prove uniqueness
More informationMATH 263 ASSIGNMENT 9 SOLUTIONS. F dv =
MAH AIGNMEN 9 OLUION ) Let F = (x yz)î + (y + xz)ĵ + (z + xy)ˆk and let be the portion of the cylinder x + y = that lies inside the sphere x + y + z = 4 be the portion of the sphere x + y + z = 4 that
More informationFinal: Solutions Math 118A, Fall 2013
Final: Solutions Math 118A, Fall 2013 1. [20 pts] For each of the following PDEs for u(x, y), give their order and say if they are nonlinear or linear. If they are linear, say if they are homogeneous or
More information