Hamilton- J acobi Equation: Weak S olution We continue the study of the Hamilton-Jacobi equation:

Similar documents
Hamilton- J acobi Equation: Explicit Formulas In this lecture we try to apply the method of characteristics to the Hamilton-Jacobi equation: u t

Hamilton Jacobi equations

Math 527 Lecture 6: Hamilton-Jacobi Equation: Explicit Formulas

KEY. Math 334 Midterm III Winter 2008 section 002 Instructor: Scott Glasgow

Chapter 2. First Order Scalar Equations

2. Nonlinear Conservation Law Equations

KEY. Math 334 Midterm III Fall 2008 sections 001 and 003 Instructor: Scott Glasgow

Finish reading Chapter 2 of Spivak, rereading earlier sections as necessary. handout and fill in some missing details!

Expert Advice for Amateurs

Math 2142 Exam 1 Review Problems. x 2 + f (0) 3! for the 3rd Taylor polynomial at x = 0. To calculate the various quantities:

t is a basis for the solution space to this system, then the matrix having these solutions as columns, t x 1 t, x 2 t,... x n t x 2 t...

MATH 4330/5330, Fourier Analysis Section 6, Proof of Fourier s Theorem for Pointwise Convergence

EXERCISES FOR SECTION 1.5

t + t sin t t cos t sin t. t cos t sin t dt t 2 = exp 2 log t log(t cos t sin t) = Multiplying by this factor and then integrating, we conclude that

Inventory Analysis and Management. Multi-Period Stochastic Models: Optimality of (s, S) Policy for K-Convex Objective Functions

Solutions from Chapter 9.1 and 9.2

Solutions to Assignment 1

THE WAVE EQUATION. part hand-in for week 9 b. Any dilation v(x, t) = u(λx, λt) of u(x, t) is also a solution (where λ is constant).

Essential Microeconomics : OPTIMAL CONTROL 1. Consider the following class of optimization problems

Vanishing Viscosity Method. There are another instructive and perhaps more natural discontinuous solutions of the conservation law

Example on p. 157

An Introduction to Malliavin calculus and its applications

arxiv: v2 [math.ap] 16 Oct 2017

MATH 128A, SUMMER 2009, FINAL EXAM SOLUTION

arxiv: v1 [math.ca] 15 Nov 2016

Homogenization of random Hamilton Jacobi Bellman Equations

on the interval (x + 1) 0! x < ", where x represents feet from the first fence post. How many square feet of fence had to be painted?

Predator - Prey Model Trajectories and the nonlinear conservation law

Physics 235 Chapter 2. Chapter 2 Newtonian Mechanics Single Particle

Math 2214 Solution Test 1A Spring 2016

Section 3.5 Nonhomogeneous Equations; Method of Undetermined Coefficients

1 Solutions to selected problems

A Necessary and Sufficient Condition for the Solutions of a Functional Differential Equation to Be Oscillatory or Tend to Zero

Math Week 14 April 16-20: sections first order systems of linear differential equations; 7.4 mass-spring systems.

Positive continuous solution of a quadratic integral equation of fractional orders

(1) (2) Differentiation of (1) and then substitution of (3) leads to. Therefore, we will simply consider the second-order linear system given by (4)

On a Fractional Stochastic Landau-Ginzburg Equation

An Introduction to Backward Stochastic Differential Equations (BSDEs) PIMS Summer School 2016 in Mathematical Finance.

Lecture 2-1 Kinematics in One Dimension Displacement, Velocity and Acceleration Everything in the world is moving. Nothing stays still.

Notes for Lecture 17-18

Math Final Exam Solutions

Optimality Conditions for Unconstrained Problems

4 Sequences of measurable functions

Boundedness and Exponential Asymptotic Stability in Dynamical Systems with Applications to Nonlinear Differential Equations with Unbounded Terms

Oscillation of an Euler Cauchy Dynamic Equation S. Huff, G. Olumolode, N. Pennington, and A. Peterson

HOMEWORK # 2: MATH 211, SPRING Note: This is the last solution set where I will describe the MATLAB I used to make my pictures.

MODULE 3 FUNCTION OF A RANDOM VARIABLE AND ITS DISTRIBUTION LECTURES PROBABILITY DISTRIBUTION OF A FUNCTION OF A RANDOM VARIABLE

