Foucult Pendulum

Download as pdf or txt
Download as pdf or txt
You are on page 1of 7

A Description of the Motion of a Foucault

Pendulum

Anthony J. Hibbs
University of Warwick

October 2, 2010

1 Introduction

In 1851, French physicist Jean Leon Foucault designed a revolutionary exper-


iment which demonstrates that the Earth is a rotating body. His apparatus
was rather simple; a 28 kg mass on a 67 m long wire which was attached to
the ceiling of the dome of the Panthon in Paris in such a way that allowed
the pendulum to swing freely in any direction. Foucault found that the plane
of oscillation rotated in a clockwise direction, as viewed from above, at a
rate of approximately 11 degrees per hour, and one full 360 degree rotation
of this plane took 32.7 hours.

2 How Does It Work?

This experiment tells us that the Earths surface is not an inertial frame of
reference, that is a frame which is either at rest or moving with a constant
velocity in a particular direction, with no external forces applied. In rotating
frames of reference, such as on the surface of the Earth, the velocity of the
frame is constantly changing direction. It is this that, in the case of the
pendulum, causes it to appear as though angular momentum is not conserved
and Newtons laws are not obeyed. As observed from a position at rest in

1
space (an inertial frame) however, one would see that angular momentum is
in fact conserved and Newtons laws are obeyed.

One may expect that the best way to treat this problem would be to
analyse the mechanics from an inertial frame of reference. In actual fact it
is somewhat less complicated to instead use a non-inertial rotating frame
fixed to the Earths surface. In order to do this, we must find a way of
translating from the inertial to the non-inertial frames and care must be
taken when considering vectors, in particular their time derivatives, as will
become apparent.

When considering mechanics in non-inertial frames, additional forces


must be considered which vanish in inertial frames. For this reason these
forces are termed fictitious forces and those which we concern ourselves with
here are the centrifugal and Coriolis forces. The details of the centrifugal
force are not essential in understanding why the plane of oscillation rotates,
however its effects are interesting and the description emerges naturally from
the mathematics which follows this introduction.

3 Rotating Frames

In a body rotating with angular velocity , a


particular point with position vector r has a
velocity v given by

v = r. (1)

In fact for any vector x (not necessarily a posi-


tion vector) fixed in a rotating body,

dx
= x. (2)
dt Figure 1: The directions
of the x, y and z axes at
Thus for unit vectors i, j and k, directed along
a point on the Earths sur-
the x (east), y (north) and z (upwards) axes
face with latitude .
fixed on the Earths surface (see figure 1), we
have
di dj dk
= i, = j and = k. (3)
dt dt dt

2
Now consider the position vector of a particle in this coordinate system,

r = rx i + ry j + rz k. (4)

Note that the origin of the coordinate system here (where r = 0) is at the
centre of the Earth. The rate of change of this vector will be measured
differently by an observer in an inertial frame, as compared to an observer
standing on the Earths surface. We shall use V to denote the rate of change
of r as measured by the inertial observer and v and r to denote the rate of
change of r as measured by the non-inertial observer on the Earths surface.
The rates of change of scalar quantities, i. e. the derivatives of rx , ry and
rz , will always be measured to be the same by both observers. According to
the non-inertial observer on the Earths surface, the directions of the unit
vectors are constant in time, so
drx dry drz
v = r = i+ j+ k. (5)
dt dt dt
For the inertial observer, however, the directions of the unit vectors are
varying in time, so we have
drx dry drz di dj dk
V= i+ j+ k + rx + ry + rz
dt dt dt dt dt dt
= r + r, (6)

from equations 3 and 4. Applying the same formula, the rate of change of
V as measured by the inertial observer is
dV
A= = V + V, (7)
dt
where V is the rate of change of V as measured by the non-inertial observer.
From equation 6,

V = r + r + r
and V = r + ( r).

Hence,
A = r + 2 r + ( r) + r. (8)
Newtons second law for a particle moving under the influence of gravity
plus an additional force F is

mA = mg + F, (9)

3
where we must consider the absolute acceleration A = dV/dt, since equation
9 holds only in inertial frames. (If one were to use V in place of A then
equation 9 would not be valid). g is a vector pointing towards the centre
of the Earth, with a magnitude equal to g ( 9.81 ms1 ). Substituting
equation 8 into 9 and rearranging, we obtain

mr = mg + F 2m r m ( r) m r. (10)

Equation 10 is the equation of motion for a particle in a non-inertial rotating


frame. The third and fourth terms on the right hand side are the Coriolis
and centrifugal forces respectively. The fifth term is the Euler force, which
can be neglected since, for the Earth, 0.

4 The Centrifugal Force

Let us delay our discussion of the Coriolis force and first focus on the cen-
trifugal force. This is a position dependent force, whereas the Coriolis force
is velocity dependent, which means that when any measurement is taken of
the acceleration due to gravity, what is actually measured is the effective
acceleration due to gravity,

geff = g ( r). (11)

As one may expect, the direction of ( r) is perpendicular to the


axis of rotation of the Earth and directed outwards, thus it has components
in the y and z directions. For a point with latitude ,

| ( r)| = | r| = 2 r sin(90 ) = 2 r cos . (12)

Thus the y and z components of geff are


