Pde Final PDF
Pde Final PDF
Pde Final PDF
EQUATIONS
January 4, 1980
Course Outline
1. Surfaces and Curves in three dimensions
3. -Methods of solution of dx
P
= dy
Q
= dz
R
8. Applications to Engineering
Reference Books
1. Elements of Partial dierential equations by I,N, Sneddon, Mc Graw
Hill(1957)
1
2. p.d.e-Schaum Outline Series
3. p.d.e-Garabedian (1964)
2
Contents
TABLE OF CONTENTS 3
1 INTRODUCTION 5
2 Partial Dierential Equation 6
2.1 Denition . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6
2.2 Formation of Partial Dierential Equations. . . . . . . . . . . 6
2.2.1 Formation of P.D.E by eliminating arbitrary constant . 7
2.2.2 Formation of P.D.E by eliminating arbitrary function . 9
2.3 Derivation of P.D.E equations by eliminating arbitrary func-
tions φ (u, v) = 0 . . . . . . . . . . . . . . . . . . . . . . . . . 10
3
10 HEAT EQUATION 62
11 DERIVATION OF THE WAVE EQUATION 72
12 DERIVATION OF HEAT EQUATION 74
13 OTHER BOUNDARY VALUE PROBLEMS 76
14 METHOD OF LAPLACE TRANSFORMS AND FOURIER
TRANSFORMS: 77
14.1 SINE AND COSINE FOURIER TRANSFORMS . . . . . . . 78
4
1 INTRODUCTION
Partial dierential equations are used to formulate problems involving func-
tions of several variables and are either solved by hand, or used to create a
relevant computer model.
Denition
A dierential equation is an equation that relates the derivatives at a scalar
function depending on one or more variables.
Example:
δ4 u
+ δδxu2 +u2 = cos(x)is a dierential equation of the function u(x)depending
2
δx4
on x, while δu = δδxu2 + δδyu2 − uis a dierential equation involving a function
2 2
δt
u(x, y, t) depending on x, y, t. This equation can be written in subscript
notation as ut = uxx + uyy − u.
Partial dierential operators
(i) δxδ = δx
(ii) δxδ 2 = δx2 or δxx
2
5
In R2 the equation x = 3tells us to graph all the points that are in the form
(3, y).This is a vertical line in a 2-D coordinate system
In R3 the equation x = 3tells us to graph all the points onf the form (3, y, z).This
is similar to the coordinates for the yz plane except this time we have
x = 3instead of x = 0. So in a 3-D system this is a plane that will be
parallel to the yz plane and pass through the x-axis at x = 3.
6
2.2.1 Formation of P.D.E by eliminating arbitrary constant
Example
(i) z = ax + by + a2 + b2 .............(i)
Solution
δz
= a.................................(ii)
δx
δz
= b.................................(iii)
δy
Replacinng (ii) and (iii) in equation (i), wehave:
2 2
δz δz δz δz
z= x+ y+ +
δx δy δx δy
z = px + qy + p2 + q 2
7
δz
= Apept sin(px)
δt
δ2z
= Ap2 ept sin(px)......................(ii)
δt2
Ading, (i) and (ii), we have:
δ2z δ2z
2
+ 2 = −Ap2 ept sin(px) + Ap2 ept sin(px) = 0
δx δt
b) xa2 + yb2 + zc2 = 1..................(i)
2 2 2
Solution
Dierentiting (i) with respect to x,
2x 2z δz δz
2
+ 2 = 0 =⇒ c2 x + a2 z = 0..............(ii)
a c δx δx
Dierentiting (i) with respect to y,
2y 2z δz δz
2
+ 2 = 0 =⇒ c2 y + b2 z = 0..............(iii)
b c δy δy
Dierentiting (ii) with respect to x,
2
δ2z
2 2 δz 2
c +a + a z 2 = 0..................(iv)
δx δx
Dierentiting (iii) with respect to y,
2
δ2z
2 2 δz
c +b + b2 z = 0..................(v)
δy δy 2
From (ii)
a2 z
2 δz
c =− .....................(vi)
x δx
Putting (vi) into(iv) 2
a2 z δ2z
δz 2 δz 2
− +a + a z 2 = 0..................(vii)
x δx δx δx
Dividing (vii) by a2 andmultiplying by x
2
δ2z
δz δz
zx 2 + x −z = 0....................(viii)
δx δx δx
Similarly from (iii),
8
b2 z
2 δz
c =− .....................(ix)
y δy
Putting (ix) into (v) 2
2
δ2z
bz δz 2 δz 2
− +b + b z 2 = 0..................(x)
y δy δy δy
Dividing (x) by a2 and multiplying by y
2
δ2z
δz δz
zy 2 + y −z = 0....................(xi)
δy δy δy
δ2z 2
00 2 00 2
0 2 0 2
= 4x f x − y + g x + y +2 g x + y + f x + y ........(iii)
δx2
9
δz
= −f 0 x2 − y + g 0 x2 + y .......(iv)
δy
δ2z 00 2
00 2
= f x − y + g x + y .............(v)
δy 2
1 δz
= f 0 x2 − y + g 0 x2 + y ............(vi)
2x δx
Substituting (vi) and (v) into (iii),
δ2z 2
2 δ z 1 δz
= 4x +2 ........(vii)
δx2 δy 2 2x δx
δ2z
2
3 δ z δz
x 2 = 4x 2
+ ........(viii)
δx δy δx
10
δv δv δu δu δv δv δu δu
+p +q = +q +p
δx δz δy δz δy δz δx δz
δv δu δv δu δv δu δv δu δv δu δv δu δv δu δv δu
· +q · +p · +pq = · +p · +q · +pq ·
δx δy δx δz δz δy δz δz δy δx δy δz δz δx δz δz
δv δu δv δu δv δu δv δu δv δu δv δu δv δu δv δu
p · − · +q · − · = − · −pq + · +pq ·
δz δy δy δz δx δz δz δx δx δy δz δz δy δx δz δz
δv δu δv δu δv δu δv δu δv δu δv δu
p · − · +q · − · = · − ·
δz δy δy δz δx δz δz δx δy δx δx δy
P p + Qq = R
Example
a) Form a partial dierential equation whose solution is φ (x2 ez , yez ) = 0
Solution
φ x2 ez , yez = 0.................................(i)
0 + pyez ez + qyez
=
2xez + px2 ez 0 + qx2 ez
11
pyez ez + qyez
=
2xez + px2 ez qx2 ez
py 1 + qy
2
=
2x + px qx2
pyqx2 = 2x + px2 (1 + qy)
x2 p + 2xyq = −2x
xp + 2yq = −2
From this, p = x, q = 2y , R = −2
b) Form a p.d.e equation whose solution is φ (x + y + z, x2 + y 2 − z 2 ) = 0
Solution
Here u = x + y + z and v = x2 + y 2 − z 2
φ (u, v) = 0
δu
δx
= 1, δu
δz
= 1, δu
δy
δv
= 1, δx = 2x, δv
δz
δv
= −2z, δy = 2y
Substituting this in the equation,
δv
δv
δx
+ p δv
δz δy
+ q δv
δz
δu
=
δx
+ p δu
δz
δu
δy
+ q δu
δz
2x − 2pz 2y − 2qz
=
1+p 1+q
(2x − 2pz) (1 + q) = (1 + p) (2y − 2qz)
x + xq − zp = y − zq + py
(x + z) q − (z + y) p = y − x
wherep = −z − y , q = x + z , R = y − x
(b) Form a partial dierential equation whose solution is φ (x2 + y 2 + z 2 , z 2 − 2xy) =
0
12
Solution
φ x2 ez , yez = 0.................................(i)
z (y + x) p − z (y + x) q = (y + x) (y − x)
zp − zq = (y − x)
wherep = z , q = −z , R = y − x
13
3 LANGRAGE METHOD OF SOLVING THE
QUASI LINEAR PARTIAL DIFFERENTIAL
EQUATIONS OF ORDER 1 NAMELY; P p+
Qq = R
The general solution of the linear partial dierential equation are P p + Qq =
R where p,q, R are functions of x,y,z is φ (u, v) = 0where φis an arbitrary
function and u (x, y, z) = c1 )andv (x, y, z) = c2 ) form a solution of equations
dx dy dz
= =
p q R
Example:
1. Solve yxz p + x2 q = y 2
2
Solution
This dierential equation is of the form P p + Qq = R
wherep = yxz , q = xz , R = y 2
2
14
x2 z 2 c2
− =
2 2 2
x2 − z 2 = c2 =⇒ v
3x − y = c
y − 3x = c1
15
thus,
v (x, y, z) = φ(u)
Solution
The equation (x − y) p + (y − x − z) q = z is of the form P p + Qq = R
wherep = (x − y), q = (y − x − z), R = z
Using langrage auxillary equation,
dx dy dz
= = ......(i)
(x − y) (y − x − z) z
Using all the fractions
dx + dy + dz dx + dy + dz
=
x−y+y−x−z+z 0
dx + dy + dz = 0
x + y + z = c1 ...................(ii)
16
(2y − c1 ) = c2 z 2
2y − x − y − z
= c2
z2
y−x−z
= c2 .......................(iii)
z2
Put z =1 in equation (ii) and (iii)
x + y = c1 − 1...................(iv)
y − x = c2 + 1...................(v)
But,
2 x2 + y 2 = (x + y)2 + (y − x)2 ........................(vi)
But x2 + y 2 = 1
2 = (c1 − 1)2 + (c2 + 1)2 ........................(vii)
17
dx dz dx dz 1
= =⇒ = − =⇒ dx = −zdz
4yz −2y 2z 1 2
1 z2
x+ = c1
2 2
x + z 2 = c1 ..................(ii)
(x + z) + z 2 + y 2 = c1 + c2
2 + 1 = c1 + c2
3 = c1 + c2
(x + z) + z 2 + y 2 = 3
18
1 1
+ = c1
x z
Taking second and third fractions
dy dz 1 1
= =⇒ − − =c
y2 −z 2 y z
1 1
+ = c2
y z
Adding (i) and (ii)
1 1 2
+ + = c1 + c2
x y z
y+x 2
+ = c1 + c2
yx z
2 + 1 = c1 + c2
1 1 2
+ + =3
x y z
yz + xz + 2xy = 3xyz
ASSIGNMENT
Find the integral surface of equation 4yzp − q + 2y = 0which passes through
the hyperbola y 2 = 9 + z and z = 5 − x.
Solution
The dierential equation 4yzp − q + 2y = 0can be written in the form P p +
Qq = R
wherep = 4yz , q = −1, R = −2y
Using langrage auxillary equation,
dx dy dz
= = ......(i)
4yz −1 −2y
Taking rst and third fractions
dx dz dx dz 1 zdz 1 z2
= =⇒ = =⇒ dx = − =⇒ x + = c1
4yz −2y 2z −1 2 1 2 2
x + z 2 = c1 ...........................(ii)
19
dy dz dz y2 1
= =⇒ ydy = =⇒ − z = c2
−1 −2y 2 2 2
y 2 − z = c2 ..................................(iii)
9 + z 2 = z + c2
But z = 5 − xand z 2 = c1 − x
9 + c1 − x = 5 − x + c2
4 = c2 − c1
4 = y2 − z − x − z2
20
The dierential equation of the family of the surface is orthogonal to the
given family is given by pdx + qdy + Rdz = 0where (y + z) dx + (z + x) dy +
(x + y) dz = 0
Rearranging,
(y + z) dx + (z + x) dy + (x + y) dz = 0
Rearranging,
(ydx + xdy) + (zdy + ydz) + (xdz + zdx) = 0
Integrate,
xy + yz + xz = c
Solution
Solution
φ z (x + y)2 , x2 − y 2 = 0.................................(i)
2x −2y
2 2 =
(x + y) + p (x + y) 2z (x + y) + q (x + y)2
x [2z + q (x + y)] = −y [2z + p (x + y)]
21
2xz + qx2 + qxy = −2zy − xyp − y 2 p
xy + y 2 p + x2 + xy q = −2z (y + x)
py (x + y) + qx (x + y) = −2z (y + x)
py + qx = −2z
xy − z 2 = c
22
4 SPECIAL TYPES OF FIRST ORDER PAR-
TIAL DIFFERENTAL EQUATIONS
We discuss special types of rst order partial dierential equations whose
solution can be determined by the Charpit's method.
Case 1:
A rst order partial dierential equation is said to be of the form Charpit's
form, if it can be written in the form:
z = px + qy + f (p, q) ................(i)
Solution
Let F (x, y, z, p, q) = z − px − q(y) − f (p, q).................(ii)
Therefore the Charpit's auxillary equation is:
dp dq dz dx dy
δF
= = = δF = δF
δx
+ p δF
δZ
δF
δy
δF
+ q δZ δF δF
−p δp − q δq − δp − δq
Using the rst and the second fractions,
dp dq
=
−p + p −q + q
This implies that
dp = 0 =⇒ p = aand dq = 0 =⇒ q = b.............................(iii)
Substituting equation(iii) into equation (i)
z = ax + by + f (a, b) ................(iv)
z = ax + by + log(a) + log(b)......................(ii)
23
To obtain the singular integral we dierentiate (ii) partially with respect to
a and b to obtain (iii) and (iv),
1 1
0=x+ =⇒ a = − ..................(iii)
a x
1 1
0 = y + =⇒ b = − ..................(iv)
b y
We then substitute equations (iii) and (iv) in (ii) to obtain the singular
integral.
1 1 1 1
z= − x+ − y + log − + log −
x y x y
1 1
z = −1 − 1 + log − + log −
x y
1
z = −2 + log
xy
√
(ii) z = px + qy − 2 pq..................(i)
Solution
Comparing this with the given equation of Chaurat's form, the complete
integral is given by
√
z = ax + by − 2 ab..................(ii)
24
r r
b a
(x − z) (y − z) = × =1
a b
(x − z) (y − z) = 1
25
2x2 y δz 2xy 2 δz
1 δz 1 δz
z= · · + +
2x δx 2y δy 4xy δx 4xy δy
1 δz 1 δz 2 1 δz 2 1 δz
z= · · +x · +y ·
2x δx 2y δy 2x δx 2y δy
Let 2xδx = δX and 2yδy = δY
let x2 = X and y 2 = Y
δz δz δz δz
z= +X +Y
δX δY δX δY
δz δz
p= δX
,q = δY
Z = pq + pX + qY
But x2 = X and y 2 = Y
Replacing yields;
z = ax2 + by 2 + ab.................(iv)
For the singular integral,we dierentiate (iv) partially with respect to a and
b to obtain equations (v) and (vi)
0 = x2 + b =⇒ b = −x2 ........................(v)
0 = y 2 + a =⇒ a = −y 2 ...........................(vi)
z = −y 2 x2 − x2 y 2 + x2 y 2 .................(vii)
z = −y 2 x2 ........(viii)
1
(v) Show that the complete integral of the equation z = px+qy+(p2 + q 2 + 1) 2 ..........(i)represents
all planes at a unit distance from the origin.
Solution
This equation is of Chaurit's form, thus
26
12
z = ax + by + a2 + b2 + 1 ............................(ii)
x y z
a + b
+ 1 = 1............(v)
1 1 1
−(a2 +b2 +1) 2 −(a2 +b2 +1) 2 (a2 +b2 +1) 2
p = qa or q = pa.................................(ii)
dz = pdx + qdy
Letu = x + ay
27
dz = pdu
dz
p= ............................(iv)
du
Thus by equation (ii),
dz
q=a ............................(v)
du
Replacing back in equation
(i)
dz dz
f , a , 0 = 0......................(vi)
du du
Equation (vi) is a p.d.e of rst order. Solving (vi), we get z as a function of
U. The Complete Integral is then obtained by replacing u with x+ay.
Example:
a) Find the Complete Integral and the Singular Integral of q 2 = z 2 p2 (1 − p2 ) .....................(i)
Solution
Step 1:
We take p = du
dz
,q = a du
dz
and u = x + ay and substitute in
! (i)
2 2 2
dz dz
2 dz
a =z 1−
du du du
" 2 #
dz
a2 = z 2 1 −
du
2
2 2 2 dz
a =z −z
du
2
z 2 − a2
dz
=
du z2
r √
dz z 2 − a2 z 2 − a2
= =
du z2 z
ˆ ˆ
z
du = √ dz
z − a2
2
28
Squaring both sides,
(u + b)2 = z 2 − a2
Replacing u = x + ay
(x + ay + b)2 = z 2 − a2 ..............................(ii)
29
5 SURFACES AND CURVES IN THREE DI-
MENSIONS:
In rectangular co-ordinates, the equation of a surface is expressed in the form
f (x, y, z) = 0..........................(i). For example, the equation x2 + y 2 + z 2 =
a2 is the equation of a sphere centered at the origin (0, 0, 0)radius r. From
equation (i) the total dierential of F is δf δx
dx + δf
δy
dy + δf
δz
dz = 0showing that
the triad δx , δy , δz are the directional rations normal to the surface at a
δf δf δf
given point while (dx, dy, dz)are the direction ratios tangent to the surface
of the same point. Consider a curve which is the intersection of the surfaces
f (x, y, z)and g (x, y, z) = 0. With the two equations of surfaces satisfying
their total dierentials i.e. df = 0and dg = 0or
δf δf δf
dx + dy + dz = 0
δx δy δz
δg δg δg
dx + dy + dz = 0
δx δy δz
which on solving for dx, we have;
δf δf δf
dx + dy = − dz
δx δy δz
δg δg δg
dx + dy = − dz
δx δy δz
δf δf δf δg δf δg
f, g
δx δy
4 = δg δg = − =J
δx δy δx δy δy δx x, y
− δf dz δf δf δf
δz δy
δz δy δz
δg δg δg δg
− δz dz δy δy δz
dx = =
δf δf f,g
δx δy J x,y
δg δg
δx δy
30
or
f,g
δzJ dx
x,y dz
dx = or = .....................(A)
f,g
J x,y J f,g
y,z
f,g
J x,y
From A and B,
dx dy dz
= =
f,g f,g f,g
J y,z
J z,x
J x,y
31
δf δf δf δf δf δf
δy δz δz δx δx δy
p= or q = or R =
δg δg δg δg δg δg
δy δz δz δx δx δy
The reverse problem is that of nding the curves of intersection for these
surface i.e. We show that the surface f (x, y, z) = c1 will involve only one
arbitrary constant c1 and f (x, y, z) = c2 will involve only one arbitrary con-
stant c2 . Therefore, the curve of intersection will involve two parameters and
we call them a two parameters and we call them a two parameter family of
curves.
32
6 METHODS OF SOLVING dx
p = dy
q = dz
R =λ
(i) The Spot Method
Suppose that U (x, y, z)is a solution of dxp = dyq = dz R
..........................(i)
then δu
δx
dx + δu
δy
dy + δu
δz
dz = 0hence λ δu
δx
dx + δu
δy
dy + δu
δz
dz = 0where λ 6=
0hence p du
dx
+ q du
dy
+ R du
dz
= 0. We then try to spot functions p0 , q 0 and R0 such
that pp0 + qq 0 + RR0 = 0implying that p0 dx + q 0 dy + R0 dz = 0will be an exact
dierential equation with the solution U (x, y, z) = c1 . Repeats the same
process with the function p00 , q 00 and R00 to get V (x, y, z) = c2 which is the
required two parameter family of curves c1 and c2 .
Example:
1. Find the integral curve of the equations y(x+y)+az
dx dy
= x(x+y)−az dz
= z(x+y)
Solution
p = y (x + y) + az , q = x (x + y) − az and R = z (x + y)
By intution take p0 = x, q 0 = y thus pp0 + qq 0 + RR0 will become exact;
xy (x + y) + xaz + yx (x + y) − yaz + R0 z (x + y) = 0
az (x + y) + R0 z (x + y) = 0 =⇒ az (x + y) = −R0 z (x + y) = 0
Thus
R0 = −a
Then pp0 +qq 0 +RR0 = 0is now exact implying that p0 dx+q 0 dy +R0 dz = 0will
be an exact dierential equation with the solution U (x, y, z) = c1
xdx − ydy − adz = 0
On integration,
x2 y 2
− − az = c or x2 − y 2 − 2az = c1
2 2
where c1 = 2c
Again by intuition, take p00 = 1, q 00 = i thus pp00 + qq 00 + RR00 = 0
33
y (x + y) + az + x (x + y) − az + R00 z (x + y) = 0
(x + y)2 + R00 z (x + y) = 0
− (x + y)
R00 =
z
Thus p00 dx + q 00 dy + R00 dz = 0will be an exact dierential equation with the
solution V (x, y, z) = c2
(x + y)
dx + dy − dz = 0
z
dx + dy 1
− dz = 0
x+y z
d (x + y) 1
− dz = 0
x+y z
ln |x + y| − ln |z| = lnc2
x+y
= c2
z
Thus the integral curves are given by U (x2 − y 2 − 2az) = c1 and V x+y
z
= c2
Note: This method can be simplied and shortened by the following prop-
erties:
so that if pp0 + qq 0 + RR0 = 0 then
0 dx+q 0 dy+R0 dz
dx
p
= dyq
= dzR
= p pp 0 +qq 0 +RR0 = dx−dy−2dz
p−q−2R
p0 dx + q 0 dy + R0 dz = 0is exact.
(ii) Method 2
Given that dx
p
= dy
q
= dz
R
, It should be noted that dx
p
= dy
q
= dz
R
= p0 dx+q 0 dy
pp0 +qq 0
=
dx+3dy−dz
pp0 +3qq 0 −RR0
e.t.c
Example:
Find the integral curves of y(x+y)+az
dx dy
= x(x+y)−az dz
= z(x+y)
Solution
Adding the rst two and equating to the third equation,
dx + dy dz
2 =
(x + y) z (x + y)
dx + dy dz
=
x+y z
34
ln |x + y| = ln |z| + lnc1
x+y
= c1
z
Also,
xdx ydy dz
= =
xy (x + y) + azx xy (x + y) − azy z (x + y)
xdx − ydy xdx − ydy
=
x2 y + xy 2 2 2
+ azx − xy − x y + azy az (x + y)
Therefore,
xdx − ydy dz
=
az (x + y) z (x + y)
xdx − ydy = adz
x2 y 2
− = az + c or x2 − y 2 − 2az = c2
2 2
where c2 = 2c
(iii) Method 3:
It is possible to use the rst solution to obtain the second solution by elimi-
nating one of the variables in one of the fractions:
Example: Obtain integral curve of y(x+y)+az
dx dy
= x(x+y)−az dz
= z(x+y) ...............(i)
Solution
Adding the rst two and equating to the third equation,
dx + dy dz
2 =
(x + y) z (x + y)
dx + dy dz
=
x+y z
ln |x + y| = ln |z| + lnc1
x+y
= c1
z
Eliminating z from the rst and second fraction using the rst solution, we
have;
x+y
z=
c1
35
Replacing in equation (i),
dx dy dz
= = x+y
y (x + y) + a x+y
x (x + y) − a x+y
c1
(x + y)
c1 c1
Solution
x is absent from the equation dyy = z+y
dz
2 and this equation is linear with z and
´ 1
the dependent variable with I.F = e − y dy = e−lny = y1 . Therefore
dz z + y2 dz z
= or − =y
dy y dy y
1 dz z
− 2 =1
y dy y
z
d =1
y
z
= y + c1 or z = c1 y + y 2
y
36
Comparing the rst and the second fractions,
dx dy
2
=
x + (c1 y + y ) y
This is a linear equation with x as the dependent variable
dx x + c1 y + y 2
=
dy y
dx x
= + c1 + y
dy y
dx x
− = c1 + y
dy y
´
− y1 dy
I.F = e = e−lny = y1 . Therefore
x 1
d = (y + c1 )
y y
x c1
d =1+
y y
x
= y + c1 lny + c2
y
Replacing c1 = z−y 2
y
x z − y2
=y+ lny + c2
y y
Example:
Find the integral curves of the equation
a) xy
dx
= dy
y2
dz
= zxy−2x 2
37
dy dz
2
=
y zxy − 2x2
But x = c1 y
dy dz
=
y 2 zc1 y − 2c21 y 2
2
dz
dy =
zc1 − 2c21
1
y = ln zc1 − 2c21 + c2
c1
y xz 2x2
y = ln − 2 + c2
x y y
b) dx
xz(z 2 +xy)
= dy
−yz(z 2 +xy)
= dz
z4
dx −dy
=
x y
lnx = −lny + lnc =⇒ xy = c1
Also,
dx dz
2
= 4
xz (z + xy) z
dx z (z 2 + xy)
= dz
x z4
dx 1 c1
= + dz
x z z3
c1
lnx = lnz − 2 + c2
2z
xy
lnx − lnz + 2 = c2
2z
x xy
+ 2 = c2
z 2z
c) x(y−z)
dx dy
= y(z−x) dz
= z(x−y)
Solution
Adding the rst two fractions and equating with the third,
38
dx + dy dz
=
−z (x + y) z (x − y)
dx + dy + dz = 0
Integrating,
x + y + z = c1
Also,
ydx + xdy dz
=
yx (y − z) + xy (z − x) z (x − y)
ydx + xdy dz
=
yx (x − y) z (x − y)
dx dy dz
+ + =0
x y z
lnx + lny + lnz = lnc2
xyz = c2
d) x2 (ydx
3 −z 3 ) =
dy
y 2 (z 3 −x3 )
= dz
z 2 (x3 −y 3 )
Solution
xdx + ydy zdz
= 3 3
x3 (y 3 3 3 3 3
− z ) + y (z − x ) z (x − y 3 )
Considering the denominator,
x3 y 3 − x3 z 3 + y 3 z 3 − x3 y 3 = z 3 x3 − z 3 y 3
z 3 y 3 − x3 = −z 3 x3 − y 3
Therefore,
xdx + ydy zdz
3 3 3
= 3 3
−z (x − y ) z (x − y 3 )
xdx + ydy + zdz = 0
x2 y 2 z 2
+ + = c or x2 + y 2 + z 2 = c1
2 2 2
dx dy dz
x2
+ y2 z2
=
y3 − x3 + z3 − x3 x3 − y 3
39
dx dy dz
x2
+ y2 z2
=
y3 + z3 x3 − y 3
dx dy dz
+ 2 + 2 =0
x2 y z
1 1 1
− − − =c
x y z
1 1 1
+ + = c2
x y z
40
7 ORTHOGONAL TRAJECTORIES OF CURVES:
δf δf δf δf δf δf
p = δy
δg
δz
δg
,q = δg
δz δx
δg
,R = δx
δg
δy
δg
δy δz
δz δx
δx δy
In three dimensions, given a surface f (x, y, z) = 0..............................(i)and
a system of curves, we have to nd a system of curves each of which lies on
the surface (i) and cuts every curve of the given system at right angles. The
new system of curve is called the system of orthogonal trajectories on the
surface of the given system of curves. The original system of curves may be
thought of as the intersection of the surface (i) with the parameter family of
curves g (x, y, z) = c1 ..............................(ii)
41
In general, the tangential direction (dx, dy, dz)to the given curve through the
points x,y,z on the surface (i) satises the equations δfδx
dx + δf
δy
dy + δf δz
dz = 0
and δx
δg δg
dx+ δy dy+ δg
δz
dz = 0. Which on solving gives dxp
= dy
q
= dzR
.............(iii)
where
δf δf f f δf δf f f δf δf f f
,q = δg δg = ,R = δx
y z δz δx z x x y
p = δy δz
= δy
=
δg
δy δg gy gz δz δx gz gx δg δg
δx δy gx gy
δz
The curve through point (x, y, z)of the orthogonal system has tangential
direction (dx0 , dy 0 , dz 0 )which lies on the surface (i) such that δfδx
dx0 + δf
δy
dy 0 +
δf
δz
dz 0 = 0...................(iv)
42
This implies that dxdx0 + dydy 0 + dzdz 0 = 0
Hence from (iii), we have pdx0 + qdy 0 + Rdz 0 = 0...........................(v)
Solving (iv) and (v) we have,
.........................(vi)where p0 = Rfy − qfz , q 0 = pfz − Rfx and
0
dx0 0
p0
= dy
q0
= dz
R0
R0 = qfx − pfy .........................(vii)
Note that:
δf
δx
+ δf
dx0 δy
dy 0 = − δf
δz
dz 0
pdx + qdy 0 = −Rdz 0
0
f δz 0 f f f
z y 0 z y
δz
0
−Rδz q −R q
dx0 = =
fx fy fx fy
p q p q
dx0 dz 0
=
f f f f
x y x y
p q p q
43
f 0
x f z δz f
0 x
f
z
δz
p −Rδz 0
p −R
dy 0 = =
fx fy
f f
x y
p q p q
dy 0 dz 0
=
f f f f
x y x y
p q p q
Example:
Find the orthogonal trajectories of the sphere x2 + y 2 + z 2 = a2 at its inter-
section with a parabola xy = cz where c is a parameter.
Solution
Let f (x, y, z)be identiacal to x2 + y 2 + z 2 = a2
Let g (x, y, z) = xyz = c
δf δf f f
Then (dx, dy, dz)satises dx
p
= dy
q
= dz
R
where p = δy
δg
δz
δg
y z
= =
δy δz
gy gz
2y 2z
2xy 2 −2xy 2 −2xz 2
= − − 2x =
z 2 z2
z − xy
x
z2
δf δf f f 2z 2x
2y (z 2 +x2 )
,q = δg z x 2 2 2y
δz δx
= = xy y = 2yz + 2xz2 y = 2yz z+2x =
δg z 2 z2
δz δx gz gx − z 2 z
and
δf δf f f 2x 2y
x y 2 2 2(x2 −y 2 )
R = δx δy
= = y x = 2x −2y =
δg δg
δx δy gx gy z z z
z
dx dy dz
2 2 = 2y(z2 +x2 ) = 2(x2 −y2 )
−2x y z+z 2 z2 z
dx dy dz
2 2
= 2 2
=
−x (y + z ) y (z + x ) z (x − y 2 )
2
44
−2yz (y 2 + z 2 )
f f 2z 2x
z x
q0 = = = −2zx (y 2 + z 2 )−2xz (x2 − y 2 ) =
R p z (x2 − y 2 ) −x (y 2 + z 2 )
−2xz(x2 + z2 )
f f 2x 2y
0 x y
R = = = 2xy (x2 + z 2 )+2xy (y 2 + z 2 ) =
2 2 2
p q −x (y + z ) y (x + z ) 2
2xy (x2 + 2z 2 + y 2 )
Hence,
dx0 dy 0 dz 0
= =
−2yz (y 2 + z 2 ) −2xz (x2 + z 2 ) 2xy (x2 + 2z 2 + y 2 )
dx0 dy 0 dz 0
= =
−yz (y 2 + z 2 ) −xz (x2 + z 2 ) xy (x2 + 2z 2 + y 2 )
For the rst two fractions,
dx0 dy 0
=
−yz (y 2 + z 2 ) −xz (x2 + z 2 )
dx0 dy 0
=
y (y 2 + z 2 ) x (x2 + z 2 )
xdx0 ydy 0
= 2
(y 2 + z 2 ) (x + z 2 )
From;x2 + y 2 + z 2 = a2 ,x2 + z 2 = a2 − y 2 and y 2 + z 2 = a2 − x2
xdx0 ydy 0
=
(a2 − x2 ) (a2 − y 2 )
1 1
− ln a2 − x2 = − ln a2 − y 2 + lnk
2 2
a2 − y 2
= k −2 = c
a2 − x 2
Thus aa2 −x 2 = cand the given surface equation is x + y + z = a2 are the
2
−y 2
2 2 2
45
The triad (dx, dy, dz)satises the equation dx
= dy = dz
p q R
δf δf f fz 2y −2ztan α
2
y
p = δy δz
= = = 2y
δg δg gy
δy δz
gz 0 1
δf δf f fx −2ztan2 α 2x
,q = δg
δz δx
δg =
z
= = −2x and
δz δx
gz gx 1 0
δf δf f fy 2x 2y
x
R = δx δy
= = =0
δg δg
gx
δx δy
gy 0 0
dx dy dz dx dy dz
= = =⇒ = =
2y −2x 0 y −x 0
Then the triad
(dx0 , dy 0 , dz 0 )satises
dx0
p0
= dy 0
q0
= dz 0
R0
where
f
y fz 2y −2ztan2 α
p0 = = = −4xztan2 α
q R −2x 0
f
z fx −2ztan2 α 2x
q0 = = = −4yztan2 α
R p 0 2y
f fy 2x 2y
0 x
R = = = −4 (x2 + y 2 )
p q 2y −2x
dx0 dy 0 dz 0
= =
−4xztan2 α −4yztan2 α −4 (x2 + y 2 )
dx0 dy 0 dz 0
= =
xztan2 α yztan2 α (x2 + y 2 )
Comparing the rst two fractions,
dx0 dy 0
=
x y
lnx = lny + lnc
x
=c
y
Therefore the required system of orthogonal trajectories are given by x2 +y 2 =
z 2 tan2 αand xy = c.
c) Find the orthogonal trajectories on the conicoid (x + y) z = 1of the conics
in which its cut by the system of planes x − y + z = k where k is a parameter.
46
Solution
The triad (dx, dy, dz)satises the equation dx
p
= dy
q
= dz
R
f = xz + yz − 1 = 0
g = x − y + z − k
δf δf f f z x + y
y z
p = δy δz
= = =z+x+y
δg δg
δy δz
gy gz −1 1
δf δf f f x + y z
,q = δg δg =
δz δx z x
= = x + y − z and
δz δx gz gx 1 1
δf δf f f z z
δx δy x y
R = δg δg = = = −z − z = −2z
δx δy gx gy 1 −1
dx dy dz dx dy dz
= = =⇒ = =
x+y+z x+y−z −2z y −x 0
Then the triad
(dx0 , dy 0 , dz 0 )satises
dx0
p0
= dy 0
q0
= dz 0
R0
where
f f z y + x
0 y z
p = = = −2z 2 − (x + y) (x + y − z)
q R x + y − z −2z
f f x + y z
z x
q0 = = = (x + y) (x + y + z) + 2z 2
R p −2z x + y + z
f f z z
x y
R0 = = = z (x + y − z) − z (x + y + z) =
p q x + y + z x + y − z
−2z 2
dx0 dy 0 dz 0
= =
−2z 2 − (x + y) (x + y − z) (x + y) (x + y + z) + 2z 2 −2z 2
Adding the rst two fractions,
dx0 + dy 0 dz 0
=
−2z 2 − (x + y) (x + y − z) + (x + y) (x + y + z) + 2z 2 −2z 2
dx0 + dy 0 dz 0
=
2z (x + y) −2z 2
dx0 + dy 0 dz 0
=
(x + y) −z
ln |x + y| = −ln |z| + lnc
47
(x + y) z = c
If c = 1,
(x + y) z = 1
48
8 METHOD OF SEPARATION OF VARIABLES:
Example 1:
Displacement of vibrating spring or wave equation.
The partial dierential equation governing the motion of the vibrating string
is given as δδt2u = c2 δδxu2 where c2 = TM
2 2 ension
ass
Therefore our task is to determine u (x, t)under the boundary conditions
u (0, t) = 0,u (1, t) = 0 and the initial conditions u (x, 0) = φ(x)and ut (x, 0) =
ϕ(x). Here u (x, t)represents displacement of the string at any position and
at any time t. u(t)represents the velocity of the wave produced by the
string at time t. By the method of separation of variables, the depen-
dent variable u (x, t)may be expressed as u (x, t) = X(x)T (t)where X(x)is
a function of x alone and T (t)is a function of t alone. Then it follows that
uxx = X 00 (x)T (t)andutt = X 00 (x)T (t)
Given,
δ2u 2
2δ u
= c ....................(i)
δt2 δx2
49
Let U = XT........................(ii)
δu δX δT δX
= T +X = T = X 0 T.........................(iii)
δx δx δx δx
δ2u δ 0 δX 0 δT δ2X
2
= (X T ) = T + X = 2
T = X 00 T.............................(iv)
δx δx δx δx δx
Similarly,
δ2u
= T 00 X.........................(v)
δt2
Substituting (v) and (iv) into (i)
XT 00 = c2 X 00 T
T 00 c2 X 00
= =k
T X
These two functions could only be equal to a constant, say k as shown above,
T 00
= k =⇒ T 00 − kT = 0
T
and
c2 X 00 k
= k =⇒ X 00 − 2 X = 0
X c
Assuming that we need to consider only the physical solution, there are three
cases to be be investigated.
Case 1: k=0
If k=0, then T 00 = 0, T 0 = A,T = At + B and X 00 = 0, X 0 = C ,X = Cx + D
Therefore, U (x, t) = X(x)T (t) =⇒ U (x, t) = (Cx + D) (At + B)
These solutions cannot describe the vibtrating of a system because it is not
periodic.
Applying U (0, t) = 0,U (1, t) = 0
D (At + B) = 0
Dt = 0, t 6= 0
hence D = 0
U (x, t) = Cx (At + B)
0 = Ct
50
hence C = 0
This shows that U (x, t) = 0has trivial solutions.
Case 2:k = λ2
Then the two dierential equations are
T 00
= λ2 =⇒ T 00 − λ2 T = 0
T
and
c2 X 00 λ2
= λ2 =⇒ X 00 − 2 X = 0
X c
T = Aeλt + Be−λt
λ λ
X = Ce c x + De− c x
Thus
λ
− λc x
x
Aeλt + Be−λt
U (x, t) = Ce c + De
51
This solution is periodic with a period of 2π
λ
. However, we need to determine
λ, A, B, C, andD from the rst boundary condition.
U (0, t) = 0 =⇒ 0 = C (Acosλt + Bsinλt) =⇒ 0 = CT, T 6= 0 =⇒ C = 0
Hence U (x, t) = Dsin λc x [Acosλt + Bsinλt]
This series can be organized as the half range since the expansion of a peri-
odic function φ(x)is denes in the range (0, l). Now An can be obtained by
multiplying the series by sin nπl xand integrating with respect to x from 0 to
´
l which gives An = π2 0l φ(x)sin nπl xdx
∞
a0 X
F (x) = + (an cos(nx) + bn sin(nx))
2 n=1
52
∞
X nπc nπc
ψ(x, t) = Bn sin x
n=1
l l
Then multiplying both sides of this equation by sin nπl xand integrating with
respect to x, we get ˆ l
2 nπ
Bn = ψ(x)sin xdx
nπc 0 l
Therefore the equation :
nπ nπc nπc
Un (x, t) = sin
x An cos t + Bn sin t
l l l
is the required solution of the wave equation with An and Bn as shown above.
Example:
a) A string is stretched and fastened to two points. The motion is started
by displacing the string in the form y (x, t) = asin πxl from which it is
53
δ2u 2
2δ u
= c ....................(i)
δt2 δx2
Boundary conditions are y (0, t) = 0and y (l, t) = 0........................(ii)
The string is in the shape of y (x, t) = asin πxl .........................(iii)at rest
Therefore a1 = 0since the solution must be true for all values of t, hence
equation (iv) reduces to
Y (x, t) = a2 sinkx (a3 coskct + a4 sinkct)
Let a2 a3 = b1 and a2 a4 = b2
Y (x, t) = sinkx (b1 coskct + b2 sinkct) .....................(v)
sinkl = 0
nπ
k= .........................(vi)
l
Using (vi) in (v),
nπ nπ nπ
Y (x, t) = sin x b1 cos ct + b2 sin ct .....................(vii)
l l l
Using = 0 at t= 0
dy
dt
nπ nπc nπ nπc nπ
0 = sin x −b1 sin ct + b2 cos ct .....................(viii)
l l l l l
54
nπ nπc
0 = sinx 0 + b2 ...................(ix)
l l
Since the answer must be true for all values of x, b2 = 0 and equation (vii)
becomes;
nπ nπ
Y (x, t) = sin x b1 cos ct
l l
nπ nπ
Y (x, t) = b1 sin xcos ct....................(x)
l l
Applying the initial conditiont = 0 in (x),
nπ nπ
Y (x, 0) = b1 sin xcos(0) = b1 sin x.......................(xi)
l l
Comparing equation (xi) with te givenequation a = b1 , n = 1
πx πct
Y (x, t) = asin cos ....................(xii)
l l
b) A tightly stretched spring with xed ends at x = 0and x = l is initially
in equilibrium position. It is vibrating by giving each of the poins an initial
velocity, Y (x, 0) = 0, ∂y 3 πx
. Find the dispacement Y (x, t).
∂t t=0
= bsin l
Solution
55
Equation of a vibrating string is
δ2u 2
2δ u
= c ....................(i)
δt2 δx2
We solve equation (i) subject to the bondary conditionsy (0, t) = 0................(ii)and
y (l, t) = 0........................(iii)
Initial conditions are y (x, 0) = 0................(iv)and ∂y πx
....................(v).
∂t t=0
= bsin3 l
The solution to equation (i) is of the form
Y (x, t) = (a1 coskx + a2 sinkx) (a3 coskct + a4 sinkct) .........................(vi)
Therefore a1 = 0since the solution must be true for all values of t, hence
equation (vii) reduces to
Y (x, t) = a2 sinkx (a3 coskct + a4 sinkct)
Let a2 a3 = b1 and a2 a4 = b2
Y (x, t) = sinkx (b1 coskct + b2 sinkct) .....................(viii)
sinkl = 0
nπ
k= .........................(ix)
l
Using (ix) in (viii),
nπ nπ nπ
Y (x, t) = sin x b1 cos ct + b2 sin ct .....................(x)
l l l
Using initial condition, y (x, 0) = 0
nπ
0 = sin x (b1 cos0 + b2 sin0) .....................(xi)
l
56
nπ
x (b1 ) ...................(xii)
0 = sin
l
Since the answer must be true for all values of x, b1 = 0 and equation (x)
becomes;
nπ nπ
Y (x, t) = sin x b2 sin ct
l l
nπ nπ
Y (x, t) = b2 sin xsin ct....................(xiii)
l l
where n=1,2,3,...
Therefore the sum of this solution is also a solution by superposition principle
of linear solutions ∞ X nπ nπ
Y (x, t) = An sin xsin ct....................(xiv)
n=0
l l
57
Equating coecient of sin πx
l
3 πc 3bl
b = A1 =⇒ A1 =
4 l 4πc
Equating coecient of sin 2πl x
2πc
0 = A2 =⇒ A2 = 0
l
Equating coecient of sin 3πx
l
1 3πc −bl
− b = A3 =⇒ A3 =
4 l 12πc
Equating coecient of sin 4πl x
4πc
0 = A4 =⇒ A4 = 0
l
Substitututing A1 , A2 , A3 , A4 into equation (xvii)
3bl πx πct bl 3πx 3πct
Y (x, t) = sin sin − sin sin
4πc l l 12πc l l
Equation (vii) represents displacement of the wave.
58
9 DERIVATION OF WAVE EQUATIONS
The wave equation is derived with respect to travelling transverse waves on
a taut string. The following assumptions are made:
(i) The amplitude of the wave is always small.
(ii) the string has a constant linear mass density
(iii) The string is under a constant tension of magnitude T.
(iv) The string is not made of atoms or molecules. It is compresses of con-
tinuous non-quantized matter.
(v) There is no damping(no loss of energy due to friction)
(vi) Gravity is ignored in the treatment.
INTRODUCTION
Through out the derivation, we will focus on a tiny bity of strings, so small
that we can pretend that it is almost straight line. The bit of a string which
will have a tension pulling at each end(with magnitude T ). Length of the
bit of string will be 4s.
Assuming that the amptitude of the wave is extremely small, then the ap-
59
proximation of the string length is 4x. Dening a linear mass density as the
mass of the string in kilograms per meter of string to be µ, then the mass of
bit of string is given by: m = µ4x.........................................(i).
Since we are looking at transverse travelling waves, we worry about the net
force in the vertical direction. X
Fnet = Fy
Note that: Deing vertical coordinates above, onm the string to be positive.
The net force on the bit of the string is:
Fnet = T sin (θ2 ) − T sin (θ1 )
60
Since we are dealing only with small amplitudes, we use this approximation
4x = 4s. Therefore, sinθ = tanθ
Then the new net force is expressed as
Fnet = T tan (θ2 ) − T tan (θ1 )
61
derivative of vertical displacements
δ2y
Fnet = m ............................(iv)
δt2
If 4xgets smaller
and smaller,
equation
(vi) becomes:
δ2y δ2y
T =µ ............................(vii)
δx2 δt2
10 HEAT EQUATION
Consider the heat conduction equation in a thin metal bar of legth l with
insulated sides. Let us assume that the ends are x = 0and x = l and
62
are hels at temperature 00 cfor all timet > 0. In addition, let us assume
that the temperature distribution at t=0 is U (x, 0) = f (x). Determine the
temperature distribution in a bar of some subsequent time t > 0.
Solution
The heat equation is given by Ut = αUxx i.e.
δu δ2u
=α 2
δt δx
where U (x, t)is the heat distribution subject to boundary conditions U (0, t) =
0and U (1, t) = 0and the initial conditions U (x, 0) = f (x).
By method of separation of variables,
Let U = XT........................(ii)
δu δX δT δX
= T +X = T = X 0 T.........................(iii)
δx δx δx δx
δ2u δ 0 δX 0 δT δ2X
= (X T ) = T + X = T = X 00 T.............................(iv)
δx2 δx δx δx δx2
Similarly,
δu
= T 0 X.........................(v)
δt
Substituting (v) and (iv) into (i)
T 0 X = αX 00 T
T0 X 00
= =k
αT X
These two functions could only be equal to a constant, say k as shown above,
T0
= k =⇒ T 0 − αkT = 0
αT
and
X 00
= k =⇒ X 00 − kX = 0
X
Assuming that we need to consider only the physical solution, there are three
cases to be be investigated.
Case 1: k=0
If k=0, then T 0 = 0, T = Aand X 00 = 0, X 0 = C ,X = Cx + D
63
Therefore, U (x, t) = X(x)T (t) =⇒ U (x, t) = (Cx + D) (A) = CAx + DA =
Ex + F
These solutions cannot describe the vibtrating of a system because it is not
periodic.
Applying U (0, t) = 0,U (1, t) = 0
F = 0 T 6= 0
hence F = 0
U (x, t) = Ex X 6= 0
hence E = 0
This shows that U (x, t) = 0has trivial solutions.
Case 2:k = λ2
Then the two dierential equations are
T 0 X = αX 00 T
T0 X 00
= = λ2
αT X
These two functions could only be equal to a constant, say k as shown above,
T0
= λ2 =⇒ T 0 − αλ2 T = 0
αT
δT δT 2
= αλ2 T =⇒ = αλ2 δt =⇒ lnT = αλ2 t + A ⇒ T = Aeαλ t
δt T
and
X 00
= k =⇒ X 00 − λ2 X = 0
X
2t
T = Aeαλ
X = Ceλx + De−λx
Thus
2
U (x, t) = Ceλx + De−λx Aeαλ t
2
U (x, t) = Beλx + F e−λx eαλ t
but when t increases the temperature in this case increases too much instead
64
of decreasing i.e. physical condition does not allow. Therefore we reject the
solution;
Case 3
k = −λ2
Then the two dierential equations are
T 0 X = αX 00 T
T0 X 00
= = −λ2
αT X
These two functions could only be equal to a constant, say k as shown above,
T0
= λ2 =⇒ T 0 + αλ2 T = 0
αT
δT δT 2
= −αλ2 T =⇒ = −αλ2 δt =⇒ lnT = −αλ2 t + A ⇒ T = Ae−αλ t
δt T
and
X 00
= λ2 =⇒ X 00 + λ2 X = 0
X
2t
T = Ae−αλ
X = Ccosλx + Dsinλx
Thus
2t
U (x, t) = (Ccosλx + Dsinλx) Ae−αλ
2t
U (x, t) = (Bcosλx + F sinλx) e−αλ
B=0
2t
U (x, t) = (F sinλx) e−αλ
b) U (l, t) = 0
65
2t
0 = (F sinλl) e−αλ
Here, F 6= 0, e−αλ t =
6 0. Thus,
2
sinλl = 0
nπ
λl = nπ =⇒ λ =
l
Then nπ −α( nπl )2 t
U (x, t) = F sin x e
l
λn are dened to be the eigen values and the function Un (x, t)corresponds to
any λn is called the eigen function. Thus the general solution becomes
nπ nπ 2
Un (x, t) = Fn sin x e−α( l ) t
l
From the superposition principle all the linear solutions, we get the general
solutions as ∞ ∞
X X nπ −α( nπl )2 t
U (x, t) = Un (x, t) = Fn sin x e
n=1 n=1
l
This can be recognized as the fourier half range sine series and Fn is the
fourier coecient which can be determined from the initial condition as fol-
lows: ∞ X nπ
U (x, t) = Fn sin x
n=1
l
Assignment
1. Use the method of separation of variables to solve the heat equation.
66
2. By stating the necessary assumption show that the one dimensional wave
equation is given by δδt2y = c2 δx 2 where c = m = .
2 δ2 y µ
2 T ension
M ass
3. Derive the heat equation.
Example 2:
A rod whose surface is measured has a length of 3 units. The end of the rod
is kept at 00 cand its initial temperature at any point x, 0 < x < 3is given
by U (x, 0) = 5sin4πx − 3sin8πx + 2sin10πx. Find the temperature at any
given time t.
Solution
Let U (x, t)be temperature at any given time t.
The boundary conditions are x = 3and x = 0
Initial conditions are U (x, 0) = 5sin4πx − 3sin8πx + 2sin10πx.
The heat equation is one dimension
δu δ2u
=α 2
δt δx
Using case 3 wherek = −λ2 since case 1 and case 2 give trivial solutions;
Then the two dierential equations are
T 0 X = αX 00 T
T0 X 00
= = −λ2
αT X
These two functions could only be equal to a constant, say k as shown above,
T0
= λ2 =⇒ T 0 + αλ2 T = 0
αT
δT δT 2
= −αλ2 T =⇒ = −αλ2 δt =⇒ lnT = −αλ2 t + A ⇒ T = Ae−αλ t
δt T
and
X 00
= λ2 =⇒ X 00 + λ2 X = 0
X
2t
T = Ae−αλ
X = Ccosλx + Dsinλx
67
Thus
2t
U (x, t) = (Ccosλx + Dsinλx) Ae−αλ
2t
U (x, t) = (Bcosλx + F sinλx) e−αλ
B=0
2t
U (x, t) = (F sinλx) e−αλ
b) U (3, t) = 0
2t
0 = (F sinλ3) e−αλ
Here, F 6= 0, e−αλ t =
6 0. Thus,
2
sinλ3 = 0
nπ
λl = nπ =⇒ λ =
3
Then nπ −α( nπ 3 )
2
t
U (x, t) = F sin x e
3
λn are dened to be the eigen values and the function Un (x, t)corresponds to
any λn is called the eigen function. Thus the general solution becomes
nπ nπ 2
Un (x, t) = Fn sin x e−α( 3 ) t
3
From the superposition principle all the linear solutions, we get the general
solutions as ∞ ∞
X X nπ −α( nπ
3 )
2
t
U (x, t) = Un (x, t) = Fn sin x e
n=1 n=1
3
68
∞ ∞
X X nπ π 2π 3π 4π 5
U (x, 0) = Un (x, 0) = Fn sin x = F1 sin x+F2 sin x+F3 sin x+F4 sin x+F5 sin
n=1 n=1
3 3 3 3 3 3
Example:
A rod of length l and with insulated side is initially at a uniform temper-
ature U0 . It ends are suddenly cooled at 00 cand kept at that tempera-
ture. Find the temperature distribution function U (x, t). Ans: U (x, t) =
c 2 n2 π 2
nπ
xe−
t
U0 sin l
l
Example:
A bar whose surface is insulated and has a length of 3 units and diusivity
2 units. If the ends are kept at temperature 0 units and initial tempera-
ture is U (x, 0) = 5sin4πx − 3sin8πx + 2sin10πx. Find the sistribution of
temperature at any point x at any time t.
Solution
69
Let the bar be placed on the x-ais such that one end is at the origin and the
other end is at x = 3. The distribution of heat is generally given by
Let U (x, t)be temperature at any given time t.
The boundary conditions are x = 3and x = 0
Initial conditions are U (x, 0) = 5sin4πx − 3sin8πx + 2sin10πx.
The heat equation is one dimension
δu δ2u
=2 2
δt δx
Using case 3 wherek = −λ2 since case 1 and case 2 give trivial solutions;
Then the two dierential equations are
T 0 X = 2X 00 T
T0 X 00
= = −λ2
2T X
These two functions could only be equal to a constant, say k as shown above,
T0
= λ2 =⇒ T 0 + 2λ2 T = 0
αT
δT δT 2
= −2λ2 T =⇒ = −2λ2 δt =⇒ lnT = −2λ2 t + A ⇒ T = Ae−2λ t
δt T
70
and
X 00
= λ2 =⇒ X 00 + λ2 X = 0
X
2t
T = Ae−2λ
X = Ccosλx + Dsinλx
Thus
2t
U (x, t) = (Ccosλx + Dsinλx) Ae−2λ
2t
U (x, t) = (Bcosλx + F sinλx) e−2λ
B=0
2t
U (x, t) = (F sinλx) e−2λ
b) U (3, t) = 0
2t
0 = (F sinλ3) e−2λ
Here, F 6= 0, e−αλ t =
6 0. Thus,
2
sinλ3 = 0
nπ
λl = nπ =⇒ λ =
3
Then nπ nπ 2
U (x, t) = F sin x e−2( 3 ) t
3
λn are dened to be the eigen values and the function Un (x, t)corresponds to
71
any λn is called the eigen function. Thus the general solution becomes
nπ nπ 2
Un (x, t) = Fn sin x e−2( 3 ) t
3
From the superposition principle all the linear solutions, we get the general
solutions as ∞ ∞
X X nπ −2( nπ
3 )
2
t
U (x, t) = Un (x, t) = Fn sin x e
n=1 n=1
3
∞ ∞
X X nπ π 2π 3π 4π 5
U (x, 0) = Un (x, 0) = Fn sin x = F1 sin x+F2 sin x+F3 sin x+F4 sin x+F5 sin
n=1 n=1
3 3 3 3 3 3
72
Let the tension at P and θbe T1 and T2 respectively.
Horizontal T along PA is −T1 cosθ1 , Vertical T along PA is −T1 sinθ1 .
Horizontal T along θB is −T2 sinθ2 , Vertical T along PA is −T2 cosθ2 .
PQ does not move horizontally. This implies that T1 cosθ1 = T2 cosθ2 =
T.......................(i)
Now, mass of PQ is given by linear denity multiplied by length. 4s · ρ
Let the displacement of PQ be u1 . This implies that velocity = δu
δt
and accel-
eration is δδt2u . Then from Newton's second law of motion F = ma.
2
δ2u
4s · ρ = −T1 sinθ1 + T2 sinθ2
δt2
δ2u
4s · ρ = T2 sinθ2 − T1 sinθ1
δt2
Dividing left hand side and right hand side by T
2
4s · ρ δδt2u T2 sinθ2 T1 sinθ1
= −
T T T
73
2
4s · ρ δδt2u T2 sinθ2 T1 sinθ1
= −
T T2 cosθ2 T1 cosθ1
2
4s · ρ δδt2u
= T anθ2 − T anθ1
T
Now T anθ2 = δu |at Q (slope or gradient)
δx x+4x
Similarly T anθ1 = δu |at P (slope or gradient)
δx x
4s · ρ δ 2 u δu δu
2
= |x+4x − |x
T δt δx δx
δu
δ2u T |
δx x+4x
− δu |
δx x
2
=
δt ρ 4x
δ2u T δ2u
=
δt2 ρ δx2
where 4x = 4s
74
Let U (x, l)be the temperature at any point x at any time t. Let one end
be the origin and the bar be kept on x axis so that the other end is kept at
x = l. At A consider the cross-sectional area a at planepand another plane
p0 at x + 4x
Let pp0 = 4x → 0. Let temperature atpq beu and atp0 q 0 be u + 4u. The
heat ow from pq to p0 q 0 is given by −ka 4u
4x
4t |at x . a is the area while k is
the constant of proportionality(thermal conductivity).
By fouries law, the heat ux is directly proportional to the temperature
gradient.
The heat fow from p0 q 0 to the next string is given by −ka 4u 4x
4t |at x+4x .
Therefore the heat retained is given by:
4u 4u
ka 4t |x+4x −ka 4t |x
4x 4x
The amount of heat retained by the small element Also, heat gained or
lost frompqp0 q 0 is given by M sθwhere M is the mass (volume multiplied by
density), s is the specic heat capacity while θis the temperature change.
= (a4x × ρ) × s4u
75
4u 4u
ρas4u4x = ka 4t |x+4x −ka 4t |x
4x 4x
With the limit 4x → 0, 4u −→ 0 in time 4t
Dividing both sides by 4tthen by 4x
4u ka 4u |
4x x+4x
−ka 4u
4x x
|
ρas =
4t 4x
δu ka δu δu
= |x+4x − |x
δt ρas δx δx
δu k δ2u
=
δt ρs δx2
By denition of a p.d.e k
ρs
= c2 . c2 is a positive quantity and is not equal to
zero.
δu δ2u
= c2 2
δt δx
Y 0 − cY = 0
76
Integrating factor is given by ´
−4cdx
I.F = e = e−4cx
´
−cdy
I.F = e = e−cy
X = Ae4cx
Y = Becy
= Dec(4x+y)
c = −7, D = 5
Thus,
U (x, y) = 5e−7(4x+y)
Exercise:
Solve the following boundary value problems:
(i) δu
δx
= 2 δu
δy
+ u ; U (x, 0) = 6e−3x
(ii) 3 δu
δx
+ 2 δu
δy
= 0 ; U (x, 0) = 4ex
Answers
(i) U (x, t) = 6e−3x−2t
(ii)U (x, t) = 4e−x− 2 y
3
77
ˆ ∞
L [U (x, t)] = U (x, t) dt
0
ˆ ∞ ˆ ∞
δu −st d
L [Ux ] = e dt = e−st U (x, t) dt
0 δx dx 0
ˆ ∞ 2 ˆ ∞ 2
δ u −st d
L [Uxx ] = 2
e dt = 2 e−st U (x, t) dt
0 δx dx 0
ˆ ∞
δu δu −st
L [Ut ] = L = e dt
δt 0 δt
δu
u = e−st =⇒ = −se−st
δt
δv δu
= =⇒ v = u
δt δt
ˆ ∞
δu −st ∞
L [Ut ] = L = U e |0 +s U e−st dt
δt 0
−
= 0 − u(0) + su
−
= u(0) + su
The rst term of the right hand side vanishes at the lower limit through the
sine term. It also vanishes at the upper limit if U (x, t)is such that Ux → ∞as
x → ∞. Thus
78
ˆ ∞
sinγxUxx dx = γ [u(0, t)] − γ 2 Fs {u}
0
´∞
where Fs {u} = 0 U (x, t)sinγxdxand U (0, t)is the value of U at x = 0.
Example
If U (x, t)is the temperature at time t and γ the thermal diusivity of the
metal bar nd U (x, t) from the partial dierential equation. Ut = αUxx ,
x > 0 and t > 0with the boundary conditions U (0, t) = U0 , t > 0 and initial
condition U (x, 0) = 0, x > 0 .
Solution
Here we use the sine transform. Thus, taking the sine transform of Ut = αUxx ,
x > 0 and t > 0, weˆhave: ˆ ∞
∞
Ut sinγxdx = α Uxx sinγxdx
0 0
or
du −
2
= α γU (0, t) − γ U
dt
− −
or ddtu + αγ 2 u = αγU0 using equation (ii) which is an ordinary dierential
− −
equation. The soluion is given by U = Ae−αγ t + γ1 U0 at t = 0, U = 0.
2
− h i
Therefore U = 1 − e−αγ t . The inversion is
U0 2
γ
ˆ ∞
2U0 2t
sinγx
U (x, t) = 1 − e−αγ dx
π 0 γ
79
Solution
This boundary value problem for Laplace's equation is dened as the dirichilet
problem. In obtaining the solution of (i), we make the following transforma-
tion. x = Rcosθand y = Rsinθ. Then the laplacee's equation can be reduced
to:
1 1
URR + UR + 2 Uθθ = 0
R R
where U = U (R, θ).
R2 = x2 + y 2 , θ = tan−1 xy
− y2
dR
dx
= x
R
= cosθ , δθ
δx
= x
y2
= −Rsinθ
R2
= − R1 sinθ
1+
x2
1
dR
dy
= y
R
= sinθ, δy
δθ
= x
2 = x
R2
= Rcosθ
R2
= cosθ
R
1+ y 2
x
We know that δ
δx
= δRδ δR
δx
+ δθδ δx
δθ δ
= cosθ δR − R1 sinθ δθδ and
δ
δy
= δ δR
δR δy
+ δθδ δy
δθ
=
δ 1 δ
sinθ δR + R cosθ δθ
δ2 u δ δu δ δ δ δ
− R1 sinθ δθδ cosθ δR 1
sinθ δθδ
δx2
= δx δx
= δx δx
u = cosθ δR − R
u =
δu 1 δu δu 1 δu
cosθ δR− R sinθ
δθ cosθ δR − R sinθ δθ
δ δu 1 1 δu δ δu 1 δu
= cosθ cosθ − sinθ
− sinθ cosθ − sinθ
δR δR R R δθ δθ δR R δθ
δ2u δ2u δ2u
1 δu 1 1 δu 1 δu 1
= cosθ cosθ 2 + 2 sinθ − sinθ − sinθ −sinθ + cosθ − cosθ − s
δR R δθ R δuδθ R δR δRδθ R δR R
δ2u 1 δu 1 δ2u 1 1 δ2u 1 δu 1
= cos2 θ 2
+ 2
cosθsinθ − sinθcosθ + sin 2
θ− sinθcosθ + 2 sinθcosθ + 2 sin
δR R δθ R δuδθ R R δRδθ R δR R
80
δ2u 2 δu 2 δ2u 1 2 δu 1 2 δ u
2
= cos2 θ + cosθsinθ − sinθcosθ + sin θ + sin θ .........................(a)
δR2 R2 δθ R δRδθ R2 δR R2 δθ2
δ2 u δ δu δ δ δ
+ R1 cosθ δθδ sinθ δR δ
+ R1 cosθ δθδ u =
δy 2
= δy δy
= δy δy
u = sinθ δR
δu 1 δu δu 1 δu
sinθ δR +
R δθcosθ sinθ δR
+ cosθ
R δθ
δ δu 1 δu 1 δ δu 1 δu
= sinθ sinθ + cosθ + cosθ sinθ + cosθ
δR δR R δθ R δθ δR R δθ
δ2u 2 δu 2 δ2u 1 2 δu 1 2
2 δ u
= sin2 θ + sinθcosθ − sinθcosθ + cos θ + cos θ ..........................(b)
δR2 R2 δθ R2 δRδθ R δR R2 δθ2
Adding (a) and (b),
θ4 + n2 θ = 0 ; n = 0, 1, 2, .........................(iv)
Since
1 0 1
R00 θ + R θ + 2 θ00 R = 0
R R
1 0 1
R + R θ + 2 θ00 R = 0
00
R R
1 0 1
00
R + R θ = − 2 θ00 R
R R
81
or
R00 + R1 R0 1 θ00
=− 2 = n2
R R θ
Equation (iii) is a cauchy euler type of equation and its solution can be
obtained by changing the independet variables. R = ez or logR = z
Now
dR dR dz
= ·
dR dz dR
1 dz
= ·
R dR
1
= DR
R
where D = d
dt
d2 R
d 1
2
= DR
dR dR R
1 dR 1 d dR
=− 2 +
R dz x dz dz
1 dR 1 dz d dR
=− 2 + · ·
R dz R dR dz dz
1 dy 1 d2 R
=− + ·
R2 dt R2 dz 2
1 1 2 1
= 2 D2 − D R
=− 2
DR + 2 DR
R R R
1
= D (D − 1) R
R2
82
Thus,
dR dR
R = DR =
dR dz
d2 R d2 R dR
R2 = D (D − 1)R = −
dR2 dz 2 dz
On substituting these values in equation (iii),
d2 R dR dR
− + − n2 R = 0
dz 2 dz dz
The solution is:
83
of variables for the solution of equation (i).
Example 3:
A transmission line of negligible resistance and capacitance has its sending
end at x = 0and receiving end at x = l. A constant voltage E0 is applied to
the sending end while an open circuit is maintained at the receiving end so
that the current there is zero. Assuming that the initial voltage and current
are zero, determine the voltage at any position x at any time t > 0
Solution
The telephone equation can be written as
d2 E d2 E dE
2
= LC 2
+ (RC + LG) + RGE
dx dt dT
d2 E dE
2
= LG ...................................(i)
dx dT
with boundary conditions given as:
x = 0, E (0, t) = E0 and x = L, dX
dE
(L, t) = 0........................................(ii)
The initial conditions are given as
t = 0, U (x, 0) = 0
Taking Laplace transforms of equation (i), we have
d2 E
−
L 2 = LG S E − E (0, 0)
dx
d2 E
− LGSE = 0
dx2
√ √
E = Ae LGSx
+ Be− LGSx
84
x = 0, L E (0, t) = E (0, t) = Es0
x = L, L E (L, t) = E (0) = 0
The application of the above boundary conditions yields
x = 0, A + B = Es0
√ √
x = l, Ae LGSl + Be− LGSl = 0
Solving the two equations, we have;
√
E0 E0 e2L Las
A= √ and B = √
s 1 + e2L Las s 1 + e2L Las
Then substituting A and B we have;
√
E0 √ E0 e2L Las √
E= √ e LGSx
+ √ e− LGSx
s 1 + e2L Las s 1 + e2L Las
√
E0 cosh Las (L − x)
E= √
scosh Las · L
Using the inverse laplace transforms, we have;
" ∞
#
4 X (−1)n (2n−1)2 π2 (2n − 1) (L − x) π
E (x, L) = E0 1+ y e 4LGR cos
π n=1 2n − 1 2L
85
c) y −→ ∞, U (∞, t) = 0
Its initial conditions are
a) t = 0, U (y, 0) = 0
The equation of motion U given by (i) which on taking laplace transforms
d2 u s d2 u s
= u or − u=0
dy 2 γ dy 2 γ
√s √s
y − y
u = Ae γ + Be γ
ˆ +∞
δ d
F {Ux } = U (x, t) e−ist dt = F [U (x, t)]
−∞ δx dx
ˆ +∞
δ2 −ist d2
F {Uxx } = U (x, t) e dt = F [U (x, t)]
−∞ δx2 dx2
86
and
ˆ +∞
δ −iγt
F {Ut } = e dt
−∞ δt
ˆ +∞
= e−iγt |+∞
−∞ +ir U e−iγt dt
−∞
If U (±∞) = 0
87