Applications of the Maximum Principle

Similar documents
u xx + u yy = 0. (5.1)

Wave Equation With Homogeneous Boundary Conditions

Strauss PDEs 2e: Section Exercise 2 Page 1 of 6. Solve the completely inhomogeneous diffusion problem on the half-line

Math 4263 Homework Set 1

GREEN S IDENTITIES AND GREEN S FUNCTIONS

Diffusion on the half-line. The Dirichlet problem

Math 342 Partial Differential Equations «Viktor Grigoryan

In this lecture we shall learn how to solve the inhomogeneous heat equation. u t α 2 u xx = h(x, t)

Suggested Solution to Assignment 7

MATH 425, FINAL EXAM SOLUTIONS

MA 201, Mathematics III, July-November 2016, Partial Differential Equations: 1D wave equation (contd.) and 1D heat conduction equation

Jim Lambers MAT 169 Fall Semester Lecture 6 Notes. a n. n=1. S = lim s k = lim. n=1. n=1

Module 7: The Laplace Equation

Electrodynamics PHY712. Lecture 3 Electrostatic potentials and fields. Reference: Chap. 1 in J. D. Jackson s textbook.

In this chapter we study elliptical PDEs. That is, PDEs of the form. 2 u = lots,

Math 124B January 31, 2012

Homework for Math , Fall 2016

Lecture 24. Scott Pauls 5/21/07

MATH 124B Solution Key HW 05

6 Non-homogeneous Heat Problems

Tangent Planes, Linear Approximations and Differentiability

Functional Analysis HW 2

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

MATH 131P: PRACTICE FINAL SOLUTIONS DECEMBER 12, 2012

Sobolev Spaces. Chapter 10

UNIVERSITY of LIMERICK OLLSCOIL LUIMNIGH

Vectors in Function Spaces

HOMEWORK 5. Proof. This is the diffusion equation (1) with the function φ(x) = e x. By the solution formula (6), 1. e (x y)2.

There are five problems. Solve four of the five problems. Each problem is worth 25 points. A sheet of convenient formulae is provided.

Jim Lambers MAT 285 Spring Semester Practice Exam 2 Solution. y(t) = 5 2 e t 1 2 e 3t.

Math 220A - Fall 2002 Homework 5 Solutions

Mathematical Methods - Lecture 9

Boundary-value Problems in Rectangular Coordinates

MA 201: Method of Separation of Variables Finite Vibrating String Problem Lecture - 11 MA201(2016): PDE

UNIVERSITY OF MANITOBA

Laplace s Equation. Chapter Mean Value Formulas

PH.D. PRELIMINARY EXAMINATION MATHEMATICS

MATH 251 Final Examination December 19, 2012 FORM A. Name: Student Number: Section:

SECTION (See Exercise 1 for verification when both boundary conditions are Robin.) The formal solution of problem 6.53 is

M545 Homework 8. Mark Blumstein. December 9, [f (x)] 2 dx = 1. Find such a function if you can. If it cannot be found, explain why not.

Jim Lambers MAT 460 Fall Semester Lecture 2 Notes

CHAPTER 3. Analytic Functions. Dr. Pulak Sahoo

13 PDEs on spatially bounded domains: initial boundary value problems (IBVPs)

Lecture Notes on Metric Spaces

Chain Rule. MATH 311, Calculus III. J. Robert Buchanan. Spring Department of Mathematics

MA 441 Advanced Engineering Mathematics I Assignments - Spring 2014

M.Sc. in Meteorology. Numerical Weather Prediction

The Maximum and Minimum Principle

Spotlight on Laplace s Equation

Analysis III Solutions - Serie 12

INTEGRATION: THE FUNDAMENTAL THEOREM OF CALCULUS MR. VELAZQUEZ AP CALCULUS

Final: Solutions Math 118A, Fall 2013

Analysis III (BAUG) Assignment 3 Prof. Dr. Alessandro Sisto Due 13th October 2017

Green s Functions and Distributions

18.02 Multivariable Calculus Fall 2007

Math 342 Partial Differential Equations «Viktor Grigoryan

MA 201: Partial Differential Equations Lecture - 10

MATH 124B Solution Key HW 03

1. Differential Equations (ODE and PDE)

Math 220a - Fall 2002 Homework 6 Solutions

Math 212-Lecture 8. The chain rule with one independent variable

CLASSIFICATION AND PRINCIPLE OF SUPERPOSITION FOR SECOND ORDER LINEAR PDE

Lecture 2 - First order linear PDEs and PDEs from physics

Section 12.6: Non-homogeneous Problems

UNIVERSITY OF MANITOBA

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

MATH 124A Solution Key HW 05

G: Uniform Convergence of Fourier Series