(y)
g = 2 r cos sin (13)
eff
(z)
g = g 2 r cos2 , (14)
eff
and the effect of the centrifugal force is to decrease the magnitude of the
acceleration due to gravity and to alter its direction slightly; in the north-
ern hemisphere a plumb line is deflected to the south and in the southern
hemisphere it is deflected to the north.

4
At the equator ( = 0 ), the y component vanishes and the z component
has its maximum value, whereas at the poles ( = 90 ), both components
vanish. Thus the acceleration due to gravity will be measured to be less
at the equator than at the poles. This is compounded by the fact that the
radius of the Earth is less at the poles than at the equator, which incidentally
is due to the centrifugal force. The overall effect is that one weighs less at
the equator than at the poles, although for those readers interested in weight
loss, there are simpler and more cost-effective methods.

5 The Coriolis Force

Before delving into the mathematics of the Coriolis force it is beneficial to


attempt to visualise its effect. Neglecting friction, consider a disk rotating
anticlockwise with a particle being fired across the diameter. Standing sta-
tionary (an inertial frame) and watching the disk rotate in front of you, if
you were to roll a ball across the diameter of the disk, you would see the ball
travel in a straight line through the centre of the disk. If you were rotating
with the disk (a non-inertial frame) however, and rolled the ball towards
the centre in the same manner, you would see the particle being deflected
to your right as if subject to a force acting perpendicular to the particles
velocity. This force is the Coriolis force and a good visual supplement to this
description is located at http://www.youtube.com/watch?v=49JwbrXcPjc.
We turn our attention back to the Earth now and imagine standing still at
the North Pole. Neglecting air resistance, a bird flying south would appear
to be following a slightly curved trajectory to the west, whereas an observer
watching from at rest in space would see the bird flying in a straight line
while the Earth rotates underneath.

With our pendulum, we assume that the distance through which the bob
moves is small enough that the gravitational and centrifugal forces can be
treated as constant and we combine them together into a constant effective
acceleration due to gravity, geff . From equation 10, the equation of motion
now becomes
mr = mgeff + F 2m r. (15)
In our coordinate system with the unit vectors i, j and k, pointing east, north
and upwards respectively, the angular velocity of a point on the Earths

5
surface with latitude is

= cos j + sin k. (16)

Hence we obtain for the Coriolis force

FC = 2m r = 2m[(ry sin rz cos )i rx sin j + rx cos k]. (17)

We concern ourselves primarily with the horizontal motion in the x, y plane.


The vertical component of the Coriolis force is negligible since it provides
a small correction to geff , with the sign alternating on each half-period.
Assuming that the length l of the pendulum is large and that the amplitude
of the oscillations is comparatively small, the velocity in the z-direction can
be neglected (rz 0) and we can utilise the small angle approximation
(sin ) for the angle the pendulum makes with the vertical (its rest
position). Thus from equations 15 and 17, assuming there are no additional
external forces (F = 0), we obtain the equations of motion for the pendulum
in the x and y directions,
geff
rx = rx + 2ry sin (18)
l
geff
ry = ry 2rx sin . (19)
l
The first terms on the right hand sides are the familiar restoring forces for
a pendulum exhibiting simple harmonic motion, and the second terms are
the contributions from the Coriolis force. In order to obtain an explicit
solution
to these equations, we can multiply equation 19 by the imaginary
unit i = 1, and add it to equation 18, giving

+ 2i sin + 02 = 0, (20)
q
where = rx + iry and 0 = geff /l is the natural angular frequency of the
oscillations, i. e. the frequency of simple harmonic oscillations in the absence
of the Coriolis force. Equation 20 resembles the equation of motion for a
damped pendulum, except for the imaginary damping term. We look for
solutions to this equation of the form = Aeit and substituting this into
equation 20 gives
2 + 2 sin 02 = 0
q
= sin 2 sin2 + 02 = sin 1 , (21)

6
q
where 1 = 2 sin2 + 02 0 since since the frequency of the rotation
of the Earth is much less than the frequency of oscillation of the pendulum
(  0 ). Hence the solution is

= ei sin t (Aei1 t + Bei1 t ). (22)

The terms in parentheses describe simple harmonic motion with frequency


1 , and A, B are constants fixed by the initial conditions. We would there-
fore expect that the prefactor ei sin t is the term which will tell us about
the rotation of the plane of oscillation. Dropping the SHO terms and taking
the real part of , we obtain

<[] = rx = cos( sin t). (23)

Defining T as the period of rotation of the plane of oscillation, clearly


rx (t = 0) = rx (t = T ), so

cos(0) = 1 = cos( sin T )


sin T = 2
2
T = . (24)
sin

The value of for the Earth in convenient units is = 2/(24 hours).


Hence we find that the period of rotation of the plane of oscillation for a
Foucault pendulum is
24
T = hours. (25)
sin
From equation 23, we can also find the angle through which the plane of
oscillation rotates in one hour by finding the value of the argument of the
cosine for t = 1 hour, thus
2 360
= sin rad. = sin deg. (26)
24 24
The latitude of Warwick University is 52 22 48 (52 degrees, 22 minutes
and 48 seconds). Inserting this into equations 25 and 26 gives values of
T = 30.3 hours and = 11.9 . The yellow lines on the floor display are
separated by this angle.

You might also like