Problems from Sections 2.2 - 2.5 in Logan

Problems from Sections 2.2 - 2.5 in Logan Problems from Sections 2.2 - 2.5 in Logan

math.utah.edu
from math.utah.edu More from this publisher
28.04.2014 Views

one family of such curves that carry signals and provide a distinguished these roles. Even first order PDEs, which are actually wave-like have perbolic problems there is always a set of characteristic curves that play influence and domain of dependence. —c2U, = 0, x > 0, — < t < oo; u(0, t) = s(t), — < t < 3. Solve the outgoing signal problem tions F and Gin the general solution (2.11) using the initial conditions (2.13). impulse. Contrast the solution with the case when f 0 and g = 0. Region of the solution surface and discuss the effect of giving a string at rest an initial Figure 2.3. displacement isf(x) = 0, and the initial velocity isg(x) = 1/(1 + 0.25x 2). Plot 2. Calculate the exact solution to the Cauchy problem when c = 2, the initial 1. Derive d’Alembert’s formula (2.14) by determining the two arbitrary func So the signal breaks into two pieces, and they propagate in opposite direc So the signal stays at the same place, but it spreads out and decreases in ill + 4kt u(x, t) = e2’(1+41t) 1 via u(x, t) = u(x, t) = O5(e_(_Ct) 2 + ex4ct) 2). Think of this signal as being a bit of information. case. For example, suppose the initial signal is a Gaussian function or bellshaped curve exp(—x gate information is to determine how a signal is propagated in a special is retained. Parabolic, or diffusion, equations propagate signals at infinite herency in the wave form as it propagates, and therefore information parabolic and hyperbolic problems. Hyperbolic, or wave-like, equations The convection equation Ut + CU = 0, which is a wave-like equation, amplitude. Any information in the signal is eventually lost. :ion at (xo, to). This rm u 0. In hy Finally, we point out again the important differences between speed; because the signals diffuse or smear out, there is a gradual loss of tions at speed c. The diffusion equation Ut = ku propagates the signal Section 1.2). coordinate sytem where the problem simplifies (recall the examples in propagate signals at a finite speed along characteristics; there is co information. A good way to understand how different equations propa propagates this signal via That is, it moves it at speed c without distortion. The wave equation =c2u moves it via 2). )lution outside this the wave equation Exercises fine a special coor ward in time to the regarding the charsignals forward in )und by tracing the y the initial values one. Looking at the itive characteriSticS, then the region of constant are paths affect the solution —time along which ence of the interval a and x — ct always zero outside :hat interval, so the n surface. :h velocity c. These z.z. The Jnbounded Domains Exercises 57

one family of such curves that carry signals and provide a dist<strong>in</strong>guished<br />

these roles. Even first order PDEs, which are actually wave-like have<br />

perbolic problems there is always a set of characteristic curves that play<br />

<strong>in</strong>fluence and<br />

doma<strong>in</strong> of<br />

dependence.<br />

—c2U, = 0, x > 0, —<br />

< t < oo; u(0, t) = s(t), —<br />

< t <<br />

3. Solve the outgo<strong>in</strong>g signal problem<br />

tions F and G<strong>in</strong> the general solution (2.11) us<strong>in</strong>g the <strong>in</strong>itial conditions (2.13).<br />

impulse. Contrast the solution with the case when f 0 and g = 0.<br />

Region of<br />

the solution surface and discuss the effect of giv<strong>in</strong>g a str<strong>in</strong>g at rest an <strong>in</strong>itial<br />

Figure 2.3.<br />

displacement isf(x) = 0, and the <strong>in</strong>itial velocity isg(x) = 1/(1 + 0.25x<br />

2). Plot<br />

2. Calculate the exact solution to the Cauchy problem when c = 2, the <strong>in</strong>itial<br />

1. Derive d’Alembert’s formula (2.14) by determ<strong>in</strong><strong>in</strong>g the two arbitrary func<br />

So the signal breaks <strong>in</strong>to two pieces, and they propagate <strong>in</strong> opposite direc<br />

So the signal stays at the same place, but it spreads out and decreases <strong>in</strong><br />

ill + 4kt<br />

u(x, t) = e2’(1+41t)<br />

1<br />

via<br />

u(x, t) =<br />

u(x, t) = O5(e_(_Ct)<br />

2 + ex4ct)<br />

