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

Similar documents
Solutions from Chapter 9.1 and 9.2

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

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

Predator - Prey Model Trajectories and the nonlinear conservation law

Differential Equations

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

Solutions to Assignment 1

556: MATHEMATICAL STATISTICS I

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

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

Class Meeting # 10: Introduction to the Wave Equation

Concourse Math Spring 2012 Worked Examples: Matrix Methods for Solving Systems of 1st Order Linear Differential Equations

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

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

Hamilton Jacobi equations

Section 3.5 Nonhomogeneous Equations; Method of Undetermined Coefficients

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.

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

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

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

An Introduction to Malliavin calculus and its applications

15. Vector Valued Functions

Chapter 2. First Order Scalar Equations

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

4.6 One Dimensional Kinematics and Integration

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

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

2. Nonlinear Conservation Law Equations

Lecture 10: Wave equation, solution by spherical means

EXERCISES FOR SECTION 1.5

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

arxiv: v1 [math.ca] 15 Nov 2016

L p -L q -Time decay estimate for solution of the Cauchy problem for hyperbolic partial differential equations of linear thermoelasticity

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

SOLUTIONS TO ECE 3084

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

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

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.

Let ( α, β be the eigenvector associated with the eigenvalue λ i

Chapter 6. Systems of First Order Linear Differential Equations

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

LAPLACE TRANSFORM AND TRANSFER FUNCTION

Chapter 3 Boundary Value Problem

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

Chapter #1 EEE8013 EEE3001. Linear Controller Design and State Space Analysis

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

Some Basic Information about M-S-D Systems

Math 36. Rumbos Spring Solutions to Assignment #6. 1. Suppose the growth of a population is governed by the differential equation.

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

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

23.5. Half-Range Series. Introduction. Prerequisites. Learning Outcomes

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

THE SINE INTEGRAL. x dt t

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

Math 221: Mathematical Notation

Evans PDE Solutions, Chapter 2

k B 2 Radiofrequency pulses and hardware

Second Order Linear Differential Equations

Two Coupled Oscillators / Normal Modes

System of Linear Differential Equations

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

4. The multiple use forestry maximum principle. This principle will be derived as in Heaps (1984) by considering perturbations

Math 23 Spring Differential Equations. Final Exam Due Date: Tuesday, June 6, 5pm

MATH 128A, SUMMER 2009, FINAL EXAM SOLUTION

7 Wave Equation in Higher Dimensions

Physics 127b: Statistical Mechanics. Fokker-Planck Equation. Time Evolution

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

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

Echocardiography Project and Finite Fourier Series

10. State Space Methods

6.2 Transforms of Derivatives and Integrals.

MA 366 Review - Test # 1

A Sharp Existence and Uniqueness Theorem for Linear Fuchsian Partial Differential Equations

NONSMOOTHING IN A SINGLE CONSERVATION LAW WITH MEMORY

Example on p. 157

Homework 10 (Stats 620, Winter 2017) Due Tuesday April 18, in class Questions are derived from problems in Stochastic Processes by S. Ross.

Physics 235 Chapter 2. Chapter 2 Newtonian Mechanics Single Particle

Lecture 20: Riccati Equations and Least Squares Feedback Control

Question 1: Question 2: Topology Exercise Sheet 3

( ) 2. Review Exercise 2. cos θ 2 3 = = 2 tan. cos. 2 x = = x a. Since π π, = 2. sin = = 2+ = = cotx. 2 sin θ 2+

Math 333 Problem Set #2 Solution 14 February 2003

Announcements: Warm-up Exercise:

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

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

Expert Advice for Amateurs

2.7. Some common engineering functions. Introduction. Prerequisites. Learning Outcomes

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

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

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.

An random variable is a quantity that assumes different values with certain probabilities.

Chapter #1 EEE8013 EEE3001. Linear Controller Design and State Space Analysis

Math 115 Final Exam December 14, 2017

Stochastic models and their distributions

CHAPTER 2 Signals And Spectra

where the coordinate X (t) describes the system motion. X has its origin at the system static equilibrium position (SEP).

Chapter 2. Motion in One-Dimension I

1 Solutions to selected problems

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

Math 2214 Solution Test 1A Spring 2016

Homogenization of random Hamilton Jacobi Bellman Equations

Transcription:

THE WAVE EQUATION 43. (S) Le u(x, ) be a soluion of he wave equaion u u xx = 0. Show ha Q43(a) (c) is a. Any ranslaion v(x, ) = u(x + x 0, + 0 ) of u(x, ) is also a soluion (where x 0, 0 are consans). par hand-in for week 9 b. Any dilaion v(x, ) = u(λx, λ) of u(x, ) is also a soluion (where λ is consan). c. Any derivaive of u(x, ) is also a soluion (assume u is as smooh as you like). d. Define he energy densiy e(x, ) by e = 2 u2 + 2 u2 x and he momenum densiy p(x, ) by p = u u x. Soluion: i. Show ha e = p x and e x = p. ii. Show ha boh e and p solve he wave equaion, i.e. e e xx = 0 and p p xx = 0. a. Clearly v xx (x, ) = u xx (x + x 0, + 0 ) and v (x, ) = u (x + x 0, + 0 ), herefore ( v)(x, ) = ( u)(x + x 0, + 0 ) = 0. (Noaion: u(x, ) = u (x, ) u xx (x, ) for -dimensional problems, and in general u(x, ) = u (x, ) 2 u(x, ) in higher dimensions.) b. v xx (x, ) = λ 2 u xx (λx, λ) and v (x, ) = λ 2 u (λx, λ), herefore ( v)(x, ) = λ 2 ( u)(λx, λ) = 0. c. To see for example ha u x is a soluion differeniae u u xx = 0 wih respec o x o ge u x u xxx = 0 and hen change he order of differeniaion (we can do his if our funcion is sufficienly smooh) o ge (u x ) (u x ) xx = 0. d. i. Assume as usual ha u is a wice coninuously differeniable soluion of he wave equaion u u xx = 0. Then ( e = 2 u2 + ) 2 u2 x = u u + u x u x ( e x = 2 u2 + ) 2 u2 x = u u x + u x u xx p = (u u x ) = u u x + u u x x and p x = (u u x ) x = u x u x + u u xx e p x = u (u u xx ) = 0 e x p = u x (u u xx ) = 0 ii. Assume now ha u is hree imes coninuously differeniable so ha we can ake second derivaives of e and p. We have e e xx = (e ) (e x ) x = p x p x = 0 and p p xx = (p ) (p x ) x = e x e x = 0

44. (S) Le u(x, ) be a real soluion of he damped wave equaion Q44 is u u xx + u = 0. (9) Show ha he energy E() = ( 2 u (x, ) 2 + 2 u x(x, ) 2) dx decreases (make any convenenien assumpions of regulariy and decay a infiniy). Soluion: Clearly we have o assume ha u is C 2 o make sense of he equaion. Moreover, we need o make sure ha all he generalized Riemann inegrals ha come in he proof exis. This will be rue if all inegrands below decay sufficienly fas a infiniy. We work exacly as in he case of he wave equaion: muliply boh sides of (9) by u o ge u u u u xx + u 2 = 0. Now u u = ( ) 2 u2 and u u xx = (u u x ) x u x u x = (u u x ) x ( ) 2 u2 x. ( ) ( ) 2 u2 (u u x ) x + 2 u2 x + u 2 = 0. Inegrae wih respec o x over R o ge (( ) 2 u2 ( ) ) (u u x ) x + 2 u2 x + u 2 dx = 0. Working as in Problem 85 we can show ha if lim u (x, )u x (x, ) = 0 for all (add his o our lis of decay x ± a infiniy assumpions, if i s no already here), hen (u u x ) x dx = 0. We are hus lef wih This can be wrien as In oher words, he energy is decreasing. (( ) ( ) 2 u2 + 2 u2 x + u 2 ) dx = 0. d ( d 2 u2 + ) 2 u2 x dx + u 2 dx = 0. de d () = u 2 dx 0. 45. (E) For a (coninuously) differeniable funcion u(x, ) define D + u and D u by compleion of hand-in for Week 9 D + u = u + u x, D u = u u x Show ha if u is wice coninuously differeniable hen D + D u = D D + u = u u xx Soluion: A sraighforward calculaion. 46. (S) Le u(x, ) be a C 2 soluion of he wave equaion in one space dimension. Le ABCD be a parallelogram on he (x, ) plane whose sides are characerisics. Using he general soluion of he wave equaion, or oherwise, show ha u(a) + u(c) = u(b) + u(d) (0) Soluion: In characerisic coordinaes ξ := x +, η := x he soluion is of he form u(ξ, η) = f(ξ) + g(η). The characerisics become ξ = cons and η = cons so he parallelogram looks like his:

y C ξ = ξ2 η = η2 B D ξ = ξ η = η A Then O u(a) = u(ξ, η ) = f(ξ ) + g(η ), u(b) = u(ξ 2, η ) = f(ξ 2 ) + g(η ) u(c) = u(ξ 2, η 2 ) = f(ξ 2 ) + g(η 2 ), u(d) = u(ξ, η 2 ) = f(ξ ) + g(η 2 ) x so boh sides of (0) are equal o f(ξ ) + f(ξ 2 ) + g(η ) + g(η 2 ). 47. (E) Le f : R 2 R and a, b : R R be coninuously differeniable funcions. Show ha d d b() f(, s)ds = b() f (, s)ds + f(, b())b () f(, )a () (Hin: Consider he funcion of hree variables F(x, y, ) = y x f(, s)ds. Then b() f(, s)ds = F(, b(), ). Use he chain rule from mulivariable calculus.) Soluion: We have Bu d d b() b() f(, s)ds = d d( F(, b(), )) f(, s)ds = F(, b(), ) = F x (, b(), )a () + F y (, b(), )b () + F (, b(), ) d d b() F x (x, y, ) = f(, x) F y (x, y, ) = f(, y) F z (x, y, ) = y x f (, s)ds b() f(, s)ds = f(, )a () + f(, b())b () + f (, s)ds

48. Le u, v : R R R be wo soluions of he wave equaion in one space dimension wih iniial daa u(x, 0) = f (x), u (x, 0) = 0, v(x, 0) = f 2 (x), v (x, 0) = 0 where f, f 2 : R R are C 2 funcions vanishing ouside some bounded inerval [ M, M]. Define Q 0, Q : R R R by Q 0 = u v u x v x, Q = u v x u x v a. (S) Express Q 0 and Q in erms of he iniial daa f and f 2. (Hin: Use d Alember s formula). b. (P) Show ha Q 0 L 2 (R R) C f L 2 (R) f 2 L 2 (R) Q L2 (R R) C f L 2 (R) f 2 L 2 (R) where C is some posiive consan independen of he iniial daa. (Hin: Use he expression you found in (a) and change o characerisic coordinaes). We have used he following noaion: For a funcion g : R R we define ( + /2 g L2 (R) = g(x) dx) 2, and similarly for a funcion G : R R R we define G L 2 (R) = ( provided ha all hese generalised Riemann inegrals exis. ) /2 G(x, ) 2 dxd [Ans (a) 2 [f (x + )f 2(x ) + f (x )f 2(x + )], 2 [f (x + )f 2(x ) f (x )f 2(x + )],] Soluion: D Alember s formula gives Using hese we ge u(x, ) = 2 (f (x + ) + f (x )), v(x, ) = 2 (f 2(x + ) + f 2 (x )) u x (x, ) = 2 (f (x + ) + f (x )) u (x, ) = 2 (f (x + ) f (x )) v x (x, ) = 2 (f 2(x + ) + f 2(x )) v (x, ) = 2 (f 2(x + ) f 2(x )) Q 0 (x, ) = u (x, )v (x, ) u x (x, )v x (x, ) = 4 [(f (x + ) f (x ))(f 2(x + ) f 2(x )) (f (x + ) + f (x )) (f 2(x + ) + f 2(x ))] = 2 [f (x + )f 2(x ) + f (x )f 2(x + )] Q 0 (x, ) 2 dxd = R 4 (x + )f 2(x ) + f (x )f 2(x + ) 2 dxd 2 R 2 2 (x + )f 2(x ) 2 dxd () R 2 + 2 (x )f 2(x + ) 2 dxd (2) R 2

where we have used (a + b) 2 2a 2 + 2b 2. We change o characerisic coordinaes ξ = x +, η = x. For he inegral in () we have: (x + )f 2(x ) 2 dxd = R 2 (ξ) 2 2(η) 2 dξdη 2 R R = ( ) ( ) 2 (ξ) 2 dξ 2(η) 2 dη R R Similarly he inegral in (2) is: herefore = 2 f L 2 (R) 2 f 2 L 2 (R) 2 R 2 (x )f 2 (x + ) 2 dxd = 2 f L2 (R) 2 f 2 L2 (R) 2 R 2 Q 0 (x, ) 2 dxd 2 f L2 (R) 2 f 2 L2 (R) 2 The proof for Q is similar. 49. (S) Le g : R 3 R be given by g(x) = g 0 (r) = r 2. We wish o solve he following IVP for he wave equaion in hree space dimensions wih radial daa: We sar by guessing ha u u = 0, u(x, 0) = 0, u (x, 0) = g(x) = g 0 (r) u(x, ) = u 0 (r, ) where u 0 : R R R is even in r i.e. u 0 ( r, ) = u 0 (r, ) (why is his a good guess?). a. Show ha he wave equaion becomes r(u 0 ) r(u 0 ) rr 2(u 0 ) r = 0 b. Show ha his means ha w(r, ) = ru 0 (r, ) solves he wave equaion in one space dimension wih iniial daa w(r, 0) = 0 and w (r, 0) = rg 0 (r). c. Use d Alember s formula o find w. w w rr = 0 d. Conclude ha u 0 (r, ) = 2r r+ r ρg 0 (ρ)dρ e. Find u. [Ans: (c) r 3 + r 3, (e) x 2 + 3 ] Soluion: Since he iniial daa are radial funcions (hey depend only on r = x ) we guess ha he soluion will also be radial (his can acually be proved), say u(x, ) = u 0 ( x, ) = u 0 (r, ), for some funcion u 0 (r, ) defined for r 0, R. Since he iniial daa are smooh even funcions of r we guess ha u 0 could be exended o all r s as a smooh even funcion of r, i.e. u 0 (r, ) = u 0 ( r, ). Using problem 32(a) wih n = 3 we have u u = (u 0 ) [(u 0 ) rr + 2r ] (u 0) r herefore, for all r 0 and R we have r(u 0 ) r(u 0 ) rr 2(u 0 ) r = 0 This is saisfied for r < 0 oo. Indeed, using u 0 ( r, ) = u 0 (r, ), we have (u 0 ) r (r, ) = (u 0 ) r ( r, ) (u 0 ) rr (r, ) = (u 0 ) rr ( r, ) (u 0 ) (r, ) = (u 0 ) ( r, )

for r < 0, R we have r(u 0 ) (r, ) r(u 0 ) rr (r, ) 2(u 0 ) r (r, ) = = r(u 0 ) ( r, ) r(u 0 ) rr ( r, ) + 2(u 0 ) r ( r, ) = [( r)(u 0 ) ( r, ) ( r)(u 0 ) rr ( r, ) (u 0 ) r ( r, )] = 0 Define w : R R R by w(r, ) = ru 0 (r, ). Then w w rr = r(u 0 ) r(u 0 ) rr 2(u 0 ) r = 0 i.e. w solves he wave equaion in one space dimension wih iniial daa w 0 (r, 0) = ru 0 (r, 0) = 0, (w 0 ) (r, 0) = r(u 0 ) (r, 0) = r 3 By d Alember s formula hence w(r, ) = 2 u 0 (r, ) = r+ r w(r, ) r ρ 3 dρ = r 3 + r 3 = r 2 + 3 u(x, ) = x 2 + 3 50. (S) Le u(x, ) be a soluion of he hree-dimensional wave equaion, and le v(, x) = u(, x + x 0 ). Verify ha v also saisfies he wave equaion. If v(x, 0) = 0 and v (x, 0) = g(x), x R 3 i is shown in class ha v(0, ) = 4π y = g(y)ds. Deduce ha u(x 0, ) = u (y, 0)dS. 4π y x 0 = Soluion: Wih v i (x, ) = x i v(x, ), ec, we have Hence resul. Nex, v (x, 0) = u (x + x 0, 0) = g(x). Hence v i (x, ) = u i (x + x 0, ) v ii (x, ) = u ii (x + x 0, ) v(x, ) = u(x + x 0, ) +v (x, ) = u (x, ). u(x 0, ) = v(0, ) = 4π = 4π = 4π y = y = g(y)ds y x 0 = u (y + x 0, 0)dS u (y, 0)dS where y = y + x 0.