Hybrid Control and Switched Systems. Lecture #3 What can go wrong? Trajectories of hybrid systems

The Asymptotic Behavior of Nonoscillatory Solutions of Some Nonlinear Dynamic Equations on Time Scales

dt = C exp (3 ln t 4 ). t 4 W = C exp ( ln(4 t) 3) = C(4 t) 3.

Homework sheet Exercises done during the lecture of March 12, 2014

dy dx = xey (a) y(0) = 2 (b) y(1) = 2.5 SOLUTION: See next page

Matrix Versions of Some Refinements of the Arithmetic-Geometric Mean Inequality

Stability and Bifurcation in a Neural Network Model with Two Delays

A universal ordinary differential equation

Module 2 F c i k c s la l w a s o s f dif di fusi s o i n

Course Notes for EE227C (Spring 2018): Convex Optimization and Approximation

Math 334 Fall 2011 Homework 11 Solutions

14 Autoregressive Moving Average Models

GMM - Generalized Method of Moments

E β t log (C t ) + M t M t 1. = Y t + B t 1 P t. B t 0 (3) v t = P tc t M t Question 1. Find the FOC s for an optimum in the agent s problem.

t 2 B F x,t n dsdt t u x,t dxdt

23.2. Representing Periodic Functions by Fourier Series. Introduction. Prerequisites. Learning Outcomes

SMT 2014 Calculus Test Solutions February 15, 2014 = 3 5 = 15.

Theory of! Partial Differential Equations-I!

Mon Apr 9 EP 7.6 Convolutions and Laplace transforms. Announcements: Warm-up Exercise:

f(s)dw Solution 1. Approximate f by piece-wise constant left-continuous non-random functions f n such that (f(s) f n (s)) 2 ds 0.

Solutions of Sample Problems for Third In-Class Exam Math 246, Spring 2011, Professor David Levermore

Some Basic Information about M-S-D Systems

15. Vector Valued Functions

Question 1: Question 2: Topology Exercise Sheet 3

ODEs II, Lecture 1: Homogeneous Linear Systems - I. Mike Raugh 1. March 8, 2004

Approximation Algorithms for Unique Games via Orthogonal Separators

TO our knowledge, most exciting results on the existence

INSTANTANEOUS VELOCITY

The Contradiction within Equations of Motion with Constant Acceleration

MA 214 Calculus IV (Spring 2016) Section 2. Homework Assignment 1 Solutions

Properties Of Solutions To A Generalized Liénard Equation With Forcing Term

ASYMPTOTIC FORMS OF WEAKLY INCREASING POSITIVE SOLUTIONS FOR QUASILINEAR ORDINARY DIFFERENTIAL EQUATIONS

1 1 + x 2 dx. tan 1 (2) = ] ] x 3. Solution: Recall that the given integral is improper because. x 3. 1 x 3. dx = lim dx.

Uniqueness of solutions to quadratic BSDEs. BSDEs with convex generators and unbounded terminal conditions

Challenge Problems. DIS 203 and 210. March 6, (e 2) k. k(k + 2). k=1. f(x) = k(k + 2) = 1 x k

Two Coupled Oscillators / Normal Modes

1. Introduction The present paper is concerned with the scalar conservation law in one space dimension: u t + (A(u)) x = 0 (1.0.1)

18 Biological models with discrete time

Time discretization of quadratic and superquadratic Markovian BSDEs with unbounded terminal conditions

KEY. Math 334 Midterm I Fall 2008 sections 001 and 003 Instructor: Scott Glasgow

Chapter 6. Systems of First Order Linear Differential Equations

Ann. Funct. Anal. 2 (2011), no. 2, A nnals of F unctional A nalysis ISSN: (electronic) URL:

4. Advanced Stability Theory

Solutions Problem Set 3 Macro II (14.452)

Lecture 20: Riccati Equations and Least Squares Feedback Control

SOLUTIONS TO ASSIGNMENT 2 - MATH 355. with c > 3. m(n c ) < δ. f(t) t. g(x)dx =

Guest Lectures for Dr. MacFarlane s EE3350 Part Deux

WEEK-3 Recitation PHYS 131. of the projectile s velocity remains constant throughout the motion, since the acceleration a x

An introduction to evolution PDEs November 16, 2018 CHAPTER 5 - MARKOV SEMIGROUP

