Class Meeting # 9: Poisson s Formula, Harnack s Inequality, and Liouville s Theorem

Similar documents
Class Meeting # 9: Poisson s Formula, Harnack s Inequality, and Liouville s Theorem

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

MATH COURSE NOTES - CLASS MEETING # Introduction to PDEs, Fall 2011 Professor: Jared Speck

Class Meeting # 2: The Diffusion (aka Heat) Equation

Midterm Exam, Thursday, October 27

Class Meeting # 13: Geometric Energy Estimates

6.3 Fundamental solutions in d

MATH 425, FINAL EXAM SOLUTIONS

Class Meeting # 12: Kirchhoff s Formula and Minkowskian Geometry

Class Meeting # 1: Introduction to PDEs

Harmonic Functions and Brownian Motion in Several Dimensions

Sobolev spaces. May 18

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

Lecture 5: Hodge theorem

Green s Functions and Distributions

(z 0 ) = lim. = lim. = f. Similarly along a vertical line, we fix x = x 0 and vary y. Setting z = x 0 + iy, we get. = lim. = i f

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:

ELLIPTIC PARTIAL DIFFERENTIAL EQUATIONS EXERCISES I (HARMONIC FUNCTIONS)

Harmonic Functions and Brownian motion

B ɛ (P ) B. n N } R. Q P

Course Description for Real Analysis, Math 156

VISCOSITY SOLUTIONS. We follow Han and Lin, Elliptic Partial Differential Equations, 5.

THE GREEN FUNCTION. Contents

Math The Laplacian. 1 Green s Identities, Fundamental Solution

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

Chapter 6: Rational Expr., Eq., and Functions Lecture notes Math 1010

u xx + u yy = 0. (5.1)

x = π m (a 0 + a 1 π + a 2 π ) where a i R, a 0 = 0, m Z.

Elliptic PDEs of 2nd Order, Gilbarg and Trudinger

Lecture Notes for LG s Diff. Analysis

Laplace s Equation. Chapter Mean Value Formulas

HARMONIC ANALYSIS. Date:

MATH 425, HOMEWORK 5, SOLUTIONS

THE DIRICHLET PROBLEM FOR THE LAPLACE OPERATOR

Partial regularity for fully nonlinear PDE

Math 342 Partial Differential Equations «Viktor Grigoryan

MATH 425, HOMEWORK 3 SOLUTIONS

Lecture Notes March 11, 2015

Mean value properties on Sierpinski type fractals

NOTES ON SCHAUDER ESTIMATES. r 2 x y 2

MATH 6337: Homework 8 Solutions

Fourier Transform & Sobolev Spaces

A RELATIONSHIP BETWEEN THE DIRICHLET AND REGULARITY PROBLEMS FOR ELLIPTIC EQUATIONS. Zhongwei Shen

Sobolev Spaces. Chapter 10

ALGEBRAIC STRUCTURE AND DEGREE REDUCTION

Lecture 4 (Sep. 18, 2017)

18.02 Multivariable Calculus Fall 2007

Boot camp - Problem set

Leplace s Equations. Analyzing the Analyticity of Analytic Analysis DEPARTMENT OF ELECTRICAL AND COMPUTER ENGINEERING. Engineering Math EECE

The Helmholtz Equation

17 The functional equation

WHY POLYNOMIALS? PART 1

MINIMAL GRAPHS PART I: EXISTENCE OF LIPSCHITZ WEAK SOLUTIONS TO THE DIRICHLET PROBLEM WITH C 2 BOUNDARY DATA

Handlebody Decomposition of a Manifold

Perron method for the Dirichlet problem.

CHAPTER 2. Laplace s equation

18 Green s function for the Poisson equation

SOLUTION OF POISSON S EQUATION. Contents

Distance in the Plane

ON THE REGULARITY OF SAMPLE PATHS OF SUB-ELLIPTIC DIFFUSIONS ON MANIFOLDS

A TALE OF TWO CONFORMALLY INVARIANT METRICS

arxiv: v1 [math.ap] 25 Jul 2012

THE HODGE DECOMPOSITION

ELECTROMAGNETIC SCATTERING FROM PERTURBED SURFACES. Katharine Ott Advisor: Irina Mitrea Department of Mathematics University of Virginia

= 2 x x 2 n It was introduced by P.-S.Laplace in a 1784 paper in connection with gravitational

is a weak solution with the a ij,b i,c2 C 1 ( )

Math 361: Homework 1 Solutions

