Heat Equation: Maximum Principles

Similar documents
u t u 0 ( 7) Intuitively, the maximum principles can be explained by the following observation. Recall

We will look for series solutions to (1) around (at most) regular singular points, which without

EECE 301 Signals & Systems Prof. Mark Fowler

STRONG DEVIATION THEOREMS FOR THE SEQUENCE OF CONTINUOUS RANDOM VARIABLES AND THE APPROACH OF LAPLACE TRANSFORM

State space systems analysis

Explicit scheme. Fully implicit scheme Notes. Fully implicit scheme Notes. Fully implicit scheme Notes. Notes

MASSACHUSETTS INSTITUTE OF TECHNOLOGY 6.265/15.070J Fall 2013 Lecture 16 11/04/2013. Ito integral. Properties

PRELIM PROBLEM SOLUTIONS

Math 213b (Spring 2005) Yum-Tong Siu 1. Explicit Formula for Logarithmic Derivative of Riemann Zeta Function

SOLUTION: The 95% confidence interval for the population mean µ is x ± t 0.025; 49

ELEC 372 LECTURE NOTES, WEEK 4 Dr. Amir G. Aghdam Concordia University

Comments on Discussion Sheet 18 and Worksheet 18 ( ) An Introduction to Hypothesis Testing

Société de Calcul Mathématique, S. A. Algorithmes et Optimisation

10-716: Advanced Machine Learning Spring Lecture 13: March 5

Seunghee Ye Ma 8: Week 5 Oct 28

STA 4032 Final Exam Formula Sheet

Math 220B Final Exam Solutions March 18, 2002

LECTURE 13 SIMULTANEOUS EQUATIONS

Last time: Completed solution to the optimum linear filter in real-time operation

f(x) dx as we do. 2x dx x also diverges. Solution: We compute 2x dx lim

Applied Mathematical Sciences, Vol. 9, 2015, no. 3, HIKARI Ltd,

1. (a) If u (I : R J), there exists c 0 in R such that for all q 0, cu q (I : R J) [q] I [q] : R J [q]. Hence, if j J, for all q 0, j q (cu q ) =

a 1 = 1 a a a a n n s f() s = Σ log a 1 + a a n log n sup log a n+1 + a n+2 + a n+3 log n sup () s = an /n s s = + t i

REVIEW OF SIMPLE LINEAR REGRESSION SIMPLE LINEAR REGRESSION

(A sequence also can be thought of as the list of function values attained for a function f :ℵ X, where f (n) = x n for n 1.) x 1 x N +k x N +4 x 3

Fig. 1: Streamline coordinates

MASSACHUSETTS INSTITUTE OF TECHNOLOGY 6.436J/15.085J Fall 2008 Lecture 19 11/17/2008 LAWS OF LARGE NUMBERS II THE STRONG LAW OF LARGE NUMBERS

Sequences and Series of Functions

Math 5C Discussion Problems 2 Selected Solutions

Lecture Notes for Analysis Class

b i u x i U a i j u x i u x j

The Performance of Feedback Control Systems

Convergence of random variables. (telegram style notes) P.J.C. Spreij

Chapter 9. Key Ideas Hypothesis Test (Two Populations)

5. Matrix exponentials and Von Neumann s theorem The matrix exponential. For an n n matrix X we define

Chapter 4. Fourier Series

A Tail Bound For Sums Of Independent Random Variables And Application To The Pareto Distribution

Math 113 Exam 3 Practice

Chapter 6 Infinite Series

M227 Chapter 9 Section 1 Testing Two Parameters: Means, Variances, Proportions

Appendix: The Laplace Transform

Sequences A sequence of numbers is a function whose domain is the positive integers. We can see that the sequence

j=1 dz Res(f, z j ) = 1 d k 1 dz k 1 (z c)k f(z) Res(f, c) = lim z c (k 1)! Res g, c = f(c) g (c)

Isolated Word Recogniser

1.3 Convergence Theorems of Fourier Series. k k k k. N N k 1. With this in mind, we state (without proof) the convergence of Fourier series.