2). Th<strong>in</strong>k of this signal as be<strong>in</strong>g a bit of <strong>in</strong>formation.<br />

case. For example, suppose the <strong>in</strong>itial signal is a Gaussian function or bellshaped<br />

curve exp(—x<br />

gate <strong>in</strong>formation is to determ<strong>in</strong>e how a signal is propagated <strong>in</strong> a special<br />

is reta<strong>in</strong>ed. Parabolic, or diffusion, equations propagate signals at <strong>in</strong>f<strong>in</strong>ite<br />

herency <strong>in</strong> the wave form as it propagates, and therefore <strong>in</strong>formation<br />

parabolic and hyperbolic problems. Hyperbolic, or wave-like, equations<br />

The convection equation Ut + CU = 0, which is a wave-like equation,<br />

amplitude. Any <strong>in</strong>formation <strong>in</strong> the signal is eventually lost.<br />

:ion at (xo, to). This<br />

rm u 0. In hy<br />

F<strong>in</strong>ally, we po<strong>in</strong>t out aga<strong>in</strong> the important differences between<br />

speed; because the signals diffuse or smear out, there is a gradual loss of<br />

tions at speed c. The diffusion equation Ut = ku propagates the signal<br />

Section 1.2).<br />

coord<strong>in</strong>ate sytem where the problem simplifies (recall the examples <strong>in</strong><br />

propagate signals at a f<strong>in</strong>ite speed along characteristics; there is co<br />

<strong>in</strong>formation. A good way to understand how different equations propa<br />

propagates this signal via<br />

That is, it moves it at speed c without distortion. The wave equation<br />

=c2u moves it via<br />

2).<br />

)lution outside this<br />

the wave equation<br />

Exercises<br />

f<strong>in</strong>e a special coor<br />

ward <strong>in</strong> time to the<br />

regard<strong>in</strong>g the charsignals<br />

forward <strong>in</strong><br />

)und by trac<strong>in</strong>g the<br />

y the <strong>in</strong>itial values<br />

one. Look<strong>in</strong>g at the<br />

itive characteriSticS,<br />

then the region of<br />

constant are paths<br />

affect the solution<br />

—time along which<br />

ence of the <strong>in</strong>terval<br />

a and x —<br />

ct<br />

always zero outside<br />

:hat <strong>in</strong>terval, so the<br />

n surface.<br />

:h velocity c. These<br />

z.z. The<br />

Jnbounded Doma<strong>in</strong>s Exercises 57


ct)/p<br />

68<br />

2. Partial Differcntial Equations on Unbounded Doma<strong>in</strong>s<br />

4.<br />

where s(t) is a known signal. H<strong>in</strong>t: Look for a right-travel<strong>in</strong>g wave solution.<br />

The three-dimensional wave equation is<br />

5. Solve the Cauchy problem<br />

—c2Au = 0,<br />

where is = u(x, y, z, t) and A is the Laplacian operator. For waves with spher<br />

