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
|
|
- Rafe Fields
- 5 years ago
- Views:
Transcription
1 Chapter 6 Laplace equation In this chapter we consider Laplace equation in d-dimensions given by u x1 x 1 + u x2 x u xd x d =. (6.1) We study Laplace equation in d = 2 throughout this chapter (excepting Section 6.2), and most of the ideas can be generalized to general space dimensions d > 2. When d = 2, the independent variables x 1, x 2 are denoted by x, y, and write x = (x, y). Thus Laplace equation in two independent variables is The non-homogeneous problem u x x + u yy =. (6.2) u xx + u yy = f, (6.3) where f is a function of the independent variables x, y only is called the Poisson equation. 6.1 Green s identities Green s Identities play an important role in the analysis of Laplace equation. They are derived from divergence theorem. Let us recall the divergence theorem now. Theorem 6.1 (Divergence theorem). Let 2 be a bounded piecewise smooth domain. Let Ψ : 2 be a function. Let Ψ = (ψ 1,ψ 2 ) where ψ i C 1 () C () for i = 1,2. Then.Ψ(x, y) d x d y = Ψ.n dσ, (6.4) where n is the unit outward normal to. Theorem 6.2 (Green s identities). Let 2 be a bounded piecewise smooth domain, and let n denote the unit outward normal to. For u, v C 2 () C 1 (), we have (i) (Green s Identity-I) u(x, y) d x d y = n u dσ, (6.5) 147
2 148 Chapter 6. Laplace equation where n u := u.n is called the normal derivative of u; as it represents the directional derivative of u in the dirction of the unit outward normal n. (ii) (Green s Identity-II) v u u v (x, y) d x d y = v n u u n v dσ (6.6) (iii) (Green s Identity-III) u. v d x d y = v n u dσ v u d x d y. (6.7) Proof. Applying (6.4) with Ψ = u, Ψ = v u u v, and Ψ = v u yield Green s identities I, II, and III respectively. 6.2 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 such symmetric properties. Let us fix ξ d, and look for solutions to u = having the form v ξ (x) = ψ(r ), where r = x ξ = d (x i ξ i ) 2. i=1 Substituting the formula for v ξ in the equation, we get Solving the last equation, we get v ξ (x) = ψ (r ) + d 1 ψ (r ) =. r ψ (r ) = C r 1 d (6.8) Integrating we get C ln r if d = 2, ψ(r ) = C 2 d r 2 d if d 3. (6.9) Theorem 6.3. Let K(x, ξ ) denote the function defined by K(x, ξ ) := ψ( x ξ ). (i) Let ξ be a point of. For u C 2 () the following identity holds u(ξ ) = K(x,ξ ) u d x (K(x,ξ ) n u u n K(x,ξ )) dσ. (6.1) (ii) The following equality holds in the sense of distributions on d : K(x,ξ ) = δ ξ. i.e., for every φ C (d ) the following equality holds. φ(ξ ) = K(x,ξ ) φ(x) d x. (6.11) d IIT Bombay
3 6.2. Fundamental solutions in d 149 (iii) If u C 2 () and harmonic in (i.e., u = in ), then for ξ we get u(ξ ) = (K(x,ξ ) n u u n K(x,ξ )) dσ. (6.12) Proof. Step 1: Proof of (i): Let u C 2 () and ξ be a point of. Note that we cannot apply Green s identity II (6.6) with v = v ξ (x) = ψ( x ξ ) since v ξ is singular at x = ξ. Thus we cut out a ball B(ξ,ρ) from along with its boundary, and then apply Green s identity II. Note that the domain ρ ;= \ B[ξ,ρ] is bounded by and S(ξ,ρ). Since v ξ = in ρ, we have v ξ u d x = vξ n u u n v ξ dσ + vξ n u u n v ξ dσ. (6.13) ρ S(ξ,ρ) Let us now compute the second term on RHS of the equation (6.13). vξ n u u n v ξ dσ = v ξ n u dσ u n v ξ dσ (6.14) S(ξ,ρ) S(ξ,ρ) S(ξ,ρ) Note that on the circle S(ξ,ρ), we have v ξ (x) = ψ( x ξ ) = ψ(ρ). Using this information, and divergence theorem, we get v ξ n u dσ = ψ(ρ) n u dσ = ψ(ρ) u d x. (6.15) S(ξ,ρ) S(ξ,ρ) B(ξ,ρ) Note that the outward unit normal n on the circle S(ξ,ρ) points towards its center ξ. Also n v ξ = ψ (ρ) holds at points on the circle S(ξ,ρ). Thus we get S(ξ,ρ) u n v ξ dσ = C ρ 1 d S(ξ,ρ) u dσ. (6.16) Since both u and u are continuous at ξ, we have the following convergences as ρ : ψ(ρ) u d x, ρ 1 d u dσ ω d u(ξ ), (6.17) B(ξ,ρ) S(ξ,ρ) where ω d denotes the surface area of the unit sphere in d. Now passing to the limit as ρ in the equation (6.13), we get v ξ u d x = vξ n u u n v ξ dσ + C ωd u(ξ ) (6.18) in view of (6.15), (6.16), (6.17). Choosing C = 1 ω d in the formulae (6.9), the equation (6.18) yields (6.1). Step 2: Proof of (ii): In the equation (6.1) if we take u = φ C (6.11). (), then we get Step 3: Proof of (iii): Equation (6.12) is an immediate consequence of (6.1). October 4, 215
4 15 Chapter 6. Laplace equation 6.3 Boundary value problems associated to Laplace equation We can study Laplace equation or Poisson equation on any open domain 2. These problems will always have an infinite number of solutions. Thus it is interesting to know if there are solutions to these equations if some constraints are placed on the unknown function. When is a bounded domain with piecewise smooth boundary, there are atleast three different kinds of restrictions are placed on the unknown function and its derivatives along the boundary of the domain. Such restrictions are called boundary conditions. Let us now describe three boundary value problems (BVP) for Poisson equation. They differ in the nature of boundary conditions. We also define the notions of solutions to these boundary value problems. 1. For given functions f, g, Dirichlet Problem consists of solving the boundary value problem u = f (x) in, (6.19a) u = g on. (6.19b) Definition 6.4 (Solution to Dirichlet BVP). Let f C (), and g C ( ). A function φ C 2 () C () is said to be a solution to Dirichlet BVP if (i) the function φ is a solution to Poisson equation (6.19a), i.e., for each x, φ(x) = f (x) holds, and (ii) for each x, the equality φ(x) = g(x) holds. 2. For given functions f, g, Neumann Problem consists of solving the boundary value problem u = f (x) in, (6.2a) n u = g on. (6.2b) Definition 6.5 (Solution to Neumann BVP). Let f C (), and g C ( ). A function φ C 2 () C 1 () is said to be a solution to Neumann BVP if (i) the function φ is a solution to Poisson equation (6.2a), i.e., for each x, φ(x) = f (x) holds, and (ii) for each x, the equality n φ(x) = g(x) holds. 3. For given functions f, g, and α, Third Boundary Value Problem a.k.a. Robin Problem consists of solving the boundary value problem u = f (x) in, (6.21a) u + α n u = g on. (6.21b) Definition 6.6 (Solution to Robin BVP). Let f C (), and g C ( ). A function φ C 2 () C 1 () is said to be a solution to Robin BVP if IIT Bombay
5 6.3. Boundary value problems associated to Laplace equation 151 (i) the function φ is a solution to Poisson equation (6.21a), i.e., for each x, φ(x) = f (x) holds, and (ii) for each x, the equality φ(x) + α n φ(x) = g(x) holds. Remark 6.7. (i) If Dirichlet BVP is posed on a bounded domain, then the corresponding Dirichlet problem is called an interior Dirichlet problem. If is the complement of a bounded domain, then corresponding Dirichlet problem is called exterior Dirichlet problem. (ii) 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. (iii) 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 wellposed. The following result says that for Neumann BVP to admit a solution, the data f, g must be compatible with each other. Lemma 6.8 (Compatibility of data for Neumann problem). Let f C (). If u C 2 () is a solution to Neumann BVP on, then f (x) d x = g(y) dσ(y). (6.22) Proof. Integrating both sides of the equation u = f on yields f (x) d x = u(x) d x. (6.23) Applying Green s identity-i (6.5), the integral on the right hand side of the equation (6.23) becomes u(x) d x = n u(y) dσ(y). (6.24) This completes the proof. In particular, when u is a solution to Neumann BVP with f =, we get g(y) dσ(y) =. Theorem 6.9 (Uniqueness theorem). Let be a smooth domain. Then (a) The Dirichlet problem has at most one solution. October 4, 215
6 152 Chapter 6. Laplace equation (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 2 be solutions of the Robin problem. Define w := u 1 u 2. Observe that w is a harmonic function(why?), and w satisfies the boundary condtion Using the Green s identity-iii (6.7) with u = v = w, we get w + α n w =. (6.25) w 2 d x d y = α ( n w) 2 d s. (6.26) Since the left hand side of (6.26) 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 (6.25), 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 2 be solutions of the Neumann problem. Define w := u 1 u 2. Observe that w is a harmonic function(why?), and w satisfies the boundary condtion Using the Green s identity-iii (6.7) with u = v = w, we get n w =. (6.27) w 2 d x d y =. (6.28) This implies that w is a constant function. Thus u 1 = u 2 + c for some constant c. Remark 6.1 (On the formula (6.12)). (i) If Dirichlet problem (6.19) has a solution, then by Theorem 6.9 Dirichlet problem has exactly one solution. Let the solution to Dirichlet problem be denoted by u. The formula (6.12) gives a representation of the solution in terms of the boundary values of u and its normal derivative n u on. (ii) Note however that in Dirichlet problem, only the values of u are prescribed on, and n u is an unknown function on, and thus the formula (6.12) is not useful for computing the solution. (iii) Note that boundary values of u already determine a solution to Dirichlet problem, and thus the quantity n u is already determined. (iv) A related question concerning Laplace equation is whether Cauchy problem is wellposed for Laplace equation. The answer to this question is negative, and is made precise in Theorem 6.31 IIT Bombay
7 6.4. Poisson s formula and its applications Poisson s formula and its applications In this section, we consider Dirichlet boundary value problem for Laplace equation on the domain = B(,1) = {(x 1, x 2 ) : x x2 2 < 1} in 2. For a given function f C (S(,1)), we consider the boundary value problem We are going to prove the following result. u = in B(,1), (6.29a) u = f on S(,1). (6.29b) Theorem 6.11 (Poisson s formula). The solution to the BVP (6.29) is given by the Poisson s formula u(x) = 1 x 2 f (y) d y. (6.3) 2π x y 2 S(,1) Defining v(r,θ) := u(r cosθ, r sinθ), and F (θ) := f (cosθ,sinθ), Poisson s formula in polar coordinates is given by v(r,θ) = 1 r 2 2π F (τ) dτ (6.31) 2π 1 2r cos(τ θ) + r 2 under the condition that F C 1 [,2π] and is a 2π-periodic function. Proof. Due to the geometry of disk, it is convenient to consider Laplace equation expressed in polar coordinates (see Exercise 6.1). Thus we consider the Dirichlet problem given by (r,θ) v v r r + 1 r v r + 1 r v 2 θθ = for < r < 1, θ 2π, (6.32a) v(1,θ) = F (θ) for θ 2π, (6.32b) The proof of theorem is divided into three main steps. In Step 1, we derive a formal solution to Dirichlet BVP. In Step 2, we justify that the formal solution derived in the first step is indded a solution. In Step 3, Poisson s formula will be obtained. Step 1: Derivation of a formal solution We obtain a solution of the equation (6.32a) using the method of separation of variables, in which solution is assumed to be of the form v(r,θ) = h(r )g(θ), (6.33) where g is assumed to be periodic of period 2π i.e., g() = g(2π) and g () = g (2π). Substituting the expression for v into (6.32a), and then dividing the resultant equation with h(r )g(θ) we get h (r ) h(r ) + 1 r Re-arranging terms in (6.34) yields h (r ) h(r ) + 1 g (θ) r 2 g(θ) = (6.34) r 2 h (r ) h(r ) + r h (r ) h(r ) = g (θ) g(θ) = λ (6.35) October 4, 215
8 154 Chapter 6. Laplace equation Note that the LHS of the equation (6.35) is a function of r only, while the RHS is a function of θ only. Thus both of them must be equal to a constant, and let us denote it by λ. Thus we have r 2 h (r ) h(r ) + r h (r ) h(r ) = g (θ) g(θ) = λ (6.36) Thus we get two ODEs from (6.36) which are given by g (θ) + λg(θ) =, r 2 h (r ) + r h (r ) = λh(r ) (6.37a) (6.37b) Let us look for solutions of (6.37a) satisfying the periodic boundary conditions g() = g(2π), g () = g (2π) (6.38) We are interested in finding non-trivial solutions to the boundary value problem (6.32), existence of which depend on the value of λ. Any λ for which the BVP (6.32) admits a non-trivial solution is called an eigenvalue and the corresponding non-trivial solutions are called eigenfunctions. Thus we are interested in finding non-trivial solutions to ODEs (6.37a) with boundary conditions (6.38). (i) (λ = ) General solution of the ODE (6.37a) is given by g(θ) = aθ + b. Applying the boundary conditions (6.38), we get a =, and b arbitrary. Thus λ = is a simple eigenvalue, and g(θ) = 1 is an eigenfunction. We need to solve ODE (6.37b) with µ =, which is given by The general solution of (6.39) is given by r 2 h (r ) + r h (r ) =. (6.39) h (r ) = Alog r + B (6.4) where A,B are real numbers. Since we are interested in a solution of Laplace equation in the unit disk, we are looking for bounded solutions, and thus A =. Thus we obtain the following solution of Laplace equation: v (r,θ) = h (r )g (θ) = 1 (6.41) (ii) (λ > ) When λ >, we may write λ = µ 2 where µ >. The ODE (6.37a) then becomes g (θ) + µ 2 g(θ) =, whose general solution is given by g(θ) = a cos(µθ) + b sin(µθ). Applying the boundary conditions (6.38), we get a = a cos(2µπ) + b sin(2µπ) bµ = aµsin(2µπ) + bµcos(2µπ). The above system of linear equations may be written as 1 cos(2µπ) sin(2µπ) a = sin(2µπ) 1 cos(2µπ) b. (6.42) IIT Bombay
9 6.4. Poisson s formula and its applications 155 Note that the system (6.42) possesses non-trivial solutions if and only if cos(2µπ) = 1. This means that µ. Thus we get a sequence of eigenvalues λ n = n 2 indexed by n, and a basis for the corresponding eigenspace is given by {cos(nθ), sin(nθ)}. For each n, we need to solve ODE (6.37b) with µ = n 2, which is given by r 2 h (r ) + r h (r ) n 2 h(r ) =. (6.43) Note that this is an ODE of Cauchy-Euler type, which reduces to a constant coefficient ODE by a change of variable s = log r, and hence can be easily solved, which is left as an exercise to the reader. The solutions are: For each n, the general solution of (6.43) is given by h n (r ) = Ar n + B r n (6.44) where A,B are real numbers. Since we are interested in a solution of Laplace equation in the unit disk, we are looking for solutions that remain bounded as r, and thus B =. For each n, we obtain the following solution of Laplace equation: v n (r,θ) = h n (r )g n (θ) = r n (a n cos(nθ) + b n sin(nθ)) (6.45) (iii) (λ < ) When λ <, we may write λ = µ 2 where µ >. The ODE (6.37a) then becomes g (θ) µ 2 g(θ) =, whose general solution is given by g(θ) = ae µθ + b e µθ. Applying the boundary conditions (6.38), we get a + b = ae 2µπ + b e 2µπ aµ bµ = aµe 2µπ bµe 2µπ. The above system of linear equations may be written as 1 e 2µπ 1 e 2µπ a 1 e 2µπ 1 + e 2µπ = b. (6.46) Note that the determinant of the coefficient matrix in (6.46) is zero if and only if e 2µπ + e 2µπ = 2. Note that this relation is of the form α + 1 = 2 for which α = 1 α is the only solutoin. Since µ >, the equation e 2µπ + e 2µπ = 2 has no positive solution for µ. Thus any λ with λ < is not an eigenvalue. Consider the function v defined by superposition of v and v n for n (which were obtained above) as follows: v(r,θ) = a 2 + r n (a n cos(nθ) + b n sin(nθ)), (6.47) n=1 where a n, b n are real numbers. We choose these constants so that v satisfies the boundary condition v(1,θ) = F (θ). Thus we are led to the relation v(1,θ) = F (θ) = a 2 + (a n cos(nθ) + b n sin(nθ)). (6.48) n=1 October 4, 215
10 156 Chapter 6. Laplace equation Let the constants a n, b n be chosen as the Fourier coefficients of F (θ), i.e., a n = 1 π 2π F (θ)cos(nθ) dθ, b n = 1 π 2π F (θ)sin(nθ) dθ. (6.49) This finishes the construction of a formal solution of the Dirichlet boundary value problem. The solution is given by (6.47), where the coefficients a n, b n are given by (6.49). Step 2: Formal solution is indeed a solution We are going to show that v defined by (6.47) is indeed a solution of the Dirichlet boundary value problem (6.32). Substituting the values of a n, b n from (6.49) into (6.47), we get v(r,θ) = 1 2π 2π F (τ) dτ + 1 π 2π r n n=1 F (τ)cos(n(τ θ)) dτ (6.5) As in the proof of Theorem 4.21, from the definition of a n, b n it follows that there exists an M > such that for all n we have a n M, n a n M; and b n M, n b n M, which follow from the formulae (6.49) and using integration by parts once in them. As a consequence of the above estimates on a n, b n, the Fourier series for F (θ) converges uniformly, and hence the formal series in (6.47) for all r R < 1 where < R < 1 is an arbitrary but fixed number. It also follows that series in (6.47) can be differentiated term-by-term and the resultant series also converges uniformly for all r R < 1. Since each term in the series is a solution to Laplace equation, so will be v(r,θ). This proves that the formal solution is indeed a solution. Step 3: Poisson s formula Since the series in (6.47) is uniformly convergent for r R < 1, the summation and integral in (6.5) can be interchanged. Thus we get v(r,θ) = 1 π = 1 π = 1 2π 2π 2π 2π 1 F (τ) 2 + r n cos(n(τ θ)) dτ n=1 1 e F (τ) 2 + r n i n(τ θ) + e in(τ θ) dτ n=1 2 F (τ) 1 + r e i(τ θ) r e i(τ θ) + dτ 1 r e i(τ θ) 1 r e i(τ θ) Simplifying the last equation, we get the following Poisson s summation formula for the solution of Laplace equation for r < 1: v(r,θ) = 1 r 2 2π 2π F (τ) dτ. (6.51) 1 2r cos(τ θ) + r 2 IIT Bombay
11 6.4. Poisson s formula and its applications 157 In cartesian coordinates, the formula (6.51) becomes u(x) = 1 x 2 f (y) d y. (6.52) 2π x y 2 S(,1) Checking that u(x) given by (6.52) satisfies the boundary condition (6.29b) is left as an exercise (see Exercise 6.13) Harnack s inequality Theorem 6.12 (Harnack s inequality). Let B(, R) d be the open ball with center at and having radius R. Let u : B(, R) [, ) be a harmonic function. Then for any x B(, R), the following inequalities hold: R d 2 (R x ) (R + x ) u() u(x) Rd 2 (R + x ) u(). (6.53) d 1 (R x ) d 1 Proof. Let us prove this for d = 2. Recall from Poisson s formula (6.52) u(x) = 1 x 2 u(y) d y. (6.54) 2π x y 2 S(,1) As a consequence of triangle inequality, for y S(,1) the following inequalities hold: 1 x x y 1 + x. Using these inequalities in the Poisson s formula (6.54) yields u(x) 1 + x 1 1 x 2π S(,1) u(y) d y = 1 + x 1 x u(), and u(x) 1 x x 2π S(,1) u(y) d y = 1 x 1 + x u(). This completes the proof of Harnack s inequalities for d = 2. Theorem 6.13 (Liouville s theorem). If u : d is a harmonic function that is bounded below, then u must be a constant function. Proof. Let m be such that u(x) m for all x d. Then the function v(x) := u(x) m is a non-negative harmonic function. Applying Harnack s inequality to the function v yields for every R >, and for every x B(, R) R d 2 (R x ) (R + x ) v() v(x) Rd 2 (R + x ) v(). (6.55) d 1 (R x ) d 1 October 4, 215
12 158 Chapter 6. Laplace equation Fix an x d, and R > such that R > x. Passing to the limit in the inequalities (6.55), we get v() v(x) v(). Thus we get v(x) = v(). Since x is arbitrary, we conclude that v is a constant function. As a consequence, u is a constant function. 6.5 Dirichlet problem on a rectangle In this section, Dirichlet problem for Laplace equation is considered on a rectangle. A solution is obtained using the method of separation of variables. u x x + u yy = for < x < π, < y < A, (6.56a) u(, y) = u(π, y) = for y A, (6.56b) u(x,a) = for x π, (6.56c) u(x,) = f (x) for x π. (6.56d) We use the method of separation of variables to solve the boundary value problem (6.56). Substituting u(x, y) = X (x)y (y) in the equation (6.56a), and re-arranging the terms, we get X (x) X (x) = Y (y) Y (y) (6.57) Since the LHS and RHS of the equation (6.57) are functions of the variables x and y respectively, each of them must be a constant function. Thus we have, for some constant λ, From the equation (6.58), we get the following two ODEs X (x) X (x) = Y (y) = λ (6.58) Y (y) X (x) + λx (x) =, Y (y) λy (y) =. (6.59a) (6.59b) Since we are interested in finding a non-trivial solution, the boundary conditions (6.56b)- (6.56c) give rise to the following conditions on the functions X and Y : X () =, X (π) =, (6.6a) Y (A) =. (6.6b) The ODE (6.59a) has non-trivial solutions satisfying the boundary conditions (6.6a) if and only if λ = n 2 for some n. Thus λ n = n 2 are eigenvalues and the corresponding eigenfunctions are given by X n (x) = sin nx. IIT Bombay
13 6.6. Weak maximum principle 159 For each n, the solution of ODE (6.59b) with λ = λ n = n 2 satisfying the condition (6.6b) (upto a constant multiple) is given by Y n (y) = sinh n(a y). Since for each n, the function u n (x, y) = X n (x)y n (y) is a solution (6.56a), we propose a formal solution of (6.56a) by superposition of the sequence of solutions (u n ) as follows: u(x, y) = b n sin nx sinh n(a y). (6.61) n=1 The coefficients b n in the formula (6.61) will be determined using the condition (6.56d). Thus we get f (x) = b n sin nx sinh na. (6.62) n=1 Choosing c n to be the Fourier sine coefficients of f, we get b n = c n. That is, if f is sinh na given by then c n are given by n=1 f (x) = c n = 2 π π c n sin nx, (6.63) n=1 f (s)sin ns d s. (6.64) Thus the formal solution of (6.56) is given by 2 π u(x, y) f (s)sin ns d s sin nx sinh n(a y). (6.65) π sinh na The formal solution (6.65) is indeed a solution to (6.56) under some conditions, and is proved in Theorem Weak maximum principle Let 2 be a bounded domain. Let u : 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 6.14 (Weak maximum principle). Let 2 be a bounded domain. Let u C 2 () 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, 2 u and 2 u. As a consequence, u at an interior maximum x 2 y 2 point. Thus if v is a function such that v > in, then the maximum value of v on October 4, 215
14 16 Chapter 6. Laplace equation 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 2: Define the function v ε by v ε (x, y) := u(x, y) + ε(x 2 + y 2 ). (6.66) Then v ε C 2 () C (). Note that v ε > in and thus v ε attains its maximum only on the boundary. Denoting M := max u, L := max(x 2 + y 2 ), (6.67) we have v ε (x, y) M + εl, (x, y). (6.68) Since u(x, y) v ε (x, y) for all (x, y), we have u(x, y) M + εl, (x, y). (6.69) Note that last inequality holds for every ε >. Thus taking limits as ε, we get u(x, y) M, (x, y). (6.7) This finishes the proof. Corollary 6.15 (Weak minimum principle). Let 2 be a bounded domain. Let u C 2 () C () be a harmonic function in. Then the minimum value of u in is achieved on the boundary. Proof. Define 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 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. IIT Bombay
15 6.6. Weak maximum principle 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 6.17 (Uniqueness of solutions to Dirichlet problem). bounded domain and consider the Dirichlet problem on : Let 2 be a u = f on, u = g on, where g is a continuous function on. Then the Dirichlet problem has at most one solution in the class C 2 () C (). Proof. Assume that u and v belonging to the class C 2 () C () solve the Dirichlet problem. Define w := u v. Then w belongs to the class C 2 () C () and solves the homogeneous Dirichlet 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 (i) Note that uniqueness of solutions to Dirichlet problem was already proved in Theorem 6.9. Recall that its proof required that n u is defined on the boundary, and thus the uniqueness of solutions to Dirichlet problem could be proved only for u C 1 (). Thanks to maximum principle, uniqueness result is valid for any harmonic function which belongs to the space C (). (ii) 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 Dirichlet problem on the upper half plane: u(x, y) = for x, < y <, u(x,) = 1 for all x. This problem has at least two solutions, namely u 1 (x, y) = xy, and, u 2 (x, y) =. Theorem Let be a bounded domain. For i = 1,2, let u i C 2 () C () solve the following Dirichlet problems on : u i = f on, u = g i on, October 4, 215
16 162 Chapter 6. Laplace equation where g i is a continuous function on for each i = 1,2. Then u 1 and u 2 satisfy max u 1 (x, y) u 2 (x, y) max g 1 (x, y) g 2 (x, y). (6.71) Proof. Define w := u 1 u 2. Then w solves the following Dirichlet problem w = on, w = g 1 g 2 on. Applying the weak maximum principle and the weak minimum principle, we get min(g 1 g 2 ) w(x, y) max(g 1 g 2 ) (x, y), (6.72) which finishes the proof of the theorem. Theorem 6.2 (Harnack). Let be a bounded domain in 2. Let (u n ) be a sequence of functions in C 2 () C () that are harmonic on, such that u n = f n on. Let f n f in C ( ) i.e., f n f uniformly on. Then (i) The sequence (u n ) converges uniformly on. (ii) Let u denote the uniform limit of the sequence (u n ). Then u is harmonic on, and u C ( ). (iii) u = f on. Proof. For each (m, n), the function u n u m is harmonic on, and belongs to C (). On the boundary, u n u m = f n f m Since ( f n ) converges uniformly on, it is a uniformly Cauchy sequence. Thus for a given ε >, there exists an N such that for all (x, y), and for all m N and n N, the following inequality holds. Applying Theorem 6.19, we get f m (x, y) f n (x, y) < ε. (6.73) max u m (x, y) u n (x, y) max f m (x, y) f n (x, y). (6.74) From the inequalities (6.73) and (6.74), we get for all m N and n N the following inequality max u m (x, y) u n (x, y) max f m (x, y) f n (x, y) < ε. (6.75) Thus the sequence (u n ) is a Cauchy sequence in C (), and hence converges uniformly, to u C (). Since each member of the sequence u n satisfies mean value property, u also satisfies mean value property being the uniform limit of the sequence (u n ). Since u is a continuous function and satisfies mean value property, it follows from Theorem 6.28 that u is harmonic in. It is clear that u = f on. IIT Bombay
17 6.6. Weak maximum principle 163 Theorem The formal solution (6.65) for the boundary value problem (6.56) is indeed a solution if the function f : [,π] is continuous, and satisfies f () = f (π) =, π ( f ) 2 (s) d s <. Proof. Recall the formal solution of (6.56) given by where u(x, y) = n=1 b n = 2 π sinh n(a y) b n sin nx sinh na, (6.76) π f (s)sin ns d s In order to establish that the formal solution given by (6.76) is indeed a solution to the BVP (6.56), we need to show that (i) The series in (6.76) is twice continuously differentiable w.r.t. both the variables x and y, thus making sure that u x x and u yy are meaningful. Then we need to show that equation (6.56a) is satisfied. (ii) We need to show that u(x, y) given by (6.76) is continuous upto the boundary of (,π) (,A) and that the boundary conditions (6.56b) - (6.56d) are satisfied. Step 1: Proof of (i) Note that sinh n(a y) sinh na = e n(a y) e n(a y) e na e na (6.77) ny 1 e 2n(A y) = e (6.78) 1 e 2nA e ny. (6.79) 1 e 2A In view of the following inequalities b sinh n(a y) n sin nx sinh na 1 1 e 2A e ny, and the fact that the series n=1 e ny converges uniformly for y > y, it follows that the series in (6.76) converges uniformly y > y. Since y > is arbitrary, it follows that u(x, y) is a continuous function in y >. Thus boundary values at x =, x = π, y = A are satisfied. By a similar reasoning, it follows that the series in (6.76) can be differentiated termby-term w.r.t. x and y twice and the resultant series define continuous functions. Step 2: Proof of (ii) In view of the assumptions on f, the fourier sine series of f converges to f uniformly. Denoting the k th partial sum of the series (6.76) by s k, we have s k (x, y) = k n=1 sinh n(a y) b n sin nx sinh na. (6.8) We wish to show that s k (x,) f uniformly on [,π], by showing that the sequence (s k ) is a Cauchy sequence in C [,π]. Note that for q > p, the function s p (x, y) s q (x, y) solves October 4, 215
18 164 Chapter 6. Laplace equation the IBVP (6.56) with s p (x,) s q (x,) = q n=p b n sin nx. Since the fourier sine series of f converges uniformly to f, the corresponding sequence of partial sums q n=1 b k sin nx is Cauchy in C [,π]. The required conclusion follows on applying maximum principle. Theorem Mean value property and its consequences Definition Let 2 be an open domain. A continuous function u : 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 2 ξ x 2 + η y 2 r 2 u(ξ,η) dξ dη. (6.81) Definition Let 2 be a domain. A continuous function u : 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 2πr where C r denotes the circle centered at P having radius r. C r u(s) d s (6.82) On expanding the line integral in the equation (6.82), we get the following form of the equation (6.82): u(x, y) = 1 2π 2π u(x + r cosθ, y + r sinθ) dθ. (6.83) Lemma Let u : be a continuous function. Then the following statements are equivalent. (i) u has the mean value property-i on. (ii) u has the mean value property-ii on. 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 2π 2π u(x + τ cosθ, y + τ sinθ) dθ. (6.84) Multiply the equation (6.84) with τ and then integrate w.r.t. τ over the interval [, r ]. This yields r τ u(x, y) dτ = 1 r 2π u(x + τ cosθ, y + τ sinθ) dθ τ dτ. (6.85) 2π IIT Bombay
19 6.7. Mean value property and its consequences 165 The last equation simplifies to Thus we get r r 2π u(x, y) = u(x + τ cosθ, y + τ sinθ)τ dθ dτ 2π = 1 u(ξ,η) dξ dη. 2π ξ x 2 + η y 2 r 2 u(x, y) = 1 πr 2 which proves that u has mean value property-i as well. ξ x 2 + η y 2 r 2 u(ξ,η) dξ dη, (6.86) Remark In view of Lemma 6.24, 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 property. This nomenclature is not standard. Theorem 6.26 (Mean value principle). Let 2 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 B(P, R) B[P, R]. We would like to prove that the function u has the surface mean value property, namely, for < r < R u(x, y ) = 1 2π 2π u(x + r cosθ, y + r sinθ) dθ. (6.87) Note that the left hand side of (6.87) is a real number, whereas the right hand side of (6.87) depends on r. Thus to prove (6.87), the strategy is to prove that the right hand side is a constant function of the variable r. To this end, define a function V : (, R) by V (r ) = 1 2π 2π u(x + r cosθ, y + r sinθ) dθ. (6.88) Let us compute dv for < r < R and prove that it is zero. As a consequence, we will d r then have V (r ) = lim V (ρ) = u(x, y ), (6.89) ρ which completes the proof of the mean value principle. It remain to show that dv = d r on (, R). dv d r (r ) = 1 2π 2π r u(x + r cosθ, y + r sinθ) dθ = ν u dσ S(,r ) = u d x =. (6.9) B(,r ) October 4, 215
20 166 Chapter 6. Laplace equation In proving the last equality, we have used the equation (6.24) and the fact that u is a harmonic function on. Theorem 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 2π 2π which holds as u satisfies mean value property on to obtain u(x + r cosθ, y + r sinθ) dθ (6.91) u x (x, y) = 1 2π u y (x, y) = 1 2π 2π 2π 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 2: Let us now prove that first order partial derivatives of u exist and are continuous. Since u has mean value property -I (6.81), we have u(x, y) = 1 u(ξ,η) dξ dη. (6.92) πr 2 ξ x 2 + η y 2 r 2 Writing the double integral on the right hand side as iterated integrals, the equation (6.92) takes the form u(x, y) = 1 1 u(ξ,η) dξ dη = πr 2 πr 2 ξ x 2 + η y 2 r 2 y+r y r x+ r 2 (η y) 2 dη x r 2 (η y) 2 u(ξ,η) dξ 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 y+r u(x + r 2 (η y) 2,η) u(x r 2 (η y) 2,η) dη. (6.93) πr 2 y r Since the right hand side of the equation (6.93) 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 If u : is continuous and u has the mean value property on, then u is harmonic in. IIT Bombay
21 6.7. Mean value property and its consequences 167 Proof. If u has the mean value property, then by Theorem 6.27 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 = (x, y ) 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 will be using the following identity in the computation that follows. d d r u(x + r cosθ, y + r sinθ) = x u(x + r cosθ, y + r sinθ).(cosθ,sinθ) Also note that the outward unit normal to the circle S(P, r ) is given by 1 r (x x, y ). Thus for < r < ε, we have < u(x, y) d x d y = ν u dσ B(P,r ) = S(P,r ) 2π r u(x + r cosθ, y + r sinθ) r dθ 2π = r u(x r + r cosθ, y + r sinθ) dθ = r 2πu(x, y r ) =. This contradiction means that our assumption that u at some point is non-zero is wrong. Thus u is a harmonic function. Theorem 6.29 (Derivative estimate). Let u : be a harmonic function, where 2 is an open set. Let k be a non-negative integer. There exists a C k > such that for each multi-index α with α = k, and for each ball B(x, r ), D α u(x ) (16k)k u(x, y) d x d y. (6.94) πr k+2 B(x,r ) Proof. Since the inequality (6.94) involves a non-negative integer k, we prove it by induction. For k =, the inequality (6.94) follows from mean value property of harmonic functions. We will prove the case k = 1 as the arguments used to prove that the required inequalities hold for k = l + 1, assuming that they hold for k = l are similar, and when k = 1 the proof can be presented cleanly without extra clutter of notations. Step 1: Case of k = 1: We prove the estimate for u x, as the proof is similar for u y. By Theorem 6.27, the partial derivative u x possesses mean value property. As a consequence, we have 1 4 u x (x ) = u π r 2 x (x, y) d x d y B(x, r 2 ) 1 4 = u.n π r 2 x dσ S(x, r 2 ) 4 r max S(x, r 2 ) u. (6.95) October 4, 215
22 168 Chapter 6. Laplace equation Let us estimate the quantity max S(x, r 2 ) u. Let x S(x, r ). Then 2 B x, r B(x 2, r ). Using (solid) mean value property, we get the estimate u(x) 4 u(x, y) d x d y. (6.96) πr 2 B(x,r ) Thus we get max u 4 u(x, y) d x d y. (6.97) S(x, r 2 ) πr 2 B(x,r ) The desired estimate (6.94) for k = 1 follows from (6.95) and (6.97). Step 2: Induction step: Let us assume that the result is true for k = l. Let us prove that the result is true for k = l + 1. Let α be a multi-index with α = l + 1. Consider a ball B(x, r ). Note that D α u will be a first order partial derivative of D β u for some multi-index β such that β = l. Proceeding as in Step 1 after writing down the (solid) mean value property for the function D α u on the ball B(x, Note that for x S(x, D α 2(l + 1) u(x ) r max r S(x, r ), we have l +1 l B x, l + 1 r B(x, r ). r l +1 ), we get l +1 ) D β u. (6.98) Applying induction hypothesis, we have following estimate for D β u(x) for x S(x, D β (16l )l u(x) πr l +2 B(x,r ) u(x, y) d x d y. (6.99) The inequalities (6.98) and (6.99) together imply the estimate (6.94) for k = l + 1. Recall that harmonic functions are infinitely differentiable (see Theorem 6.27). In fact harmonic functions are nicer than being infinitely differentiable, they have convergent Taylor series expansions at every point of, and such functions are known as real-analytic functions. Theorem 6.3 (Harmonic functions are analytic). Let u : be a harmonic function, where 2 is an open set. Then u is analytic in. Proof. Let x be fixed. Since u is an infinitely differentiable function, we may write down the Taylor series for u around the point x which is given by u(x) Dα u(x ) (x x α! ) α (6.1) α ( {}) 2 r l +1 ). IIT Bombay
23 6.7. Mean value property and its consequences 169 In order to show that u is analytic, we need to show that the series in (6.1) converges to u. This is shown by expressing u(x) = T N (x)+r N (x), where T N is the Taylor polynomial of degree N associated with u and R N is the remainder. We then show that R N as N. Note that T N (x is given by T N (x) = α ( {}) 2, α N D α u(x ) (x x α! ) α. (6.11) By mean value theorem, there exists z on the line joining x and x such that R N (x) has the following expression: R N (x) = u(x) T N (x) = α ( {}) 2, α =N+1 D α u(z) (x x α! ) α. (6.12) Let r = 1 4 d(x 1, ), and M = πr B(x u(x, y) d x d y. For each z B(x 2,2r ), r ), we have B(z, r ) B(x,2r ). Applying Theorem 6.29 with B(x, r ), and α = N + 1, we get D α (16(N + 1))N+1 u(z) πr N+1+2 (16(N + 1))N+1 πr N+1+2 (16(N + 1))N+1 M r N+1 B(z,r ) B(x,2r ) u(x, y) d x d y u(x, y) d x d y From the last string of inequalities, we conclude that 16 N+1 sup D α u(z) M (N + 1) N+1 (6.13) z B(x,r ) r We now state two general inequalities. They are (i) From Stirling s formula [3] which is concerned with growth rate of factorials, there exists a C > such that for all n the following inequality holds: n n C e n n! (6.14) (ii) This is an estimate on n!. By multinomial theorem, we have Hence we get d k = ( ) k = α =k α! α!. α! d α α! (6.15) In view of (6.14) and (6.15), from the estimate (6.13) we get 32e N+1 sup D α u(z) C M α! (6.16) z B(x,r ) r October 4, 215
24 17 Chapter 6. Laplace equation Let us now estimate the remainder term R N given by (6.12). R N (x) α ( {}) 2, α =N+1 α ( {}) 2, α =N+1 D α u(z) x x α! α 32e N+1 C M x x r N+1 (6.17) Let s > (to be chosen later) be such that x x s. Then the estimate (6.17) yields R N (x) α ( {}) 2, α =N+1 Choose an s such that 32e s = 1 r. In other words, choose s = r 4 the estimate (6.18) takes the form R N (x) α ( {}) 2, α =N+1 C M 32e s N+1 C M (6.18) r 128e. With this choice of s, 1 N+1 1 N+1 C M 2 N+2 C M (6.19) N Thus the remainder term R N as N. Thus u is an analytic function on. As a consequence of the last result, it follows that Cauchy problem for Laplace equation need not have solutions. This is the content of the next result. Theorem Let Γ 2 denote a segment of x-axis. Let g be a function defined on Γ. Then the Cauchy problem for u = where Cauchy data is given by u(x,) =, u (x,) = g(x) for (x,) Γ. (6.11) y has no solutions unless x g(x) is a real-analytic function. Proof. Let P(ξ,) Γ. Let B be an open disk with center at P such that B x-axis is contained in Γ. Let B + denote the part of B which lies above x-axis, i.e., B + = {(x, y) B : y }. Let u C 2 (B + ) be a solution of Laplace equation corresponding to the Cauchy data (6.11). Let u be extended to whole of B by reflection by setting for (x, y) B, u(x, y) := u(x, y) for y <. Verifying that the extended function u : B has the following properties (i) u C 2 (B), and (ii) u = on B is left as an exercise to the reader. By Theorem (6.3), the harmonic function u is analytic in B. In particular, u is analytic at P, and so is its derivative u. Thus g, which is the y IIT Bombay
25 6.7. Mean value property and its consequences 171 restriction of an analytic function to x-axis is itself real-analytic. This completes the proof of the theorem. Theorem 6.32 (Harnack s inequality). Let be a connected open subset of such that. There exists a constant C > such that holds for all non-negative harmonic functions u : [, ). sup u C inf u (6.111) Proof. We will first prove the result on each of the balls contained in the domain, and then deduce the result for using connectedness of and compactness of. Step 1: Let r < d(, ). Let x and y be such that x y < r. By (solid) mean value property, we have u(x) = 1 u(z) d z. (6.112) 4πr 2 B(x,2r ) Since B(y, r ) B(x,2r ), from the equation (6.112) we get u(x) = 1 u(z) d z 1 4πr 2 4πr 2 u(z) d z = 1 u(y). 4 (6.113) B(x,2r ) B(y,r ) Thus we have 1 4 u(y) u(x). Interchanging the roles of x and y we get 1 u(x) u(y). 4 Combining these two estimates, we get the following estimates, valid for all x and y satisfying x y < r : 1 u(y) u(x) 4u(y). (6.114) 4 Step 2: Since is connected and its closure is compact, there exists finitely many balls B 1,B 2,,B N having radius r, and B i B i+1 for i = 1,2,,N 1. Let x and y. Then there exists an i and a j such that x B j and y B k. Without loss of generality, we may assume that j k. For i = 1,2,,N 1, let z i B i B i+1 be fixed. u(x) 4u(z j ) 4 2 u(z j +1 ) 4 k j u(z k 1 ) 4 k+1 j u(y) 4 N u(y) (6.115) Using the first inequality in (6.114), and following a similar argument, we get Combining the inequalities (6.115) and (6.116), we get u(x) 1 u(y) (6.116) 4N 1 4 N u(y) u(x) 4N u(y), (6.117) which is valid for every x and y. The required conclusion (6.111) now follows. October 4, 215
26 172 Chapter 6. Laplace equation Strong maximum principle Let 2 be a bounded domain and u : 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 6.33 (Strong maximum principle). Let 2 be a domain (domain need not be bounded) and u : 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 2: Join P to P by a smooth curve γ 12. 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 radius d γ 2 with center at P denoted by B d γ (P ). Since u is a harmonic function in, it 2 has the mean value property in. Applying the mean value property 13 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 2 arguments of Step 2 with the disk B d γ (P 1 ) and then find a P 2 in a similar way. Continuing 2 this process till we get a k with the property that P B d γ (P k ) 14. We will then have 2 u(p) = M. This finishes the proof of the theorem. Remark This is a continuation of Remark 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 Let be the domain exterior to the unit disk. The function u(x, y) = log(x 2 +y 2 ) is a harmonic function on. It has neither a maximum value nor a minimum value in. 12 This is possible since is open and connected and hence it is path-connected. Prove this. 13 Mean Value Property-I or Mean Value Property-II? 14 Why is it possible to find such a k? IIT Bombay
27 6.8. Dirichlet principle Dirichlet principle Dirichlet principle characterises the solution of Dirichlet BVP in terms of minimizer of a functional. An analogous result in the context of systems of linear equations says that the solutoin of the linear system Ax = b (where A is a symmetric positive definite matrix) is the minimizer of the functional 1 2 x T Ax x T b. Theorem 6.36 (Dirichlet principle). Let 2 be a bounded domain with smooth boundary. Let f C (), g C ( ), and u C 2 (). The following statements are equivalent. (i) The function u is a solution to Dirichlet BVP u = f (x) in, (6.118a) u = g on. (6.118b) (ii) u is a minimizer of the functional J : S where S := {v C 2 () : v = g on }, and J (v) := 1 v 2 d x f (x)v(x) d x (6.119) 2 Proof. Proof of (i)= (ii): Let u be a solution to Dirichlet BVP (6.118). Let v S. Multiplying the equation (6.118a) with u v and then integrating on yields u(x)(u(x) v(x)) d x = f (x)(u(x) v(x)) d x Integrating by parts on the LHS of the last equation gives u(x).( u v)(x) d x = f (x)(u(x) v(x)) d x (6.12) Note that the boundary terms resulting from integration by parts do no appear as u = v on. Re-arranging terms in the equation (6.12), we get u(x) 2 d x f (x) u(x) d x = u(x). v(x) d x f (x) v(x) d x (6.121) Note that u. v d x 1 2 u v d x u 2 d x 1 2 u 2 d x + v 2 d x 1 2 v 2 d x (6.122) Using the inequality (6.122) in the equation (6.121), we get October 4, 215
28 174 Chapter 6. Laplace equation u(x) 2 d x f (x) u(x) d x 1 u 2 d x + v 2 d x f (x) v(x) d x 2 (6.123) Re-arranging the terms in (6.123), we get 1 u(x) 2 d x f (x) u(x) d x 1 v 2 d x f (x) v(x) d x 2 2 (6.124) Note that the equation (6.123) is nothing but J (u) J (v) Thus u is the minimizer of the functional J over the set S. Proof of (ii)= (i): Let v C (), and t. We have J (u + t v) = J (u) + t ( u. v f v) d x + t 2 v 2 d x 2 Rewriting the last equality, we get J (u + t v) J (u) = t ( u. v f v) d x + t 2 v 2 d x. (6.125) 2 Since the functional achieves its minimum at u, the real valued function h : defined by h(t) = J (u + t v) J (u) achieves its minimum at t =. Thus h () =, which yields u(x). v(x) d x = f (x) v(x) d x (6.126) Integrating by parts on the LHS of the equation (6.126) gives u(x) v(x) d x = f (x) v(x) d x (6.127) Since the equality (6.127) holds for every v C (), we conclude that for each x, u(x) = f (x). Note that u satisfies the boundary condition (6.118b) as u S. This proves that u is a solution to the BVP (6.118). IIT Bombay
29 Exercises 175 General Exercises 6.1. [16] Show that the laplace operator in two space dimensions takes the form u = 1 r r u u r r r 2 θ 2 in the polar coordinates x 1 = r cosθ, x 2 = r sinθ [16] Show that the laplace operator in three space dimensions takes the form u = 1 r r u u r r r 2 θ + 2 u 2 z 2 in the cylindrical coordinates x 1 = r cosθ, x 2 = r sinθ, x 3 = z [16] Show that the laplace operator in three space dimensions takes the form u = 2 u r + 2 u 2 r r u r 2 sin 2 ϕ θ u 2 r 2 ϕ + +cotϕ 2 r 2 u ϕ in the spherical coordinates x 1 = r cosθ sinϕ, x 2 = r sinθ sinϕ, x 3 = r cosϕ Let 2 be an open set and let u :. (a) Let (a, b) 2 and let T : 2 2 be given by T (x, y) = (x + a, y + b). Define be given by := {(x a, y b) : (x, y) }. Let T denote the restriction of T to. Show that u C 2 () if and only if u T C 2 ( ) and that (u T ) = ( u) T. In particular, u = (u T ) =. (b) Laplace equation is invariant under any real orthogonal transformation. Let U : 2 2 be any real orthogonal transformation i.e., U (x, y) = M U x y (x, y) 2, where M U is a real 2 2 invertible matrix whose inverse is equal to its transpose. Define := U 1 () and let Ũ be the restriction of U to. Show that u C 2 () if and only if u Ũ C 2 ( ) and that (u Ũ ) = ( u) Ũ. In particular, u = (u Ũ ) = Let u : 2 be a polynomial function, which is also radial and harmonic. Show that u is a constant function. October 4, 215
u 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 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 informationPartial 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
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 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 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 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 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 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 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 informationLAPLACE 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 informationUNIVERSITY 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
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 informationLEGENDRE POLYNOMIALS AND APPLICATIONS. We construct Legendre polynomials and apply them to solve Dirichlet problems in spherical coordinates.
LEGENDRE POLYNOMIALS AND APPLICATIONS We construct Legendre polynomials and apply them to solve Dirichlet problems in spherical coordinates.. Legendre equation: series solutions The Legendre equation is
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 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
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 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 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 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 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 informationHartogs Theorem: separate analyticity implies joint Paul Garrett garrett/
(February 9, 25) Hartogs Theorem: separate analyticity implies joint Paul Garrett garrett@math.umn.edu http://www.math.umn.edu/ garrett/ (The present proof of this old result roughly follows the proof
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 informationA 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
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 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 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 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 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 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 informationIntroduction 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
More informationNotes on Complex Analysis
Michael Papadimitrakis Notes on Complex Analysis Department of Mathematics University of Crete Contents The complex plane.. The complex plane...................................2 Argument and polar representation.........................
More informationMA 201: Method of Separation of Variables Finite Vibrating String Problem Lecture - 11 MA201(2016): PDE
MA 201: Method of Separation of Variables Finite Vibrating String Problem ecture - 11 IBVP for Vibrating string with no external forces We consider the problem in a computational domain (x,t) [0,] [0,
More information4.10 Dirichlet problem in the circle and the Poisson kernel
220 CHAPTER 4. FOURIER SERIES AND PDES 4.10 Dirichlet problem in the circle and the Poisson kernel Note: 2 lectures,, 9.7 in [EP], 10.8 in [BD] 4.10.1 Laplace in polar coordinates Perhaps a more natural
More informationLecture notes on Ordinary Differential Equations. S. Sivaji Ganesh Department of Mathematics Indian Institute of Technology Bombay
Lecture notes on Ordinary Differential Equations S. Ganesh Department of Mathematics Indian Institute of Technology Bombay May 20, 2016 ii IIT Bombay Contents I Ordinary Differential Equations 1 1 Initial
More informationMetric 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
More informationHarmonic Functions and Brownian motion
Harmonic Functions and Brownian motion Steven P. Lalley April 25, 211 1 Dynkin s Formula Denote by W t = (W 1 t, W 2 t,..., W d t ) a standard d dimensional Wiener process on (Ω, F, P ), and let F = (F
More informationIntroduction to Sturm-Liouville Theory and the Theory of Generalized Fourier Series
CHAPTER 5 Introduction to Sturm-Liouville Theory and the Theory of Generalized Fourier Series We start with some introductory examples. 5.. Cauchy s equation The homogeneous Euler-Cauchy equation (Leonhard
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 informationGeorgia Tech PHYS 6124 Mathematical Methods of Physics I
Georgia Tech PHYS 612 Mathematical Methods of Physics I Instructor: Predrag Cvitanović Fall semester 2012 Homework Set #5 due October 2, 2012 == show all your work for maximum credit, == put labels, title,
More informationNATIONAL 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;
More informationSolutions to Laplace s Equation in Cylindrical Coordinates and Numerical solutions. ρ + (1/ρ) 2 V
Solutions to Laplace s Equation in Cylindrical Coordinates and Numerical solutions Lecture 8 1 Introduction Solutions to Laplace s equation can be obtained using separation of variables in Cartesian and
More informationSAMPLE FINAL EXAM SOLUTIONS
LAST (family) NAME: FIRST (given) NAME: ID # : MATHEMATICS 3FF3 McMaster University Final Examination Day Class Duration of Examination: 3 hours Dr. J.-P. Gabardo THIS EXAMINATION PAPER INCLUDES 22 PAGES
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 informationTaylor and Laurent Series
Chapter 4 Taylor and Laurent Series 4.. Taylor Series 4... Taylor Series for Holomorphic Functions. In Real Analysis, the Taylor series of a given function f : R R is given by: f (x + f (x (x x + f (x
More information1 Directional Derivatives and Differentiability
Wednesday, January 18, 2012 1 Directional Derivatives and Differentiability Let E R N, let f : E R and let x 0 E. Given a direction v R N, let L be the line through x 0 in the direction v, that is, L :=
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 informationMath The Laplacian. 1 Green s Identities, Fundamental Solution
Math. 209 The Laplacian Green s Identities, Fundamental Solution Let be a bounded open set in R n, n 2, with smooth boundary. The fact that the boundary is smooth means that at each point x the external
More informationMath 2930 Worksheet Final Exam Review
Math 293 Worksheet Final Exam Review Week 14 November 3th, 217 Question 1. (* Solve the initial value problem y y = 2xe x, y( = 1 Question 2. (* Consider the differential equation: y = y y 3. (a Find the
More informationSPECTRAL PROPERTIES OF THE LAPLACIAN ON BOUNDED DOMAINS
SPECTRAL PROPERTIES OF THE LAPLACIAN ON BOUNDED DOMAINS TSOGTGEREL GANTUMUR Abstract. After establishing discrete spectra for a large class of elliptic operators, we present some fundamental spectral properties
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 informationApplied Mathematics Masters Examination Fall 2016, August 18, 1 4 pm.
Applied Mathematics Masters Examination Fall 16, August 18, 1 4 pm. Each of the fifteen numbered questions is worth points. All questions will be graded, but your score for the examination will be the
More informationPreliminary 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
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 informationPerron method for the Dirichlet problem.
Introduzione alle equazioni alle derivate parziali, Laurea Magistrale in Matematica Perron method for the Dirichlet problem. We approach the question of existence of solution u C (Ω) C(Ω) of the Dirichlet
More informationMultivariable Calculus
2 Multivariable Calculus 2.1 Limits and Continuity Problem 2.1.1 (Fa94) Let the function f : R n R n satisfy the following two conditions: (i) f (K ) is compact whenever K is a compact subset of R n. (ii)
More informationf (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
More informationHarmonic Functions and Brownian Motion in Several Dimensions
Harmonic Functions and Brownian Motion in Several Dimensions Steven P. Lalley October 11, 2016 1 d -Dimensional Brownian Motion Definition 1. A standard d dimensional Brownian motion is an R d valued continuous-time
More informationWaves on 2 and 3 dimensional domains
Chapter 14 Waves on 2 and 3 dimensional domains We now turn to the studying the initial boundary value problem for the wave equation in two and three dimensions. In this chapter we focus on the situation
More information2 A Model, Harmonic Map, Problem
ELLIPTIC SYSTEMS JOHN E. HUTCHINSON Department of Mathematics School of Mathematical Sciences, A.N.U. 1 Introduction Elliptic equations model the behaviour of scalar quantities u, such as temperature or
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 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 informationCauchy 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
More informationξ,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
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 informationCourse 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.........
More informationSOLUTION 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.
More informationAnalysis 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
More informationThe Helmholtz Equation
The Helmholtz Equation Seminar BEM on Wave Scattering Rene Rühr ETH Zürich October 28, 2010 Outline Steklov-Poincare Operator Helmholtz Equation: From the Wave equation to Radiation condition Uniqueness
More information2 Sequences, Continuity, and Limits
2 Sequences, Continuity, and Limits In this chapter, we introduce the fundamental notions of continuity and limit of a real-valued function of two variables. As in ACICARA, the definitions as well as proofs
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 informationu( 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
More informationLast Update: April 7, 201 0
M ath E S W inter Last Update: April 7, Introduction to Partial Differential Equations Disclaimer: his lecture note tries to provide an alternative approach to the material in Sections.. 5 in the textbook.
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 informationi=1 α i. Given an m-times continuously
1 Fundamentals 1.1 Classification and characteristics Let Ω R d, d N, d 2, be an open set and α = (α 1,, α d ) T N d 0, N 0 := N {0}, a multiindex with α := d i=1 α i. Given an m-times continuously differentiable
More informationElliptic PDEs of 2nd Order, Gilbarg and Trudinger
Elliptic PDEs of 2nd Order, Gilbarg and Trudinger Chapter 2 Laplace Equation Yung-Hsiang Huang 207.07.07. Mimic the proof for Theorem 3.. 2. Proof. I think we should assume u C 2 (Ω Γ). Let W be an open
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 informationSome lecture notes for Math 6050E: PDEs, Fall 2016
Some lecture notes for Math 65E: PDEs, Fall 216 Tianling Jin December 1, 216 1 Variational methods We discuss an example of the use of variational methods in obtaining existence of solutions. Theorem 1.1.
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 informationConsequences of Orthogonality
Consequences of Orthogonality Philippe B. Laval KSU Today Philippe B. Laval (KSU) Consequences of Orthogonality Today 1 / 23 Introduction The three kind of examples we did above involved Dirichlet, Neumann
More informationENGI 9420 Lecture Notes 8 - PDEs Page 8.01
ENGI 940 ecture Notes 8 - PDEs Page 8.0 8. Partial Differential Equations Partial differential equations (PDEs) are equations involving functions of more than one variable and their partial derivatives
More information= 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
More informationSobolev Spaces. Chapter Hölder spaces
Chapter 2 Sobolev Spaces Sobolev spaces turn out often to be the proper setting in which to apply ideas of functional analysis to get information concerning partial differential equations. Here, we collect
More informationDepartment of Mathematics
INDIAN INSTITUTE OF TECHNOLOGY, BOMBAY Department of Mathematics MA 04 - Complex Analysis & PDE s Solutions to Tutorial No.13 Q. 1 (T) Assuming that term-wise differentiation is permissible, show that
More informationSolutions to Laplace s Equations- II
Solutions to Laplace s Equations- II Lecture 15: Electromagnetic Theory Professor D. K. Ghosh, Physics Department, I.I.T., Bombay Laplace s Equation in Spherical Coordinates : In spherical coordinates
More informationTraces, extensions and co-normal derivatives for elliptic systems on Lipschitz domains
Traces, extensions and co-normal derivatives for elliptic systems on Lipschitz domains Sergey E. Mikhailov Brunel University West London, Department of Mathematics, Uxbridge, UB8 3PH, UK J. Math. Analysis
More informationMA3111S COMPLEX ANALYSIS I
MA3111S COMPLEX ANALYSIS I 1. The Algebra of Complex Numbers A complex number is an expression of the form a + ib, where a and b are real numbers. a is called the real part of a + ib and b the imaginary
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 informationBOUNDARY VALUE PROBLEMS ON A HALF SIERPINSKI GASKET
BOUNDARY VALUE PROBLEMS ON A HALF SIERPINSKI GASKET WEILIN LI AND ROBERT S. STRICHARTZ Abstract. We study boundary value problems for the Laplacian on a domain Ω consisting of the left half of the Sierpinski
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 informationu(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.
More informationPhysics 6303 Lecture 15 October 10, Reminder of general solution in 3-dimensional cylindrical coordinates. sinh. sin
Physics 6303 Lecture 15 October 10, 2018 LAST TIME: Spherical harmonics and Bessel functions Reminder of general solution in 3-dimensional cylindrical coordinates,, sin sinh cos cosh, sin sin cos cos,
More informationDiffusion 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
More informationMATH 205C: STATIONARY PHASE LEMMA
MATH 205C: STATIONARY PHASE LEMMA For ω, consider an integral of the form I(ω) = e iωf(x) u(x) dx, where u Cc (R n ) complex valued, with support in a compact set K, and f C (R n ) real valued. Thus, I(ω)
More informationSequences and Series of Functions
Chapter 13 Sequences and Series of Functions These notes are based on the notes A Teacher s Guide to Calculus by Dr. Louis Talman. The treatment of power series that we find in most of today s elementary
More informationHomework for Math , Fall 2016
Homework for Math 5440 1, Fall 2016 A. Treibergs, Instructor November 22, 2016 Our text is by Walter A. Strauss, Introduction to Partial Differential Equations 2nd ed., Wiley, 2007. Please read the relevant
More informationELEMENTARY APPLICATIONS OF FOURIER ANALYSIS
ELEMENTARY APPLICATIONS OF FOURIER ANALYSIS COURTOIS Abstract. This paper is intended as a brief introduction to one of the very first applications of Fourier analysis: the study of heat conduction. We
More informationExistence 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
More informationSubharmonic functions
Subharmonic functions N.A. 07 Upper semicontinuous functions. Let (X, d) be a metric space. A function f : X R { } is upper semicontinuous (u.s.c.) if lim inf f(y) f(x) x X. y x A function g is lower semicontinuous
More informationMTH4101 CALCULUS II REVISION NOTES. 1. COMPLEX NUMBERS (Thomas Appendix 7 + lecture notes) ax 2 + bx + c = 0. x = b ± b 2 4ac 2a. i = 1.
MTH4101 CALCULUS II REVISION NOTES 1. COMPLEX NUMBERS (Thomas Appendix 7 + lecture notes) 1.1 Introduction Types of numbers (natural, integers, rationals, reals) The need to solve quadratic equations:
More information