PAPER : IIT-JAM 2010

This section is optional.

A sequence of numbers is a function whose domain is the positive integers. We can see that the sequence

Ma 530 Infinite Series I

On the Signed Domination Number of the Cartesian Product of Two Directed Cycles

Notes 19 : Martingale CLT

REGRESSION WITH QUADRATIC LOSS

Infinite Sequences and Series

Application to Random Graphs

Math 113, Calculus II Winter 2007 Final Exam Solutions

MAT1026 Calculus II Basic Convergence Tests for Series

Integrable Functions. { f n } is called a determining sequence for f. If f is integrable with respect to, then f d does exist as a finite real number

20. CONFIDENCE INTERVALS FOR THE MEAN, UNKNOWN VARIANCE

EULER-MACLAURIN SUM FORMULA AND ITS GENERALIZATIONS AND APPLICATIONS

Lecture 19: Convergence

We are mainly going to be concerned with power series in x, such as. (x)} converges - that is, lims N n

1+x 1 + α+x. x = 2(α x2 ) 1+x

Quiz. Use either the RATIO or ROOT TEST to determine whether the series is convergent or not.

MIDTERM 3 CALCULUS 2. Monday, December 3, :15 PM to 6:45 PM. Name PRACTICE EXAM SOLUTIONS

Stat 3411 Spring 2011 Assignment 6 Answers

6.3 Testing Series With Positive Terms

Math Solutions to homework 6

1 Approximating Integrals using Taylor Polynomials

If a subset E of R contains no open interval, is it of zero measure? For instance, is the set of irrationals in [0, 1] is of measure zero?

Solutions to HW Assignment 1

Positive solutions of singular (k,n-k) conjugate boundary value problem

Singular Continuous Measures by Michael Pejic 5/14/10

Lesson 10: Limits and Continuity

Math 61CM - Solutions to homework 3

Math 113 Exam 3 Practice

On Elementary Methods to Evaluate Values of the Riemann Zeta Function and another Closely Related Infinite Series at Natural Numbers

Assignment 1 - Solutions. ECSE 420 Parallel Computing Fall November 2, 2014

Probability 2 - Notes 10. Lemma. If X is a random variable and g(x) 0 for all x in the support of f X, then P(g(X) 1) E[g(X)].

Fall 2013 MTH431/531 Real analysis Section Notes

f(x)g(x) dx is an inner product on D.

Chapter 3. Strong convergence. 3.1 Definition of almost sure convergence

Recitation 4: Lagrange Multipliers and Integration

Generalized Likelihood Functions and Random Measures

MA1200 Exercise for Chapter 7 Techniques of Differentiation Solutions. First Principle 1. a) To simplify the calculation, note. Then. lim h.

7.1 Convergence of sequences of random variables

lim za n n = z lim a n n.

2 Banach spaces and Hilbert spaces

Taylor expansion: Show that the TE of f(x)= sin(x) around. sin(x) = x - + 3! 5! L 7 & 8: MHD/ZAH

Expectation of the Ratio of a Sum of Squares to the Square of the Sum : Exact and Asymptotic results

Numerical Method for Blasius Equation on an infinite Interval

Regression with quadratic loss

MATH4822E FOURIER ANALYSIS AND ITS APPLICATIONS

n outcome is (+1,+1, 1,..., 1). Let the r.v. X denote our position (relative to our starting point 0) after n moves. Thus X = X 1 + X 2 + +X n,

Sequences, Series, and All That

Carleton College, Winter 2017 Math 121, Practice Final Prof. Jones. Note: the exam will have a section of true-false questions, like the one below.

1 Review and Overview

Series III. Chapter Alternating Series

5.1. The Rayleigh s quotient. Definition 49. Let A = A be a self-adjoint matrix. quotient is the function. R(x) = x,ax, for x = 0.

Solutions to home assignments (sketches)

On The Computation Of Weighted Shapley Values For Cooperative TU Games

Transcription:

Heat Equatio: Maximum Priciple Nov. 9, 0 I thi lecture we will dicu the maximum priciple ad uiquee of olutio for the heat equatio.. Maximum priciple. The heat equatio alo ejoy maximum priciple a the Laplace equatio, but the detail are lightly differet. Recall that the domai uder coideratio i Ω T Ω [0, T); Ω T (Ω 0}) ( Ω [0, T]). () We have the followig trog maximum priciple. Theorem. (Maximum priciple of the heat equatio) Aume u C (Ω T ) C(Ω T) olve i Ω T. i. (Weak maximum priciple) The u t u=0 () maxu=maxu. (3) Ω T Ω T ii. (Strog maximum priciple) Furthermore, if Ω i coected ad there exit a poit (x 0,t 0 ) Ω T uch that the u i cotat i Ω t 0. u(x 0, t 0 ) = maxu, (4) Ω T Remark. Ituitively, the maximum priciple ca be explaied by the followig obervatio. Recall that u(x, t) = Φ(x y, t ) u(x, )dy (5) R for ay < t, ad the oegative fuctio Φ atifie R Φ(x y, t ) dy =, (6) it i clear that u(x, t) at time t i obtaied from averagig u(y, ) at a earlier time, ad therefore the maximum hould be decreaig. We preet two proof of thi theorem. The firt i through clever dicuio of ig, the ecod i through etablihig a mea value property. Note that the firt method oly yield the weak maximum priciple, that i the maximum iide i bouded by that o the boudary, itead of the trog maximum priciple, that i the maximum ca oly be attaied at the boudary, ule the fuctio i a cotat. Proof. Theorem 3. (Weak maximum priciple) Let Ω R be ope ad bouded. Let u C (Ω T ) C 0 (Ω T), ad atify i Ω T, we the have (If T <, the upreme ca be replaced by maximum). u t u 0 (7) upu= up u. (8) Ω T Ω T. Note that Ω t0 caot be replaced by Ω T, a ca be ee by the followig coideratio. Let u olve the heat equatio with iitial value M (cotat) ad (ide) boudary value g(x) M. Now if u t 0 a t ց0, we ca exted u to t<0 by ettig u M.

Heat Equatio: Maximum Priciple Remark 4. The coectio to the theory of Laplace equatio ca be ee by writig the iequality ito u u t. (9) However it i ot clear how to really ue thi fact a u t 0 oly hold at (x 0, t 0 ). Proof. If uffice to prove the cocluio for Ω T for all T < T. Fix oe uch T. By compacte of Ω T we kow that the maximum i attaied at ome poit (x 0, t 0 ).. Firt coider the cae Now if (x 0, t 0 ) Ω T, at (x 0, t 0 ), we have u > u t i Ω T. (0) a) u t (t 0 ) 0 becaue otherwie u(x 0, t 0 ε) >u(x 0, t 0 ) for ε mall eough; b) u(x 0, t 0 ) 0 becaue x 0 maximize u o Ω t 0 }. But the two cocluio cotradict each other.. For the geeral cae we coider a auxiliary fuctio We eaily check ad ca apply. Fially lettig ε ց 0 we obtai the cocluio. u u t, () v(x, t) =u(x, t) ε t. () v v t (3) From weak maximum priciple oe immediately obtai the uiquee for ohomogeeou iitial/boudary-value heat equatio whe the domai Ω i bouded. The ituatio i more complicated whe Ω i ot bouded. Theorem 5. (Weak maximum priciple for Ω = R ) Suppoe the u t u 0 i Ω T ; u(x, t) Me λ x i Ω T for M,λ>0; u(x, 0)= g(x), (4) upu upg. (5) Ω T R Remark 6. The reao why we eed a boud o the growth rate at ifiity i the followig. Recall that the effect of the heat operator i to obtai olutio at later time from averagig olutio at earlier time via a kerel of the form e x /4. Note that the kerel decay at ifiity a e c x. Therefore, if the olutio grow fater tha it ivere rate e λ x, the it i poible to iput extra tuff from ifiity ad thu violate uiquee. Here i a example of how oe ca tore eergy at ± ad the iput it through averagig by the heat kerel, ad thu obtai a ozero olutio with zero iitial data. Set u(x, t) g () (t) ()! x (6) 0 with g(t) e t k t > 0, k >. (7) 0 t = 0. Oe ca verify that u atifie the equatio by formally differetiatig iide the ummatio.. The reao we ue T < T here i to coclude u t 0 at the maximizer. Sice the equatio i olved o Ω T oly, u t at t=t i ot kow.