ical symmetry, is = u(p, ,<br />

where p = .Jx2 + y2 + z2. In this special case<br />

the Laplacian is given by (Section 1.8) Au = u + u. By <strong>in</strong>troduc<strong>in</strong>g a<br />

change of dependent variable U = pu, show that the general solution for the<br />

spherically symmetric wave equation<br />

is<br />

1<br />

u= —(F(p—ct)+G(p+cfl).<br />

p<br />

Why do you th<strong>in</strong>k an outward-mov<strong>in</strong>g wave is = F(p —<br />

amplitude?<br />

u1—c<br />

2u=0, xcR, t>0,<br />

2<br />

=c2(u + u)<br />

decays <strong>in</strong><br />

-i<br />

u(x, 0) =e1, Ut(x, 0) = cosx, x € R.<br />

Use a computer algebra program to graph the wave profile at t = 1, 2, 3. Take<br />

c= 1<br />

6. In Section 1.7 we showed that any solution to Laplace’s equation has the<br />

property that its value at a po<strong>in</strong>t is approximately the average of four nearby<br />

values surround<strong>in</strong>g the po<strong>in</strong>t. Can we make a statement about solutions to<br />

the wave equation? Consider any characteristic parallelogram (see Figure<br />

2.4) whose sides are positive and negative characteristics, and let A, B, C, D<br />

be the vertices as shown. Show that any solution to the wave equation satisfies<br />

the relation<br />

u(A) + u(C) = u(B) + u(D).<br />

Figure 2.4. Characteristic<br />

parallelogram.


oo.<br />

Therefore, for each t > 0,<br />

w(x, t) = f(()<br />

lfr(y))G(x<br />

-<br />

y,<br />

t)dy.<br />

x = 0 the solution is<br />

Exercises 71<br />

tisfies the Cauchy<br />

(x, t) are close. Let<br />

her. We would like<br />

i for Laplace’s equa<br />

I close <strong>in</strong> the sense<br />

b this end, consider<br />

easily observe that<br />

ther problems such<br />

)n the <strong>in</strong>itial and/or<br />

(ii) the solution is<br />

ions for various PDE<br />

m<strong>in</strong>g. We say that a<br />

an <strong>in</strong>itial boundary<br />

) solve the problem,<br />

eed only speci1y the<br />

m (2J5)—(2.16) does<br />

ty of cont<strong>in</strong>uous de<br />

e get the boundary<br />

<strong>in</strong> the boundary data<br />

lution should depend<br />

ehavior is disturb<strong>in</strong>g,<br />

er all, we want to be<br />

ation (2.15)—(2.16) is<br />

lation models steady<br />

:he data on the bound<br />

on Unbounded Doma<strong>in</strong>s<br />

this mean with regard to stability?<br />

showthatl u’(x,t)—u<br />

2(x,t) I 8 +8Tforallx c R, 0 < t < T.Whatdoes<br />

s<strong>in</strong>ce f G(x —<br />

f6G(x —y,t)dy = 6,<br />

Iu(x, t) —<br />

v(x,<br />

t)I f<br />

(y) — fr(y) II G(x<br />

—<br />

y,<br />

y,<br />

1. Show that the Cauchy problem for the backward diffusion equation,<br />

2. Let u = u(x, y). Is the problem<br />

uxy = 0, 0 < x,y < 1,<br />

3. Consider two Cauchy problems for the wave equation with different <strong>in</strong>itial<br />

u’(x, 0) = f’(x), u(x, 0) = g’(x), x<br />

I f’(x) —f<br />

2(x) &, I g’(x) —g<br />

2(x) IS 2,<br />

<strong>in</strong>dices and not exponents). If for all x E R we have<br />

for i = 1, 2, where f’ ,f2, g’, and g2 are given functions (the superscripts are<br />

= x E R, 0 < t < T,<br />

data:<br />

square, a well-posed problem? Discuss.<br />

on the unit square, where the value of u is prescribed on the boundary of the<br />

for large n.<br />

u(x, t) = 1 + _en2ts<strong>in</strong>nx<br />

is unstable by consider<strong>in</strong>g the solutions<br />

u(x, 0) = ft), x € R,<br />

u+u=0, xeRt>O,<br />

Exercises cJ<br />

Lc<br />

closeness of the <strong>in</strong>itial data implies closeness of the solution.<br />

t)dy = 1. Therefore, <strong>in</strong> the sense <strong>in</strong>terpreted above,<br />

t) I<br />

The solution formula (2.8) gives<br />


led Doma<strong>in</strong>s Exercises 75<br />

of an odd solution to the positive real axis is the solution to the given <strong>in</strong>i<br />

tial boundary value problem. If this <strong>in</strong>tuitive reason<strong>in</strong>g leaves the reader<br />

perplexed, then one can always verify analytically that the solutions we<br />

ace—time<br />

have obta<strong>in</strong>ed by this reflection method are, <strong>in</strong> fact, solutions to the given<br />

where the<br />

problems.<br />

<strong>2.2</strong>4) is If the boundary condition (2.18) along x = 0 <strong>in</strong> the heat flow problem<br />

ion x > Ct<br />

(2.17)—(219) is replaced by a Neumann condition<br />

y the<br />

g and can U(O, t) = 0, t > 0,<br />

embert’s<br />

then the problem can be solved by extend<strong>in</strong>g the <strong>in</strong>itial data to an<br />

even function. The same is true for the wave equation. We leave these<br />

calculations as exercises.<br />

(<strong>2.2</strong>7)<br />

0,<br />

(<strong>2.2</strong>8)<br />

Exercises<br />

1. Solve the problem<br />

< 0. = ku, x > 0, t > 0,<br />