HOMEWORK 4 1. P45. # 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:

2 Two-Point Boundary Value Problems

PH.D. PRELIMINARY EXAMINATION MATHEMATICS

MATH-UA 263 Partial Differential Equations Recitation Summary

Strauss PDEs 2e: Section Exercise 4 Page 1 of 6

the convolution of f and g) given by

Lucio Demeio Dipartimento di Ingegneria Industriale e delle Scienze Matematiche

Complex Variables. Chapter 2. Analytic Functions Section Harmonic Functions Proofs of Theorems. March 19, 2017

23 Elements of analytic ODE theory. Bessel s functions

INTRODUCTION TO PDEs

Summer Jump-Start Program for Analysis, 2012 Song-Ying Li

Brownian Motion and Stochastic Calculus

On semilinear elliptic equations with measure data

Math 240 Calculus III

17 Source Problems for Heat and Wave IB- VPs

Jim Lambers ENERGY 281 Spring Quarter Lecture 5 Notes

Numerical Solutions to Partial Differential Equations

A First Course of Partial Differential Equations in Physical Sciences and Engineering

[2] (a) Develop and describe the piecewise linear Galerkin finite element approximation of,

PDEs, Homework #3 Solutions. 1. Use Hölder s inequality to show that the solution of the heat equation

Hermite Interpolation

MATH COURSE NOTES - CLASS MEETING # Introduction to PDEs, Spring 2018 Professor: Jared Speck

MATH 2400: Calculus III, Fall 2013 FINAL EXAM

Partial Differential Equations

MATH 173: PRACTICE MIDTERM SOLUTIONS

1 From Fourier Series to Fourier transform

MATH 425, HOMEWORK 3 SOLUTIONS

Strauss PDEs 2e: Section Exercise 1 Page 1 of 6

Additional Homework Problems

A global solution curve for a class of free boundary value problems arising in plasma physics

DIFFERENTIAL EQUATIONS

Transcription:

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 very useful for proving uniqueness of solutions of various problems. Theorem 1 The irichlet problem in a bounded domain has at most one solution in C 2 () C(). Proof: Assume that u 1 and u 2 are solutions of the irichlet problem on. Then the difference v u 1 u 2 is also a solution of the irichlet problem, with v 0 on. By the weak maximum principle, the maximum and minimum of v are assumed on the boundary, so v 0 in, and therefore u 1 = u 2. Another application is stability of solutions. Theorem 2 Let u 1, u 2 C 2 () C() be solutions of u = F in a bounded domain with irichlet boundary conditions u 1 = g 1, u 2 = g 2 on. Then max u 1(x, y) u 2 (x, y) max g 1(x, y) g 2 (x, y). Proof: Let v u 1 u 2. Then v is a harmonic function with irichlet boundary data g g 1 g 2. By the weak maximum principle, from which the theorem follows. Green s Identities min g(x, y) v(x, y) max g(x, y), (x, y), The following fundamental theorem from vector calculus can be used to prove several theorems regarding uniqueness of solutions. Theorem 3 (ivergence Theorem) Let be a bounded domain with piecewise smooth boundary, and let ψ C 1 () C(). Then ψ(x, y) dx dy = ψ(x(s), y(s)) n ds. From this theorem, we obtain the three Green s identities: 1

From ψ = u, we obtain u dx dy = n u ds. From ψ = v u u v and the Product Rule for differentiation, we obtain v n u u n v ds = (v u u v) dx dy = v u + u v u v u v dx dy = v u u v dx dy. From ψ = v u we obtain v n u ds = v u + u v dx dy. We can now prove the following uniqueness results. Theorem 4 Let be a bounded domain with smooth boundary. 1. The irichlet problem u = F on, u = g on has at most one solution. 2. If α 0, the Robin problem u = F on, u + α n u = g on has at most one solution. 3. Any two solutions of the Neumann problem u = f on, n u = g on differ by a constant. Proof: For each of these three problems, let u 1, u 2 be two solutions and let v u 1 u 2. Then v is harmonic in, and satisfies a homogeneous boundary condition. 1. By the weak maximum principle, v = 0 on. 2. From the third Green s identity, and the fact that v = 0 in, we have v n v ds = v v dx dy. From the boundary condition, v = α n v on, which yields α( n v) 2 ds = v 2 dx dy. 2