Nov. 9, 0 3 With ome calculatio 3, oe ca further how that the erie o the RHS i majorized by the Taylor expaio of [ ( x exp U(x, t) = t θ )] t k t > 0. (8) 0 t 0 for ome θ > 0. Note that U(x, t) caot be bouded by M e λx for ay cotat M,λ. Proof. Firt we divide (0, T) ito ubiterval with ize r <. It uffice to prove the claim o a ubiterval. From ow o we aume T < 4 λ 4 λ Coider the auxiliary fuctio ( x v(x, t) u(x, t) δ (4 π (T +ε t)) e 4(T +ε t) / ). (9) where ε i choe uch that T + ε < 4 λ. It i eay to verify that The trategy i the how firt v t v 0. (0) v(x, t) upg () R for ay (x, t) ad the lettig δ ց 0. We firt otice that, o the phere x y = R, we have v(x, t) Me λ( x +R) δ (4 π (T + ε t)) / e ( R 4(T +ε t) ) up R g () oce R i big eough. Now apply the weak maximum priciple o the domai B R (0, T), we obtai v(x, t) up g. Proof Mea value property. The mea value property for the heat equatio tur out to be much more complicated tha the oe for the Laplace equatio. We firt defie the heat ball : N Defiitio 7. (Heat ball) For fixed x R, t R, r > 0, we defie E(x, t; r) (y, ) R + t, Φ(x y, t ) r }. (3) Remark 8. We try to gai ome idea of what E(x, t; r) look like. It i defied by Thu i the time directio, we have (4π (t )) / e x y t t r 4 π where the lower boud i from the fact that e x y 4(t ). Furthermore 4(t ) r. (4) (5) E(x, t; r) = 0}=E(x, t; r) =t}=x}. (6) Next, to fid out the correct formula, we eed to fid a kerel K(x y, t ) uch that K(x y, t )dy d (7) E(x,t;r) 3. See pp. 3 i F. Joh Partial Differetial Equatio, 4ed. The cotat θ i determied i Problem 3 o p. 73 of the ame book.

4 Heat Equatio: Maximum Priciple i idepedet of r. Notice that where Thu we have Thi implie whe (4 π (t )) E(x,t;r) / e x y 4(t ) r Oe ca omehow 4 verify that 5 4 E(0,0;) Thu fially we have the correct formulatio. (4 π (t )) / e x y 4(t ) (8) t = t/r, =/r, x = x/r, y = y/r. (9) K(x y, t )dy d = K(x y, t )r + dy d. (30) E(x,t ;) K(x, t )=K(x, t) r + (3) x =x/r, t = t/r. (3) x dxdt =. (36) t Theorem 9. (Mea value property for the heat equatio) Let u C (Ω T ) olve the heat equatio, the for each E(x, t; r) Ω T. u(x, t)= x y 4 r E(x,t;r) u(y, ) dy d. (37) (t ) Proof. Without lo of geerality we ca et (x, t)=(0, 0) ad deote E(0, 0; r) by. Defie φ(r) r u(y, ) y dy d= u(ry, r ) y E dy d. (38) I the followig we will omit the prime ad imply ue y ad y. thu we compute φ (r) = d [ ] u(ry, r ) y dy d dr E y = ( yi u) y i + r ( u) y dy d i= = y r + ( yi u) y i + r + ( u) y i= A+B. (39) 4. I have o idea how to come up with thi particular fuctio. 5. A the computatio i tediou ad log, we put it i a a footote ad furthermore oly preet the mai tep. Firt otice x E(0,0;) e 4t / }= 0 t (4 π t) 4 π ; x ( t) log }. (33) 4 π t Thu the itegral become [ /4π t 0 x tlog 4πt x dx ] dt= π/ + (+)/ (+)/ Γ ( 4π + ) 0 ( ) t log 4 π t + dt. (34) where we have ued polar coordiate ad the formula α() = π/ Γ ( for the volume of -dimeio ball. Now ettig + ) =4πt ad uig the formula ( λ z Γ(z)= t λ log ) z dt; Γ(z + ) = z Γ(z), (35) 0 t we will ee after careful calculatio that mot term cacel out ad what remai i 4.

