# u xx + u yy = 0. (5.1)

Size: px
Start display at page:

## Transcription

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 of the independent variables x, y only is called the Poisson equation. efinition 5.1 (Solution) (i) A twice continuously differentiable function φ : R is said to be a solution of the Laplace equation (5.1) if φ xx (x, y) + φ yy (x, y) = (x, y). (5.3) (ii) Let F be a given continuous function on. A twice continuously differentiable function φ : R is said to be a solution of the Poisson equation (5.) if φ xx (x, y) + φ yy (x, y) = F (x, y) (x, y). (5.4) Boundary Value Problems associated to Laplace Equation For domains other than the whole XY -plane, there is always a boundary and we can impose restrictions on the unknown function u or its derivatives or a mix of both on the boundary. These kinds of problems are called boundary value problems. There are at least three boundary value problems associated Laplace and Poisson equations. For simplicity, we describe them only for Laplace equation and the corresponding notions can be easily extended to Poisson equation. (1) irichlet Problem Let be a domain in R with piecewise smooth boundary and f be a continuous function on the closure of. To find a solution u of (5.1) in the domain such that (i) u is continuous on closure of and (ii) u = f If is the unit disk, then corresponding irichlet problem is called interior irichlet problem. If is the complement of the closed unit disk, then corresponding irichlet problem is called exterior irichlet problem. () Second Boundary Value Problem a.k.a. Neumann Problem To find a solution u of (5.1) in the domain such that (i) u, u x, u y are continuous on closure of and 41

2 4 5. (ii) at every point of, the directional derivative of u in the direction of normal (denoted by u n or nu) satisfies n u = f. (3) Third Boundary Value Problem a.k.a. Robin Problem To find a solution u of (5.1) in the domain such that (i) u, u x, u y are continuous on closure of and (ii) at every point of, u + α u n = f, where α is a given constant. Remark 5. (i) There are other kinds of boundary value problems possible. For example, on a part of the boundary one of the three boundary value problems described above is imposed and on the remaining part another one of the above three. We do not consider such boundary value problems here. (ii) Cauchy-Kowalewski theorem guarantees that a solution of an analytic Cauchy problem for an elliptic equation exists and is unique (locally), but is not always well-posed. Hadamard s Example We proved that Cauchy problem for the Wave equation is well-posed for x R and t [, T ] for every T >. Let us consider a similar problem for the Laplace equation now. More precisely, we consider the following Cauchy problem for Laplace equation in the upper half-plane: u = on R (, ), u(x, ) = f(x) for all x R, u y (x, ) = g(x) for all x R. (5.5a) (5.5b) (5.5c) Note that we proved the well-posedness of Cauchy problem for Wave equation, the initial conditions were exactly same as (5.5b)-(5.5c), except that the Laplace equation (5.5b) was replaced with Wave equation. Hadamard proved that Cauchy problem (5.5) is not well-posed on R [, T ] for any T >, by proving that stability estimate does not hold. For the moment, let us accept that solution to Cauchy problem exists and is unique 1. If the problem were to be well-posed, the following stability estimate is expected to be satisfied: Given ɛ >, there exists a δ > such that for every f 1, f, g 1, g satisfying f 1 (x) f (x) + g 1 (x) g (x) < δ x R, (5.6) the corresponding solutions u 1, u of the Cauchy problem satisfy u 1 (x, y) u (x, y) < ɛ x R, y. (5.7) What should we do to prove that solutions to Cauchy problem do not have stability? We produce a sequence (f n ), (g n ) which are close to the zero function but the corresponding solutions are far from the function u(x, y) which is the solution with zero Cauchy data. Consider f n, g n given by sin nx f n (x), g n (x) = n, (5.8) 1 This assumption says that if we somehow find a solution of the Cauchy problem then that is the only solution. Formulate the negation and convince yourself that what we are dong here is indeed equivalent.