u(0,t)=0, t >0,<br />

u(x, 0) = (x), x > 0,<br />

with an <strong>in</strong>sulated boundary condition by extend<strong>in</strong>g<br />

an even function. The solution is<br />

to all of the real axis as<br />

u(x, y) = f[G(x y, t) + G(x + y, ](y)dy.<br />

—<br />

s)ds<br />

0<br />

2. F<strong>in</strong>d a formula for the solution to the problem<br />

u=ku, x>0,t>0,<br />

u(0,t)=0, t >0,<br />

eplaced by<br />

may write u(x, 0) = 1, x > 0.<br />

(s)ds,<br />

Sketch a graph of several solution profiles with k = 0.5.<br />

3. F<strong>in</strong>d the solution to the problem<br />

(<strong>2.2</strong>9) 2<br />

Utt = C U, x > 0, t > 0,<br />

.24) is given u(0, t) = 0, t > 0,<br />

)rx 0.<br />

rn and wave<br />

hr these two Pick c = 0.5 and sketch several time snapshots of the solution surface to<br />

e restriction<br />

observe the reflection of the wave <strong>from</strong> the boundary.


—<br />

ku<br />

— =<br />

r)<br />

G(X<br />

r))f(y,<br />

80 2. Partial Differential Equations on Unbounded Doma<strong>in</strong>s 2.6. i<br />

Therefore, transformation of the dependent variable has changed the 2.6<br />

problem <strong>in</strong>to one with a homogeneous boundary condition, but a price<br />

was paid—an <strong>in</strong>homogeneity, or source term —g’(t), was <strong>in</strong>troduced <strong>in</strong>to<br />

Laplac<br />

the PDE. In general, we can always homogenize the boundary condi-<br />

equati<br />

tions <strong>in</strong> a l<strong>in</strong>ear problem, but the result is an <strong>in</strong>homogeneous PDE; SO<br />

ord<strong>in</strong>a<br />

<strong>in</strong>homogeneous boundary conditions can be traded for <strong>in</strong>homogeneous<br />

an alg<br />

PDEs.<br />

operat<br />

We can solve (2.49)—(<strong>2.5</strong>1) for v(x, t) by formulat<strong>in</strong>g a Duhamel’s pr<strong>in</strong>-<br />

Let<br />

ciple. However, <strong>in</strong> Section <strong>2.5</strong> we will observe that Laplace transform<br />

grow t<br />

methods can also be applied to f<strong>in</strong>d the solution.<br />

means<br />

the La;<br />

Exercises<br />

The La<br />

1. Write a formula<br />

given f<br />

for the solution to the problem<br />

U(s), ir<br />

s<strong>in</strong> x, x E R, t > o, variabi<br />

u(X, 0)<br />

where<br />

Ut(X, 0) = 0, X E R.<br />

is calle<br />

Graph the solution surface when c = 1.<br />

Laplac<br />

2. Write a formula for<br />

many<br />

the solution<br />

s<br />

to the problem<br />

shortt<br />

SrnX, XE R, t > 0, system<br />

u(x,0)=0,<br />

XER.<br />

3. Us<strong>in</strong>g Duhamel’s pr<strong>in</strong>ciple, f<strong>in</strong>d a formula for the solution to the <strong>in</strong>itial value returns<br />

problem for the convection equation<br />

Ut + cu = f(x, t), X E R, t > 0; u(X, 0) = 0, X E R.<br />

H<strong>in</strong>t: Look at the problem<br />

w(x,t;r)+cw(X,t;r)=O, XER,t >0; w(x,0;r)=flX,j, xeR.<br />

Solve the problem<br />

Ut + 2u = xet, X E R, t > 0; U(x, 0) = 0, x E R.<br />

which i<br />

use<br />

which r<br />

4. Formulate Duhamels pr<strong>in</strong>ciple and solve the <strong>in</strong>itial boundary value problem The<br />

that it c<br />

Ut = kU +f(x, t), X > 0, t > 0,<br />

transfoi<br />

U(X,0)=0, X >0,<br />

u(0,t)0, t >0.<br />

The solution is<br />

u(x, t) =<br />

f f (G(X y,<br />

—<br />

t<br />

—<br />

—<br />

+ y t<br />

—<br />

r)dydr.<br />

Formuh<br />

parts, ai<br />

equatlo;

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!