Math 334 Test 1 KEY Spring 2010 Section: 001. Instructor: Scott Glasgow Dates: May 10 and 11.

Ordinary Differential Equations

Comparison between the Discrete and Continuous Time Models

Undetermined coefficients for local fractional differential equations

Transcription:

M ah 5 7 Fall 9 L ecure O c. 4, 9 ) Hamilon- J acobi Equaion: Weak S oluion We coninue he sudy of he Hamilon-Jacobi equaion: We have shown ha u + H D u) = R n, ) ; u = g R n { = }. ). In general we canno expec he exisence of classical soluions ha is u C R n, ) ) saisfying he equaion everywhere) ;. The Hopf-Lax formula saisfies he equaion almos everywhere. { x y u x, ) = min L y ) } + g y) We see ha exisence is guaraneed as soon as we only require he equaion o be saisfied almos everywhere. Thus i is naural o propose he following as he definiion for weak soluions : u Lipschiz is a weak soluion if u saisfies he equaion a. e. and akes g as is iniial value. However such definiion leads o non-uniqueness, as shown in he following example. Example. Consider he D H-J equaion u + u x = R, ) ; u = R { = }. 3) I is clear ha u is a classical soluion. However, we can ry o find oher soluions using separaion of variables. Se u = T ) + X x). Then he equaion gives ) which gives T ) + X x) ) = 4) T ) = λ, X x) ) = λ. 5) Obviously if we require he above o hold everywhere, hen necessarily λ = Recall Exercise 5 in Homework ). However, if we only require he equaion o be saisfied almos everywhere, hen here are infiniely many soluions. More specifically, we can pick any X such ha X = ± λ piecewise, and hen solve T = λ. We can even ake differen λ for differen x,. One example of almos everywhere soluions obained his way is x u x, ) = x x 6) x x { x x Noe ha u x, ) is in fac assembled using wo soluions : and. I is easy o see ha x x one can assemble any soluions in anyway o ge a funcion saisfying he equaion almos everywhere, he only hing o be careful is o make he resul be consisen wih he iniial value. For example: + gives a almos everywhere soluion noe ha T so he wedge is moving downward, hus saisfying he iniial value). In conras,

is no a good consrucion, as i does no saisfy he iniial value. I is clear ha he soluion is no unique. How do we regain uniqueness? Or more specifically, how o do we ge rid of all non-zero soluions as u is obviously a reasonable soluion for his problem). One way is o add viscosiy. Consider u ε + u ε x ε = ε u x x 7) and hen le ε. In many cases u ε will hen converge o one of he many almos everywhere soluions of he inviscid problem. In oher words, u ε is a smoohed version of he correc soluion. Smoohing he wedge soluion, we see ha x = is a minimizer of u ε. A his poin, we have u <, u x =, ε u x x > 8) a conradicion! Therefore he limi of u ε canno be he wedge soluion! From he above example we see ha wedges should no appear in he correc weak soluion. In oher words, a correc soluion should no have u = + anywhere. Moivaed by his, we inroduce he following noion. Definiion. Semi-concaviy) f is called semi- concave if here is C > such ha for any x, z. f x + z) f x) + f x z) C z 9) I is easy o see ha he condiion is equivalen o f x) C x ) is concave. Hence he name semi-concave. Now we can ry he following definiion of weak soluions for he H-J equaion: Definiion 3. A Lipschiz funcion u is said o be a weak soluion if. u akes g as iniial value;. u saisfies he equaion a. e. ; 3. u is semi- concave: u x + z, ) u x, ) + u x z, ) C + ) z. ) Now we ry o esablish well-posedness. More specifically, we will show ha under cerain assumpions, he soluion given by he Hopf-Lax formula is he unique weak soluion. Theorem 4. Exisence) The funcion u x, ) given by he Hopf-Lax formula is a weak soluion if eiher one of he following is rue. a ) g is semi-concave; b ) H is uniformly convex, ha is for all p, ξ. ξ T D H ξ θ ξ. Of course, for his problem, one can direcly show ha u ε. However he argumen used here is more appropriae for moivaing he general heory. )

