Basic notions of
Celestial Mechanics
Mohamad Magdy
Abd El Rasoul
Egypt-kafr El
Shiekh -Balteem
Don’t worry
about your mathematical difficulties;
I assure you
that mine are greater !
Albert
Einstein
1. Law of
Universal Gravitation
As his friend Halley
was right in stating "Nec fas
est propius mortali attingere divos" - "No
mortal nearer the Gods could reach", we should feel when
studying the great genius of Sir
Isaac Newton.
His great contributions to
science, astronomy, optics, physics, the invention of Calculus and
marvelous geometric perception of the Universe make us remember the Alexander
Pope (XVIII century) poem, where in some place says "Nature
and nature’s laws lay hid
in night. God said 'let Newton be' and all was light".
The law of gravitation formulated
by Newton says "Each particle of
matter attracts any other particle of matter with a force proportional
to the product of their masses and inversely proportional to the square
of the distance between them". In the universal language
of mathematics, this can be stated as:
Where m and m1 are the masses of both particles, r is the distance between them, F is
the attractive gravitational force and G is the so-called constant of
universal gravitation.
The numeric value of
is 6.672 x 10-11 m3 kg-1s-2.
Before going any deeper in the
application of the law of universal gravitation and its relationship
with the planetary orbits, let us see an example of said law in the
circular orbits of the terrestrial artificial satellites
Figure 1.1
Figure 1.1 depicts a mass ms attached to the tip of a string which is being swung around a fixed point
with an angular velocity
. According to the dynamics of our
example, the mass experiences a direct action of the centrifugal force
which is proportional to the mass, its distance to the rotation point
and to the square of the angular velocity.
The resulting centrifugal force
is Fc = msr
2 and it is said to be in equilibrium with
the force made at the point of rotation, should not this be the case the
mass would be shot in a tangential trajectory to the circle of rotation.
Applying this example to a
satellite in rotation around the Earth we have that the force that keeps
the satellite in orbit is the force of gravity, which replaces the
string in our above example, and is proportional to the mass of the
Earth mt, the mass of the
satellite ms and inversely
proportional to the square of the distance between the center of the
Earth and the satellite.
Therefore, this force
Fg is expressed as:
1.1.
If we set this equation equal to
the corresponding centrifugal force, we have:
1.2
Equation 1.2
is the basis of the circular motion of a satellite in orbit around the
Earth. It is important to note that the mass of the satellite, in this
equation, can be cancelled as it appears in both sides and therefore if
we were to calculate the mass of a body, for instance that of the Earth,
the reference central body should be the Earth itself and not the
observed body.
From equation
1.2 we have:
1.3
Let us
introduce the concept of mean
movement n, instead of the angular velocity
. From equation 1.3 we have that said
mean movement is:
1.4
The
instantaneous circular speed of the satellite is obtained by multiplying
each member of the equation by
r.
1.5
Based in the
satellite speed around the Earth or the value of its mean movement we
can calculate the period of the satellite as the longitude of the arc
for a complete revolution is 2
r. Denoting the speed or our satellite
as Vc , the duration or a
revolution is:
1.6
where P is the
period, ie.
1.7
Now, equation
1.5 allows us to deduce that the speed of a satellite can be determined
from the mass of the Earth, its distance and the constant of universal
gravitation.
As a matter of
fact, we can express the distance to the center of the Earth as the sum
of the height of the satellite and the radius of the Earth, as.
1.8
where Rt
is the equatorial Earth’s radius and h is the height or
elevation of the satellite above the surface of the Earth.
Based on the
above, we can transform equation 1.5 in:
1.9
As a practical
application we can study the case of the geostationary satellites, whose
orbits are placed about the plane of the equator and which allow the
transmission and reception of radio communication through very long
distances by means of parabolic antennas. The main parameter of these
satellites is that they always stay on the same spot in the sky, that is
why the satellites dishes are always fixed !. In this sense, their
period with respect to the fixed stars is 23h 56m., a concept closely
related to that of sidereal time.
Because this
period is 23h 56m or 86,160 seconds, their mean movement is:
From equations
1.4 and 1.7 we have:
; where we have replaced r by Rt +
h
From this last
equation we can deduce that the height h of the satellite is:
The numerical
values used are
for G,
for the Earth’s mass and 6378.14 km. for its radius.
The speed of
these satellites is obtained as:
But going back
to our study if we square both member of equation 1.4 and multiply them
by r3, we obtain:
1.10
This equation
is a very famous one in celestial mechanics and its a form of one of the
Kepler´s Laws which will be discussed below. In fact, equation 1.1 is
an approximation to the exact relationship as it does not consider the
mass of the satellite in respect to that of the central body. This
approximation, however, for the case or artificial satellites is
perfectly valid for all practical purposes.
Figure 1.2
Figure 1.2
allows us to derive the mentioned exact relationship if we suppose that
it represents the Earth and a satellite, two asteroids of similar mass
or a binary star system. In this figure r1 and r2
represent de distances of the bodies with masses m1 and m2
to their mutual center or mass.
For the case
of the Earth and an artificial satellite, we have m1= mt
y m2 = ms.
The system is
in static equilibrium if
For the
satellite:
1.11
and for the
Earth:
1.12
If we cancel ms
from equation 1.11 and mt from
equation 1.12, and equal them we have:
1.13
By adding
equations 1.11 and 1.12 we
obtain:
|