good PDE note 1
PDF · 279 pages · 1.0 MB
Open PDF file
A first-course textbook, not Phil's own writing, by Marcel B. Finan, dated August 2009. It reviews calculus and ODE background, then covers first-order PDEs and characteristics, the wave, heat and Laplace equations, Fourier series, separation of variables, and Laplace and Fourier transform methods. Chapters include practice problems and an answers section.
AI-written summary; may contain errors.
Extracted text (machine-read; may contain errors)
Undergraduate Notes in Mathematics
Arkansas Tech University
Department of Mathematics
A First Course of Partial Dierential
Equations
in Physical Sciences and Engineering
Marcel B. Finan
Arkansas Tech University
c
All Rights Reserved
2
Preface
Partial dierential equations are often used to construct models of the most
basic theories underlying physics and engineering. The goal of this book is to
develop the most basic ideas from the theory of partial dierential equations,
and apply them to the simplest models arising from the above mentioned
elds.
It is not easy to master the theory of partial dierential equations. Unlike
the theory of ordinary dierential equations, which relies on the fundamental
existence and uniqueness theorem, there is no single theorem which is central
to the subject. Instead, there are separate theories used for each of the major
types of partial dierential equations that commonly arise.
It is worth pointing out that the preponderance of dierential equations aris-
ing in applications, in science, in engineering, and within mathematics itself,
are of either rst or second order, with the latter being by far the most preva-
lent. We will mainly cover these two classes of PDEs.
This book is intended for a rst course in partial dierential equations at
the advanced undergraduate level for students in engineering and physical
sciences. It is assumed that the student has had the standard three semester
calculus sequence, and a course in ordinary dierential equations.
Marcel B Finan
August 2009
3
4 PREFACE
Contents
Preface 3
Preliminaries 7
1 Some Results of Calculus . . . . . . . . . . . . . . . . . . . . . . . 7
2 Sequences of Functions: Pointwise and Uniform Convergence . . . 14
Review of Some ODEs Results 23
3 The Method of Integrating Factor . . . . . . . . . . . . . . . . . . 23
4 The Method of Separation of Variables for ODEs . . . . . . . . . 28
5 Second Order Linear ODEs . . . . . . . . . . . . . . . . . . . . . 33
Introduction to PDEs 43
6 The Basic Concepts . . . . . . . . . . . . . . . . . . . . . . . . . 43
7 Solutions and Related Topics . . . . . . . . . . . . . . . . . . . . 52
First Order Partial Dierential Equations 65
8 Classication of First Order PDEs . . . . . . . . . . . . . . . . . 65
9 The One Dimensional Spatial Transport Equations . . . . . . . . 71
10 The Method of Characteristics . . . . . . . . . . . . . . . . . . . 79
11 The Cauchy Problem for First Order Quasilinear Equations . . . 86
Second Order Linear Partial Dierential Equations 99
12 Second Order PDEs in Two Variables . . . . . . . . . . . . . . . 99
13 Hyperbolic Type: The Wave equation . . . . . . . . . . . . . . . 104
14 Parabolic Type: The Heat Equation in One-Dimensional Space . 112
15 An Introduction to Fourier Series . . . . . . . . . . . . . . . . . 120
16 Fourier Sines Series and Fourier Cosines Series . . . . . . . . . . 132
17 Separation of Variables for PDEs . . . . . . . . . . . . . . . . . 139
5
6 CONTENTS
18 Solutions of the Heat Equation by the Separation of Variables
Method . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 146
19 Elliptic Type: Laplace's Equations in Rectangular Domains . . . 154
20 Laplace's Equations in Circular Regions . . . . . . . . . . . . . . 165
The Laplace Transform Solutions for PDEs 177
21 Essentials of the Laplace Transform . . . . . . . . . . . . . . . . 177
22 Solving PDEs Using Laplace Transform . . . . . . . . . . . . . . 192
The Fourier Transform Solutions for PDEs 199
23 Complex Version of Fourier Series . . . . . . . . . . . . . . . . . 199
24 The One Dimensional Fourier Transform . . . . . . . . . . . . . 205
25 Applications of Fourier Transforms to PDEs . . . . . . . . . . . 213
Answers and Solutions 221
Preliminaries
In this chapter we include some results from calculus which we will use often
in the study of partial dierential equations. Details and proof of these results
can be found in most calculus books.
1 Some Results of Calculus
The rst result provides a mean of showing when a function is zero on an
interval.
Theorem 1.1
(a) Suppose that fis continuous on an interval IRsuch thatRb
af(x)dx= 0
for all subintervals [ a;b]I:Thenf(x) = 0 for all x2I:
(b) Suppose that f: [a;b]!Ris continuous and non-negative. IfRb
af(x)dx=
0 thenf(x) = 0 on [a;b]:
(c) Suppose that f: [a;b]!Ris continuous such thatRb
af(x)g(x)dx= 0
for all continuous functions gon [a;b]:Thenf(x) = 0 on [a;b]:
Proof.
(a) Fixa2I:Letx2I:By the Fundamental Theorem of Calculus we have
0 =d
dxZx
af(t)dt=f(x):
Sincexwas arbitrary, we have f(x) = 0 for all x2I:
(b) Suppose the contrary. That is, suppose that x02[a;b] such that f(x0)>
0:By the continuity of f(x) atx0;there is a >0 such thatjx x0j<
impliesjf(x) f(x0)j<f(x0)
2:That is,jx x0j<impliesf(x)>f(x0)
2>0:
7
8 PRELIMINARIES
In words, there exists an open interval I[a;b] centered at x0such that
f(x)>0 for allx2I:Hence, because f(x)0 we must have
Zb
af(x)dxZ
If(x)dx> 0
which contradicts our assumption that the integral is zero. We conclude that
f(x) = 0 on [a;b]:
(c) This follows from (b) by taking g(x) =f(x)
Remark 1.1
The above theorem remains valid for functions in two variables. For example,
iff(x;y) is dened for xin an interval Iandyin an interval Jsuch that
Zb
aZd
cf(x;y)dxdy = 0
for all [a;b]Jand [c;d]Ithenf(x;y) = 0 over the rectangle IJ:
Example 1.1
Letf;g: [a;b]!Rbe continuous and such that f(x)g(x) for allxin
[a;b]:Show that ifRb
a(g(x) f(x))dx= 0 thenf(x)g(x) on [a;b]:
Solution.
Apply part (b) of previous theorem to the function h(x) =g(x) f(x)
Partial Derivatives
For multivariable functions, there are two common notations for partial
derivatives, and we shall employ them interchangeably. The rst is the Leib-
nitz notation that employs the symbol @to denote partial derivative. The
second, a more compact notation, is to use subscripts to indicate partial
derivatives. For example, utrepresents@u
@t;whileuxxrepresents@2u
@x2;anduxxt
becomes@3u
@2x@t:
An important formula of dierentiation is the so-called chain rule. If
u=u(x;y) wherex=x(s;t) andy=y(s;t) then
@u
@s=@u
@x@x
@s+@u
@y@y
@s:
1 SOME RESULTS OF CALCULUS 9
Likewise,
@u
@t=@u
@x@x
@t+@u
@y@y
@t:
Example 1.2
Compute the partial derivatives indicated:
(a)@
@y(y2sinxy)
(b)@2
@x2[ex+y]2
Solution.
(a) We have@
@y(y2sinxy) = sinxy@
@y(y2)+y2@
@y(sinxy) = 2ysinxy+xy2cosxy:
(b) We have@
@x[ex+y]2=@
@xe2(x+y)= 2e2(x+y):Thus,@2
@x2[ex+y]2=@
@x2e2(x+y)=
4e2(x+y)
Example 1.3
Supposeu(x;y) = sin (x2+y2);wherex=tesandy=s+t:Findusandut:
Solution.
We have
us=uxxs+uyys= 2xcos (x2+y2)tes+ 2ycos (x2+y2)
=[2tes+ 2(s+t)] cos [t2e2s+ (s+t)2]
Likewise,
ut=uxxt+uyyt= 2xcos (x2+y2)es+ 2ycos (x2+y2)
=[2tes+ 2(s+t)] cos [t2e2s+ (s+t)2]
Often we must dierentiate an integral with respect to a parameter which
may appear in the limits of integration, or in the integrand.
Letf(x;t) be a continuous function in the rectangle faxbgfct
dg:Assume that@f
@tis continuous on this rectangle. Dene the function
J(t) =Zb(t)
a(t)f(x;t)dx
wherea(t) andb(t) are continuously dierentiable functions of tsuch that
aa(t)b(t)b:Recall that a function f(x) is said to be continnu-
ously dierentiable if the derivative f0(x) exists, and is itself a continuous
function.
10 PRELIMINARIES
Theorem 1.2
dJ
dt=d
dtZb(t)
a(t)f(x;t)dx
=f(b(t);t)b0(t) f(a(t);t)a0(t) +Zb(t)
a(t)@f
@t(x;t)dx
Example 1.4
Consider the heat problem
ut=kuxx u; > 0; k> 0;0<x<L; t> 0
with boundary conditions ux(0;t) = 0 =ux(L;t) and initial condition u(x;0) =
f(x):LetE(t) =1
2RL
0u2dx:
(a) Show that E0(t)0:
(b) Show that E(t)RL
01
2jf(x)j2dx:
Solution.
(a) We have
dE
dt=1
2ZL
0@
@tu2(x;t)dx
=ZL
0u(x;t)ut(x;t)dx=kZL
0u(x;t)uxx(x;t)dx ZL
0u2(x;t)dx
=ku(x;t)ux(x;t)jL
0 kZL
0u2
x(x;t)dx ZL
0u2(x;t)dx
= kZL
0u2
x(x;t)dx ZL
0u2(x;t)dx0:
(b) From (a) we conclude that E(t) is a decreasing function of t>0:Thus,
E(t)E(0) =1
2ZL
0u2(x;0)dx=ZL
01
2jf(x)j2dx
The Least Upper Bound
A function f:D!Ris said to be bounded from above inDif there is a
constantMsuch thatf(x)Mfor allx2D:We callMan upper bound
1 SOME RESULTS OF CALCULUS 11
off:Note that the numbers, M+ 1;M+ 2;are also upper bounds of
f:The smallest upper bound of fis called the least upper bound or the
supremum. IfMis the supremum of finDwe write
M= supff(x) :x2Dg:
Note that if Nis any upper bound of finDthenMN:
Example 1.5
Find the supremum of f(x) = sinx:
Solution.
The graph of fis bounded between 1 and 1:Thus, supff(x) :x2Rg= 1
Example 1.6
Find
sup2sinx
sint
:x2R; t> 0
Solution.
The answer is
sup2sinx
sint
:x2R;t> 0
=2
12 PRELIMINARIES
Practice Problems
Exercise 1.1
Compute the partial derivatives indicated:
(a)@
@x(y2sinxy)
(b)@2
@x2(ex2y)
(c)@4
@x@y2@z
zln
x2
y
:
Exercise 1.2
Find all the rst partial derivatives of the functions:
(a)f(x;y) =x4+ 6py
(b)f(x;y;z ) =x2y 10y2z3+ 43x 7 tan (4y)
(c)f(s;t) =t7ln (s2) +9
t3 7p
s4
(d)f(x;y) = cos 4
x
ex2y 5y3
(e)f(u;v) =9u
u2+5v
(f)f(x;y;z ) =xsiny
z2
(g)f(x;y) =p
x2+ ln (5x 3y2)
Exercise 1.3
Letf(x;y) =e3xcosy:Computefx(0;2):
Exercise 1.4
Ifz=exsiny; x =st2;andy=s2t, nd@z
@sand@z
@t:
Exercise 1.5
In the equation
@u
@x @u
@y=x 2y
identify the independent variable(s) and the dependent variable.
Exercise 1.6
Letfbe an odd function, that is, f( x) = f(x) for allx2R:Show that
for alla2Rwe haveZa
af(x)dx= 0:
1 SOME RESULTS OF CALCULUS 13
Exercise 1.7
Letfbe an even function, that is, f( x) =f(x) for allx2R:Show that
for alla2Rwe have
Za
af(x)dx= 2Za
0f(x)dx:
Exercise 1.8
Use the product rule of derivatives to derive the formula of integration by
parts Z
uv0dx=uv Z
u0vdx:
Exercise 1.9
Letu(x;t) =2sin x
sin t
:Finduttanduxx:
Exercise 1.10
Letu(x;t) =2sin x
sinh t
;where
sinhx=ex e x
2:
Finduttanduxx:
Exercise 1.11
Find
sup2sinht
sinx
:x2R
:
Exercise 1.12
Letun(x;t) = 1 +en2t
nsinnx.
(a) Find supfjun(x;0) 1j:x2Rg:
(b) Find supfjun(x;t) 1j:x2Rg:
14 PRELIMINARIES
2 Sequences of Functions: Pointwise and Uni-
form Convergence
Later in this book we will be constructing solutions to PDEs involving innite
sums of sines and cosines. These innite sums or series are called Fourier
series. Fourier series are examples of series of functions. Convergence of
series of functions is dened in terms of convergence of a sequence of func-
tions. In this section we study the two types of convergence of sequences of
functions.
Recall that a sequence of numbers fang1
n=1is said to converge to a number
Lif and only if for every given >0 there is a positive integer N=N()
such that for all nNwe havejan Lj<:
What is the analogue concept of convergence when the terms of the sequence
are variables? Let DRand for each n2Nconsider a function fn:D!R:
Thus, we obtain a sequence of functions ffng1
n=1:For such a sequence, there
are two types of convergenve that we consider in this section: pointwise con-
vergence and uniform convergence.
We say thatffng1
n=1converges pointwise onDto a function f:D!Rif
and only if for a given a2Dand>0 there is a positive integer N=N(a;)
such that if nNthenjfn(a) f(a)j<:In symbol, we write
lim
n!1fn(a) =f(a):
It is important to note that Nis a function of both aand:
Example 2.1
Denefn: [0;1)!Rbyfn(x) =nx
1+n2x2:Show that the sequence ffng1
n=1
converges pointwise to the function f(x) = 0 for all x0:
Solution.
For allx0;
lim
n!1fn(x) = lim
n!1nx
1 +n2x2= 0
Example 2.2
For each positive integer nletfn: (0;1)!Rbe given by fn(x) =nx:Show
thatffng1
n=1does not converge pointwise on D= (0;1):
2 SEQUENCES OF FUNCTIONS: POINTWISE AND UNIFORM CONVERGENCE 15
Solution.
This follows from the fact that lim n!1nx=1for allx2D
As pointed out above, for pointwise convergence, the positive integer Nde-
pends on both the given xand:A stronger convergence concept can be
dened where Ndepends only on :
LetDbe a subset of Rand letffng1
n=1be a sequence of functions dened on
D:We say thatffng1
n=1converges uniformly onDto a function f:D!R
if and only if for all >0 there is a positive integer N=N() such that if
nNthenjfn(x) f(x)j<for allx2D:
This denition says that the integer Ndepends only on the given so that
fornN, the graph of fn(x) is bounded above by the graph of f(x) +and
below by the graph of f(x) :
Example 2.3
For each positive integer nletfn: [0;1]!Rbe given by fn(x) =x
n:Show
thatffng1
n=1converges uniformly to the zero function.
Solution.
Let >0 be given. Let Nbe a positive integer such that N >1
:Then for
nNwe have
jfn(x) f(x)j=jxj
n1
n1
N<
for allx2[0;1]
Clearly, uniform convergence implies pointwise convergence to the same limit
function. However, the converse is not true in general.
Example 2.4
Denefn: [0;1)!Rbyfn(x) =nx
1+n2x2:By Example 2.1, this sequence
converges pointwise to f(x) = 0:Let=1
3:Show that there is no positive
integerNwith the property nNimpliesjfn(x) f(x)j<for allx0:
Hence, the given sequence does not converge uniformly to f(x):
Solution.
For any positive integer Nand fornNwe have
fn1
n
f1
n=1
2>
16 PRELIMINARIES
Exercise 2.1 below shows a sequence of continuous functions converging point-
wise to a discontinuous function. That is, pointwise convergence does not
preserve the property of continuity. One of the interesting features of uniform
convergence is that it preserves continuity as shown in the next example.
Example 2.5
Suppose that for each n1 the function fn:D!Ris continuous in D:
Suppose thatffng1
n=1converges uniformly to f:Leta2D:
(a) Let>0 be given. Show that there is a positive integer Nsuch that if
nNthenjfn(x) f(x)j<
3for allx2D:
(b) Show that there is a >0 such that for all jx aj<we havejfN(x)
fN(a)j<
3:
(c) Using (a) and (b) show that for jx aj< we havejf(x) f(a)j< :
Hence,fis continuous in Dsinceawas arbitrary. Symbolically we write
lim
x!alim
n!1fn(x) = lim
n!1lim
x!afn(x):
Solution.
(a) This follows from the denition of uniform convergence.
(b) This follows from the fact that fNis continuous at a2D:
(c) Forjx aj<we havejf(x) f(a)j=jf(a) fN(a) +fN(a) fN(x) +
fN(x) f(x)jjfN(a) f(a)j+jfN(a) fN(x)j+jfN(x) f(x)j<
3+
3+
3=
Does pointwise convergenvce preserve integration? In real analysis, it is
proven that pointwise convergence does not preserve integrability. That is,
the pointwise limit of a sequence of integrable functions need not be inte-
grable. Even when a sequence of functions converges pointwise, the process
of interchanging limits and integration is not true in general.
Contrary to pointwise convergence, uniform convergence preserves integra-
tion. Moreover, limits and integration can be interchanged. That is, if
ffng1
n=1converges uniformly to fon a closed interval [ a;b] then
lim
n!1Zb
afn(x)dx=Zb
alim
n!1fn(x)dx:
Now, what about dierentiablility? Again, pointwise convergence fails in
general to conserve the dierentiability property. See Exercise 2.1. Does
uniform convergence preserve dierentiability? The answer is still no as
shown in the next example.
2 SEQUENCES OF FUNCTIONS: POINTWISE AND UNIFORM CONVERGENCE 17
Example 2.6
Consider the family of functions fn: [ 1;1] given by fn(x) =q
x2+1
n:
(a) Show that fnis dierentiable for each n1:
(b) Show that for all x2[ 1;1] we have
jfn(x) f(x)j1pn
wheref(x) =jxj:Hint: Note thatq
x2+1
n+p
x21pn:
(c) Let>0 be given. Show that there is a positive integer Nsuch that for
nNwe have
jfn(x) f(x)j<for allx2[ 1;1]:
Thus,ffng1
n=1converges uniformly to the non-dierentiable function f(x) =
jxj:
Solution.
(a)fnis the composition of two dierentiable functions so it is dierentiable
with derivative
f0
n(x) =x
x2+1
n 1
2
:
(b) We have
jfn(x) f(x)j=r
x2+1
n p
x2=(q
x2+1
n p
x2)(q
x2+1
n+p
x2)
q
x2+1
n+p
x2
=1
nq
x2+1
n+p
x2
1
n
1pn=1pn
(c) Let >0 be given. Since lim n!11pn= 0 we can nd a positive integer
Nsuch that for all nNwe have1pn<: Now the answer to the question
follows from this and part (b)
Even when uniform convergence occurs, the process of interchanging lim-
its and dierentiation may fail as shown in the next example.
18 PRELIMINARIES
Example 2.7
Consider the functions fn:R!Rdened byfn(x) =sinnx
n:
(a) Show thatffng1
n=1converges uniformly to the function f(x) = 0:
(b) Note thatffng1
n=1andfare dierentiable functions. Show that
lim
n!1f0
n(x)6=f0(x) =h
lim
n!1fn(x)i0
:
That is, one cannot, in general, interchange limits and derivatives.
Solution.
(a) Let>0 be given. Let Nbe a positive integer such that N >1
:Then
fornNwe have
jfn(x) f(x)j=sinnx
n1
n<
and this is true for all x2R:Hence,ffng1
n=1converges uniformly to the
functionf(x) = 0:
(b) We have lim n!1f0
n() = limn!1cosn= limn!1( 1)nwhich does not
converge. However, f0() = 0
Pointwise convergence was not enough to preserve dierentiability, and nei-
ther was uniform convergence by itself. Even with uniform convergence the
process of interchanging limits with derivatives is not true in general. How-
ever, if we combine pointwise convergence with uniform convergence we can
indeed preserve dierentiability and also switch the limit process with the
process of dierentiation.
Theorem 2.3
Letffng1
n=1be a sequence of dierentiable functions on [ a;b] that converges
pointwise to some function fdened on [ a;b]:Ifff0
ng1
n=1converges uniformly
on [a;b] to a function g;then the function fis dierentiable with derivative
equals tog:Thus,
lim
n!1f0
n(x) =g(x) =f0(x) =h
lim
n!1fn(x)i0
:
Finally, we conclude this section with the following important result that is
useful when a given sequence is bounded.
2 SEQUENCES OF FUNCTIONS: POINTWISE AND UNIFORM CONVERGENCE 19
Theorem 2.4
Consider a sequence fn:D!R:Then this sequence converges uniformly to
f:D!Rif and only if
lim
n!1supfjfn(x) f(x)j:x2Dg= 0:
Example 2.8
Show that the sequence dened by fn(x) =cosx
nconverges uniformly to the
zero function.
Solution.
We have
0supfjcosx
nj:x2Rg1
n:
Now apply the squeeze rule for sequences we nd that
lim
n!1supfjcosx
nj:x2Rg= 0
which implies that the given sequence converges uniformly to the zero func-
tion on R
20 PRELIMINARIES
Practice Problems
Exercise 2.1
Denefn: [0;1]!Rbyfn(x) =xn:Denef: [0;1]!Rby
f(x) =0 if 0x<1
1 ifx= 1:
(a) Show that the sequence ffng1
n=1converges pointwise to f:
(b) Show that the sequence ffng1
n=1does not converge uniformly to f:Hint:
Suppose otherwise. Let = 0:5 and get a contradiction by using a point
(0:5)1
N<x< 1:
Exercise 2.2
Consider the sequence of functions
fn(x) =nx+x2
n2
dened for all xinR:Show that this sequence converges pointwise to a
functionfto be determined.
Exercise 2.3
Consider the sequence of functions
fn(x) =sin (nx+ 3)pn+ 1
dened for all xinR:Show that this sequence converges pointwise to a
functionfto be determined.
Exercise 2.4
Consider the sequence of functions dened by fn(x) =n2xnfor all 0x1:
Show that this sequence does not converge pointwise to any function.
Exercise 2.5
Consider the sequence of functions dened by fn(x) = (cosx)nfor all
2
x
2:Show that this sequence converges pointwise to a noncontinuous
function to be determined.
2 SEQUENCES OF FUNCTIONS: POINTWISE AND UNIFORM CONVERGENCE 21
Exercise 2.6
Consider the sequence of functions fn(x) =x xn
ndened on [0 ;1):
(a) Doesffng1
n=1converge to some limit function? If so, nd the limit func-
tion and show whether the convergence is pointwise or uniform.
(b) Doesff0
ng1
n=1converge to some limit function? If so, nd the limit func-
tion and show whether the convergence is pointwise or uniform.
Exercise 2.7
Letfn(x) =xn
1+xnforx2[0;2]:
(a) Find the pointwise limit f(x) = limn!1fn(x) on [0;2]:
(b) Doesfn!funiformly on [0 ;2]?
Exercise 2.8
For eachn2Ndenefn:R!Rbyfn(x) =n+cosx
2n+sin2x:
(a) Show that fn!1
2uniformly.
(b) Find lim n!1R7
2fn(x)dx:
Exercise 2.9
Show that the sequence dened by fn(x) = (cosx)ndoes not converge uni-
formly on [
2;
2]:
Exercise 2.10
Letffng1
n=1be a sequence of functions such that
supfjfn(x)j: 2x5g2n
1 + 4n:
(a) Show that this sequence converges uniformly to a function fto be found.
(b) What is the value of the limit lim n!1R5
2fn(x)dx?
22 PRELIMINARIES
Review of Some ODEs Results
Later on in this book, we will encounter problems where a given partial
dierential is reduced to an ordinary dierential function by means of a given
change of variables. Then techniques from the theory of ODE are required in
solving the transformed ODE. In this chapter, we include some of the results
from ODE theory that will be needed in our future discussions.
3 The Method of Integrating Factor
In this section, we discuss a technique for solving the rst order linear non-
homogeneous equation
y0+p(t)y=g(t) (3.1)
wherep(t) andg(t) are continuous on the open interval a<t<b:
Sincep(t) is continuous, it has an antiderivative namelyR
p(t)dt:Let(t) =
eR
p(t)dt:Multiply Equation (3.1) by (t) and notice that the left hand side of
the resulting equation is the derivative of a product. Indeed,
d
dt((t)y) =(t)g(t):
Integrate both sides of the last equation with respect to tto obtain
(t)y=Z
(t)g(t)dt+C
Hence,
y(t) =1
(t)Z
(t)g(t)dt+C
(t)
or
y(t) =e R
p(t)dtZ
eR
p(t)dtg(t)dt+Ce R
p(t)dt
23
24 REVIEW OF SOME ODES RESULTS
Notice that the second term of the previous expression is just the general
solution for the homogeneous equation
y0+p(t)y= 0
whereas the rst term is a solution to the nonhomogeneous equation. That
is, the general solution to Equation (3.1) is the sum of a particular solution of
the nonhomogeneous equation and the general solution of the homogeneous
equation.
Example 3.1
Solve the initial value problem
y0 y
t= 4t; y(1) = 5:
Solution.
We havep(t) = 1
tso that(t) =1
t:Multiplying the given equation by the
integrating factor and using the product rule we notice that
1
ty0
= 4:
Integrating with respect to tand then solving for ywe nd that the general
solution is given by
y(t) =tZ
4dt+Ct= 4t2+Ct:
Sincey(1) = 5;we ndC= 1 and hence the unique solution to the IVP is
y(t) = 4t2+t;0<t<1
Example 3.2
Find the general solution to the equation
y0+2
ty= lnt; t> 0:
Solution.
The integrating factor is (t) =eR2
tdt=t2:Multiplying the given equation
byt2to obtain
(t2y)0=t2lnt:
3 THE METHOD OF INTEGRATING FACTOR 25
Integrating with respect to twe nd
t2y=Z
t2lntdt+C:
The integral on the right-hand side is evaluated using integration by parts
withu= lnt;dv =t2dt;du =dt
t;v=t3
3obtaining
t2y=t3
3lnt t3
9+C
Thus,
y=t
3lnt t
9+C
t2
26 REVIEW OF SOME ODES RESULTS
Practice Problems
Exercise 3.1
Solve the IVP: y0+ 2ty=t; y(0) = 0:
Exercise 3.2
Find the general solution: y0+ 3y=t+e 2t:
Exercise 3.3
Find the general solution: y0+1
ty= 3 cost; t> 0:
Exercise 3.4
Find the general solution: y0+ 2y= cos (3t):
Exercise 3.5
Find the general solution: y0+ (cost)y= 3 cost:
Exercise 3.6
Given that the solution to the IVP ty0+ 4y=t2; y(1) = 1
3exists on the
interval 1<t<1:What is the value of the constant ?
Exercise 3.7
Suppose that y(t) =Ce 2t+t+ 1 is the general solution to the equation
y0+p(t)y=g(t):Determine the functions p(t) andg(t):
Exercise 3.8
Suppose that y(t) = 2e t+et+ sintis the unique solution to the IVP
y0+y=g(t); y(0) =y0:Determine the constant y0and the function g(t):
Exercise 3.9
Find the value (if any) of the unique solution to the IVP y0+ (1 + cost)y=
1 + cost; y(0) = 3 in the long run?
Exercise 3.10
Solve
aux+buy+cu= 0
by using the change of variables s=ax+byandt=bx ay:
3 THE METHOD OF INTEGRATING FACTOR 27
Sample Exam Questions
Exercise 3.11
Solve the initial value problem ty0=y+t; y(1) = 7:
Exercise 3.12
Show that if aandare positive constants, and bis any real number, then
every solution of the equation
y0+ay=be t
has the property that y!0 ast!1 . Hint: Consider the cases a=and
a6=separately.
Exercise 3.13
Solve the initial-value problem y0+y=ety2;y(0) = 1 using the substitution
u(t) =1
y(t)
Exercise 3.14
Solve the initial-value problem ty0+ 2y=t2 t+ 1; y(1) =1
2
Exercise 3.15
Solvey0 1
ty= sint; y (1) = 3:Express your answer in terms of the sine
integral ,Si(t) =Rt
0sins
sds:
28 REVIEW OF SOME ODES RESULTS
4 The Method of Separation of Variables for
ODEs
The method of separation of variables that you have seen in the theory of
ordinary dierential equations has an analogue in the theory of partial dif-
ferential equations (Section 17). In this section, we review the method for
ordinary dierentiable equations.
A rst order dierential equation is separable if it can be written with one
variable only on the left and the other variable only on the right:
f(y)y0=g(t)
To solve this equation, we proceed as follows. Let F(t) be an antiderivative
off(t) andG(t) be an antiderivative of g(t):Then by the Chain Rule
d
dtF(y) =dF
dydy
dt=f(y)y0
Thus,
f(y)y0 g(t) =d
dtF(y) d
dtG(t) =d
dt[F(y) G(t)] = 0
It follows that
F(y) G(t) =C
which is equivalent to
Z
f(y)y0dt=Z
g(t)dt+C
As you can see, the result is generally an implicit equation involving a func-
tion ofyand a function of t:It may or may not be possible to solve this to
getyexplicitly as a function of t:For an initial value problem, substitute the
values oftandybyt0andy0to get the value of C:
Remark 4.2
IfFis a dierentiable function of yandyis a dierentiable function of tand
bothFandyare given then the chain rule allows us to nddF
dtgiven by
dF
dt=dF
dydy
dt
For separable equations, we are given f(y)y0=dF
dtand we are asked to nd
F(y):This process is referred to as \reversing the chain rule."
4 THE METHOD OF SEPARATION OF VARIABLES FOR ODES 29
Example 4.1
Solve the initial value problem y0= 6ty2; y(1) =1
25:
Solution.
Separating the variables and integrating both sides we obtain
Zy0
y2dt=Z
6tdt
or
Zd
dt1
y
dt=Z
6tdt
Thus,
1
y(t)= 3t2+C
Sincey(1) =1
25;we ndC= 28:The unique solution to the IVP is then
given explicitly by
y(t) =1
28 3t2
Example 4.2
Solve the IVP yy0= 4 sin (2t); y(0) = 1:
Solution.
This is a separable dierential equation. Integrating both sides we nd
Zd
dty2
2
dt= 4Z
sin (2t)dt
Thus,
y2= 4 cos (2t) +C
Sincey(0) = 1;we ndC= 5:Now, solving explicitly for y(t) we nd
y(t) =p
4 cost+ 5
Sincey(0) = 1;we havey(t) =p 4 cost+ 5:The interval of existence of
the solution is the interval 1<t<1
30 REVIEW OF SOME ODES RESULTS
Practice Problems
Exercise 4.1
Solve the (separable) dierential equation
y0=tet2 lny2:
Exercise 4.2
Solve the (separable) dierential equation
y0=t2y 4y
t+ 2:
Exercise 4.3
Solve the (separable) dierential equation
ty0= 2(y 4):
Exercise 4.4
Solve the (separable) dierential equation
y0= 2y(2 y):
Exercise 4.5
Solve the IVP
y0=4 sin (2t)
y; y(0) = 1:
Exercise 4.6
Solve the IVP:
yy0= sint; y(
2) = 2:
Exercise 4.7
Solve the IVP:
y0+y+ 1 = 0; y(1) = 0:
Exercise 4.8
Solve the IVP:
y0 ty3= 0; y(0) = 2:
4 THE METHOD OF SEPARATION OF VARIABLES FOR ODES 31
Exercise 4.9
Solve the IVP:
y0= 1 +y2; y(
4) = 1:
Exercise 4.10
Solve the IVP:
y0=t ty2; y(0) =1
2:
32 REVIEW OF SOME ODES RESULTS
Sample Exam Questions
Exercise 4.11
For what values of the constants ;y 0;and integer nis the function y(t) =
(4 +t) 1
2a solution of the initial value problem?
y0+yn= 0; y(0) =y0:
Exercise 4.12
Solve the equation 3 uy+uxy= 0 by using the substitution v=uy:
Exercise 4.13
Solve the IVP
(2y siny)y0= sint t; y(0) = 0:
Exercise 4.14
State an initial value problem, with initial condition imposed at t0= 2;
having implicit solution y3+t2+ siny= 4:
Exercise 4.15
Can the dierential equation
dy
dx=x2 xy
be solved by the method of separation of variables? Explain.
5 SECOND ORDER LINEAR ODES 33
5 Second Order Linear ODEs
When solving second order partial dierential equations such as the heat,
wave, and Laplace's equations using the method of separation of variables
for PDEs one ends up confronting second order linear ODEs. Thus, it is
deemed necessary to review some of the techniques used in solving second
order linear ordinary dierential equations which we do in this section.
We start rst by considering the second order linear ODE with constant
coecients given by
ay00+by0+cy= 0 (5.1)
wherea;bandcare constants with a6= 0:
Notice rst that for b= 0 andc6= 0 the function y00is a constant multiple
ofy:So it makes sense to look for a function with such property. One such
function is y(t) =ert:Substituting this function into (5.1) leads to
ay00+by0+cy=ar2ert+brert+cert= (ar2+br+c)ert= 0
Sinceert>0 for allt, the previous equation leads to
ar2+br+c= 0 (5.2)
Thus, a function y(t) =ertis a solution to (5.1) when rsatises equation
(5.2). We call (5.2) the characteristic equation for (5.1) and the polyno-
mialC(r) =ar2+br+cis called the characteristic polynomial .
The characteristic equation is a quadratic equation. Thus, this equation can
have two distinct real solutions, two equal solutions, or two conjugate com-
plex solutions depending on the sign of the expression b2 4ac:Hence, we
consider the following three cases:
Case 1:b2 4ac> 0:
In this case, equation (5.2) have two distinct real roots r1= b p
b2 4ac
4aand
r2= b+p
b2 4ac
4a:The general solution to (5.1) is given by
y(t) =c1er1t+c2er2t
wherec1andc2are arbitrary constants.
Example 5.1
Solve the initial value problem
y00 y0 6y= 0; y(0) = 1; y0(0) = 2:
Describe the behavior of the solution y(t) ast! 1 andt!1:
34 REVIEW OF SOME ODES RESULTS
Solution.
The characteristic polynomial is C(r) =r2 r 6 = (r 3)(r+2) so that the
characteristic equation r2 r 6 = 0 has the solutions r1= 3 andr2= 2:
The general solution is then given by
y(t) =c1e3t+c2e 2t:
Taking the derivative to obtain
y0(t) = 3c1e3t 2c2e 2t:
The conditions y(0) = 1 and y0(0) = 2 lead to the system
c1+c2= 1
3c1 2c2= 2:
Solving this system by the method of elimination we nd c1=4
5andc2=1
5:
Hence, the unique solution to the initial value problem is
y(t) =1
5(4e3t+e 2t):
Ast! 1; e3t!0 ande 2t!1:Thus,y(t)!1:Similarly,y(t)!1
ast!1
Case 2:b2 4ac= 0:
In this case, the characteristic equation has the single root r= b
2a:The
general solution to (5.1) is given by
y(t) =c1e b
2at+c2te b
2at
wherec1andc2are arbitrary constants.
Example 5.2
Solve the initial value problem: y00+ 2y0+y= 0; y(0) = 1; y0
1(0) = 1:
Solution.
The characteristic equation r2+2r+1 = 0 has a repeated root: r1=r2= 1:
Thus, the general solution is given by
y(t) =c1e t+c2te t:
5 SECOND ORDER LINEAR ODES 35
The two conditions y(0) = 1 and y0(0) = 1 lead toc1= 1 andc2= 0:
Hence, the unique solution is y(t) =e t
Case 3:b2 4ac< 0:
In this case, the complex roots of equation (5.1) are given by
r1;2= bip
4ac b2
2a
wherei=p 1:The general solution is given by
y(t) =et(c1cost+c2sint)
where= b
2a; =p
4ac b2
2a;andc1andc2are real numbers.
Example 5.3
Solve the initial value problem
y00 10y0+ 29y= 0; y(0) = 1; y0(0) = 3:
Solution.
The characteristic equation r2 10r+ 29 = 0 has the complex roots r1;2=
52i:Thus, the general solution is given by the expression
y(t) =e5t(c1cos 2t+c2sin 2t):
Findingy0we obtain
y0(t) =e5t[(5c1+ 2c2) cos 2t+ (5c2 2c1) sin 2t]:
The initial conditions yield c1= 1 andc2= 1:Thus, the unique solution
to the initial value problem is
y(t) =e5t(cos 2t sin 2t)
An Eigenvalue Problem
Consider the question of nding a nontrivial twice dierentiable function u
satisfying the ordinary dierential equation
d2u
dx2=u;0<x< 1:
36 REVIEW OF SOME ODES RESULTS
subject to the boundary conditions u(0) =u(1) = 0:This problem is referred
to as the eigenvalue problem for the following reason: Dene the function
Ld2
dx2:Then the given equation can be written as Lu=u:In linear
algebra,is called an eigenvalue ofLwith corresponding eigenvector u:
Dierent solutions to the eigenvalue problem are obtained depending on the
sign of:Suppose rst that = 0:Thenu(x) =C1x+C2for arbitrary
constantsC1andC2:Using the boundary conditions we nd C1=C2= 0:
Hence,u0:
Suppose that > 0:Thenu(x) =Aep
x+Be p
x:Again, the boundary
conditions imply that u0:
Now, suppose that < 0:Thenu(x) =Acosp
x+Bsinp
x:Using
the condition u(0) = 0 to obtain A= 0:Using the condition u(1) = 0
and assuming we are looking for non-trivial solution uwe expect to have
sinp
= 0:This happens when =n= (n)2wheren2N:We calln
an eigenvalue with corresponding eigenfunction un(x) = sinnx:
Finally, using the principle of superposition we nd that the general solution
to the eigenvalue problem is given by
u(x) =1X
n=1Ansinnx
where the convergence is pointwise convergence (See Section 2).
Euler Equations
A second order linear dierential equations of the form
ax2y00+bxy0+cy= 0
wherea;b;c are constants is called an Euler equation.
To solve Euler equation, one starts with solutions of the form y=xr(with
x > 0) whereris to be determined. Plugging this into the dierential
equation to get
ax2r(r 1)xr 2+bxrxr 1+cxr=0
(ar2 ar+br+c)xr=0
ar2 (a b)r+c=0
This last equation is a quadratic equation in rand so we will have three cases
to look at : Real distinct roots, double roots, and complex conjugate roots.
5 SECOND ORDER LINEAR ODES 37
If the quadratic equation has two distinct real roots r1andr2then the general
solution is given by
y(x) =Axr1+Bxr2:
If the quadratic equation has two equal roots r1=r2=rthen the general
solution is given by
y(x) =xr(A+Blnx):
If the quadratic equation has two complex conjugate solutions r1;2=i
then the general solution is given by
y(x) =x(Acos (lnx) +Bsin (lnx)):
Example 5.4
Solve the initial value problem
2x2y00+ 3xy0 15y= 0
y(1) = 0; y0(1) = 1:
Solution.
Lettingy=xrwe obtain the quadratic equation 2 r2+r 15 = 0 whose
roots arer1=5
2andr2= 3:Hence, the general solution is given by
y(x) =Ax5
2+Bx 3:
The condition y(1) = 0 implies A+B= 0:The condition y0(1) = 1 implies
5
2A 3B= 1:Solving this system of two unknowns we nd A=2
11and
B= 2
11:Hence, the unique solution is given by
y=2
11x5
2 2
11x 3
Second Order Linear nonhomogeneous ODE: The Method of Un-
determined Coecients
We consider the nonhomogeneous second order
ay00+by0+cy=g(t); a<t<b:
We know that the general solution has the structure
y(t) =c1y1(t) +c2y2(t) +yp(t)
38 REVIEW OF SOME ODES RESULTS
whereyp(t) is a particular solution to the nonhomogeneous equation. We
will writey(t) =yh(t) +yp(t) whereyh(t) =c1y1(t) +c2y2(t):
One way to nding ypis by using the method of undetermined coecients.
The idea behind the method of undetermined coecients is to look for yp(t)
which is of a form like that of g(t):This is possible only for special functions
g(t);but these special cases arise quite frequently in applications.
We will assume that g(t) being simple means it is some combination of terms
likeert;cos (kt);sin (kt);and polynomials antn+an 1tn 1+a1t+a0:Based
on those terms we will put together a candidate ypthat has some constants in
it we need to solve for: Those are the undetermined coecients this method
is named for.
In the following table we list examples of g(t) along with the corresponding
form of the particular solution.
Form ofg(t) Form ofyp(t)
antn+an 1tn 1++a1t+a0 tr[Antn+An 1tn 1++A1t+A0
[antn+an 1tn 1++a1t+a0]ettr[Antn+An 1tn 1++A1t+A0]et
[antn+an 1tn 1++a1t+a0] costtr[(Antn+An 1tn 1++A1t+A0) cost
or +(Bntn+Bn 1tn 1++B1t+B0) sint]
[antn+an 1tn 1++a1t+a0] sint
et[antn+an 1tn 1++a1t+a0] sinttr[(Antn+An 1tn 1++A1t+A0)etcost
or +(Bntn+Bn 1tn 1++B1t+B0)etsint]
et[antn+an 1tn 1++a1t+a0] cost
The number ris chosen to be the smallest nonnegative integer such that
no term in the assumed form is a solution of the homogeneous equation
ay00+by0+cy= 0:The value of rwill be 0, 1, or 2.
Example 5.5
List an appropriate form for a particular solution of
(a)y00+ 4y=t2e3t:
(b)y00+ 4y=te2tcost:
(c)y00+ 4y= 2t2+ 5 sin 2t+e3t:
(d)y00+ 4y=t2cos 2t:
Solution.
The general solution to the homogeneous equation is yh(t) =c1cos 2t+
c2sin 2t:
(a) Forg(t) =t2e3t, an appropriate particular solution has the form yp(t) =
5 SECOND ORDER LINEAR ODES 39
tr(A2t2+A1t+A0)e3t:We taker= 0 since no term in the assumed form for
ypis present in the expression of yh(t):Thus
yp(t) = (A2t2+A1t+A0)e3t
(b) An appropriate form is
yp(t) =tr[(A1t+A0)e2tcost+ (B1t+B0)e2tsint]
We taker= 0 since no term in the assumed form for ypis present in the
expression of yh(t):Thus
yp(t) = (A1t+A0)e2tcost+ (B1t+B0)e2tsint
(c)
yp(t) =A2t2+A1t+A0+B0tcos 2t+C0tsin 2t+D0e3t
(d)
yp(t) =t(A2t2+A1t+A0) cos 2t+t(B2t2+B1t+B0) sin 2t
Example 5.6
Find the general solution of
y00 2y0 3y= 4t 5 + 6te2t
Solution.
The characteristic equation of the homogeneous equation is r2 2r 3 = 0
with rootsr1= 1 andr2= 3:Thus,
yh(t) =c1e t+c2e3t
A guess for the particular solution is yp(t) =At+B+Cte2t+De2t:Inserting
this into the dierential equation leads to
3At 2A 3B 3Cte2t+ (2C 3D)e2t= 4t 5 + 6te2t
From this identity we obtain 3A= 4 so that A= 4
3:Also, 2A 3B= 5
so thatB=23
9:Since 3C= 6 we nd C= 2:From 2C 3D= 0 we nd
D= 4
3:It follows that
y(t) =c1e t+c2e3t 4
3t+23
9
2t+4
3
e2t
40 REVIEW OF SOME ODES RESULTS
Practice Problems
Exercise 5.1
Solve the initial value problem
y00 4y0+ 3y= 0; y(0) = 1; y0(0) = 1
Describe the behavior of the solution y(t) ast! 1 andt!1:
Exercise 5.2
Solve the initial value problem
y00+ 4y0+ 2y= 0; y(0) = 0; y0(0) = 4
Describe the behavior of the solution y(t) ast! 1 andt!1:
Exercise 5.3
Solve the initial value problem
2y00 y= 0; y(0) = 2; y0(0) =p
2
Describe the behavior of the solution y(t) ast! 1 andt!1:
Exercise 5.4
Find a homogeneous second-order linear ordinary dierential equation whose
general solution is y(t) =c1e2t+c2e t:
Exercise 5.5
Solve the IVP
9y00 6y0+y= 0; y(3) = 2; y0(3) = 5
3
Exercise 5.6
Solve the IVP
25y00+ 20y0+ 4y= 0; y(5) = 4e 2; y0(5) = 3
5e 2
Exercise 5.7
The graph of a solution y(t) of the dierential equation 4 y00+ 4y0+y= 0
passes through the points (1 ;e 1
2) and (2;0):Determiney(0) andy0(0):
5 SECOND ORDER LINEAR ODES 41
Exercise 5.8
Find the general solution of y00 6y0+ 9y= 0:
Exercise 5.9
Solve the IVP
y00+ 2y0+ 2y= 0; y(0) = 3; y0(0) = 1
Exercise 5.10
Solve the IVP
2y00 2y0+y= 0; y( ) = 1; y0( ) = 1
Exercise 5.11
Find the general solution of
y00 y0+y= 2 sin 3t
Exercise 5.12
Find the general solution of
y00+ 4y0 2y= 2t2 3t+ 6
42 REVIEW OF SOME ODES RESULTS
Sample Exam Questions
Exercise 5.13
Find the general solution to the following dierential equation.
x2y00 7xy0+ 16y= 0:
Exercise 5.14
Find the general solution to the following dierential equation.
x2y00+ 3xy0+ 4y= 0:
Exercise 5.15
Consider the dierential equation
d2y
dx2+y= 0:
Determine the eigenvalues and the corresponding eigenfunctions if ysatis-
es the following boundary conditions:
(a)y(0) =y() = 0
(b)y(0) =y0(L) = 0
(c)y0(0) =y(1) = 0:
Exercise 5.16
Show by direct computation that the eigenvalue problems
(ky0(x))0+y(x) = 0; k> 0
with the following boundary conditions have no negative eigenvalues :
(a)y(0) =y(L) = 0
(b)y0(0) =y0(L) = 0
(c)y(L) =y( L); y0(L) =y0( L):
Exercise 5.17
Solve the initial-value problem: 2 y00+ 5y0 3y= 0; y(0) = 2; y0(0) = 1:
Exercise 5.18
Find the general solution of
y00 y0= 5et sin 2t
Exercise 5.19
Solve using undetermined coecients:
y00+y0 2y=t+ sin 2t;y(0) = 1;y0(0) = 0
Introduction to PDEs
Many elds in engineering and the physical sciences require the study of ODE
and PDE. Examples of those elds are acoustics, aerodynamics, elasticity,
electrodynamics,
uid dynamics, geophysics (seismic wave propagation), heat
transfer, meteorology, oceanography, optics, petroleum engineering, plasma
physics (ionized liquids and gases), quantum mechanics.
So the study of partial dierential equation is of great importance to the
above mentioned elds. The purpose of this chapter is to introduce the
reader to the basic terms of partial dierential equations.
6 The Basic Concepts
The goal of this section is to introduce the reader to the basic concepts and
notations that will be used in the remainder of this book.
Adierential equation is an equation that involves an unknown scalar
function (the dependent variable) and one or more of its derivatives. For
example,
d2y
dx2 5dy
dx+ 3y= 3 (6.1)
or
@u
@t @2u
@x2 @2u
@y2+u= 0: (6.2)
If the unknown function is a function in one single variable then the dier-
ential equation is called an ordinary dierential equation. An example
of an ordinary dierential equation is Equation (6.1). In contrast, when the
unknown function is a function of two or more independent variables then
the dierential equation is called a partial dierential equation , in short
PDE. Equation (6.2) is an example of a partial dierential equation. In this
book we will be focusing on partial dierential equations.
43
44 INTRODUCTION TO PDES
Example 6.1
Identify which variables are dependent variable or independent variable(s)
for the following dierential equations.
(a)d4y
dx4 x2+y= 0
(b)utt+xutx= 0:
(c)xdx
dt= 4:
(d)@y
@u 4@y
@v=u+ 3y:
Solution.
(a) Independent variable is xand the dependent variable is y:
(b) Independent variables are xandtand the dependent variable is u:
(c) Independent variable is tand the dependent variable is x:
(d) Independent variables are uandvand the dependent variable is y
Example 6.2
Classify the following as either ODE or PDE.
(a)ut=c2uxx:
(b)y00 4y0+ 5y= 0:
(c)ut+cux= 5:
Solution.
(a) PDE (b) ODE (c) PDE
Theorder of a partial dierential equation is the highest order derivative
occurring in the equation. Thus, (6.2) is a second order partial dierential
equation.
Example 6.3
Find the order of each of the following partial dierential equations:
(a)xux+yuy=x2+y2
(b)uux+uy= 2
(c)utt c2uxx=f(x;t)
(d)ut+uux+uxxx= 0
(e)utt+uxxxx= 0:
Solution.
(a) First order (b) First order (c) Second order (d) Third order (e) Fourth
order
6 THE BASIC CONCEPTS 45
A partial dierential equation is called linear if it is linear in the unknown
function and all its derivatives with coecients depend only on the indepen-
dent variables. For example, a rst order linear partial dierential equation
has the form
A(x;y)ux+B(x;y)uy+C(x;y)u=D(x;y)
whereas a second order linear partial dierential equation has the form
A(x;y)uxx+B(x;y)uxy+C(x;y)uyy+D(x;y)ux+E(x;y)uy+F(x;y)u=G(x;y):
A partial dierential equation is called quasi-linear if the highest-order
derivatives which appear in the equation are of degree 1(regardless of the
manner in which lower-order derivatives and unknown functions occur in the
equation). For example, a rst order quasi-linear partial dierential equation
has the form
A(x;y;u )ux+B(x;y;u )uy=C(x;y;u )
whereas a second order quasi-linear partial dierential equation has the form
A(x;y;u;ux;uy)uxx+B(x;y;u;ux;uy)uxy+C(x;y;u;ux;uy)uyy=D(x;y;u;ux;uy):
A partial dierential equation is semi-linear if it is quasi-linear and the
coecients of the highest-order derivatives are functions of independent vari-
ables only. For example, a rst order semi-linear partial dierential equation
has the form
A(x;y)ux+B(x;y)uy=C(x;y;u )
whereas a second order semi-linear partial dierential equation has the form
A(x;y)uxx+B(x;y)uxy+C(x;y)uyy=D(x;y;u;u x;uy):
Note that linear and semi-linear partial dierential equations are special cases
of quasi-linear equations.
A partial dierential equation that is not linear is called nonlinear . For
example,u2
x+ 2uxy= 0:
As for ODEs, linear PDEs are usually simpler to analyze/solve than nonlinear
PDEs.
46 INTRODUCTION TO PDES
Example 6.4
Determine whether the given PDE is linear, quasilinear, semilinear, or non-
linear:
(a)xux+yuy=x2+y2
(b)uux+uy= 2
(c)utt c2uxx=f(x;t)
(d)ut+uux+uxxx= 0
(e)u2
tt+uxxxx= 0:
Solution.
(a) Linear, quasilinear, semilinear.
(b) Quasilinear, nonlinear.
(c) Linear, quasilinear, semilinear.
(d) Quasilinear, semilinear, nonlinear.
(e) Quasilinear, semilinear, nonlinear
A more precise denition of a linear dierential equation begins with the
concept of a linear dierential operator L:The operator Lis assembled
by summing the basic partial derivative operators, with coecients depend-
ing on the independent variables. The operator acts on suciently smooth
functions depending on the relevant independent variables. Linearity im-
poses two key requirements:
L[u+v] =L[u] +L[v] andL[u] =L[u];
for any two (suciently smooth) functions u; vand any constant :
Example 6.5
Dene a linear dierential operator for the PDE
ut=c2uxx:
Solution.
LetL[u] =ut c2uxx:Then one can easily check that L[u+v] =L[u] +L[v]
andL[u] =L[u]
A linear partial dierential equation is called homogeneous if every term
of the equation involves the unknown function or its partial derivatives. A
linear partial dierential equation that is not homogeneous is called nonho-
mogeneous. In this case, there is a term in the equation that involves only
6 THE BASIC CONCEPTS 47
the independent variables.
A homogeneous linear partial dierential equation has the form
L[u] = 0
whereLis a linear dierential operator.
Example 6.6
Determine whether the equation is homogeneous or nonhomogeneous:
(a)xux+yuy=x2+y2:
(b)utt=c2uxx:
(c)uxx+uyy= 0:
Solution.
(a) Nonhomogeneous because of x2+y2:
(b) Homogeneous.
(c) Homogeneous
Finally, we shall be employing a few basic notational conventions regard-
ing the variables that appear in our dierential equations. We always use
tto denote time, while x;y;z will represent (Cartesian) space coordinates.
Polar coordinates r;will also be used when needed, and our notational con-
ventions appear at the appropriate places in the exposition.
Anequilibrium equation models an unchanging physical system, and so
only involves the space variables. The time variable tappears when mod-
eling dynamical , meaning time-varying, processes. Both time and space
coordinates are independent variables.
48 INTRODUCTION TO PDES
Practice Problems
Exercise 6.1
Classify the following equations as either ODE or PDE.
(a) (y000)4+t2
(y0)2+4= 0
(b)@u
@x+y@u
@y=y x
y+x
(c)y00 4y= 0
Exercise 6.2
Write the equation
uxx+ 2uxy+uyy= 0
in the coordinates s=x; t=x y:
Exercise 6.3
Write the equation
uxx 2uxy+ 5uyy= 0
in the coordinates s=x+y; t= 2x:
Exercise 6.4
For each of the following PDEs, state its order and whether it is linear or
nonlinear. If it is linear, also state whether it is homogeneous or nonhomo-
geneous:
(a)uux+x2uyyy+ sinx= 0
(b)ux+ex2uy= 0
(c)utt+ (siny)uyy etcosy= 0:
Exercise 6.5
For each of the following PDEs, determine its order and whether it is linear
or not. For linear PDEs, state also whether the equation is homogeneous or
not. For nonlinear PDEs, circle all term(s) that are not linear.
(a)x2uxx+exu=xuxyy
(b)eyuxxx+exu= siny+ 10xuy
(c)y2uxx+exuux= 2xuy+u
(d)uxuxxy+exuuy= 5x2ux
(e)ut=k2(uxx+uyy) +f(x;y;t ):
6 THE BASIC CONCEPTS 49
Exercise 6.6
Which of the following PDEs are linear?
(a)Laplace's equation: uxx+uyy= 0:
(b)Convection (transport) equation: ut+cux= 0:
(c)Minimal surface equation: (1+Z2
y)Zxx 2ZxZyZxy+(1+Z2
x)Zyy= 0:
(d)Korteweg-Vries equation: ut+ 6uux=uxxx:
Exercise 6.7
Classify the following dierential equations as ODEs or PDEs, linear or
nonlinear, and determine their order. For the linear equations, determine
whether or not they are homogeneous.
(a) The diusion equation foru(x;t) :
ut=kuxx:
(b) The wave equation forw(x;t) :
wtt=c2wxx:
(c) The thin lm equation forh(x;t) :
ht= (hhxxx)x:
(d) The forced harmonic oscillator fory(t) :
ytt+!2y=Fcos (!t):
(e) The Poisson Equation for the electric potential ( x;y;z ) :
xx+ yy+ zz= 4(x;y;z ):
where(x;y;z ) is a known charge density.
(f)Burger's equation forh(x;t) :
ht+hhx=hxx:
Exercise 6.8
Write down the general form of a linear second order dierential equation of
a function in three variables.
50 INTRODUCTION TO PDES
Exercise 6.9
Give the orders of the following PDEs, and classify them as linear or nonlin-
ear. If the PDE is linear, specify whether it is homogeneous or nonhomoge-
neous.
(a)x2uxxy+y2uyy log (1 +y2)u= 0
(b)ux+u3= 1
(c)uxxyy+exux=y
(d)uuxx+uyy u= 0
(e)uxx+ut= 3u:
Exercise 6.10
Consider the second-order PDE
uxx+ 4uxy+ 4uyy= 0:
Use the change of variables v(x;y) =y 2xandw(x;y) =xto show that
uww= 0:
6 THE BASIC CONCEPTS 51
Sample Exam Questions
Exercise 6.11
Write the one dimensional wave equation utt=c2uxxin the coordinates
v=x+ctandw=x ct:
Exercise 6.12
Write the PDE
uxx+ 2uxy 3uyy= 0
in the coordinates v(x;y) =y 3xandw(x;y) =x+y:
Exercise 6.13
Write the PDE
aux+buy= 0
in the coordinates s(x;y) =ax+byandt(x;y) =bx ay:Assumea2+b2>0:
Exercise 6.14
Write the PDE
ux+uy= 1
in the coordinates s=x+yandt=x y:
Exercise 6.15
Write the PDE
aut+bux=u; a;b6= 0
in the coordinates v=ax btandw=1
at:
52 INTRODUCTION TO PDES
7 Solutions and Related Topics
By aclassical solution orstrong solution to a partial dierential equation
we mean a function that satises the equation. To solve a PDE is to nd all
its classical solutions. In the case of only two independent variables xandy;
a solutionu(x;y) is visualized geometrically as a surface, called a solution
surface or an integral surface in the (x;y;u ) space.
A formula that expresses all the solutions of a PDE is called the general
solution of the equation.
Example 7.1
Show that u(x;t) =e 22t(cosx sinx) is a solution to the equation
ut 2uxx= 0:
Solution.
Since
ut 2uxx= 22e 22t(cosx
sinx) 2e 22t( 2cosx+2sinx) = 0
the given function is a solution to the given equation
Example 7.2
Find the general solution of uxy= 0:
Solution.
Integrating rst we respect to ywe ndux(x;y) =f(x);wherefis an
arbitrary dierentiable function. Integrating uxwith respect to xwe nd
u(x;y) =R
f(x)dx+g(y);wheregis an arbitrary dierentiable function
Note that the general solution in the previous example involves two arbitrary
functions. In general, the general solution of a partial dierential equation
is an expression that involves arbitrary functions. This is in contrast to the
general solution of an ordinary dierential equation which involves arbitrary
constants.
Usually, a classical solution enjoys properties such as smootheness (i.e. a
function that has continuous derivatives up to some desired order over some
domain.) and continuity. However, in the theory of nonlinear pdes, there are
solutions that do not require the smoothness property. Such solutions are
7 SOLUTIONS AND RELATED TOPICS 53
called weak solutions orgeneralized solutions. We illustrate this con-
cept using equations rather than pdes. Consider the equation x2 y2= 0:
The function y=xis a classical solution of this equation. This solution is
innitely dierentiable function. On the other hand, the function y=jxjis
also a solution to the given equation. However, this solution is not dieren-
tiable at 0. We call such a solution a weak solution. In this book, the word
solution will refer to a classical solution.
Example 7.3
Show thatu(x;t) =t+1
2x2is a classical solution to the PDE
ut=uxx: (7.1)
Solution.
Assume that the domain of denition of uisDR2:Sinceu;ut;ux;utx;uxx
exist and are continuous in D(i.e.,uis smooth in D) andusatises equation
(7.1), we conclude that uis a classical solution to the given PDE
Now, consider the linear dierential operator Las dened in the previous
section. The dening properties of linearity immediately imply the key facts
concerning homogeneous linear (dierential) equations.
Theorem 7.1
The sum of two solutions to a homogeneous linear dierential equation is
again a solution, as is the product of a solution by any constant.
Proof.
Letu1;u2be solutions, meaning that L[u1] = 0 andL[u2] = 0:Then, thanks
to linearity,
L[u1+u2] =L[u1] +L[u2] = 0;
and hence their sum u1+u2is a solution. Similarly, if is any constant, and
uany solution, then
L[u] =L[u] =0 = 0;
and so the scalar multiple uis also a solution
The following result is known as the superposition principle for homo-
geneous linear equations.
54 INTRODUCTION TO PDES
Theorem 7.2
Ifu1;;unare solutions to a common homogeneous linear partial dieren-
tial equation L[u] = 0;then the linear combination u=c1u1++cnunis
a solution for any choice of constants c1;;cn:
Proof.
The key fact is that, thanks to the linearity of L, for any suciently smooth
functionsu1;;unand any constants c1;;cn;
L[u] =L[c1u1++cnun] =L[c1u1++cn 1un 1] +L[cnun]
==L[c1u1] ++L[cnun] =c1L[u1] ++cnL[un]:
In particular, if the functions are solutions, so L[u1] = 0;;L[un] = 0;then
the right hand side of the above equation vanishes, proving that uis also a
solution to the homogeneous equation L[u] = 0
In physical applications, homogeneous linear equations model unforced sys-
tems that are subject to their own internal constraints. External forcing
is represented by an additional term that does not involve the dependent
variable. This results in the nonhomogeneous equation
L[u] =f
whereLis a linear partial dierential operator, uis the dependent variable,
andfis a given non-zero function of the independent variables alone.
You already learned the basic philosophy for solving of nonhomogeneous
linear equations in your study of elementary ordinary dierential equations.
Step one is to determine the general solution to the homogeneous equation.
Step two is to nd a particular solution to the nonhomogeneous version. The
general solution to the nonhomogeneous equation is then obtained by adding
the two together. Here is the general version of this procedure:
Theorem 7.3
Letuibe a particular solution to the nonhomogeneous linear equation L[u] =
f:Then the general solution to L[u] =fis given by u=ui+uh, whereuhis
the general solution to the corresponding homogeneous equation L[u] = 0:
Proof.
Let us rst show that u=ui+uhis also a solution to L[u] =f:By linearity,
L[u] =L[ui+uh] =L[ui] +L[uh] =f+ 0 =f:
7 SOLUTIONS AND RELATED TOPICS 55
To show that every solution to the nonhomogeneous equation can be ex-
pressed in this manner, suppose usatisesL[u] =f:Setuh=u ui:Then,
by linearity,
L[uh] =L[u ui] =L[u] L[ui] = 0;
and henceuhis a solution to the homogeneous dierential equation. Thus,
u=ui+uhhas the required form
In physical applications, one can interpret the particular solution uias a
response of the system to the external forcing function, while the solution
uhto the homogeneous equation represents the system's internal, unforced
motion. The general solution to a linear nonhomogeneous equation is thus a
combination of the external and internal responses.
As you have noticed by now, one solution of a linear PDE leads to the cre-
ation of lots of solutions. In contrast, nonlinear equations are much tougher
to deal with, for example, knowledge of several solutions does not necessarily
help in constructing others. Indeed, even nding one solution to a nonlinear
partial dierential equation can be quite a challenge.
In this introductory course, we will primarily but not exclusively con-
centrate on analyzing the most basic linear partial dierential equations. But
we will have occasion to brie
y foray into the nonlinear realm, to appreciate
some recent developments in this fascinating area of contemporary research
and applications.
As observed above, a general solution of a partial dierential equation has
innitely many solutions. In almost all cases, this general solution is of little
use since it has to satisfy other supplementary conditions, usually called ini-
tial or boundary conditions. These conditions determine the unique solution
of interest.
Aboundary value problem is a partial dierential equation where either
the unknown function or its derivatives have values assigned on the physical
boundary of the domain in which the problem is specied. These conditions
are called boundary conditions . For example,
uxx+uyy=0 if 0 <x;y< 1
u(x;0) =u(x;1) =0 if 0 <x< 1
ux(0;y) =ux(1;y) =0 if 0 <y< 1:
56 INTRODUCTION TO PDES
There are three types of boundary conditions which arise frequently in for-
mulating physical problems:
1.Dirichlet Boundary Conditions: In this case, the dependent function
uis prescribed on the boundary of the bounded domain. For example, if the
bounded domain is the rectangular plate 0 < x < L 1and 0< y < L 2;the
boundary conditions u(0;y);u(L1;y);u(x;0);andu(x;L 2) are prescribed.
The boundary conditions are called homogeneous if the dependent variable
is zero at any point on the boundary, otherwise the boundary conditions are
called nonhomogeneous.
2.Neumann Boundary Conditions: In this case, rst partial derivatives
are prescribed on the boundary of the bounded domain. For example, the
Neuman boundary conditions for a rod of length L;where 0<x<L; are of
the formux(0;t) =andux(L;t) =;whereandare constants.
3.Robin or mixed Boundary Conditions: This occurs when the depen-
dent variable and its rst partial derivatives are prescribed on the boundary
of the bounded domain.
Aninitial valur problem (orCauchy problem ) is a partial dierential
equation together with a set of additional conditions on the solution or its
derivatives at either a given point or a given curve in the domain of the so-
lution. These conditions are called initial value conditions. For example,
thetransport equation
ut(x;t) +cux(x;t) =0
u(x;0) =f(x)
is a Cauchy problem.
It can be shown that initial conditions for a PDE are necessary and sucient
for the existence of a unique solution.
We say that an initial and/or boundary value problem associated with a PDE
iswell-posed if it has a solution which is unique and depends continuously
on the data given in the problem. The last condition, namely the continuous
dependence is important in physical problems. This condition means that
the solution changes by a small amount when the conditions change a little.
Such solutions are said to be stable .
Example 7.4
7 SOLUTIONS AND RELATED TOPICS 57
Forx2Randt>0 we consider the initial value problem
utt uxx=0
u(x;0) =ut(x;0) =0
Clearly,u(x;t) = 0 is a solution to this problem.
(a) Let 0<<< 1 be a very small number. Show that the function u(x;t) =
2sin x
sin t
is a solution to the problem
utt uxx=0
u(x;0) =0
ut(x;0) =sinx
(b) Show that sup fju(x;t) u(x;t)j:x2R;t > 0g=2:Thus, a small
change in the initial data leads to a small change in the solution. Hence, the
initial value problem is well-posed.
Solution.
(a) We have
@u
@t=sinx
cost
@2u
@t2= sinx
sint
@u
@x=cosx
sint
@2u
@x2= sinx
sint
Thus,@2u
@t2 @2u
@x2= 0:Moreover,u(x;0) = 0 and@
@tu(x;0) =sin x
:
(b) We have
supfju(x;t) u(x;t)j:x2R;t> 0g=2supfsinx
sint
:x2R;t> 0g
=2
A problem that is not well-posed is referred to as an ill-posed problem. We
illustrate this concept in the next example.
58 INTRODUCTION TO PDES
Example 7.5
Forx2Randt>0 we consider the initial value problem
utt+uxx=0
u(x;0) =ut(x;0) =0
Clearly,u(x;t) = 0 is a solution to this problem.
(a) Let 0<<< 1 be a very small number. Show that the function u(x;t) =
2sin x
sinh t
;where
sinhx=ex e x
2
is a solution to the problem
utt+uxx=0
u(x;0) =0
ut(x;0) =sinx
(b) Show that sup fj@
@tu(x;0) ut(x;0)j:x2Rg=and supfju(x;t)
u(x;t)j:x2Rg=2sinh t
:
(c) Find lim t!1supfju(x;t) u(x;t)j:x2Rg:
Solution.
(a) We have
@u
@t=sinx
cosht
@2u
@t2= sinx
sinht
@u
@x=cosx
sinht
@2u
@x2= sinx
sinht
Thus,@2u
@t2+@2u
@x2= 0:Moreover,u(x;0) = 0 and@
@tu(x;0) =sin x
:
(b) We have
supfj@
@tu(x;0) ut(x;0)j:x2Rg= supfsinx
:x2Rg
=supfsinx
:x2Rg=
7 SOLUTIONS AND RELATED TOPICS 59
and
supfju(x;t) u(x;t)j:x2Rg=2supfsinht
sinx
:x2Rg
=2sinht
:
(c) We have
lim
t!1supfju(x;t) u(x;t)j:x2Rg= lim
t!12sinht
=1:
Thus, a small change in the initial data leads to a catastrophically change in
the solution. Hence, the given problem is ill-posed
60 INTRODUCTION TO PDES
Practice Problems
Exercise 7.1
Determineaandbso thatu(x;y) =eax+byis a solution to the equation
uxxxx+uyyyy+ 2uxxyy= 0:
Exercise 7.2
Consider the following dierential equation
tuxx ut= 0:
Supposeu(t;x) =X(x)T(t):Show that there is a constant such that
X00=XandT00=tT:
Exercise 7.3
Consider the initial value problem
xux+ (x+ 1)yuy= 0; x;y> 1
u(1;1) =e:
Show thatu(x;y) =xex
yis the solution to this problem.
Exercise 7.4
Show thatu(x;y) =e 2ysin (x y) is the solution to the initial value prob-
lem ux+uy+ 2u= 0 forx;y> 0
u(x;0) = sinx
Exercise 7.5
Solve each of the following dierential equations:
(a)du
dx= 0 whereu=u(x):
(b)@u
@x= 0 whereu=u(x;y):
Exercise 7.6
Solve each of the following dierential equations:
(a)d2u
dx2= 0 whereu=u(x):
(b)@2u
@x@y= 0 whereu=u(x;y):
7 SOLUTIONS AND RELATED TOPICS 61
Exercise 7.7
Show thatu(x;y) =f(y+ 2x) +xg(y+ 2x);wherefandgare two arbitrary
twice dierentiable functions, satisfy the equation
uxx 4uxy+ 4uyy= 0:
Exercise 7.8
Find the dierential equation whose general solution is given by u(x;t) =
f(x ct)+g(x+ct);wherefandgare arbitrary twice dierentiable functions
in one variable.
Exercise 7.9
Letp:R!Rbe a dierentiable function in one variable. Prove that
ut=p(u)ux
has a solution satisfying u(x;t) =f(x+p(u)t);wherefis an arbitrary
dierentiable function. Then nd the general solution to ut= (sinu)ux:
Exercise 7.10
Find the general solution to the pde
uxx+ 2uxy+uyy= 0:
Hint: See Exercise 6.2.
62 INTRODUCTION TO PDES
Sample Exam Questions
Exercise 7.11
Letu(x;t) be a function such that uxxexists andu(0;t) =u(L;t) = 0 for all
t2R:Prove thatZL
0uxx(x;t)u(x;t)dx0:
Exercise 7.12
Consider the initial value problem
ut+uxx= 0; x2R; t> 0
u(x;0) = 1:
(a) Show that u(x;t)1 is a solution to this problem.
(b) Show that un(x;t) = 1 +en2t
nsinnxis a solution to the initial value
problem
ut+uxx= 0; x2R; t> 0
u(x;0) = 1 +sinnx
n:
(c) Find supfjun(x;0) 1j:x2Rg:
(d) Find supfjun(x;t) 1j:x2Rg:
(e) Show that the problem is ill-posed.
Exercise 7.13
Find the general solution of each of the following PDEs by means of direct
integration.
(a)ux= 3x2+y2; u=u(x;y):
(b)uxy=x2y; u =u(x;y):
(c)uxyz= 0; u=u(x;y;z ):
(d)uxtt=e2x+3t; u=u(x;t):
Exercise 7.14
Consider the second-order PDE
uxx+ 4uxy+ 4uyy= 0:
(a) Use the change of variables v(x;y) =y 2xandw(x;y) =xto show
thatuww= 0:
(b) Find the general solution to the given PDE.
7 SOLUTIONS AND RELATED TOPICS 63
Exercise 7.15
Derive the general solution to the PDE
utt=c2uxx
by using the change of variables v=x+ctandw=x ct:
64 INTRODUCTION TO PDES
First Order Partial Dierential
Equations
Many problems in the mathematical, physical, and engineering sciences deal
with the formulation and the solution of rst order partial dierential equa-
tions. Our rst task is to understand simple rst order equations. In ap-
plications, rst order partial dierential equations are most commonly used
to describe dynamical processes, and so time, t;is one of the independent
variables. Most of our discussion will focus on dynamical models in a single
space dimension, bearing in mind that most of the methods can be readily
extended to higher dimensional situations. First order partial dierential
equations and systems model a wide variety of wave phenomena, including
transport of solvents in
uids,
ood waves, acoustics, gas dynamics, glacier
motion, trac
ow, and also a variety of biological and ecological systems.
From a mathematical point of view, rst order partial dierential equations
have the advantage of providing conceptual basis that can be utilized in the
study of higher order partial dierential equations.
In this chapter we introduce the basic denitions of rst order partial dif-
ferential equations. We then derive the one dimensional spatial transport
eqution and discuss some methods of solutions. One general method of solv-
ability for quasilinear rst order partial dierential equation, known as the
method of characteristics, is analyzed.
8 Classication of First Order PDEs
In this section, we present the basic denitions pertained to rst order PDE.
By a rst order dierential equation in two variables xandywe mean
65
66 FIRST ORDER PARTIAL DIFFERENTIAL EQUATIONS
any equation of the form
F(x;y;u;u x;uy) = 0: (8.1)
In what follows the functions a;b; andcare assumed to be continuously
dierentiable functions. If Equation (8.1) can be written in the form
a(x;y;u )ux+b(x;y;u )uy=c(x;y;u ) (8.2)
then we say that the equation is quasilinear . The following are examples
of quasilinear equations:
uux+uy+cu2= 0
x(y2+u)ux y(x2+u)uy= (x2 y2)u:
If Equation (8.1) can be written in the form
a(x;y)ux+b(x;y)uy=c(x;y;u ) (8.3)
then we say that the equation is semilinear . The following are examples of
semilinear equations:
xux+yuy=u2+x2
(x+ 1)2ux+ (y 1)2uy= (x+y)u2:
If Equation (8.1) can be written in the form
a(x;y)ux+b(x;y)uy+c(x;y)u=d(x;y) (8.4)
then we say that the equation is linear . Examples of linear equations are:
xux+yuy=cu
(y z)yx+ (z x)uy+ (x y)uz= 0:
A rst order pde that is not linear is said to be nonlinear. Examples of
nonlinear equations are:
ux+cu2
y=xy
u2
x+u2
y=c:
First order partial dierential equations are classied as either linear or non-
linear. Clearly, linear equations are a special kind of quasilinear equation
8 CLASSIFICATION OF FIRST ORDER PDES 67
(8.2) ifaandbare functions of xandyonly andcis a linear function of u:
Likewise, semilinear equations are quasilinear equations if aandbare func-
tions ofxandyonly. Also, semilinear equations (8.4) reduces to a linear
equation if cis linear in u:
A linear equation is called homogeneous ifd(x;y)0 and nonhomoge-
neous ifd(x;y)6= 0:Examples of linear homogeneous equations are:
xux+yuy=cu
(y z)ux+ (z x)uy+ (x y)uz= 0:
Examples of nonhomogeneous equations are:
ux+ (x+y)uy u=ex
yux+xuy=xy:
Recall that for an ordinary linear dierential equation, the general solution
depends mainly on arbitrary constants. Unlike ODEs, in linear partial dif-
ferential equations, the general solution depends on arbitrary functions.
Example 8.1
Solve the equation ut(x;t) = 0:
Solution.
The general solution is given by u(x;t) =f(x) wherefis an arbitrary dif-
ferentiable function of x
Example 8.2
Consider the transport equation
aut(x;t) +bux(x;t) = 0
whereaandbare constants. Show that u(x;t) =f(bt ax) is a solution
to the given equation, where fis an arbitrary dierentiable function in one
variable.
Solution.
Letv(x;t) =bt ax:Using the chain rule we see that ut(x;t) =bfv(v) and
ux(x;t) = afv(v):Hence,aut(x;t) +bux(x;t) =abfv(v) abfv(v) = 0
68 FIRST ORDER PARTIAL DIFFERENTIAL EQUATIONS
Practice Problems
Exercise 8.1
Classify each of the following PDE as linear, quasilinear, semi-linear, or non-
linear.
(a)xux+yuy= sin (xy):
(b)ut+uux= 0
(c)u2
x+u3u4
y= 0:
(d) (x+ 3)ux+xy2uy=u3:
Exercise 8.2
Show thatu(x;y) =exf(2x y);wherefis a dierentiable function of one
variable, is a solution to the equation
ux+ 2uy u= 0:
Exercise 8.3
Show thatu(x;y) =xpxysatises the equation
xux yuy=u
subject to
u(y;y) =y2; y0:
Exercise 8.4
Show thatu(x;y) = cos (x2+y2) satises the equation
yux+xuy= 0
subject to
u(0;y) = cosy2:
Exercise 8.5
Show thatu(x;y) =y 1
2(x2 y2) satises the equation
1
xux+1
yuy=1
y
subject tou(x;1) =1
2(3 x2):
8 CLASSIFICATION OF FIRST ORDER PDES 69
Exercise 8.6
Find a relationship between aandbifu(x;y) =f(ax+by) is a solution to
the equation 3 ux 7uy= 0 for any dierentiable function f:
Exercise 8.7
SupposeLis a linear operator, that is, L(u+v) =L(u)+L(v):Consider
the homogeneous and nonhomogeneous linear equations
Lu= 0
Lu=f
wherefis some function. Suppose vis a solution to the homogeneous equa-
tion, andwis a solution to the nonhomogeneous equation. Show u=av+w
is a solution to the nonhomogeneous equation for any constant a:
Exercise 8.8
Reduce the partial dierential equation
aux+buy+cu= 0
to a rst order ODE by introducing the change of variables s=ax+byand
t=bx ay:
Exercise 8.9
Solve the partial dierential equation
ux+uy= 1
by introducing the change of variables s=x+yandt=x y:
70 FIRST ORDER PARTIAL DIFFERENTIAL EQUATIONS
Sample Exam Questions
Exercise 8.10
Show thatu(x;y) =e 4xf(2x 3y) is a solution to the rst-order PDE
3ux+ 2uy+ 12u= 0:
Exercise 8.11
Derive the general solution of the PDE
aut+bux=u; a;b6= 0
by using the change of variables v=ax btandw=1
at:
Exercise 8.12
Derive the general solution of the PDE
aux+buy= 0; a;b6= 0
by using the change of variables s(x;y) =ax+byandt(x;y) =bx ay:
Assumea2+b2>0:
Exercise 8.13
Write the equation
ut+cux+u=f(x;y)
in the coordinates v=x ct; w =t:
Exercise 8.14
Suppose that u(x;t) =w(x ct) is a solution to the PDE
xux+tut=Au
whereAandcare constants. Let v=x ct:Write the dierential equation
with unknown function w(v):
9 THE ONE DIMENSIONAL SPATIAL TRANSPORT EQUATIONS 71
9 The One Dimensional Spatial Transport Equa-
tions
Modeling is the process of writing a dierential equation to describe a physi-
cal situation. In this section we discuss the one-dimensional transport equa-
tion and discuss an analytical method for solving it.
Linear Transport Equation for Fluid Flows
We shall describe the transport of a dissolved chemical by water that is trav-
eling with uniform velocity cthrough a long thin tube Gwith uniform cross
sectionA:(The very same discussion applies to the description of the trans-
port of gas by air moving through a pipe.) We identify Gwith the open
interval (a;b);and the velocity c>0 is in the (rightward) positive direction
of thex axis. We will assume that the concentration of the chemical is con-
stant across the cross section Aat each point xso that the chemical changes
in thex direction and thus the term one-dimensional spatial equation. See
Figure 9.1
Figure 9.1
Letu(x;t) be a continuously dierentiable function denoting the concentra-
tion of the chemical (i.e. amount of chemical/area) at position xat time
t:Then at time t;the amount of chemical stored in a section of the tube
between positions aandxis given by the denite integral
Zx
aAu(s;t)ds:
72 FIRST ORDER PARTIAL DIFFERENTIAL EQUATIONS
Since the water is
owing at a speed c;so at timeh+tthe same quantity of
chemical will be
Zx
aAu(s;t)ds=Zx+ch
a+chAu(s;t+h)ds:
Taking the derivative of both sides with respect to xwe nd
u(x;t) =u(x+ch;t+h):
Now taking the derivative of this last equation with respect to hwe nd
0 =ut(x+ch;t+h) +cux(x+ch;t+h):
Taking the limit of this last equation as happroaches 0 we nd
ut(x;t) +cux(x;t) = 0 (9.1)
for all (x;t):This equation is called the transport equation in one-dimensional
space. It is a linear, homogeneous rst order partial dierential equation.
Example 9.1
Show thatu(x;t) =f(x ct) is a solution to (9.1), where fis an arbitrary
dierentiable function in one variable.
Solution.
Using the chain rule we nd
ut= cf0(x ct) andux=f0(x ct):
Hence, by substituting these results into the equation we nd
ut+cux= cf0(x ct) +cf0(x ct) = 0:
The solution u(x;t) =f(x ct) is called the right traveling wave, since
the graph of the function f(x ct) at a given time tis the graph of f(x)
shifted to the right by the value ct:Thus, with growing time, the function
f(x) is moving without changes to the right at the speed c
An initial value condition determines a unique solution to the transport equa-
tion as stated in the next theorem.
9 THE ONE DIMENSIONAL SPATIAL TRANSPORT EQUATIONS 73
Theorem 9.4
Letgbe a continuously dierentiable function. Then there is a unique con-
tinuously dierentiable solution u(x;t) to the IVP
aux(x;t) +but(x;t) = 0
u(x;0) =g(x):
Indeed,uis given explicitly by the formula
u(x;t) =f(bx at); g(x) =f(bx):
Method of Solutions: The Coordinate Method
We will solve (9.1) by solving the more general equation
aux+buy= 0 (9.2)
wherea2+b2>0:
We introduce a new rectangular system by the substitution
s=ax+by; t =bx ay
According to the chain rule for the derivative of a composite function, we
have
ux=ussx+uttx=aus+but
uy=ussy+utty=bus aut
Substituting these into (9.2) to obtain
a2us+abut+b2us abut= 0
or
(a2+b2)us= 0
and sincea2+b2>0 we obtain
us= 0:
Solving this equation, we nd
u(s;t) =f(t)
wherefis an arbitrary dierentiable function of one variable. Now, in terms
ofxandywe nd
u(x;y) =f(bx ay):
74 FIRST ORDER PARTIAL DIFFERENTIAL EQUATIONS
Example 9.2
Use the coordinate method to nd the solution to ut 3ux= 0; u(x;0) =e x2:
Solution.
Letv= 3x+tandw=x+ 3t:Thenux= 3uv+uwandut=uv+ 3uw:
Substituting these into the given equation we nd 10 uv= 0 oruv= 0:
Hence,u(v;w) =f(w) oru(x;t) =f(x+ 3t) wherefis a dierentiable
function in one variable. Since u(x;0) =e x2;we nde x2=f(x):Hence,
u(x;t) =e (x+3t)2
Transport Equation with Decay: The Method of Characteristic
Coordinates
Atransport equation with decay is an equation given by
ut+cux+u=f(x;t)
whereandcare constants and fis a given function representing external
resources. Note that the decay is characterized by the term u:
To solve this equation, we introduce the characteristic coordinates given
by
v=x ct; w =t:
Using the chain rule, we nd
ut=uvvt+uwwt= cuv+uw
ux=uvvx+uwwx=uv:
Substituting these into the original equation we obtain the equation
uw+u=f(v+cw;w )
which can be solved by the method of integrating factor. We illustrate this
approach in the next example.
Example 9.3
Find the general solution of the transport equation
ut+ux u=t:
9 THE ONE DIMENSIONAL SPATIAL TRANSPORT EQUATIONS 75
Solution.
The characteristic coordinates are
v=x t; w =t:
These transform the original equation to the rst order ODE
uw u=w:
Using the method of integrating factor, we nd
d
dw(e wu) =we w
and solving this equation we nd
u(v;w) = (1 +w) +ewf(v)
and in terms of xandtwe nd
u(x;t) =f(x t)et (1 +t)
A more general method for solving quasilinear rst order partial dierential
equations, known as the method of characteristics, will be discussed in the
next section.
76 FIRST ORDER PARTIAL DIFFERENTIAL EQUATIONS
Practice Problems
Exercise 9.1
Use the coordinate method to nd the solution to ut+ 3ux= 0; u(x;0) =
sinx:
Exercise 9.2
Use the coordinate method, solve the equation aux+buy+cu= 0:
Exercise 9.3
Use the coordinate method, solve the equation ux+ 2uy= cos (y 2x) with
the initial condition u(0;y) =f(y);wheref:R!Ris a given function.
Exercise 9.4
Show that the initial value problem ut+ux=x; u(x;x) = 1 has no solution.
Exercise 9.5
Solve the transport equation ut+ 2ux= 3uwith initial condition u(x;0) =
1
1+x2:
Exercise 9.6
Solveut+ux 3u=twith initial condition u(x;0) =x2:
Exercise 9.7
Show that the decay term uin the transport equation with decay
ut+cux+u= 0
can be eliminated by the substitution w=uet:
Exercise 9.8
Use the coordinate method to solve
ux+uy=u2
u(x;0) =h(x)
9 THE ONE DIMENSIONAL SPATIAL TRANSPORT EQUATIONS 77
Exercise 9.9 (Well-Posed)
Letube the unique solution to the IVP
ut+cux= 0
u(x;0) =f(x)
andvbe the unique solution to the IVP
ut+cux= 0
u(x;0) =g(x)
wherefandgare continuously dierentiable functions.
(a) Show that w(x;t) =u(x;t) v(x;t) is the unique solution to the IVP
ut+cux= 0
u(x;0) =f(x) g(x)
(b) Write an explicit formula for win terms of fandg:
(c) Use (b) to conclude that the transport problem is well-posed. That is, a
small change in the initial data leads to a small change in the solution.
Exercise 9.10
Solve the initial boundary value problem
ut+cux= u; x> 0; t> 0
u(x;0) = 0; u(0;t) =g(t); t> 0:
78 FIRST ORDER PARTIAL DIFFERENTIAL EQUATIONS
Sample Exam Questions
Exercise 9.11
Solve the rst-order equation 2 ut+3ux= 0 with the initial condition u(x;0) =
sinx:
Exercise 9.12
Solve the PDE
ux+uy= 1
using the coordinate method.
Exercise 9.13
Consider the rst order linear homogeneous PDE
Aux+Buy+Cu= 0
whereA;B; andCare constants with A6= 0:
(a) Determine a;b;c;d in terms of A;B;C such thatad bc6= 0 and so that
the change of variables v=ax+byandw=cx+dywill reduce the given
PDE to a rst order PDE of the form uv+u= 0:
(b) Use (a) to nd the general solution of the given PDE.
Exercise 9.14
Use the result of the previous problem to solve the PDE
ux+uy+u= 0:
10 THE METHOD OF CHARACTERISTICS 79
10 The Method of Characteristics
In this section we develop a method for nding the general solution of a
quasilinear rst order partial dierential equation. This method is called
themethod of characteristics orLagrange's method. This method of
solution can be described by the following result.
Theorem 10.1
The general solution of the quasilinear rst order PDE
a(x;y;u )ux+b(x;y;u )uy=c(x;y;u ) (10.1)
is given by
f(v;w) = 0 (10.2)
wherefis an arbitrary dierentiable function of v(x;y;u ) andw(x;y;u ) and
v=constant= c1; w=constant= c2are solutions to the ODE system
dx
a=dy
b=du
c: (10.3)
Equations (10.3) are called the characteristic equations in non-parametric
forms. The corresponding parametric forms are given by the system of ODEs
dx
ds=a
dy
ds=b
du
ds=c
Remark 10.1
Sometimes (10.2) is written explicitly as v=g(w) orw=g(v) wheregis an
arbitrary dierentiable function.
Example 10.1
Find the general solution of the PDE x2ux+y2uy= (x+y)u:
Solution.
The characteristic equations for this PDE aredx
x2=dy
y2=du
(x+y)u:Using the
80 FIRST ORDER PARTIAL DIFFERENTIAL EQUATIONS
rst two fractions, we haveRdx
x2=Rdy
y2and this impliesx y
xy=c1:Also,
we can solve for xobtainingx=1
1
y c1andx+y=y
1 c1y+y:Using the
last two fractions we nddy
y2=du
u(x+y)=)1
y2
y
1 c1y+y
dy=du
u=)
1
y c1y2+1
y
dy=du
u=)
1 c1y
y c1y2+c1
1 c1y+1
y
dy=du
u=)R
2
y+c1
1 c1y
dy=
Rdu
u=)2 lnjyj lnj1 c1yj= lnjuj+c2=)y2
1 c1y=cu=)yy
1 c1y=
cu=)y1
1
y c1=cu=)xy=cu:Hence, the general solution is
fx y
xy;xy
u
= 0
wherefis an arbitrary dierentiable function
Example 10.2
Find the general solution of the PDE yuux+xuuy=xy:
Solution.
The characteristic equations aredx
yu=dy
xu=du
xy:Using the rst two fractions
we ndx2 y2=c1:Using the last two fractions we nd u2 y2=c2:Hence,
the general solution is f(x2 y2;u2 y2) = 0 oru2=y2+g(x2 y2);where
fandgare arbitrary dierentiable functions
Example 10.3
Find the general solution of the PDE x(y2 u2)ux y(u2+x2)yy= (x2+y2)u:
Solution.
The characteristic equations aredx
x(y2 u2)=dy
y(u2+x2)=du
(x2+y2)u:Using a
property of proportions we can write
xdx+ydy+udu
x2(y2 u2) y2(u2+x2) +u2(x2+y2)=du
(x2+y2)u:
That is
xdx+ydy+udu
0=du
(x2+y2)u
or
xdx+ydy+udu= 0:
10 THE METHOD OF CHARACTERISTICS 81
Hence, we nd x2+y2+u2=c1:Also,
dx
x dy
y
y2 u2+u2+x2=du
(x2+y2)u
or
dx
x dy
y=du
u:
Hence, we nd lnyu
x=constant oryu
x=c2:The general solution is given by
f
x2+y2+u2;yu
x
= 0
or
u=x
yg(x2+y2+u2)
wherefandgare arbitrary dierentiable functions
Example 10.4
Solve the transport equation using the method of characteristics
ut+cux= 0:
Solution.
The characteristic equations are given by
dt
1=dx
c=du
0:
Solving the rst two fractions we nd x ct=k:The last fraction implies
u=k0:The general solution is given by f(x ct;u) = 0 oru=g(x ct)
Solution curves to the ODE
dy
dx=b
a
are called characteristic curves or simply characteristics. These are
curves in the xy plane.
Example 10.5
Find the characteristics of cos yux+uy+xu= 0:
82 FIRST ORDER PARTIAL DIFFERENTIAL EQUATIONS
Solution.
Solving the equationdy
dx=1
cosyby the separation of variable method we nd
siny x=k
Example 10.6
Find the characteristics of ux+ 2uy u= 0:
Solution.
We havea= 1 andb= 2:Thus,dy
dx= 2 so that the characteristics are given
by 2x y=k
10 THE METHOD OF CHARACTERISTICS 83
Practice Problem
Exercise 10.1
Find the characteristics of the PDE
xux yuy=u:
Exercise 10.2
Find the characteristics of the PDE
yux+xuy= 0:
Exercise 10.3
Find the characteristics of the PDE
(x+y)(ux+uy) =u 1:
Exercise 10.4
Find the general solution of the PDE xux+yuy= 1 +u2:
Exercise 10.5
Find the general solution of the PDE ln ( y+u)ux+uy= 1:
Exercise 10.6
Find the general solution of the PDE xux+yuy=u:
Exercise 10.7
Find the general solution of the PDE xux+yuy=nu:
Exercise 10.8
Find the general solution of the PDE x(y u)ux+y(u x)uy=u(x y):
Exercise 10.9
Find the general solution of the PDE u(u2+xy)(xux yuy) =x4:
Exercise 10.10
Find the general solution of the PDE ( y+xu)ux (x+yu)uy=x2 y2:
Exercise 10.11
Find the general solution of the PDE ( y2+u2)ux xyuy+xu= 0:
84 FIRST ORDER PARTIAL DIFFERENTIAL EQUATIONS
Exercise 10.12
Find the general form of solutions to
ux+ 2uy=u
and sketch some of the characteristics. Hint: dene a new variable v=e xu:
What equation does vsatisfy?
Exercise 10.13
Find the general form of solutions to
(1 +x2)ux+uy= 0
and sketch some of the characteristics.
10 THE METHOD OF CHARACTERISTICS 85
Sample Exam Questions
Exercise 10.14
Find the general solution of the equation
ux+yuy=u:
Exercise 10.15
Find the characteristics associated with the PDE
ux+xuy+ 3u= 2:
Exercise 10.16
Find the general solution of the rst order PDE
ux+yuy+xu= 0:
Exercise 10.17
Find the characteristics of the PDE
1
xux+1
yuy= 0:
Exercise 10.18
Find the characteristics of the PDE
1
xux+1
yuy=1
y:
86 FIRST ORDER PARTIAL DIFFERENTIAL EQUATIONS
11 The Cauchy Problem for First Order Quasi-
linear Equations
When solving a partial dierential equation, it is seldom the case that one
tries to study the properties of the general solution of such equations. In
general, one deals with those partial dierential equations whose solutions
satisfy certain supplementary conditions. In the case of a rst order partial
dierential equation, we determine the particular solution by formulating an
initial value porblem also known as a Cauchy problem.
In this section, we discuss the Cauchy problem for the rst order quasilinear
partial dierential equation
a(x;y;u )ux+b(x;y;u )uy=c(x;y;u ): (11.1)
Recall that the initial value problem of a rst order ordinary dierential
equation asks for a solution of the equation which has a given value at a
given point in R:The Cauchy problem for the PDE (11.1) asks for a solution
of (11.1) which has given values on a given curve in R2:A precise statement
of the problem is given next.
Initial Value Problem or Cauchy Problem
LetCbe a given curve in R2dened parametrically by the equations
x=x0(t); y=y0(t)
wherex0;y0are continuously dierentiable functions on some interval I:Let
u0(t) be a given continuously dierentiable function on I:The Cauchy prob-
lem for (11.1) asks for a continuously dierentiable function u=u(x;y)
dened in a domain
R2containing the curve Cand such that:
(1)u=u(x;y) is a solution of (11.1) in
.
(2) On the curve C; u equals the given function u0(t);i.e.
u(x0(t);y0(t)) =u0(t); t2I: (11.2)
We callCtheinitial curve of the problem, u0(t) the initial data , and
(11.2) the initial condition of the problem. See Figure 11.1.
11 THE CAUCHY PROBLEM FOR FIRST ORDER QUASILINEAR EQUATIONS 87
Figure 11.1
If we view a solution u=u(x;y) of (11.1) as an integral surface of (11.1),
we can give a simple geometrical statement of the problem: Find a solu-
tion surface of (11.1) containing the curve described parametrically by the
equations
:x=x0(t); y=y0(t); u=u0(t); t2I:
Note that the projection of this curve in the xy plane is just the curve C:
The following theorem asserts that under certain conditions the Cauchy prob-
lem (11.1) - (11.2) has a unique solution.
Theorem 11.1
Suppose that x0(t);y0(t);andu0(t) are continuously dierentiable functions
oftin an interval I;and thata;b;andcare functions of x;y; anduwith
continuous rst order partial derivatives with respect to their argument in
some domain Dof (x;y;u ) space containing the initial curve
:x=x0(t); y=y0(t); u=u0(t)
wheret2I:Then for each point ( x0(t);y0(t);u0(t)) on that satises the
condition
a(x0(t);y0(t);u0(t))dy0
dt(t) b(x0(t);y0(t);u0(t))dx0
dt(t)6= 0: (11.3)
there exists a unique solution u=u(x;y) of (11.1) in a neighborhood Uof
(x0(t);y0(t)) such that the initial condition (11.2) is satised for every point
onCcontained in U:See Figure 11.2.
88 FIRST ORDER PARTIAL DIFFERENTIAL EQUATIONS
Figure 11.2
Note that condition (11.3) implies that
dy0(t)
dx0(t)6=b(x0;y0;u0)
a(x0;y0;u0)
which means that the vector ( a(x0;y0;u0);b(x0;y0;u0);c(x0;y0;u0)) is not
tangent to :(Recall that the normal vector to Chas components
dy0(t)
dt; dx0(t)
t
so that a vector ( a;b) is tangent to Cifady0(t)
dt bdx0(t)
dt= 0:) It follows that
the Cauchy problem has a unique solution if Cis nowhere characteristic.
We construct the desired solution using the method of characteristics as fol-
lows: Pick a point ( x0(t);y0(t);u0(t))2 :Using this as the initial value
we solve the system of ODEs consisting of the characteristic equations in
parametric form
dx
ds=a(x(s);y(s);u(s))
dy
ds=b(x(s);y(s);u(s))
du
ds=c(x(s);y(s);u(s))
satisfying the initial condition
(x(0);y(0);u(0)) = (x0(t);y0(t);u0(t)):
11 THE CAUCHY PROBLEM FOR FIRST ORDER QUASILINEAR EQUATIONS 89
The solution depends on the parameter sso it consists of a triples of functions
x=x(s;t); y=y(s;t); u=u(s;t): (11.4)
This system represents the parametric representation of the integral surface
of the problem in which the curve corresponds to s= 0:The solution uis
recovered by solving the rst two equations in (11.4) for
t=t(x;y); s=s(x;y)
and substituting these into the third equation to obtain u(x;y) =u(s(x;y);t(x;y)):
Example 11.1
Solve the Cauchy problem
ux+uy=1
u(x;0) =f(x):
Solution.
The initial curve in R3can be given parametrically as
:x0(t) =t; y 0(t) = 0; u0(t) =f(t):
We have
a(x0(t);y0(t);u0(t))dy0
dt(t) b(x0(t);y0(t);u0(t))dx0
dt(t) = 16= 0
so by the above theorem the given Cauchy problem has a unique solution.
To nd this solution, we solve the system of ODEs
dx
ds=1
dy
ds=1
du
ds=1:
Solving this system we nd
x(s;t) =s+(t); y(s;t) =s+(t); u(s) =s+
(t):
90 FIRST ORDER PARTIAL DIFFERENTIAL EQUATIONS
Butx(0;t) =tso that(t) =t:Similarly,y(0;t) = 0 so that (t) = 0
andu(0;t) =f(t) implies
(t) =f(t):Hence, the unique solution is given
parametrically by the equations
x(s;t) =t+s; y (s;t) =s; u (s;t) =s+f(t):
Solving the rst two equations for sandtwe nd
s=y; t =x y
and substituting these into the third equation we nd
u(x;y) =y+f(x y):
Alternative Computation
We can apply the results of the previous section to nd the unique solution.
If we solve the characteristic equations in non-parametric form
dx
1=dy
1=du
1
we ndx y=c1andu x=c2:Thus, the general solution of the PDE
is given by u=x+F(x y):Using the Cauchy data u(x;0) =f(x) we nd
f(x) =x+F(x) which implies that F(x) =f(x) x:Hence, the unique
solution is given by
u(x;y) =x+f(x y) (x y) =y+f(x y)
If condition (11.3) is not satised than Cis a characteristic curve. If the
curve satises the characteristic equations than the problem has innitely
many solutions. To see this, pick an arbitrary point P0= (x0;y0;u0) on :
Pick a new initial curve 0passing through P0which is not tangent to at
P0:In this case, condition (11.3) is satised and the new Cauchy problem has
a unique solution. Since there are innitely many ways of selecting 0;we
obtain innitely many solutions. We illustrate this case in the next example.
Example 11.2
Solve the Cauchy problem
ux+uy=1
u(x;x) =x:
11 THE CAUCHY PROBLEM FOR FIRST ORDER QUASILINEAR EQUATIONS 91
Solution.
The initial curve in R3can be given parametrically as
:x0(t) =t; y 0(t) =t; u 0(t) =t:
We have
a(x0(t);y0(t);u0(t))dy0
dt(t) b(x0(t);y0(t);u0(t))dx0
dt(t) = 0:
As in Example 11.1, the general solution of the PDE is u(x;y) =y+f(x
y) wherefis an arbitrary dierentiable function. Using the Cauchy data
u(x;x) =xwe ndf(0) = 0:Thus, the solution is given by
u(x;y) =y+f(x y)
wherefis an arbitrary function such that f(0) = 0:There are innitely
many choices for f:Hence, the problem has innitely many solutions. Note
that satises the characteristic equations
If condition (11.3) is not satised and if does not satisfy the characteristic
equations then it can be shown that the Cauchy problem has no solutions.
We illustrate this case next.
Example 11.3
Solve the Cauchy problem
ux+uy=1
u(x;x) =1:
Solution.
The initial curve in R3can be given parametrically as
:x0(t) =t; y 0(t) =t; u 0(t) = 1:
We have
a(x0(t);y0(t);u0(t))dy0
dt(t) b(x0(t);y0(t);u0(t))dx0
dt(t) = 0:
Solving the characteristic equations in parametric form we nd
x(s;t) =s+(t); y(s;t) =s+(t); u(s;t) =s+
(t):
92 FIRST ORDER PARTIAL DIFFERENTIAL EQUATIONS
Clearly, does not satisfy the characteristic equations. Now, the general
solution to the PDE is given by u=y+f(x y):Using the Cauchy data
u(x;x) = 1 we nd f(0) = 1 x;which is not possible since the LHS is a
xed number whereas the RHS is a variable expression. Hence, the problem
has no solutions
Example 11.4
Solve the Cauchy problem
ux uy=1
u(x;0) =x2: (11.5)
Solution.
The initial curve is given parametrically by
:x0(t) =t; y 0(t) = 0; u 0(t) =t2:
We have
a(x0(t);y0(t);u0(t))dy0
dt(t) b(x0(t);y0(t);u0(t))dx0
dt(t) = 16= 0
so the Cauchy problem has a unique solution.
The characteristic equations are
dx
1=dy
1=du
1:
Using the rst two fractions we nd x+y=c1:Using the rst and the third
fractions we nd u x=c2:Thus, the general solution can be represented
by
u=x+f(x+y)
wherefis an arbitrary dierentiable function. Using the Cauchy data
u(x;0) =x2we ndx2 x=f(x):Hence, the unique solution is given
by
u=x+ (x+y)2 (x+y) = (x+y)2 y
Example 11.5
Solve the initial value problem
ut+uux=x; u (x;0) = 1:
11 THE CAUCHY PROBLEM FOR FIRST ORDER QUASILINEAR EQUATIONS 93
Solution.
The initial curve is given parametrically by
:x0(t) =t; y 0(t) = 0; u 0(t) = 1:
We have
a(x0(t);y0(t);u0(t))dy0
dt(t) b(x0(t);y0(t);u0(t))dx0
dt(t) = 16= 0
so the Cauchy problem has a unique solution.
The characteristic equations are
dt
1=dx
u=du
x:
Since
dt
1=d(x+u)
x+u
we nd that ( x+u)e t=c1:Now, using the last two fractions we nd
u2 x2=c2:Hence, the general solution is given by
f((x+u)e t;u2 x2) = 0
wherefis an arbitrary dierentiable function. Using the Cauchy data we
ndc1= 1 +xandc2= 1 x2= 2(1 +x) (1 +x)2= 2c1 c2
1:Thus,
u2 x2= 2(x+u)e t (x+u)2e 2t
or
u x= 2e t (x+u)e 2t:
This can be reduced further as follows: u+ue 2t=x+ 2e t xe 2t=
2e t+x(1 e 2t) =)u=2e t
1+e 2t+x1 e 2t
1+e 2t= sech(t) +xtanh(t)
Example 11.6
Solve the initial value problem
uux+uy= 1
with the initial curve
:x0(t) = 2t2; y 0(t) = 2t; u 0(t) = 0; t> 0:
94 FIRST ORDER PARTIAL DIFFERENTIAL EQUATIONS
Solution.
We have
a(x0(t);y0(t);u0(t))dy0
dt(t) b(x0(t);y0(t);u0(t))dx0
dt(t) = 4t6= 0; t> 0
so the Cauchy problem has a unique solution.
The characteristic equations in parametric form are given by the system of
ODEs
dx
ds=u
dy
ds=1
du
ds=1:
Thus, the solution of this system depends on two parameters sandt:Solving
the last two equations we nd
y(s;t) =s+(t); u(s;t) =s+
(t):
Solving the rst equation with ubeing replaced by s+
(t) we nd
x(s;t) =1
2s2+
(t)s+(t):
Using the initial conditions
x(0;t) = 2t2; y(0;t) = 2t; u(0;t) = 0
we nd
x(s;t) =1
2s2+ 2t2; y(s;t) =s+ 2t; u(s;t) =s:
Eliminating sandtwe nd
(u y)2+u2= 2x:
Solving this quadratic equation in uto nd
2u=y(4x y2)1
2:
The solution surface satisfying u= 0 ony2= 2xis given by
2u=y (4x y2)1
2:
This represents a solution surface only when y2<4x:The solution does not
exist fory2>4x
11 THE CAUCHY PROBLEM FOR FIRST ORDER QUASILINEAR EQUATIONS 95
Practice Problems
Exercise 11.1
Solve
(y u)ux+ (u x)uy=x y
with the condition u(x;1
x) = 0:
Exercise 11.2
Solve the linear equation
yux+xuy=u
with the Cauchy data u(x;0) =x3:
Exercise 11.3
Solve
x(y2+u)ux y(x2+u)uy= (x2 y2)u
with the Cauchy data u(x; x) = 1:
Exercise 11.4
Solve
xux+yuy=xe u
with the Cauchy data u(x;x2) = 0:
Exercise 11.5
Solve the initial value problem
xux+uy= 0; u(x;0) =f(x)
using the characteristic equations in parametric form.
Exercise 11.6
Solve the initial value problem
ut+aux= 0; u(x;0) =f(x):
Exercise 11.7
Solve the initial value problem
aux+uy=u2; u(x;0) = cosx
96 FIRST ORDER PARTIAL DIFFERENTIAL EQUATIONS
Exercise 11.8
Solve the initial value problem
ux+xuy=u; u (1;y) =h(y):
Exercise 11.9
Solve the initial value problem
uux+uy= 0; u(x;0) =f(x):
Exercise 11.10
Solve the initial value problem
p
1 x2ux+uy= 0; u(0;y) =y:
11 THE CAUCHY PROBLEM FOR FIRST ORDER QUASILINEAR EQUATIONS 97
Sample Exam Questions
Exercise 11.11
Consider
xux+ 2yuy= 0:
(i) Find and sketch the characteristics.
(ii) Find the solution with u(1;y) =ey:
(iii) What happens if you try to nd the solution satisfying either u(0;y) =
g(y) oru(x;0) =h(x) for given functions gandh?
(iv) Explain, using your picture of the characteristics, what goes wrong at
(x;y) = (0;0):
Exercise 11.12
Solve the equation ux+uy=usubject to the condition u(x;0) = cosx:
Exercise 11.13
(a) Find the general solution of the equation
ux+yuy=u:
(b) Find the solution satisfying the Cauchy data u(x;3ex) = 2:
(c) Find the solution satisfying the Cauchy data u(x;ex) =ex:
Exercise 11.14
Solve the Cauchy problem
ux+ 4uy=x(u+ 1)
u(x;5x) = 1:
Exercise 11.15
Solve the Cauchy problem
ux uy=u
u(x; x) = sinx:
98 FIRST ORDER PARTIAL DIFFERENTIAL EQUATIONS
Exercise 11.16
(a) Find the characteristics of the equation
yux+xuy= 0:
(b) Sketch some of the characteristics.
(c) Find the solution satisfying the boundary condition u(0;y) =e y2:
(d) In which region of the plane is the solution uniquely determined?
Exercise 11.17
Consider the equation ux+yuy= 0:Is there a solution satisfying the extra
condition
(a)u(x;0) = 1
(b)u(x;0) =x?
If yes, give a formula; if no, explain why.
Second Order Linear Partial
Dierential Equations
In this chapter we consider the three fundamental second order linear partial
dierential equations of parabolic, hyperbolic, and elliptic type. These types
arise in many applications such as the wave equation, the heat equation
and the Laplace's equation. We will study the solvability of each of these
equations.
12 Second Order PDEs in Two Variables
In this section we will brie
y review second order partial dierential equa-
tions.
Asecond order partial dierential equation in the variables xandy
is an equation of the form
F(x;y;u;u x;uy;uxx;uyy;uxy) = 0: (12.1)
If Equation (12.1) can be written in the form
A(x;y;u;ux;uy)uxx+B(x;y;u;ux;uy)uxy+C(x;y;u;ux;uy)uyy=D(x;y;u;ux;uy)
(12.2)
then we say that the equation is quasilinear .
If Equation (12.1) can be written in the form
A(x;y)uxx+B(x;y)uxy+C(x;y)uyy=D(x;y;u;u x;uy) (12.3)
then we say that the equation is semilinear .
If Equation (12.1) can be written in the form
A(x;y)uxx+B(x;y)uxy+C(x;y)uyy+D(x;y)ux+E(x;y)uy+F(x;y)u=G(x;y)
(12.4)
99
100SECOND ORDER LINEAR PARTIAL DIFFERENTIAL EQUATIONS
then we say that the equation is linear .
A linear equation is said to be homogeneous whenG(x;y)0 and non-
homogeneous otherwise.
Equation (12.4) resembles the general equation of a conic section
Ax2+Bxy +Cy2+Dx+Ey+F= 0
which is classied as either parabolic, hyperbolic, or elliptic based on the sign
of the discriminant B2 4AC: We do the same for a second order linear
partial dierential equation:
Hyperbolic: This occurs if B2 4AC > 0 at a given point in the domain
ofu:
Parabolic: This occurs if B2 4AC= 0 at a given point in the domain
ofu:
Elliptic: This occurs if B2 4AC < 0 at a given point in the domain of
u:
Example 12.1
Determine whether the equation uxx+xuyy= 0 is hyperbolic, parabolic or
elliptic.
Solution.
Here we are given A= 1;B= 0;andC=x:SinceB2 4AC= 4x;the
given equation is hyperbolic if x<0, parabolic if x= 0 and elliptic if x>0
Second order partial dierential equations arise in many areas of scientic
applications. In what follows we list some of the well-known models that are
of great interest:
1. The heat equation in one-dimensional space is given by
ut=kuxx
wherekis a constant.
2. The wave equation in one-dimensional space is given by
utt=c2uxx
wherecis a constant.
3. The Laplace equation is given by
u=uxx+uyy= 0:
12 SECOND ORDER PDES IN TWO VARIABLES 101
Practice Problems
Exercise 12.1
Classify each of the following equation as hyperbolic, parabolic, or elliptic:
(a) Wave propagation: utt=c2uxx; c> 0:
(b) Heat conduction: ut=cuxx; c> 0:
(c) Laplace's equation: u=uxx+uyy= 0:
Exercise 12.2
Classify the following linear scalar PDE with constant coecents as hyper-
bolic, parabolic or elliptic.
(a)uxx+ 4uxy+ 5uyy+ux+ 2uy= 0:
(b)uxx 4uxy+ 4uyy+ 3ux+ 4u= 0:
(c)uxx+ 2uxy 3uyy+ 2ux+ 6uy= 0:
Exercise 12.3
Find the region(s) in the xy plane where the equation
(1 +x)uxx+ 2xyuxy y2uyy= 0
is elliptic, hyperbolic, or parabolic. Sketch these regions.
Exercise 12.4
Show thatu(x;t) = cosxsintis a solution to the problem
utt=uxx
u(x;0) = 0
ut(x;0) = cos x
ux(0;t) = 0
for allx;t> 0:
Exercise 12.5
Classify each of the following PDE as linear, quasilinear, semi-linear, or non-
linear.
(a)ut+uux=uuxx
(b)xutt+tuxx+u3u2
x=t+ 1
(c)utt=c2uxx
(d)u2
tt+ux= 0:
102SECOND ORDER LINEAR PARTIAL DIFFERENTIAL EQUATIONS
Exercise 12.6
Show that, for all ( x;y)6= (0;0); u(x;y) = ln (x2+y2) is a solution of
uxx+uyy= 0;
and that, for all ( x;y;z )6= (0;0;0); u(x;y;z ) =1p
x2+y2+z2is a solution of
uxx+uyy+uzz= 0:
Exercise 12.7
Consider the eigenvalue problem
uxx=u; 0<x<L
ux(0) =k0u(0)
ux(L) = kLu(L)
with Robin boundary conditions, where k0andkLare given positive numbers
andu=u(x):Can this system have a nontrivial solution u60 for>0?
Hint: Multiply the rst equation by uand integrate over x2[0;L]:
Exercise 12.8
Show that u(x;y) =f(x)g(y), wherefandgare arbitrary dierentiable
functions, is a solution to the PDE
uuxy=uxuy:
Exercise 12.9
Show that for any n2N;the function un(x;y) = sinnxsinhnyis a solution
to the Laplace equation
u=uxx+uyy= 0:
Exercise 12.10
Solve
uxy=xy:
12 SECOND ORDER PDES IN TWO VARIABLES 103
Sample Exam Questions
Exercise 12.11
Classify each of the following second-oder PDEs according to whether they
are hyperbolic, parabolic, or elliptic:
(a) 2uxx 4uxy+ 7uyy u= 0:
(b)uxx 2 cosxuxy sin2xuyy= 0:
(c)yuxx+ 2(x 1)uxy (y+ 2)uyy= 0:
Exercise 12.12
Letc>0:By computing ux;uxx;ut;anduttshow that
u(x;t) =1
2(f(x+ct) +f(x ct)) +1
2cZx+ct
x ctg(s)ds
is a solution to the PDE
utt=c2uxx
wherefis twice dierentiable function and gis a dierentiable function.
Then compute and simplify u(x;0) andut(x;0):
Exercise 12.13
Consider the second-order PDE
yuxx+uxy x2uyy ux u= 0:
Determine the region DinR2;if such a region exists, that makes this PDE:
(a) hyperbolic, (b) parabolic, (c) elliptic.
Exercise 12.14
Consider the second-order hyperbolic PDE
uxx+ 2uxy 3uyy= 0:
Use the change of variables v(x;y) =y 3xandw(x;y) =x+yto solve the
given equation.
Exercise 12.15
Solve the Cauchy problem
uxx+ 2uxy 3uyy= 0:
u(x;2x) = 1; ux(x;2x) =x:
104SECOND ORDER LINEAR PARTIAL DIFFERENTIAL EQUATIONS
13 Hyperbolic Type: The Wave equation
The wave equation has many physical applications from sound waves in air
to magnetic waves in the Sun's atmosphere. However, the simplest systems
to visualize and describe are waves on a stretched elastic string.
Initially the string is horizontal with two xed ends say a left end Land a
right endR:Then from end Lwe shake the string and we notice a wave
propogate through the string. The aim is to try and determine the vertical
displacement from the x axis of the string, u(x;t);as a function of position
xand timet:A displacement of a tiny piece of the string between points P
andQis shown in Figure 13.1.
Figure 13.1
where
(x;t) is the angle between the string and a horizontal line at position x
and timet;
T(x;t) is the tension in the string at position xand timet;
(x) is the mass density of the string at position x:
To derive the wave equation we need to make some simplifying assumptions:
(1) The density of the string, ;is constant so that the mass of the string
betweenPandQis simplytimes the length of the string between Pand
13 HYPERBOLIC TYPE: THE WAVE EQUATION 105
Q, where the length of the string is sgiven by
s=p
(x)2+ (u)2= xs
1 +u
x2
xs
1 +@u
@x2
(2) The displacement, u(x;t), and its derivatives are assumed small so that
sx
and the mass of the portion of the string is
x:
(3) The only forces acting on this portion of the string are the tensions
T(x;t) atPandT(x+ x;t) atQ:(In physics, tension is the magnitude of
the pulling force exerted by a string). The gravitational force is neglected.
(4) Our tiny string element moves only vertically. Then the net horizontal
force on it must be zero.
Next, we consider the forces acting on the typical string portion shown in
Figure 13.1. These forces are:
(i) tension pulling to the right, which has magnitude T(x+ x;t);and acts
at an angle (x+ x;t) above the horizontal.
(ii) tension pulling to the left, which has magnitude T(x;t);and acts at an
angle(x;t) above the horizontal.
Now we resolve the forces into their horizontal and vertical components.
Horizontal: The net horizontal force of the tiny string is T(x+x;t) cos(x+ x;t)
T(x;t) cos(x;t):Since there is no horizontal motion, we must have
T(x;t) cos(x;t) =T(x+ x;t) cos(x+ x;t) =T: (13.1)
Vertical: AtPthe tension force is T(x;t) sin(x;t) whereas at Qthe
force isT(x+ x;t) sin(x+ x;t). Then Newton's Law of motion
massacceleration = Applied Forces
gives
x@2u
@t2=T(x+ x;t) sin(x+ x;t) T(x;t) sin(x;t):
106SECOND ORDER LINEAR PARTIAL DIFFERENTIAL EQUATIONS
Dividing by Tand using (13.1) we obtain
Tx@2u
@t2=T(x+ x;t) sin(x+ x;t)
T(x+ x;t) cos(x+ x;t) T(x;t) sin(x;t)
T(x;t) cos(x;t)
= tan(x+ x;t) tan(x;t):
But
tan(x;t) = lim
x!0u
x=ux(x;t):
Likewise,
tan(x+ x;t) =ux(x+ x;t):
Hence, we get
Txutt(x;t) =ux(x+ x;t) ux(x;t):
Dividing by xand letting x!0 we obtain
Tutt(x;t) =uxx(x;t)
or
utt(x;t) =c2uxx(x;t) (13.2)
wherec2=T
:We callcthewave speed.
D'Alembert Solution of (13.2)
Letv=x+ctandw=x ct:Then by application of the chain rule we nd
ut=c(uv uw)
ux=uv+uw
utt=c2(uvv 2uvw+uww)
uxx=uvv+ 2uvw+uww:
Substituting into (13.2) we obtain
c2(uvv+ 2uvw+uww) =c2(uvv 2uvw+uww)
and this simplies to
4c2uvw= 0 oruvw= 0:
13 HYPERBOLIC TYPE: THE WAVE EQUATION 107
It follows that
u(v;w) =f(v) +g(w)
wherefandgare arbitrary dierentiable functions. Now, writing uin terms
ofxandywe nd the general solution
u(x;y) =f(x+ct) +g(x ct):
D'Alembert's solution involves two arbitrary functions that are determined
(normally) by two initial conditions.
Example 13.1
Find the solution to the Cauchy problem
utt=c2uxx
u(x;0) =v(x)
ut(x;0) =w(x):
Solution.
We have
u(x;0) =f(x) +g(x) =v(x)
and
ut(x;0) =cf0(x) cg0(x) =w(x)
which implies that
f(x) g(x) =1
cW(x) =1
cZ
w(x)dx:
Therefore,
g(x) =1
2(v(x) 1
cW(x)):
Hence,
f(x) =1
2(v(x) +1
cW(x)):
Finally,
u(x;t) =1
2[v(x ct) +v(x+ct) +1
c(W(x+ct) W(x ct))]
=1
2[v(x ct) +v(x+ct) +1
cZx+ct
x ctw(s)ds]
108SECOND ORDER LINEAR PARTIAL DIFFERENTIAL EQUATIONS
Practice Problems
Exercise 13.1
Show that if v(x;t) andw(x;t) satisfy equation (13.2) then v+wis also
a solution to (13.2), where andare constants.
Exercise 13.2
Show that any linear time independent function u(x;t) =ax+bis a solution
to equation (13.2).
Exercise 13.3
Find a solution to (13.2) that satises the homogeneous conditions u(x;0) =
u(0;t) =u(L;t) = 0:
Exercise 13.4
Solve the initial value problem
utt=9uxx
u(x;0) = cosx
ut(x;0) =0:
Exercise 13.5
Solve the initial value problem
utt=uxx
u(x;0) =1
1 +x2
ut(x;0) =0:
Exercise 13.6
Solve the initial value problem
utt=4uxx
u(x;0) =1
ut(x;0) = cos (2x):
13 HYPERBOLIC TYPE: THE WAVE EQUATION 109
Exercise 13.7
Solve the initial value problem
utt=25uxx
u(x;0) =v(x)
ut(x;0) =0
where
v(x) =1 ifx<0
0 ifx0:
Exercise 13.8
Solve the initial value problem
utt=c2uxx
u(x;0) =e x2
ut(x;0) = cos2x:
Exercise 13.9
Prove that the wave equation, utt=c2uxxsatises the following properties,
which are known as invariance properties. If u(x;t) is a solution, then
(i) Any translate, u(x y;t) whereyis a xed constant, is also a solution.
(ii) Any derivative, say ux(x;t);is also a solution.
(iii) Any dilation, u(ax;at );is a solution, for any xed constant a.
Exercise 13.10
Findv(r) ifu(r;t) =v(r)
rcosntis a solution to the PDE
urr+2
rur=utt:
110SECOND ORDER LINEAR PARTIAL DIFFERENTIAL EQUATIONS
Sample Exam Questions
Exercise 13.11
Find the solution of the wave equation on the real line ( 1< x < +1)
with the initial conditions
u(x;0) =ex; ut(x;0) = sinx:
Exercise 13.12
The total energy of the string (the sum of the kinetic and potential energies)
is dened as
E(t) =1
2ZL
0(u2
t+c2u2
x)dx:
(a) Using the wave equation derive the equation of conservation of energy
dE(t)
dt=c2(ut(L;t)ux(L;t) ut(0;t)ux(0;t)):
(b) Assuming xed ends boundary conditions, that is the ends of the string
are xed so that u(0;t) =u(L;t) = 0;for allt>0;show that the energy is
constant.
(c) Assuming free ends boundary conditions for both x= 0 andx=L;that
is bothu(0;t) andu(L;t) vary with t;show that the energy is constant.
Exercise 13.13
For a wave equation with damping
utt c2uxx+dut= 0; d> 0;0<x<L
with the xed ends boundary conditions show that the total energy decreases.
Exercise 13.14
(a) Verify that for any twice dierentiable R(x) the function
u(x;t) =R(x ct)
is a solution of the wave equation utt=c2uxx:Such solutions are called
traveling waves.
(b) Show that the potential and kinetic energies (see Exercise 13.12) are
equal for the traveling wave solution in (a).
13 HYPERBOLIC TYPE: THE WAVE EQUATION 111
Exercise 13.15
Find the solution of the Cauchy wave equation
utt= 4uxx
u(x;0) =x2; ut(x;0) = sin 2x:
Simplify your answer as much as possible.
112SECOND ORDER LINEAR PARTIAL DIFFERENTIAL EQUATIONS
14 Parabolic Type: The Heat Equation in One-
Dimensional Space
In this section, We will look at a model for describing the distribution of
temperature in a solid material as a function of time and space.
Before we begin our discussion of the mathematics of the heat equation, we
must rst determine what is meant by the term heat ? Heat is type of energy
known as thermal energy . Heat travels in waves like other forms of energy,
and can change the matter it touches. It can heat it up and cause chemical
reactions like burning to occur.
Heat can be released through a chemical reaction (such as the nuclear re-
actions that make the Sun \burn") or can be trapped for a limited time by
insulators. It is often released along with other kinds of energy such as light
waves or sound waves. For example, a burning candle releases light and heat
waves. On the other hand, an explosion releases light, heat, and sound waves.
The most common units of heat are BTU (British Thermal Unit), Calorie
and Joule.
Consider now a rod made of homogeneous heat conducting material (i.e. it
is composed of the exact same material and no foreign bodies are in it) of
uniform density and constant cross section A;placed along the x axis from
x= 0 tox=Las shown in Figure 14.1.
Figure 14.1
Assume the heat
ows only in the x direction, with the lateral sides well
insulated, and the only way heat can enter or leave the rod is at either end.
Also we assume that the temperature of the rod is constant at any point of
the cross section. In other words, temperature will only vary in xand we
can hence consider the rod to be a one spatial dimensional rod. We will also
assume that heat energy in any piece of the rod is conserved.
Letu(x;t) be the temperature of the cross section at the point xand the
timet:Consider a portion Uof the rod from xtox+ xof length xas
14 PARABOLIC TYPE: THE HEAT EQUATION IN ONE-DIMENSIONAL SPACE 113
shown in Figure 14.2.
Figure 14.2
Consider the portion SofUof height s:From the theory of heat conduction,
the quantity of heat Qfromxtox+ sat timetis given by
Q=cu(x;t)V
where Vis the volume of Sandcis the specic heat , that is, the amount
of heat energy that it takes to raise one unit of mass of the material by one
unit of temperature.
ButSis a cylinder of height sand area of base Aso that V=As:
Hence,
Q=cAu (x;t)s:
The quantity of heat in the portion Uis given by
Q(t) =Zx+x
xcAu (s;t)ds:
By dierentiating we take the partial of uto nd the change in heat with
respect to time.
dQ
dt=Zx+x
xcAut(s;t)ds:
Assuming that uis continuously dierentiable, we can apply the mean value
theorem for integrals and nd xx+ xsuch that
Zx+x
xut(s;t)ds= xut(;t):
Thus, the rate of change of heat in Uis given by
dQ
dt=cAxut(;t):
114SECOND ORDER LINEAR PARTIAL DIFFERENTIAL EQUATIONS
On the other hand, by Fourier (or Fick's) law of heat conduction, the rate
of heat
ow through any cross section is proportional to the area Aand
the negative gradient of the temperature normal to the cross section, and
heat
ows in the direction of decreasing temperature. Thus, the rate of heat
owing inUthrough the cross section at xis KAux(x;t) and the rate of
heat
owing out of Uthrough the cross section at x+xis KAux(x+x;t);
whereKis the thermal conductivity of the rod.
Now, the conservation of energy law states
rate of change of heat in U= rate of heat
owing in rate of heat
owing
out
or mathematically written as,
cAxut(;t) = KAux(x;t) +KAux(x+ x;t)
or
cAxut(;t) =KA[ux(x+ x;t) ux(x;t)]:
Dividing this last equation by cAxand letting x!0 we obtain
ut(x;t) =kuxx(x;t) (14.1)
wherek=K
cis called the diusivity constant.
Equation (14.1) is the one dimensional heat equation which is second order,
linear, homogeneous, and of parabolic type.
The non-homogeneous heat equation
ut=kuxx+f(x)
is known as the heat equation with an external heat source f(x):An ex-
ample of an exterenal heat source is the heat generated from a candle placed
under the bar.
The function
E(t) =ZL
0cu(x;t)dx
is called the total thermal energy at timetof the entire rod.
Example 14.1
The two ends of a uniform rod of length Lare insulated. There is a con-
stant source of thermal energy q06= 0 and the temperature is initially
14 PARABOLIC TYPE: THE HEAT EQUATION IN ONE-DIMENSIONAL SPACE 115
u(x;0) =f(x):
(a) Write the equation and the boundary conditions for this model.
(b) Calculate the total thermal energy of the entire rod.
Solution.
(a) The model is given by the PDE
cut(x;t) =Kuxx+q0
with boundary conditions
ux(0;t) =ux(L;t) = 0:
(b) First note that
d
dtZL
0cu(x;t)dx=ZL
0cut(x;t)dx=ZL
0Kuxxdx+ZL
0q0dx
=KuxjL
0+q0L=q0L
sinceux(0;t) =ux(L;t) = 0:Integrating in time from 0 to twe nd
E(t) =q0Lt+C:
ButC=E(0) =RL
0cu(x;0)dx=RL
0cf(x)dx:Hence, the total thermal
energy is given by
E(t) =ZL
0cf(x)dx+q0Lt
Initial Boundary Value Problems
In order to solve the heat equation we must give the problem some initial
conditions. If you recall from the theory of ODE, the number of conditions
required for solving initial value problems always matched the highest order
of the derivative in the equation.
In partial dierential equations the same idea holds except now we have to
pay attention to the variable we are dierentiating with respect to as well.
So, for the heat equation we have got a rst order time derivative and so we
will need one initial condition and a second order spatial derivative and so
we will need two boundary conditions.
116SECOND ORDER LINEAR PARTIAL DIFFERENTIAL EQUATIONS
For the initial condition, we dene the temperature of every point along the
rod at time t= 0 by
u(x;0) =f(x)
wherefis a given (prescribed) function of x:This function is known as the
initial temperature distribution.
The boundary conditions will tell us something about what the temperature
is doing at the ends of the bar. The conditions are given by
u(0;t) =T0andu(L;t) =TL:
and they are called as the Dirichlet conditions. In this case, the general
form of the heat equation initial boundary value problem is to nd u(x;t)
satisfying
ut(x;t) =kuxx(x;t);0xL; t> 0
u(x;0) =f(x);0xL
u(0;t) =T0;u(L;t) =TL; t> 0:
In the case of insulated endpoints, i.e. there is no heat
ow out of them, we
use the boundary conditions
ux(0;t) =ux(L;t) = 0:
These conditions are examples of what is known as Neumann boundary
conditions. In this case, the general form of the heat equation initial bound-
ary value problem is to nd u(x;t) satisfying
ut(x;t) =kuxx(x;t);0xL; t> 0
u(x;0) =f(x);0xL
ux(0;t) =ux(L;t) = 0; t> 0:
14 PARABOLIC TYPE: THE HEAT EQUATION IN ONE-DIMENSIONAL SPACE 117
Practice Problems
Exercise 14.1
Show that if u(x;t) andv(x;t) satisfy equation (14.1) then u+vis also a
solution to (14.1), where andare constants.
Exercise 14.2
Show that any linear time independent function u(x;t) =ax+bis a solution
to equation (14.1).
Exercise 14.3
Find a linear time independent solution uto (14.1) that satises u(0;t) =T0
andu(L;T) =TL:
Exercise 14.4
Show that to solve (14.1) with the boundary conditions u(0;t) =T0and
u(L;t) =TLit suces to solve (14.1) with the homogeneous boundary
conditions u(0;t) =u(L;t) = 0:
Exercise 14.5
Find a solution to (14.1) that satises the conditions u(x;0) =u(0;t) =
u(L;t) = 0:
Exercise 14.6
Let (I) denote equation (14.1) together with intial condition u(x;0) =f(x),
wherefis not the zero function, and the homogeneous boundary conditions
u(0;t) =u(L;t) = 0:Suppose a nontrivial solution to (I) can be written in
the formu(x;t) =X(x)T(t):Show thatXandTsatisfy the ODE
X00
kX= 0 andT0 T= 0
for some constant :
Exercise 14.7
Consider again the solution u(x;t) =X(x)T(t):Clearly,T(t) =T(0)et:
Suppose that >0:
(a) Show that X(x) =Aexp+Be xp;where=
kandAandBare
arbitrary constants.
(b) Show that AandBsatisfy the two equations A+B= 0 andA(eLp
e Lp) = 0:
(c) Show that A= 0 leads to a contradiction.
(d) Using (b) and (c) show that eLp=e Lp:Show that this equality leads
to a contradiction. We conclude that <0:
118SECOND ORDER LINEAR PARTIAL DIFFERENTIAL EQUATIONS
Exercise 14.8
Consider the results of the previous exercise.
(a) Show that X(x) =c1cosx+c2sinxwhere=q
k:
(b) Show that =n= kn22
L2;wherenis an integer.
Exercise 14.9
Show thatu(x;t) =Pn
k=1uk(x;t);whereun(x;t) =cne kn22
L2tsin n
L
xsat-
ises (14.1) and the homogeneous boundary conditions.
Exercise 14.10
Suppose that a wire is stretched between 0 and a:Describe the boundary
conditions for the temperature u(x;t) when
(i) the left end is kept at 0 degrees and the right end is kept at 100 degrees;
and
(ii) when both ends are insulated.
Exercise 14.11
Letut=uxxfor 0< x < andt >0 with boundary conditions u(0;t) =
0 =u(;t) and initial condition u(x;0) =f(x):LetE(t) =R
0(u2
t+u2
x)dx:
Show thatE0(t)<0:
Exercise 14.12
Suppose
ut=uxx+ 4; ux(0;t) = 5; ux(L;t) = 6; u(x;0) =f(x):
Calculate the total thermal energy of the one-dimensional rod (as a function
of time).
14 PARABOLIC TYPE: THE HEAT EQUATION IN ONE-DIMENSIONAL SPACE 119
Sample Exam Questions
Exercise 14.13
Consider the heat equation
ut=kuxx
forx2(0;1) andt>0;with boundary conditions u(0;t) = 2 andu(1;t) = 3
fort >0 and initial conditions u(x;0) =xforx2(0;1):A function v(x)
that satises the equation v00(x) = 0;with conditions v(0) = 2 and v(1) = 3
is called a steady-state solution. That is, the steady-state solutions of the
heat equation are those solutions that don't depend on time. Find v(x):
Exercise 14.14
Consider the equation for the one-dimensional rod of length Lwith given
heat energy source:
ut=uxx+q(x):
Assume that the initial temperature distribution is given by u(x;0) =f(x):
Find the equilibrium (steady state) temperature distribution in the following
cases.
(a)q(x) = 0;u(0) = 0;u(L) =T:
(b)q(x) = 0;ux(0) = 0;u(L) =T:
(c)q(x) = 0;u(0) =T;ux(L) =:
Exercise 14.15
Consider the equation for the one-dimensional rod of length Lwith insulated
ends:
cut=Kuxx; ux(0;t) =ux(L;t) = 0:
(a) Give the expression for the total thermal energy of the rod.
(b) Show using the equation and the boundary conditions that the total
thermal energy is constant.
Exercise 14.16
Suppose
ut=uxx+x; u (x;0) =f(x); ux(0;t) =; ux(L;t) = 7:
(a) Calculate the total thermal energy of the one-dimensional rod (as a func-
tion of time).
(b) From part (a) nd the value of for which a steady-state solution exist.
(c) For the above value of nd the steady state solution.
120SECOND ORDER LINEAR PARTIAL DIFFERENTIAL EQUATIONS
15 An Introduction to Fourier Series
In this and the next section we will have a brief look to the subject of Fourier
series. The point here is to do just enough to allow us to do some basic so-
lutions to partial dierential equations later in the book.
Motivation: In Calculus we have seen that certain functions may be repre-
sented as power series by means of the Taylor expansions. These functions
must have innitely many derivatives, and the series provide a good approx-
imation only in some (often small) vicinity of a reference point.
Fourier series constructed of trigonometric rather than power functions, and
can be used for functions not only not dierentiable, but even discontinuous
at some points. The main limitation of Fourier series is that the underlying
function should be periodic.
Recall from calculus that a function series is a series where the summands
are functions. Examples of function series include power series, Laurent se-
ries, Fourier series, etc.
Unlike series of numbers, there exist many types of convergence of series of
functions, namely, pointwise, uniform, etc. We say that a series of functionsP1
n=1fn(x)converges pointwise to a function fif and only if the sequence
of partial sums
Sn(x) =f1(x) +f2(x) ++fn(x)
converges pointwise to f:We write
1X
n=1fn(x) = lim
n!1Sn(x) =f(x):
Likewise, we say that a series of functionsP1
n=1fn(x)converges uniformly
to a function fif and only if the sequence of partial sums fSng1
n=1converges
uniformly to f:
In this section we introduce a type of series of functions known as Fourier
series . They are given by
f(x) =a0
2+1X
n=1h
ancosn
Lx
+bnsinn
Lxi
; LxL(15.1)
whereanandbnare called the Fourier coecients. The expression on the
right is called a trigonometric series. Note that we begin the series with
a0
2as opposed to simply a0to simplify the coecient formula for anthat we
15 AN INTRODUCTION TO FOURIER SERIES 121
will derive later in this section.
The main questions we want to consider next are the questions of determin-
ing which functions can be represented by Fourier series and if so how to
compute the coecients anandbn:
Before answering these questions, we look at some of the properties of Fourier
series.
Periodicity Property
Recall that a function fis said to be periodic with period T > 0 if
f(x+T) =f(x) for allx;x+Tin the domain of f:The smallest value
ofTfor whichfis periodic is called the fundamental period. A graph of
aT periodic function is shown in Figure 15.1.
Figure 15.1
For aT periodic function we have
f(x) =f(x+T) =f(x+ 2T) =:
Note that the denite integral of a T periodic function is the same over any
interval of length T:By Exercise 15.1 below, if fandgare two periodic func-
tions with common period T;then the product fgand an arbitrary linear
combination c1f+c2gare also periodic with period T:It is an easy exercise
to show that the Fourier series (15.1) is periodic with fundamental period 2 L:
Orthogonality Property
Recall from Calculus that for each pair of vectors ~ uand~ vwe associate a
scalar quantity ~ u~ vcalled the dot product of~ uand~ v:We say that ~ uand~ v
areorthogonal if and only if ~ u~ v= 0:We want to dene a similar concept
for functions.
Letfandgbe two functions with domain the closed interval [ a;b]:We dene
122SECOND ORDER LINEAR PARTIAL DIFFERENTIAL EQUATIONS
a function that takes a pair of functions to a scalar. Symbolically, we write
<f;g> =Zb
af(x)g(x)dx:
We call< f;g > theinner product offandg:We say that fandg
areorthogonal if and only if < f;g > = 0:A set of functions is said to
bemutually orthogonal if each distinct pair of functions in the set is
orthogonal.
Example 15.1
Show that the set
1;cos n
Lx
;sin n
Lx
:n2N
is mutually orthogonal in
[ L;L]:
Solution.
We have ZL
L1cosn
Lx
dx=L
nh
sinn
LxiL
L= 0
andZL
L1sinn
Lx
dx= L
nh
cosn
LxiL
L= 0:
Now, forn6=mwe have
ZL
Lcosm
Lx
cosn
Lx
dx=1
2ZL
L
cos(m+n)
Lx
+ cos(m n)
Lx
dx
=1
2L
(m+n)sin(m+n)
Lx
+L
(m n)sin(m n)
LxL
L= 0
where we used the trigonometric identity
cosacosb=1
2[cos (a+b) + cos (a b)]:
In the exercises below, we show that
ZL
Lsinm
Lx
sinn
Lx
dx= 0
15 AN INTRODUCTION TO FOURIER SERIES 123
and ZL
Lcosm
Lx
sinn
Lx
dx= 0
The reason we care about these functions being orthogonal is because we will
exploit this fact to develop a formula for the coecients in our Fourier series.
Now, in order to answer the rst question mentioned earlier, that is, which
functions can be expressed as a Fourier series expansion, we need to intro-
duce some mathematical concepts.
A function f(x) is said to be piecewise continuous on [a;b] if it is contin-
uous in [a;b] execept possibly at nitely many points of discontinuity within
the interval [ a;b];and at each point of discontinuity, the right- and left-
handed limits of fexist. An example of a piecewise continuous function is
the function
f(x) =x 0x<1
x2 x1x2:
We will say that fispiecewise smooth in [a;b] if and only if f(x) as well
as its derivatives are piecewise continuous.
The following theorem, proven in more advanced books, ensures that a
Fourier decomposition can be found for any function which is piecewise
smooth.
Theorem 15.1
Letfbe a 2L-periodic function. If fis a piecewise smooth on [ L;L] then
for all points of discontinuity x2( L;L) we have
f(x ) +f(x+)
2=a0
2+1X
n=1h
ancosn
Lx
+bnsinn
Lxi
:
where as for points of continuity x2( L;L) we have
f(x) =a0
2+1X
n=1h
ancosn
Lx
+bnsinn
Lxi
:
Remark 15.1
(1) Almost all functions occurring in practice are piecewise smooth functions.
(2) Given a non-periodic function fon [ L;L]. The above theorem applies
to the periodic extension FoffwhereF(x+ 2nL) =f(x) (n2Z) and
F(x) =f(x) on [ L;L]:
124SECOND ORDER LINEAR PARTIAL DIFFERENTIAL EQUATIONS
Convergence Results of Fourier Series
We list few of the results regarding the convergence of Fourier series:
(1) The type of convergence in the above theorem is pointwise convergence.
(2) The convergence is uniform for a continuous function fon [ L;L] such
thatf( L) =f(L):
(3) The convergence is uniform wheneverP1
n=1(janj2+jbnj2) is convergent.
(4) Iff(x) is periodic, continuous, and has a piecewise continuous derivative,
then the Fourier Series corresponding to fconverges uniformly to f(x) for
the entire real line.
(5) The convergence is uniform on any closed interval that does not contain
a point of discontinuity.
Euler-Fourier Formulas
Next, we will answer the second question mentioned earlier, that is, the ques-
tion of nding formulas for the coecients anandbn:These formulas for an
andbnare called Euler-Fourier formulas which we derive next. We will as-
sume that the RHS in (15.1) converges uniformly to f(x) on the interval
[ L;L]:Integrating both sides of (15.1) we obtain
ZL
Lf(x)dx=ZL
La0
2dx+ZL
L1X
n=1h
ancosn
Lx
+bnsinn
Lxi
dx:
Since the trigonometric series is assumed to be uniformly convergent, from
Section 2, we can interchange the order of integration and summation to
obtain
ZL
Lf(x)dx=ZL
La0
2dx+1X
n=1ZL
Lh
ancosn
Lx
+bnsinn
Lxi
dx:
But ZL
Lcosn
Lx
dx=L
nsinn
LxiL
L= 0
and likewise
ZL
Lsinn
Lx
dx= L
ncosn
LxiL
L= 0:
Thus,
a0=1
LZL
Lf(x)dx:
15 AN INTRODUCTION TO FOURIER SERIES 125
To nd the other Fourier coecients, we recall the results of Exercises 15.2
- 15.3 below.
ZL
Lcosn
Lx
cosm
Lx
dx=Lifm=n
0 ifm6=n
ZL
Lsinn
Lx
sinm
Lx
dx=Lifm=n
0 ifm6=n
ZL
Lsinn
Lx
cosm
Lx
dx= 0;8m;n:
Now, to nd the formula for the Fourier coecients amform> 0;we multiply
both sides of (15.1) by cos m
Lx
and integrate from LtoLto otbain
ZL
Lf(x) cosm
Lx
=ZL
La0
2cosm
Lx
dx+1X
n=1
anZL
Lcosn
Lx
cosm
Lx
dx
+bnZL
Lsinn
Lx
cosm
Lx
dx:
Hence,ZL
Lf(x) cosm
Lx
dx=amL
and therefore
am=1
LZL
Lf(x) cosm
Lx
dx:
Likewise, we can show that
bm=1
LZL
Lf(x) sinm
Lx
dx:
Example 15.2
Find the Fourier series expansion of
f(x) =0; x0
x; x> 0
on the interval [ ;]:
126SECOND ORDER LINEAR PARTIAL DIFFERENTIAL EQUATIONS
Solution.
We have
a0=1
Z
f(x)dx=1
Z
0xdx=
2
an=1
Z
0xcosnxdx =1
xsinnx
n+cosnx
n2
0=( 1)n 1
n2
bn=1
Z
0xsinnxdx =1
xcosnx
n+sinnx
n2
0=( 1)n+1
n
Hence,
f(x) =
4+1X
n=1( 1)n 1
n2cos (nx) +( 1)n+1
nsin (nx)
Example 15.3
Apply Theorem 15.1 to the function in Example 15.2.
Solution.
LetFbe a periodic extension of fof period 2:Thus,f(x) =F(x) on the
interval [ ;]:Clearly,Fis a piecewise smooth function so that by the
previous thereom we can write
4+1X
n=1( 1)n 1
n2cos (nx) +( 1)n+1
nsin (nx)
=8
<
:
2; ifx=
f(x);if <x<
2; ifx=
Takingx=we have the identity
4+1X
n=1( 1)n 1
n2( 1)n=
2
which can be simplied to
1X
n=11
(2n 1)2=2
8:
This provides a method for computing an approximate value of
15 AN INTRODUCTION TO FOURIER SERIES 127
Remark 15.2
An example of a function that does not have a Fourier series representation
is the function f(x) =1
x2on [ L;L]:For example, the coecient a0for this
function does not exist. Thus, not every function can be written as a Fourier
series expansion.
The nal topic of discussion here is the topic of dierentiation and integration
of Fourier series. In particular we want to know if we can dierentiate a
Fourier series term by term and have the result be the Fourier series of the
derivative of the function. Likewise we want to know if we can integrate a
Fourier series term by term and arrive at the Fourier series of the integral of
the function. Answers to these questions are provided next.
Theorem 15.2
A Fourier series of a piecewise smooth functionfcan always be integrated
term by term and the result is a convergent innite series that always con-
verges toRL
Lf(x)dxeven if the original series has jumps.
Theorem 15.3
A Fourier series of a continuous function f(x) can be dierentiated term by
term iff0(x) ispiecewise smooth. The result of the dierentiation is the
Fourier series of f0(x):
128SECOND ORDER LINEAR PARTIAL DIFFERENTIAL EQUATIONS
Practice Problems
Exercise 15.1
Letfandgbe two functions with common domain Dand common period
T:Show that
(a)fgis periodic of period T:
(b)c1f+c2gis periodic of period T;wherec1andc2are real numbers.
Exercise 15.2
Show that for m6=nwe have
(a)RL
Lsin m
Lx
sin n
Lx
dx= 0 and
(b)RL
Lcos m
Lx
sin n
Lx
dx= 0:
Exercise 15.3
Compute the following integrals:
(a)RL
Lcos2 n
Lx
dx:
(b)RL
Lsin2 n
Lx
dx:
(c)RL
Lcos n
Lx
sin n
Lx
dx:
Exercise 15.4
Find the Fourier coecients of
f(x) =8
<
: ; x<0
; 0<x<
0; x = 0;
on the interval [ ;]:
Exercise 15.5
Find the Fourier series of f(x) =x2 1
2on the interval [ 1;1]:
Exercise 15.6
Find the Fourier series of the function
f(x) =8
<
: 1; 2<x<