Proof. a) Since g is semi-concave, here is C > such ha g x + z) g x) + g x z) C z. 3) Now for any, le y be such ha ) x y u x, ) = L + g y). 4) Then we have 3 ) x y u x ± z, ) L + g y ± z) and he conclusion follows. 6) b) Again we choose y such ha ) x y u x, ) = L + g y). 7) Now ) x ± z y u x ± z, ) L + g y). 8) Thus [ ) ) ) ] x + z y x y x z y u x + z, ) u x, ) + u x z, ) L L + L. 9) The conclusion follows from he fac ha ) p + p H uniformly convex H ) q + q L Which is Problem 3. 5. 9. 4 Now we esablish uniqueness. H p ) H p ) θ 8 p p L q ) L q ) 8 θ q q. ) Theorem 5. Uniqueness) Assume H is C, convex, wih super-linear growh a infiniy and g Lipschiz. Then here exiss a mos one weak soluion of he H-J equaion. Proof. Suppose u, ũ are wo weak soluions. hen leing w = u ũ, we have w is differeniable a. e. and Leing w x, ) = u x, ) ũ x, ) = H D ũ x, ) ) H D u x, ) ) d = H r D u + r) D ũ) dr dr = D H r D u + r) Dũ) D w x, ) dr [ ] = D H r D u + r) Dũ) dr D w x, ). ) b x, ) D H r D u + r) Dũ) dr ). Equivalenly, H θ p is sill convex. 3. By he Hopf-Lax formula Take y = y ± z. { x ± z y u x ± z, ) = min L y 4. Hin: For any q, q, ake p, p such ha H p i) + L q i) = p i q i. ) } + g y ) 5)

we conclude ha w saisfies w + b D w = a. e. ; w = =. 3) If w is a classical soluion o he above equaion, hen obviously w and we are done. Bu his is no he case. Neverheless, if b is Lipschiz, we sill can show w as follows. Le x ) be such ha x = b. Then w saisfies d w x ), ) = a. e. 4) d As w is Lipschiz in x, and x ) Lipschiz in, W ) w x ), ) is also Lipschiz and hus absoluely coninuous. As a consequence we have W ) W ) = W s) ds = 5) as W s) = almos everywhere. Now since W ) = we conclude ha W. Then i s clear ha w = almos everywhere. However he problem is ha b x, ) D H r D u + r) Dũ) dr 6) is definiely no Lipschiz. The idea now is o approximae b by a Lipschiz funcion b ε, ge some esimae and hen le ε. Then we can wrie w + b ε D w = b ε b) D w a. e. 7) If we ry o do energy-ype esimae of he above, we quickly realize ha a uniform upper bound of b ε is needed. This means ha he naïve way b ε = η ε b 8) will no work as his way b ε /ε. The correc smoohing here is we use subscrip o emphasize ha his is no he usual mollificaion) b ε D H r D u ε + r) D ũ ε ) dr 9) where u ε, ũ ε are usual mollificaion f ε = η ε f. This way we have [ b ε = D H: r D u ε + r) D ũ ] ε 3) where A: B ) i, j A i j B i j. Now using he fac ha boh u, ũ are semi-concave, we have D u ε, D ũ ε C + I. Noe ha as he mollificaion is aken over space-ime, such bound only holds for > C ε. Using his we have b ε C + ). 3 ) Now le v φ w) o be fixed laer. Then v saisfies he same equaion as w. Se R max { D H p) p max { Lip u), Lip ũ) } }. 3) We ry o show ha v in he cone 5 Se C { x, ), x x R ) }. 33) e ) B x, R ) ) v x, ) dx. 34) 5. Recall our proof of he uniqueness for he wave equaion! Here he slope R is an upper bound of he speed of characerisics.

We compue e ) = v R v ds B x, R ) ) = b ε v) + b ε ) v + b ε b) D v dx R v ds B B x, R ) ) = v b ε n + R) ds + b ε ) v + b ε b) D v. 35) B B ) The firs erm is due o he definiions of b ε and R. The second erm is bounded by C + e ), he hird erm vanishes as ε due o dominaed convergence. Thus we have e ) C + ) e ) 36) for a. e. < <. Finally, aking φ such ha φ z) = for z ε [ Lip u) + Lip ũ) ] and posiive oherwise, we have v = φ w) = 37) for ε. Now for > ε, + + and we have ε e ) C + ) e ) a. e., > ε; e ) = ε. 38) ε which gives e ) =. This means u ũ = w ε [ Lip u) + Lip ũ) ] 39) in he cone C. By he arbirariness of ε we conclude ha u ũ = almos everywhere and ends he proof.