Hamilton Jacobi equations

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- J acobi Equation: Weak S olution We continue the study of the Hamilton-Jacobi equation:

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

An Introduction to Malliavin calculus and its applications

Math 334 Fall 2011 Homework 11 Solutions

arxiv: v2 [math.ap] 16 Oct 2017

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...

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

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).

Chapter 2. First Order Scalar Equations

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

Predator - Prey Model Trajectories and the nonlinear conservation law

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

1 Solutions to selected problems

2. Nonlinear Conservation Law Equations

Example on p. 157

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

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

Differential Equations

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

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

556: MATHEMATICAL STATISTICS I

Solutions to Assignment 1

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)

ENGI 9420 Engineering Analysis Assignment 2 Solutions

LAPLACE TRANSFORM AND TRANSFER FUNCTION

Algorithmic Trading: Optimal Control PIMS Summer School

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

Expert Advice for Amateurs

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

Homogenization of random Hamilton Jacobi Bellman Equations

Theory of! Partial Differential Equations-I!

Laplace Transforms. Examples. Is this equation differential? y 2 2y + 1 = 0, y 2 2y + 1 = 0, (y ) 2 2y + 1 = cos x,

Lecture 20: Riccati Equations and Least Squares Feedback Control

Theory of! Partial Differential Equations!

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

Class Meeting # 10: Introduction to the Wave Equation

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

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

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

Math 333 Problem Set #2 Solution 14 February 2003

Lecture 10: The Poincaré Inequality in Euclidean space

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

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

EXERCISES FOR SECTION 1.5

Section 3.5 Nonhomogeneous Equations; Method of Undetermined Coefficients

Convergence of the Neumann series in higher norms

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.

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

Hamilton-Jacobi Equations

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.

Y 0.4Y 0.45Y Y to a proper ARMA specification.

Math Final Exam Solutions

Solutions Problem Set 3 Macro II (14.452)

SOLUTIONS TO ECE 3084

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

u(x) = e x 2 y + 2 ) Integrate and solve for x (1 + x)y + y = cos x Answer: Divide both sides by 1 + x and solve for y. y = x y + cos x

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

THE 2-BODY PROBLEM. FIGURE 1. A pair of ellipses sharing a common focus. (c,b) c+a ROBERT J. VANDERBEI

Basic Circuit Elements Professor J R Lucas November 2001

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

4.5 Constant Acceleration

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

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

STABILITY OF NONLINEAR NEUTRAL DELAY DIFFERENTIAL EQUATIONS WITH VARIABLE DELAYS

15. Vector Valued Functions

THE BERNOULLI NUMBERS. t k. = lim. = lim = 1, d t B 1 = lim. 1+e t te t = lim t 0 (e t 1) 2. = lim = 1 2.

THE MYSTERY OF STOCHASTIC MECHANICS. Edward Nelson Department of Mathematics Princeton University