Exercise 2. Prove that [ 1, 1] is the set of all the limit points of ( 1, 1] = {x R : 1 <

Asymptotic behavior of infinity harmonic functions near an isolated singularity

Now, the above series has all non-negative terms, and hence is an upper bound for any fixed term in the series. That is to say, for fixed n 0 N,

Regularity for Poisson Equation

CS 542G: The Poisson Problem, Finite Differences

MATH 51H Section 4. October 16, Recall what it means for a function between metric spaces to be continuous:

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

Local Asymmetry and the Inner Radius of Nodal Domains

A LOWER BOUND FOR THE GRADIENT OF -HARMONIC FUNCTIONS Edi Rosset. 1. Introduction. u xi u xj u xi x j

Probability inequalities 11

BOUNDARY VALUE PROBLEMS ON A HALF SIERPINSKI GASKET

Subharmonic functions

ESTIMATES FOR THE MONGE-AMPERE EQUATION

2.4 The Extreme Value Theorem and Some of its Consequences

APPLICATIONS OF DIFFERENTIABILITY IN R n.

THE INVERSE FUNCTION THEOREM

Final Exam May 4, 2016

An overview of key ideas

The Mollifier Theorem

4.1. Dirichlet, Poisson and Neumann boundary value problems

Power series and Taylor series

Partial Differential Equations

Geometry and the Dirichlet Problem in Any Co-dimension

Elementary Theory and Methods for Elliptic Partial Differential Equations. John Villavert

Algebraic Geometry Spring 2009

Cauchy Integral Formula Consequences

Math 40510, Algebraic Geometry

Metric spaces and metrizability

The shortest route between two truths in the real domain passes through the complex domain. J. Hadamard

Regularity estimates for the solution and the free boundary of the obstacle problem for the fractional Laplacian

HARMONIC FUNCTIONS TSOGTGEREL GANTUMUR

Recall that any inner product space V has an associated norm defined by

Exercises - Chapter 1 - Chapter 2 (Correction)

Transcription:

MATH 85 COUSE NOTES - CLASS MEETING # 9 85 Introduction to PDEs, Fall 0 Professor: Jared Speck Class Meeting # 9: Poisson s Formula, Harnack s Inequality, and Liouville s Theorem epresentation Formula for Solutions to Poisson s Equation We now derive our main representation formula for solution s to Poisson s equation on a domain Ω Theorem (epresentation formula for solutions to the boundary value Poisson equation) Let Ω be a domain with a smooth boundary, and assume that f C ( Ω) and g C( Ω ) Then the unique solution u C (Ω) C( Ω) to (0) u( x) = f( x ), x Ω n, can be represented as u( x) = g( x ), x Ω (0) u(x) = f(y)g(x, y Ω where G(x, y) is the Green function for Ω ) d n y + g( ˆN( G(x, dσ, Ω Poisson kernel Proof Applying the epresentation formula for u Proposition, we have that (0) u(x) = Φ(x y)f(y) d n y Φ(x ˆN( u( dσ + g( ˆN ( Φ( x dσ ecall also that Ω (04) G(x, y) = Φ(x y) φ(x, y), where (05) y φ(x, y) = 0, x Ω, and (06) G( x, = 0 when x Ω and σ Ω Ω The expression (0) is not very useful since don t know the value of ˆN ( u( along Ω To fix this, we will use Green s identity Applying Green s identity to the functions u(y) and φ(x, y), and recalling that y φ( x, y) = 0 for each fixed x Ω, we have that Ω

MATH 85 COUSE NOTES - CLASS MEETING # 9 u(y) Φ(x u( (07) 0 = φ(x, y) f(y) d n y φ(x, ˆNu( dσ + g( ˆNφ(x, dσ Ω Ω Ω Subtracting (07) from (0), and using (04), we deduce the formula (0) Poisson s Formula Let s compute the Green function G( x, y) and Poisson kernel P ( x, = ˆNG(x, from (0) in the case that Ω = B ( 0) is a ball of radius centered at the origin We ll use a technique called the method of images that works for special domains Warning 0 Brace yourself for a bunch of tedious computations that at the end of the day will lead to a very nice expression The basic idea is to hope that φ( x, y) from the decomposition G( x, y) = Φ( x y) φ( x, y ), where φ( x, y) is viewed as a function of x that depends on the parameter y, is equal to the Newtonian potential generated by some imaginary charge q placed at a point x B c ( 0 ) To ensure that G(x, = 0 when σ B ( 0 ), q and x have to be chosen so that along the boundary { y y = }, φ( x, y) = In a nutshell, we guess that (08) 4π x y G(x, y) = 4π x y + q 4 π x y, φ( x,y)? and we try to solve for q and x so that G( x, y) vanishes when y = q emark 0 Note that y 4π x y = 0, which is one of the conditions necessary for constructing G( x, y) By the inition of G(x, y), we must have G(x, y) = 0 when y =, which implies that (09) Simple algebra then leads to 4π x y = q 4π x y (00) x y = q x y When y =, we use (00) to compute that (0) x x y + = x y = q x y = q ( x x y + ), where denotes the Euclidean dot product Then performing simple algebra, it follows from (0) that (0) x + q ( + x ) = y ( x q x ) Now since the left-hand side of (0) does not depend on y, it must be the case that the right-hand side is always 0 This implies that x = q x, and also leads to the equation

MATH 85 COUSE NOTES - CLASS MEETING # 9 (0) q 4 x q ( + x ) + = 0 Solving (0) for q, we finally have that (04) q = x, (05) Therefore, x = x x (06) φ(x, y) = 4π x x x y (07) φ( 0, y) =,, where we took a limit as x 0 in (06) to derive (07) Next, using (08), we have (08) G(x, y) = 4π x y +, 4π x x x y x 0, (09) G( 0, y ) = + 4π y For future use, we also compute that (00) y G(x, y) = x y 4π x y Now when σ B ( 0 ), (00) and (04) imply that x y + 4π x x y (0) x σ = x σ x Therefore, using (00) and (0), we compute that ( ) = x σ x x σ x σ x x x σ (0) σg x, σ + = 4π x σ 4π x σ 4π x σ + 4π x σ σ x = ( ) 4π x σ Using (0) and the fact that Nˆ ( σ ) = σ, w e deduce x (0) ˆ ˆN ( G( x, = σg( x, N( = x σ

4 MATH 85 COUSE NOTES - CLASS MEETING # 9 emark 0 If the ball were centered at the point p instead of the origin, then the formula (0) would be replaced with ˆ x p (04) ˆN( G( x, = σg( x, N( σ ) = x σ Let s summarize this by stating a lemma Lemma 0 The Green function for a ball B ( ) is (05a) (05b) p ( ) = + G x, y, x p, 4π x y 4π x p (x p) (y p) G( p, y ) = 4π y p + Furthermore, if x B (p) and σ B (p), then ( ) = x p (05c) ˆN( G x, σ x σ We can now easily derive a representation formula for solutions to the Laplace equation on a ball Theorem (Poisson s formula) Let B (p) be a ball of radius centered at p = ( p, p, p ), and let x = ( x x ) ( ( )) ( ( )) (,, (x denote a point in Let g C B p Then the unique solution u C B p C B p)) of the PDE { u(x) = 0, x B ( p), (06) u( x) = g( x ), x B (p), can be represented using the Poisson formula: ( ) = x p (07) u x x p B (p) g( dσ x σ emark 0 In n dimensions, the formula (07) gets replaced with x p g( (08) u(x) = ω n B ( p) x σ where as usual, ωn is the surface area of the unit ball in Proof The identity (07) follows immediately from Theorem and Lemma 0 n n dσ, Harnack s inequality We will now use some of our tools to prove a famous inequality for Harmonic functions The theorem provides some estimates that place limitations on how slow/fast harmonic functions are allowed to grow

MATH 85 COUSE NOTES - CLASS MEETING # 9 5 Theorem (Harnack s inequality) Let B ( 0) n be the ball of radius centered at the origin, and let u C ( B( 0)) C( B ( 0)) be the unique solution to (06) Assume that u is non-negative on B ( 0 ) Then for any x B ( 0 ), we have that (09) n ( x ) + x u( 0) u(x) n ( ) u ( 0 ) ( + x ) n ( x ) n Proof We ll do the proof for n = The basic idea is to combine the Poisson representation formula with simple inequalities and the mean value property By Theorem, we have that (00) ( ) = x u x B (0) g( x σ By the triangle inequality, for σ B ( 0) (ie σ = ), we have that x x σ x + Applying the first inequality to (00), and using the non-negativity of g, we deduce that (0) u(x) + x x B (0) Now recall that by the mean value property, we have that dσ g( dσ (0) u( 0) = g σ dσ ( ) (0) B Thus, combining (0) and (0), we have that ( + x ) (0) u(x) u(0), ( x ) which implies one of the inequalities in (09) The other one can be proved similarly using the remaining triangle inequality We now prove a famous consequence of Harnack s inequality The statement is also often proved in introductory courses in complex analysis, and it plays a central role in some proofs of the fundamental theorem of algebra Corollary 0 (Liouville s theorem) Suppose that u C ( n ) is harmonic on n Assume that there exists a constant M such that u(x) M for all x n, or such that u( x) M for all x n Then u is a constant-valued function Proof We first consider the case that u( x) M Let v = u + M Observe that v 0 is harmonic and verifies the hypotheses of Theorem Thus, by (09), if x n and is sufficiently large, we have that n ( x ) n ( + x ) (04) ( + ) v(0) v(x) v(0 ) x n ( x ) n Allowing in (04), we conclude that v( x) = v( 0) Thus, v is a constant-valued function (and therefore u is too)

6 MATH 85 COUSE NOTES - CLASS MEETING # 9 To handle the case u(x) M, we simply consider the function w(x) = u(x) + M in place of v( x ), and we argue as above

MIT OpenCourseWare http://ocwmitedu 85 Introduction to Partial Differential Equations Fall 0 For information about citing these materials or our Terms of Use, visit: http://ocwmitedu/terms