3 Chapter 5 : Laplace Equation 43 note that they are very close to the function zero. The solution to Cauchy problem with the Cauchy data (5.8) is given by u n (x, y) = 1 sin(nx) sinh(ny). (5.9) n For n large, the initial conditions are very close to zero and hence can be thought of as a perturbation of the zero initial state. However the corresponding solution of Cauchy problem u n (x, y), due to the presence of hyperbolic sine function, is unbounded even on [R [, T ] for any T >. Summary : This example shows the difference between the Cauchy problems for Wave and Laplace equations posed on the upper half-plane. Even though the nature of the Cauchy data imposed is the same, changing the equation from Wave to Laplace changes the stability property drastically. 5.1 Green s identities Green s Identities form an important tool in the analysis of Laplace equation. They all are derived from divergence theorem, equivalently from Integration by parts formula. Let us recall the divergence theorem now. Theorem 5.3 (ivergence theorem) Let R be a bounded piecewise smooth domain. Let Ψ : R be a function. Let Ψ = (ψ 1, ψ ) where ψ i C 1 () C() for i = 1,. Then Ψ(x, y) dx dy = Ψ(x(s), y(s)) n(x(s), y(s)) ds, (5.1) where n(x(s), y(s)) is the unit outward normal at the point (x(s), y(s)). Green s Identity-I Let u : R be such that u C () C 1 (). Now apply (5.1) with Ψ = u to get Green s identity-i: u(x, y) dx dy = u(x(s), y(s)) n(x(s), y(s)) ds (5.11) With the notation n u := u n, Green s identity-i takes the form u(x, y) dx dy = n u ds (5.1) The quantity n u is called the normal derivative of u; as it represents the directional derivative of u in the dirction of the outward unit normal n. Green s Identity-II Let u : R be such that u C () C 1 (). Now apply (5.1) with Ψ = v u u v to get Green s identity-ii: ( ) ( v u u v (x, y) dx dy = v n u u n v ) ds (5.13) Green s Identity-III Let u : R be such that u C () C 1 (). Now apply (5.1) with Ψ = v u to get Green s identity-iii: u v dx dy = v n u ds v u dx dy. (5.14)

4 Green s identities Consequences of Green s identities Lemma 5.4 Let F C(). If u solves u = F on a domain, then F (x, y) dx dy = ν u dσ. (5.15) Proof: Integrating both sides of the equation u = F on yields F (x, y) dx dy = u(x, y) dx dy (5.16) Applying Green s identity-i (5.1), the integral on the right hand side of the equation (5.16) becomes u(x, y) dx dy = n u ds. (5.17) This finishes the proof. In particular, when u solves u =, we have n u dσ =. (5.18) Theorem 5.5 (Uniqueness theorem) Let be a smooth domain. Then (a) The irichlet problem has at most one solution. (b) If u solves the Neumann problem, then any other solution is of the form v = u + c for some real number c. (c) If α, then the Robin problem has at most one solution. Proof: Proof of (c): Note that (a) is a special case of (c), hence it is enough to prove (c). Let u 1, u be solutions of the Robin problem. efine w := u 1 u. Observe that w is a harmonic function(why?), and w satisfies the boundary condtion w + α n w =. (5.19) Using the Green s identity-iii (5.14) with u = v = w, we get w dx dy = α ( n w) ds. (5.) Since the left hand side of (5.) is non-negative while the right hand side is non-positive, we conclude that both sides must be zero. This implies that w = in and also n w = on. Since w = in, w must be a constant function. Since n w =, in view of the equality (5.19), we conclude that w = on the boundary. Since w C(), it follows that w in. This finishes the proof of (c). Proof of (b): Let u 1, u be solutions of the Neumann problem. efine w := u 1 u. Observe that w is a harmonic function(why?), and w satisfies the boundary condtion n w =. (5.1) Using the Green s identity-iii (5.14) with u = v = w, we get w dx dy =. (5.) This implies that w is a constant function. Thus u 1 = u + c for some constant c R.

5 Chapter 5 : Laplace Equation Weak Maximum Principle Let R be a bounded domain. Let u : R be a continuous function such that u can be extended to the closure of as a continuous function. Such class of functions is denoted by C(). When is bounded every function in C() attains both its maximum and minimum values, somewhere in. Weak Maximum principle says that if u happens to be a harmonic function, then the maximum and minimum values of u are definitely attained on the boundary of the domain whether or not they are attained in. Theorem 5.6 (The Weak Maxmimum Principle) Let R be a bounded domain. Let u C () C() be a harmonic function in. Then the maximum value of u in is achieved on the boundary. Proof: Step 1: Recall from differential calculus of two variables that at a point of interior maximum, u x and u y. As a consequence, u at an interior maximum point. Thus if v is a function such that v > in, then the maximum value of v on cannot be attained in. Hence v attains its maximum value only on the boundary. The idea to prove Weak maximum principle is to find such a function v starting from the given harmonic function u. Step : efine the function v ɛ by v ɛ (x, y) := u(x, y) + ɛ(x + y ). (5.3) Then v ɛ C () C(). Note that v ɛ > in and thus v ɛ attains its maximum only on the boundary. enoting M := Max u, L := Max (x + y ), (5.4) we have Since u(x, y) v ɛ (x, y) for all (x, y), we have v ɛ (x, y) M + ɛl, (x, y). (5.5) u(x, y) M + ɛl, (x, y). (5.6) Note that last inequality holds for every ɛ >. Thus taking limits as ɛ, we get This finishes the proof. u(x, y) M, (x, y). (5.7) Corollary 5.7 (The weak minimum principle) Let R be a bounded domain. Let u C () C() be a harmonic function in. Then the minimum value of u in is achieved on the boundary. Proof: efine by v(x, y) = u(x, y) on. Then v is a harmonic function. Now apply the weak maximum principle to the harmonic function v and conclude. Remark 5.8 We have to understand clearly what the Weak maximum principle says and what it does not talk about. The weak maximum principle says that on bounded domains, any harmonic function takes its maximum value on the boundary surely. The weak maximum principle is silent on whether the harmonic function will take or will not take the maximumum value in the domain. For example, consider to be the unit disk centered at the origin. Then the constant function u 1 is definitely a harmonic function in. It attains its maximum in also, apart from attaining on. There is another maximum principle, called Strong maximum priniciple, which says that on bounded domains non-constant harmonic functions will never attain maxima in and maxima are attained only on the boundary.

6 Weak Maximum Principle 5..1 Consequences of Weak Maximum Principle We already discussed uniqueness of solutions to all the three boundary value problems for Laplace equation posed on bounded domains using Green s identities. For Wave and Heat equations, we discussed uniqueness via Energy method and same analysis can be doen for Laplace equation as well. We are now going to discuss uniqueness questions concerning the boundary value problems using another tool, namely, the weak maximum principle for harmonic functions on a bounded domain. Theorem 5.9 (Uniqueness of solutions to irichlet problem) Let R be a bounded domain and consider the irichlet problem on : u = f on, u = g on, where g is a continuous function on. Then the irichlet problem has at most one solution in the class C () C(). Proof: Assume that u and v belonging to the class C () C() solve the irichlet problem. efine w := u v. Then w belongs to the class C () C() and solves the homogeneous irichlet problem: w = on, w = on. By weak maximum principle w and by weak minimum principle w. Thus w. This finishes the proof of the theorem. Remark 5.1 In the above theorem, it is essential that is a bounded domain. For unbounded domains the theorem may not hold good. For example, consider the following irichlet problem on the upper half plane: u(x, y) = for x R, < y <, u(x, ) = 1 for all x R. This problem has at least two solutions, namely u 1 (x, y) = xy, and, u (x, y) =. Theorem 5.11 Let R be a bounded domain. For i = 1,, let u i C () C() solve the following irichlet problems on : u i = f on, u = g i on, where g i is a continuous function on for each i = 1,. Then u 1 and u satisfy Max u 1 (x, y) u (x, y) Max g 1 (x, y) g (x, y). (5.8) Proof: efine w := u 1 u. Then w solves the following irichlet problem w = on, w = g 1 g on. Applying the weak maximum principle and the weak minimum principle, we get Min (g 1 g ) w(x, y) Max (g 1 g ) (x, y), (5.9) which finishes the proof of the theorem.

7 Chapter 5 : Laplace Equation Mean Value Property efinition 5.1 Let R be a domain. A continuous function u : R is said to possess the Mean Value Property-I on if for every point P = (x, y) and for every r > such that the open disk B r (P ), the function u satisfies the relation u(x, y) = 1 πr ξ x + η y r u(ξ, η) dξ dη. (5.3) efinition 5.13 Let R be a domain. A continuous function u : R is said to possess the Mean Value Property-II on if for every point P = (x, y) and for every r > such that the open disk B r (P ), the function u satisfies the relation u(x, y) = 1 πr where C r denotes the circle of radius r that is centered at P. C r u(s) ds (5.31) On expanding the line integral in the equation (5.31), we get the following form of the equation (5.31): u(x, y) = 1 π u(x + r cos θ, y + r sin θ) dθ. (5.3) π Lemma 5.14 Let u : R be a continuous function. equivalent. (i) u has the mean value property-i on. (ii) u has the mean value property-ii on. Then the following statements are Proof: We will show (ii) implies (i) and all the steps in this proof can be reversed thereby proving that (i) implies (ii). Let P = (x, y) be an arbitrary point. Further let R > be such that the open disk B R (P ). Then if u satisfies mean value property-ii, then for every < τ R we have u(x, y) = 1 π u(x + τ cos θ, y + τ sin θ) dθ. (5.33) π Multiply the equation (5.33) with τ and then integrate w.r.t. τ over the interval [, r]. This yields r τ u(x, y) dτ = 1 r ( π ) u(x + τ cos θ, y + τ sin θ) dθ τ dτ. (5.34) π The last equation simplifies to Thus we get r 1 u(x, y) = π = 1 π r π u(x, y) = 1 πr u(x + τ cos θ, y + τ sin θ) τ dθ dτ ξ x + η y r u(ξ, η) dξ dη. which proves that u has mean value property-i as well. ξ x + η y r u(ξ, η) dξ dη, (5.35) Remark 5.15 In view of Lemma 5.14, we say that a continuous function u has the mean value property on a domain if either of the mean value properties holds on and as a consequence both the mean value properties hold on. The mean value property-i may also be called the Solid Mean Value Property (as the average is taken over the entire disk) and the mean value property-ii may also be called Surface/Circle mean value theorem. This nomenclature is not standard.

8 Mean Value Property Theorem 5.16 (Mean Value Principle) Let R be a domain and u be a harmonic function in. Then u has the mean value property on. Proof: Let P = (x, y ) and let R > be such that the disk, denoted by B R (P ), of radius R centered at (x, y ) with the property B R (P ). We have to prove that, the function u satisfies the relation u(x, y ) = 1 π π u(x + R cos θ, y + R sin θ) dθ., (5.36) Note that the left hand side of (5.36) is a real number, whereas the right hand side of (5.36) depends on r. Thus to prove (5.36), the strategy is to prove that the right hand side is a constant function of the variable r. So, let us consider V (r) defined for r (, R] by Let us compute dv dr V (r) = 1 π π u(x + r cos θ, y + r sin θ) dθ. (5.37) for < r < R and prove that it is zero. As a consequence, we will then have V (r) = lim ρ V (ρ) = u(x, y ), (5.38) which completes the proof of the mean value principle. It remain to show that dv dr = on (, R). dv dr (r) = 1 π u(x + r cos θ, y + r sin θ) dθ = ν u dσ =. (5.39) π r In proving the last equality, we have used the equation (5.17) and the fact that u is a harmonic function on. Theorem 5.17 If u is continuous and has the mean value property on a domain, then it has continuous derivatives of all orders and all of them have the mean value property on. Proof: Step 1: It is sufficient to show that the first order partial derivatives of u exist and are continuous. Once it is known that first order partial derivatives of u exist, it follows by computing partial derivatives from the equation u(x, y) = 1 π π which holds as u satisfies mean value property on to obtain u(x + r cos θ, y + r sin θ) dθ (5.4) u x (x, y) = 1 π u y (x, y) = 1 π π π u x (x + r cos θ, y + r sin θ) dθ, u y (x + r cos θ, y + r sin θ) dθ. By repeating the arguments by replacing u with its partial derivatives, we conclude that u has partial derivatives of all orders and all of them satisfy the mean value relation. Step : Let us now prove that first order partial derivatives of u exist and are continuous. Let us recall the mean value relation (5.3) now. u(x, y) = 1 πr u(ξ, η) dξ dη. (5.41) ξ x + η y r

9 Chapter 5 : Laplace Equation 49 Writing the double integral on the right hand side as iterated integrals, the equation (5.41) takes the form u(x, y) = 1 πr ξ x + η y r u(ξ, η) dξ dη = 1 πr y+r y r x+ r (η y) dη u(ξ, η) dξ x r (η y) The last equation tells us that u x (x, y) exists, in view of fundamental theorem of integral calculus. In fact u x (x, y) is given by u x (x, y) = 1 πr y+r y r [ u(x + r (η y), η) u(x ] r (η y), η) dη. (5.4) Since the right hand side of the equation (5.4) is a continuous function of (x, y), we conclude that u x is a continuous function on. Similarly one can deduce the existence of u y at each point and continuity of u y on. Theorem 5.18 If u : R is continuous and u has the mean value property on, then u is harmonic in. Proof: If u has the mean value property, then by Theorem 5.17 we know that all partial derivatives of u exist and all of them are continuous on. We want to show that u =. On the contrary, suppose that there exists a point P in such that u(p ) ; and without loss of generality assume that u(p ) >. By continuity of the second order partial derivatives, there exists a disk of radius ɛ > centered at P (denote it by B ɛ (P ) on which u >. We now have < u(x, y) dx dy = ν u(s) ds B ɛ (P ) = ɛ ɛ C ɛ (P ) π = ɛ [πu(x, y)] =. ɛ u(x + ɛ cos θ, y + ɛ sin θ) dθ This contradiction means that our assumption that u at some point is non-zero is wrong. Thus u is a harmonic function Consequences of Mean Value Property Let R be a bounded domain and u : R be a harmonic function. The Weak maximum principle asserted that maximum value of u on the closure of is attained on the boundary of, and was silent about attaining the said maximum value in. Strong maximum principle that we are going to prove says that if such a maximum is also attained inside, then necessarily the harmonic function u is a constant function. Theorem 5.19 (Strong Maximum Principle) Let R be a domain (domain need not be bounded) and u : R be a harmonic function. If u attains its maximum in, then u is constant. Proof: Step 1: Let u assume its maximum at P and let this maximum value be denoted by M. We want to prove that u is the constant function that takes the value M everywhere in. Let P be an arbitrary point in. We wil show that u(p ) = M. Step : Join P to P by a smooth curve γ 3. Since γ is a compact set, it maintains a positive distance from (in case is non-empty), let us denote this distance by d γ. Take a disk of 3 This is possible since is open and connected and hence it is path-connected. Prove this.

10 Mean Value Property radius d γ with center at P denoted by B dγ (P ). Since u is a harmonic function in, it has the mean value property in. Applying the mean value property 4 on this disk, we conclude that u equals M on the entire disk. Step 3:(Continuation argument) We now take a point P 1 on γ which lies in the disk B dγ (P ) which is at least at a distance d γ 4 from P. Note that u(p 1 ) = M. Now repeat the arguments of Step with the disk B dγ get a k N with the property that P B dγ the proof of the theorem. (P 1 ) and then find a P in a similar way. Continuing this process till we (P k ) 5. We will then have u(p ) = M. This finishes Remark 5. This is a continuation of Remark 5.8. It is important to note that Strong maximum principle says that if a harmonic function attains its maximum in a domain (bounded domain or otherwise), then it is necessarily a constant function. In particular Strong maximum principle does not talk about where the maximum will be attained; sometimes a harmonic function may not have a maximum value also. This happens when the domain is not bounded, because in that case though u is continuous on we cannot assert the existence of a maximum value, since is not compact. See the next example in this context. Example 5.1 Let be the domain exterior to the unit disk. The function u(x, y) = log(x +y ) is a harmonic function on. It has neither a maximum value nor a minimum value in. Exercise 5. Formulate the strong minimum principle and prove it by following the method of proof of the Strong maximum principle. 4 Mean Value Property-I or Mean Value Property-II? 5 Why is it possible to find such a k?

### Laplace equation. In this chapter we consider Laplace equation in d-dimensions given by. + u x2 x u xd x d. u x1 x 1

Chapter 6 Laplace equation In this chapter we consider Laplace equation in d-dimensions given by u x1 x 1 + u x2 x 2 + + u xd x d =. (6.1) We study Laplace equation in d = 2 throughout this chapter (excepting

### Applications of the Maximum Principle

Jim Lambers MAT 606 Spring Semester 2015-16 Lecture 26 Notes These notes correspond to Sections 7.4-7.6 in the text. Applications of the Maximum Principle The maximum principle for Laplace s equation is

### Math 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

### GREEN 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

### In 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

### MATH 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

### MATH 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

### Partial differential equation for temperature u(x, t) in a heat conducting insulated rod along the x-axis is given by the Heat equation:

Chapter 7 Heat Equation Partial differential equation for temperature u(x, t) in a heat conducting insulated rod along the x-axis is given by the Heat equation: u t = ku x x, x, t > (7.1) Here k is a constant

### Laplace 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

### MATH 124B Solution Key HW 05

7.1 GREEN S FIRST IENTITY MATH 14B Solution Key HW 05 7.1 GREEN S FIRST IENTITY 1. erive the 3-dimensional maximum principle from the mean value property. SOLUTION. We aim to prove that if u is harmonic

### The 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

### Module 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

### Math 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

### Spotlight 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

### SOLUTIONS 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

### NATIONAL UNIVERSITY OF SINGAPORE Department of Mathematics MA4247 Complex Analysis II Lecture Notes Part II

NATIONAL UNIVERSITY OF SINGAPORE Department of Mathematics MA4247 Complex Analysis II Lecture Notes Part II Chapter 2 Further properties of analytic functions 21 Local/Global behavior of analytic functions;

### 21 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

### MATH 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

### Analysis III Solutions - Serie 12

.. Necessary condition Let us consider the following problem for < x, y < π, u =, for < x, y < π, u y (x, π) = x a, for < x < π, u y (x, ) = a x, for < x < π, u x (, y) = u x (π, y) =, for < y < π. Find

### u( x) = g( y) ds y ( 1 ) U solves u = 0 in U; u = 0 on U. ( 3)

M ath 5 2 7 Fall 2 0 0 9 L ecture 4 ( S ep. 6, 2 0 0 9 ) Properties and Estimates of Laplace s and Poisson s Equations In our last lecture we derived the formulas for the solutions of Poisson s equation

### MATH 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

### Suggested 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

### A proof for the full Fourier series on [ π, π] is given here.

niform convergence of Fourier series A smooth function on an interval [a, b] may be represented by a full, sine, or cosine Fourier series, and pointwise convergence can be achieved, except possibly at

### Metric Spaces and Topology

Chapter 2 Metric Spaces and Topology From an engineering perspective, the most important way to construct a topology on a set is to define the topology in terms of a metric on the set. This approach underlies

### General 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

### Partial 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,

### = 2 x y 2. (1)

COMPLEX ANALYSIS PART 5: HARMONIC FUNCTIONS A Let me start by asking you a question. Suppose that f is an analytic function so that the CR-equation f/ z = 0 is satisfied. Let us write u and v for the real

### Calculus - II Multivariable Calculus. M.Thamban Nair. Department of Mathematics Indian Institute of Technology Madras

Calculus - II Multivariable Calculus M.Thamban Nair epartment of Mathematics Indian Institute of Technology Madras February 27 Revised: January 29 Contents Preface v 1 Functions of everal Variables 1 1.1

### u(0) = u 0, u(1) = u 1. To prove what we want we introduce a new function, where c = sup x [0,1] a(x) and ɛ 0:

6. Maximum Principles Goal: gives properties of a solution of a PDE without solving it. For the two-point boundary problem we shall show that the extreme values of the solution are attained on the boundary.

### 6.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

### Final: 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

### The 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

### Preliminary Exam 2018 Solutions to Morning Exam

Preliminary Exam 28 Solutions to Morning Exam Part I. Solve four of the following five problems. Problem. Consider the series n 2 (n log n) and n 2 (n(log n)2 ). Show that one converges and one diverges

### Simple 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

### APPM 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

### Partial 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

### Existence and Uniqueness

Chapter 3 Existence and Uniqueness An intellect which at a certain moment would know all forces that set nature in motion, and all positions of all items of which nature is composed, if this intellect

### Minimal Surfaces: Nonparametric Theory. Andrejs Treibergs. January, 2016

USAC Colloquium Minimal Surfaces: Nonparametric Theory Andrejs Treibergs University of Utah January, 2016 2. USAC Lecture: Minimal Surfaces The URL for these Beamer Slides: Minimal Surfaces: Nonparametric

### Course 212: Academic Year Section 1: Metric Spaces

Course 212: Academic Year 1991-2 Section 1: Metric Spaces D. R. Wilkins Contents 1 Metric Spaces 3 1.1 Distance Functions and Metric Spaces............. 3 1.2 Convergence and Continuity in Metric Spaces.........

### Math 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

### n 2 xi = x i. x i 2. r r ; i r 2 + V ( r) V ( r) = 0 r > 0. ( 1 1 ) a r n 1 ( 1 2) V( r) = b ln r + c n = 2 b r n 2 + c n 3 ( 1 3)

Sep. 7 The L aplace/ P oisson Equations: Explicit Formulas In this lecture we study the properties of the Laplace equation and the Poisson equation with Dirichlet boundary conditions through explicit representations

### Math 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

### z x = f x (x, y, a, b), z y = f y (x, y, a, b). F(x, y, z, z x, z y ) = 0. This is a PDE for the unknown function of two independent variables.

Chapter 2 First order PDE 2.1 How and Why First order PDE appear? 2.1.1 Physical origins Conservation laws form one of the two fundamental parts of any mathematical model of Continuum Mechanics. These

### Introduction 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.

### SINC PACK, and Separation of Variables

SINC PACK, and Separation of Variables Frank Stenger Abstract This talk consists of a proof of part of Stenger s SINC-PACK computer package (an approx. 400-page tutorial + about 250 Matlab programs) that

### SOLUTION OF THE DIRICHLET PROBLEM WITH A VARIATIONAL METHOD. 1. Dirichlet integral

SOLUTION OF THE DIRICHLET PROBLEM WITH A VARIATIONAL METHOD CRISTIAN E. GUTIÉRREZ FEBRUARY 3, 29. Dirichlet integral Let f C( ) with open and bounded. Let H = {u C ( ) : u = f on } and D(u) = Du(x) 2 dx.

### α u x α1 1 xα2 2 x αn n

Chapter 1 Introduction In this chapter we introduce the notions of a Partial Differential Equation (PDE) and its solution. We attempt classifying all partial differential equations in at least three different

### Cauchy Integral Formula Consequences

Cauchy Integral Formula Consequences Monday, October 28, 2013 1:59 PM Homework 3 due November 15, 2013 at 5 PM. Last time we derived Cauchy's Integral Formula, which we will present in somewhat generalized

### Green 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

### 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

### Introduction to Topology

Introduction to Topology Randall R. Holmes Auburn University Typeset by AMS-TEX Chapter 1. Metric Spaces 1. Definition and Examples. As the course progresses we will need to review some basic notions about

### CLASSIFICATION 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

### Lucio Demeio Dipartimento di Ingegneria Industriale e delle Scienze Matematiche

Scuola di Dottorato THE WAVE EQUATION Lucio Demeio Dipartimento di Ingegneria Industriale e delle Scienze Matematiche Lucio Demeio - DIISM wave equation 1 / 44 1 The Vibrating String Equation 2 Second

### CHAPTER 2. Laplace s equation

18 CHAPTER 2 Laplace s equation There can be but one option as to the beauty and utility of this analysis by Laplace; but the manner in which it has hitherto been presented has seemed repulsive to the

### Sequences. Chapter 3. n + 1 3n + 2 sin n n. 3. lim (ln(n + 1) ln n) 1. lim. 2. lim. 4. lim (1 + n)1/n. Answers: 1. 1/3; 2. 0; 3. 0; 4. 1.

Chapter 3 Sequences Both the main elements of calculus (differentiation and integration) require the notion of a limit. Sequences will play a central role when we work with limits. Definition 3.. A Sequence

### THE WAVE EQUATION. d = 1: D Alembert s formula We begin with the initial value problem in 1 space dimension { u = utt u xx = 0, in (0, ) R, (2)

THE WAVE EQUATION () The free wave equation takes the form u := ( t x )u = 0, u : R t R d x R In the literature, the operator := t x is called the D Alembertian on R +d. Later we shall also consider the

### Lecture 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

### Diffusion on the half-line. The Dirichlet problem

Diffusion on the half-line The Dirichlet problem Consider the initial boundary value problem (IBVP) on the half line (, ): v t kv xx = v(x, ) = φ(x) v(, t) =. The solution will be obtained by the reflection

### FOURIER METHODS AND DISTRIBUTIONS: SOLUTIONS

Centre for Mathematical Sciences Mathematics, Faculty of Science FOURIER METHODS AND DISTRIBUTIONS: SOLUTIONS. We make the Ansatz u(x, y) = ϕ(x)ψ(y) and look for a solution which satisfies the boundary

### Numerical 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

### Aero 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

### An Introduction to Complex Analysis and Geometry John P. D Angelo, Pure and Applied Undergraduate Texts Volume 12, American Mathematical Society, 2010

An Introduction to Complex Analysis and Geometry John P. D Angelo, Pure and Applied Undergraduate Texts Volume 12, American Mathematical Society, 2010 John P. D Angelo, Univ. of Illinois, Urbana IL 61801.

### Integral 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

### LECTURE 15: COMPLETENESS AND CONVEXITY

LECTURE 15: COMPLETENESS AND CONVEXITY 1. The Hopf-Rinow Theorem Recall that a Riemannian manifold (M, g) is called geodesically complete if the maximal defining interval of any geodesic is R. On the other

### The Divergence Theorem Stokes Theorem Applications of Vector Calculus. Calculus. Vector Calculus (III)

Calculus Vector Calculus (III) Outline 1 The Divergence Theorem 2 Stokes Theorem 3 Applications of Vector Calculus The Divergence Theorem (I) Recall that at the end of section 12.5, we had rewritten Green

### Chapter 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,

### MATH 126 FINAL EXAM. Name:

MATH 126 FINAL EXAM Name: Exam policies: Closed book, closed notes, no external resources, individual work. Please write your name on the exam and on each page you detach. Unless stated otherwise, you

### UNIVERSITY OF MANITOBA

Question Points Score INSTRUCTIONS TO STUDENTS: This is a 6 hour examination. No extra time will be given. No texts, notes, or other aids are permitted. There are no calculators, cellphones or electronic

### Lecture Introduction

Lecture 1 1.1 Introduction The theory of Partial Differential Equations (PDEs) is central to mathematics, both pure and applied. The main difference between the theory of PDEs and the theory of Ordinary

### 4.1 Analysis of functions I: Increase, decrease and concavity

4.1 Analysis of functions I: Increase, decrease and concavity Definition Let f be defined on an interval and let x 1 and x 2 denote points in that interval. a) f is said to be increasing on the interval

### Analysis II - few selective results

Analysis II - few selective results Michael Ruzhansky December 15, 2008 1 Analysis on the real line 1.1 Chapter: Functions continuous on a closed interval 1.1.1 Intermediate Value Theorem (IVT) Theorem

### Numerical Solutions to Partial Differential Equations

Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Discretization of Boundary Conditions Discretization of Boundary Conditions On

### M.Sc. in Meteorology. Numerical Weather Prediction

M.Sc. in Meteorology UCD Numerical Weather Prediction Prof Peter Lynch Meteorology & Climate Centre School of Mathematical Sciences University College Dublin Second Semester, 2005 2006. In this section

### Theory 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

### ELLIPTIC 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.

### 4.1. Dirichlet, Poisson and Neumann boundary value problems

Chapter 4 Laplace s equation 4.. Dirichlet, Poisson and Neumann boundary value problems The most commonly occurring form of problem that is associated with Laplace s equation is a boundary value problem,

### converges as well if x < 1. 1 x n x n 1 1 = 2 a nx n

Solve the following 6 problems. 1. Prove that if series n=1 a nx n converges for all x such that x < 1, then the series n=1 a n xn 1 x converges as well if x < 1. n For x < 1, x n 0 as n, so there exists

### MATH 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,

### Metric spaces and metrizability

1 Motivation Metric spaces and metrizability By this point in the course, this section should not need much in the way of motivation. From the very beginning, we have talked about R n usual and how relatively

### PROBLEM SET 6. E [w] = 1 2 D. is the smallest for w = u, where u is the solution of the Neumann problem. u = 0 in D u = h (x) on D,

PROBLEM SET 6 UE ATE: - APR 25 Chap 7, 9. Questions are either directl from the text or a small variation of a problem in the text. Collaboration is oka, but final submission must be written individuall.

### f (n) (z 0 ) Theorem [Morera s Theorem] Suppose f is continuous on a domain U, and satisfies that for any closed curve γ in U, γ

Remarks. 1. So far we have seen that holomorphic is equivalent to analytic. Thus, if f is complex differentiable in an open set, then it is infinitely many times complex differentiable in that set. This

### ξ,i = x nx i x 3 + δ ni + x n x = 0. x Dξ = x i ξ,i = x nx i x i x 3 Du = λ x λ 2 xh + x λ h Dξ,

1 PDE, HW 3 solutions Problem 1. No. If a sequence of harmonic polynomials on [ 1,1] n converges uniformly to a limit f then f is harmonic. Problem 2. By definition U r U for every r >. Suppose w is a

### Method 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

### McGill University Math 354: Honors Analysis 3

Practice problems McGill University Math 354: Honors Analysis 3 not for credit Problem 1. Determine whether the family of F = {f n } functions f n (x) = x n is uniformly equicontinuous. 1st Solution: The

### Various 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

### Functional Analysis. Franck Sueur Metric spaces Definitions Completeness Compactness Separability...

Functional Analysis Franck Sueur 2018-2019 Contents 1 Metric spaces 1 1.1 Definitions........................................ 1 1.2 Completeness...................................... 3 1.3 Compactness......................................

### INTRODUCTION 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

### Notions such as convergent sequence and Cauchy sequence make sense for any metric space. Convergent Sequences are Cauchy

Banach Spaces These notes provide an introduction to Banach spaces, which are complete normed vector spaces. For the purposes of these notes, all vector spaces are assumed to be over the real numbers.

### (x x 0 ) 2 + (y y 0 ) 2 = ε 2, (2.11)

2.2 Limits and continuity In order to introduce the concepts of limit and continuity for functions of more than one variable we need first to generalise the concept of neighbourhood of a point from R to

### 18 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

### Math 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

### FOURIER TRANSFORMS. 1. Fourier series 1.1. The trigonometric system. The sequence of functions

FOURIER TRANSFORMS. Fourier series.. The trigonometric system. The sequence of functions, cos x, sin x,..., cos nx, sin nx,... is called the trigonometric system. These functions have period π. The trigonometric

### Topological properties

CHAPTER 4 Topological properties 1. Connectedness Definitions and examples Basic properties Connected components Connected versus path connected, again 2. Compactness Definition and first examples Topological

### McGill 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

### INTEGRABILITY OF SUPERHARMONIC FUNCTIONS IN A JOHN DOMAIN. Hiroaki Aikawa

INTEGRABILITY OF SUPERHARMONIC FUNCTIONS IN A OHN OMAIN Hiroaki Aikawa Abstract. The integrability of positive erharmonic functions on a bounded fat ohn domain is established. No exterior conditions are

### [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,

### 1. (a) (5 points) Find the unit tangent and unit normal vectors T and N to the curve. r (t) = 3 cos t, 0, 3 sin t, r ( 3π

1. a) 5 points) Find the unit tangent and unit normal vectors T and N to the curve at the point P 3, 3π, r t) 3 cos t, 4t, 3 sin t 3 ). b) 5 points) Find curvature of the curve at the point P. olution:

### Introduction to Real Analysis Alternative Chapter 1

Christopher Heil Introduction to Real Analysis Alternative Chapter 1 A Primer on Norms and Banach Spaces Last Updated: March 10, 2018 c 2018 by Christopher Heil Chapter 1 A Primer on Norms and Banach Spaces