Nov. 9, 0 5 Now we itroduce ad ote that ψ = 0 o. Now we have ψ B = r + = r + log ( 4π)+ y 4 ( u) i= + log r, (40) y i ( yi ψ)dy d 4 u ψ +4 u yi y i ψ dy d. (4) i= Now itegrate by part w.r.t., we have B = r + 4u ψ + 4 u yi y i ψ dy d i= = r + 4u ψ + 4 ( u yi y i ) y 4 dy d = r + uyi y i dy d A. (4) Coequetly, uig the equatio, we have 4u ψ φ (r) = A+B = r + 4 uψ uyi y i dy d = r + 4 u yi ψ yi u y i y i dy d = 0. (43) Thu we have ( ) y φ(r)= limφ(t)=u(0, 0) lim t 0 t 0 t dy d = 4 u(0, 0). (44) E t From the mea value property we immediately have the trog maximum priciple: Proof. (of the trog maximum priciple) Suppoe there i (x 0, t 0 ) Ω T uch that u(x 0, t 0 ) = max Ω Tu, the by pickig r mall eough o that E(x 0, t 0 ; r) Ω T, ad uig the mea value property, we coclude that u i cotat iide E(x 0, t 0 ; r). Next for ay (y 0, 0 ) Ω T uch that the lie egmet coectig x 0, y 0 i i Ω T, we ca how that u(y 0, 0 )=u(x 0,t 0 ) wheever 0 <t 0 by coverig the lie egmet coectig (y 0, 0 ) ad (x 0,t 0 ) with the heat ball. Fially, ice Ω i coected, ay y 0 ca be coected from x 0 via fiitely may lie egmet. Ad therefore u(y, )=u(x 0, t 0 ) for all y Ω, < t 0. Remark 0. It tur out there i alo a verio of mea value property ivolvig oly the urface itegral over E(x, t; r). The idea i to compute the itegral of [Φ(x y, T t) r ] (u t u) over E(x, t; r) via itegratio by part. Thi will give a repreetatio formula for u(x, t), which tur out to be u(x, t)= x y r E(x,t;r) u(y, ) ds. (45) (t ) See J. Jot Partial Differetial Equatio, pp. 8 83 for detail. The trog maximum priciple ca alo be proved uig thi formula. See Theorem 4..3 o p.86 of J. Jot Partial Differetial Equatio.. Uiquee. A i the elliptic cae, uiquee i etablihed through maximum priciple.

6 Heat Equatio: Maximum Priciple Theorem. (Uiquee o bouded domai) Let g C( Ω T ), f C(Ω T ). The there exit at mot oe olutio u C (Ω T ) C(Ω T) of the iitial/boudary value problem I the cae Ω = R, we have to add extra coditio. u t u = f i Ω T ; u = g o Ω T. (46) Theorem. (Uiquee whe Ω = R ) Suppoe g C(R ), f C(R [0, T]). The there exit at mot oe olutio u C (R (0, T]) C(R [0, T]) of the iitial value problem atifyig the growth etimate for cotat M,λ>0. u t u= f i R (0, T); u= g o R t = 0} (47) u(x, t) M e λ x (48) Remark 3. There are i fact ifiitely may phyically icorrect olutio, which do ot atify the growth boud, to the zero iitial value problem. See Chapter 7 of F. Joh Partial Differetial Equatio. Exercie. Exercie. Let u C (Ω T ) C(Ω T) be a olutio to u t u=0 i Ω T ; u = 0 o Ω [0, T]; u = g o Ω t =0}, (49) where g 0. The u i poitive everywhere withi Ω T if g i poitive omewhere o Ω.