= ( ) ) or a system of differential equations with continuous parametrization (T = R

Admissible Solution for Hyperbolic Conservation Laws

Lie Derivatives operator vector field flow push back Lie derivative of

( ) a system of differential equations with continuous parametrization ( T = R + These look like, respectively:

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

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

QUANTITATIVE DECAY FOR NONLINEAR WAVE EQUATIONS

Simulation-Solving Dynamic Models ABE 5646 Week 2, Spring 2010

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

Linear Surface Gravity Waves 3., Dispersion, Group Velocity, and Energy Propagation

DISCRETE GRONWALL LEMMA AND APPLICATIONS

4.6 One Dimensional Kinematics and Integration

Semiconcave Functions, Hamilton Jacobi Equations, and Optimal Control. Piermarco Cannarsa & Carlo Sinestrari

Fréchet derivatives and Gâteaux derivatives

Online Appendix to Solution Methods for Models with Rare Disasters

arxiv: v1 [math.ca] 15 Nov 2016

Chapter 7 Response of First-order RL and RC Circuits

6.2 Transforms of Derivatives and Integrals.

Backward stochastic dynamics on a filtered probability space

NONSMOOTHING IN A SINGLE CONSERVATION LAW WITH MEMORY

A Note on Superlinear Ambrosetti-Prodi Type Problem in a Ball

(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)

ME 391 Mechanical Engineering Analysis

Evans PDE Solutions, Chapter 2

8. Basic RL and RC Circuits

Physics 235 Chapter 2. Chapter 2 Newtonian Mechanics Single Particle

Viscosity solutions of Hamilton-Jacobi equations - summer course at Universidade Federal Fluminense, Niterói. Diogo Aguiar Gomes

FEEDBACK NULL CONTROLLABILITY OF THE SEMILINEAR HEAT EQUATION

Math 10B: Mock Mid II. April 13, 2016

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

System of Linear Differential Equations

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

Transcription:

Hamilon Jacobi equaions Inoducion o PDE The rigorous suff from Evans, mosly. We discuss firs u + H( u = 0, (1 where H(p is convex, and superlinear a infiniy, H(p lim p p = + This by comes by inegraion from special hyperbolic sysems of he form (n = m v + F j (v j v = 0 when here exiss a ponenal for F j, i.e. F j = j H(v, and when we seek soluions as gradiens v = u. Typical example (n = 1 is Burgers equaion v + vv x = 0 which by inegraion v = u x gives, omiing consans, u + 1 2 u2 x = 0. We know soluions of Burgers or conservaion laws conserve L norms and produce shocks. So we expec soluions of HJ equaions o be Lipschiz, and loose second derivaives. Tha hey do. We recall ha he Legendre ransform of a convex superlinear-a-infiniy funcion L is L (p = sup q R n [p q L(q] Noe ha L is convex and superlinear a infiniy and (L = L. 1

Le H = L. Noe ha, if L is smooh (and i is, if H is sricly convex, hen a a maximum (hey are aained becuse of superlineariy q L(q = p, so, solving for q gives a funcion of p. (This is unique if L is sricly convex. So, H(p = p q(p L(q(p Differeniaing his, we have p H(p = q(p + p p q q L(q(p p q = q(p This can be undersood: he Legendre ransform of L is given by H(p = p q(p L(q(p where q solves q L = p, or by he equaion H(p = p p H L( p H(p which looks horrible, and iself i is a seady Hamilon-Jacobi equaion. Similarly, H (q = L(q is compued by solving p H(p = q and seing L(q = p(q q H(p(q. Of course, if L is no sricly convex we do no have necessarily a unique soluion, and he minimizaion is necessary. Now we have a magical soluion of (1, u(x, = min y { L ( x y } + u 0 (y where u 0 is he value of u a = 0. This is ermed he Hopf-Lax formula in Evans. The magic is hen enveloped in mysery: {ˆ } u(x, = min L(ẇ(sds + u 0 (w(0 A where 0 A = {w C 1 ([0, ], R n w( = x} and ẇ(s = dw. This is easier o prove han i looks (see Evans, and makes ds conac wih minimum acion principles, or conrol, bu does no illuminae (2 2

he Hopf Lax formula by even a single lumen. By wha divinaion process could hey have arrived a such an amazing formula? Le us differeniae our equaion (1, moivaed by where i came from. We obain v + ( p H(v v = 0 where v = u. So, le us inroduce characerisics, dx d = ph(v(x,, x(0 = x 0. Denoe for a second x(, x 0 he characerisic issued from x 0. On i, v is consan, v(x(, x 0, = v(x 0, 0. Bu v = u, which is he sole argumen of p H. So he characerisics are sraigh lines and x(, x 0 = x 0 + q, u(x, = p, q = p H(p for x = x 0 + p H(p. Now le us look a u on he characerisic. We differeniae d ds u(x 0 + sq, s = s u(x 0 + sq, s + q u(x 0 + sq, s Bu we are on he characerisic, so u = p, and q = p H(p. Using he equaion (1 we also have s u(x 0 + sq, s = H(p. The righ-hand side is a consan, p q H(p. Tha we know how o inegrae in ime: u(x, u(x 0, 0 = (p q H(p We recall ha q = p H(p. Wai a minue, ha means ha p q H(p = L(q. So u(x, = L(q + u 0 (x 0 Now we express q in erms of x and x 0 : From x = x 0 + q, i follows ha q = x x 0. So, wihou magic, we arrived a ( x x0 u(x, = L + u 0 (x 0. Now, we assumed sric convexiy, and smooh iniial daa, and his would work only for shor ime. Remarkably, he variaional formula is rue for all ime and gives he good weak soluion even afer u develops shocks. When hey do, x 0 is no uniquely deermined on characerisics by u(x,, x an. We sar wih a semigroup propery: 3

Lemma 1. If H is convex and superlinear a infiniy and if u is given by (2, hen { ( } x y u(x, = min ( sl + u(y, s (3 y s for 0 s <. Proof. Take z so ha u(y, s = sl ( y z s + u0 (z and use convexiy of L o show ( x z ( L 1 s ( x y L + s ( y z s L s Then u(x, L ( x z = ( sl ( x y + u0 (z ( s L ( x y s ( + sl y z s + u0 (z s + u(y, s. On he oher hand, pick w so ha u(x, = L ( x w + u0 (w and ake y = sx + (1 sw. Then x y s = x w = y w s and consequenly ( sl ( ( x y s + u(y, s ( sl x w ( + sl y w s + u0 (w = L ( x w + u0 (w = u(x,. Proposiion 1. If u 0 is Lipschiz, hen he funcion u(x, given in (2 is Lipschiz coninuous and coninuous up o = 0. Proof. Fix > 0. Pick x, x 1. Find y so ha ( x y u(x, = L + u 0 (y Then { ( u(x 1, u(x, = min z L x1 z + u0 (z } L ( x y u0 (y g(x 1 x + y g(y C x x 1. Then we swich he roles of x 1 and x. For he coninuiy a = 0, on one hand we have u(x, L(0 + u 0 (x 4

direcly from (2 by aking y = x. On he oher hand, using u 0 (y u 0 (x C x y in he definiion (2 we have { ( u(x, u 0 (x + min y C x y + L x y } = u 0 (x max w {C w L(w} = u 0 (x max z C max w {z w L(w} = u 0 (x max z C H(z We also have u(x, u(x, s C s using he semigroup propery and he bounds above. Proposiion 2. Le u 0 be Lipschiz and ake a poin (x, where he funcion u(x, defined by (2 is differeniable. (This happens a.e. by Rademacher s heorem and he resul above. Then u(x, + H( u(x, = 0. Proof. le q R n, h > 0. By he semigroup propery { ( u(x + hq, + h = min y hl x+hq y } h + u(y, hl(q + u(x,, which implies This is valid for all q, so q u(x, + u(x, L(q u(x, + H( u(x, = u(x, + min {q u(x, L(q} 0. q On he oher hand, choose z so ha u(x, = L ( x z + u0 (z. Take h > 0, pu s = h, y = sx + (1 s x z z. Then = y z and so, s u(x, u(y, s L ( x z + u0 (z { sl ( y z s + u0 (z } = ( sl ( x z. This means 1 h Leing h 0 ( (( u(x, u 1 h x + h z, h L u(x, + x z ( x z u(x, L 5 ( x z

Then u(x, + H( u(x, = u(x, + max w {q u(x, L(q} u(x, + x z u(x, L ( x z 0. This concludes his verificaion. Now i urns ou ha an addiional propery holds, semi-concaviy. Lemma 2. If here exiss a consan such ha u 0 (x + z + u 0 (x z 2u 0 (x C z 2 holds for all x, z, hen he soluion (2 saisfies, wih he same C u(x + z, + u(x z, 2u(x, C z 2 Proof. Le y so ha u(x, = L ( x y + u0 (y and use his y + z and y z in he definiions of u(x + z,, u(x z,. This will cancel he L erms and give he resul. Semi-concaviy of u implies ha u ɛ = φ ɛ u where φ ɛ is a sandard mollifier saisfies u ɛ CI I urns ou ha if H is sricly convex, hen for > 0 soluions become semiconcave, even if u 0 was no, wih bounds ha explode a = 0. Lemma 3. If H is sricly convex ( p1 + p 2 H 1 2 2 H(p 1 + 1 2 H(p 2 c p 1 p 2 2 hen he soluion (2 saisfies u(x + z, + u(x z, 2u(x, C z 2. Proof. Using definiions, i urns ou ha 1 2 L(q 1 + 1 ( 2 L(q q1 + q 2 2 L + C q 1 q 2 2 2 Then, choosing y o have u(x, = L ( x y + u0 (y and using he same y for u(x + z, and u(x z, he u 0 erms cancel ou, and he resul emerges. 6

Theorem 1. Le H C 2, sricly convex, superlinear a infiniy. Assume u 0 is Lipschiz. Then here is only one Lipschiz coninuous funcion u saisfying u(x, 0 = u 0 (x and u + H( u(x, = 0, a.e. u(x + z, + u(x z, 2u(x C(1 + 1 z 2 Proof. Taking he difference u of wo soluions u 1, u 2 we arrive a where v is bounded, and given by v(x, = ˆ 1 0 u + v u = 0 We mollify he soluions u 1, u 2 and define v ɛ (x, = ˆ 1 0 p H(λ u 1 (x, + (1 λu 2 (x, dλ p H(λ u 1ɛ (x, + (1 λu 2ɛ (x, dλ where u jɛ = u j φ ɛ. The divergence saisfies a one-sided bound divv ɛ C(1 + 1. This used he semiconcaviy of boh soluions. Now because boh soluions are Lipschiz u(x, c, u C. We wrie u + v ɛ u = (v ɛ v u We ake a funcion w = f(u wih f smooh and nonnegaive. This sill solves w + v ɛ w = (v ɛ v w We inegrae on a fasly shrinking ball ˆ e( = x x 0 V ( 0 7 wdx

wih V > b L. Then de d C(1 + 1 e + ˆ x x 0 V ( 0 (v v ɛ w We le ɛ 0: de d C(1 + 1 e Now we selec f o vanish for u cɛ, posiive oherwise. Then e = 0 for ɛ and from hen on, by Gronwall e( = 0. I follows ha u cɛ, bu ɛ was arbirary. 8