PHY203
PHY203
136
Course Code PHY 203
Course Title OSCILLATIONS AND WAVES
Course Developer DR. A. B. Adeloye
PHYSICS DEPARTMENT
UNIVERSITY OF LAGOS
137
OSCILLATIONS AND WAVES
Course Introduction
The phenomena we normally observe in nature can be broadly classified into two categories:
those concerned with matter and those concerned with waves. Physics courses usually begin with
discussion of phenomena dealing with matter mechanics and properties of matter. Next comes the
phenomena of waves. Out of our five senses (touch, taste, smell, hearing and seeing), two deal
with the waves—hearing and seeing. Our contact with the external world is mainly through these
two senses. Sound and light, though of entirely different nature, have many properties in
common. In this course, you will learn about waves in general. This unified approach to wave
motion is meant to bring out the underlying similarity between apparently widely differing
phenomena. Even our understanding of modern physics, particularly quantum mechanics,
depends upon a clear understanding of this course.
In Unit 1 of this course we have developed the mathematical theory of simple harmonic motion.
It is used to study oscillations, by analogy, of entirely different systems from different branches
of physics. Unit 2 deals with the superposition of two or more collinear or orthogonal harmonic
oscillations of same/different frequencies.
In nature, most oscillations left to themselves die down gradually. This happens because of
damping. The effect of damping on harmonic oscillations is discussed in Unit 3. In Unit 4 you
will learn about the motion of a damped harmonic oscillator on which a periodic harmonic force
is acting. This leads to the spectacular phenomenon of resonance. Unit 5 deals with the analysis
of coupled oscillations. You will learn that in the limit of their number becoming very large, we
are lead to the phenomenon of wave motion.
In Unit 6, you will learn the basic concepts and vocabulary of wave motion as well as wave
propagation in one and more dimensions. The wave equation for one dimensional progressive
waves in a stretched string as well as fluids (gases and liquids) are established. Its connection
with wave impedance presented by a medium is also discussed for transverse as well as
longitudinal waves.
In Unit 7, you will learn the changes a wave undergoes at the interface of two different media,
using Huygen’s construction and the concept of wave impedance. Expressions for reflection and
transmission amplitude and energy coefficients are derived. In this unit we have also discussed
Doppler’s effect and shock waves.
In Units 8 and 9, you will learn about the superposition of waves. You will study superposition of
waves under different conditions. You will see that the superposition of two waves which have
the same amplitude, frequency and wavelength but are moving in opposite directions result in the
formation of stationary waves. These are basically responsible for production of music. Two
waves having slightly different frequencies but traveling in the same direction give rise to a wave
group and beats. In Unit 9
138
In each unit, we have given many SAQ's and TQ's to' fix-up your ideas. If you are not able to
solve them yourself, you can look for solutions at the end of each unit.
We hope that after studying this Block you will realise the wide applicability of simple harmonic
motion and its connection with wave motion. You are, therefore, expected to master the
mathematical technique needed to study SHM under different conditions.
139
UNIT 1 SIMPLE HARMONIC MOTION
Structure
1.1 Introduction Objectives
1.2 Simple Harmonic Motion (SHM): Basic Characteristics
Oscillations of Spring-mass System ;
1.3 Differential Equation of SHM
1.4 Solution of the Differential Equation for SHM Amplitude and Phase Time Period and
Frequency Velocity and Acceleration
1.5 Transformation of-Energy in Oscillating Systems: Potential and Kinetic Energies
1.6 Calculation of Average Values of Quantities Associated with SH./l
1.7 Examples of Physical Systems Executing SHM
Simple Pendulum
Compound Pendulum
Torsional Systems
An L-C Circuit
An Acoustic Oscillator
A Diatomic Molecule: Two-body Oscillations
1.8 Summary
1.9 Terminal Questions
1.10 Solutions
1.1 INTRODUCTION
In your school science courses you must have learnt about different types of motions. You are
familiar with the motion of falling bodies, planets and satellites. A body 1 released from rest and
falling freely (under the action of gravity) moves along a straight line. But an object dropped
from an aeroplane or a ball thrown up in the air follows a curved path (except when it is thrown
exactly vertically). You must have also observed the motion of the pendulum of a wall clock and
vibrating string of a violin or some other string instrument. These are examples of oscillatory
motion. The simplest kind of oscillatory motion which can be analysed mathematically is the
Simple Harmonic Motion (SHM). We can analyse oscillatory motions of systems of entirely
different physical nature in terms of SHM. For example, the equation of motion that we derive for
a pendulum will be similar to the equation of motion of a charge in a circuit containing an
inductor and a capacitor. The form of solutions of these equations and the time variation of
energy in these systems show remarkable similarities. However, there are many important
phenomena which arise due to superposition of two or more harmonic oscillations. For example,
our eardrum vibrates under a complex combination of harmonic vibrations. But we shall discuss
this aspect in the next unit.
In this unit we will study oscillatory systems using simple mathematical techniques. Our
emphasis would be on highlighting the similarities between different systems.
Objectives
After studying this unit you should be able to:
state the basic criteria for the simple harmonic motion of a system
establish the differential equation for a system executing SHM and solve it
define the terms amplitude, phase and time period
compute potential, kinetic and total energies of a body executing SHM
deduce expressions for average potential and average kinetic energies and discuss their
significance
140
write down the equation of motion and expressions for displacement, time period and
energy of
simple physical systems executing SHM
identify similarities between different oscillating systems.
define wave motion an state its characteristics
distinguish between longitudinal and transverse waves
represent graphically waves at a fixed position or at a fixed time
relate wavelength, frequency and speed of a wave
establish wave equations for longitudinal and transverse waves
compute the energy transported by a progressive wave
derive expressions for velocities of longitudinal and transverse waves
derive expressions for characteristic impedance and acoustic impedance
write two and three dimensional wave equations
141
If the spring is stretched by pulling the mass longitudinally, due to elasticity a restoring force
comes into play which tends to bring the mass back towards the equilibrium position (Fig l.lb). If
the spring were compressed the restoring force would tend to extend the spring and restore the
mass to its equilibrium position (Fig. l.lc).
Fig. 1.1 A Spring-mass System as an ideal oscillator (a) The equilibrium configuration, (b) An
extended configuration, (c) A compressed configuration.
The more you stretch/compress the spring, the more will be the restoring force. So the direction of
the restoring force is always opposite to the displacement. If total change in the length is small
compared to the original length, then the magnitude of restoring force is linearly proportional to
the displacement. Mathematically, we can write
F kx (1.1)
The negative sign signifies that the restoring force opposes the displacement. The quantity k is
called the spring constant or the force constant of the spring. It is numerically equal to the
magnitude of restoring force exerted by the spring for unit extension. Its SI unit is Nm -1.
SAQ 1
The spring in Fig. 1.1 a is stretched by 5 cm when a force of 2 N is applied. Calculate the spring
constant. How much will this spring be compressed by a force of 2.5 N?
How does the spring-mass system oscillate? To answer this question, we note that when we pull
the mass, the spring is stretched. The restoring force tends to bring the mass back to its
equilibrium position ( x = 0). Therefore, on being released, the mass moves towards the
equilibrium position. In this process it acquires kinetic energy and overshoots the equilibrium
position. Do you know why? It is because of inertia. Once it overshoots and moves to the other
142
side, the spring is compressed and the mass is acted upon by a restoring force but in the opposite
direction. Thus we can say that oscillatory motion results from two intrinsic properties of the
system: (i) elasticity and (ii) inertia.
What is the direction of the restoring force vis-a-vis the equilibrium position of an oscillating
body?
The restoring force is always directed towards the equilibrium position of the oscillating body.
In discussing the spring-mass system we observed two important points:
(i) The restoring force is linearly proportional to the displacement of mass from its
equilibrium position.
(ii) The restoring force is always directed towards the equilibrium position.
Any oscillatory motion which satisfies both these conditions is called simple harmonic motion.
The study of SHM is important because, as you will see, oscillatory motion of systems of entirely
different physical nature can be analysed in terms of it.
Let us now study the effect of gravity on oscillations of spring-mass system. Consider a spring of
negligible mass suspended from a rigid support with a mass m attached to its lower end (Fig.1.2).
Fig 1.2: A vertically hanging spring-mats system, (a) The sprint with no object suspended from
it, (b) The spring in equilibrium with mass m suspended, (c) Spring-mass system
displaced from equilibrium position.
Let us choose the X-axis along the length of the spring. We take the bottom of the spring as our
reference point, X = 0, when no wejght is attached to it (Fig. l.2a). When a mass m is suspended
from the spring, let the reference point move to X=X0 (Fig. 1.2b). At equilibrium, the weight, mg,
balances the spring force, kX0.. Since the net force is zero, we have
mg kX 0 0
or
mg kX 0 (1.2)
143
Now if the mass is pulled downwards so that the reference mark shifts to X1 (Fig 1.2c), then the
total restoring force will be kX1 and points in the upward direction. The net downward force will
therefore be (using Eq. (1.2)),
mg kX1 k ( X 0 X 1 ) kx
where x X 1 X 0 .
F kx
where x is its displacement from the equilibrium position, X 0 . This result is of the same form as
Eq. (1.1) for the horizontal arrangement. It is thus clear that gravity has no effect on the frequency
of oscillations of a mass hanging vertically from a spring; it only displaces the equilibrium.
Let us now find the differential equation which describes the oscillatory motion of a spring-mass
system. The equation of motion of such a system is given by equating the two forces acting on the
mass:
It is important to note that in this equation, the equilibrium position of the body is taken as the
origin, x = 0.
You will note that the quantity k/m has units of Nm 1kg 1 (kg.mg 2 )kg 1m 1 s 2 . Hence
we can replace k/m by 0 where 0 is called the angular frequency of the oscillatory motion.
2
It may be remarked here that Eq. (1.3) is the differentia] form of Eq. (1.1) and describes simple
harmonic motion in one dimension.
A differential equation having terms involving only the first power of the variable and its
derivatives is known as a linear differential equation. If such an equation contains no term
independent of the variable it is said to be homogeneous. We may, therefore, say that Eq. (1.3) is
a second order linear homogeneous equation. Its solution will contain two arbitrary constants.
144
1.4 SOLUTION OF THE DIFFERENTIAL EQUATION
FOR SHM
To find the displacement of the mass at any time t, we have to solve Eq. (1.3) subject to given
initial conditions. A close inspection of Eq. (1.3) shows that x should be such a function that its
second derivative with respect to time is the negative of the function itself, except for a
multiplying factor o>o. From elementary calculus, we know that sine and cosine functions have
this property.
You can check that this property does not change even if sine and cosine functions have a
constant multiplying factor.
A general solution for x (t) can thus be expressed as a linear combination of both sine and cosine
terms, i.e.,
d (sin t ) d 2 (sin t )
cos t and 2
2 sin t
dt dt
Similarly,
d (cos t ) d 2 (cos t )
sin t and 2
2 cos t
dt dt
Differentiating this equation twice with respect to time and comparing the resultant expression
with Eq. (1.3), we obtain 0 . The negative sign is dropped as it gives negative frequency
which is a physically absurd quantity.
The constants A and occurring in Eq. (1.5) are determined using the initial conditions on
dx
displacement (x) and velocity .
dt
145
Let us assume that the mass is held steady at some distance a from the equilibrium position and
then released at t = 0. Thus the initial conditions are: at t = 0, x = a
dx
and = 0. Then, from Eq. (1.5) we would have
dt
x (at t = 0) = A cos a
and
dx
(at t = 0) = A 0 sin 0
dt
These conditions are sufficient to fix A and . The second condition tells us that is either zero
or n (n = 1, 2, ...). We reject the second option because the first condition requires cos to be
positive. Thus with the above initial conditions, Eq. (1.5) has the simple form
x a cos t (1.6)
SAQ2
Take A1 B sin and A2 B cos in Eq.(1.4). In this case show that the solution is
We therefore observe that both cosine and sine forms are valid solutions of Eq. (1.3). If you plot
Eq. (1.5), the graph will be a cosine curve with a definite initial phase (Fig. 1.3).
Fig. 1.3: Displacement-time graph of simple harmonic motion with an initial phase
146
x0 a cos
We know that the value of the sine and cosine .functions lies between 1 and 1 . When
cos( 0 t ) = 1 or 1, the displacement has the maximum value. Let us denote it by a or a .
The quantity a is called the amplitude of oscillation.
The displacement-time graphs for = 0, / 2 and are shown in Fig.1.4. In all the cases, the
graphs have exactly the same shape if we shift the origin along the time axis. When the phase
difference is two oscillations are said to be in opposite phase or out of phase by .
The mass in Fig .1.1 oscillates with an amplitude a. If the time is measured from the instant when
it is at (i) x = + a, (ii) x = a , and (iii) x a / 2 , calculate the phase constant for the
equations (a) x a sin( 0 t ) and (b) x a cos( 0 t ) .
That is, the displacement of the particle repeats itself after an interval of time 2 / 0 . In other
words, the oscillating particle completes one vibration in time 2 / 0 . This time is called the
period of vibration or the time period. We denote it by T :
147
T 2 / 0 (1.8)
T 2 m / k (1.9)
The number of vibrations executed bv the oscillator per second is called the frequency. The unit
of frequency is Hertz (Hz). Denoting it 0 , we have for a spring-mass system
1 1 k
0 (1.10)
T 2 m
This means that stiffer the spring, higher will be the frequency of vibration.
Therefore, the instantaneous velocity, which is the lust time derivative of displacement, is given
by
dx
v 0 a sin( 0 t ) (1.11)
dt
We can rewrite it as
v 0 a cos( / 2 0 t ) (1.12a)
You may also like to know the value of v at any point .v. To this end, we rewrite Eq. (1.11) as
We also know that acceleration is the first time derivative of velocity. From Eq. (1.11)
dv
0 a cos( 0 t )
2
dt
0 a cos( 0 t )
2
(1.13a)
148
If you compare Eqs. (1.7), (1.12 a) and (1.13 a), you will note that (i) 0 a is the velocity
amplitude and 0 a is the acceleration amplitude, and (ii) velocity is ahead of displacement by
2
If you plot displacement, velocity, and acceleration as functions of time, you will get graphs as
shown in Fig.1.5.
Fig. 1.5: Time variation of displacement, velocity and acceleration of a body executing SHM
( = 0)
SAQ4
The displacement of a particle executing simple harmonic motion is given hy
x 0.01cos 4 (t 0.0625) metre. Deduce (i) the amplitude, (ii) the time-period, (iii) maximum
speed, (iv) maximum acceleration and (v) initial displacement.
Consider the spring-mass system shown in Fig 1.1. When the mass is pulled, the spring is
elongated. The amount of energy required to elongate the spring through a distance dx is equal to
the work done in bringing about this change. It is given by dW dU F0 dx , where F0 is the
applied force (such as by hand). This force is balanced by the restoring force. That is, its
magnitude is same as that of F and we can write F0 kx . Therefore, the energy required to
elongate the spring through a distance x is
x x
1 2
U F0 dx k xdx kx (1.14)
0 0 2
This energy is stored in the spring in the form of potential energy and is responsible for
oscillations of the spring-mass system.
On substituting for the displacement from Eq. (1.7) in Eq. (1.14), we get
149
1 2
U ka cos( 0 t ) (1.15)
2
1 2
U0 ka cos 2 (1.16)
2
As the mass is released, it moves towards the equilibrium position and the potential energy starts
1 2
changing into kinetic energy (K.E). The kinetic energy at any time / is given by K .E. mv .
2
Using Eq. (1.11), we get
1
K .E. m 0 sin 2 ( 0 t )
2
2
1
= ka 2 sin 2 ( 0 t ) (1.17)
2
since 0 k /m.
2
One can also express K.E. in terms of the displacement by writing K .E. in terms of the
displacement by writing
1 2
K .E. ka [1 cos 2 ( 0 t )]
2
1 1
ka2 ka2 cos 2 ( 0 t ) (1.17)
2 2
1 1 1
ka2 kx2 k (a 2 x 2 ) (1.18)
2 2 2
This shows that when an oscillating body passes through the equilibrium position ( x 0) , its
1 2
kinetic energy is maximum and equal to ka .
2
SAQ5
Show that the periods of potential and kinetic energies are one-half of the period of vibration.
It is thus clear from the explicit time dependence of Eqs. (1.15) and (1.17) that in a spring-mass
system the mass and the spring alternately exchange energy. Let us consider that the initial phase
= 0. At t = 0, the potential energy stored in the spring is maximum and K.E. of the mass is
zero. At t = T/4, the potential energy is zero and K. E. is maximum. As the mass oscillates,
energy oscillates from kinetic form to potential form and vice versa. At any instant, the total
energy, E, of the oscillator will be sum of both these energies. Hence, from Eqs. (1.15) and
(1.17), we can write
1 2 1 1
E U K .E. ka cos 2 ( 0 t ) ka2 sin 2 ( 0 t ) ka2 (1.19)
2 2 2
150
This means that the total energy remains constant (independent of time) and is proportional to the
square of the amplitude. As long as there are no dissipative forces like friction, the total
mechanical energy will be conserved.
Fig. 1.6 Variation of potential energy (U), kinetic energy (K.E) and total energy (E) with
displacement according to Eqs. (1.14), (1.18) and (1.19)
The plots of U and K.E. as a function of x as obtained from Eqs. (1.14) and (1.18) are shown in
Fig 1.6. You will note that
(i) the shape of these curves is parabolic, (ii) the shape is symmetric about the origin, and (iii) the
potential and kinetic energy curves are inverted with respect to one another. Why? This is due to
the phase difference of / 2 between the displacement and velocity of a harmonic oscillator. At
1 2
any value of x, the total energy is the sum of kinetic and potential energies and is equal to ka .
x 2
This is represented by the horizontal line.
The points where this horizontal line intersects the potential energy curve are called the 'turning
points.' The oscillating particle cannot go beyond these and turns back towards the equilibrium
1 2
position At these points, the total energy of the oscillator is entirely potential ( E U ka )
2
and K.E. is zero. At the equilibrium position (x = 0) the energy is entirely kinetic
1 2 1 2
( K .E. E ka ) so that the maximum speed, v max is given by the relation mvmax E,
2 2
i.e., vmax 2E / m .
At any intermediate position, energy is partly kinetic and partly potential, but the total energy
always remains the same. The transformation of energy in a spring-mass system is shown in Fig
1.7.
151
Fig 1.7 Energy transformation in a spring-mass system at various times. The bars indicate
potential and kinetic energies are shown at intervals of t T / 8 .
Do you know why the point of minimum potential energy is regarded as the position of stable
equilibrium? This is because there is no net force acting on the system in this position.
SAQ 6
A body of mass m fell from height h onto the pan of a spring balance. The masses of the pan and
the spring are negligible. The stiffness constant of the spring is k. Having stuck to the pan, the
body executes harmonic oscillations in the vertical direction. Find the amplitude and the energy
of oscillation.
The time average of kinetic energy over one complete cycle is defined as
152
T
K .E.dt
K .E. 0
(l.20a)
T
ka2
sin ( t )dt
T
K .E. 2
0 (1.20b)
2T 0
On solving the integral in Eq. (l.20b) you will find that its value is T / 2 . So, the expression for
average kinetic energy reduces to
ka 2
K .E. (1.21)
4
Similarly, one can show that the average value of potential energy over one cycle is
ka 2
U (1.22)
4
That is, the average kinetic energy of a harmonic oscillator is equal to the average potential
energy over one complete period.
Thus the sum of average kinetic and average potential energies is equal to the total energy :
1 2 1 2
K.E. U = ka ka
4 4
1 2
= ka E
2
F mg sin
153
Fig. 1.8 A simple pendulum
The bob is moving along the arc whose length at any instant is given by x. If the corresponding
angular displacement from the equilibrium position is , then the length of arc is
x l (1.24)
where l is length of the string by which the bob is suspended. Differentiating Eq. (1.24) twice
with respect to t, and substituting the result in Eq. (1.23), we get
d 2 g
2
sin (1.25)
dt l
For small angular displacements, sin may be approximated to . In this approximation, Eq.
(1.25) takes the form
d 2
2
02 0 (1.26)
dt
where 0 g /l .
Eq. (1.26) is exactly of the standard form (1.3) showing that pendulum executes simple harmonic
motion. The time period of oscillation is given by
2
T 2 l / g (1.27)
0
By analogy, we can write the general solution of the Eq. (1.26) as
m cos( 0 t ) (1.28)
154
where is the maximum angular displacement.
From Eq. (1.27) you will note that for small angular displacements, the frequency of oscillation of
a simple pendulum depends on g and but not on the mass of the bob. The appearance of the
factor g in Eq. (1.27) implies that a pendulum clock will move slower near the equator than at the
poles. Do you know why? This is because the value of g varies with latitude. For the same reason,
the period of a pendulum will be different on moons and planets.
When the amplitude of oscillation is not small, we are required to solve the general Eq. (1.25).
The time period, which can be expressed in the form of a series involving the maximum angular
displacement m is given by
l 1 1 32
T 2 1 2 sin 2 m 2 2 sin 4 m ...
g 2 2 2 4 2
You can check the accuracy of Eq. (1.27) by comparing the value of T obtained from Eq. (1.29).
For example, you will find that when m is 15° (corresponding to a total to and fro angular
displacement of 30°), the actual value of time period differs from that given by Eq. (1.27) by less
than 0.5%.
SAQ7
Use the principle of conservation of energy to show that the angular speed of a simple pendulum
is given by,
1/ 2
If I is the moment of inertia of the body about the horizontal axis passing through S, the restoring
torque equals Id 2 / dt 2 . Hence, the equation of motion can be written as
155
Fig. 1.9 A rigid body oscillating about a horizontal axis: Compound pendulum.
d 2
I mgl sin (1.30)
dt 2
For small angular displacement, sin and Eq. (1.30) takes the form
d 2 mgl
0 (1.31)
dt 2 I
This equation shows that a compound pendulum executes SHM and the time period is given by
2
T 2 I / mgl (1.32)
0
There is a very useful and important theorem of parallel axes in the study of moment of inertia.
According to this theorem, the moment of inertia I of a body about any axis and its inertia I g
about a parallel axis passing through its center of gravity are connected by the relation
I I g ml 2 (1.33)
where l is the distance between the two axes and I g mkr2 . The quantity kr is the radius of
gyration of the body about the axis passing through G. It is the radial distance at which the whole
mass of the body could be placed without any change in the moment of inertia of the body about
that axis.
On substituting the expression for I from Eq. (1.33) in Eq. (1.32), we obtain
k r2 l 2
T 2 (1.34)
gl
156
On comparing this expression for T with that given by Eq. (1.27) for a simple pendulum, you will
note that two periods become equal if l in Eq. (1.27) is replaced by L (k r2 / l ) l . This is
called the length of an equivalent simple pendulum. If we produce the line SG and take a point O
k r2
on it such that SO = l , then O is called the centre of oscillation.
l
When a torsional system is twisted and then left free, it executes torsional oscillations in a
horizontal plane. For an angular displacement , the restoring torque is k t . Here, k t is a
constant which depends on the properties of the wire.
If I is the moment of inertia of the system about the axis of rotation and d 2 / dt 2 is the angular
acceleration, then the equation of motion is
d 2
I k t
dt 2
d 2
or 2
02 (1.35)
dt
157
where 0 k t / l . This is exactly of the standard form (1.3). Hence, the motion is SHM and
the period of oscillation is
T 2 I / k t (1.35a)
You will note that this expression for T contains no approximation. This means that the time
period for large amplitude oscillations will also remain the same, provided the elastic limit of the
suspension wire is not exceeded. The solution of Eq. (1.35) is given by Eq. (1.28).
SAQ8
A solid sphere of mass 3 kg and radius 0.01 m suspended from a wire. Find the period of
oscillations, if the torque required to twist the wire is 0.04 N-m rad-1. The moment of inertia of a
sphere about an axis passing through its centre is given by
2 2
I mr
5
In a pendulum, the mean position is taken as the equilibrium state. What is the equilibrium state
in an L-C circuit? It corresponds to the state when there is no current in the circuit. It may be
disturbed by charging or discharging the capacitor. Let the capacitor be given a charge Q0
coulomb. Then the voltage across the capacitor plates will be Qo/C. Now if the circuit is
disconnected, the capacitor discharges through the inductor. As a result current starts building up
in the circuit gradually and the charge on the plates of the capacitor decreases. At any time t , let
the current in the circuit be q and the charge on capacitor plates be q. Then the voltage drop
across the inductor will be
158
dI
VL L
dt
This must be equal to the voltage Vc q / C across the capacitor plates at that time. Thus, we
can write
Vc VL
q dI
or L (1.36)
C dt
dI d 2 q
Since I dq / dt and , Eq. (1.36) takes the form
dt dt 2
d 2q
2
02 q 0 (1.37)
dt
1
where 02 . This means that one can have a wide range of frequencies by changing the
LC
values of L and C. That is how you tune different stations in your radio sets.
q q0 cos( 0 t ) (1.38)
T 2 LC (1.39)
I Q0 0 sin( 0 t )
I 0 cos( 0 t / 2)
where I 0 Q0 0 .
Thus the current leads the charge in phase by / 2 . In practice, you will always find that an
inductor offers some resistance in an L-C circuit. Its effect on charge oscillations will be
discussed in Unit 3.
Let us now calculate the energy stored in the inductor L and the capacitor C at any instant t. As
the current rises from zero to I in time t, the energy stored in the inductor, EL, is obtained by
integrating the instantaneous power with respect to time,
i.e.
159
t
E L IVL dt
0
The negative sign implies that work is done against, rather than by the emf. On substituting for
VL, we get
I
dI 1
EL L Idt LI 2
0 dt 2
EC q 2 / 2C
1 2 1 q2
E E L EC LI (1.40)
2 2 C
This expression, for total energy is similar to the one for mechanical oscillator
1 2 1 2
(E mv kx ) . As q and I vary with time, the inductor and capacitor exchange energy
2 2
periodically. This is similar to the energy exchange in the spring-mass system. Further, the mass
and inductor play analogous roles in mechanical and electrical systems, respectively.
SAQ 9
Calculate the frequency of electrical oscillations when an inductor of 20 mH is connected with a
capacitor of 1 F . If the maximum potential difference across the capacitor is 10 V, calculate the
energy of oscillation.
Fig. 1.12 (a) An Acoustic oscillator, (b) As air in the neck is pushed, air in the flask is
compressed, and (c) Due to elasticity, air in the flask exerts a restoring force on the air
in the neck.
160
We consider free vibrations of air in the neck of the flask. As the air in the neck moves in, the air
in the flask is compressed. If air in the neck goes out, the air in the flask is rarefied. So the air in
the neck behaves like the mass and the air in the flask behaves like the spring in a mechanical
oscillator.
Suppose that the air in the neck moves inward through a distance x from the equilibrium position.
The change in the volume of the air in the bulb V = x A. Let the increase in pressure over the
atmospheric pressure be p . We know that the volume of a gas depends on the pressure as well
as the temperature. Therefore, the pressure changes in acoustic vibrations should alternately heat
and cool the air in the flask as it gets compressed and rarefied. We assume that the pressure
changes are so rapid that they do not permit any exchange of heat. That is, the process is
adiabatic. Hence, we can write
V Ax
p E E (1.41)
V V
where E is the adiabatic elasticity of the gas. It is defined as the ratio of the stress to volume
strain. Numerically, stress is same as pressure. So we can write
p
E
(V / V )
The negative sign signifies the fact that as pressure increases, volume decreases and vice-versa.
This excess pressure p of air inside the bulb provides the restoring force F, which acts upward.
We can therefore write
E A 2
F pA x
V
If is the density of air, the mass of the air in the neck m lA . Hence, the equation of motion
of air in the neck can be written as
d 2x E A 2
lA 2 x
dt V
d 2 x E A
or x0 (1.42)
dt 2 Vl
This equation has the standard form for simple harmonic motion. Hence, the frequency of
oscillation of air in the neck is
1 E A vs A
0
2 vl 2 vl
161
where v s = E / is the speed of sound. We know that v s is proportional to the square root of
temperature. So the frequency of vibration of air in a flask is also proportional to the square root
of the temperature.
SAQ 10
A flask has a neck of radius 1 cm and length 10 cm. If the capacity of the flask is 2 litres,
determine the frequency at which the system will resonate (speed of sound in air = 350 ms-1).
Suppose that two masses m1 and m 2 are connected by a spring of force constant k. The masses
are constrained to oscillate along the axis of the spring (Fig. 1.13a). Let r0 be the normal length
of the spring. We choose X -axis along the line joining the two masses. If X 1 and X 2 are the
coordinates of the two ends of the spring at time t, the change in length is given by
x ( X 2 X 1 ) r0 (1.44)
(b) Fig. 1.13 (a) A two-body oscillator (b) An equivalent one-body oscillator.
For x > 0, x = 0 and x < 0, the spring is extended, normal and compressed respectively. Suppose
that at a given instant of time the spring is extended, i.e. x > 0. Though the spring exerts the same
force (kx) on the two masses, the force F1 (= kx) acting on m1 opposes the force F2 ( kx) on
m2 , i.e.,
F1 kx and F2 kx
162
d 2 X1
m1 kx
dt 2
and
d2X2
m2 kx
dt 2
Since r0 denotes a constant length of the spring, Eq. (1.44) tells us that
d 2x k
x0 (1.46)
dt 2
1
1 1 m1 m2
where is called the reduced mass of the system.
m1 m2 m1 m2
1
0 k/ (1.47)
2
This means that a diatomic molecule behaves as a single object of mass , connected by a spring
of force constant k (Fig. 1.13b).
SAQ 11
For an HCl molecule, r0 = 1.3Å. Find the value of the force constant and the frequency of
1
oscillation. Given: m H = 1.67 10 27 kg and mcl 35mH . Use 9 10 9 Nm 2 C 2 .
4 0
1.7 SUMMARY
1. Simple Harmonic Motion: An oscillatory motion is said to be simple harmonic when
the acceleration is proportional to the displacement and is always directed against it.
163
We can also say that in SHM the restoring force is linearly proportional to the
displacement and acts against it.
2. Differential Equation of SHM is
d 2x
2
02 x 0 , where 0 = k/m
dt
4. The period and frequency characterizing a SHM are represented by the relations:
2 0 1
T and 0
0 2 T
1 2
The time averaged kinetic energy and potential energy are the same; equal to ka .
4
164
(a) 2 m /( k1 k 2 )
(b) 2 m /( k1 k 2 )
(c) 2 m /(1 / k1 1 / k 2 )
2. A smooth tunnel is bored through the earth along one of its diameters and a ball is
dropped into it. Show that the ball will execute simple harmonic motion with period
T 2 R / g where R is the radius of the earth and g is the acceleration due to
gravity at the surface of the earth. Assume the earth to be a homogeneous sphere of
uniform density.
3. Find the angular frequency and the amplitude of harmonic oscillations of a particle if
at distances x1 and x 2 from the equilibrium position its velocity equals v1 and v 2
respectively.
4. Show that the centers of suspension and oscillations in a compound pendulum are
mutually interchangeable.
e2 c
U (r ) 9
4 r 0 r
The first term represents the attractive part and the second term represents the
repulsive part. Show that the force constant is 2e 2 / 0 r0 , where r0 is the
3
equilibrium separation.
dU
Hint: F kr and F .
dr
165
1.9 SOLUTIONS
SAQ 1
Force 2.0 N
Force constant k = 2
40 Nm 1
Displacement 5.0 10 m
2.5 N
Compressed length = Force / Force constant = 1
6.3 10 3 m
40 Nm
SAQ 2
Putting A1 B sin and A2 B cos in equation (1.4), we get x(t ) B sin(t )
since sin( A B) sin A cos B sin B cos A
SAQ 3
(a) x a sin( 0 t )
Since time is measured from the instant x a , we get
(i) x a a sin , i.e., sin = 1 or / 2
Similarly,
and
a 1
(iii) x a sin , i.e., sin = or / 4
2 2
(b) x a cos( 0 t )
At t = 0,
(i) x a a cos , i.e., cos 1 or 0
(ii) x a a cos , i.e., cos 1 or
a 1
(iii) x a cos , i.e., cos or / 4
2 2
SAQ 4
x 0.01cos 4 (t 0.0625)
Compare it with the standard equation x a cos( 0 t ) . We can write
(i) Amplitude a = 0.01 m
166
2 2
(ii) Period T = 0.5 s
0 4
(iii) Maximum speed = 0 a 4 s 1 0.01m 0.13 ms 1
(iv) Maximum acceleration = 02 a (4 ) 2 s 2 0.01m 1.6m s 2
(v) Displacement at t = 0 is x0 0.01cos 4 0.0625 m
1
= 0.01 m 7.1 10 3 m
2
SAQ 5
The graphs for the variation of U, K.E and E with time are shown in Fig. 1.6.
Since U (t / 0 ) (1 / 2)ka2 cos 2 [ 0 (t / 0 ) ]
= (1 / 2)ka2 cos 2 ( 0 t )
= (1 / 2)ka2 cos 2 ( 0 t )
= U(t)
This means that the period of oscillation of potential energy is / 0 , i.e., one half of
that of vibration.
Similarly, you can show that the period of oscillation of kinetic energy is / 0 . This
means that in each cycle, fixed amount of energy is transferred from the mass to the
spring and back again twice.
SAQ 6
Potential energy (P.E.) = mgh and maximum kinetic energy ( K .E.) max (1 / 2)ka2 . By
equating these expressions we can calculate the amplitude of the oscillation, i.e.,
1/ 2
2mgh
a
k
SAQ 7
U mg(l l cos ) ; K .E. (1 / 2)m(l) 2
E K .E. U mgl(1 cos ) (1 / 2)m(l) 2
[2 / ml 2 {E mgl(1 cos )}]1 / 2
SAQ 8
T 2 I / k t ; k t is the torque producing unit angular displacement.
For a sphere, I (2 / 5)mR 2
2mR 2
T 2
5k t
167
2 3 kg (0.01) 2 m 2
= 2
5 0.04 Nm rad 1
= 0.34 s
SAQ 9
1
0 (2 20 10 3 mH 10 6 F ) 1 = 1125 Hz
2 LC
E (1 / 2)CV 2 (1 / 2) 10 6 F 10 2 V 2 5 10 5 J
SAQ 10
vs A
0
2 Vl
Since V = 2 litres = 2000 c.c. = 2000(10 2 ) 3 m 3 2 10 3 m 3 , we get
350ms 1 (.01) 2 m 2
0
2 2 10 3 m 10 1 m
= 69.8 Hz
An audible note of about this frequency can be heard when an empty flask of this size is
suddenly uncorked.
SAQ 11
The force constant of a diatomic molecule is given by
2e 2 2(1.6 10 19 C ) 2 9 10 9 Nm 2 C 2
k
4 0 r 3 (1.3 10 10 ) 3 m 3
= 209.7 Nm 1
Since
mH mCl 35 1.67 10 27
kg
mH mCl 36
1 1 209.7 Nm 1 36
0 k/
2 2 35 1.67 10 27 kg
= 5.7 1013 Hz
Terminal Questions
1a. In this arrangement, both springs will be extended by the same length x . The restoring
force
F k1 x k 2 x
168
or
d 2x
m (k1 k 2 ) x 0
dt 2
Hence, the period of the system is given by
m
T 2
k1 k 2
b. In this arrangement, if the mass is displaced up or down by x , the restoring forces are
F1 k1 x and F2 k 2 x
Hence,
F k1 x k 2 x
i.e.,
d 2x
m (k1 k 2 ) x 0
dt 2
m
so that the period T 2 .
k1 k 2
c. Here the two springs are connected in series. When the mass is displaced by x , the same
restoring force will act in the springs, extending them by x1 and x 2 due to their different
force constants. Thus,
F k1 x1 k 2 x2
and
x x1 x2 F / k1 F / k 2
i.e.,
x (1 / k1 1 / k 2 ) F
1
or F x
1 1
k1 k 2
m
Therefore time period T = 2 .
1 1
k1 k 2
169
Dividing Eq. (iii) by Eq. (ii), we get
g' / g (R d ) / R (iv)
If the distance to be measured from the center of the earth, let us put R d x . Then,
Eq. (iv) can be rewritten as
d 2x g
g' 2 x
dt R
where the negative sign shows that acceleration is directed towards the center of the
earth. Thus,
d 2x g
x0
dt 2 R
This equation represents SHM whose period is
T 2 R / g
3. x a cos( 0 t )
dx
v a 0 sin( 0 t )
dt
x2
= a 0 1 2
a
Hence, we can write
v1 x2
1 12
a 0 a
v2 x 22
and 1 2
a 0 a
On squaring these expressions, we get
2
v1 x2
1 12
a 0 a
2
v2 x2
1 22
a 0 a
These equations may be combined to give
v12 v 22
0
x 22 x12
Using this result in Eq. (i), we get on simplification
v12 x 22 v 22 x12
a
v12 v 22
170
Suppose T ' is the time period, when the pendulum is suspended from the center of
oscillation. Then we can write
(k r2 / l ' ) l '
T ' 2 (ii)
g
The distance between the center of oscillation and centre of gravity l ' k r2 / l , i.e.,
k r2 ll ' . Using this in Eq. (ii), we get
(k r2 / l ' ) l '
T ' 2 T
g
That is, the periods about the centre of suspension and centre of oscillation are equal.
This property of a compound pendulum is called mutual interchangeability of the centres
of suspension and oscillation. The mutual interchangeability of centres of suspension and
oscillation of a compound pendulum arises because the periods of oscillation about S and
O are equal. In all, there are four points (S, O´, O, S´) on the line SS´(so that GS = GS´
and GO = GO´) about which the periods of oscillation are the same.
e2 c
5. U (r )
4 0 r r9
dU e 2 1 9c
F
dr 4 0 r 2 r 10
At the equilibrium separation, r r0 , force vanishes, i.e.,
e 2 1 9c
0
4 0 r02 r010
e 2 r08
or c
36 0
d 2U e2 2 90c
Now,
dr 2 r r0
4 0 r 3 r r0 r 11 r r0
e2 1 5 e2
=
2 0 r03 r03 2 0
2e 2
=
0 r03
which is positive. Hence, r r0 is the separation at stable equilibrium, and
Force constant
2e 2
k
0 r03
171
SUPERPOSITION OF HARMONIC OSCILLATIONS
Structure
2.1 Introduction Objectives
2.2 Principle of Superposition
2.3 Superposition of Two Harmonic Oscillations of the Same Frequency along the Same
Line
2.4 Superposition of Two Collinear Harmonic Oscillations of Different Frequencies
2.5 Superposition of Many Harmonic Oscillations of the Same Frequency
Method of Vector Addition Method of Complex Numbers '
2.6 Oscillations in Two Dimensions
Superposition of Two Mutually Perpendicular Harmonic Oscillations of the Same
Frequency Superposition of Two Rectangular Harmonic Oscillations with Nearly Equal
Frequencies: Lissajous Figures
2.7 Summary
2.8 Terminal Questions
2.9 Solutions
2.1 INTRODUCTION
In Unit 1, we studied simple harmonic motion and considered a number of examples from
different areas of physics. We found that in each case the motion is governed by a homogeneous
second order differential equation. The solution of this equation gives us information regarding
displacement of the body as a function of time. In many situations, one has to deal with a
combination of two or more simple harmonic oscillations. Do you know that our eardrums vibrate
under a complex combination of harmonic vibrations? The resultant effect is given by the
principle of superposition. You must have observed that oscillations of a swing gradually die out,
when left to itself. This is due to factors like friction and air resistance. The system loses energy
and its motion is said to be damped. We will discuss damped harmonic oscillations in the next
unit.
In this unit we first discuss the principle of superposition. Then you will learn to apply this
principle to situations where two (or more) harmonic oscillations are superposed, either along the
same line or in perpendicular directions.
Objectives
After studying this unit you should be able to
state the principle of superposition
apply the principle of superposition to two harmonic oscillations of (a) the same
frequency and (b) different frequencies along the same line
apply the methods of vector addition and complex numbers for superposition of many
simple harmonic oscillations, and
apply the principle of superposition to two mutually perpendicular harmonic oscillations
of different frequencies/phases and describe the formation of Lissajous figures.
172
initial displacement to be the sum of the earlier displacements, viz., a1 a 2 , then according to
the superposition principle,, the displacement x 3 after the same interval of time t will be
x3 x1 x2
You can perform this activity by taking three identical simple pendulums. Release all three bobs
simultaneously such that their initial velocities are zero and initial displacements of the first,
second and the third pendulum are a1 , a 2 and a1 a 2 , respectively. You will find that at any
time the displacement x 3 of the third pendulum will be the algebraic sum of the displacements of
the other two. In general, the initial velocities may be non-zero. Thus, the principle of
superposition can be stated as follows:
You will note that the principle of superposition holds for any number of simple harmonic
oscillations. These may be in the same or mutually perpendicular directions, i.e. in two
dimensions.
d 2x
02 x (2.1)
dt 2
Such an equation has an important property that the sum of its two linearly independent solutions
is itself a solution. We have already used this property in Unit I while writing Eq. (1.4).
d 2 x1
2
02 x1 (2.2)
dt
d 2 x2
2
02 x2 (2.3)
dt
173
SAQ 1
For a simple pendulum we know that the equation of motion is
d 2
0 sin
dt 2
3 5 7
sin ...
3! 5! 7!
will it remain linear in ? If you retain the first two terms and consider the resulting equation for
the two displacements 1 and 2 , will the principle of superposition still hold? If not, why?
You will find that in the case of a simple pendulum you can apply the principle of superposition
only for small oscillations, i.e. when sin . Here we shall study only those oscillations for
which the displacement satisfies linear homogeneous differential equations.
Let us superpose two collinear (along the same line) harmonic oscillations of amplitudes a1 and
a 2 having frequency 0 and a phase difference of . The displacements of these oscillations
are given by
x1 a1 cos 0 t (2.6)
and x2 a2 cos( 0 t )
= a2 cos 0 t (2.7)
According to the principle of superposition, the resultant displacement is given by
x(t ) = x1 (t ) x2 (t )
= a1 cos 0 t a2 cos 0 t
= (a1 a2 ) cos 0 t (2.8)
174
This represents a simple harmonic motion of amplitude (a1 a2 ) . In particular, if two
amplitudes are equal, i.e. a1 a 2 , the resultant displacement will be zero at all times. The
displacement-time graph depicting this situation is shown in Fig. 2.1.
Fig. 2.1 Superposition of two collinear harmonic oscillations of equal amplitude but out of
phase by
SAQ 2
Two harmonic oscillations of amplitudes a1 and a 2 have the same frequency 0 and are in
phase. Show that their superposition gives a harmonic oscillation of amplitude a1 a 2 .
We will now discuss the general case of superposition of two harmonic oscillations. Let one of
these be characterised by amplitude a1 and inital phase 1 and the other with amplitude a 2 and
phase 2 . Both oscillations have frequency 0 and are collinear, i.e. they are along the same
line. Then, we can write
x1 a1 cos( 0 t 1 ) (2.9)
and x2 a2 cos( 0 t 2 ) (2.10)
According to the principle of superposition, the resultant displacement is the sum of x1 and x 2
and we have
Using the expression for the cosine of the sum of two angles, this can be written as
175
where a and have to be determined. Then, we can rewrite Eq. (2.11) in the form
It has the form of the cosine of the sum of two angles and can be expressed as
This equation has the same form as either of our original equations for separate harmonic
oscillations. Hence, we have the important result that the sum of two collinear harmonic
oscillations of the same frequency is also a harmonic oscillation of the same frequency and along
the same line. But it has a new amplitude and a new phase constant . The amplitude can easily
be calculated by squaring Eqs. (2.12) and (2.13) and adding the resultant expressions. On
simplification we have
a1 sin 1 a 2 sin 2
tan 1 (2.16)
a1 cos 1 a 2 cos 2
SAQ 3
Two harmonic oscillations of frequency 0 have initial phases 1 and 2 and amplitudes a1
and a 2 . Their resultant has the phase
(a) 1 2 2n
and (b) 1 2 (2n 1)
where n is an integer. Using Eq. (2.15), show that the amplitudes of the resultant oscillations are
equal to (a1 a2 ) and (a1 a2 ) , respectively.
SAQ 4
Two harmonic oscillations of frequency 0 having an amplitude 1 cm and initial phases zero and
/ 2 , respectively, are superposed. Calculate the amplitude and the phase of the resultant
vibration.
176
having the same amplitude a but slightly different frequencies
1 and 2 such that 1 2 :
x1 a cos(1t 1 )
x2 a cos( 0 t 2 )
(1 2 )t (1 2 )
x1 (t ) a cos 1t
and x2 (t ) a cos 2 t
(2.17)
2 1 2
x(t ) 2a cos 1 t cos t
2 2
(2.19)
1 2
av
2
(2.20a)
177
and a modulated angular frequency
mod (1 2 ) / 2
(2.20b)
mod 1 2
varies with a frequency
2 4
time.
178
The modulated amplitude and phase constant are
respectively given by
(2.23)
(a a ) sin mod t
and mod 1 2
(a1 a 2 ) cos mod t
(2.24)
or (1 2 )t 2n n = 0, 1, 2, …
1 2 n
or t 0 , , , …,
( 1 2 ) ( 1 2 ) ( 1 2 )
cos 2 mod t 1
1 3 5
t , , ,…
2( 1 2 ) 2( 1 2 ) 2( 1 2 )
179
2.5 SUPERPOSITION OF MANY HARMONIC OSCILLATIONS OF THE SAME
FREQUENCY
In the preceding sections, we considered superposition of two collinear harmonic oscillations. How
will you calculate the resultant of a number of harmonic oscillations of the same frequency? You may
suggest that an obvious way is to extend the procedure outlined in Sec. 2.3. But we find that the
mathematical analysis, though simple, becomes unwieldy. A convenient way out in such a case is to use
either the method of vector analysis or complex numbers. We will now discuss these in turn.
Fig. 2.3 Uniform circular motion and its connection with SHM
length of arc Vt
radius of the circle R
Thus,
dx d
R sin
dt dt
180
= 0 R sin
d V
since 0
dt R
dy
0 R cos
dt
d 2x
02 x
dt 2
d2y
and 2
02 y
dt
181
Fig. 2.4a Projections of a rotating vector along the diameter of
a circle
If you project OP' on the y -axis, you will find that the point corresponding to the foot of the normal
satisfies the equation
y a0 sin( 0 t ) (2.28b)
This means that a rotating vector can, in general, be resolved into two orthogonal components, and we
can write
r xx y y (2.29)
where x x and y y are unit vectors along the x - and the y -axes, respectively.
Vector Addition
Let us now consider the superposition of n harmonic oscillations, all having the same amplitude
a 0 and angular frequency 0 . The initial phases of successive oscillations differ by 0 . Let the
first of these oscillations be described by the equation
x1 (t ) a0 cos 0 t
x2 (t ) a0 cos( 0 t 0 )
x3 (t ) a0 cos( 0 t 2 0 )
.
xn (t ) a0 cos[ 0 t (n 1) 0 ] (2.30)
From the principle of superposition, the resultant oscillation is given by
Let us denote the harmonic oscillations given in Eq. (2.31) as projections of rotating vectors
OP1' , OP2' , OP3' , … (Fig. 2.4b).
182
Fig. 2.4b Projections of rotating vectors OP1' , OP2' , … on the diameter of a circle
To find the resultant of these vectors, we translate them parallel to themselves so that the head of
the first coincides with the tail of the second and so on. You will observe that
(i) combining vectors form successive sides of an incomplete n -sided polygon (Fig.
2.4c)
(ii) OP1' || OP1' , P1' P2' || OP2' , P2' P3' || OP3' and so on.
Let us now project each of these vectors along the x -axis. Then, we get
x1 = Proj (OP1' ) x = a0 cos 0 t
x 2 = Proj (P1' P2' ) x = a0 cos( 0 t 0 )
183
x3 = Proj (P2' P3' ) x = a0 cos( 0 t 2 0 )
.
.
x n = Proj (Pn' 1Pn' ) x = a0 cos( 0 t (n 1) 0 )
The law of vector addition implies that the resultant of OP1' , P1' P2' , P2' P3' , … is given by the
vector OPn' , i.e.,
Let the length of OPn' be a and its phase with respect to the first vector be . Then the
projection of OPn' along the x -axis is given by
x(t ) Proj (OPn' ) x a cos( 0 t )
Hence, the sum in Eq. (2.31) reduces to calculating a and characterising the resultant vector
OPn' . To this end, we recall that any regular polygon will lie on a circle of radius r , as shown in
Fig. 2.4b. The angle subtended at the centre C of the circle by individual vectors will be equal to
0 . Hence, the total angle subtended at C by the resultant vector OPn' will be n 0 .
A circle is an infinite-sided regular polygon
Using the trigonometric relation cos 2 1 2 sin 2 and simplifying the resultant expression,
we get
a 2r sin(n 0 / 2) (2.33a)
184
In the isosceles COP1' , OCP1' = 0 and OP1' C / 2 .
Therefore, 2COPn' n 0
0
Or COPn' n (2.36b)
2 2
That is, the initial phase of the resultant oscillation is equal to half the phase difference between
the n th and the first oscillations. Hence,
sin n 0
x(t ) a0 2 cos t (n 1) 0 (2.38)
0 0 2
sin
2
We shall obtain the same result in the next subsection using the method of complex numbers. For
the moment, let us examine the behaviour of the amplitude of the resultant oscillation defined by
Eq. (2.38):
sin n 0
a a0 2
sin 0
2
You will notice that the value of a depends on the value of 0 . When n is very large, 0
becomes very small. Then, using Eq. (2.37), we can write
0 n 0
(n 1)
2 2
0 0
so that sin
2 2 n
Hence, for large n , we have
185
sin sin
a a0 na0
( / n)
sin
= A
A sin
a A na0
2A
2
2
A sin
Fig. 2.5 Plot of as a function of
That is, the polygon becomes an arc of the circle with centre at O and length na0 A with a as
sin
chord. The plot of A for different values of is shown in Fig. 2.5.
The pattern is symmetric about = 0 and is zero for sin n ( n = 1, 2,...) except at = 0.
When = 0, 0 = 0 and the resultant of n oscillations (vectors) is a straight line of length A. As
increases, A becomes the arc of the circle until at = / 2 the last and first contributions are
out of phase and the arc A becomes a semi-circle whose diameter is the resultant a. A further
increase in curls the length A into the circumference of a circle ( ) with a zero resultant
and so on.
SAQ5
Three collinear harmonic oscillations, represented by x1 4 cos 20 t ,
x2 4 cos(20 t / 3) , x3 4 cos(20 t 2 / 3) are superposed. Determine the amplitude
and phase of the resultant vibration.
186
sin Im[exp(i )] =
Let us now see how this technique of complex numbers is used to obtain the resultant of n
harmonic oscillations given by Eq. (2.30). In the complex exponential notation, we can write
Z1 a0 exp(i 0 t )
Z 2 a0 exp[i( 0 t 0 )
Z 3 a0 exp[i( 0 2 0 )] (2.40)
This series is in geometric progression with common ratio e i0 . Its sum is given by
1 e in0
Z a0 exp(i 0 t )
1 e i 0
e in0 / 2 e in0 / 2 e in0 / 2
= a0 exp( i 0 t ) i / 2 i / 2
e 0 e 0 e i0 / 2
Using the relation
e i e i
sin
2i
we get
sin n 0
in 0 i 0
Z a0 exp(i 0 t )
2
exp exp
0 2 2
sin
2
sin(n 0 / 2) i(n 1) 0
= a0 exp(i 0 t ) exp (2.42)
sin( 0 / 2) 2
Since Z a exp[i( 0 t )] , we find that the amplitude and phase of the resultant vibration are
the same as given by Eqs. (2.34) and (2.37), respectively.
The cosine form of the resultant oscillation is obtained by taking the real part of Eq. (2.42).
187
velocity is maximum and vice versa. Remember that since the time period of the pendulum
depends only on the acceleration due to gravity and the length of the cord, the frequency of the
superposed SHM's will be the same. The result is a curved path, in general, an ellipse.
We now apply the principle of superposition to the case where two harmonic oscillations are
perpendicular to each other.
x1 a1 cos 0 t (2.43)
and x2 a2 cos( 0 t ) (2.44)
Here we have taken the initial phase of the vibrations along the x and the y-axes to be zero and
respectively. That is, is the phase difference between the two vibrations.
We shall first find out the resultant oscillation for a few particular values of phase difference .
Case 1. = 0 or
For = 0
x a1 cos 0 t
and
y a2 cos 0t
Therefore
y / x a2 / a1
or
y (a2 / a1 ) x (2.45)
Similarly, for ,
x a1 cos 0 t
y a2 cos 0 t
So that
y (a2 / a1 ) x (2.46)
Eqs. (2.45) and (2.46) describe straight lines passing through the origin. This means that the
resultant motion of the particle is along a straight line. However, for = 0 the motion is along
one diagonal (PR in Fig. 2.6a) but when the motion is along the other diagonal (QS in Fig.
2.6b).
188
Case II. = / 2
In this case the two vibrations are given by
x a1 cos 0 t
y a2 cos( 0 t / 2) = a2 sin 0 t
Fig. 2.6 Superposition of two mutually perpendicular harmonic oscillations having the same
frequency different phases
This is the equation of an ellipse. Thus the resultant motion of the particle is along an ellipse
whose principal axes lie along the x- and the y-axes. The semi-major and semi-minor axes of the
ellipse are a1 and a 2 . Note that as time increases x decreases from its maximum positive value
but y becomes more and more negative. Thus, the ellipse is described in the clockwise direction
as shown in Fig. 2.6c. If you analyse the case when 3 / 2 or / 2 , you will obtain the
same ellipse, but the motion will be in the anticlockwise direction (Fig. 2.6d).
x2 y2 a2
This equation represents a circle of radius a . This means that the ellipse reduces to a circle.
General case
We will now consider the general case for any arbitrary value of . Let the two SHM’s given by
Eqs. (2.43) and (2.44) be superposed. To find the resultant oscillation, we write Eq. (2.44) as
189
y
cos( 0 t ) cos 0 t cos sin 0 t sin (2.48)
a2
x2 y2 xy
2
2 2 cos sin 2 (2.49)
a1 a 2 a1a 2
as the equation of the resultant path. This describes an ellipse whose axes are inclined to the
coordinate axes.
For some typical values of lying between 0 and , the resultant paths traced out by the
resultant oscillation when two mutually perpendicular SHM’s of equal frequency are superposed
are shown in Fig. 2.7. These can be most easily demonstrated on a cathode ray oscilloscope.
Fig. 2.7 Superposition of two mutually perpendicular harmonic oscillations of the same
frequency and having values of lying between 0 and 2
190
We may thus conclude that an elliptical motion is a combination of two mutually perpendicular
linear harmonic oscillations of unequal amplitudes and having a phase difference . A circular
motion is a combination of two harmonic oscillations of equal amplitudes
SAQ6
In a cathode ray oscilloscope, the deflection of electrons by two mutually perpendicular electric
fields is given by
x 4 cos 2 t
y 4 cos(2 t / 6)
We now know that when two orthonormal vibrations have exactly the same frequency, the shape
of the curve traced out by the resultant oscillation depends on the phase difference between
component vibrations. For a few values of the phase difference in the range 0 to 2 radians,
these curves are shown in Fig. 2.6. When the two individual rectangular vibrations are of slightly
different frequencies, the resulting motion is more complex. This is because the relative phase
[ 2 t 0 1t ( 2 1 )t 0 ] of the two vibrations gradually changes with time. This
makes the shape of the figure to undergo a slow change. If the amplitudes of vibrations are a1
and a 2 , respectively, then the resulting figure always lies in a rectangle of sides 2a1 and 2a 2 .
The patterns which are traced out are called Lissajous figures. When the two vibrations are in the
same phase, i.e. = 0, the Lissajous figure reduces to a straight line and coincides with the
diagonal y = (a2 / a1 ) x of the rectangle. As changes from 0 to / 2 , the Lissajous figure is an
ellipse and passes through oblique positions in the rectangle. When increases from / 2 to ,
the ellipse closes into a straight line which coincides with the (other) diagonal y (a2 / a1 ) x of
the rectangle. Further, as changes from to 2 , the series of changes mentioned above take
place in the reverse order. In general, the shape of curve depends on the amplitudes, frequencies
and the phase difference. All these changes are shown'in Fig 2.7. The phase changes by 2 in
the time interval 2 /( 2 1 ) . Therefore, the period of the complete cycle of changes is
2 1
2 /( 2 1 ) and its frequency is 1 2 , i.e., equal to the difference of the
2
frequencies of the individual vibrations.
Lissajous figures can be illustrated easily by means of a cathode ray oscilloscope (CRO).
Different alternating sinusoidal voltages are applied at XX and YY deflection plates of the CRO.
The electron beam traces the resultant effect on the fluorescent screen. When the applied voltages
have the same frequency, we can obtain various curves of Fig. 2.7 by adjusting the phases and
amplitudes.
If the frequencies of individual perpendicular vibrations are in the ratio 2:1, the Lissajous figures
are relatively complex. It has the shape of parabola for = 0 or and for = / 2 its shape is
that of figure '8'. To clarify this let us study the following example:
191
Two rectangular harmonic vibrations having frequencies in the ratio 2:1 are represented as
follows:
x a1 cos(2 0 t ) (2.50)
and y a2 cos 0 t (2.51)
(ii) When
2
x a1 sin 2 0 t
x
or 2 sin 0 t cos 0 t
a1
and y a2 cos 0 t
Since we can write
y
cos 0 t
a2
y2
and sin 0 t 1 2
a2
x 2y y2
1 2
a1 a 2 a2
192
4y2 y2 x2
2 1 2 0 (2.8b)
a 22 a2 a1
(iii) When
x a1 cos 2 0 t
x
or 2 cos 2 0 1
a1
and y a2 cos 0 t
0 /2
Fig 2.8 Superposition of two harmonic oscillations having frequencies in the ratio 2:1 and
phase difference = (i) 0, (ii) / 2 and , respectively
2y2 x
2
1
a2 a1
a 22
or y
2
( x a1 ) (2.54)
2a1
2.7 SUMMARY
The principle of superposition states that when we superpose the initial conditions corresponding
to velocities and amplitudes, the resultant displacement of two (or more) harmonic displacements
will be simply the algebraic sum of the individual displacements at all subsequent times.
x(t ) x1 (t ) x2 (t )
x1 a1 cos( 0 t 1 )
and x2 a2 cos( 0 t 2 )
193
are superposed, the resultant is given by
x(t ) a cos( 0 t )
where a [a12 a22 2a1a2 cos(1 2 )]1 / 2
a1 sin 1 a 2 sin 2
and tan 1
a1 cos 1 a 2 cos 2
When two collinear harmonic oscillations of different frequencies are superposed, the modulated
oscillation is represented by
x(t ) amod (t ) cos avt
1 2
with mod
2
1 2
and av
2
Superposition of n harmonic collinear oscillations of the same amplitude (a0 ) and frequency ( 0 )
but having a constant phase difference ( 0 ) between successive oscillations yields a harmonic
oscillation. It is given by
x(t ) a cos( 0 t )
n
a 0 sin 0
where a
2
sin 0
2
0
and (n 1)
2
When two mutually perpendicular harmonic oscillations are superpose, the resultant form traces
out different curves. If the oscillations have equal frequencies, the shape of the curve depends on
the phase difference. In general, the curve is elliptical but for certain phases, it closes into a
straight line. When the frequencies are nearly equal, the curves are termed Lissajous figures.
194
Show that for a new displacement x3 x1 x2 , the initial conditions on the bob are the
superposition of the initial conditions of x1 and x 2 .
2.9 SOLUTIONS
SAQ
1. On using the given expansion, we get
d 2 3 5
2
0 ... 0
dt 3! 5!
Since this equation contains terms of power higher than , it is not linear.
Even if we retain the first two terms in the expansion, the resulting equation will not be
linear and hence the principle of superposition will not hold.
2. x1 a1 cos 0 t
x2 a2 cos 0 t
According to the principle of superposition
x x1 x2 (a1 a2 ) cos 0 t
Since the cosine function varies between + 1 and 1 , the amplitude of the resultant
oscillation is | a1 a 2 | .
3. The resultant of two harmonic oscillations having amplitudes a1 and a 2 and initial
phases 1 and 2 is given by
195
a 2 a12 a22 2a1a2 cos(1 2 ) (i)
(a) When 1 2 = 2n , cos(1 2 ) 1 and Eq. (i) reduces to
5. Here, n = 3, a 0 = 4 and 0 / 3 rad. From Eq. (2.34) we note that the amplitude
of the resultant oscillation is given by
n
sin 0
a a0 2
sin 0
2
3
sin
= a0
23
sin
23
sin
= a0
2
sin
6
1
Since sin 1 and sin , we get
2 6 2
196
a 2a0 = 8 units
0
(n 1)
2
= 2
6
= rad.
3
Terminal Questions
d 2x
1. 4x 0 (i)
dt 2
d 2x
0 x 0 ,
2
Comparing it with the standard equation for SHM 2
dt
We find that the solution of this equation is
x a cos(2t ) (ii)
Differentiating Eq. (ii) with respect to t , we get
dx
2a sin(2t ) (iii)
dt
dx
(1) Since at t = 0, x 3 cm and 0 , from Eqs. (ii) and (iii), we obtain
dt
3 cm = a cos
and 0 2a sin
The latter of these two relations implies that = 0. Using this in the former, we
get
a = 3 cm
Therefore, the complete solution is
x(t ) 3 cos 2t cm (iv)
dx
(2) Again if at t 0, x = 2 cm and 4 cm s 1 , we find
dt
197
2 cm = a cos
and 4 cm s 1 = 2a sin
or 2 cm s 1 = a sin
2
Hence, tan
5
2
sin
29
5
cos
29
and a 29 cm
198
2. x1 3 cos(20 t / 6 / 2) cm
and x2 4 cos(20 t / 3 / 2) cm
or x1 3 cos(20 t / 3) cm
and x2 4 cos(20 t / 6) cm
= (9 16 12 3 )1 / 2 cm
= 6.77 cm
3 sin / 3 4 sin / 6 1 3 3 4
and tan 1 = tan
3 cos / 3 4 cos / 6 3 4 3
= 0.24
3. z a1 exp(i1t ) a2 exp(i 2 t )
a 2 ( zz*) (a1e i1t a2 e i2t ) (a1e i1t a2 e i2t )
= a12 a22 a1a2 exp[ i(1 2 )t ] a1a2 exp[ i(1 2 )t ]
4. A B 0.05 Hz
Now on loading the prong of the tuning fork A with wax, the frequency of A will decrease. However, now the cycle of chang es
of figures is completed in 10 s and hence the frequency difference increases to 0.1. This means that the frequencies of A before
and after loading are, respectively, (300-0.05) Hz = 299.95 Hz and (300 – 0.1) Hz = 299.9 Hz.
199
3.1 Introduction
Objectives
3.2 Differential Equation of a Damped Oscillator
3.3 Solutions of the Differential Equation
Heavy Damping
Critical Damping
Weak or Light Damping
3.4 Average Energy of a Weakly Damped Oscillator Average Power Dissipated Over One
Cycle
3.5 Methods of Describing Damping
Logarithmic Decrement
Relaxation Time
The Quality Factor
3.6 Examples of Weakly Damped Systems
An LCR Circuit
A Suspension Type Galvanometer
3.7 Summary
3.8 Terminal Questions
3.9 Solutions
3.1 INTRODUCTION
In Unit 1 you learnt that SHM is a universal phenomenon. Now you also know that in the ideal
case the total energy of a harmonic oscillator remains constant in time and the displacement
follows a sine curve. This implies that once such a system is set in motion it will continue to
oscillate forever. Such oscillations are said to be free or undamped. Do you know of any physical
system in the real world which experiences no damping? Probably there is none. You must have
observed that oscillations of a swing, a simple or torsional pendulum and a spring-mass system
when left to themselves, die down gradually. Similarly, the amplitude of oscillation of charge in
an LCR circuit or of the coil in a suspended type galvanometer becomes smaller and smaller. This
implies that every oscillating system loses some energy as time elapses. The question now arises:
Where does this energy go? To answer this, we note that when a body oscillates in a medium it
experiences resistance to its motion. This means that damping force comes into play. Damping
force can arise within the body itself, as well as due to the surrounding medium (air or liquid).
The work done by the oscillating system against the damping forces leads to dissipation of energy
of the system. That is, the energy of an oscillating body is used up in overcoming damping. But in
some engineering systems we knowingly introduce damping. A familiar example is that of brakes
– we increase friction to reduce the speed of a vehicle in a short time. In general, damping causes
wasteful loss of energy. Therefore, we invariably try to minimise it.
Many a time it is desirable to maintain the oscillations of a system. For this we have to feed
energy from an outside agency to make up for the energy losses due to damping. Such
oscillations are called forced oscillations. You will learn various aspects of such oscillations in
the next unit.
In this unit you will learn to establish and solve the equation of motion of a damped harmonic
oscillator. Damping may be quantified in terms of logarithmic decrement, relaxation time and
quality factor. You will also learn to compute expressions for the logarithmic decrement, power
dissipated in one cycle and the quality factor.
Objectives
After going through this unit you will be able to:
200
establish the differential equation for a damped harmonic oscillator and solve it
analyse the effect of damping on amplitude, energy and period of oscillation
highlight differences between weakly damped, critically damped and over-damped
systems
derive expressions foe power dissipated in one oscillation
compute relaxation time and quality factor of a damped oscillator, and
draw analogies between different physical- systems.
Fd v (3.1)
The force experienced by a body falling freely in a viscous medium is given by Fd 6rv . This is
known as Stokes' law. Here is the coefficient of viscosity of the medium and r is radius of body –
assumed to be spherical, and v is its velocity.
The negative sign signifies that the damping force opposes motion. The constant of
proportionality is called the damping coefficient. Numerically, it is equal to force per unit
N kgms2
velocity and is measured in 1
1
kg s 1 .
ms ms
We will now establish the differential equation which describes the oscillatory motion of a
damped harmonic oscillator. Let us take the x -axis to be along the length of the spring. We
define the origin of the axis (x = 0) as the equilibrium position of the mass. Imagine that the mass
201
(in the spring-mass system) is pulled longitudinally and then released. It gets displaced from its
equilibrium position. At any instant, the forces acting on the spring-mass system are:
(i) a restoring force: kx where k is the spring factor, and
dx
(ii) a damping force: v . where v is the instantaneous velocity of the oscillator.
dt
This means that for a damped harmonic oscillator, the equation of motion must include the
restoring force as well as the damping force. Hence, in this case Eq. (1.2) is modified to .
d 2x dx
m 2
kx (3.2)
dt dt
After rearranging terms and dividing throughout by m, the equation of motion of a damped
oscillator takes the form
d 2x dx
2
2b 02 x 0 (3.3)
dt dt
where 02 k / m . (You will note that a factor of 2 has been introduced in the damping term as it
helps us to obtain a neat expression for the solution of this equation.) The constant b has the
dimensions of
force MLT 2
1
T 1
velocity mass LT M
You will note that like Eq. (1.3), Eq. (3.3) is a linear second order homogeneous differential
equation with constant coefficients. If there were no damping, the second term in Eq. (3.3) will be
zero and the general solution of the resulting equation will be given by Eq. (1.5), i.e.
x A cos( 0 t ) . On the other hand, if there is damping and no restoring force, the third term
in Eq.(3.3) will be zero. Then the general solution of the resulting equation is given by
x(t ) Ce 2bt D , where C and D are constants. (You can show this by substituting the assumed
solution in Eqs. (3.3).) This means that the displacement will decrease exponentially in the
absence of any restoring force. Thus we expect that the general solution of Eq. (3.3) will
represent an oscillatory motion whose amplitude decreases with time.
202
Differentiating Eq.(3.4) twice with respect to time, we get
dx
a exp( t )
dt
d 2x
and 2
a 2 exp( t )
dt
a =0
which is trivial, or
2 2b 02 0 (3.6)
b b 2 4ac
x
The roots of the equation ax bx c 0 are given by
2
2a
1 b (b 2 02 )1 / 2 (3.7a)
and 2 b (b 2 02 )1 / 2 (3.7b)
These roots determine the motion of the oscillator. Obviously has dimensions of inverse time.
Did you not expect it from the form of exp( t ) ?
x1 (t ) a1 exp[ {b (b 2 02 )1 / 2 }t ]
and x2 (t ) a2 exp[ {b (b 2 02 )1 / 2 }t ] (3.8)
Since Eq.(3.3) is linear, the principle .of superposition is applicable. Hence, the general solution
is obtained by the superposition of x1 and x 2 :
203
You will note the quantity (b 2 02 ) can be negative, zero or positive respectively, depending on
whether b is less than, equal, to or greater than 0 respectively. We, therefore, have three
possibilities:
(i) If b 0 , we say that the system is over damped,
(ii) If b 0 , we have a critically damped system,
(iii) If b 0 , we have an under-damped system.
Each of these conditions gives a different solution, which describes a particular behaviour.
b 2 02
The general solution for damped oscillator given by Eq. (3.9) reduces to
This represents non-oscillatory behaviour. Such a motion is called dead-beat. The actual
displacement will, however, be determined by the initial conditions. Let us suppose that to begin
with, the oscillator is at its equilibrium position, i.e x 0 at t 0 . Then we give it a sudden
kick so that it acquires a velocity v 0 , i.e v v 0 at t 0 . Then from Eq. (3.10) we have
a1 a2 0
and b(a1 a2 ) (a1 a2 ) v0
v0
a1 a 2
2
On substituting these results in Eq. (3.10), we can write the solution in compact form:
v0
x(t ) exp( bt )[exp( t ) exp( t )]
2
204
v0
= exp( b t ) sinh t (3.11)
where sinh t (1 / 2)[exp( t ) exp( t )] is the hyperbolic sine function. From Eq. (3.11) it
is clear that x(t ) will be determined by the interplay of an increasing hyperbolic function and a
decaying exponential. These are plotted separately in Fig. 3.2(a). Fig. 3.2(b) shows the plot of Eq.
(3.11) for a heavily damped system when it is suddenly disturbed from its equilibrium position.
You will note that initially the displacement increases with time. But soon the exponential term
becomes important and displacement begins to decrease gradually.
y y
y sinh x y ex
0 x 0 x
Fig. 3.2a Plot of sinh x and exp( x)
x(t )
0 t
Fig. 3.2b Plot of Eq. (3.11) for a heavily damped system
Mathematically, we say that a system is critically damped if b is equal to the natural frequency,
0 , of the system. This means that b 2 02 = 0, so that Eq (3.9) reduces to
205
where a a1 a2 .
Let us pause for a minute and recall that the solution of the differential equation for SHM
involves two arbitrary constants which are fixed by giving the initial conditions. But Eq. (3.12)
has only one constant. Does this mean that it is not a complete solution? It is important to
understand how this happens. The reason is simple: the quadratic equation for (Eq. 3.6) has
equal roots. So, the two terms in Eq. (3.9) give the same time dependence and reduce to one term.
It can be easily verified that in this case the general solution of Eq. (3.3) is
where p and q are constants, p has the dimensions of length and q those of velocity. These can be
determined easily from the initial conditions.
Let us assume that the system is disturbed from its mean equilibrium position by a sudden
impulse. (The coil of a suspended type galvanometer receives some electric charge at t = 0.) That
dx
is, at t = 0, x(0) = 0 and v0 . This gives p = 0 and q v0 , so that the complete solution
dt t 0
is
Fig. 3.3 illustrates the displacement-time graph of a critically damped system described
dx d 2x
by Eq.(3.13 b). At maximum displacement, 0 and 0.
dt x xmax dt 2 x xmax
v0 mv
xmax v0 te 1 0.368 0.736 0
b
Fig. 3.3 Displacement-time graph for a critically damped system described by Eq. (3.13b)
206
3.3.3 Weak or light damping
When b 0 , we refer to it as a case of weak damping. This implies that (b 2 02 ) is a
negative quantity, i.e. (b 2 02 )1 / 2 is imaginary. Let us rewrite it as
(b 2 02 )1 / 2 1( 02 b 2 )1 / 2
= i d
where i 1 and
1/ 2
k 2
d ( b )
2
0
2 1/ 2
2
(3.14)
m 4m
is a real positive quantity. You will note that for no damping (b 0) , d reduces to 0 , the
natural frequency of the oscillator.
On combining Eqs. (3.9) and (3.14) we find that the displacement now has the form
To compare the behaviour of a damped oscillator with that of a free oscillator, we should recast
Eq.(3.15) so that the displacement varies sinusoidally. To do this, we write the complex
exponential in terms of sine and cosine functions. This gives
a1 a2 a0 cos
and i(a1 a2 ) a0 sin (3.17)
207
From the second of these results we note that tan 0 is a complex quantity. Does this moan that
is also complex? How can we interpret a complex angle? To know this, we use the identity
sec 2 1 tan 2
a1 a 2
cos
2 a1 a 2
with d as given by Eq. (3.14). You will note that the solution given by Eq. (3.19) describes
sinusoidal motion with frequency d which remains the same throughout the motion. This
property is crucial for the use of oscillators in accurate timepieces. How is the amplitude modified
vis-a-vis an ideal SHM? You will note that the amplitude decreases exponentially with time at a
rate governed by b. So we can say that the motion of a weakly damped system is not simple
harmonic.
The damped oscillatory behaviour described by Eq. (3.19) is plotted in Fig. 3.4 for the particular
case of = 0. Since the cosine function varies between +1 and 1 , we observe that the
displacement-time curve lies between a0 exp( bt ) and a0 exp( bt ) . Thus, we may conclude
that damping results in decrease of amplitude and angular frequency.
x(t )
a0 exp( bt )
t
0
a0 exp( bt )
How does damping influence the period of oscillation? You can discover this effect by noting that
the period of oscillation is given by
208
2 2 2
T
d ( b 2 )1 / 2
2
k 2
0
m 4m 2
If b 0 , d 0 . This means that the period of vibration of a damped oscillator is more than
that of an ideal oscillator. Did you not expect it since damping forces resist motion?
SAQ 1
The amplitude of vibration of a damped spring-mass system decreases from 10 cm to 2.5 cm in
200 s. If this oscillator performs 100 oscillations in this time, compare the periods with and
without damping.
We have discussed solutions of the differential equation for a damped oscillator for heavy, critical
and weak dampings. In the following discussion we shall concentrate only on weakly damped
systems.
E(t ) K .E.(t ) U (t )
2
1 dx 1
= m kx2 (3.20)
2 dt 2
For a weakly damped harmonic oscillator, the instantaneous displacement is given by Eq.(3.19):
x(t ) a0 exp( bt ) cos( d t ) . By differentiating it with respect to time, we get instantaneous
velocity:
dx(t )
v a0 exp( bt )[b cos( d t ) d sin( d t )]2 (3.21)
dt
2
1 dx 1
K .E. m ma02 exp( 2bt )[b cos( d t ) d sin( d t )]2
2 dt 2
1
= ma02 exp( 2bt )[b 2 cos 2 ( d t ) d2 sin 2 ( d t )
2
209
+ b d sin 2( d t )] (3.22a)
1 2 1
U kx m 02 x 2
2 2
since k m 02 .
1
U ma02 02 exp( 2bt ) cos 2 ( d t ) (3.22b)
2
1
E (t ) ma02 exp( 2bt )[(b 2 02 ) cos 2 ( d t ) d2 sin 2 ( d t )
2
+ b d sin 2( d t )] (3.23)
When damping is small, the amplitude of oscillation does not change much over one oscillation.
So we may take the factor exp( 2bt ) as essentially constant. Further, since
sin 2 ( d t ) cos 2 ( d t ) 1 / 2 and sin( d t ) 0 , the energy of a weakly
damped oscillator when averaged over one cycle is given by
1
E ma02 exp( 2bt ) [(b 2 02 ) cos 2 ( d t ) d2 sin 2 ( d t )
2
+ b d sin 2( d t )]
1 b 2 02 d2
= ma02 exp( 2bt )
2 2 2
1
= ma02 02 exp( 2bt ) (3.24a)
2
1
From Unit 1 we recall that ma02 02 is the total energy of an undamped oscillator. Hence, we
2
can write
E E0 exp( 2bt ) (3.24b)
This shows that the average energy of a weakly damped oscillator decreases exponentially with
time. This is illustrated in Fig. 3.5. From Eq.(3.24 b) you will also observe that the rate of decay
of energy depends on the value of b; larger the value of b, faster will be the decay.
210
E
E0
E0 exp( bt )
0 t
Fig. 3.5 Time variation of average energy for a weakly damped system
dE d 2x dx
P(t ) m 2 kx
dt dt dt
On combining this result with Eq. (3.2) we find that power dissipated by a damped oscillator is
given by
2
dx
P(t )
dt
This relation shows that the rate of doing work against the frictional force is directly proportional
dx
to the square of instantaneous velocity. On substituting for from Eq. (3.21), we obtain
dt
1
P a02 02 exp( 2bt )
2
= E
m
= 2b E
211
3.5 METHODS OF CHARACTERISING DAMPED SYSTEMS
We now know that in the viscous damping model, a damped oscillator is characterised by and
0 . We also know that this model applies to vastly different physical systems. Therefore, you
may ask: Are there other ways of characterizing damped oscillations? Experience tells us that in
certain cases it is more convenient to use other parameters to characterise damped motion. In all
cases we can relate these to y and coo. We will now discuss these briefly.
Fig. 3.6 A damped oscillation. The first two amplitudes are a 0 and a1
You will note that these amplitudes lie in the same direction/quadrant. If T is the period of
oscillation, then using Eq. (3.19) for a weakly damped oscillator, we can write
a1 a0 exp( bT )
a0
so that exp(bT ) exp( T / 2m) (3.26)
a1
You will note that in the ratio a0 / a1 , the larger amplitude is in the numerator. That is why this
ratio is called the decrement. It is denoted by the symbol d. You may now ask: Is the decrement
same for any two consecutive amplitudes? The answer is: yes, it is. To show this let us consider
the ratio of the second and the third amplitudes. These are observed for t = Tand t = 2T,
respectively in Eq. (3.19). Then, we can write
212
a1 a exp( bT )
0 exp(bT )
a 2 a0 exp( 2bT )
So, we may conclude that for any two consecutive amplitudes separated by one period, we have
a n 1
exp(bT ) d (3.27)
an
That is, decrement is the same for two successive amplitudes and we can write
a0 a1 a 2 a
... n 1 d (3.28)
a1 a 2 a3 an
The logarithm of the ratio of successive amplitudes of oscillation separated by one period is
called the logarithmic decrement. It is usually denoted by the symbol :
a T
ln n1 (3.29a)
a n 2m
This equation shows that we can measure by knowing two successive amplitudes. But from an
experimental point of view it is more convenient and accurate to compare amplitudes of
oscillations separated by n periods. That is, we measure a0 / a n . To compute this ratio, we first
invert Eq. (3.29 a) to write
a n 1
exp( ) (3.29b)
an
The ratio a0 / a n can now be written as
a0 a0 a1 a 2 a n1
... [exp( )]n
a n a1 a 2 a3 a n
(3.30)
1 a
ln 0 (3.31)
n an
This shows that if we plot ln( a0 / an ) versus n for different values of n, we will obtain a straight
line. The slope of the line gives us .
SAQ2
213
A damped harmonic oscillator has the first amplitude of 20 cm. It reduces to 2 cm after 100
oscillations, each of period 4.6 s. Calculate the logarithmic decrement and damping constant.
Compute the number of oscillations in which the amplitude drops
by 50%.
If we denote the amplitude of oscillation after an interval of time r by a (t + T), we can write
a(t )
exp( b )
a(t )
1
= for b = 1 (3.32)
e
This shows that for b 1 the amplitude drops to 1 / e = 0.368 of its initial value. Using this
result in Eq. (3.25), we get
2E
P
The relaxation time, , is therefore a measure of the rapidity with which motion is damped. (You
will note that the negative sign occurring in Eq. (3.25) has been dropped here.)
d m d d
Q (3.33)
2b 2
214
You will note that Q is only a number and has no dimensions. In general, is small so that Q is
very large. A tuning fork has Q of a thousand or so, whereas a rubber band exhibits a much lower
(~10) Q. This is due to the internal friction generated by the coiling of the long chain of
molecules in a rubber band. An undamped oscillator ( 0) has an infinite quality factor.
For a weakly damped mechanical oscillator, the quality factor can be expressed in terms of the
spring factor and damping constant. For weak damping,
d 0 k / m
Hence, Q km / 2
That is, the quality factor of a weakly damped oscillator is directly proportional to the square root
of k and inversely proportional to .
We can rewrite Eq. (3.33) in a more physically meaningful form using Eq. (3.25):
d 2 E
Q
2b Td P
average energy stored in the systemin one cycle
= 2 (3.34)
average energy lost in one cycle
The quality factor is related to the fractional change in the frequency of an undamped oscillator.
To establish this relation, we note that,
d 02 b 2
d2 b2
or 1
02 02
1
= 1
4Q 2
where we have used Eq. (3.33). This result can be rewritten as
d
1/ 2
1
1
0 4Q 2
1
= 1
8Q 2
where in the binomial expansion we have retained terms up to first order in Q 2 . Hence, the
fractional change in 0 is 1 /(8 Q 2 ) .
215
SAQ4
The quality factor of a tuning fork of frequency 256 Hz is 10 3. Calculate the time in which its
energy becomes 10% of its initial value.
q dI dq
L
C dt dt
d 2q dq q
L 2
R 0 (3.36)
dt dt C
Comparing it with Eq. (3.2) we find that L, R and 1/C are respectively analogous to m, and k.
This means that a resistor in an electric circuit has an exactly analogous effect as that of the
viscous force in a mechanical system.
d 2 q R dq 1
2
q0 (3.37)
dt L dt LC
216
In this form, Eq. (3.37) is analogous to Eq. (3.3) and the two may be compared directly. This
gives
1
02
LC
R
and b (3.38)
2L
We know that b has dimensions of time inverse. This means that R/L has the unit of s 1 , same as
that of 0 . That is why 0 L is measured in ohm.
With these analogies all the results of Section 3.3 apply to Eq. (3.37). For a weakly damped
circuit, the charge on the capacitor plates at time t is
R
q(t ) q0 exp t cos( d t ) (3.39a)
2L
with angular frequency
1 R2
d 2 (3.39b)
LC 4 L
R
Eq. (3.39 a) shows that the charge amplitude q0 exp t will decay at a rate which depends
2L
on the resistance. Thus in an LCR circuit, resistance is the only dissipative element; an increase in
R increases the rate of decay of the charge and decreases the frequency of oscillations.
For R 0 , Eq (3.39 a) reduces to Eq. (1.38) and d 0 . The Q value of a weakly damped
LCR circuit is
d L 1 L
Q 0 (3.40)
2b R R C
This equation shows that for a purely inductive circuit ( R 0 ), the quality factor will be infinite.
SAQ 5
In an -LCR circuit, L= 2mH and C= 5 F. If R =1 , 40 and I00 , calculate the
frequency of oscillation and the quality factor when the discharge is oscillatory.
217
always along the magnetic lines of force. To ensure uniform strength, an iron cylinder is
suspended between the poles of the magnet, as shown in Fig. (3.8). When we pass Charge
through the galvanometer coil, it rotates through some angle . Since the coil is mechanically a
d
torsional pendulum, it experiences a restoring couple k t and a damping couple . Do
dt
you know how damping creeps in, in this case? It has origin in air friction and electromagnetic
induction.
Part of the damping arises from the viscous drag of air. In general, it is small. As the
galvanometer coil rotates in the magnetic field, an induced e.m.f. is produced, which
opposes its motion in accordance with Lenz's law. This so-called electromagnetic
damping controls the motion of the coil when galvanometer is in use.
d 2 d
I 2
k t (3.41)
dt dt
where I is the moment of inertia of the coil about the axis of suspension. Comparing it with Eq.
(3.2) we find that I and k t are analogous to m and k respectively.
218
2b / I
and 02 k t / I (3.42)
d 2 d
we get 2
2b 02 0 (3.43)
dt dt
This equation is of the same form as Eq. (3.3). Hence, all results deduced earlier will apply to the
motion of the coil described by Eq. (3.43).
where 0 exp( bt ) is the amplitude of oscillation. Eq. (3.44) describes oscillatory motion with
the period of oscillation T given by
2 2 2
T (3.45)
0 ( b )
2
0
2 1/ 2
kt 2
1/ 2
I 4I 2
This explains why a weakly damped suspension type galvanometer is called a ballistic
galvanometer. You will note that for damping to be small, we must decrease and increase I .
The question now arises: How can we reduce ? As mentioned earlier, air damping is usually
small. Nevertheless, it will always be present. To reduce electromagnetic damping, we must
minimise the induced emf. To ensure this, we wind the coil over a non-conducting bamboo or
ivory frame. If the frame is metallic, it is cut at one place, so that no current can flow through it.
The quality factor of a ballistic galvanometer is
d I kt 2
Q (3.46a)
2b I 4I 2
kt 2
If 2 , this expression reduces to
I 4I
kt I
Q 2
This relation shows that a lightly damped suspension type galvanometer will have high quality
factor.
SAQ6
The period of vibration of a galvanometer coil is 4 s. The amplitude of its vibration decreases to
one-tenth of its original value in 46 s. Calculate the damping constant and the quality factor.
219
3.7 SUMMARY
1. The differential equation of a damped harmonic oscillator is
d 2x dx
2
2b 02 x 0
dt dt
where 2b / m and 02 k / m . The solution of the equation for heavy damping is
x(t ) exp(bt )[a1 exp( t ) a2 exp( t )]
where b 2 02
2. The amplitude and average energy of a weakly damped oscillator decrease exponentially
with time:
a a0 e bt
and E E0 exp( 2bt )
where a 0 is the initial amplitude and E 0 is the total initial energy.
0
2
m 4m
5. The rate of loss of energy or power dissipated by a weakly damped harmonic oscillator
over one cycle is
P 2 E /
d 2 q R dq 1
2
q0
dt L dt LC
220
The effect of L, R and 1/C in an LCR circuit is respectively analogous to those of m ,
and k in a mechanical oscillator. In a weakly damped circuit, the charge oscillates
harmonically:
R
q(t ) q0 exp t cos( d t )
2L
and the frequency of oscillation is given by
1 1 R2
d 2
2 LC 4 L
For low R circuit
1
Q 0 L / R L/C
R
2. The quality factor of a sonometer wire is 4,000. The wire vibrates at a frequency of
300 Hz. Find the time in which the amplitude decreases to half of its original value.
3. A box of mass 0.2 kg is attached to one end of a spring whose other end is fixed to a
rigid support. When a mass of 0.8 kg is placed inside the box, the system performs 4
oscillations per second and the amplitude falls from 2 cm to 1 cm in 30 s. Calculate
(i) the force constant, (ii) the relaxation time and (iii) the Q-factor.
5. The quality factor of a tuning fork of frequency 512 Hz is 6 10 4 . Calculate the time
in which its energy is reduced to e 1 of its energy in the absence of damping. How
many oscillations will the tuning fork make in this time?
3.9 SOLUTIONS
SAQ’s
221
200s
1. Td =2 s
100
2
Now, Td 2s
( b 2 )1 / 2
2
0
so that 02 2 b 2
2 2
Hence, T0
0 ( b 2 )1 / 2
2
2
T0 2 s Td
[ (0.015942) 2 ]1 / 2
2
2. We know that
1 a
ln 0
n an
1
= ln 10
100
2.3
= log 10 10 2.3 10 2
100
Since b
T
we get
2.3 10 2
b
4.6s
= 5.0 10 3 s 1
Further, to calculate n for which the amplitude drops by 50%, we invert (i) to write
1 a
n ln 0
an
222
ln 2 2.3 log 10 2
= 2
2.3 10 2.3 10 2
= 30
3. Since
1 a
bT ln 0
n an
we can write
1
b ln( a0 / a n )
nT
1
= ln 4
200 s
2.3 0.6010 1
= s
200
Hence,
1 200
145 s
b 2.3 0.6010 s 1
4. Q = 10 3 and = 256 Hz
We know that for weak damping, Q is given by
0
Q
2
On inverting this relation, we get
Q 10 3
= 1.24 s
256 s 1
Since E E0 exp( 2bt ) E0 exp( 2t / ) , we get for E / E0 1 / 10
1 2t
exp
10
Hence,
t ln 10
2
1.24
= 2.3
2
= 1.4 s
5. L 2 10 3 H and C = 5 10 6 F
1 1
3 4
108 s 2
LC 2 10 H 5 10 F
Case I: R = 1
223
1 2 4
2
R2
= 6 10
4 L2 4 (2 10 3 ) 2 H 2 H2
Thus,
1 R2
2 so that the discharge is oscillatory. The frequency of oscillation
LC 4L
1 1 R2
2 LC 4 L2
= 1.6 kHz
Case II: R = 40
In this case,
R2 40 40 2 8
2
10
4 L2 4 (2 10 3 ) 2 H 2 H2
1 R2
Hence, and this is the case of critical damping.
LC 4 L2
You will note that increasing resistance in the circuit increases damping.
2
6. T =4 s
02 b 2
or
2
02 b 2
4
Also,
a
ln 0 ln 10 bt
an
or
1
b ln 10
t
224
2.3
= log 10 10 = 0.05 s 1
46 s
Hence,
02 (0.0025 2.4649) s 2
= 2.467 s 2
0 = 1.57 s 1
and
0 0 1.57
Q 15.7
2 2b 0.1
Terminal Questions
1. Since 0 e bt , we can write
1
b ln 0
t
Substituting the given data, we get
1 5
b ln
40 4
= 5.57 10 3 s 1
and
1
179.5 s 1
b
0 2Q 2 4000
2. Since Q , we can write = = 4.24 s
2 0 2 300
bt t /
Now, a a0 e a0 e
a
t ln 0 4.24 s ln 2 2.94 s
a
0
(iii) For a weakly damped system, Q 25 43.5 = 1088
2
225
1 1
Here, 3 6
10 8 s 2
LC 5 10 H 2 10
R2 (0.2) 2 2 2
And 400 2
4 L2 4 (5 10 3 ) 2 H 2 H
1 R2
Since , the discharge is oscillatory and has frequency
LC 4 L2
1 1 R2
2 = 1.59 10 3 Hz
2 LC 4 L
E E0 exp( 2bt )
= E0 exp( 2t / )
where b 1
is the relaxation time.
E0
When t / 2 , E
e
d
Also, Q
2
Hence,
2Q 2 6 10 4
= = 37.3 s
d 2 512 s 1
Thus, energy will reduce to 1 / e of its initial value in 18.7 s.
The number of oscillations made by the tuning fork in this time is given by
n d t
= 512 18.7
= 95.7 10 2
226
UNIT 4 FORCED OSCILLATIONS AND RESONANCE
Structure
4.1 Introduction Objectives
4.2 Differential Equation for a Weakly Damped Forced Oscillator
4.3 Solutions of the Differential Equation
Steady-state Solution
4.4 Effect of the Frequency of the Driving Force on the Amplitude and Phase of Steady-state
Forced Oscillations
Low Driving Frequency
Resonance Frequency
High Driving Frequency
4.5 Power Absorbed by a Forced Oscillator
4.6 Quality Factor
Q in Terms of Band Width: Sharpness of a Resonance
4.7 An LCR Circuit
4.8 Summary
4.9 Terminal Questions
4.10 Solutions
4.1 INTRODUCTION
In the previous unit we studied how the presence of damping affects the amplitude and the
frequency of oscillation of a system. However, in systems, such as a wall clock or an ideal LC
circuit, oscillations do not seem to die out. To maintain oscillations we have to feed energy to the
system from an external agent called a driver. In general, the frequencies of the driver and the
driven system may not match. But in steady-state, irrespective of its natural frequency, the system
oscillates with the frequency of the applied periodic force. Such oscillations are called forced
oscillations. However, when the frequency of the driving force exactly matches the natural
frequency of the vibrating system a spectacular effect is observed; the amplitude of forced
oscillations becomes very large and we say that resonance occurs. Do you know that Galileo was
the first physicist who understood how and why resonance occurs?
Resonances are desirable in many mechanical and molecular phenomena. But resonance can be
disastrous also; it can literally break an oscillating system apart. For instance, fast blowing wind
may set a suspension bridge in oscillation. If the frequency of the fluctuating force produced by
the wind matches the natural frequency of the bridge, it gains in amplitude and may ultimately
collapse. In 1940, the Tacoma Narrows bridge in Washington State collapsed within 4 months of
its being opened. Similarly, when the army marches on a suspension bridge, soldiers are
instructed to break step to avoid resonant vibrations. In practice, isolated systems are rare. In solid
state and molecular physics, two or more systems are coupled through interatomic forces. In an
electric circuit we have inductive and capacitative couplings. The oscillations of such systems
will be studied in the next unit.
227
In this unit we shall study, in detail, the response of a system when it is driven by an external
harmonic force.
Objectives
After studying this unit you should be able to:
establish the differential equation of a system driven by a harmonic force and solve it
analyse the response of the oscillator at different frequencies
compute resonance width and the quality factor of a forced oscillator, and
establish the differential equation for an LCR circuit under the influence of harmonic emf
and write its solution by drawing an analogy with a mechanical system.
F (t ) F0 cos t (4.1)
where F0 is a constant.
Let the mass be displaced from its equilibrium position and then released. At any instant, it is
dx
subject to (i) a restoring force, kx , (ii) a damping force, and (iii) a driving force,
dt
F (t ) F0 cos t .
d 2x dx
m 2
kx F0 cos t (4.2)
dt dt
Dividing by m and rearranging terms, the equation of motion of a forced oscillator takes the form
d 2x dx
2
2b 02 x f 0 cos t (4.3)
dt dt
228
where 2b / m , 02 k / m and f 0 F0 / m is a measure of the driving force.
You may now ask: Does this equation apply only to a mass on a spring? No, it applies to any
oscillator whose natural frequency is 0 and is subject to a harmonic driving force.
You will note that Eq. (4.3) is an inhomogeneous second order linear differential equation with
constant coefficients. We will now solve this equation to learn about the motion of a forced
oscillator.
x(t ) x1 (t ) x2 (t )
where x1 (t ) is a solution of the equation obtained by replacing the RHS of Eq. (4.3) by zero.
On substituting this result in Eq. (4.3), you will find that x 2 (t ) satisfies the equation
d 2 x2 dx
2
2b 2 02 x2 f 0 cos t
dt dt
It is thus clear that ( x1 x2 ) is the complete solution of Eq. (4.3). In your course on differential
equations you must have learnt that x1 is called the complementary function and x 2 (t ) is called
the particular integral.
You may recall that when there is no driving force, the displacement of a weakly damped
( b 0 ) system at any instant is given by Eq. (3.19):
x1 (t ) a0 e bt cos( d t )
Obviously this complementary function decays exponentially and after some time it will
disappear. That is why it is also referred to as the transient solution. In the transient state, the
system oscillates with some frequency which is other than its natural frequency or the frequency
of the driving force.
After a sufficiently long time (t ) , natural oscillations of the spring-mass system will
disappear due to damping. However, we know that the general solution of Eq. (4.3) will not decay
229
with time. That is, the system will oscillate with the frequency of the driving force. The system is
then said to be in the steady-state. We will now obtain the steady-state solution of Eq. (4.3).
where a and are unknown constants. By comparing Eqs. (4. 1) and (4.4) you will note that the
driving force leads the displacement in phase by an angle .
To determine a and we differentiate Eq. (4.4) twice with respect to time. This gives,
dx2
a sin( t )
dt
d 2 x2
and 2
a 2 cos( t )
dt
We know that both cos t and sin t never simultaneously become zero; when one vanishes,
the other takes a maximum value. Therefore, Eq. (4.5) can be satisfied only when both terms
within the square brackets become zero separately, i.e.
( 02 2 )a cos 2ab sin f 0 (4.6a)
( 02 2 )a sin 2ab cos 0 (4.6b)
Eq. (4.6b) readily gives the phase by which the driving force leads the displacement:
2b
tan 1 (4.7a)
2
2
0
230
The amplitude of steady-state displacement can be determined from Eq. (4.6a) once we know the
values of sin and cos . To get these values we construct the so-called acoustic impedance
triangle, as shown in Fig. 4.2. We can readily write
2b
sin
[( ) 4b 2 2 ]1 / 2
2
0
2 2
( 02 2 )
and cos
[( 02 2 ) 2 4b 2 2 ]1 / 2
Using these values of sin and cos in Eq. (4.6a) and rearranging terms, we get
f0 F0
a (4.7b)
[( ) 4b ]
2
0
2 2
2 2 1/ 2
m[( 0 ) 4b 2 2 ]1 / 2
2 2 2
Thus, we find that the steady-state amplitude of forced oscillations depends on (i) amplitude and
angular frequency of the driving force, (ii) mass and the natural angular frequency of the
oscillating system and (iii) the damping constant.
Putting this value of a in Eq. (4.4) we can write the steady-state solution of Eq.(4.3) as
F
x2 (t ) 0
cos( t ) (4.8)
m[( ) 4b 2 2 ]1 / 2
2
0
2 2
The important point to note here is that the steady-state solution has the frequency of the driving
force and its amplitude is constant. Moreover, its phase is also defined completely with respect to
the driving force. Therefore, it does not depend on the initial conditions. In other words, the
motion of a driven system in steady-state is independent of the way we start the oscillation.
F0 cos( t )
x(t ) a0 e bt cos( d t )
m[( 2 ) 2 4b 2 2 ]1 / 2
2
0
231
i.e. the complete general solution of Eq. (4.3) are shown in Fig. 4.3. The contribution of the
transient part diminishes with time and ultimately disappears completely. The time for which
transients persist is determined by b and hence by the damping factor . The greater the value of
b, the more quickly do the transients die out.
For an undamped system, the steady-state solution is obtained by putting b = 0 in Eqs. (4.7a) and
(4.8). This gives
0
F0
and x2 (t ) cos t (4.9)
m( 2 )
2
0
That is, the driving force and the displacement are in phase ( 0) . From this we may conclude
that phase lag is essentially a consequence of damping. We further note that if the frequency of
the driving force equals the frequency of the undamped oscillator, its amplitude will become
infinitely large. Then resonance is said to occur.
You may now ask: Do we observe infinitely large amplitude in practice? No, the amplitude is
finite since some damping is always present in every system.
Fig. 4.3 Time variation of the transient solution, steady-state solution and the general solution of
Eq. (4.3) for a weakly damped system.
232
4.4.1 Low Driving Frequency ( 0 )
To know the behaviour of a( ) at low driving frequencies, we first rewrite
Eq. (4.7 b) as
f0
a( ) 1/ 2
2
2
4b 2 2
02 1 2
0 02
For 0 , we note that the ratio 2 / 02 will be much less than one. So, we neglect terms
containing 2 / 02 . This gives
f0 F0 F
a( ) 0 (4.10a)
2
0 m 02
k
Thus, at very low driving frequencies, the steady-state amplitude of the oscillation is controlled
by the stiffness constant and the magnitude of the driving force.
2b
tan 0 for 1 (4.10b)
2
0
2
0
That is, the driving force and the steady-state displacement are in phase.
f0
a( 0 ) (4.11a)
2b 0
From this we note that at resonance the amplitude depends upon the damping; it is inversely
proportional to b. That is why in actual practice the amplitude never becomes infinite.
tan
so that
/2 (4.11b)
This means that the driving force and the displacement are out of phase by / 2 . You may be
thinking that the value of a( ) given by Eq. (4.11a) is maximum. This however is not true.
Why? To answer this, let us maximise a( ) . That is, differentiate Eq. (4.7b) with respect to
233
and set the resulting expression equal to zero. The frequency at which the first derivative becomes
zero and the second derivative is negative gives the correct answer:
da( ) d f0
2
d d [( 0 2 ) 2 4b 2 2 ]1 / 2
f 0 [4 ( 02 2 ) 8b 2
=
2[( 02 2 ) 2 4b 2 2 ]1 / 2
f 0 [4 ( 02 2 ) 8b 2 ]
= 0 for r
2[( 02 2 ) 2 4b 2 2 ]1 / 2
r2 02 2b 2 0
r ( 02 2b 2 )1 / 2 (4.12)
The root corresponding to the negative sign is physically meaningless and is ignored.
For a( ) to be maximum, its second derivative with respect to should be negative. You can
d 2a
easily verify that at r ( 02 2b 2 )1 / 2 is negative. Thus, we can conclude that the
d 2 r
peak value of amplitude is attained at a frequency slightly below 0 . The shift is caused due to
damping. We can visualize it as follows: When the driver imparts maximum push, the driven
system does not accept it instantly due to a finite phase difference between x(t ) and F (t ) .
On substituting for 0 from Eq. (4.12) in Eq. (4.7b) and simplifying the resulting expression, we
get the peak value of steady-state amplitude:
f0
a max (4.13)
2b 02 b 2
When at a particular frequency, the amplitude of the driven system becomes maximum, we say
that amplitude resonance occurs. The frequency r is referred to as the resonance frequency. It
is instructive to note that r is less than 0 as well as d 02 b 2 .
234
For 0 we rewrite Eq. (4.7b) as
f0
a( ) 1/ 2
2
2
4b 2
2 1 02 2
and neglect terms containing 02 / 2 as well as (2b / ) 2 , as they are both much smaller than
unity. Then the amplitude of resulting vibration is given by
f0
a( ) (4.14a)
2
That is, at high frequencies the amplitude decreases as 1 / 2 and ultimately becomes zero.
2b 2b
tan 0
( )
2
2
0
or (4.14b)
This means that at high frequencies the driving force and displacement are out of phase by .
235
Fig. 4.4 Frequency variation of (a) steady-state amplitude, and (b) phase of a forced oscillator
4.5 POWER ABSORBED BY A FORCED OSCILLATOR
You now know that every oscillating system loses energy in doing work against damping. But
oscillations of a forced oscillator are maintained by the energy supplied by the driving force. It is,
therefore, important to know the average rate at which energy must be supplied to the system to
sustain steady-state oscillations. So, we now calculate the average power absorbed by the
oscillating system.
dx2 (t ) F0
v
dt m[( 0 ) 4b 2 2 ]1 / 2
2 2 2
is the phase difference between velocity and the applied force. On substituting for F (t ) and v
from Eqs. (4.1) and (4.15), respectively we find that the instantaneous power absorbed by the
oscillator is given by
Since cos( t ) cos t cos sin t sin , we can rewrite the expression for
instantaneous power as
From this we can easily calculate the average power absorbed over one cycle:
1
P F0 v0 sin sin 2 t F0 v0 cos cos 2 t (4.16)
2
From Unit 1 you may recall that sin 2 t = 0 so that the first term on the RHS of Eq. (4.16)
drops out. Also cos 2 t 1 / 2 . Then Eq. (4.16) reduces to
236
1 1
P F0 v0 cos F0 v0 sin (4.17)
2 2
On substituting for sin from Fig. 4.2 and v 0 from Eq. (4.15a) in Eq. (4.17), we get
bF02 2
P (4.18)
m [( 0 ) 4b ]
2 2 2 2 2
From Eq. (4.17) we note that the average power absorbed by a forced oscillator will be maximum
when sin 1 cos , i.e., / 2 ( 0) . This happens for 0 . This happens for
0 . Using this result in Eq. (4.18), we get
1
P max F02 (4.19)
4bm
That is, the peak value of average power absorbed by a maintained system is determined by
damping and the amplitude of the driving force. The frequency variation of P is shown in
Fig. 4.5
Fig. 4.5. Frequency variation of average power, a., and a>: correspond to hair-power points.
It is important to note that unlike the case of amplitude resonance, maximum average Forced
Oscillations and Resonance power is transferred at the natural frequency of the system. This
arises because velocity and driving force are in phase.
237
You can use the same definition to calculate Q of a forced oscillator once you know E and
P .
SAQ1
Show that the average energy of a forced oscillator is
1
E m( 2 02 )a 2
4
2 02
Q
4b
Another equivalent and more useful interpretation of the quality factor is in terms of amplitudes.
The Q factor is defined as the ratio of the amplitude at resonance to the amplitude at low
frequencies ( 0) . Using this definition, the value of the quality factor can be calculated
rather easily on dividing Eq. (4.13) by Eq. (4. l0a).
a max f0 02 02
Q (4.20a)
a( 0) 2b( 02 b 2 )1 / 2 f 0 2b( 02 b 2 )1 / 2
If damping is small, b 2 02 and the expression for the quality factor reduces to
0 0
Q (4.20b)
2b 2
SAQ 2
Using Eq. (4.20b), show that the amplitude and phase of a weakly damped forced oscillator can
be expressed as
0 /
a( ) a0 1/ 2
0 2 1
2
0 Q
1/ Q
and tan
0
0
where a0 F0 / m 02 .
238
For different values of Q frequency variation of a( ) and ( ) based on these equations is
shown in Fig.4.6. We observe that as Q increases (i.e., damping decreases), the value of
a( ) increases.
Fig. 4.6 (a) Amplitude as a function of driving frequency for different values of Q,
(b) Phase difference as a function of driving frequency for different values of Q.
SAQ3
1 F02 0
Express P in terms of Q and show that P Q
2 k
To calculate the frequency at which average power drops to half its maximum value we can write
from SAQ 3
239
On simplification we can write
2 02
( 02 2 )
Q2
so that
2 02
( )
2
0
2
Q2
This equation has 4 roots. Of these, two roots correspond to negative frequencies and are
physically unacceptable. The other two acceptable roots are
0
1/ 2
1
1 0 1
2
2Q 4Q
0
1/ 2
1
and 2 0 1
2
(4.22)
2Q 4Q
Obviously, the second of these roots is greater than 0 and the other root is smaller than 0 .
This is illustrated in Fig. 4.5.
0
2 1 2
Q
From Eq. (4.23) it is clear that a high Q system has small bandwidth and the resonance is said to
be sharp. On the other hand, a low Q system has a large bandwidth and the resonance is said to
be flat. This is illustrated in Fig. 4.6. Thus, the sharpness of resonance refers to the rapid rate of
the fall of power with frequency on either side of resonance. We measure it in terms of the Q -
value of the system. The Q factor has its greatest importance in reference to electrical circuits
which we will discuss now.
SAQ 4
Calculate the energy stored in a mass of 0.1 kg attached to a spring. The mass is oscillating with
an amplitude of 5 cm and is in resonance with a driving force of frequency 30 Hz. If the Q factor
is 100, calculate the power loss.
240
From Unit 3 we know that in an LCR circuit, charge oscillations die out because of power losses
in the resistance. What changes do you expect in this behaviour when a source of alternating
e.m.f. of frequency is introduced? To answer this question, let us consider Fig, 4.7. Let I be
the current in the circuit at a given time. Then, the applied emf is equal to the sum of the potential
differences across the capacitor, resistor and the inductor. Then Eq. (3.35) modifies to
q dI
RI L E0 cos t (4.24)
C dt
dq
Since I , this equation can be rewritten as
dt
d 2q dq q
L 2 R E0 cos t (4.25)
dt dt C
d 2 q R dq q E
2
0 cos t (4.26)
dt L dt LC L
In this form Eq. (4.26) is similar to Eq. (4.3). Hence its steady-state solution can be written by
analogy. For a weakly damped system, the charge on capacitor plates at any instant of time is
given by
E0 / L
q 1/ 2
cos( t ) (4.27)
1 2 R
2
LC L
1 R2
where 2 is the angular frequency of oscillation and
LC 4 L
(R / L)
tan (4.28)
1
2
LC
241
defines the phase with respect to applied emf.
The current in the circuit is obtained by differentiating Eq. (4.27) with respect to t . The result is
E0
I cos(t ) (4.29)
2
1
R 2 L
C
From Eq. (4.29) we note that current in a LCR circuit is a function of the frequency.
1
When L , the circuit is capacitive in nature and we can write
C
2
1 1
L 2 2
C C
Thus, if we are working at low frequencies and R is also small, the current amplitude will be
small. What will be its magnitude for 0? In this limit I 0 and leads the applied emf by
/ 2.
1
As the driving frequency increases, the reactance L decreases and current amplitude
C
increases. When
1
L (4.31)
C
the term under the radical sign in Eq. (4.29) becomes minimum; equal to R . Then the current
attains its peak value I 0 E0 / R and the circuit is said to resonate with frequency
1
r (4.32)
2 LC
242
At resonance, the current and applied emf are in phase. When the driving frequency is high, the
circuit will be inductive and the current lags behind emf by / 2 .
For different values of R , the frequency variation of peak current and phase is shown in Fig. 4.8.
You will observe that the lower the resistance, the higher the peak value of the current and the
sharper the resonance.
Fig. 4.8 Frequency variation of peak current and phase for different values of R in
a driven LCR circuit
The power in an electric circuit is defined as the product of current and the emf. For an LCR
circuit, we can write
P EI E0 I 0 cos t cos( t )
2
1
where I 0 E0 / R L
2
C
Using the formula 2 cos A cos B cos( A B) cos( A B) , we can rewrite the above
expression for power as
E0 I 0
P [cos cos(2 t )]
2
Power averaged over one complete cycle is obtained by noting that cos(2 t ) 0 .
Hence
E0 I 0
P cos
2
= Erms I rms cos
243
where Erms E0 / 2 and I rms I 0 / 2 are, respectively the root mean square values of emf
and current. Since P varies with cos , it is customary to call cos the power factor.
0 L
A (4.34)
R
where 0 1 / LC
You can verify that the bandwidth of power resonance curve for an LCR circuit is given by
2 R 0
2 1 (4.35)
L Q
so that
Frequency at resonance
Q
Full width at half power points
The Q of a circuit determines its ability to select a narrow band of frequencies from a wide range
of input frequencies. This, therefore, acquires particular importance in relation to radio receivers.
Signals of various frequencies from all stations are present around the antenna. But the receiver
selects just one particular station to which we wish to tune and discards others. Normally, radio
receivers operating in MHz region have Q values of the order of 10 2 to 10 3 . Microwave
cavities have Q values of the order of 10 5 .
4.8 SUMMARY
When a harmonic force F F0 cos t is impressed upon a damped harmonic oscillator,
the oscillator executes forced oscillations. The differential equation of motion of a driven
oscillator is
d 2x dx
2
2b 02 x F0 cos t
dt dt
k F
where 2b = , 0 and f 0 0
m m m
f0
a
[( ) 4b 2 2 ]1 / 2
2
0
2 2
and
2b
tan 1
2
2
0
244
are steady-state amplitude and phase respectively.
LC L
R / L
with tan
1
2
LC
245
3. An alternating potential of frequency 10 5 Hz and amplitude 1.2 V is applied to a series
LCR circuit. If L = 0.5 mH and R = 40 , find the value of the capacitance C to get
resonance. Also calculate the rms value of this current.
4.10 SOLUTIONS
SAQ 1
K .E. (1 / 2)m v 2
(1 / 2)m 2 a 2 sin 2 ( t ) (1 / 4)m 2 a 2
U (1 / 2)k x 2
(1 / 2)m 02 a 2 cos 2 ( t ) (1 / 4)m 02 a 2
Therefore, the time-averaged energy is
E (1 / 4)m( 2 02 )a 2
The average energy dissipated per second is v 2 = mb 2 a 2 .
By definition,
average energy stored in one cycle
Q 2
average energy dissipated in one cycle
average energy stored in one cycle
2
time period average energy dissipated per second
m( 02 2 ) a 2
2
4Tmb 2 a 2
2 2
0
4b
SAQ 2
From Eq. (4.7b) we recall that the amplitude of a weakly damped forced oscillator is given by
F
a 0
m[( ) 4b 2 2 ]1 / 2
2
0
2 2
F0 / m
1/ 2
2
4b 2
0 0 2
0 0
F0
If we put a 0 and use Eq. (4.20), we get the required result:
m 02
0 /
a( ) a0 1/ 2
0 2 1
2
0 Q
246
Similarly, from Eq. (4.7a), we recall that
2b
tan
2
2
0
2b
0 0
0
1/ Q
or tan
0
0
SAQ 3
From Eq. (4.18)
bF02 2
P
m [( 02 2 ) 2 4b 2 2 ]
0
Putting Q , we get
2b
0 F02 2
P
2mQ [( 02 2 ) 2 ( 2 02 / Q 2 )]
At 0 , the denominator in the parenthesis will become minimum and the average power
absorbed by the oscillator becomes maximum:
F02 Q
P max
2m 0
1 F02 0 Q
=
2 m 02
1 F02 0 Q
=
2 k
SAQ 4
From SAQ 1, you would recall that the average energy of a weakly damped oscillator is given by
m 2
E ( 0 2 ) a 2
4
At resonance, 0 , and the expression for average energy reduces to
E (1 / 2)m 02 a 2
= (1 / 2) 0.1 kg (2 30 s 1 ) 2 (5 10 2 ) m) 2
= 4.44 J
247
2 4.44
Average energy dissipated in (1 / 30) s
100
= 0.28 J
Average energy dissipated per second = 30 0.28 J = 8.4 J.
Terminal Questions
1. For free oscillations, the differential equation is
mx kx x
Substituting m = 0.1 kg, k = 100 Nm 1 and = 5 N m s 1 we get
d 2x dx
0.1 2
100 x 5
dt dt
2
d x dx
or 2
50 1000 x 0
dt dt
2 2
Period T s 0.32 s .
[ 0 b ]
2 2 1/ 2
[1000 625]1 / 2
a max
2. Let r be the value of the angular frequency when a .
2
Using Eq. (4.7b) and (4.13), we get
f0 f0
[( ) 4b ]
2
0
2 2
2 2 1/ 2
4b( 0 b 2 )1 / 2
2
( 02 2 ) 2 4b 2 2 16b 2 ( 02 b 2 )
248
( r2 2 2b 2 ) 2 4b 2 2 16b 2 ( r2 b 2 )
= 5.1 10 9 F
E rms E 1.2V
I rms 0
R R 2 4 2
= 0.21 A (at resonance Z R )
= 95.5 V
249
UNIT 5 COUPLED OSCILLATIONS
Structure
5.1 Introduction Objectives
5.2 Oscillations of Two Coupled Masses
Differential Equation
Normal Co-ordinates and Normal Modes
Modulation of Coupled Oscillations
Energy of Two Coupled Masses
General Procedure for Calculating Normal Mode Frequencies
5.3 Normal Mode Analysis of other Coupled Systems
Two Coupled Pendulums
Inductively Coupled LC-circuits
5.4 Longitudinal Oscillations of N Coupled Masses: The Wave Equation
5.5 Summary
5.6 Terminal Questions
5.7 Solutions
5.1 INTRODUCTION
In this block so far you have studied isolated (single) oscillating systems such as a spring-mass
system, a pendulum or a torsional oscillator. In nature we also come across many examples of
coupled oscillators. We know that atoms in a solid are coupled by interatomic forces. In
molecules, say the water molecule, two hydrogen Moms are coupled to an oxygen atom while in a
carbon dioxide molecule oxygen atoms are coupled to one carbon atom. In all these cases,
oscillations of one atom are affected by the presence of other atom(s). In radio and TV
transmission, we use electrical circuits with inductive/ capacitative couplings. Therefore, it is
important to extend our study of preceding units to cases where such simple systems are coupled.
We begin this unit with a study of longitudinal oscillations of coupled masses. Do you expect the
motion of these masses to be simple harmonic? You will learn that their motion is not simple
harmonic. But it is possible to analyse it in terms of normal modes, each of which has a definite
frequency and represents SHM. The presence of coupling leads to exchange of energy between
two masses. To illustrate this further, we will determine, by analogy, normal mode frequencies of
two coupled pendulums and two inductively coupled LC circuits. This analysis will then be
extended to N coupled oscillators. When N becomes very large, i.e. we have a homogeneous
medium, exchange of energy leads to the phenomenon of wave«motion.
In the next unit you will learn the details of wave propagation with particular reference to waves
in strings, liquids and gases.
Objectives
After studying this unit you should be able to:
250
5.2 OSCILLATIONS OF TWO COUPLED MASSES
To analyse the effect of coupling we start again with the model spring-mass system. We consider
two such identical systems connected (coupled) by a spring, as shown in Fig. 5.la. In this system
we have two equal masses attached to springs of stiffness constant k' and coupled to each other by
a spring of stiffness constant k. In the equilibrium position, springs do not exert any force on
either mass. The motion of this system will depend on the initial conditions. That is, the motion
may be transverse or longitudinal depending on how the masses are disturbed. For simplicity, we
first consider longitudinal motion of these two coupled masses.
We pull one of the masses longitudinally and then release H. The restoring force will tend to
bring it back to its equilibrium position. As it overshoots the equilibrium . mark, the coupling
spring will pull the other mass. As a result both masses start oscillating longitudinally. This
means that motion imparted to one of the two coupled masses is not confined to it only; it is
transmitted to the other mass as well. We now establish the equation of motion of these masses.
Let X A and X B be the coordinates of the centre of the masses A and B respectively. When mass
B is displaced towards the right and then released, mass A will also get pulled towards the right
due to the coupling spring. The coupled system would then start oscillating. Suppose X A and
X B are the instantaneous positions of masses A and B respectively. Then their displacements
from their respective equilibrium positions are given by
x2 x B X B and x1 x A X A
Now at any instant of time during oscillation, the forces acting on mass A are
(i) restoring force : k ' ( x A X A ) k ' x1 ; and
(ii) a coupling force : k ( x B x A ) ( X B X A ) k ( x2 x1 )
251
We are here assuming that the masses are moving on a frictionless surface. By Newton's second
law, the equation of motion of mass A is thus given by
d 2 xA
m k ' ( x A X A ) k[ x B X B x A X A ]
dt 2
d 2 (xA X A ) d 2 x1
or m m 2 k ' x1 k ( x2 x1 ) (5.1)
dt 2 dt
dX A
since 0.
dt
d 2 x1
2
02 x1 s2 ( x2 x1 ) 0 (5.2)
dt
k' k
where 02 and s2
m m
d 2 x2
m k ' x2 k ( x2 x1 ) (5.3)
dt 2
Let us pause for a minute and ask: Do Eqs. (5.2) and (5.4) represent simple harmonic motion?
No, we cannot, in general, identify the motion described by these equations as simple harmonic
because of the presence of the coupling term s2 ( x2 x1 ) . This means that the analysis of
previous units will not work since these equations are coupled in x1 and x 2 . The question now
arises: How to solve these equations? These equations will have to be solved simultaneously. For
this purpose we first add Eqs. (5.2) and (5.4) to obtain
d2
2
( x1 x2 ) 02 ( x1 x2 ) 0 (5.5a)
dt
Next we subtract Eq. (5.4) from Eq. (5.2) and rearrange terms. This gives
d2
( x1 x2 ) ( 02 2 s2 )( x1 x2 ) 0 (5.5b)
dt 2
By looking at Eqs. (5.5a) and (5.5b) you will recognise that these are standard equations for
SHM. This suggests that if we introduce two new variables defined as
252
1 x1 x2 (5.6a)
and
2 x1 x2 (5.6b)
the motion of a coupled system can be described in terms of two uncoupled and independent
equations:
d 21
2
121 0 (5.7)
dt
and
d 2 2
2
22 2 0 (5.8)
dt
We therefore find that new co-ordinates 1 and 2 have decoupled Eqs, (5.2) and (5.4) into two
independent equations which describe simple harmonic motions of frequencies 1 and 2 and
2 1 . The new coordinates are referred to as normal coordinates and simple harmonic
motion associated with each coordinate is called a normal mode. Each normal mode has its own
characteristic frequency called the normal mode frequency.
1 (t ) a1 cos(1t 1 ) (5.11)
and
2 (t ) a2 cos( 2 t 2 ) (5.12)
where a1 and a 2 are the amplitudes of normal modes and, 1 and 2 are their initial phases.
Since x1 (t ) = (1 2 ) / 2 , we can write the displacement of mass A as
1
x1 (t ) [a1 cos(1t 1 ) a 2 cos( 2 t 2 )] (5.13)
2
253
1
x2 (t ) [a1 cos(1t 1 ) a 2 cos( 2 t 2 )] (5.13)
2
The constants a1 , a 2 , 1 , 2 are fixed by the initial conditions. Once we know these, we can
completely determine the motion of the coupled masses.
SAQ 1
Solve Eqs. (5.13) and (5.14) subject to the following initial conditions:
dx1 dx
(A) x1 (0) a , x2 (0) a , 0 and 2 0
dt t 0 dt t 0
dx dx
(B) x1 (0) a , x2 (0) a , 1 0 and 2 0
dt t 0 dt t 0
On solving this SAQ you will observe that when both masses are initially given the same
displacement to the right and then released, their displacements are equal, i.e. x1 (t ) x2 (t ) , or
2 = 0 at all times. The motion is described by Eq. (5.7) and the normal mode frequency is the
same as that of the uncoupled masses. This means that coupling has no influence and both masses
oscillate in phase. In this mode of vibration, the coupling spring is neither stretched nor
compressed (and is as good as not being there), as shown in Fig. 5.2a.
Fig. 5.2. Normal modes (a) When two coupled masses are given equal displacement in the same
direction, (b)When two coupled masses are pulled together equally.
When two coupled masses are initially pulled, equally towards each other and then released, the
displacements are equal but out of phase by , i.e., x1 x2 or 1 = 0 (Fig 5.2b). The normal
mode frequency will be higher than that of the uncoupled masses ( 2 1 ) . This means that the
coupling spring is either compressed or stretched and we say that coupling is effective. We thus
conclude that normal coordinates allow us to write the equation of motion of a coupled system
into a set of linear differential equations with constant coefficients. Each equation contains only
one dependent variable. Moreover, the motion of a coupled system may be regarded as a
superposition of its possible normal modes.
254
5.2.3 Modulation of Coupled Oscillations
In the above discussion we assumed that the two coupled masses are pulled equally in the same
direction or in opposite directions. What will happen if only one of them is pulled and then
released? To understand this we have to solve Eqs. (5.13) and (5.14). Suppose the initial
condition is as follows:
dx1 dx 2
x1 (0) 2a , x2 (0) 0 , 0 and 0 (5.15)
dt t 0 dt t 0
You will find that the displacements of two coupled masses are given by
Expressing the sum (difference) of two cosine functions into their product, these equations can be
rewritten in a physically more familiar form:
1 2 1
x1 (t ) 2a cos 2 t cos t (5.17)
2 2
and
1 2 1
x2 (t ) 2a sin 2 t sin t (5.18)
2 2
You would recall that Eq. (5.17) is essentially the same as Eq. (2.19) obtained for modulated
oscillations. As before, we define av (1 2 ) / 2 as the average angular frequency and
mod ( 2 1 ) / 2 as the modulated angular frequency.
Then Eqs. (5.17) and (5.18) represent modulated oscillations respectively defined by
What is the phase difference between the displacements of the two masses? Since sine and cosine
functions differ by / 2 , the phase difference between the displacements of the coupled masses
is / 2 . The same is true of modulated amplitudes as well.
The displacement-time graphs for the two masses are shown in Fig. 5.3. We observe at t = 0, the
amplitude of A is maximum while that of the mass at B is zero. With time, the amplitude of A
decreases and becomes zero at t T / 4 while that of B becomes a maximum. After t T / 4 ,
255
this trend is reversed for the succeeding quarter of the period. This process will repeat itself
indefinitely provided damping is not present.
1
E1 m av
2 2
a mod (t ) 2ma 2 av
2
cos 2 mod t
2
and
1
E2 m av
2 2
bmod (t ) 2ma 2 av
2
cos 2 mod t
2
The total energy of two masses coupled through a spring which stores almost no energy is given
by
E E1 E2 2ma 2 av
2
(5.24)
E
E1 [1 cos( 2 1 ) t ] (5.25a)
2
and
E
E2 [1 cos( 2 1 ) t ] (5.25b)
2
256
These equations show that at t 0 , E1 E and E2 0 . That is, to begin with, mass at A
possesses all the energy. As time passes, energy of mass at A starts decreasing. But mass at B
begins to gain energy such that the total energy of the system remains constant.
T 2 /( 2 1 )
We now illustrate this procedure for two unequal masses m1 and m 2 coupled through a spring of
force constant k. The equations of motion of two coupled masses are
x1 a1 cos( t )
and
x2 a2 cos( t )
2 k k
0 2 x1 x2 (5.27a)
m1 m
257
and
2 k k
0 2 x2 x1 (5.27b)
m2 m
x1 k / m1
x2 2 k
0 2
m1
2 k
0 2
k / m1
=
m2
2 k k / m2
0 2
m1
On cross-multiplying, we get
k k k k
4 2 02 2 02 02 0
m1 m2 m1 m2
12 02 (5.28a)
and
1 1
22 02 k (5.28b)
m1 m2
You will note that for m1 m2 , these results reproduce Eqs. (5.9) and (5.10).
258
5.3.1 Two Coupled Pendulums
Let us consider two identical simple pendulums A and B having bobs of equal mass, m , and
suspended by strings of equal length l , as shown in Fig. 5.4 (a). The bobs of the two pendulums
are connected by a weightless, spring of force constant k. In the equilibrium position, the distance
between the bobs is equal to the length of the unstretched spring.
Suppose that both bobs are displaced to the right from their respective equilibrium positions. Let
x1 (t ) and x2 (t ) be the displacements of these bobs at time t , as shown in Fig. 5.4 (b). The
tension in the coupling spring will be k ( x1 x2 ) . It opposes the acceleration of A but will
support the acceleration of B. For small amplitude approximation, we recall from Unit 1 that the
equation of motion of a simple pendulum is
d 2x mg
m 2
x
dt l
Fig. 5.4 Two identical pendulums simple together (a) Equilibrium configuration (b) Instantaneous
configuration
d 2 x1 mg
m x1 k ( x1 x2 )
l
2
dt
and
d 2 x2 mg
m x2 k ( x1 x2 )
l
2
dt
The term k ( x1 x2 ) arises due to the presence of coupling. Dividing throughout by m and
rearranging terms, we get
259
d 2 x1
2
02 x1 s2 ( x1 x2 ) 0 (5.29a)
dt
and
d 2 x2
2
02 x2 s2 ( x1 x2 ) 0 (5.29b)
dt
You will recognize that these equations are respectively identical to Eqs. (5.2) and (5.4). Thus the
entire analysis of preceding sections applies and we can describe the motion of coupled
pendulums by drawing analogies. The normal modes of this system are shown in Fig. 5.5. In
mode 1 ( x1 x2 ) , the bobs are in phase and oscillate with frequency 1 0 g / l . But in
mode 2 ( x1 x2 or x2 x1 ) , the bobs are in opposite phase and oscillate with frequency
2 [ 02 2 s2 ]1 / 2 [( g / l ) 2(k / m)]1 / 2 .
Fig. 5.5 Normal modes of a coupled pendulum (a) In-phase normal mode (b) Out-of-phase
normal mode.
SAQ 2
The kinetic and potential energies of two identical coupled pendulums are given by
m
K .E. [( x1 ) 2 ( x 2 ) 2 ]
2
and
1 mg 2 1
U ( x1 x2 ) k ( x1 x2 )
2 2
2 l 2
dxi
where x i ( i = 1,2). Express these in terms of normal coordinates.
dt
260
On solving this SAQ, you will observe that
m m
K .E. (1 ) 2 ( 2 ) 2
4 4
and
1
U (m1212 m 22 22 )
4
We can rewrite these expressions in a more elegant form by defining normal coordinates as
m m
1 ( x1 x2 ) and 2 ( x1 x2 ) (5.30)
2 2
SAQ 3
Using the definition given in Eq. (5.30), calculate the total energy of a system of two coupled
pendulums in terms of normal coordinates. At any instant, 1 1.5 10 3 m kg and
2 0.5 10 3 m kg . Calculate x1 (t ) and x2 (t ) at the same instant. Given m = 0.1 kg.
In the above discussion you have learnt to calculate normal mode frequencies of two pendulums
whose bobs are coupled. Will these frequencies remain the same if the strings of these pendulums
were coupled, as shown in Fig. 5.6 (a). To discover the answer to this question, we consider the
configuration shown in Fig 5.6 (b). At any time t , let the change in the length of the spring be
x2 ' x1 ' (d / l )( x2 x1 ) , where x1 and x 2 are displacements of the bob from their equilibrium
positions and d is the distance between the points of suspension and coupling. Hence, the
restoring force in the spring is given by
kd
F k ( x2 ' x1 ' ) ( x2 x1 )
l
Fig. 5.6 Two identical simple pendulums whose strings are coupled by a spring (a) Equilibrium
configuration (b) Instantaneous configuration
261
SAQ 4
Write down the equations of motion of two pendulums coupled at a distance d from the point of
suspension. Compute the normal mode frequencies by analogy. You will find that the frequencies
of the normal modes are given by
g
1
l
and
g 2kd 2
2
l ml 2
This shows that 2 depends on the distance between the points of suspension and coupling.
Obviously, for d l , the expression for 2 reduces to that for coupled bobs.
Two electrical circuits are said to be inductively coupled when a change in the
magnetic flux linked with one circuit induces emf (and hence gives rise to a
current) in the other circuit. The coupling coefficient is given by
M / L1 L2 , where M is mutual inductance and L1 and L2 are self-
inductances of two coupled circuits.
Let I and I 2 be the instantaneous values of currents in the two circuits. The equation of motion
of charge in circuit 1 is obtained by modifying Eq. (1.36) as
262
q1 dI dI
L1 1 M 2
C1 dt dt
where MdI 2 / dt is the emf produced in circuit 1 due to current I 2 in the second circuit.
Obviously, it tends to increase I 1 . Similarly, for circuit 2 we can write
q2 dI dI
L2 2 M 1
C2 dt dt
dq dI d 2 q
Since I and , we can rewrite these equations as
dt dt dt 2
d 2 q1 M d 2 q2
p q1
2
(5.31a)
dt 2 L1 dt 2
and
d 2 q2 M d 2 q1
2
q
s 2 (5.31b)
dt 2 L2 dt 2
Eqs. (5.31a) and (5.3Ib) are two coupled equations. To find normal mode frequencies, we write
q1 A cos( t )
M 2
( 2p 2 )q1 q2
L1
and
M 2
( s2 2 )q 2 q1
L2
M 2 L2 ( s2 2 )
L1 ( 2p 2 ) M 2
M2 4
( 2p 2 ) ( s2 2 ) 2 2 (5.32)
L1 L2
263
Eq. (5.32) is quadratic in 2 ; its roots give us normal mode frequencies. For simplicity, we
assume that the circuits are identical so that their natural frequencies are equal, i.e.
p s 0 , say, then
( 02 2 ) 2 2 4
or
02 2 2
so that
0
1
The acceptable normal mode frequencies are those values of which correspond to positive
roots and are given by
0
(5.33a)
1
and
0
(5.33b)
1
SAQ 5
Two identical inductively coupled circuits, each having a natural frequency of 600 Hz, have
coupling coefficient 0.44. Calculate the two normal mode frequencies.
264
d 2 n
m k ( n n1 ) k ( n n1 ) (5.34)
dt 2
From Unit 1 you may recall that the spring constant k is defined as restoring force per unit
extension. So we can write k = F/d, where d is extension in the spring. Using this relation in Eq.
(5.34), we get
d 2 n F n1 n n n1
(5.35)
dt 2 m d d
Now let us assume that N . The separation between any two consecutive masses will
become very small. That is, we will have a continuous distribution of masses. Then we can
replace d by x and Eq. (5.35) becomes
d 2 n F n1 n n n1
dt 2 m x x
or
d 2 n F
m x n1 x n
2
(5.36)
dt
If the n th mass is located at a distance x from the origin, then in the limit x 0 , we have
d 2 n F d d
m dx n1 dx n
2
(5.37)
dt
This means that is now a function of t a well as x . We know that any continuous function
f ( x x) can be expressed in terms of the function defined at x and its derivatives using the
following expansion:
f ( x, t ) 1 2 f ( x, t )
f ( x x, t ) f ( x, t ) x (x) 2 ... (5.38)
x 2! x 2
Taking d / dx as f and retaining terms only up to first order in Ax, we can write
2
x ... (5.39)
x x x x x x 2
so that terms within the square brackets in Eq. (5.37) can be rewritten as
2
2 x (5.40)
x x x x x x
265
2 F 2
t 2 l x 2
where l = m / x . You will observe that Eq. (5.41) is a partial differential equation. Moreover,
the quantity F / l has the dimensions of the square of velocity. For this reason, this equation is
referred to as the wave equation.
We thus find that longitudinal motion of a large number of coupled masses results in the
phenomenon of wave propagation. We obtain a similar equation for a system of large number of
coupled masses executing transverse oscillations. These will be discussed in detail later.
5.5 SUMMARY
The longitudinal oscillations of two (identical) coupled masses are nut simple harmonic.
The resultant amplitudes of coupled masses resemble a modulated pattern.
The displacement of either of the two (identical) coupled masses can be regarded as
superposition of normal modes of the system. Each mode represents .an independent
SHM. The normal mode frequencies are given by
k' k '2k
1 and 2
m m
where k and k ' are spring constants.
The normal mode frequencies of a system of two coupled pendulums are given by
g g 2k
1 and 2
l l m
266
(a) (b)
Fig. 5.9
2. Two equal masses (m) are connected to each other with the help of a spring of force
constant K and then the upper mass is connected to a rigid support by an identical spring
as shown in Fig. 5.10. The system is made to oscillate in the vertical direction. Show that
the two normal frequencies are given by
K
2 (3 5 )
2m
Fig. 5.10
3. Consider two capacitively coupled circuits shown in Fig. 5.11. Write down the equations
of motion for current and compute normal mode frequencies.
Fig. 5.11
267
5.7 SOLUTIONS
SAQ 1
1 2 1
x1 [a1 cos(1t 1 ) a 2 cos( 2 t 2 )] (i)
2 2
and
1 2 1
x2 [a1 cos(1t 1 ) a 2 cos( 2 t 2 )] (ii)
2 2
dx1 1
[a11 sin(1t 1 ) a 2 2 sin( 2 t 2 )] (iii)
dt 2
dx2 1
[a11 sin(1t 1 ) a 2 2 sin( 2 t 2 )] (iv)
dt 2
SAQ 2
At the displacements x1 and x 2 , the speeds of pendulums A and B will be dx1 / dt x1 and
dx2 / dt x 2 respectively. Therefore,
Kinetic energy K .E. (1 / 2)m( x1 ) 2 (1 / 2)m( x 2 ) 2 (i)
Potential energy P.E. = (1 / 2)m 02 ( x12 x22 ) (1 / 2)m s2 ( x1 x2 ) 2 (ii)
g k
Where 02 and s2
l m
From Eqs. (5.6a) and (5.6b), we have
1 2 1 2
x1 and x 2
2 2
Hence,
1 2 1 2
x1 and x 2
2 2
1 1
K .E. = m(1 2 ) 2 m(1 2 ) 2
8 8
m
= [(1 ) 2 (2 ) 2 ] (iv)
4
[ (a b) 2 (a b) 2 2a(a 2 b 2 ) ]
268
and
1 1
P.E. m 02 ( x12 x22 ) m 02 ( x12 x22 )
2 2
1 1
= m 02 (12 22 ) m s2 22
4 2
m 2 2
= ( 0 1 02 22 ) (v)
4
where 2 0 2 s2
2 2
SAQ 3
If normal coordinates are defined as
m m
1 ( x1 x2 ) and 2 ( x1 x2 )
2 2
we can write
1 1
x1 (1 2 ) and x2 (1 2 )
2m 2m
Hence, the total energy of the coupled pendulum in terms of normal coordinates is
1
E [(1 ) 2 (2 ) 2 1212 22 22 ]
2
1
x1 (1.50 0.50) 10 3 kg1 / 2 m 4.5 10 3 m
0.2kg
1
x2 (1.50 0.50) 10 3 kg1 / 2 m 2.2 10 3 m
0.2kg
SAQ 4
The equations of motion of the coupled pendulums A and B are
g ka2
x1 x1 2 ( x2 x1 )
l ml
and
g ka2
x2 x2 2 ( x2 x1 )
l ml
By comparing these equations with Eqs. (5.29a) and (5.29b), we get
g
12 02
l
and
269
g 2ka2
22 02 2 s2
l ml 2
SAQ 5
From Eqs. (5.33a) and (5.33b), we have
0
1
1
and
0
2
1
Here, 0 = 600 Hz and = 0.44 so that
600
1 Hz = 500 Hz
1.44
and
600
2 Hz = 802 Hz
0.56
Terminal Questions
1 K 0
1. 0 = = 2 Hz
2 m 2
k
02 4 2 02 4 2 2 2 16 2 (Hz)2 (i)
m
When A is clamped, the equation of motion of B will be given by
d 2 xB
m 2 ( K K ' ) x B
dt
2
d xB K K'
or 2
( ) xB 0
dt m
The frequency of this simple harmonic motion is given by
1 K K'
B = 2.5 Hz
2 m m
K K'
4 2 B2 4 2 (2.5) 2 (Hz)2 = 25 2 (Hz)2 (ii)
m m
Subtracting Eq. (i) from Eq. (ii), we get
K'
9 2 (Hz)2
m
Now, the angular frequencies of two normal modes of vibration are given by
K
12 4 2 12 02 16 2 (Hz)2
m
or 1 = 2 Hz
K 2K '
22 4 2 22 = (16 2 18 2 ) (Hz)2 = 34 2 (Hz)2
m m
270
17
2 Hz = 2.9 Hz
2
K K
2
m m
K 2K
2
m m
Hence,
3K 2 K 2
4 2 0
m m
3K 5K
2
2m 2m
271
K
Thus, 2 (3 5 ) for the slower mode.
2m
and
K
2 (3 5 ) for the faster mode.
2m
3. Refer to Fig. 5.11. The equations governing the balance of emf’s in two circuits are
di1 q q
L 1 3
dt C C
and
di2 q q
L 2 3
dt C C
Differentiating these equations with respect to time and using the relation i dq / dt , we
get
d 2 i1 1 1
L 2
i1 (i2 i1 )
dt C C
and
d 2 i1 1 1
L 2
i2 (i2 i1 )
dt C C
If we replace L by m , 1 / C by k ' k and i by x , then these equations become
identical to Eqs. (5.1) and (5.3). Hence, the two normal frequencies of the system are:
1
1
LC
and
3
2 .
LC
272
UNIT 6 WAVE MOTION
Structure
6.1 Introduction
Objectives
6.2 Basic Concepts of Wave Motion
Types of Waves
Propagation of Waves
Graphical Representation of Wave Motion
Relation between Phase Velocity, Frequency and Wavelength
6.3 Mathematical Description of Wave Motion
Phase and Phase Difference
Energy Transported by Progressive Waves
6.4 One-Dimensional Progressive Waves: Wave Equation
Waves on a Stretched String
Waves in a Fluid
Waves in a Uniform Rod
6.5 Wave Motion and Impedance
Impedance offered by Strings: Transverse Waves 1
Impedance offered by Gases: Sound Waves
6.6 Waves in Two and Three Dimensions
6.7 Summary
6.8 Terminal Questions
6.9 Solutions
6.1 INTRODUCTION
In Unit 5 you have learnt that when one mass in a system of N coupled masses is disturbed, the
disturbance is gradually felt by all other masses. You can think of many other similar situations in
which oscillations at one place are transmitted to some other place through the intervening
medium. When we talk, our vocal cord inside the throat vibrates. It causes air molecules to
vibrate and the effect-speech is transmitted. When it makes our ear drum to vibrate, it is heard.
Do you know what carries the audio information? The information is carried by a (sound) wave
which propagates through the medium (air). If you have ever stood at a seashore, you would need
no description of waves.
In addition to sound and water waves, other familiar types of waves are: ultrasound waves and
electromagnetic waves, which include visible light, radio waves, microwaves and x-rays. Matter
waves, shock waves, and seismic waves are other less familiar but important types of waves. You
will note that all our communications depend on transmission of signals through waves. The use
of x-rays in medical diagnosis is so very well known. Nowadays we also use ultrasound waves –
sound waves of frequency greater than 20 kHz – to make images of soft tissues in the interior of
humans. Sound waves are used in sound ranging; sonars and prospecting for mineral deposits and
oil – commodities governing the economy of nations these days. This means that understanding
of the physics of wave motion is of fundamental importance to us. In this unit we will confine to
mechanical waves with particular reference to sound waves.
When a progressive wave reaches the boundary of a finite medium or an interface between two
media, waves undergo reflection and/or refraction. These will be discussed in detail in the next
unit.
You would recall that our discussion of oscillations was simplified because of some basic
similarities between different physical systems. Once we understood the behaviour of a model
273
spring-mass system, we could easily draw analogies for others. Exactly the same simplification
occurs in the study of waves. The basic description of a wave and the parameters required to
quantify this description remain the same when we deal with a one-dimensional (1-D) wave
travelling along a string, a 2-D wave on the surface of a liquid or a 3-D sound wave. For this
reason, in this unit we shall first consider basic characteristics of wave motion. Then we would
calculate the energy transported by progressive waves. The vocabulary language and ideas
developed here will then be applied to waves on strings, liquids and gases.
Objectives
After going through this unit, you should be able to:
Fig. 6.1 The motion of a disturbed mass in the coupled spring-mass system. The disturbance is
eventually communicated to adjacent masses. This results in wave propagation. You would note
that regions of compression and elongation move along the system, which is shown here at two
different times.
A vivid demonstration of the energy carried by water waves is in damage caused in coastal areas
by tidal waves in stormy weather. Do you know that 3m high oceanic wave can lift 30 bags of
wheat by about 10 ft.?
274
Another important characteristic of mechanical waves is their velocity of propagation, referred to
as wave velocity. It is defined as the distance covered by a wave in unit time. It is different from
the particle velocity, i.e. the velocity with which the particles of the medium vibrate to-and-fro
about their respective equilibrium positions. Moreover, the wave velocity depends on the nature
of the medium in which a wave propagates. A wave has a characteristic amplitude, wavelength
and frequency. You must have learnt about these in your earlier classes. We will, however
recapitulate these in sub-section 6.2.4.
You can see with unaided eyes the actual propagation of a disturbance in water. Can you see a
sound wave propagating in air? Obviously, you cannot. Then you may like to know as to how we
detect sound waves. We observe the motion at the source (like a guitar string or drum membrane)
or at the receiver (microphone membrane). Another question that comes to our mind is: Are
sound waves and water waves similar? If not, how are waves classified? Let us now proceed to
know the answer to this question.
Waves set up by a single, isolated disturbance are called pulses. The dropping of a stone in still
water of a pond, the sound produced by clapping of hands, a single word of greeting or command
shouted from one person to another belong to this category. When an engine joins the
compartments, the jerk produces a disturbance which is carried through as a pulse. But
continuous and regular oscillations produce periodic waves. This, along with waveforms for
sound produced by a violin and a piano, is shown in Fig. 6.2. The simplest type of a periodic
wave is a harmonic wave.
Fig. 6.2 Waveform for (a) harmonic wave (b) the violin and (c) piano
275
When the motion of particles of the medium is perpendicular to the direction in which the wave
propagates, it is called a transverse wave. Waves on a string under tension are transverse, as in a
violin. You can generate transverse waves on a coupled spring-mass system of Fig. 6.1 by
displacing a mass at right angles to the spring as shown in Fig. 6.3a. (Electromagnetic waves are
also transverse in nature. But they do not require medium for propagation.)
Fig. 6.3 (a) A transverse wave and (b) a longitudinal wave on a coiled spring-mass system.
The broad arrow Indicates the direction of wave propagation.
When the motion of particles of the medium is along the direction in which wave propagates, the
wave is called a longitudinal wave. Sound waves in air are the most familiar example of
longitudinal waves. You can generate a longitudinal wave on a coupled spring-mass system of
Fig. 6.1 by displacing a mass along the length of the spring (Fig. 6.3b).
In your school you may have been told that water waves, produced by winds or otherwise, are
transverse and the motion is confined to the surface. But this is not correct. Strictly speaking, the
motion gradually extends with diminishing amplitude to deeper layers. Moreover, oscillations
have longitudinal as well as transverse components. That is, water waves are composite; partly
transverse and partly longitudinal. This is illustrated in Fig. 6.4. Similar waves can occur at the
surfaces of elastic solids. Such waves are called Rayleigh waves.
Fig. 6.4 Water waves are composite. The arrows at different points indicate instantaneous
velocity of water whereas dotted circles are paths traced out by parcels of water as the
wave passes. The direction of wave propagation is shown by the broad arrow.
In general, only longitudinal waves can propagate in gases and liquids but in solids both
transverse and longitudinal waves can propagate.
Earlier in this course, you have read about torsional oscillations. When such a disturbance
propagates in a medium, we have a torsional wave. In the following sections, you will learn about
mechanical waves in general and sound waves in particular. But before that you should answer
the following SAQ.
276
SAQ1
(i) The frequency of ultrasound wave is more than ..........
(ii) Water waves are .......... waves
(iii) Waves transfer.......... not ..........
(iv) Light waves require .......... medium
(v) Waves on sitar strings are ..........
Take a long elastic string and fix its one end to a distant wall. Hold the other end tightly. Move
your hand up and down. What do you observe? A disturbance travels along the string. This
disturbance is due to the up and down motion of the particles of the string about their respective
mean positions. When the motion of the arm (hand) is periodic, the disturbance on the string is a
wave with a sinusoidal profile. The shape of a portion of the string at intervals of T / 8 is shown
in Fig. 6.5. You will observe that the waveform moves to the right, as shown by the broad arrow.
You may ask: Why is the whole string not displaced simultaneously? The time lag between
different parts is due to gradual transfer of disturbance between successive particles.
Fig. 6.5 A periodic motion of the hand generates waves with a sinusoidal profile. In parts (b)-(j)
we have depicted wave profile on the string at intervals of T / 8
277
In this connection, it is important to distinguish between the motion of the waveform and the
motion of a particle of the string. While the waveform moves with a constant speed, the particles
of the string execute SHM. To illustrate this difference clearly, let us mark nine equidistant points
on the initial portion of the string. We assume that this string oscillates with a period T. Let us tie
one end of this string (at mark 1) to a vertically oscillating spring-mass system as shown in Fig.
6.6. As the mass m on the spring moves up and down, the particles at the marked positions begin
to oscillate one after the other. In time /; the disturbance initiated at the first particle will reach the
ninth particle. This means that in the interval T / 8 , the disturbance propagates from the particle at
mark 1 to the particle at mark 2. Similarly, in the next 778 interval, the disturbance travels from
the panicle at mark 2 to the particle at mark 3 and so on. In parts (a) - (i) of Fig. 6.7 we have
shown the instantaneous positions of particles at all nine marked positions at intervals of 778.
(The arrows indicate the directions of motion along with which particles at various marks are
about to move.) You will observe that
Fig. 6.6 A string fastened to an oscillating spring mass system: Illustration of the difference
between motion of the waveform and the motion of particles
The important point to note here is that while the wave moves along the string all panicles of the
string are oscillating up and down about their respective equilibrium positions with the same
period (D and amplitude (a). This wave remains progressive till it reaches the fixed end.
We would now like you to know to represent a wave graphically as well as mathematically. This
forms the subject of our discussion in the following sections.
278
Fig. 6.7 Instantaneous profiles at intervals of T / 8 when a transverse wave is
generated on a string.
t x
y (t ) a sin 2 and y ( x) a sin 2
T
The argument of the sine function ensures that the function repeats itself regularly.
We can draw an analogy between the wavelength and the period. The wavelength is separation in
space between successive in-phase points on the wave. On the other hand, period is separation in
time between equivalent instants in successive cycles of vibration. This means that the
wavelength and the period are respectively the spatial and the temporal properties of a wave.
279
Fig. 6.8 (a) The profile of vibrations at a given position in the path of a wave, and (b) the profile
of a wave at a particular Instant. It is snapshot of the wave travelling along the string.
It is important to note that the scales for y (x) and y (t ) are different, For sound waves, the
displacement amplitudes are a small fraction of 1 mm whereas x extends to several metres.
Human ears can hear sound of 1000 Hz quite clearly. The amplitude of the wave corresponding to
the faintest sound that a normal human ear can hear is approximately 10 11 m. This is smaller
than the radius of the atom (= 10 10 m). Yet our ears respond to such a small displacement!
Another point about graphical representation is that it can be used for both transverse and
longitudinal waves.
In longitudinal waves, the displacement of particles is along the direction of wave propagation. In
Fig. 6.9, the hollow circles represent the mean positions of equidistant particles in a medium. The
arrows show their (rather magnified) longitudinal displacements at a given time. You will observe
that the arrows are neither equal in length nor in the same direction. This is clear from the
positions of solid circles, which describe instantaneous positions of the particles corresponding to
the heads of the arrows. The displacements to the right are shown in the graph towards + y axis
and displacements to the left towards the – y axis.
For every arrow directed to the right, we draw a proportionate line upward. Similarly, for every
arrow directed to the left, a proportionate line is drawn downward. On drawing a smooth curve
through the heads of these lines, we find that the graph resembles the displacement-time curve for
a transverse wave. If we look at the solid circles, we note that around the positions A and B, the
280
particles have crowded together while around the position C, they have separated farther. These
represent regions of compression and rarefaction. That is, there are alternate regions where
density (pressure) are higher and lower than average. A sound wave propagating in air is very
similar to the longitudinal waves that you can generate on your string/spring. This similarity is
clearly illustrated in Fig. 6.10. A sound wave may be considered either as a pressure wave or as a
displacement wave. However, the pressure wave is 90° put of phase with the displacement wave.
That is, when displacement from equilibrium at a point is maximum, the excess pressure (over the
normal) is zero and vice versa. The variations of pressure and density are represented graphically
in. Fig. 6.11. This means that in longitudinal waves, alternate high and low pressures propagate
along the wave.
Fig. 6.11 The sound waves can be viewed in terms of changes in pressure or density
wavelength
v (6.3)
period T
v f
That is, the speed of any wave is equal to the product of its frequency and the wavelength. This
equation predicts that in a given medium, the speed of a wave of given frequency is constant,
This is a very important relation.
You will note that we have derived Eq. (6.4) with reference to a transverse wave in a string. But it
holds for all other media like air, water, glass etc. as well as longitudinal waves. At STP, the
speeds of sound waves in air, water and steel are 332 m s 1 , 1500 m s 1 and 5100 m s 1 ,
respectively. (This explains why the whistle of an approaching train may be heard twice - first as
the sound travels through the railroad track and again as it (ravels through the air.) Ripples on the
surface of a pond move with a speed of about 0.2 m s 1 . The seismic waves move with a speed of
the order of 6 10 3 m s 1 in the earth's outer crust and light moves with a speed of 3 108 m s 1 .
That is why light that originates on or near the earth reaches us almost instantly.
281
6.3 MATHEMATICAL DESCRIPTION OF WAVE MOTION
In the preceding section you have learnt that at a particular time, a wave is described by Eq. (6.1).
As time passes, the wave propagates along the + x direction. So at a given value of x, the
displacement of medium particles must change with time. This information is not contained in
Eq. (6.1). This means that it is not a complete equation for the wave. You would like to know as
to how we can modify Eq. (6.1). To answer this question, let us consider Fig. 6.12, which shows a
'snapshot' of a wave moving with speed v along x -axis. Now imagine two particles, say at A
( x = 0) and at B separated by a distance x. You can easily convince yourself that a disturbance
created at A will reach B in time t / v . This means that the particle at B will have the same
displacement as particle at A at time t ' t x / v . Mathematically, we can express this as
We can obtain y( x 0) on replacing t in Eq. (6.2) by t ' t x / v . Then it readily follows that
when a wave moves along + x direction with speed v , the displacement of medium particles as a
function of x and t is given by the equation
2 x
y ( x, t ) a sin t (6.6a)
T v
or
2
y ( x, t ) a sin (vt x) (6.6b)
since v / T . We can also rewrite this equation as
t x
y ( x, t ) a sin 2 (6.6c)
T
At t = 0 (as also when t/T is an integer), this equation reduces to Eq. (6.1). [The negative sign
implies slight re-adjustment in the phase of the wave at t = 0.] Further, you can check that this
equation displays periodicity by calculating the displacement at a distance x + from the origin.
To ensure that you have grasped these ideas we would like you to solve the following SAQ.
SAQ3
A 275 Hz sound wave travels with a speed of 340 ms along the x-axis. Each point of the medium
moves up and down through 5.0 mm. Write down the equation for the wave. Calculate (i) the
wavelength of the wave, and (ii) velocity and acceleration of medium particles.
282
Eqs. (6.6a), (6.6b) and (6.6c) show that if we sit at a fixed position x = 0 say, then the
displacement varies sinusoidally in time: y( x 0, t ) a sin(2 t / T ) . We know that this
relation expresses SHM with angular frequency 0 2 / T .
Another equivalent and convenient form for y( x, t ) is written in terms of the wave number, k
defined as the number of radians of wave cycle per unit distance:
2
k
(You should not confuse this k with spring constant used earlier.)
You will recognize that the wave number is the spatial analog of the angular frequency. In
terms of 0 and k , we can also write the wave equation as
The simple way in which k and 0 enter this description of the wave explains why these
quantities are so often used.
On comparing Eqs. (6.6a) and (6.6d) you will observe that v and k are connected through the
relation
0
v (6.7)
k
Eqs. (6.6a - d) provide us equivalent description of a travelling wave moving in the + x direction.
But the choice of the particular form to be used for a specific problem is a matter of convenience.
How would you describe a wave propagating in the negative x -direction? You can easily
convince yourself that to describe a wave moving in the negative x -direction, we should replace
x by x in Eq. (6.6).
SAQ4
A sinusoidal water wave having a maximum height of 7.4 cm above the equilibrium water level is
propagating in the x direction with a speed of 93 cm s 1 . The distance between two successive
crests is 55 cm. Write the wave equation in terms of angular frequency and wave number. Also
calculate the particle velocity.
0 t kx (6.8)
283
You will note that the phase changes both with time and the space coordinate. With time, it
changes according to
The minus sign in this equation signifies that in a wave moving along x direction, the forward
points lag in phase. That is, they reach the successive stages of vibration later.
d 0 dt kdx
dx
vp 0 (6.10)
dt k
On comparing Eqs. (6.7) and (6.10) you will observe that the expression deduced earlier is
actually the phase velocity.
Let us consider a thin layer of thickness and cross-sectional area A at a distance x from the source.
If is the density of the medium, the mass of the layer is x A . Therefore, kinetic energy of
the layer
y ( x, t )
2
1 1
K .E.( x, t ) = mv 2 xA
2 2 t
1
= xA 02 a 2 cos 2 ( 0 t kx)
2
284
= 2 2 f 02 Axa 2 cos 2 ( 0 t kx) (6.11)
1
This expression implies that kinetic energy oscillates between zero and xA 02 a 2 . This is
2
because the value of the function cos 2 ( 0 t kx) varies between 0 and 1.
1
Over one full cycle, the average value of cos 2 is . So, the average kinetic energy over a time
2
period is
1
K .E. 02 a 2 Ax 2 a 2 f 02 A x (6.12)
4
What about the potential energy? From Unit 1 you would recall that in SHM, the average kinetic
energy and average potential energy are equal. Is the same true for a harmonic wave as well?
Physically, we except so. Let us now compute potential energy analytically.
2 y( x, t )
F Ax
t 2
= 4 2 f 02 xAy ( x, t )
We know that the work done by this force, when the layer of interest is displaced through y from
its equilibrium position, is stored in the layer as its potential energy. So we can write
A4
y
U ( x, t ) = 2
f 02 xy ' dy '
0
= 2 2 f 02 xAy 2
= 2 2 f 02 Axa 2 sin 2 ( 0 t kx) (6.13)
The minus sign tells us that the work is done on the layer. (This is of no consequence when we
calculate total energy of the wave.) The time-averaged potential energy of the wave is
U 2 f 02 Aa 2 x K .E. (6.14)
On combining Eqs. (6.11) and (6.13), we find that the total energy of the wave is
E K.E. U
= 2 2 a 2 f 02 Ax
= K .E. U (6.15)
This shows that half the energy of the wave is kinetic and the other half is potential.
285
What happens to this energy? As the layer moves, it pushes the next layer. In the process, it
transmits its energy. Now you may like to know: How long does this layer take to give up its
energy? Or, what is the average rate of energy flow? To calculate this, we note that if the wave is
moving at speed v the energy passes the layer in time t x / v . Hence
E 2 2 a 2 f 02 xA
P
t x / v
= 2 a f 0 vA
2 2 2
(6.16a)
This shows that average rate of energy flow, or what we call power, is proportional to wave
Speed and the square of the amplitude.
I 2 2 a 2 f 02 v (1 / 2)vVA2 (6.16b)
where V A2 2 f 0 a .
The SI units of intensity are Jm 2 s 1 or Wm 2 . From Unit 1, you will recall that total energy is
proportional to the square of the amplitude of oscillation. In the same way, the intensity of a wave
at a given position is proportional to the square of amplitude at that position. For a second wave,
we can write
I p02 (6.18 a)
where p 0 is the maximum change in pressure over normal pressure. Note that when we express
intensity in terms of p 0 , the frequency does not appear explicitly in the expression. This means
that 100 Hz sound wave has the same intensity as a 10 kHz sound wave; both have the same
amplitude.
How does the intensity of a wave at a point vary with distance from the source?
We know that the area crossed by a wave increases as it spreads out. If it originates from a point
source or the distance from the source is much greater than the size of the source, the area will be
286
almost spherical ( r 2 ) . Then principle of conservation of energy demands that E 4Ir 2 be
constant. So as r increases, intensity decreases as 1 / r 2 :
1
I (6.18b)
r2
1
p0
r
Since a p0 , this relation implies that amplitude of a wave is inversely proportional to the
distance from the source. This explains why we can be heard up to a certain distance. (Beyond
this the amplitude becomes too small to affect our eardrums.) You must note that these results
will hold if the wave is not absorbed or obstructed. In Table 6.1 we have listed intensities of
waves generated from different sources.
Source/Wave Intensity
Sound
Threshold of hearing 10-12
Rustle of leaves 10-11
Whisper, intensity at eardrum 10-10
Ordinary conversation 3.2 10-6
Street traffic 10-5
Bursting cracker, at 1 m 8 10-5
Jet Taking off, at 30m 5
Electromagnetic Waves
Radio in home 10-8
TV signal, 5.0 km from 50 kW transmitter 1.6 10-4
Sunlight intensity at earth's orbit 1368
1m from typical camera flash 4000
Inside microwave oven 6000
Target of laser fusion experiment 1018
Seismic wave
5 km from Richter 7.0 quake 4 104
SAQ5
At a distance of 1m from a bursting cracker, the intensity of sound is 8 10 5 Wm 2 . The
threshold of human hearing is about 10 12Wm 2 . If sound waves spread out evenly in all
directions, how far from the source could such a sound be heard?
Our ear is sensitive to an extremely large range of intensities. So we can define a logarithmic
intensity scale. The intensity level of a sound wave is defined by the equation
287
10 ln( I / I 0 )
where I 0 (= 10 12Wm 2 ) denotes the threshold of hearing. Intensity levels are expressed in
decibels, abbreviated db. Our ear can tolerate intensity up to 120 db.
Fig. 6.13: A magnified view of a small element of the stretched string. The net force acting on it
is non-zero. The vertical distortion is exaggerated (for clarity)
SAQ6
Using dimensional analysis, show that
v K F /m
We will now analyse the problem mechanically by considering a small element along the string.
288
Suppose that the string is distorted slightly so that the magnitude of tension on the string
essentially remains unchanged. Fig. 6.13 shows (magnified view of) a small element of the
distorted string. You will observe that the direction of the (tension) force varies along the element
under consideration. Why? This is because the string is curved. This means that the tension forces
pulling at opposite ends of the element, although of the same magnitude, do not exactly cancel
out. To calculate the net force along x and y-axis, we resolve F in rectangular components. The
difference in the x and y components of tension between the right and the left ends of the element
is respectively given by
Fx F cos 2 F cos 1
and Fy F sin 2 F sin 1
Where 1 and 2 are angles which the tangents drawn at the ends make with the horizontal. For
small oscillations ( 4°).
cos 1 cos 2 1
This means that there is no net force in the x -direction, Fx = 0, and the string will be very nearly
horizontal. This implies that the sines of the angles are very nearly the same as their tangents, i.e.
sin 1 tan 1
and sin 2 tan 2
But the tangents of the two angles are just the slopes of the string - or the derivatives dy / dt at
the ends of the element under consideration. Then, the -component of the force on the element is
approximately
Fy F (tan 2 tan1 )
dy ( x, t ) dy ( x, t )
= F (6.20)
dx x x dx x
From the previous unit you may recall that the quantity in parentheses is just the change in the
first derivative from one end of the interval x to the other. Dividing that change by x gives, in
the limit x 0 , the rate of change of the first derivative. But we know that the displacement of
the string is a function of position as well as time. If either of these variables changes, the
displacement also changes. You will recognise that Eq. (6.20) is valid for the configuration of the
string at a particular instant of time. Therefore, the derivative in this equation is to be taken with
the time fixed. We call a derivative taken with respect to one variable while other(s) is (are) kept
constant a partial derivative. We denote partial derivatives with the symbol in place of the
usual symbol 'd'. Then, Eq. (6.20) can be rewritten as
289
2 y ( x, t )
Fy F x
x 2
This equation gives the net force on the segment x . By Newton's second law of motion, we can
equate this force to the product of mass and acceleration of the segment. The mass of the segment
of length x is mx . Then, we can write
2 y( x, t ) 2 y ( x, t )
mx F x
t 2 t 2
2 y( x, t ) m 2 y( x, t )
(6.21)
x 2 F t 2
Eq. (6.21) has been obtained by applying Newton's second law to a small element of a stretched
string. Since there is nothing special about this particular element on the string, this equation
applies to the whole of it.
Let us pause for a minute and ask: What goals we set for ourselves and how Eq. (6.21) helps us in
attaining them? We wish to know what determines the speed of a wave. To know this let us
assume that a harmonic wave described by
y( x, t ) a sin( 0 t kx)
moves on the string. If this mathematical form is consistent with Newton's law, then you can be
sure that such waves can move on the string. To see this, you should calculate the second partial
derivatives of the particle displacement:
2 y
k 2 a sin( 0 t kx)
x 2
and
2 y
02 a sin( 0 t kx)
t 2
m 2
k2 0
F
or
0
2
F
k m
290
What is implied by this equality? We know that it has followed from Newton's law of motion
applied to a stretched string when a harmonic wave is travelling along it. So, the above relation
tells us that only those waves can propagate on the string for which wave properties 0 and k are
related to F and m through the relation
0 F
k m
0 F
v (6.22)
k m
This relation tells us that the velocity of a (transverse) wave on a stretched string depends on
tension as well as mass per unit length of the string, not on wavelength or time period. This
means that v is not a property of the material of the string. It involves an external factor – tension
– which can be adjusted for fine-tuning. This explains why musicians are seen adjusting tension
in their stringed instruments. However, no such thing is done in case of a flute or a harmonium.
2 y ( x, t ) 1 2 y
2 2 (6.23)
x 2 v t
This result is identical to Eq. (5.41) and expresses one dimensional wave equation. It holds so
long as the oscillations of the string have small amplitude. You may now ask: Will Eq. (6.23)
hold for large amplitude disturbances as well? The answer to this question is that large amplitude
oscillations result in a more complicated equation and the wave speed tend to depend on
wavelength as well. Before you proceed further, you may like to solve the following
SAQ 7
A 1 m long string having mass 1 g is sketched with a force of 10 N. Calculate the speed of
transverse waves.
We now know that the speed of a wave is determined by the interplay of elasticity and inertia of
the medium. Elasticity gives rise to the restoring force and inertia tells us how the medium
responds to them. Since a fluid (gas or liquid) lacks rigidity, transverse waves can propagate only
in solids. However, longitudinal waves can propagate in all phases of matter – plasmas, gases,
liquids and solids – in the form of condensations and rarefactions. We will now consider wave
propagation in a fluid.
291
tuning fork at its one end or displacing the fluid to the right using a piston. As the wave passes
through the column under consideration, its pressure, density and volume change. Let us assume
that in time t planes PQ and SR move to P'Q' and S'R', respectively, as shown in Fig. 6.14b. If
the planes PQ and SR are displaced through (x) and ( x x) the change in thickness, l , is
l ( x x) ( x) x
x
In writing this expression we have used Taylor series expansion for ( x x) about (x) . This
means that the change in volume V is
V Al Ax
x
Fig, 6.14: (a) Equilibrium state of the column PQRS of a fluid contained in a long tube of cross
sectional area A. (b) Displaced position of column under pressure difference
The volume strain, which is defined as the change in volume per unit volume, is given by
Ax
V
x (6.24)
V Ax x
The minus sign signifies that the column is compressed. This happens because the pressure on its
two sides does not balance. Let the difference of pressures be p( x x) p( x) .
The net force acting on the column is, therefore A[ p( x x) p( x)] . Using Taylor series and
retaining only first order term in x , we can write
p( x)
F= A x
x
292
( p0 p)
= A x
x
(p)
= A x
x
where p 0 is the equilibrium pressure and p is the excess pressure due to the wave.
Hence, the equation of motion for the column under consideration, according to Newton's second
law, is
2 (p)
xA Ax (6.25)
t 2
x
To express this result in a familiar form, we note that p and E , the bulk modulus of elasticity,
are connected by the relation
stress p
E
volume strain V / V
The negative sign is included to account for the fact that when pressure increases, volume
decreases. This ensures that E is positive.
We can rewrite it as
V
p E
V
p E
x
2
2 = E
x x x
2
= E
x 2
or
2 E 2 ( x, t )
(6.26)
t 2 x 2
If we identify
v E/ (6.26a)
293
as the speed of longitudinal waves, Eq. (6.26) becomes identical to Eq. (5.41).
You will note that the wave speed is determined only by F and ; two properties of the
medium through which the wave is propagating. Let us now consider the propagation of sound
waves in a gas.
v p0 / (6.27)
For air at STP, = 1.29 kg m 3 and p 0 = 1.01 105 N m 2 . Then the velocity of sound in air
using Newton's formula comes out to be
1.01 10 5 Nm 2
v = 280 m
1.29 kg m 3
But experiments show that the velocity of sound in air at STP is 332 m s 1 . This gives rise to an
interesting question: How could Newton come so close to the correct answer and yet miss it by
about 15%? It means that something is wrong with his derivation. You may now ask: How to
explain the discrepancy? The problem lies with the use of Boyle's law, which holds only at
constant temperature. The discrepancy was solved when Laplace pointed out that when sound
waves move in a medium, the particles oscillate very rapidly. In the process, regions of
compression are heated up while the regions of rarefaction get cooled. That is, local changes in
temperature do occur when sound propagates in air. This effect produces a larger phase velocity.
However, the total energy of the system is conserved. This means that adiabatic changes occur in
air when sound propagates.
For an adiabatic change, E is times the pressure, where is the ratio of the specific heats at
constant pressure and at constant volume, i.e. E p0 . Then Eq. (6.27) becomes
p0
v (6.28)
For air, = 1.4. So the velocity of sound in air at STP works out to be 331 m s 1 , which is an
excellent agreement with the measured value. This shows that Laplace's agreement is correct.
At a given temperature, p0 / is constant for a gas. So Eq, (6.28) shows that the velocity of a
longitudinal wave is independent of pressure.
294
A process is said to be isothermal if temperature remains constant during the
process. In an adiabatic process, the total energy of the system remains constant.
For an isothermal process, Boyle’s law states that pV = constant. Any change in p
p
and/or V is related as Vp pV 0 or E p.
V / V
You will now like to know why the heat does not have time to flow from a compression to a
rarefaction and equalize the temperature everywhere. To attain this condition, heat has to flow a
distance of one-half wavelength in a time much shorter than one-half of the period of oscillation.
This means that we would need
Since heat flow is mostly due to conduction, the speed of air molecules is given by
kHT
v rms (6.30)
M
where M is the mass of air molecules and T is the absolute temperature. We can similarly write
kHT
vsound (6.31)
M
Thus, even if air molecules travel a distance of / 2 , they will not be able to transfer heat in
time. In practice they move randomly in zig-zag paths of the order of 10 5 cm and as long as
105 cm , the adiabatic flow is a very good approximation. The shortest wavelength for
audible sound (1.6cm) corresponds to 20 kHz.
The ability to measure the speed of sound has been put to many uses in defence. During World
War I, a technique called sound ranging was developed to locate the position of enemy guns by
using the sound of cannon in action.
295
water than in air. This has led to the development of Sonar. High frequency sound waves are used
in Sonar which can measure the depth of sea bed, detect submarines and enemy torpedoes.
stress p
Y
longitudin al strain l / l
2 Y 2
(6.32a)
t 2 x 2
which is the equation of wave motion in the rod. The velocity of longitudinal waves is given
Y
v (6.32b)
SAQ8
For a steel rod, Y 2 1011 Nm 2 and 7800 kgm3 . Compute the speed of sound. On
working out this SAQ you will find that v 5 10 3 ms 1 , which shows that longitudinal waves
travel faster in a solid than in a gas or a liquid. This means that you can know about a coming
train by putting your ear on the rails. However, you are advised never to do so!
F Zv
296
The constant of proportionality Z is called the wave impedance. From this equation it is clear
that impedance is numerically equal to the driving force which imparts unit velocity to a particle.
We will now consider some specific examples.
Let T be the tension in the string. The vertical component of tension (T ) along the negative y-
direction is equal to the applied transverse force (to give zero resultant force at the x = 0 end
of the string):
F0 cos 0 t T sin
The tangent (or slope) is defined at the x = 0 end of the string. Using Eq. (6.6d) we can relate
(dy / dx) to (dy / dt ) :
dy k dy
dx 0 dt
Tk
F0 cos 0 t a 0 cos 0 t
0
297
T
= a 0 cos 0 t
v
where v 0 / k .
Writing a 0 = v 0 as the velocity amplitude of the wave, the above equation reduces to
Tv0
F0 cos 0 t cos 0 t
v
giving
Tv0
F0
v
F0 T
(6.34)
v0 v
This result specifies the ratio of the amplitude of the applied force to the amplitude of particle
velocity for transverse waves in terms of the tension in the string and particle velocity. This result
can be used to get an expression for the characteristic impedance (Z) of the string which is
defined as:
F0 T
Z (6.35)
v0 v
This result shows that the characteristic impedance has units of Nm 1 s but its dimensions are
MT 1 .
From Eq. (6.26a) we recall that v T / m where m is the mass per unit length of the string.
T
Z Tm (6.36a)
v
v2m
Z mv (6.36b)
v
298
From Eq. (6.36a) we find that the characteristic impedance is governed by the mass per unit
length of the string and the tension in it. This means that a sonometer wire will offer different
impedance when it is loaded by different weights. Eq. (6.36b) tells us that since Z is related to the
velocity of the wave, it depends on the inertia as well as the elasticity of the medium.
SAQ9
Calculate the characteristic impedance offered by a thin wire of steel stretched by a force of
80 N. It weights 2 g per metre.
It means that dimensionally Z is the ratio of force per unit area to velocity.
The excess pressure p experienced by the medium when a longitudinal wave propagates
through it is given by
p E (6.38)
x
where E is the Bulk modulus of elasticity of the medium. This means that to know Z, we must
compute and . To do so, we recall that the particle displacement for a longitudinal wave
t x
travelling in the +ve x -direction is written as
2
( x, t ) a sin
(vt x)
2 2
a cos (vt x) (6.39a)
x
and
2 v 2
a cos (vt x) (6.39b)
t
2 2
p Ea cos (vt x) (6.40)
299
On substituting for p and / t from Eqs. (6.40) and (6.39b) in Eq. (6.37) we find that the
acoustic impedance is given by
where v is the wave velocity. This result shows that the units of acoustic impedance are Nm 3 s
and the dimensions are ML2T 1 (You should verify these before proceeding further.)
E
v
where is the density of the medium. Hence the acoustic impedance Z can also be expressed as:
E
Z E v (6.42)
v
This result shows that the acoustic impedance Z is given by the product of the density of the
medium and the wave velocity. This means that the denser the medium, the greater will be the
impedance offered. Yet we know that sound moves faster in solids than gases.
In the next unit, you will apply these results to compute reflection and transmission amplitude and
energy coefficients for a wave incident on a boundary separating two media.
SAQ10
Calculate the acoustic impedance of air at standard temperature and pressure. Use =
1.29 kg m 3 and v = 332 m s 1 . Will this value be more for air or water? Justify your answer.
We will not go into mathematical details to answer these questions. However, from physical
consideration, extension of Eq. (6.23) for 2-D wave is a straightforward exercise. Since forces
300
along x and y axes act independently, each one will contribute analogous term to the wave
equation so that Eq. (6.23) modifies to
2 ( x, y, t ) F 2 2
( x, y, t )
x 2 y 2
(6.43)
t 2
Let us pause for a minute and ask: Do sound and light waves emanate radially from a small two-
dimensional source? How can we describe seismic waves or a wave propagating in an elastic
solid? These are three-dimensional waves. To analyse 3-D waves we have to extend the
preceding arguments. The result is a 3-D wave equation.
SAQ 11
Generalize Eq. (6.43) in three dimensions.
6.7 SUMMARY
Mechanical (elastic) waves can be transverse as well as longitudinal. In a transverse
wave, particles of the medium vibrate normal to the direction in which a wave moves, in
a longitudinal wave vibrations of the particles of the medium are along the direction of
wave propagation.
The wave velocity, frequency and wavelength are connected by the relation: v f .
0
We can also express v as ratio of the angular frequency and wave number: v .
k
t x
y ( x, t ) a sin 2 =
T
2
= a sin (vt x)
= a sin( 0 t kx)
where T is the time period. The phase of a wave, = 0 t kx varies both with time
and space.
Waves carry energy. The total energy carried by a wave is half kinetic and half potential:
E 2 2 a 2 f 02 Ax
= K .E. U
The average rate of energy flow or average power is proportional to the wave speed and
to the square of the wave amplitude:
P 2a 2 f 02 A
301
For waves propagating in space, intensity is a more useful measure of energy carried by
waves. The intensity of a plane wave remains constant as the wave propagates. But for
spherical waves, the intensity decreases as the inverse square of the distance from the
source.
A wave propagating along a string (1-D wave) is described by the equation
2 ( x, t ) 2 ( x, t )
2
v
t 2 x 2
where k r k x x k y y and k 2 k x2 k y2 .
302
2. In normal conversation, the intensity (energy flux/area) of sound is 5 10 6 Wm 2 . The
frequency of normal human voice is about 1000 Hz. Compute the amplitude of sound
waves, given that the density of air at STP is 1.29 kg m 3 . Take the velocity of sound at
STP as 332 m s 1 .
6.9 SOLUTIONS
SAQs
1. (i) 20 kHz (ii) composite (iii) no (v) transverse
303
The angular frequency is connected to v and by the relation
2 v
0
Here v = 93 cm s 1 and = 55 cm. Hence,
2 93 cm s 1
0
55 cm
372 1
= s
55
Similarly, the wave number is related to wavelength as
2 2 22 4
k cm 1
7 55 cm 35
Hence,
372 4
y( x, t ) (7.4 cm) sin t x
35 35
1
= 78.6 cm s cos(10.6t 0.1x)
5. I 8 10 5 Wm 2
I 0 10 12Wm 2
The intensity falls off as the inverse square of the distance. If r denotes the distance at
which this sound could just be heard, then
8 10 5 Wm 2 r2
10 12Wm 2 (1 m) 2
r 9 10 3 m = 9 km
But this is not observed in practice. It is because of absorption of energy by the medium.
6. v kf ( F , m)
[v] [ K ][ F ]a [m]b
[ LT 1 ] [MLT 2 ] A [ML1 ]b
or
[ LT 1 ] [M a b ][ Lab ][T 2a ]
On comparing the powers of T , we get
1 2a
or
a 1/ 2
On comparing the powers of L , we have
a b 1
or
b a 1
= 1/ 2
Hence,
v k F /m
304
7. We know that the wave velocity on a stretched string is given by
v F /m
On substituting the given values, we get
10 N
v
10 kgm1
3
= 100 m s 1
From Eq. (6.32b),
v /
On inserting the given values, you will get
2 1011 Nm 2
v
7800 kg m 3
= 5 10 3 ms 1
TQs
1. We know that a harmonic wave in 1-D is described by
2
( x, t ) a sin
(vt x)
where a = 1 cm
(a) At x 10 cm , ( x, t ) 0.5 cm
305
2
0.5 cm (1cm) sin (vt 10)
or
2
sin (vt 10) 1 / 2 sin
6
This equality implies that
2 7
(vt 10)
6
or
7
vt 10 (i)
12
(b) At x 20 cm , ( x, t ) 0.5 cm
2
0.5 cm (1cm) sin (vt 20) sin
6
so that
vt 20 (ii)
12
From (i) and (ii), we get
10
2
or
20 cm
= 0.2 m
We know that
v f
v (500 Hz) (0.2 m) 100 ms 1
Hence,
2
( x, t ) (0.01cm) sin
(100t x) (0.01 m) sin[10 (100t x)]
0.2
306
3. Here we use the expression
p p0
v1 ; v0
0
and the gas equation
p0V0 pV
T0 T
307
UNIT 7 WAVES AT THE BOUNDARY OF TWO MEDIA
Structure
7.1 Introduction
Objectives
7.2 The Concept of Wavefront and Huygens' Construction
Reflection of Waves
Refraction of Waves
7.3 Reflection and Transmission Amplitude Coefficients
Transverse Wives
Longitudinal Waves
7.4 Reflection and Transmission Energy Coefficients
7.5 The Doppler Effect
Source in Motion and Observer Stationary
Source Stationary and Observer in Motion
Source and Observer both in Motion
7.6 Shock Waves
7.7 Summary
7.8 Terminal Questions
7.9 Solutions
7.1 INTRODUCTION
In Unit 6 we discussed the basic characteristics of wave motion. The propagation of waves on
strings and in fluids was discussed with particular reference to sound. You may now ask: What
happens to a wave when it encounters a rigid barrier, as for instance, in the case of a string whose
one end is tied to a rigid wall. The wave energy will not flow into the wall. But the wave cannot
stop there. Then where will its energy go? What happens is that the wave turns around and
bounces back along the string. We say that the wave has been reflected.
You must have experienced sound reflection in the form of echoes in large halls or in the
neighbourhood of hills. You must have also observed reflection of water (sea) waves from a fixed
barrier (sea shore). In the case of light, reflection from silvered surfaces, say in a looking mirror,
is the most common optical effect we know. The reflection of ultrasonic (sound) waves forms the
operating principle of sonars in depth-ranging, navigation, prospecting for oil and mineral
deposits. The reflection of e.m. waves governs the working of a radar for detection of aircrafts.
Reflection of radiowaves by the ionosphere makes signal transmission from one place to another
possible and is so crucial in the area of communications.
You may now like to know as to what would happen to the incident wave. You would agree that
the boundary is not very rigid and properties of the medium change suddenly. Now suppose that
we connect two strings of different mass per unit lengths. We observe that in such a case energy
is partly transmitted into the second string and the rest is reflected back along the first. The
phenomenon of partial reflection and transmission at a junction of strings has its analog in the
behaviour of all waves at interfaces between two different media. Shallow water waves are
partially reflected if water depth changes suddenly. Light incident on our atmosphere undergoes
partial reflection because of changes in the density of the medium. Partial reflection of ultrasound
waves at the interfaces of body tissues with different densities makes ultrasound a valuable
diagnostic tool.
Does this mean that waves never undergo complete refraction? Were this true, we could not
explain the working of lenses, which is fundamental to seeing and our contact with the
308
surroundings. You may have seen the sun before actual sunrise and after actual sunset. This is
because of refraction of light in the atmosphere.
In Sections 7.2 and 7.3 you will learn, using the concepts of Huygens' construction and the
concept of impedance, that when a wave is incident at a boundary separating two media, its
wavelength changes but frequency remains constant. But there are many situations where the
frequency of a wave also undergoes a change. This effect is known as Doppler Effect. You will
learn it in Section 7.5.
define a wavefront
construct the wave front for a given source
explain reflection and refraction of waves using Huygens' construction
compute the reflection and transmission amplitude coefficients
compute reflection and transmission energy coefficients
compute the apparent frequency of sound when the source and/or the observer (listener)
are in motion.
i) Each point on a wavefront becomes a fresh source of secondary wavelets, which move
out in all directions with the speed of the wave in that medium.
ii) The new wavefront, at any later time, is given by the forward envelope of the secondary
wavelets at that time.
iii) In an isotropic medium, the energy carried by waves is transmitted equally in all
directions.
If S is the source of sound or light (Fig. 7.la), then after an interval of time t , all particles of the
medium lying on the surface AB vibrate in the same phase. This is because all particles on the
surface AB are equidistant from the source. Any disturbance emanating from S is handed on to
them at the same time.
According to Huygen's construction, surface AB is called a primary wavefront. Each point on AB,
like a, b, c, etc., acts as secondary source (derived from the original source S). These secondary
sources give out waves (or disturbances) in all directions as demonstrated by drawing circles
around the points a, b, c, etc. The envelope of all these waves (which acts as a tangent to all of
them at any given instant), like the one at CD, forms another wavefront, called the secondary
wavefront. This, in short, means that the source S gives out wavelets in all directions. The
envelope of these wavelets acts as a primary wavefront. Each point on this primary wavefront
acts as a source for secondary wavelets. An envelope of these secondary wavelets forms a
309
secondary wavefront. Each point on this secondary wavefront gives out further wavelets to form
further secondary wavefronts. This process goes on and the wave keeps on spreading in space.
The direction SP (Fig. 7.la) in which the disturbance (originating at S) propagates is called a ray.
A ray is always normal to the expanding wavefront.
To visualise the Huygens' construction in space, you may imagine a point source to be at the
centre of a hollow sphere. The outer surface of this sphere then acts as a primary wavefront. If
this sphere is further enclosed by another hollow sphere of larger radius, the outer surface of the
second hollow sphere will then act as a secondary wavefront. If this sphere is further enclosed by
another sphere of still bigger radius, the surface of the outermost sphere becomes the secondary
wavefront. For this, the surface of the inner sphere acts as the primary wavefront. In two
dimensions, the primary and secondary wavefronts appear to be concentric circles, the parts of
which are shown in Figs. 7.la and 7.1 b.
Fig. 7.1 (a) Construction of Huygens' wavefront, (b) Depletion of • Secondary source
The formation of secondary sources as visualised by Huygens can also be understood pictorially
through a simple diagram. If we place a screen XY with a tiny hole at S' in the path of waves
emanating from the source S, S' acts as a secondary source (Fig. 7.1b). This gives out waves on
the other side of the screen. These waves spread out from S' as if S' is an original source itself.
In your school classes you have studied reflection and refraction of waves. We observe these
whenever a wave travelling in one medium, say air, meets the boundary of another medium.
Suppose we clamp one end of a string to a rigid wall and generate a pulse by moving the other
end. You will observe that the pulse is reflected at the fixed end. Similarly, you can study the
reflection of ripples in a water basin. You will be surprised to know that the same physical laws
govern the reflection (refraction) of all waves, including light. We will now consider reflection
and refraction of waves using Huygens' wave theory.
310
Can you locate the reflected wavefront? To discover this, we note that at the instant D is just
disturbed, the wavelet from A has grown for time DC / v and has travelled to E so that the
distance AE is equal to DC. We can draw a circle of radius AE (= DC) to represent this wavelet
with A as centre. Similarly, we can draw many circles from the intermediate points. The tangent
or the envelope to these circles from D defines the reflected wavefront.
From Fig, 7.2 it is clear that s ACD and DEA, are congruent. Hence
CAD ADE
or
i r (7.1)
That is, the angle of incidence is equal to the angle of reflection. Moreover, you will note that the
incident ray, the reflected ray and the normal at the point of incidence lie in the plane of the
paper.
In this connection it is important to mention here that the reflected wavefront undergoes a phase
change of . In fact, it is true for any wave travelling in a rarer medium (air) and undergoing
reflection at the interface with a denser medium (water). However, the reverse is not true.
311
1 2
v1 v2
Let v t and t respectively denote the speed and wavelength of a wave of frequency f in a rarer
medium. On being refracted at the interface of a denser medium, let its speed and wavelength be
v d and d , respectively. Mathematically, we can connect these quantities through the relation
vt vd
(7.2)
t d
This relation holds for waves in water, air and string alike.
Using Huygens’ principle, you can prove the laws of refraction as well (TQ1). But you will agree
that Huygens' method is essentially geometrical and can be used when the wave is either reflected
or refracted at the interface. You may now ask: Can we apply this method to study partial
reflection and refraction, as in the case of two strings having different mass per unit lengths? In
principle, we .can do so but it is more convenient to study partial reflection and refraction in
terms of impedance offered by a medium. To this end, we normally compute reflection and
transmission amplitude coefficients. You will now learn to compute these in the following
section.
312
7.3.1 Transverse Waves
Let us reconsider the strings AO and OB joined together at O and kept under the same tension T.
Let us assume that they offer characteristic impedances of Z 1 and Z 2 , respectively. A wave
travelling in the positive x -direction (Fig. 7.4) gets partly reflected and partly transmitted at O.
The particle displacements due to incident, reflected and transmitted waves can be written as:
yi ( x, t ) ai sin( 0 t k1 x) (7.3)
y r ( x, t ) ar sin( 0 t k1 x) (7.4)
and
yt ( x, t ) at sin( 0 t k 2 x) (7.5)
Fig. 7.4 Transverse waves on strings having different mass per unit length
where the subscripts i , r and t on displacements and the amplitudes refer to the incident,
reflected and the transmitted waves, respectively. You will note that the angular frequency of
these waves remains the same. Moreover the propagation constant for the incident and the
reflected waves is the same but differs for the transmitted wave. Do you know why? This is
because the wave speed changes as the density of the medium changes. You would also note that
for the reflected wave we have used a positive sign before k1 x . This is because it is travelling in
the negative x -direction.
To give physical meaning to the reflection and transmission coefficients, we have to consider the
boundary conditions. The boundary conditions are the conditions which must be satisfied at the
interface where the two media meet. Here the total displacement and the total transverse
component of tension on one side of the boundary are the result of the combination of incident
and reflected waves. So the boundary conditions in this case are:
1. The particle displacements immediately to the left and the right of the boundary (i.e. at
y ( x, t )
x = 0 must also be the same. This implies that the particle velocities should also
t
be the same.
y ( x, t )
2. The transverse components of tension T must also be the same immediately
x
on two sides of the boundary.
313
These conditions require:
yi ( x, t ) x0 y r ( x, t ) x0 yt ( x, t ) x0 (7.6)
and
yi y r yt
T T T (7.7)
x x 0 x x 0 x x 0
We know that
2 2f
k1T T T 2 fm1v1 = 2 f Z1
1 v1
In arriving at this result, we have used Eqs. (6.22) and (6.36b). Similarly, you can write
k 2T 2 f Z 2
2 f Z1 (ai ar ) 2 f Z 2 at
or
Z 1 ( ai a r ) Z 2 at (7.10)
Eqs. (7.8) and (7.10) enable us to calculate the ratios ar / ai and at / ai . These ratios give us the
fractions of the incident amplitude reflected and transmitted at the boundary. These ratios are
usually called the reflection and transmission amplitude coefficients. We will denote these by the
symbols R12 and T12 :
a r Z1 Z 2
R12 (7.11)
ai Z 1 Z 2
314
and
at 2Z1
T12 (7.12)
ai Z 1 Z 2
We note that the reflection and transmission amplitude coefficients depend only on the
impedances of the two media.
Let us now consider the implications of results arrived at in Eqs. (7.11) and (7.12):
(i) Assume that the string is rigidly fixed to a wall. This means that the second medium is
extremely heavy, meaning thereby that Z 2 . In such a case, R12 1 and T12 = 0.
This result implies that a r ai and a t = 0. That is, the amplitude of reflected wave is
equal to the amplitude of incident wave with just a reversal of sign and there is no
transmitted wave. This means that the incident wave suffers a change of phase of on
reflection from a denser medium.
(ii) When Z 2 Z1 , i.e. the second string (medium) is denser, R12 is still negative, implying
a phase change of on reflection. In this case, however, the incident wave is partly
reflected and partly transmitted.
(iv) When Z1 Z 2 , R12 = 0 showing no reflected wave. In this case T12 = 1, which gives
at ai . This means that the amplitude of a transmitted wave is equal to the amplitude of
the incident wave.
The points, (i), (ii) and (iii) above clearly show that if a wave travelling in a medium of lower
impedance meets the boundary of a medium of higher impedance (air to water), the reflected
wave undergoes a phase change of . If, however, a wave travelling through a medium of higher
impedance meets the boundary of a medium of lower impedance (water to air), no change of
phase takes place for the reflected wave. You may also note that T12 is always positive,
indicating that there is no change of phase for the transmitted wave in any case. These results are
depicted in Fig. 7.5.
From Eq. (6.36a, b), you will recall that for a given tension, the wave velocity will be lower in a
medium of higher impedance. Using this observation, can you now connect the above discussion
with the one given in Sec. 7.2.1? Is there not a one to one correspondence between the two cases?
This explains why we expected all waves, whether sound waves, water waves, waves on string or
light waves to follow the same laws.
Coming to the point (iv) above, we note that when Z1 Z 2 , the two strings are made up of the
same material and there effectively exists no boundary. That is why there is no reflection at all.
315
SAQ 1
Two strings of linear densities m1 and m 2 (= 4 m1 ) are joined together and stretched with the
same tension T. For transverse wave, calculate the reflection and transmission amplitude
coefficients.
Fig. 7.5 Reflected and transmitted waves when the incident wave (a) travels front a medium of
lower impedance to a medium of higher Impedance, and (b) when reverse b the case
(i) The particle displacement ( x, t ) is continuous at the boundary. That is, it has the same
value immediately to the left and right of the boundary at x 0 .
316
(ii) The excess pressure is also the same immediately on the two sides of the boundary.
Using the boundary conditions stated above you can show that the reflected and the transmitted
longitudinal waves obey the same characteristics as transverse waves (TQ5).
You have seen in Unit 6 that when a string of mass per unit length m vibrates with amplitude a
and angular frequency 0 , the total energy is given by
1
E ma 2 02 (7.13)
2
Let us assume that the wave is travelling with a speed v . Then the rate at which the energy is
carried along the string is obtained by multiplying the expression for energy with the speed v of
1
the wave and is equal to ma 2 02 v .
2
Now refer to the case of the transverse waves discussed in Section 7.3.1. The rate at which the
energy reaches the boundary alongwith the incident wave is given by
1 1
Pi m1 ai2 02 v Z1 02 ai2 (7.14)
2 2
Similarly, the rates at which the energy leaves the boundary alongwith the reflected and the
transmitted waves are
1
Pr Z1 02 a r2 (7.15)
2
and
1
Pt Z1 02 at2 (7.16)
2
Using Eqs. (7.8) and (7.10), we can write a r and a t in terms of a i . Substituting the resulting
expression in Eqs. (7.15) and (7.16) we find that
2
1 Z Z2 2
Pr Z1 02 1 ai (7.17)
2 Z1 Z 2
and
317
2
1 2Z1 2
Pt Z1 02 ai (7.18)
2 Z1 Z 2
These results can be used to obtain the reflection and transmission energy coefficients R E and
RE :
2
rate at which energy is reflected at the interface Pr Z1 Z 2
RE (7.19)
rate at which energy is incident at the interface Pi Z1 Z 2
rate at which energy is transmitted at the interface Pt 4Z1 Z 2
TE (7.20)
rate at which energy is incident at the interface Pi ( Z1 Z 2 ) 2
We note from Eq. (7.19) that if Z1 Z 2 : (which is also possible if we have T1m1 T2 m2 ,
RE 0 . That is, no energy is reflected back when impedances match. Such an impedance
matching plays a very important role in energy transmission. Long distance cables carrying
energy need to be matched accurately at all joints; otherwise a lot of energy will be wasted due to
reflection. We need impedance matching when we wish to transfer sound energy from air in a
loudspeaker to the air of the room. Similarly, when light waves travel from air into glass lens or a
slab, we wish not to have reflections (as it will reduce intensity).
SAQ 2
Show that the energy is conserved when a transverse wave meets the boundary between two
media of characteristic impedances Z 1 and Z 2 .
For longitudinal waves, it is customary to calculate energy transfer in terms of their intensity.
From Unit 6 we recall that the intensity of sound waves in a gas is given by
1
I a 2 02 v
2
= 2 2 f 2 a 2 Z (7.21)
where Z is the impedance offered by the medium to wave motion. Hence the incident, reflected
and transmitted wave intensities can be written as
I i 2 2 f 2 ai2 Z1 (7.22)
I r 2 f a Z1
2 2 2
r (7.23)
and
I t 2 2 f 2 at2 Z 2 (7.24)
Using these equations, you can easily show that the reflection and transmission energy
coefficients are given by
2
I Z Z2
RE r 1 (7.25)
I i Z1 Z 2
318
It 4Z1 Z 2
TE (7.26)
I i ( Z1 Z 2 ) 2
You will observe that these relations are the same as for transverse waves. This means that the
same conclusions hold even for longitudinal waves.
SAQ 3
Sound waves are incident on a water-steel interface. Show that 86% of the energy is reflected
back. Impedances of water and steel are respectively 1.43 10 6 Nm 3 s and 3.9 10 7 Nm 3 s .
This dialogue suggests that frequency can change with the speed of the observer or source. You
all must have heard the whistle of a moving train. What do you feel when the train approaches
you? The pitch of the whistle seems to rise. But when the engine passes by, the pitch appears to
decrease. The apparent change of frequency due to the relative motion between the source and the
observer (or the listener) is known as the Doppler Effect.
In general, when the source approaches the listener or the listener approaches the source, or both
approach each other, the apparent frequency is higher than the actual frequency of the sound
produced by the source. Similarly when the source moves away from the listener, or when the
listener moves away from the source, or when both move away from each other, the apparent
frequency is lower than the actual frequency of the sound produced by the source.
Do you know that Doppler shift in ultrasound reflected from moving body tissues allows
measurement of blood flow? It is commonly used by obstetricians to detect foetal heartbeat. Do
you know how it arises? As the heart muscle pulsates, the reflected ultrasound waves are Doppler
shifted from the incident waves. Similarly, a sonar makes use of the Doppler effect in
determining the velocity of a submarine relative to a ship.
The electromagnetic waves, including light, are also subject to the Doppler effect. In air
navigation, radar works by measuring the Doppler shift of high frequency radio waves reflected
from moving aeroplanes. The Doppler shift of starlight allows us to study stellar motion. When
we examine light from stars in a spectrograph, we observe several spectral lines. These lines are
slightly shifted as compared to the corresponding lines from the same elements on the earth. This
shift is generally towards the red-end and is attributed to stellar motion. This is illustrated in Fig.
7.6 for hydrogen atoms in a double star system. (The Doppler shift of light from distant galaxies
is an evidence that our universe is expanding.)
319
Fig. 7.6 The wavelength of light emitted by hydrogen atoms in a binary star reveals the stellar
motion
To study the Doppler effect for sound waves, we have to consider the following situations:
(i) Whether the source is in motion, or the observer is in motion, or both are in motion.
(ii) Whether the motion is along the line joining the source and the observer, or inclined (at
an angle) to it.
(iii) Whether the direction of motion of the medium is along or opposite to the direction of
propagation of sound.
(iv) Whether the speed of the source is greater or smaller than the speed of sound produced by
it.
320
Fig. 7.7 Successive wavefronts emitted by a moving source
Representing the same situation in terms of waves, as shown in Fig. 7.8a, we find that if v the
speed of sound produced, f waves occupy a length v in one second, if the source is stationary.
After one second, when the source has moved a distance u , towards the listener, the same
number of waves get crowded a length (v u s ) as shown in Fig. 7.8 (b).
v us
'
f
v v
f ' f (7.27)
' v us
321
If, however, the source moves away from the observer (in a direction opposite to sound), u s , is
negative and Eq. (7.27) becomes
v
f ' f (7.28)
v us
To fix up the ideas discussed above, you may now like to solve a SAQ.
SAQ4
A person is standing near a railway track. A train approaches him/her with a speed of 72 km h 1 .
The apparent frequency of the whistle heard by the person is 700 Hz. What is the actual
frequency of the whistle? Use the speed of sound in air as 350 ms 1 .
If the listener moves towards the source, u s is negative, and the apparent frequency is given by
v u0
f ' f (7.30)
v
322
length of block of waves received
f '
reduced wavelength
v u0
= (7.31)
v u0
You may now ask: Is there any difference in the apparent frequency when the source approaches
the listener or the listener approaches the source with the same velocity? Eq. (7.31) tells us that
the apparent frequency will be different in these cases.
For electromagnetic waves, Eq. (7.31) has to be modified. For sound, u 0 and u s are measured
relative to the medium. This is because the medium determines the wave speed. However, e.m.
waves do not require a medium for propagation so that their speed relative to source or the
observer is always the same. For these waves we have to consider only the relative motion of the
source and the observer. If u 0 is the speed of source relative to observer, and u s v , we can
rewrite Eq. (7.31) as
u
f ' f
v us
1
u
= f 1 s (7.32)
v
Using binomial expansion and retaining only first order terms in (u s / u v ) , we get
u
f ' f 1 s
v
In air navigation, we take u s to be twice the approach velocity of the aeroplane. This is because
the radar detects e.m. waves sent by it and reflected back by the aeroplane.
SAQ 5
A stationary observer notes that the spectral line of wavelength 4000 A emitted by a star is shifted
towards the red from its normal position by 100 A. Calculate the speed of the star in the line of
sight? Speed of light = 3 108 ms 1 .
Now you may ask: What happens when the speed of source exceeds the speed of sound waves as
for supersonic planes? To discover the answer to this question, let us see if we can draw wave
patterns similar to those shown in Fig. 7.7.
323
us v
Fig, 7.10 Schematic representation of piling of waves when the source moves at the wave speed
Let us suppose that the source is at point A at t = 0. After time t , the waves emitted at A are on a
sphere of radius vt . Since u s v , the distance travelled by the source AS = u s t is more than the
distance travelled by sound waves. The waves emitted at successive points, B,C,D, E.... are on the
line A'S, where the circles are most crowded. We thus see that sound waves pile up on a cone
whose half angle is given by sin v / u s as shown in Fig. 7.11a. No sound waves are present
outside this cone. The velocity of sound waves is normal to the surface of the cone. When this
cone hits an observer, he detects the sudden arrival of a large amplitude wave, known as a shock
wave. A supersonic aircraft generates shock waves, also called sonic booms, due to the formation
of two principal shock fronts; one at its nose and the other at its tail (Fig. 7.11b). A strong boom
can break window glasses or cause other damage to buildings.
324
Fig. 7.11 (a) Shock waves created by a sound moving faster than the speed of sound,
(b) Sonic booms produced by a supersonic aircraft
Shock waves are also generated in a ripple tank by a moving source for Mach numbers greater
than one. You can also observe that shock waves are formed by a boat moving faster than the
speed of water waves.
7.7 SUMMARY
The locus of points in a given phase is called a wavefront. The shape of a wavefront
depends on the nature of the source.
According to Huygens, each point on a wavefront becomes a fresh source of secondary
wavelets, which move out in all directions with the speed of the wave in that medium.
When waves travelling through one medium meet the boundary of another medium with
a different impedance, they are partly reflected and partly transmitted. The reflection and
transmission amplitude coefficients are respectively given by
Z1 Z 2
R12
Z1 Z 2
and
2Z1
T12
Z1 Z 2
When a wave travelling from a medium of lower impedance is reflected from a medium
of higher impedance, a phase change of takes place.
Due to the relative motion between the source of sound waves and the observer (listener),
the apparent frequency of the sound is different from the actual one. This is known as the
Doppler effect. The Doppler shifted frequency (when both approach each other) is given
by
v u0
f ' f
v vs
325
of the amplitude of sound wave that enters the second medium to the amplitude of the
incident wave. Use = 1.29 kgm3 . Speeds of sound in air and water are 350 ms 1 and
1500 ms 1 respectively.
3. A rope is made up of a number of identical strands twisted together. At one point, the
rope becomes frayed so that only a single strand continues (Fig. 7.12). The rope is held
under tension and a wave of amplitude 1.0 cm is sent from the single strand. The wave
reflected back along the single strand has an amplitude of 0.45 cm. How many strands are
in the rope?
7.9 SOLUTIONS
SAQs
1. We know that impedance is related to tension and mass per unit length through the
relation
Z mT
For the given strings
Z1 m1T and Z 2 m2T
Z1 m1 1
Z 2 m2 2
From Eqs. (7.11) and (7.12), we note that the reflection and transmission amplitude
coefficients are
a r Z1 Z 2 Z1 / Z 2 1 1 / 2 1 1
R12
ai Z 1 Z 2 Z 1 / Z 2 1 1 / 2 1 3
and
at 2Z1 2Z1 / Z 2 2
T12
ai Z 1 Z 2 Z 1 / Z 2 1 3
The negative sign in R12 implies a phase change of at the interface.
2. From Eq. (7.14) we know that the rate at which energy reaches the boundary is given by
1
P1 Z1 02 ai2
2
Similarly, the rate at which energy leaves the boundary with reflected and transmitted
326
waves is given by
1 1
P2 Z1 02 a r2 Z1 02 at2
2 2
On substituting for a r and a t we get
2 2
1 Z Z2 2 1 2Z1 2
P2 Z1 02 1 ai Z 2 02 ai
2 Z1 Z 2 2 Z1 Z 2
1
2 Z1 Z 2
2
4Z1 Z 2 2
Z1 0 ai
2 Z1 Z 2 ( Z1 Z 2 ) 2
1
= Z1 02 ai2
2
Since the rate at which energy arrives al the interface is equal to the rate at which energy
leaves the interface (with reflected and transmitted waves), we can say that energy is
conserved in this process.
37.57
= = 0.86
40.43
This means that when sound waves are incident on water-steel interface, only 86% of the
energy is reflected back.
5. Since the wavelength increases, we can say that the star is moving away along the line of
sight. This means that frequency decreases. Using Eq. (7.28) for the case of light you can
write
c
f ' f
c us
1
us
= f 1
c
327
us
= f 1 for u s c
c
Since v c / , you can write
1 1 us
1
' c
or
c
us ( ' )
'
Here ' = 4100 Å and c 3 108 ms 1
Hence,
3 10 8 ms 1
us 100
4000
= 7.3 10 6 ms 1
= 7.3 10 3 kms1
TQs
1. Refer to Fig 7.13. AB represents a part of a wavefront moving towards the interface, S1S2
which separates the two media say air and water. Let us assume that wave speeds in
medium 1 and medium 2 be v1 and v 2 respectively.
The wavefront will first strike at C and then at successive points towards D. The point B
on the wavefront reaches D at a time t B' D / v1 later than the point A reaches C. From
each point on S1S2, a secondary wavelet starts growing into the second medium at speed
v 2 . At the instant when D is just disturbed, the wavelet from C has grown for time
t ( B' D / v1 ) and acquired the radius
B' D
CC ' v2 (i)
v1
328
You can represent this wavelet by drawing an arc of radius CC ' with C as centre. Draw
a tangent DC" from D to this arc If you repeat this process for other intermediate points
between C and D you will observe that DC" is a common tangent to all of them. Thus,
DC" represents the refracted wavefront.
That is, the sine of the angle of incidence to the sine of angle of refraction of the
wavefront is equal to the ratio of the wave speeds and is a constant This constant is
known as the refractive index of medium 2 with respect to medium I. We denote it by the
symbol 12 .
2. 1 1.29 kgm3
2 1000 kgm3
v1 350 ms 1
v2 1500 ms 1
329
( Z1 / Z 2 ) 1
=
( Z1 / Z 2 ) 1
4. In this case, the velocity of the car is not directed towards the sound source (Fig. 7.14a),
and we have to find the component of the velocity vector directed towards the source.
Referring to Fig. 7.14b it is given by
Fig. 7.14 (a) Observer moving along a line not intersecting the line of motion of source
330
(b) The component of velocity of observer towards the source is responsible for
Doppler shift
x x
u 0 cos u 0 = 20 m s 1
x 2 20 2 x 2 20 2
Then the space dependence Doppler-shifted frequency is given by
v u 0 cos
v( x) f 0
v
x
= 500 Hz1 0.06
x 2 20 2
You can plot this as a function of x for 100 m < x < 100 m. At x = 0, the car is moving
perpendicular to the wave and at the instant when the car passes this point, the driver
hears the true frequency, 500 Hz.
5. The particle displacement for the incident, reflected and transmitted waves are
i ( x, t ) ai sin( 0 t k1 x) (i)
r ( x, t ) ar sin( 0 t k1 x) (ii)
and
t ( x, t ) at sin( 0 t k 2 x) (iii)
1. The particle displacement ( x, t ) is continuous at the boundary. That is, it has the same
value immediately to the left and the right of the boundary at x = 0.
2. The excess pressure is the same on the two sides of the boundary.
331
For a longitudinal wave, p E , where E is the elasticity. Since E p0 ,
x
cp
where and p 0 is the equilibrium pressure, we find that p 0 cancels out on both
cv
sides and the second condition implies that
i r t
(v)
x x x
332
UNIT 8 SUPERPOSITION OF WAVES-1
Structure
8.1 Introduction
Objectives
8.2 Principle of Superposition of Waves
8.3 Stationary Waves
Velocity of a Particle and Strain at any Point in a Stationary Wave
Harmonics in Stationary Waves
Properties of Stationary Waves
8.4 Wave Groups and Group Velocity
8.5 Beats
8.6 Summary
8.7 Terminal Questions
8.8 Solutions
8.1 INTRODUCTION
You have studied in Unit 2, how a particle acted upon simultaneously by two simple harmonic
oscillations gives rise to the formation of Lissajous figures.
You have also read about the general characteristics of waves in Unit 6; and of their behaviour at
the interface of two media in Unit 7. In this unit you will study the principle of superposition of
waves. Under certain conditions, the superposition of waves leads to some interesting phenomena
like the formation of stationary waves, beats, wave groups, interference, diffraction etc. In the
present unit you will study the phenomena of stationary waves, wave group and beats. The other
two topics, viz. Interference and Diffraction, will be discussed in Unit 9.
In the present unit you will study the basic features, especially the sound producing part of the
woodwind instruments. There are two basic types of pipes, viz. flute pipe and reed pipe, which
you will study in this unit. Stationary waves are formed when two waves of the same angular
frequency (i.e., same ( ), same wavelength (i.e., of the same wave vector or propagation
constant k) and of same amplitude, travelling in opposite directions superpose on each other. On
the other hand, if two sound waves of slightly different frequencies are superposed, they produce
beats.
Wave groups, sometimes also called the wave packets, are the result of superposition of waves of
slightly different frequencies. The concept of wave packet is of great importance in the study of
quantum mechanics, which we consider later.
In the next Unit you will study the superposition of two waves, which leads to the phenomena of
interference. There you will also study the necessary conditions for the interference of two waves.
Towards the end, you will learn about diffraction of waves and some typical cases of diffraction
phenomena.
Objectives
After going through this unit, you will be able to
Describe the principle of superposition of waves
Explain the ideas underlying the formation of stationary waves
Identify the positions of nodes and antinodes on a stationary wave
List the characteristics of stationary waves
Describe the formation of wave groups
333
Compute the value of group velocity knowing the dependence of wave velocity on
wavelength
Calculate the number of beats produced if the frequencies of two superposing notes are
known.
Two or more waves can traverse the same path in a given space, independent of one another. This
means that the resultant displacement of a particle at a given time is simply the algebraic sum of
the displacements that are given to the particle by the individual waves. In other words, we can
say that the resultant displacement of the particle is found simply by adding algebraically the
displacements due to the individual waves. This is known as superposition of waves.
An interesting case of superposition of waves is that of radio waves. You know that radio waves
of different frequencies are transmitted hy different radio stations to broadcast their programmes.
When they fall on the receiving antenna, the resultant electric current set up in the antenna is
quite complex because of the superposition of different waves. Nevertheless, we find that we can
still tune to a particular station. That is, out of the many, we can still choose and pick up the
particular wave we want. In other words, if we have a wave group obtained by the superposition
of a large number of individual waves, we can still separate the different waves that were
superposed. This is indicative of the individual behaviour of waves, which is the basis of the
superposition principle in waves.
Now you can demonstrate the principle of superposition by considering two pulses travelling on a
rope in opposite directions as shown in Fig. 8.1. Before and after crossing each other, they act
completely independently. At the time of crossing, the resultant displacement is the algebraic sum
of the individual displacements.
You have also studied in Unit 2, the mathematical basis for the superposition of oscillations. It
lies in the linearity of the equation. Consider two waves acting independently on a particle at any
position x. Let y1 ( x, t ) and y 2 ( x, t ) be the displacements of the particle at the instant of time t
334
due to the two waves. Then the resultant displacement Y ( x, t ) of the particle is mathematically
written as:
Y ( x, t ) y1 ( x, t ) y 2 ( x, t ) (8.1)
You have studied in Unit 6, that a wave is essentially characterised by its amplitude, angular
frequency, wave vector and phase. Depending on which of these components are the same or
different, you will study the various phenomena in Physics due to the superposition of waves. Let
us consider some of these phenomena. For this, you consider the superposition of the following
pair of waves.
You will draw the following conclusion from above combination of waves.
Eq. (8.3) implies that the resultant wave has same the frequency and phase and the resultant
amplitude is (a1 a2 ) . It is shown in Fig. 8.2.
Fig. 8.2 Superposition of two waves of same frequency, wave vector and phase, but differing
amplitudes a1 and a 2
335
(c) In case (iii), the frequency and wave vector k of the two waves differ.
Now let us consider the case when the frequencies and wave vectors of two waves differ slightly.
In such a case, irrespective of phase difference the superposition results in an interesting
phenomenon of 'Beats'. If however many waves of slightly different frequencies superpose, then
they form a waves group (or a wave packet). This gives rise to group velocity, quite distinct from
the wave velocity. You will study group velocity in detail in Section 8.4.
(d) In case (iv), the waves equations have different signs before the wave vector (k). In this case,
the first wave y1 ( x, t ) is propagating along the positive direction of x -axis, while the other
wave, y 2 ( x, t ) , is propagating in negative direction along x -axis. This implies that they are
propagating in opposite directions. When such kind of waves superpose then stationary or
standing waves are produced. You will study stationary waves in Section 8.3.
The reflection of the incident wave can take place at a fixed boundary (like that of a string fixed
to a wall, or the closed end of an organ pipe) or at a free boundary (like the free end of a string, or
the open end of an organ pipe). We have learnt in the last Unit that at a fixed boundary, the
displacement y( x, t ) stays zero, and the reflected wave changes its sign. At a free boundary,
however, the reflected wave has the same sign as the incident wave. In other words, at a fixed
boundary, a phase change of takes place, while at a free boundary, no such change of phase
takes place.
Let us consider the case where the reflection is taking place at a free boundary. In this case, the
resultant displacement is given by:
From Eq. (8.5) you see that the amplitude is given by (2a cos kx) which is not fixed. It is
dependent (or varies harmonically) on the position x of the particle. Further, the resultant motion
has the same frequency and the wavelength as the individual waves.
Looking at equation (8.4) we note that the particles distributed along the x -axis execute
vibrations perpendicular to the x -axis. The amplitudes with which they execute these vibrations
are different at different positions along the x -axis. However, the time period of vibrations of all
the particles is same.
336
We note that Eq. (8.5) does not represent a travelling wave since the argument of the sine
function is independent of the space variable x . We thus see that although we started with two
waves propagating in opposite directions, we have ended up with something that does not
propagate in space. The wave that does not travel (or propagate) is called a stationary (or a
standing) wave. Since it does not propagate, it transports no energy along with it.
From equation (8.5) it is clear that the displacement Y ( x, t ) is maximum when
2
cos kx cos x 1 (8.6)
and minimum when
2
cos kx cos x0 (8.7)
2 2
To satisfy Eq. (8.6) we require, x m . Similarly Eq, (8.7) requires x (2m 1)
2
with m = 0, 1, 2, ... These give the points of maximum displacement at x = 0, / 2 , , ...,
m / 2 ; and minimum displacement at x = / 4 , 3 / 4 , …, (2m 1) / 4 .
The points of maximum displacement are called 'Antinodes', while those of minimum
displacement are called 'Nodes'. The distance between any two consecutive nodes or antinodes is
/ 2 , while that between a node and an antinode is / 4 (Fig. 8.3).
From the above discussion you have learnt that a stationary wave results due to the superposition
of two identical progressive waves travelling in opposite directions. The result is a non-
progressive wave in which the disturbance is not handed over from one particle to the next. The
space (or the region) where the two waves superpose gets divided into compartments or segments
(Fig. 8.3). Each segment ends with points called the nodes where the displacement of the particles
is always zero.
Fig 8.3 The envelope of a standing wave showing the pattern of nodes and antinodes
337
The particles at the central points of these segments (called the antinodes) execute vibrations with
maximum amplitude. The particles lying in-between the nodes and the antinodes execute
vibrations with amplitudes lying in between zero and the maximum amplitude. This is shown in
Fig. 8.4.
The particle a, for example, is always at rest. The particle b always executes vibration with
maximum amplitude, and the particle c always with intermediate amplitude as shown in Fig. 8.4.
Fig 8.4 Stationary wave, with arrowheads indicating the amplitudes with which various particles
vibrate
SAQ 1
Derive equation for the displacement of a particle lying on a standing wave on a string fixed at
both ends. Will the fixed end of a string be a node or an antinode? If the standing wave is in an
open-ended air pipe, will there be a node or an antinode at the end? How do you explain the
absence of energy flow in a standing wave?
Y
velocity = 2a cos kx cos t (8.8)
t
The velocity is maximum when cos kx = ± 1, i.e. at points where x = 0, / 2 , , ..., m / 2 .
(see Eq. (8.6) and the discussion that follows). The velocity is minimum (zero) when cos kx = 0,
i.e. at points where x = / 4 , 3 / 4 , …, (2m 1) / 4 .It means that the velocity is maximum at
the antinodes where the displacement is also maximum. The velocity is zero at the nodes where
the displacement is zero. At points in between the antinodes and nodes, the velocity gradually
338
decreases from maximum at the antinodes to zero at the nodes. The lengths of the arrowheads in
Fig. 8.4 may also be taken to represent the velocities of the particles in a stationary wave.
The strain on a particle in a stationary wave can be calculated by differentiating the resultant
amplitude i.e. Y ( x, t ) w.r.t. x keeping t constant. If we differentiate Eq. (8.5) w.r.t. x we get
strain
y
2ak sin kx sin t (8.9)
x
You can show that the strain is maximum for the particles at the nodes where the displacement
and the velocity are zero. This can also be visualised from Fig. 8.4. The particles at the nodes are
stretched by particles moving in opposite directions. The strain is minimum at the antinodes
where the displacement and velocity are maximum. Again referring to Fig. 8.4, we can see that
the particles at the antinodes always move along with the particles at their sides, not causing
much strain on particles at the antinodes.
In the case of stationary waves, the particles get divided into segments like the P, Q and R in as
shown in Fig. 8.4. Particles in one segment always move along in the same direction. When
particles in segment P move up, those in Q move down. When those in Q move up, the ones in P
move down. That is, in any two adjacent segments, particles move in opposite directions.
All particles in a particular segment reach the extreme positions at the same time, and also pass
through the mean positions at the same time. This is shown in Fig. 8.5. All this is possible since
all the particles have the same time period T but have different velocities. The particles which
have to cover larger distances have greater velocities. Those which have to cover smaller
distances, have smaller velocities.
Fig. 8.5 Stationary waves on a string fixed at both ends. Shape of the string at different times
during a time period is shown
339
Now coming to the individual particle, we can see when its velocity is maximum, and when it is
zero. Writing Eq. (8.8) as:
Y 2
4 af cos x cos 2 f t
t
4 a 2 2
= cos x cos t
T T
you can see that the particle velocity is zero for t T / 4 and 3T / 4 , and is maximum for t 0 ,
T / 2 and T . Thus, during each time period, the particles of the medium have their maximum
velocity when they pass through the mean position, and have zero velocity when they are at the
extreme positions. Now in the next section you will study the conditions for producing different
harmonics in stationary waves.
If the length of the string is L , the wavelength of the possible stationary waves on this string,
2 L
starting from the longest wavelength are: 2 L , L , ,..., etc. (See Fig. 8.6).
3 2
These wavelengths determine the frequencies of oscillation of the string through the relation
f v . Here v the velocity of the transverse wave on the string. It is given by the relation
T
v
where T is the tension in the string, and is the linear mass density (mass per unit length) of
the string.
The lowest frequency f 0 of vibration is called the fundamental frequency. It is given by:
v 1 T
f0 (8.10)
2L
The other frequencies are called the overtones, and are integral multiples of the fundamental
frequency f 0 (See Fig. 8.6).
The fundamental frequency is also called the first harmonic. The first overtone, with frequency
f 2 f 0 , is called the second harmonic. The second overtone, with frequency f 3 f 0 , is called
the third harmonic, and so on.
340
Fig. 8.6 Allowed stationary waves on a string or length L clamped at its ends
The musical instruments based on the principle of standing wave are flute, reed, etc. The primary
elements, which determine the tone, quality and overall sound are (1) the source of noise or
vibration (2) the size and shape of the bore, and (3) the size and positions of the finger holes. The
quality of woodwind tones depends on the combination of physical and musical experience. From
the physics point of view, the air is stored under pressure in the wind chest. A large reservoir is
required to keep the pressure steady, while the various combinations of notes are played with
fingers. In the above instruments one end is open, making them open ended organ pipes. The
closed end of an organ pipe acts as a fixed boundary, while the open end as a free boundary. At
the closed end there is always a node, and at the open end there is always an antinode.
For a pipe having one end closed, the fundamental wavelength is 4L . This gives the
fundamental frequency f 0 v / 4 L . In such a pipe, the even-numbered harmonics are absent
(See Fig. 8.7). For a both ended open pipe, the fundamental wavelength is 2L , giving
fundamental frequency f 0 v / 2 L .
Fig. 8.7 Modes of vibration of longitudinal stationary waves in organ pipes with (a) one end
closed, and (b) both ends open
341
SAQ 2
(a) A piano string of length 1 m (mass per unit length 0.015 kgm1 ) fixed at both ends is to
be used to strike a fundamental note of frequency f = 220 Hz. Find the tension to be
applied to the string.
(b) Estimate the frequency of the fundamental mode in a one end closed organ pipe of length
1.0 m. Use velocity of sound v =350 m/s. What happens to frequency if the pipe is
overblown?
(i) Stationary waves are not progressive. In these the disturbance is not handed over from
one particle to the next.
(ii) The amplitude of each particle is not the same. It is maximum at the antinodes and zero at
the nodes. In between, it gradually decreases from that at the antinode to the one at the
node, i.e., zero.
(iii) The distance between two consecutive nodes or two consecutive antinodes is half the
wavelength of the stationary wave. The medium splits into segments, with the length of
each segment equal to half the wavelength.
(iv) All the particles between two consecutive nodes are in phase, i.e. they reach their
maximum and minimum displacement positions (mean positions) at the same time. The
phase of particles, in one segment is opposite to that of particles in the adjoining segment.
(v) The velocity of particles at the nodes is zero. The velocity of particles at the antinodes is
maximum. For particles in between, the velocity gradually decreases from that at the
antinodes to the one at the nodes (i.e. zero).
Further, writing
342
1 2 k1 k 2
av and k av ,
2 2
k
Y ( x, t ) 2a sin( avt k av x) cos t x (8.12)
2 2
Now let us see what the new wave form represented by Eq. (8.12) looks like. Firstly, its
amplitude is twice that of the amplitude of either wave. Secondly, it is made up of two parts. The
faster varying part (i.e. sine part) has a frequency which is the mean of the frequencies of the two
component waves. The slowly varying part (i.e. cosine part) has a frequency which is half of the
difference of the two frequencies. The propagation vector of the slowly varying part of the
superposed wave is k / 2 . It acts as an envelope over the faster varying part as shown in
Fig.8.8.
The superposition, as you can see in Fig. 8.8, results in the formation of groups (or segments)
called wave groups (or wave packets). A wave group can travel with a velocity which may be
different from that of the individual waves, or of the resultant wave. The velocity of the wave
group is called the group velocity. The ratio of angular frequency and wave vector of the slowly
moving part of the superposed wave is called group velocity. It is given by the following relation
/ 2
vg (8.13)
k / 2 k
If a group consists of a number of component waves with angular frequencies lying between 1
and 2 (with 1 2 ), and similarly in wave vector k1 and k 2 (with k1 k 2 ), the group
velocity v g is then written as :
d
vg (8.14)
k dk
343
Here d and dk represent the spreads (gaps between the maximum and the minimum) in
angular frequencies and propagation constants of the component waves that go to make a wave
group.
The velocity of the resultant superposed wave is called the phase velocity. You can obtain this
using Eq. (8. 12), i.e.
av
vp
k av
1 2
= v (say)
k1 k2
1 2 k1v k 2 v
then, vp v
k1 k 2 k1 k 2
( 2 ) / 2 vk1 vk2
and vg 1 v
(k1 k 2 ) / 2 k1 k 2
The group velocity is a more fundamental quantity in physics as the energy is transferred by the
wave with the group velocity. The relation between the phase and the group velocities is given
by:
d d
vg (vk)
dk dk
dv
vg v k
dk
If we write
k 2 / (8.15)
2
then, dk 2 d
Inserting this in Eq. (8.15), we get
2
dv
vg v
(2 / 2 )(d )
dv
= v (8.16)
d
344
This gives another relation connecting the phase and the group velocities. The wavelength of the
resultant wave is given by
2
k
2 4
e
k / 2 k
e 1 2
(8.17)
2 2 1
This gives the length (or the extent) of the wave group. We can see from Fig. 8.9 that the length
of the wave group is half of the wavelength of the enveloping wave, i.e., it is equal to e / 2 .
To illustrate the difference between phase and group velocities, we consider the striking example
of waves in deep water – called ―gravity waves." These waves are strongly dispersive. For them,
the phase velocity is found to be proportional to the square root of the wavelength, i.e.,
v p C1 / 2
or
v p C1k 1 / 2 (since k 2 / )
Here, the new constant, C1 C 2 . v p C1 k 1 / 2 , therefore, C1k 1 / 2 . Differentiating
with respect to k , we get,
d 1 1
vg C1 k 1 / 2 v p
dk 2 2
That is, the group velocity for gravity waves is just half of the phase velocity. In other words, for
these waves, the component wave crests move faster through the group as a whole.
345
SAQ 3
The phase velocity of a wave in a certain medium is represented by:
v C1 C2
8.5 BEATS
We have seen above that the superposition of two waves of slightly different angular frequencies
1 and 2 leads to the formation of wave groups. You may have noticed from Figs. 8.8 and Fig.
8.9 that we have plotted the resultant displacement Y ( x, t ) against distance x . This may be
called the superposition in space. For this we kept the time t as constant. We may now consider
another type of superposition, where we may plot Y ( x, t ) against t and call it superposition in
time. For this we may keep x as constant.
The superposition in time for sound waves leads to the interesting phenomenon of beats. The
beats are loud sounds which we hear at regular intervals of time depending on the difference in
frequencies of the two superposing waves. The beats are often used by musicians for tuning their
instruments.
Let us consider two waves of slightly different angular frequencies 1 and 2 , and of the same
amplitude a, proceeding in the same direction, as we have done in the last section. Let us fix the
spatial coordinate x in Eq. (8.10), say, at x = 0. This corresponds to an observer standing at x = 0,
and watching the waves passing by. He will observe a resultant waveform given by:
Like the earlier case discussed in Section 8.4, Eq. (8.18) indicates that the amplitude of the
resultant wave at a given point is not constant, but varies in time. This has an angular frequency
1 2
av . Its amplitude varies between 2a and zero, because of the presence of the
2
cos t term. This term acts as an envelope on the sine term.
2
If 1 and 2 are nearly equal, is small. In that case the amplitude of the resultant wave
varies slowly. The periodic rise and fall of this wave leads to the appearance of beats; or to the
hearing of loud sounds at regular intervals of time.
Beat are heard at the maxima of amplitude (See Fig. 8.10). They occur whenever
cos t 1. This is because the intensity of sound is directly proportional to the square of
2
the amplitude. The maximum amplitude occurs twice in every time period associated with the
346
angular frequency . Thus the frequency of beau is simply the difference of the two
2
component frequencies i.e. (1 2 ) .
In terms of frequencies f 1 and f 2 , the beat frequency is f f1 f 2 . The time
2
1
elapsed between any two consecutive beats, called the beat period = (see Fig. 8.10).
f
Fig. 8.10 Formation of bead due to superposition of two waves of nearly the same frequency
SAQ 4
When a certain note of a piano is sounded with a tuning fork of frequency 560 Hz, 6 beats are
heard every second. Find the frequency of the note.
8.6 SUMMARY
1. When two waves travelling through the same space superpose on each other, the resultant
displacement at any point is given by the algebraic sum of the individual displacements.
2. Stationary waves result because of the superposition of two waves of same amplitude,
frequency and wavelength travelling in opposite directions and confined between two
points.
3. On a stationary wave, nodes and antinodes are points of zero and maximum
displacement, respectively. Distance between any two consecutive nodes or antinodes is
half the wavelength of the stationary wave.
4. Both transverse and longitudinal waves can have different modes of vibration.
6. The number of beats produced per second is equal to the difference in the frequencies of
the two waves.
347
7. The velocity with which a wave group travels is called the group velocity. It is equal to
the wave velocity if the two component waves have the same velocity; otherwise, it is
different from the wave velocity.
8. The smaller the difference between the wavelengths of the component waves, the greater
the length of the wave group.
2. Fifty tuning forks are arranged in order of increasing frequency and any two successive
forks gives 5 beats per second when sounded together. If the last fork gives the octave of
the first, calculate the frequency of the latter. (A note is octave of another nole if its
frequency is double that of the other.)
3. A closed pipe, 25cm long, resounds when full of oxygen, to a given tuning fork. Find the
length of a closed pipe, full of hydrogen which will resound to the same tuning fork.
(Velocity of sound in oxygen = 320 m/s and velocity of sound in hydrogen = 1280 m/s).
4. The phase velocity (v) of transverse wave in a crystal of atomic separation d is given by
sin(kd / 2)
vC
(kd / 2)
where C is a constant. Show that the group velocity is C cos(kd / 2)
8.8 SOLUTIONS
SAQ l
Since there is a phase change of on reflection at the fixed end, the reflected wave is given by:
y 2 a sin(t kx)
This leads to the resultant displacement
Y ( x, t ) a sin( t kx) a sin( t kx)
= 2a sin kx cos t
= A cos t
with A 2a sin kx
At the fixed end, there is always a node, as the displacement is zero. In open-ended pipes, as
shown, there is always an antinode at the end.
A standing wave is formed because of a positive x directed incident wave, and a negative x
directed reflected wave. Each carry the same amount of energy in opposite directions. The net
energy flow is thus always zero.
348
SAQ2
(a) Wavelength of fundamental mode
2L 2 1 2 m
velocity of wave = 220 z 2 m 440 m / s
T
From v , T v2
= (440 m / s) 2 0.015kg / m
= 2.9 10 3 N
SAQ3
We know that
dv
vg v ,
d
dv
For the wave in question, C2
d
Inserting in the above equation,
v g C1 C2 C2 C1
SAQ4
Let the frequency of the note be f . Then
6 = [560 f ]
f 554 Hz or 566 Hz
In this case, the frequency of the note cannot be found without ambiguity. However, it is
either of the above two.
TQs
1(a) f =1.5 Hz
16
(b) m , n 1, 2, 3, … for positive moving wave.
16n 1
16
m , n 1, 2, 3, … for positive moving wave.
16n 1
(c) u 8 / 5 m / s , etc.
v 24 m / s , etc.
349
Then the frequency of the Second fork = n + 5 = n + (2 – 1)5
Frequency of the Third fork = n + 5 + 5 = n + 10 = n + (3 – 1)5
Frequency of the Fourth fork = n + 5 + 5 + 5 = n + 15 – n + (4 – 1)5
Frequency of the Fifth fork = n + 20 = n + (5 – 1)5
Therefore the frequency of the 50th fork
n (50 1) 5 n 245
Since the frequency of 50th fork is 2n then
n + 245 = 2n
So n = 245 Hz.
4. Group velocity
d
vg
dk
and
kv , we have
(kd / 2)
= C sin
(kd / 2)
Now
sin(kd / 2)
kv kC
(kd / 2)
2C
= sin(kd / 2)
d
or
d 2C kd d
vg cos
dk d s 2
or
v g C cos(kd / 2)
350
UNIT 9 SUPERPOSITION OF WAVES- II
Structure
9.1 Introduction Objectives
9.2 Interference Coherent Sources
Interference between Waves from Two Slits
Intensity Distribution in Interference Pattern
Interference in Thin Films
9.3 Diffraction
Different Types of Diffraction: Fraunhofer and Fresnel
Fraunhofer Diffraction by a Single Slit
Diffraction at a Straight Edge
9.4 Summary
9.5 Terminal Questions
9.6 Solutions
9.1 INTRODUCTION
In the last unit you have studied the principle of superposition of waves and employed it to study
the phenomena of formation of stationary waves, wave groups and beats.
You have also learnt about the superposition of two waves which have the same amplitude and
frequency but differ in phase. When such waves superimpose on each other, the phenomenon of
interference takes place. For producing interference, the sources of waves must be coherent. That
is, they must emit waves with zero or constant difference of phase. In this unit you will study how
coherent sources are produced, and how intensity varies in an interference pattern. You will also
learn about the appearance of colours in thin films of oil spread over water.
The phenomenon of diffraction which results due to the superposition of many waves of same
amplitude and frequency, but differing slightly in phase, is usually referred to as the bending of
waves round the corners. Because of this phenomenon, we are at times inclined to think as if
waves do not travel in straight lines. There are two classes of diffraction patterns, called Fresnel
and Fraunhofer classes of diffraction. You will learn that the distinction between these two types
of diffraction is related to the relative separations between the sources of waves and the obstacles
(or the apertures) producing the diffraction patterns.
Both interference and diffraction are very important phenomena in physics. They have
contributed immensely in justifying the wave nature of light. The difference between the two is
quite subtle. Interference arises because of superposition of waves originating from two (or more)
narrow sources, derived from the same source. Diffraction arises from superposition of wavelets
from different numerous parts of the same wavefront, as will be discussed later in this unit.
In Unit 6, you have studied different kinds of waves like sound and light waves. You have also
studied that sound waves are longitudinal while light waves are transverse. Both give rise to the
same phenomena when waves superpose on each other. Basically whatever is true for one kind of
wave is also true for the other; if light waves show the phenomenon, of interference and
diffraction, so do sound waves. Light wave effects have to be observed while sound wave effects
have to be heard. Since the wavelength of sound waves is much greater than the wavelength of
visible region light waves, the sound wave effects are in general on a larger scale compared to the
effects of the light waves.
351
Objectives
After going through this unit, you will be able to:
give examples of coherent sources
derive the condition connecting the path difference between waves from two coherent
sources and the wavelength of the waves used for getting maxima and minima of
intensity on a screen placed in the path of waves
outline the variation of intensity in an interference pattern
explain the principle associated with the appearance of colours in thin films
explain the phenomenon of diffraction
explain the diffraction obtained by a single slit, and
describe the intensity distribution in a diffraction pattern.
9.2 INTERFERENCE
In Unit 8, you have studied the superposition of two waves. You have seen that under certain
special condition superposition leads to the phenomenon of interference. We will now study the
phenomenon of interference in detail. Let us consider the superposition of the following two
waves:
y1 a sin( t kx)
and
y 2 a sin( t kx ) .
These two waves have the same angular frequency and the same wave vector k, and are
travelling along the same direction. They have a phase difference that remains constant with
time. Can you determine the energy distribution after these waves have superposed? If you try,
you will find that this distribution is not uniform in space. The energy is found to be maximum at
certain points and minimum (or probably zero) at others. This type of energy distribution in space
is known as an interference pattern.
In Fig. 9.1, the interference pattern of two waves in a shallow water tank is shown. Here S1 and
S 2 are two sources which produce circular waves on the water surface. The sources S1 and S 2
have to be adjusted in such a way that the waves produced by them on the surface of water are in
constant phase. The resultant of these waves will produce an interference pattern. When the crest
of one wave falls on the crest of another wave, a stronger crest (i.e., one with larger amplitude) is
produced. Similarly, when a trough of one wave falls on the trough of another wave, a shallower
trough is produced. However, when a crest of one wave meets a trough of another wave, their
effects cancel out. Thus it leads to doubling of the amplitude at some points and its reduction to
zero at others. This leads to, what we call, 'an interference pattern.'
352
In most interference experiments performed in the laboratory, the interference pattern is in the
form of fringes. The interference fringes are alternately bright and dark bands, as shown in Fig.
9.2.
A bright band appears on the screen wherever the crest of one wave falls on the crest of another,
or wherever the trough of one falls on the trough of another. Wherever the crest of one falls on
the trough of another, a dark band is produced.
From the above discussion, we can draw the conclusion that for producing an interference pattern,
you require basically two sources. Now, the question arises, what should be the nature of these
sources?
We find from experience that to have a stable and a well defined interference pattern, the two
sources must emit waves either with zero, or with constant difference of phase, say, . If the
sources emit waves with zero or constant difference of phase, they are called coherent sources.
How can we obtain such sources? The easiest way to obtain coherent sources is to obtain them
from the same original source.
One way to obtain such sources in optics is to put an opaque screen containing two slits in the
path of waves emitted by a single source, as shown in Fig. 9.3. The waves originating from the
slits have zero, or a constant difference of phase. When these waves overlap, an interference
pattern is obtained. This you will study in detail in Subsection 9.3.2.
Fig.9.3 Two coherent sources S1 and S 2 obtained from a single source S in an optics experiment
353
In sound, two coherent sources may be obtained by dividing the original longitudinal wave into
two parts as shown in Fig. 9.4. Here one part goes via path I, while the other part goes via path II.
These parts combine again to produce interference. The intensity of sound at any point can easily
be noted, though qualitatively, by listening to the sound around this point.
Fig 9.4 Coherent sources. Path I and II for producing interference of sound waves
Let us now pause for a minute and think as to what will happen if instead of two coherent sources
we use two independent sources. We find from experimentation that if we use any two
independent sources, interference pattern is not produced. This is because with independent
sources, the phase difference between the waves changes rapidly and at random, giving rise to
fast changing interference patterns. Two independent light sources just give a general
illumination on the screen.
SAQ I
Can two small bulbs of 60 W each placed behind two slits form two coherent sources for
interference purpose? If not, why?
354
Fig 9.5: Set-up for observing the interference of wave
The waves from slits S1 and S 2 , separated by a distance d , are in phase with each other.
Whatever phase difference gets created subsequently between the waves from the two sources is
due to their travelling different distances. The slits S1 and S 2 act like coherent sources and
waves of angular frequency and amplitude A. Let us consider a point P at a distance of x i
from S and x 2 from x 2 . Let these distances be sufficiently large compared to d. Let the
displacement at P due to the waves from S, be
It is clear from Eqs (9.1) and (9.2) that the path difference (i.e., the difference between the paths
covered) between the two waves at P is given by ( x1 x2 ) . This will lead to a phase difference
of
2
( x2 x1 )
This is because the phase difference is always associated with the path difference according to
the relation
y y1 y 2
355
In the expanded form Eq. (9.4) can be written as
= A[sin t cos kx1 cos t sin kx1 sin t cos kx2 cos t sin kx2 ]
= A (cos k1 cos k 2 ) sin t A(sin kx1 sin kx2 ) cos t (9.5)
where, we have used Eq. (9.3) for the phase difference . Dividing Eq. (9.7) by (9.6), and
expressing the sum of sine and cosine terms as produced, we get
356
which is the maximum intensity and may be denoted by I max .
Now we calculate the position of maxima. Let d be the distance between the centres of the two
slits, be the angle at which we observe the beams and ( x2 x1 ) the path difference between
the two waves.
Then from Fig. 9.6,
x 2 x1
sin or d sin ( x2 x1 )
d
Maxima in intensity are obtained whenever this path difference is an integral multiple of , the
wavelength of the waves used. Thus for maxima,
Interference minima occur whenever this path difference becomes an odd integral multiple of
/ 2 , i.e.,
SAQ 2
Light passes through two narrow slits with d = 0.8mm. On a screen 1.6m away the distance of the
second order maximum from the axis is 2.5 mm. What is the wavelength of the light used?
In our discussion so far, we have assumed that the two waves have equal amplitudes. If, however,
we assume they have amplitudes a1 and a 2 respectively, then one can show that the resultant
amplitude A1 will be given by
In this case the maxima and minima of intensities (as will be shown in Subsection 9.2.3) become:
I max (a1 a 2 ) 2
(9.16)
I min (a1 a 2 ) 2
357
with their ratio as
I max (a1 a 2 ) 2
(9.17)
I min (a1 a 2 ) 2
To study the variation of intensity with phase difference, , let us plot a graph of I versus .
This is shown in Fig. 9.7. When the phase difference is 0, 2 , 4 , etc., cos = 1. We then have
maximum of intensity, i.e.
On the other hand, whenever the phase differences is , 3 , 5 , etc., cos 1 . We then
have minimum of intensity, i.e.
For a particular case, when the amplitudes are equal i.e., a1 a 2 = a , the intensity varies from a
maximum of 4a 2 to a minimum of zero. In the case of light waves from a monochromatic
source, one would observe the dark fringes of zero intensity separated by bright fringes.
The intensity distribution curve shows that when the two waves arrive at a point on the screen
(exactly) out of phase, they interfere destructively and the resulting intensity (or energy flux) is
zero. Whatever amount of energy is lost from a dip in the zero intensity is, by energy
conservation, found to be redistributed in the maximum intensity peak.
We have seen earlier that for waves of equal amplitude, the intensity can be written as
358
I 2a 2 (1 cos ) 4a 2 cos 2 / 2
Since the average value of cos 2 / 2 is 1/2, the dotted line at I 2a 2 in Fig. 9.7 is the average
intensity, which is actually the sum of the separate intensities from each slit.
It can be seen from Fig. 9.5 that the path difference between the waves reaching P from S1 and
S 2 = ( x2 x1 ) d sin . If is very small, and is measured in radians we can use the
approximation
sin tan
y0
d sin d n
D
giving
nD
yn , where n = 0, 1, 2, …
d
where y n is the distance of the n th maxima from the point where the perpendicular bisector or
the line joining the two slits meets the screen.
nD
yn
d
and
(n 1) D
y n 1
d
the separation y between any two consecutive maxima (or the fringe widths B ) is
D
B y n 1 y n
d
This shows that as long as is very small, the separation between the two consecutive maxima
of intensity is independent of n, i.e., the maxima are evenly spaced. Similarly, it can be shown
D
that the separation between two adjacent minima is also equal to , and that they too are
d
equally spaced.
359
SAQ 3
(a) If the two sources S1 and S 2 in Fig. 9.3 emit waves (i) in phase, and (ii) out of phase by
, discuss the intensity of resultant wave along the perpendicular bisector of S1 and S 2 .
(b) If two waves of amplitude ratio 5:1 interfere, deduce the ratio of intensities at maxima
and minima.
SAQ 4
Two loudspeakers connected to a common amplifier are 5m apart. As one walks along a straight
path 100m away from the speakers, at what spatial period does the intensity vary? Assume that
the wavelength of sound waves = 0.3m.
The interference of light also explains the origin of beautiful colours from oil films on water or
soap bubbles. In the next section, we attempt to discuss these in brief.
Consider a ray of light AB incident on a thin film of uniform thickness t and refractive index
as shown in Fig. 9.8. A part of this is reflected along BC while the remaining part is refracted
along BD into the film. At D it is again partly reflected along DE. The ray DE partly emerges into
the air along EF, which is parallel to BC. The incident ray thus divides at B into two beams of
different amplitudes, out of which the refracted beam suffers multiple reflections at D, E, etc. EH
is perpendicular from E on BC.
The path difference between the rays BC and EF in the reflected system is
( BD DE) BH
Fig. 9.8 Interference in thin films BD is extended to I so that BD = DI. The path difference
between the reflected ray BC and EF is ( BD DE) GI 2 cos r
360
Path difference = 2 t cos r
where r is the angle of refraction in the film. We have already learnt that a phase change of
takes place on reflection at a denser medium. This is equivalent to a path difference of / 2 . The
ray BC is due to reflection at a denser medium. Hence the net path difference between the
reflected rays BC and EF is given by:
2 t cos r / 2 = n (9.22)
We thus see that with monochromatic light alternate bright and dark fringes are obtained. With
white light, which is a mixture of several colours, coloured fringes are obtained.
We have seen above that the path difference depends on and r, apart from t and . Path
difference is different for different colours as is different for different colours. Similarly, for
different angles of incidence, the angles of refraction r are different. Viewing it from different
direction shows different colours. These all lead to the appearance of colours in thin films. These
arise because of the interference of light.
9.3 DIFFRACTION
It is experimentally observed that when a beam of light passes through a small opening (a small
circular hole or a narrow slit) it spreads to some extent into the region of the geometrical shadow.
This is known as diffraction of light.
Consider a point source of monochromatic light S as shown in Fig. 9.9. Place an obstacle, say a
penny or a sharp razor blade, halfway between the source and the screen. Following the rules of
geometrical optics, we expect to see a well-defined and distinct shadow as shown in Fig. 9.9.
Now you carefully examine the shadow. If the experiment is performed in a dark room and the
wavelength of the light used is of the order of the size of the edges of the obstacle, you will find
that the edges of the shadow are not sharp. Inside the shadow, near the edges, the intensity of
light gradually decreases. Outside the shadow it gradually increases, forming alternatively bright
and dark fringes as shown for a penny in Fig. 9.10.
361
Fig 9.9 Diffraction due to a rectangular block
Consider water waves in a ripple tank. Suppose you generate plane waves on the surface of water
in a ripple by giving a periodic up and down motion to a straight vibrator such as a ruler.
Consider an obstacle, such as a slit AB placed in the path of the waves which are travelling as
shown in Fig. 9.11a. As long as the opening AB is large, the plane waves passing through it
appear nearly plane waves. The edges of the emergent plane waves roughly correspond to the
edges of the slit AB. However, when the width of the slit is made narrow, say comparable to the
wavelength of the water waves, then the parallel plane waves entering into the small opening
spread out in the form of approximately circular concentric arcs as shown in Fig. 9.11b.
Fig 9.11 Diffraction pattern obtained by the bending of waves round corners
362
Consequently, these circular waves not only travel in the straight direction but also around the
boundaries of the opening. This phenomenon of bending of waves around the edges of the
aperture (or of the obstacle) is known as diffraction.
This property of bending of waves is distinctly observed when the wavelength of the waves
involved is comparable to the size of the opening (or of the obstacle) through which the waves
pass. For sound waves of frequency 500 Hz the wavelength in air is 0.6m, while that for yellow
light it is 6 10 7 m. Clearly, in the case of sound if the entrance to the room is about 1 m, we
may observe diffraction of sound. However, in the case of light, the dimension of the opening
must be of the order of 10 6 m for a similar diffraction of light to be observed. Thus it is far more
difficult to observe diffraction of light than of water waves or sound waves. Anyway diffraction
of light can be observed in specially devised experiments and we will discuss these now.
Whatever may be the class of diffraction, the resultant distribution of energy in space, or on a
screen, is obtained due to the superposition of waves from different parts of the same wave front.
In the Fraunhofer class, the wavefront considered is plane, while it is spherical or cylindrical in
the Fresnel class. In interference, we have two or more wave sources; while in diffraction, we
have many, almost tending to infinity.
In the discussion to follow, we consider two cases in optics, one of diffraction due to a narrow
single slit, and the other of diffraction at a straight edge. The former belongs to the Fraunhofer
class, while the latter to the Fresnel class.
Fig. 9.12 represents an enlarged diagram of a narrow slit of width, d. Let us assume that as the
plane wavefront reaches the slit, all points in it emit the secondary wavelets in the same phase.
Thus if the disturbance is observed at point P on the far side of the slit at angle to the normal,
then there is a net path difference of d sin between the waves from the two edges of the slit
AB. According to Eq. (9.3), this corresponds to a phase difference of 2 d sin / .
363
Fig 9.12a Diffraction due to a single slit. Note that light is not travelling in a straight line.
sin
Fig 9.12b Graph of sin versus
Now imagine the slit AB is divided into a large number of strips of equal width, d . Each of
these strips sends secondary wavelets and the path difference between the waves arriving at a
point P from two adjacent strips is equal to d sin . The corresponding phase difference is
given by
2 d sin
(9.24)
If we divide the slit AS into a total number of N strips, then clearly d Nd , and the total phase
difference
2 d sin 2
N d sin N (9.25)
364
Suppose that the amplitude of the secondary wave from each strip is denoted by As . Then the
resultant disturbance at P is obtained by the superposition of the waves from all these strips. In
other words,
where 2 r / is the phase difference corresponding to the distance r from the first slit to
point P.
You would recall that we have already discussed this problem of superposition in detail in Unit 2.
It was shown there that the amplitude A of the resultant is obtained by the vector sum of N vectors
of length As , each of which makes an angle with its adjacent vector (see Fig. 9.13).
Fig 9.13 Vectorial addition of contributions from adjacent sources lo give the resultant amplitude.
sin( N / 2)
A As (9.27)
sin( / 2)
We must, however, remember that this subdivision of a slit into a finite number of sub-slits (or
strips) is artificial. We take the limit as N , d 0 . In this case we have a continuous
variation of phase. The vector diagram given in Fig. 9.13 then becomes a smooth circular arc of
radius R. The resultant amplitude is then
sin
A A0 (9.28)
d sin
with A0 RN , and .
From Eq. (9.28), the intensity I 0 of light at any angle with respect to the incident direction is
given by:
365
sin 2
I I0 (9.29)
2
The plot of intensity on screen is given in Fig. (9.12).
For a single slit Fraunhofer diffraction pattern, the minima of intensity are observed at angles
n from the incident direction, where
n d sin n (9.30)
Here n = 1, 2, 3, etc. is the number of the diffraction dark band, starting from the central
maximum.
From Eq. (9.29), you may also note that in the limit as 0 , I 0 approaches I 0 . This becomes
sin
the intensity of the central maximum. This is because lim 1.
0
For values of , for which sin 0 , I 0 0 . This gives us the positions of various minima
which appear for values of equal to n . The corresponding values of n may be calculated
using the earlier given relation, i.e.,
d sin n /
For finding the positions of various maxima lying in-between the different minima, we have to
sin
differentiate the function with respect to and equate that to zero. An elaborate
calculation shows that these maxima, also called the secondary maxima, appear at values of
1.429 , 2.459 , etc. (The details of this calculation will be given in our course on optics).
The heights of these secondary maxima are 1/21, 1/61, 1/120 respectively of the central maxima.
This gives us an idea of the intensity distribution in a single slit diffraction pattern which is
shown in Fig. 9.12b.
sin
d
This shows that as the wavelength increases, or the width of the slit decreases, the angular
spread increases. That is, the narrower the slit, the wider the diffraction pattern. Similarly, the
greater the wavelength, the more widely spread is the pattern.
In terms of the distance D between the slit and the screen the width of the central maximum, y ,
on the screen is given by
366
D
y (9.31)
d
The central peak in the intensity curve is called the primary maximum, while the other peaks are
called secondary maxima. The height of the primary maximum is much more than any of the
secondary maxima.
SAQ5
Calculate the angular spread of the central maximum for light of wavelength 6000Å when the
width of the slit is (i) 10 2 m and (ii) 2 10 5 m.
Fig 9.14 Diffraction due to an obstacle placed in the path of tight waves. Note that there is light
intensity in the Geometrical shadow region, showing that light is not travelling in a straight line.
The intensity curve of a diffraction pattern is always quite different from that of an interference
pattern. In the latter, the heights and widths of peaks are always equal (see Fig. 9.7) meaning
thereby that all the maxima (or minima) are of the same intensity, and are equally spaced. This,
however, is not the case in a diffraction pattern. In a diffraction pattern, the maxima (or minima)
are not of same intensity, and are not equally spaced (see Figs. 9.12 and 9.14).
SAQ6
Sound waves of frequency 1650 Hz fall normally on an opening of width 0.6m. A listener walks
parallel to the opening at a distance of 3 m, starting from a point on the perpendicular bisector of
the opening. Find the positions at which he will observe a minima of sound.
Take the speed of second in air to be 330 m s 1 .
367
9.4 SUMMARY
1. Two sources are said to be coherent if they emit waves with no or constant difference of
phase.
2. As a result of the superposition of waves from two coherent sources, the distribution of
energy in space is not uniform. It is found to alternately pass through maxima and
minima. Such a distribution of energy is called an interference pattern.
3. If two waves of the same frequency and of amplitudes a1 and a 2 , and phases 1 and 2
are acting on a particle, then according to the superposition principle, the amplitude A of
the resultant wave is represented by
A2 a12 a22 2a1a2 cos(1 2 )
4. The phase difference between two interfering waves (from two coherent sources), if they
travel in different paths, is calculated by using the relation:
phase difference = 2 / (path difference)
Maxima in intensity are observed where the path difference is an integral multiple of ,
the wavelength of light used.
5. The distance between any two adjacent maxima or minima in an interference pattern is
given by:
D
d
where is called the fringe width, is the wavelength of the light used, d is the
distance between the two coherent sources, and D is the distance between the sources and
the screen.
6. Diffraction refers to the bending of waves around corners. There are two classes of
diffraction patterns, named as the Fraunhofer and the Fresnel diffractions.
7. Fresnel diffraction phenomena are observed when the source and the screen for observing
the diffraction pattern are at a finite distance from the diffracting aperture or the obstacle.
8. In the Fraunhofer diffraction the source and the screen are at infinite distance from the
aperture causing the diffraction.
9. In a single slit diffraction pattern, the minima in intensity are observed at angles n given
by: n d sin n .
2. Young's experiment is performed with the light of the green mercury line. If the fringes
are measured with a micrometer eye-piece 80 cm behind the double slit, it is found that
20 of them occupy a distance of 10.92 mm. Find the distance between two slits, given
that the wavelength of green mercury line is 5460 Å.
368
3. Light of wavelength 5000Å is incident normally on a slit. The first minimum of the
diffraction pattern is observed to lie at a distance of 5 mm from the central maximum on
the screen placed at a distance of 2 m from the slit. Calculate the width of the slit.
9.6 SOLUTIONS
SAQs
SAQ 1
The light waves emitted by these bulbs are neither in the same phase, nor are they with a constant
difference of phase.
SAQ 2
For small angles, sin can be set equal to y/D, where y is the distance of a given maximum from
the axis and D is the distance from the slits to the screen. For a maximum of second order, we can
write:
d sin n with n = 2
which gives:
y
d n
D
dy
or
nD
where y = 2.5 mm, D = 1.6 m and d = 0.8 mm
(0.8mm)(2.5mm)(1.25 10 3 mm)
1.25 10 3 mm
1.8m
SAQ 3
(a) Path difference between waves along the perpendicular bisector is zero. If sources emit
waves in phase, there will be maximum intensity along the perpendicular bisector. In the
second case, there will be minimum of intensity along the perpendicular bisector.
a12 25a 22
(b) 5
a 22 a 22
I max (5 1) 2 36
2.25
I min (5 1) 2 16
SAQ 4
D 0.3 100m 2
y
d 5.0 m
= 6.0 m
369
SAQ 5
(i) Using Eq. 9.30
6000 10 10
sin 2
6000 10 8 m
d 10 m
SAQ 6
speed 330
0.2 m and
frequency 1650
d 0.6 m
Positions of minima appear where
370
d sin n , with n = 1, 2, 3, …
n
Positions of minima thus lie along the directions n sin 1 , i.e.,
d
0.2 1
1 sin 1 sin 1 19 0
0.6 3
2 0.2 2
2 sin 1 sin 1 42 0
0.6 3
The observer will find minima at
x1 D tan1 3 tan19 0 3 0.344
= 1.03 m
x2 D tan 2 3 tan 42 0 3 0.9
= 2.7 m
(Note that the positions of minima are not equidistant.)
TQs
1. Let be the path length between waves for the 3 rd dark fringe. Then,
(2n 1) / 2 with n 2
5896Å = 5896 10 10 m
(5 5896 10 10 ) / 2 1.4740 10 6 m
3. The angles of diffraction for minimum intensity due to Fraunhofer diffraction at a single
slit are given by:
d sin n with n = 1, 2, 3, …
For the first minimum, n = 1. We can thus write
d sin
If is small, then sin ( is in radians)
d or / d radians
Here, = 5000Å = 5000 10 8 cm and d = ?
5000 10 8 / d radians (A)
The distance between the first minimum and the central maximum is 0.5 cm, and the
distance of the screen from the slit is 2 m.
i.e., 200 cm. This gives
0.5 / 200 radians (B)
Equating Eqs. (A) and (B), we get
371
5000 10 8 cm 0.5 cm
d 200 cm
or
d 5000 10 8 cm 200 / 0.5
= 0.02 cm
372