One Dimensional Wave Equation
One Dimensional Wave Equation
CHAPTER THIRTEEN
and
Partial differential equations arise in connection with various physical
on two or more independent
geometrical problems when the functions involved depend
heat
variables. Most problems in fluid and solid mechanics (dynamics, elasticity),
areas of physics lead to
transfer, electromagnetic theory, quantum mechanics and other
be time and one or
partial differential equations. The independent variables involved maysome of the most
several co-ordinates in space. The present chapter is devoted
to
We derive
important partial differential equations occurring in engineering applications.
consider methods for obtaining
these equations as models of physical systems and
situations.
utions of those equations corresponding to the given physical
differential equation in a region R contains
Since the general solution of a partial
the unique solution of a partial differential
bitrary constants or arbitrary functions,
will satisfy certain other conditions at the
equation corresponding to a physical problem When these
are known as boundary conditions.
boundary of the region R. These
conditions are specified for the time t 0, they are known as initial conditions. A partial
=
(13.1)
Engineering Mathematics - Ill (Mech.) 13.2 Applications of Partial Differential Equations
0Z = 0
and 2 (y-b) +2 (z- c)oy
Differentiating first equation partially w.r.t. x and second equation partially w.r.t. y,
we get
1 -(Z -c)ax
and
1+ (-
i.e. a partial differential equation of order two obtained by eliminating four arbitrary
constants a, b, c and d.
z = xf
(ii)
Ox )
and
oy
+ y oy xf 2.
Engineering Mathematics Il (Mech.) 13.3 Applications of Partial Differential Equations
i.e. a partial differential equation of order one obtained by climinating one arbitrary
function f. Note that it is same as obtained in (i).
(iv) z
=f(y+ ax) + (y ax)
OZ
a f (y+ax) - a $'(y-ax)
and
f'(y+
oy ax) +®(y- ax)
af" (y+ax) +a2 ¢" (y -ax)
Ox2
i.e. a partial differential equation of order two obtained by eliminating two arbitrary
functions fand o.
(v) z ax + by + f
x a +f"|:
and b+r
x + Y ayy = ax + by
xx2 + y =a
and
Xoy x y +yb
and second by y and adding, we get
Multiplying first by x
+xy y t +y =0
two arbitrary constants a and b and
i.e. a PDE of order two obtained by eliminating
arbitrary function f.
Note: We observe that there is relationship between the order of PDE and the
some
XX+Ax
T cos a =
T2 cos ß = T = constant . (1)
n vertical direction, have two forces, the vertical components T, sin a and
we -
where, m= mass of the undeflected string per unit length, öx = the length of the portion
of the
undeflected string.
T sin B T , sin a = m ôx.
tantan a =
and tan B =
Note: Here
x)+8x
we write partial derivatives because also
u depends on t.
u (u
Ox) o'u
M Ox
As dx 0,
Let c2 T
m
Ou
L.C.:ux,0) =
f(x) ,
l g . (3)
Our problem is now to find a solution of () the conditions (2) and (3).
satisfying
We take the following important steps :
Step 3 Using Fourier series, we shall compose those solutions, in order to get a
solution of the wave equation (1) that also satisfies the given initial conditions.
Let us consider these steps one by one.
Step 1: Two ordinary differential equations: In the method of separating variables
(or product method), we determine solution of the wave equation (1) of the form
But ou Ou
O2
d'G = c2 d2 G ()
F(x).d2
FG" = c2 F" G
G"
c2G F
G"
i.e. F C2G
Here L.H.S. is a function of x alone and R.H.S. is a function of t alone, hence above
two expressions are independent of each other. Therefore we can equate to any constant,
say k.
Engineering Mathematics Ill (Mech.) 13.8 Applications of Partial Dilferential Equations
G"
Also, c2 k G"-c? kG = 0 - c2k G(t) =0
Step 2: Satisfyingthe boundaary conditions:The boundary conditions are
u (0, t) =
0,u(l,t) = 0 for allt.
But u (x, t) = F(x) G (t)
u (0. t) = 0 F(0) G () = 0
u(, t) = 0 F() G ()=0
If G () = 0 then u = 0 is a trivial solution. Thus G() + 0, and hence F(0) = F() = 0
F" 0 d'F 0
dx2=
whose solution is F(x) =
c, X
+ C2.
F(0) = 0 c2 = 0 ; F() = 0 c, = 0
GG" =m
F m2 F-m F = 0
d?F
dx2-n
F0 let D dx
DF-m2F = 0, (D-m?) F =0
A.E. D - m2 = 0 D=tm
Now F (0) = 0 0
C +C2 =
F) = 0 C em+ C2 em = 0
m2 F"+ m2 F
=
=0
-
dF +m2F 0 ddx
dx2
dx
=
let D
D'F+ m2F =
0, (D' + m?) F =0
A.E. D+ m2 0 =
D =t im
F (x) = c, cos mx +C2 Sin mx.
Now F (0) 0 =
c, =0
F (x) = C2 sin mx,
m =
We thus obtain infinitely many solutions F (x) =F, (x), where
G"
Also, c2G - m? where m=T
G" +c* m G = 0 +c m2 G =0
dt2
d
Let D = du D:G + c m2 G =0 (D'+c* m) G = 0
G (t) =
c3 cos cmt + C Sin cmt
nTct TCt
Or G(t) =
G, (t) =
Cj cos +C sin m
u (x. t) = F(x) G() =F% (x) Ga() = U, (x, 1 ) . (say)
Now,
u,x. ) =
c, sin c, cos
nTCt
+C sin nTct
nict nitX
u, (x, ) : A,
= cos +B sin Sin
B, c: Ci , n =
1,2,3,...
where, A =
C2 Ci, =
Engineering Mathematics -II (Mech.) 13.1 Applications of Partial Differential Equations
evident that no single solution of
Step 3: Solution of the entire problem It is
P.D.E. is linear, principle of
u,(x, t) can satisfy the initial conditions. However, the given
several solutions then their sum
superimposition is valid meaning thereby that if, we have
is also a solution.
nTct nTct
u(x, t) = 2
n =1|
A, cos B sin sin T (5)
A =Sro)sind
To determine B, we have
ntX
g)= B . sin
n=l
We must choose B, so that becomes the half range Fourier sine series
of g(x).
B =
Sg. sin d
B, g).sindx
Substituting values of A, and B, in (5), we get the required solution.
Engineening Mathematics -ll (Mech.) 13.11 Applications of Partial Differential Equations
Result: Our discussion shows that u (x, t) given by (5) with cocfficients A, and B,
the
is a solution of (1) that satisfies all the conditions (2) and (3) of our problem, provided
series (5) converges.
Summary: To solve the one-dimensional wave equation,
@'u
subjectsubjectto the conditions
1. u (0, t) = 0
Boundary conditions
u(1,t) = 0J
3 u (x, 0) = f(x)
Initial conditions
g (x)J
ot=0
The most general solution is given by
u (x, t) = (C Cos mx +C Ssin mx) (C, cos cmt + c sin cmt)| ..(6)
ILLUSTRATIONS
***N9**************
(ii) y1,1) = 0
and initial conditions,
Here. C* 0, sin mx # 0 = 0
The most general solution will be
y (x, t)= (C2 sin mx) (C3 cos cmt)
Cs sin
.. (1)
y (x, t) = mx cos cmt
Now, Cs0
Cs 0 (otherwise y (x, t) =0 will become trivial solution.)
COS cmt # 0
m nt n = 1,2,3
sin ml = 0 ml = n t, .
x*), 0SxSl
niTx
k (lx - x2) = 2 b, sin . 0SxSl
n = 1|
This is Fourier's half range sine series for f(x) = k{lx - x?) in 0 Sxsl.
ntX niTX
b,= sin =j f(x) dx J k(lx -x)
0
sin dx
2 - ) cos-a-2) nTX
nTx
+-2)i
T3co
COSs
2k2
I n T) cos nTX
Engineering Mathematics - I (Mech.) 13.13 Applications of Partial Differential Equations
b, kP
t)
4kP y E I ) nTct
y (x, =
sin
n=l
cos . Ans.
OO
y (x, t) =
8kP 1 sinn mX COS
(2n-1)nct
(2n 1)3
Ex. 2: A string is stretched and fastened to two points l apart. Motion is started by
displacing the string in the form u = a sin from which it is released at time t = 0.
u
Findthe displacement u (x, t) from one end. Use wave equationp
(ii) u (l, t) 0, Vt =
(Ou = 0
ii)
TTX
(iv) u
(x, 0) = a sin
Condition (i) c = 0
Condition (ii) c = 0.
Engineering Mathematics -!ll (Mech.) 13.14 Applications of Partial Differential Equations
Solution becomes
u (x, t) (C, sin mx) (C, cOs cmt)
u (x, t) = Cs sin mx cos cmt ..(1)
Applying condition (ii), we get,
0 = Cs sin ml cos cmt
.. sin ml = 0
Cs0, cos cmt # 0
ml = T m = n =
l1, 2, 3, . . .
cos nTct|
nTX
u (x, t) = 2 b, sin . (2)
nl
(x, 0)
TLX
u = a
sin
TX n+TX
a sin 2 b, sin
n= 1
. b =a, b, = 0 =
b, =
... =
ba =
( 2 ) will become
TX TCt 2TX
u (x, t) = b, sin cos + b, sin cos
TTX Tct|
u (x, t) = a sin T COS
differential equation Ox? = c? Og prove that the oscillations of the string are
given by
y = a
sec cos-2sinp. (May 2005)
Engineering Mathematics Ili (Mech.) 13.15 Applications of Partial Differential Equations
Sol.: We have the G.S..
y = (C, cos mx + C2 Sin mx) (C, cos cmt + C^ sin cmt) .. )
The condition y = a sin pt for x 0
gives =
a sin
pt =
C, (C3 cos cmt + C Sin cmnt)
This implies that c, Cz = 0 ... (ii)
CC a
(111)
and cm =
p . (iv)
From ii), (ii), c^ =0 and from (iv), m =
Substituting in (i), we get
sin px
pX
y =
c, cos+ c
c sin pt
= px
C C cos+C C, sinsinpt
y a cos CC, sinsin pt . (V)
The condition y = a sin pt for x = l gives
p
a sin pt lacos+, c, sinpsin pt|
a-cos=c;e sin C
pl
2sin 2c sin 2c
C2 C aa = a
2 sincos cos 2c
p
sin 2c
y = a cos+a sin pt
cos 2c
pX
COS
cos2sinsin
a sin pt
cOs 2c
y =a sec cos s i n pt
Applications of Partial Differential Equations
Engineering Mathematics -ll(Mech.) 13.1
Ex. 4:A rightly stretched string with fixed end points x = 0 and x = l is initially ina
(iv) y (x, 0) =
yo sin'
The most general solution is given by:
y (x, t) (C, cos mx + C2 sin mx) (C3
= cos cnmt + C4 Sin cmt)
Condition (i) c, = 0
Condition (iii) > c, = 0
The most general solution will become
m n
=1, 2, .
Solution () becomes
NTct
y (x, t) = Cs sin cos n =
1,2, ...
Combining all these solutions, we get
nTct|
y(x, t) =2 b, sinTcOs (
by using sin' 0 =
sin0- sin3e
JTX
sin sin=b,sin b, sin+b sin7
b, =0; b^=-b=0 =b, =
b, = =
b, =
Comparing we get, b,
. . .
X TCt 3Tct
y (x, t) =
4sin cos sincos
is stretched between two fixed points at a distance l apart,
Ex. 5: An elastic string
and distance from this end the string is displaced a
one end is taken at the origin at a
Slope of OB = 21/73
a 3a
Sol. 21
3ax
Equation of OB is y = i
Slope of BA = 2
X
A4,0)
(0,0)
Fig. 13.2
Equation of BA is y
-0 = -
(x-)
y -x)
We have to solve
Engineering Mathematics Ill (Mech) 13.18 Applications of Partial Differential Equations
(ii) = 0
= 0
iv) y (x, 0) = 3a X,
0SxS
-.sxsl
The most general solution is
y (x, t) = (C, cos mx + C2 sin mx) (C; cos cmt + C sin cmt)
Condition (i) C= 0
Condition (ii) C = 0
ntX nTct
y (x, t) =2 b, sin cosT . (1)
- ) . SxS
nTX
b, =7 f(x) sin dx
0
3ax
2 sin dx J(-x) sindx +
Engineering Mathematics ! (Mech.) 13.19 Applications of Partial Differential Equations
n2T2 Sin
nTx
- cos-1) sin
2nT 2nT ZnT 2nt|
P3nn COs
3 +2n T? SIn3 3n COS3+n?T?
SIn
3
nT in 9asin 2nt
Substituting in (I), we get,
ntct
y Cx,t)= 2sinsin COS
EXERCISE 13.1
1. A taut string of a length 21 is fastened at both ends. The mid point of the string is
taken to a height b and then released from rest in that position. Obtain the
displacement.
(oy y(21, t) c
Hint: B.C.: 0
=0, (ii) =
0, Use=
DX0Sx[l
L.C.: (ii) y (0, 1) = 0; (iv) y (x, 0) =
2x-x). isxs2
The most GS is y(x, t) = (Ci cos mx +C2 Sin mx) (¢3 cos cmt + C4 Sin cmt)
2
nTX 8b sin 2 .
and use formula for b, = . f(x) sin 27 dx =
0
8b nT nTx
(nTct
Ans.:
y,) = n=1
sin2 sin2 COS
2
2. If a string of length / is initially at rest in its equilibrium position and each of its
such that
point is given a velocity v (x)
CX, 0<xS
v (x) =
c-) sxsi
Obtain the displacement y (x, t) at any timet. (Nov. 2013)
Engineering Mathematics -Ill (Mech.) 13.20 Applications of Partial Diferential Equations
nmat
Ans.: y x, t) = 2 s i ninsin sin
a ns
3.A tightly stretched string with fixed ends x = 0 and x = l is initially at rest in its
cquilibrium position. Ifit is set vibrating giving each pointa velocity 3x (1- x) for
0<x<I, find the displacement.
Hint: BC. : (i) y (0, t) 0, (ii) y (1, t)
=
0; =
and use bn =
i f(x) sindx =
1 -cos n)
Ans.:: y (x, t) = (Zn 1 ) Sin
(2n 1) TX 2n-1) zct
Sin
nTX
b 7 f(x) sin dx =
n (1-cos nt)
24/ 2 . N sin (2n- 1) Tx
Ans.:
Ans. y n -
1)3 Sin COS
(2n-) rct
D=
Engineening Mathematics -1(Mech) 13.21 Applications of Partial Differential Equations
12800
Ans. y (x, t) =
T 2 (2n s1)3 , 2n-) ax
i SIn COS
(2n-) nct
40
40
7. An clastic string is stretched between two points at a distance l apart. One end is
3d sxsl
TX tat
Ans.: y (x, t) =
yo sin cos
9. The points of trisection of a tightly stretched string of length I with fixed ends are
pulled aside through a distance d on opposite sides of the position of equilibrium
and the string is released from rest. Obtain the displacement of the string and
show that mid-point always remains at rest.
9d
Ans. y (x,) =
sinsin cos 2nTat
when x =
(4d
3 7 X, xs
.C. :(iii) t 0
= 0 iv) y (x, 0) = <
NTX n7Tat
Ans. y (x, t) = sin4 sin cos
n =l
0 and x =l is given the initial
I. A uniform string stretched between the points x =
TTX
sin, 0<x<land initial velocity,
displacement y(x, 0) =
0. 0<x<
v (x) = a,
0,
12. A string of length I fixed at its ends satisfies the wave equation
Ox
Find the solution ifthe string has initial triangular deflection given by
x, 0sxs
y (x, 0) =
d-. sist
and initial velocity zero. (Dec. 2007)
n-1
(2n 1)? Sin COS
13. Find the deflection u x, ) of a vibrating string lengthl= R, ends fixed and c? =
corresponding to zero velocity and initial deflection 0.01 (7-x).
(May 2008)
Hint: c =4. c 2, bi = 0.1, b2 = 0, b3 = 0, ba = 0.01, bs = b6 = ... =0.
Ans. : y (x, t) =
0.1 sin x cos 2t + 0.01 sin 4x cos 8t
13.5 MODELING OF ONE-DIMENSIONAL HEAT FLOW
Derivation of Equation:
We make use of following experimental facts or empirical laws:
()
) Heat flows from higher temperature to lower
temperature.
(i) The rate of flow of heat through an area is proportional to the area and to the
(Ou
temperature gradient in degrees per unit distance | where u(x, t) is
temperature distribution normal to the area. Constant of proportionality is
called the thermal conductivity of the material and denoted generally by k.
ii) The amount of heat required to change the temperature through a given range is
proportional to the mass of the body and the change of temperature. The
constant of proportionality is termed as specific heat and generally denoted
by S.
Consider a homogeneous bar of uniform cross-section, sides coated with insulating
material. It is assumed that the loss of heat from the sides by conduction or radiation is
negligible. One end of the bar is treated as the origin and the direction of heat flow as
positive X-axis. Let p be the density (gm/cm*), 'S" the specific heat (cal/gm deg) and
k the thermal conductivity (cal/cm-deg.sec). The temperature at any point of the bar
depends on the distance x of the point from one end and time t ' and is denoted by u(x, t)
or u. Also the temperature distribution through a cross-section is same.
-öx
X
Fig. 13.3