Because the integral on the left side is nonpositive, and the integral on the right side is nonnegative, both integrals must be zero. Because neither integrand changes sign, it follows that both integrands must be identically zero. Therefore, α n v 0 on, but from the boundary condition, we then have v 0 on. By the weak maximum principle, v 0 in. 3. From the third Green s identity, and the fact that n v = 0 on and v = 0 in, we have v v dx dy = v 2 dx dy = 0. It follows that v 0 in, so v must be a constant. The Maximum Principle for the Heat Equation Consider the following IBVP for heat equation on a bounded three-dimensional domain, that is also bounded in time: u t = k u, (x, y, z), t > 0, (1) Q T = {(x, y, z, t) (x, y, z), 0 < t T } P Q T = { {0}} { [0, T ]} Let C H be the class of functions that are twice differentiable in Q T with respect to x, y, z, once differentiable with respect to t, and continuous in Q T. We then have the following preliminary result. Proposition 5 Let v C H satisfy v t k v < 0 in Q T. Then v has no local maximum in Q T, and v achieves its maximum on P Q T. Proof: Let v have a local maximum at q Q T. Then v t (q) = 0, and v(q) < 0, which contradicts the assumption v t k v < 0 in Q T. Because v is continuous on Q T, v has a maximum on Q T, which must be assumed on the boundary Q T. If this maximum is achieved at a point q { {T }}, we must have v t (q) 0, and because it is a local maximum in (x, y, z), v(q) < 0. By continuity, this again contradicts the assumption v t k v < 0 in Q T. The preceding proposition leads us to the following maximum principle for the heat equation. Theorem 6 Let u C H be a solution of the heat equation (1) in Q T. Then u achieves its maximum and minimum values on P Q T. Proof: Let v = u ɛt. Then v t k v = u t k u ɛ, and therefore v t k v < 0. By the preceding proposition, v has no local maximum in Q T, and assumes its maximum on P Q T. Furthermore, max v = max v M = max u, Q T P Q T P Q T because v u in Q T. By rearranging v = u ɛt M, we obtain u M + ɛt on Q T. Letting ɛ 0, we obtain u M on Q T, and therefore the maximum of u is assumed on P Q T, since M is defined to be the maximum of u on P Q T. The preceding theorem has this consequence. 3

Theorem 7 Let u 1, u 2 be solutions of the heat equation u t k u = F (x, y, z, t), (x, y, z), 0 < t < T, with initial conditions u 1 (x, y, z, 0) = f 1 (x, y, z), u 2 (x, y, z, 0) = f 2 (x, y, z), (x, y, z), and boundary conditions u 1 (x, y, z, t) = h 1 (x, y, z, t), u 2 (x, y, z, t) = h 2 (x, y, z, t), (x, y, z), 0 < t < T. Then u 1 u 2 max f 2 f 2 + max {t>0} h 1 h 2. Proof: Let v u 1 u 2. Then v solves the IBVP v t k u = 0, (x, y, z), 0 < t < T, v(x, y, z, 0) = f 1 (x, y, z) f 2 (x, y, z), (x, y, z), v(x, y, z, t) = h 1 (x, y, z, t) h 2 (x, y, z, t), (x, y, z), 0 < t < T. It follows from the preceding theorem that v assumes its maximum and minimum on P Q T. This theorem states that the solution of the heat equation is stable; by setting f 1 = f 2 and h 1 = h 2, we obtain uniqueness as well. Specializing to the 1- heat equation, we obtain the following result. Corollary 8 Let be the solution of the IBVP u(x, t) = n=1 L e k(nπ/l)2 t (2) If the series u t ku xx = 0, 0 < x < L, t > 0, u(0, t) = u(l, t) = 0, t > 0, u(x, 0) = f(x), 0 < x < L. f(x) = n=1 L converges uniformly on (0, L), then the series (2) converges uniformly on (0, L) (0, T ), and u is a classical solution. Proof: From the uniform convergence of the series for f(x), we have that for any ɛ > 0 there exists N ɛ such that for k, l N ɛ, l B n sin nπx < ɛ, 0 < x < L. L n=k 4

Because l n=k L e k(nπ/l)2 t is a classical solution of the heat equation, the maximum principle applies, and therefore this solution is also bounded by ɛ on (0, L) (0, T ). It follows that the series (2) converges uniformly on (0, L) (0, T ). Because u satisfies the initial condition, boundary conditions, and PE, it is a classical solution. We conclude with a maximum principle for the Neumann problem. Proposition 9 Let u(x, t) be the solution of the IBVP u t = ku xx, a < x < b, t > 0, u(x, 0) = f(x), a < x < b, u x (a, t) = 0, u x (b, t) = 0, 0 < t < T. Let f C 1 ([a, b]) satisfy f (a) = f (b) = 0, and let Q T = (a, b) (0, T ]. Then max u x max f (x). Q T (a,b) Proof: Let w = u x. Then w satisfies the heat equation with irichlet boundary conditions, with initial condition w(x, 0) = f (x). The proposition then follows from the maximum principle for the heat equation. Exercises Section 7.9: Exercises 1, 4, 5, 8 5