C H R

Gravitational Force and Potential

10.1 INTRODUCTION

In this chapter we shall investigate Newton's universal law of gravitation and its application. We
introduce the concepts of gravitational field intensity (or simply gravitational field) g and grav-
itational potential V. We calculate these quantities for different mass distributions by applying
Newton's law of gravitation. Gauss's laws will be applied to calculate g and V for simple sym-
metrical mass distributions. Finally, gravitational field equations will be introduced, which are
differential equations satisfied by functions such as g and V. These equations provide a more
general procedure of interest.

10.2 NEWTON'S UNIVERSAL LAW OF GRAVITATION

Newton's universal law of gravitation (which is a law of force), together with Newton's laws of
motion, has been applied by physicists to predict and calculate very precisely the motion of the
planets, moons, satellites, and other objects in the universe. In 1666, 23-year-old Isaac Newton
stated the universal law of gravitation in the following form:

Newton's Universal Law of Gravitation. The gravitational force (or interaction) of at-
traction between any two objects in the universe is directly proportional to the product of
their masses and inversely proportional to the square of the distance between them.

Thus the magnitude of the force F between any two objects of masses m, and m; separated
by a distance rtj is given by

(10.1)
1/
387
388 Gravitational Force and Potential Chap. 10

(a)

Figure 10.1 (a) Gravitational forces

between two masses mt and rrij.
(b) Gravitational force on mass m due
(b) to mass M.

where G is the gravitational constant; its presently accepted value is

G = (6.673 0.003) X 10~ n N-m 2 /kg 2 (10.2)

Referring to Fig. 10.1 (a), we may write the law in vector form as

'y mmj
(10.3)
r r r
ij ij u

where F,-,- is the gravitational force by which mass m, is attracted by mass mp rtj = r, r,- is the
distance between the two masses mi and m,-, and F,-,- is the force by which nij is attracted by mass
mt. According to Newton's third law, we have

7
V=-FJ,

(10.4)

From Fig. 10.1(b), mass m is attracted by mass M with a force F; we may write

(10.5)

where the unit vector u r is in the direction from M to m. The minus sign indicates that F is the
force of attraction with its line of action passing through a fixed point on the line joining the two
masses. Thus the force is directed toward the center of mass M, and the gravitational force is a
Sec. 10.2 Newton's Universal Law of Gravitation 389

central force. The preceding equations are applicable to the situation in which the masses may
be considered point masses. This is possible only if the dimensions of the masses are negligible
compared to the distances between them.
Let us consider a point mass m at P attracted by an extended body of mass M, as shown
in Fig. 10.2. To calculate the force onmat^ we must assume that the gravitational field is a lin-
earfield.That is, the force at P may be calculated by the vector addition of the individual forces
produced by the interactions between the point particle m and the large number of particles in
the extended body. The force dF between m and a small element of volume dV of mass dm is

dV = -G^- ur (10.6)

where dm = p(r) dV, p(r) being the density. The force F acting on m due to the extended body
of mass M may be obtained by integrating Eq. (10.6); that is,

F = - G^KdV' (10.7)

where V indicates integration ever the whole volume. If the extended body is a thin shell that
has a surface density or area density a so that dm = a dA, we may write
ma{T)
F--J G
r2 urdA (10.8)

where A indicates the integration over the whole area. If the extended body is a line source with
a linear mass density A so that dm = A dL, we may write
p _ ,mA(r)
(10.9)

Figure 10.2 Gravitational force on

mass m at P due to an extended body of
mass M and volume V.
390 Gravitational Force and Potential Chap. 10

If the extended body is replaced by a large number of discrete masses mb m2, m3, . . . , m,, the
force on mass m may be written as

mm,
(10.101

where u r is the unit vector in the direction along the line joining m, and m.
According to Eq. (10.7), the system offerees acting on different portions of the extended
body due to mass m at P has a resultant force F acting along a line through the mass m. Ac-
cording to Newton's third law, the force acting on m is F, as shown in Fig. 10.3. On this line
of action of F, we locate a point CG at a distance r from matP such that

, mM
F = G (10.11)

Under these conditions, the gravitational force between the body of mass M and the particle of
mass m is equivalent to a single resultant force F acting on M at CG and F acting on m at P.
The extended body behaves as if all its mass is concentrated at CG. The point CG is called the
center of gravity of the body of mass M relative to the point mass m at P. If the position of m at
P changes, so will the position of CG. In general, CG does not coincide with the center of mass
of M; it may not even be on the line joining the center of mass of M with P. The center of grav-
ity will coincide with the center of mass under the following conditions: (1) If the mass m is far
away from M, the gravitational field will be uniform, different parts of the body will be acted on
by the same force, and the center of gravity will coincide with the center of mass; (2) for a sym-
metrical body, such as a uniform sphere, its center of gravity coincides with its center of mass.
We will encounter another complication if the mass m is also an extended body. In such
cases Eqs. (10.6) and (10.7) must be rearranged, which will involve integrals of both m and dm.

Figure 10.3 Center of gravity CG of

an extended body of mass M relative to
M
mass m at point P.
 Sec. 10.3 Gravitational Field and Gravitational Potential 391

10.3 GRAVITATIONAL FIELD AND GRAVITATIONAL POTENTIAL

As stated before, a gravitational force is a central force; that is, it is a purely radial force pass-
ing through a given point, the center of force. Furthermore, the gravitational force is spherically
symmetric; that is, the magnitude of the force depends only on the radial distance from the cen-
ter of the force and not on its direction. We shall show that spherically symmetric central forces
are conservative; hence the sum of the kinetic energy and the potential energy is constant. Con-
versely, if a central force field is conservative, it must also be spherically symmetric. (Note of
caution: A force that is conservative may or may not be both central and spherically symmetric.)
Suppose a particle of mass m is under the action of a spherically symmetric central force
F with its center of force at O, as shown in Fig. 10.4. In this situation, the force F has only a ra-
dial component Fr, which is a function of r only and may be written as

Fr = F(r) (10.12)
The work dW done by the central force F when m undergoes a small displacement ds, as
shown, is
dW = F ds = F ds cos d (10.13)
But ds cos d dr
where dr is the change of the radial distance from O when mass m undergoes a displacement
ds. Thus
dW= Fdr (10.14)
Since the magnitude of the force F depends only on r, the total work done in going from A to B,
as shown in Fig. 10.4, will be

WAB= F(r)dr (10.15)

Figure 10.4 Work done by a central

force F when a mass m is displaced
from point A to point B.
392 Gravitational Force and Potential Chap. 10

Since this integral and hence the work done depend only on the initial and the final values of r
(not on the path itself) the spherically symmetric force must be conservative.
Once we know that the force is conservative, we can proceed to define a potential energy
function U(r) of an object in such a spherically symmetric central force field. Thus, in going
from A to B, the change in potential energy of an object is

AU=UB-UA= - f[r) dr (10.16)

From Eqs. (10.15) and (10.16), we obtain

WAB=-AU = ~(UB-UA) (10.17)

But the work done is also equal to change in kinetic energy; that is,
(10.18)

Thus, if Eis the total energy, Eq. (10.18) yields

KA + UA = KB UB = (10.19)
which is the law of the conservation of energy.
Since the gravitational force is an inverse square law force,

(10.20)

where C is a constant. Substituting this in Eq. (10.16), we get

r
UB-UA=-j f " C?
which on integration gives
1 1
uB-uA = (10.21)

As is usually done , we define UA = 0 when rA - and UB = U(r) where rB = r; thus

we get

U(r) = - (10.22)
r
which states that the potential energy of a particle in a central force field is a function of the dis-
tance r from the force center. The constant C is negative for attractive forces and positive for re-
pulsive forces. Since the gravitational force is attractive and has the general form
GMm C
F(r) = z = -r where C = GMm (10.23)
r r
Sec. 10.3 Gravitational Field and Gravitational Potential 393

the potential energy of m in the field of M at a distance r from M is

GMm
U(r) = (10.24)
r
If M is a continuous mass distribution of arbitrary shape, the potential energy of m at a dis-
tance r is

W)
H (10.25)

To make the preceding three equations independent of m (the test mass), we introduce the con-
cepts of gravitational field and gravitational potential.
The gravitational field intensity, or gravitational field vector, or simply gravitational field,
g, is defined as the force per unit mass exerted on a particle in the gravitational field of mass M.
That is,
F GM
g (10.26)
m r2

or, for an extended body of mass M, we may write

g = (10.27)

where g has the dimensions of force per unit mass, that is, acceleration. The magnitude of this
gravitational acceleration on the surface of Earth is approximately 9.8 m/s2.
Whenever there is a conservative vector field, as is the gravitational force field, we can al-
ways introduce a gravitational potential (which is a scalar quantity) to represent this field, pro-
vided certain conditions are satisfied. The condition required is that the curl of the vector field
g must be zero. Since g is proportional to 1/r2,

curl g = V X g = 0 (10.28)
(as proved in Chapter 6). This condition will also be satisfied if g is equal to the gradient of a
scalar; that is,
g = - grad V = - V V (10.29)

(remembering V X V V = 0), where V is called the gravitational potential and has the dimen-
sions of energy per unit mass. Since g is only r dependent, V will be only r dependent. Substi-
tuting for g from Eq. (10.26) into Eq. (10.29), we get

GM dV
r2 r
dr

which on integration gives

GM
V(r) = (10.30)
r
 394 Gravitational Force and Potential Chap. 10

It is not necessary to have a constant of integration in Eq. (10.30) because we assume that
V(r) > 0 as r > o.
The gravitational potential due to a continuous distribution of mass M may be written as

T)
dV (10.31)
'v
We may summarize this discussion as follows:
Force:

F= -f Gtn^urdV (10.7)
Jv '
Potential energy:
Gmp(r)
dV (10.25)
Jv r
Gravitational field:
Gp(r)
= -f^firrfV (10.27)
Jy
Jv' '
Gravitational potential:
(r)
dV (10.31)
Jv r

Also,

F = mg (10.32a)
U = mV (10.32b)
g = -VV= -grad V (10.32c)
F = -VU = -grad U (10.32d)
Whenever a mass m is placed in the field of M, it is conventional to speak of the potential en-
ergy of mass m even though such potential energy resides in thefieldand not in the mass itself.

10.4 LINES OF FORCE AND EQUIPOTENTIAL SURFACES

The lines offeree and equipotential lines in two dimensions and equipotential surfaces in three
dimensions are very helpful in visualizing a force field. Let us consider a mass M that produces
a gravitationalfieldin the surrounding space and that may be described by the gravitational field
Sec. 10.4 Lines of Force and Equipotential Surfaces 395

Field lines
Potential lines

Figure 10.5(a) Gravitational field

lines (boldface lines) and equipotential
lines due to a sphere of mass M. The
graph shows the relative value of V(r)
versus r.

vector g. We start from an arbitrary point and draw an infinitesimal line element in the direction
of the vector g at that point. At the end of this line element, we draw another line element in the
direction of g at this new point. We continue this process, and when we join these small line el-
ements, we obtain a smooth line or curve called the line of force ox force field line. We can draw
a large number of such lines in the space surrounding a mass, as shown in Fig. 10.5(a). [See also
Fig. 10.5(b).] These lines start from the surface of a mass and extend to infinity. For a single
mass point, the force lines are straight lines (or radial) extending to infinity as shown. This is
not true in all mass configurations and may be very complicated. For example, Fig. 10.6 will
396 Gravitational Force and Potential Chap. 10

Figure 10.5(b)
Below is the graph of the
gravitational field and gravitational N := 50 i:= 1..N R E := 6.368-106
potential versus the distance r from
the center of mass M of Earth
G : = 6.673-10"u M := 5.98-1024 r. := i-RE

Which value, g or V, G-M -G-M

decreases faster with a change i
in distance r and why? r.

y =-9.84 V =-6.266-107

=-0.025 V =-3.133-106 V.
20
10'
=-0.004 V =-1.253-10 6
50

-10
3.333*10

show potential curves resulting from two unequal masses. The force field lines will be perpen-
dicular to potential curves at every point.
This picture of the lines of force may be used to describe the direction and magnitude of
the field vector g. A tangent drawn at any point to the field line gives the direction of the force
field (of F or g) at that point. The density of these lines, the number of lines passing through a
unit volume (the volume being small, but including the point), gives the magnitude of the vec-
tor field g at that point. No two field lines cross each other because g is a single-valued func-
tion; that is, it has only one value at any given point. It may be pointed out that these field lines
have no real existence, but give a vivid picture depicting the properties of the force field.
We now seek to investigate the relation between the force field lines and the gravitational
potential lines. Suppose we know the gravitational potential Vin the space surrounding a mass.
Sec. 10.4 Lines of Force and Equipotential Surfaces 397

Since the gravitational potential Vis defined for each point in space and is a single-valued func-
tion, we may write

V = V(x, y, z) (10.33)

Suppose we join all the points having the same value of gravitational potential Vo. The equation
representing these points is

V = V(x, y, z) = Vo = constant (10.34)

This is the equation of a surface, called an equipotential surface. We can draw a surface for each
different value of Vo, hence resulting in a large number, or a whole family, of equipotential sur-
faces. In a two-dimensional case instead of equipotential surfaces, we get equipotential lines.
Once again, since V(x, y, z) is a single-valued function, no two equipotential surfaces or lines
will cross each other. Suppose we move a mass m from one point to another point on an equipo-
tential line. By definition, no work will be done. This leads us to the conclusion that the lines of
force are everywhere perpendicular {or orthogonal) to the equipotential lines. This is true be-
cause g = W ; it means that g cannot have a component along an equipotential surface be-
cause V is constant. Thus every line of force must be normal to the equipotential surface, as
shown in Fig. 10.5(a). We shall elaborate on this point shortly. Meanwhile, Fig. 10.6 shows
equipotential lines resulting from two masses M{ and M2. The equipotential surfaces in this case
are defined by the equation

M, Af,
V= -G = constant (10.35)

Consider a mass at point P and let it be displaced a distance ds. The change in its poten-
tial energy, which is equal to the work done, is given by

dU = - F ds = -F.ds (10.36)

where Fs is the component of the force in the direction of the displacement ds. Equation (10.36)
may be written as

dU
(10.37)

This equation states that the component of F in any direction is equal to the negative rate of
change ofpotential energy with distance in that direction. The right side of Eq. (10.37) is called
the directional derivative because its value will depend on the direction of ds relative to F. For
example, consider two equipotential energy lines Uo and Uo + At/ or two equipotential lines Vo
and Vo + AK as shown in Fig. 10.7. If we move form P to Q, which is on the same equipoten-
tial line, dUlds will be zero. But if we move from P to Rh R2, or R on a different equipotential
line, dUlds will be different for different paths, such that dUlds > dU/dsh dU/ds2,.... In this
case, dUlds is maximum when ds is the shortest and hence perpendicular to the equipotential
line at that point. The particular direction for which dUlds is maximum in the direction of the
398 Gravitational Force and Potential Chap. 10

Figure 10.6
Below is a graph of the gravitational potential lines (only three are shown drawn)
due to two nearby unequal masses located at the center of the circles (not shown).
The force field lines, not shown, are perpendicular to potential lines at every point.

j :=!.. 10 = 0..200 rl. : = 20-j r2. :=10-j xl. : = i - 80 x2. : = i - 120

-25

-i
-75

0 25 50 75 100 125 150 175 200

line of force, and the maximum magnitude of dUlds is the magnitude of the vector force at that
point. The maximum value of dUlds and its direction is called the gradient of the potential en-
ergy and is equal to the force F; that is,
F = -gradf/ (10.38)
Since F = mg and U = mV, we may write
g = -grad V = - V V (10.39)
 Sec. 10.5 Calculation of Gravitational Force and Gravitational Potential 399

UQ+AU

Figure 10.7 Gradient of the potential energy. The magnitude of the gradient is
A [//As.

10.5 CALCULATION OF GRAVITATIONAL FORCE AND

GRAVITATIONAL POTENTIAL

We shall start by calculating the gravitational force between a uniform spherical shell of mass
M and a point mass m. We shall show that any spherical shell may be treated as a point mass lo-
cated at the center of the shell. Actually, this is true for any uniform spherically symmetric dis-
tribution of matter. In any of these situations, instead of calculating the force (which is a vector
quantity), it is easier to calculate gravitational potential (which is a scalar quantity). Once the
gravitational potential is known, the gravitational force may be calculated from it. We shall elab-
orate on both these procedures.

Spherical Shell

Consider a thin uniform shell of mass M and radius R, as shown in Fig. 10.8. A particle of mass
m is placed outside the shell at point P a distance r (r > R) from the center of the shell. We
divide the shell into a large number of circular rings like the one shown shaded in the figure.
We can calculate the force between one of these rings and mass m and then sum over all the
rings. As shown in the figure, the width of the shaded ring is R dd, while the radius of the ring
is R sin 9. The circumference of the ring is 2vR sin 9, while the area dA of the circular strip or
shaded ring is
dA = (2ITR sin d)R d6 = 2TTR2 sin 9 dd (10.40)

If a is the density per unit area of the material of the shell, then the mass of the whole spherical
shell is
M
M = (4TT7? V , (10.41)
4TTR2

while the mass dM of the shaded ring is

M
dM = adA = O2TTR2 sin 6d0 = sin 6d9 (10.42)
400 Gravitational Force and Potential Chap. 10

Rdd

(a)

Figure 10.8 Gravitational force between a point mass m and spherical shell of
mass M and radius R.

Point Q, or any other point on the shaded ring, is at the same distance s from the point mass m
at P. The force dFt on m due to any small section of this ring, such as at Q, points toward that
section [see Fig. 10.8(b)]. This force can be resolved into transverse component dF{ sin tf>.
which is perpendicular to PO, and another component dFt cos <j>, which is parallel to PO. Due
to the symmetry of the situation, all the transverse components resulting from considering the
whole ring add up to zero, while the force components parallel to PO due to the whole ring add
up to give

dF = ^dF, = % cos (10.43)

or, substituting for dM, we have

GMm
dF = = sin 9 d6 cos <p (10.44)
2s
Sec. 10.5 Calculation of Gravitational Force and Gravitational Potential 401

The force due to the entire shell is

17
GMm sin 0 d8
2 cos (f>
2s
v
cos (b sin 0 dd
or F = GMm (10.45)

From triangle OPQ, using the law of cosines, we obtain

s2 = r2 + R2 - 2rR cos 6 (10.46)

Since r and R are constants, differentiation yields

2s ds = 2rR sin d dO (10.47)

and, similarly, from the same triangle OPQ, we obtain

R2 = s2 + r2 - 2srcos(t)

or cos <f> = (10.48)

2sr

Substituting for sin 6 dd and cos <f> from Eqs. (10.47) and (10.48) into Eq. (10.45) and chang-
ing the limits by using Eq. (10.46) from 0>77tor R -^ r + R, we obtain

GMm f 1 + -7 ) ds (10.49)
F= r- s
Ar2R j r _ R
which on integration yields

GMm
F= (10.50)
r
In vector notation, this may be written as

GMm .
F= ,-n^ forr>R (10.51a)
r

GM ^
and g= ^-u r , forr > R (10.51b)

where u r is the unit radial vector from the origin O. This result indicates that a uniform spheri-
cal shell acts as if the whole mass of the shell were concentrated at the center. A solid uniform
spherical body may be assumed to consist of a large number of concentric shells. Each shell may
be treated as if its mass is concentrated on the center; hence the mass of the whole sphere may
be assumed to be at the center.
402 Gravitational Force and Potential Chap. 10

To calculate the force on a point mass m placed inside the shell, all we must do is change
the lower limit r R to R r and the upper limit r + R to R + r. Integrating Eq. (10.49) with
appropriate limits is
R+r
GMm I R
2 s + = 0
4r R s R-r

Thus F = 0 and g = 0, for r < R (10.52)

It must be kept in mind that this result (for r < R) is true only for a spherical shell and not for a
solid sphere.
Using the relation given by Eq. (10.16), that is,

--J
At/ = - f B F{r) dr = - [ " F dr (10.16)
'A

and Eqs. (10.51a) and (10.52), we can calculate the potential energy to be
f / GMm\
GMm \ , GMm
U = - I I -j- I dr = , for r > R (10.53)

and U = constant = Cl, for r < R

We can evaluate the constant by substituting r = R in Eq. (10.53), that is,

for r < R (10.54)

R

while the gravitational potential V(=U/m) is

GM
V = , for r > R (10.55)

GM
V = constant = C2 = , for r < R (10.56)

The variation is g and Vfor this case is shown in Fig. 10.9.

We can obtain the preceding results by first calculating the potential energy U(r) and then
calculating F(r) from the relation F = dUldr, as shown next.
The potential energy of mass m at P due to the circular ring of mass dM given by
Eq. (10.42) at a distance s (each point of the ring is at the same distance s) is (see Fig. 10.8)

_ GmdM _ GMm sin Odd

s 2 s

while the total potential energy of m at P is

GMm fw sin Odd

U(r) = (10.58)
I Jo s
Sec. 10.5 Calculation of Gravitational Force and Gravitational Potential 403

Figure 10.9
Below is the graph of the variation
in g(r) and V(r) versus r in the case ,-H
G=6.67210 M=5.98-1024 R=6.37106
of a spherical shell.
24
6.672-10" 10 6.37-10
G: = M :=5.98 R :=
i Because of the large variation in the 10,-n .24
10' 10
values of G, M, and R, divide them with
n
appropriate numbers to make the graphs N := 100 n :=1,5..N 6 :=
50
easier to interpret.)
x : = R-cos(9 y :=fR-sinfe
x and y are used for drawing
a spherical surface or a circular
G-M GM n
surface in two dimensions. Vln :=- V2 :--
R n V n :=ifr n <R,V2 n ,Vl n

The expressions give the values

of g and V inside and outside the gi ,= - G M : = if(r n <R,g2 n ,gl n )
sphere.

/ \
min(V) =-6.264

(
min(g) =-0.916 \ o
-X-
* * * * * *
* * * * * *
max(V) =0 -5
\ /

max(g) =0
-10 -5

(a) Explain the variations in the values of g and V for the values of r given above.

(b) Since max(V) = 0 and max(g) = 0, what do the variations in V and g mean?

From the triangle OPQ in Fig. 10.8, we obtain

s2 = r2 + R2 - 2rR cos 8 (10.46)
Differentiating, while keeping in mind that r and R are constants, and rearranging, we get
sin 6 dd _ ds
s Rr
404 Gravitational Force and Potential Chap. 10

Substituting in Eq. (10.58) yields

GMm
U(r) = - I ds (10.59
2Rr
The limits of integration will depend on the position of the point mass m, as discussed next.

Case (i) r>R: That is, the point mass m at P is outside the shell. As before, the lim-
its 0 > 77 change to smin = r /? smm = r + R. Thus

GMm
ds=-^2R
/::
GMm
gives U(r) = for r > R (10.60)

That is, the potential energy varies as l/r, while

dU d GMm
F= -
dr dr
GMm
gives F= forr> R (10.61)

We may also write

GM
V(r) = - for; (10.62)
r
GM
and for r > R (10.63)
r2

Case (ii) r < R: That is, the point mass m at P is inside the shell. Hence the limits of
integration 0 > TT change to smin = R r > smax = R + r. Thus
.R+r
GMm GMm
U(r) = - ds = 2r
2Rr
GMm
gives U(r) = for r < R (10.64)
R
That is, the potential inside the shell is constant, while

dU d I GMm
F= - = 0
dr dr\ R
F = 0, for r < R (10.65)
Sec. 10.5 Calculation of Gravitational Force and Gravitational Potential 405

as expected. We may also write

GM
V(r) = forr<R (10.66)
R '
g(r) = 0, for r < R (10.67)
These results are graphed in Fig. 10.9, as already mentioned.

Solid Sphere

The results derived for a spherical shell may be extended to a solid sphere. The only require-
ment is that the distribution of matter, that is, the density, be spherically symmetric. Further-
more, the problem becomes simple if the density is uniform.

Case (i) r> R: That is, mass m is at r outside a solid sphere of mass M and radius R.
The sphere may be divided into a large number of shells, each behaving as if the mass of the
shell were concentrated at the center. Independent of the variation in density with radial distance
(that is, symmetric but not necessarily uniform), as in the case of a shell, we obtain

GMm
F= for r> R (10.68)
r1 '
GM
for r > R (10.69)
gif) -
GMm
U(r) = for r > R (10.70)
r
GM
V(r) = for r > R (10.71)

The graphs of V(r) and g(r) are shown in Fig. 10.10.

Case (ii) r<R: That is, the mass m is inside a solid sphere of mass M. Once again we
draw spherical shells. All the shells that are outside a sphere of radius r give zero contribution
to the force, while the shells inside r contribute to the force. For convenience, let us assume that
the density is uniform; that is, the sphere is homogeneous. The fraction of the mass contained
within r is
(477-/3)r3p _ r^
(4TT/3)R3P ~ R3

(see Fig. 10.11) where p is the density of the material. Thus the mass concentrated at the center
is Mr'IR3. Hence the force at r is given by
Gm GMm
forr<R (10.72)
406 Gravitational Force and Potential Chap. 10

> Figure 10.10

V G=6.67210" U M=5.98-1024 R=6.3710 6
Below
versus risdue
the to
graph of V(r)
a solid and g(r)
homogeneous
0A
sphere of radius R and mass M. r\ ,_6.672-10" u 10 _ ,_6.37-106
IN.
M
1 -" c5 -no
98 u 24
10" n io 10 6
Before graphing, we divided the G = 6.672 M =5.98 R=6.37
constants by appropriate powers
of 10 to make the graphs easier N :=100 n : =I N Tn :=- 6 :=TC
to interpret. ' 9 50

Graphing y versus x gives a circle. x :-R-cos( yn:= R-sineJ

GM
VI :=- Vn : = if r n <R,V2 n ,Vl n

_-GM _-GM
: = if(r n <R,g2 n ,gl n )
R3

V and g versus r

min(V) =-9.376

min(g) =-0.944
/ \
(

{\ \
max(V) = 0

max(g) = 0

-10 -5 5 10 15 20
x ,r ,r

(a) How do you think the plot for F will differ from these? Explain.

(b) Explain the variations in the values of g and V.

The potential energy U(r) of the mass inside the sphere may be calculated by using
Eq. (10.72). For r < R, we obtain
GMm
Sec. 10.5 Calculation of Gravitational Force and Gravitational Potential 407

Figure 10.11 Fraction of mass of a

sphere of radius r inside a homoge-
neous sphere of radius R and mass M.

that is, U(r) - U(R) = - ) (10.73)

2R'
But at r = R, from Eq. (10.70), we obtain
GMm
U(R) = -
R
Substituting this in Eq. (10.73), we get

U(r) = - for (10.74)

2R
for r R
or < (10.75)

We can calculate U(r) and V(r) at r = 0, that is, at the center

3GMm
/(0) = - atr = 0 (10.76)
2R
3GM
and atr = 0 (10.77)
2R '
The graphs of g(r) and V{r) for r > R and r < Rare shown in Fig. 10.10.

Shell of Finite Thickness

Consider a shell offinitethickness of inner radius /?] and outer radius R2, as shown in Fig. 10.12.
We want to calculate the potential at point P at a distance r from the center of the shell. By de-
finition, the potential V(r) at P is

V(r) = dV (10.78)

where dV is a small volume element at R:

dV = (R sin 9 dcf>)(R dQ) dR = R2 sin 6 dR dd d<$> (10.79)
408 Gravitational Force and Potential Chap. 10

Figure 10.12 Potential and force on a point mass m at P due to a shell of finite
thickness.

Because of the symmetry about the line connecting O with P the azimuthal angle (f> may be
eliminated by integrating over d(f>; that is, J d<f> = 2v. Also, p(R) = p = constant for a homo-
geneous sphere. Thus

V(r) = -2TTGP f R2 dR f dd (10.80)

From triangle OPQ (Fig. 10.12), we obtain

(10.81)
where R and r are constants; hence differentiating gives

2s ds = 2Rrsinddd
sin 6 dd ds
or Rr
Substituting this in Eq. (10.80) yields

V(r) = - (10.82)

From Eq. (10.81), if 6 = 0, smin = r - R, and if 0 = IT, ^raax = r + R. Therefore, for r > R2,

V(r) = -l7rC}^ [ 2 RdR ( ds (10.83)

r JR, K-R

^ P (RlR2'dR

That is, (10.84)

Sec. 10.5 Calculation of Gravitational Force and Gravitational Potential 409

4T7
Since M= 3
(10.85)

we get the following expression for the gravitational potential outside the shell:
GM
V(r) = - for r > R7 (10.86)

Thus the gravitational potential at any point outside a shell or sphere with a spherically sym-
metric mass distribution is independent of the distribution. It behaves as if the whole mass were
located at the center.
For a point inside the shell, changing the limits in Eq. (10.83) (as we did in a previous
case), we get, for r < Ru
2-nGp
V(r) = - RdR ds (10.87)

RdR

which gives
V(r) = ~ 2irGp(R\ - R2) = constant, for r < /?, (10.88)
Thus the potential inside (r < R\) the shell is constant and is independent of the position.
The potential inside (Rx < r < R2) the shell is a little bit tricky to calculate. But an easy
approach to this problem is to change the lower limit by replacing R} by r in Eq. (10.87) for V(r)
for r < Ru and to change the upper limit by replacing R2 by r in Eq. (10.83) for V{r) for r > R2.
Thus, combining the two gives the potential inside the shell:

V(Rl <r<R2)= -2TTGP(R22 - r2) -

That is,
R3
R2) = - - -- - for/?, <r<R, (10.89)

The field intensity vector g can be calculated from the relation g = dVldr for each of the three
regions by using Eqs. (10.86), (10.89), and (10.88). That is,
GM
g(r) = (10.90)

ATTGP (R]
R, < r < R? (10.91)

g(r) = 0, r<Rt (10.92)

The plots of V(r) using Eqs. (10.86), (10.88), and (10.89) and of g(r) using Eqs. (10.90), (10.91),
and (10.92) are shown in Fig. 10.13. Let us make some important observations. The potential
function V(r) plotted in Fig. 10.13 is continuous across the points r = Rx and r = R2, and its gra-
dient dV(r)/dr, which is force g(r), is also continuous, as shown. If the potential function V(r)
410 Gravitational Force and Potential Chap. 10

-V(r)

~g(r)

Figure 10.13 Variation of V(r) and

g(r) versus r for a shell of finite
thickness.

were not continuous, its derivative would be infinite. That is, the force would be infinite, which
has no meaning. Thus the potential function must be continuous for the force to have any phys-
ical meaning. (Note that the derivative of the force function is not continuous.) That the poten-
tial function is continuous may be seen mathematically as follows. In Eq. (10.89), if we substi-
tute r = R2, we get the same result as from Eq. (10.84). Similarly in Eq. (10.89), if we substitute
r = /?j, we get the same result as from Eq. (10.88).

y Example 10.1
Consider a homogeneous circular disk of radius R, thickness t, and average density p (mass per unit vol-
ume). Calculate the gravitational potential and gravitational field intensity at a point outside the disk and
on the axis of symmetry.
Solution
A circular disk of radius R is shown in Fig. Ex. 10.1 (a). Its mass is given by

M = pV = p(iTR2)t (i)
Sec. 10.5 Calculation of Gravitational Force and Gravitational Potential 411

where t is the thickness of the disk. Let us consider a ring of radius r and width dr, as shown in
Fig. Ex. 10.1 (a). As shown in Fig. Ex. 10.1(b), any small element of this ring is at the same distance s from
point P. We shall calculate potential at point P due to this ring. The point P is at the same distance from
all parts of the ring and lies on the axis of the disk. The mass of the ring is

dm = pdV = p(2irrdrt) (ii)

Hence the potential at point P due this ring is

_ G dm _ Glirptrdr
d V =
~ = - (Z2"+ ,2)' (iii)

Refer to two small elements of the ring at A and B in Fig. Ex. 10.1(b). The intensity at P due to these is
given by g, and g2 pointing along the lines PA and PB, respectively. When gj and g2 are resolved, the hor-
izontal components cancel. Because of the symmetry of the situation, all the horizontal components can-
cel and the vertical components add.

(a)

(b)

Figure Ex. 10.1

412 Gravitational Force and Potential Chap. 10

Using the value of dV from Eq. (iii),

we can calculate V by integrating it as V2-jt-pG-t- -dr
shown. Jz -t-r

To calculate g, we have two alternatives: direct integration or the definition g = dV/dz.

rR
-Gdm 2-ji-p-G-t-z-r
dr g=V g=-
dz dz
2
A2
-A/Z + R -t-z
g=2-jt-p-G-t-

2 2
g=2-;rp-G-t- z +R g=2-7tp-G-t-
z2+R2

How will the graphs of V and g versus z look?

EXERCISE 10.1 Repeat the example for the case of a circular ring of radius R and a linear
mass density X, but with the same mass as the disk. Compare the two results.

) Example 10.2
Consider a thin rod of length L and mass M. Calculate the gravitational potential and
gravitational field intensity at a point P that is at a distance r from the center of the rod and
perpendicular to the rod as shown. Make a plot of F versus r.

S2 = J + r2

M
Sec. 10.5 Calculation of Gravitational Force and Gravitational Potential 413

Solution
L
We can apply Eq. (10.31) using M/L as 2
linear density, and replacing the volume -G-M 1
dx
element dV by the length element dx and L fl 2
r by the perpendicular distance from L
the rod. The gravitational potential energy V is 2

Given the value of V, we can V=G-M-

calculate the value of g by
differentiating V. Simplifying the
expression for g gives the value of
F (= mg).
The graph of F versus r is as shown. g=8-G-M

i:=0.. 10 r.

G:=6.672-10"" M:=5 m:=l L:=1O

r.
F : = 8-G-M-m
i

How do you explain

,-26
the variation in the -1.292*10 r-*r-*
0
value of F for
-6.542 10 \

1/
(a)r = 0, -3.333'10~

(b) r very small as -3.097 JO' "

11
compared to L and -1.907 10' -6.667*10

(c) r very large as -1.302 10' n

compared to L? 12 0 2 4 6
-9.436 10"
12
-7.119 10"
Gravitational force field versus r
Graph V and g -5.54- 10"
12

versus r. 12
-4.42- 10"
2
-3.6- 0" min(F) =-6.542*10 " max(F) =0
12
-2.984 10

EXERCISE 10.2 Repeat the example for the case of a cylindrical rod of radius a, length L,
and mass M.
 414 Gravitational Force and Potential Chap. 10

10.6 GAUSS'S LAW

Gauss's law is used extensively in connection with electric fields in electrostatics. Actually,
Gauss's law is applicable to any situation that involves the inverse-square force law; it is fair
to say that Gauss's law is a compact form of the statement of the inverse-square law. Since
gravitational force is an inverse-square force, let us apply Gauss's law and see its usefulness in
calculating gravitational field intensity g in simple situations.
Let us consider a point mass M. The gravitational field g at a distance r from M is given by

GM (10.93)
g(r) = -

Draw a sphere of radius r with point mass M at the center. We define the radially outward di-
rection as positive. A quantity flux <f> of the gravitational field g through the surface of the sphere
is defined as, using Eq. (10.93),

4> = \irr2gr = -ATTGM (10.94)

where gr is the radial component of g and 4irr2 is the surface area of a sphere of radius r. We
shall show that the total flux due to any mass is independent of the distance r.
Let us consider a mass M that is completely enclosed by an arbitrarily shaped surface, as
shown in Fig. 10.14. Such an arbitrary surface is called a Gaussian surface (GS). Let us con-
sider a point P on this surface where the outward normal to the surface makes an angle 6 with
r from M to P. Resolve g into two components, a radial (or normal) component and a transverse
component (component parallel to the surface). It is only the normal component that contributes
to the flux 4>. Since the normal component of g is -(GM/r2) cos d, the flux dcj> through an ele-
ment of area dA is

GM (10.95)
dcj> = cos 6 dA = g JA - n g dA

GS

Figure 10.14 Mass M enclosed by a

Gaussian surface (GS).
Sec. 10.6 Gauss's Law 415

The projection of dA perpendicular to r is dA cos 6 = n dA, while dA cos d/r2 = dfl, where
dfl is the solid angle subtended by dA at M, as shown. Thus Eq. (10.95) may be written as
d4> = -GMdQ. (10.96)

The total flux $> due to g is obtained by integrating Eq. (10.96) in which contributions from all
the solid angle elements are taken into consideration. Remembering that the complete solid
angle is 4TT, we get

<[> = f d<f>= -GM ( dH= -ATTGM (10.97)

This is the same result we obtained using a spherical surface, as in Eq. (10.94).
From Fig. 10.14 and Eqs. (10.95) and (10.96), we may conclude that the flux d(f> is a scalar
product of g and dA; that is,

dcf> = g dA = n g dA (10.98)

where n dA = dA cos 6 = projection of area dA perpendicular to g or r.

Let us extend our discussion to a large number of masses M1? M2, A/ 3 ,. . . , inside an arbi-
trary closed surface. At any point P on this surface, the total gravitational field is

g = gj + g2 + g3 + (10.99)

while the total gravitational flux through the enclosed surface is

<f>=\g-dA=\n-gdA (10.100)

Combining Eqs. (10.99) and (10.100),

<j> = { g , d A + I" g 2 d A + = - 4 i r G ( M l + M2 + )

If M total = M, + M2+ , then

(10.101)

Equation (10.101) is a statement of Gauss's law, and its validity is based on the fact that the force
is an inverse-square law. Once the flux is calculated by using Eq. (10.101), we can use
Eq. (10.100) in simple symmetrical situations to calculate g, as illustrated in Examplel0.3.

/Example 10.3
By using Gauss's law, calculate the gravitational field intensity at a distance x from an infinite plane sheet
having a surface mass density a.
416 Gravitational Force and Potential Chap. 10

Figure Ex. 10.3

Solution
Draw a Gaussian surface in the form of a right circular cylinder, as shown in Fig. Ex. 10.3. The two end
caps of the cylinder each have a surface area A. The total flux through the cylinder is due to the three sur-
face areas (right, left, and curved) and is given by

in-gdA=jg-dA=i g dA + I g dA + ( g dA (i)
J * -'right -'left 'curved

g will always be perpendicular to the sheet. One the curved surface of the cylinder, an element like dA3
(as shown in Fig. Ex. 10.3) will be perpendicular to g; hence there is no flux through the curved surface;
that is,

f g dA = f g dA cos 90 = 0 (ii)
^curved *

At the end caps, g is always antiparallel to the areas as shown; hence

f g dA = f g dA cos 180 = -gA (iii)

'right >

\ g dA = \ g dA cos 180 = -gA (iv)

 Sec. 10.7 Gravitational Field Equations 417

Thus the total flux is

<& = J g dA = - 2gA (v)

J total

where A is the area of each cap, which is also equal to the area of the sheet enclosed.
The mass of the enclosed sheet is Mtota] = aA. Hence the total flux according to Gauss's law is

= -47rGM total = -ATTG(O-A) (vi)

Combining Eqs. (v) and (vi) gives

-2gA = -4TTG(OJ4)

or g = 277-Gcr (vii)

Thus g is independent of the distance from the plane sheet; that is, it is the same everywhere and is di-
rected toward and perpendicular to the sheet.
EXERCISE 10.3 By using Gauss's law, calculate the gravitational field intensity just outside a spherical
shell having a mass M, radius R, and surface density a.

10.7 GRAVITATIONAL FIELD EQUATIONS

We have briefly outlined the procedure for calculating g and cf> by using Gauss's law for sym-
metrical mass distributions and also by the direct application of the inverse-square gravitational
force law. A more general procedure of interest will be to find differential equations that are sat-
isfied by the gravitational field intensity g(r) and the gravitational potential V(r). Thus, from
Eq. (10.29), we know the relation between g and V to be

g(r) = -VV(r) (10.102)

Taking the curl on both sides, and noting that the curl of a gradient is zero, we get

V x g = - V x V V = 0 (10.103)
That is,
V x g(r) = 0 (10.104)

This vector equation is a set of three differential equations giving relations between the com-
ponents {gx, gy, gz) of g; that is,

^ - ^ = 0, ^ - ^ = 0, ^ - ^ = o (10.105)
dy dz dz dx dx dy
These equations are satisfied by any gravitational field. For a unique determination of g, we need
a relation between g and the distribution of mass. This can be achieved as follows.
418 Gravitational Force and Potential Chap. 10

Figure 10.15 Mass m enclosed by a

V"
volume V of surface area A.

As shown in Fig. 10.15, let us consider a mass m enclosed by a volume V having a sur-
face area A. Thus the flux <j> through this area is

<f>= [ [ n g dA (10.106)
A

But 4> is also given by Eq. (10.97); that is,

<f) = -ATTGM

or for a continuous mass distribution

4>= ~\\ f 4 T T G dm= - [ \ UirGp dV (10.107)

V V

where dm = p dV, p being the mass density. Equating the preceding two equations gives

[ {n g dA = - [ \ UirGp dV (10.108)
A V

Gauss's divergence theorem applied to any vector B is (see Chapter 5)

f f fi B dA = - f f [ V B dV (10.109)
A V

Applying this to the left side of Eq. (10.108), we obtain

JJn -gdA = ~\\\V (10.110)

A V

Substituting this in Eq. (10.108) gives

/ / / = - / / /
Problems 419

fff ( V - g + 4vGp)dV =0 (10.111)

Since this holds for any arbitrary volume V , we may write

V g + A-rrGp = 0
or V-g=-477Gp (10.112)
which is the required relation between g and the mass distribution described by the density
p(x, y, z). In Cartesian coordinates, Eq. (10.112) may be written as

(10.113)
dx dy dz
Knowledge of p(x, y, z), using Eqs. (10.105) and (10.113) and the boundary condition that g >
0 as r > , will uniquely determine g.
Substituting g = - W i n Eq. (10.112) yields
V g = V (-VV) = -AirGp
or V2V = AirGp (10.114)
Rewriting in Cartesian coordinates gives
d2V d2V d2V
2W + 2r- + T = 477^0 (10.115)
dx dy dz2 r

This equation is called Poisson's equation and uniquely determines the value of V(r) with the
boundary conditions that V(r) > 0 as r > . The general solution of Eq. (10.114) is

=
V(r) =
\\rr (10.116)

where p(f) is the density of the volume element dV. This is in agreement with the value of V(r)
given earlier by Eq. (10.31).
In short, Newton's theory of gravitation may be completely summarized by a set of three
equations [(10.26), (10.29), (10.114)]; that is,
F
g(r) = , g = - VV, and V2V = -4-nGp (10.117)
m
or, alternatively, Eqs. (10.26), (10.104), and (10.112); that is,

g(r) = , V X g = 0, and V g = -AirGp (10.118)

m

PROBLEMS
10.1. Starting with Kepler's laws of planetary motion and Newton's laws of motion, derive Newton's
universal law of gravitation.
10.2. Consider a planet of mass M and radius R and with a uniform density p. A tunnel is bored through
this planet to connect any two points on its surface. An object of mass m is thrown in this tunnel.
420 Gravitational Force and Potential Chap. 10

(a) Show that the motion of this object is simple harmonic and calculate the time period of this
motion.
(b) Calculate the time period if the planet were Earth.
(c) Calculate the time period if the hole were drilled through the Moon.
10.3. In Fig. 10.3, if there are a large number of forces acting on M, the result of combining them should
be a single force and a torque. Why is there no torque? (Hint: Take an arbitrary point P and show
that the forces pass through it.)
10.4. An object of mass m released at a very large distance from Earth falls toward Earth's center. Cal-
culate the time it will take (a) to reach halfway to Earth's center, and (b) to reach from halfway to
Earth's center. Compare the two time intervals. Assume Earth to be a point mass.
10.5. Using the expression U(r) = GMm Ir for r > R, show that as a mass m is moved from the sur-
face of Earth to a height h the change in the potential energy is = mgh.
10.6. A particle of mass m in a certain force field given by F = k/x3 is moving toward the center of
the force. Calculate the time it will take the particle to move from a point at a distance D from the
center to the center of the force.
10.7. Suppose an object is dropped from a height h(h< RE, where RE is the radius of Earth). Show that
the speed with which it will hit the ground is

10.8. An object has a free fall in the gravitational field of Earth from infinity to Earth's surface, while
another object falls from a height h = RE with constant acceleration g. Show that they both arrive
at Earth's surface with the same speed.
10.9. Draw gravitational field lines and equipotential lines for a thin rod of finite length. What can you
say about the equipotential surfaces of the rod?
10.10. The gravitational potential at any point P due to two masses Mx and M2 is given by (see Fig. P10.10)

GM, GM,
V(r) = -

Suppose Mi = nM2, where n = 2 or 3. Outline a method for drawing equipotential lines, and draw
them for these two cases.

Figure P10.10

10.11. In Problem 10.10, if the masses are not point masses but are spheres of finite sizes, what changes
will take place in the equipotential lines?
10.12. Draw lines of force and equipotential lines due to two masses Mt and M2 when (a) M, = M2, and
(b) Mj > M2.
Problems 421

10.13. Explain the steps necessary to arrive at Eq. (10.89).

10.14. Consider a uniform hemispherical shell of radius R and mass M with its center at z = 0. Let the Z-
axis be its symmetry axis. Calculate the gravitational potential and field intensity at any point on
the Z-axis. Graph V(r) and g{r). How do these compare with those due to a full (solid) shell?
10.15. Consider a uniform solid hemisphere of radius R and mass M with its center at z = 0. Let the Z-
axis be its symmetry axis. Calculate the gravitational potential and field intensity at any point on
the Z-axis and graph the results.
10.16. Consider a planet of radius R{ and mass M that is surrounded by a cloud of mixed gases with an
average density p. Calculate the gravitational potential and the gravitational field intensity in re-
gions I, II, and III (see Fig. P10.16). Graph V(r) and g(r) for different regions.

Figure P10.16
10.17. Consider a sphere of radius R having a variable density given by p = poe~ar for r < R, where a is
a constant. Calculate the gravitational potential and intensity at some point r inside and outside the
sphere.
10.18. In a tunnel at a distance h below Earth's surface, the density of Earth's material is ph. What will be
the change in the time period of a clock pendulum, that is, AT/T, at this depth? Will the clock run
fast or slow? Calculate your answer in terms of h, M, RE, and ph. If we measure ATIT experimen-
tally, can we determine Earth's mass, knowing the other variables?
10.19. Consider a thin cylindrical rod of length L, radius a, and mass M. Calculate gravitational potential
and the gravitationalfieldintensity at a distance r from the center of the rod and in a direction per-
pendicular to the rod. Graph the results.
10.20. Consider a thin rod of length L and mass M. Calculate the gravitational potential and gravitational
field intensity at a point P that is at a distance r (>L) from the center of the rod and making an an-
gle 6 with the rod, as shown in Fig. P10.20. Calculate only to the second order in Llr. Graph the
results and compare them with the results derived in the text.

Figure P10.20
422 Gravitational Force and Potential Chap. 10

10.21. Calculate the gravitational potential and the gravitational field intensity due to a thin circular ring
of mass M and radius R at a point in the plane of the ring. For large distances from the ring, expand
the expression for the potential and find the first-order correction term. Draw the corresponding
graphs.
10.22. Calculate the gravitation potential and the gravitational field intensity due to a thin circular ring of
mass M and radius R for a point P on the Z-axis at right angles to the plane of the disk, as shown
in Fig. P10.22. Assume z > R, and expand the expression for the potential, keeping terms only in
the second order in Rlz. Graph V(z) and g(z).

Figure P10.22

10.23. Calculate the gravitational potential and the gravitational field intensity due to a thin circular ring
of mass M and radius R at a point P at a distance r from the center of the disk and making an angle 6
with the Z-axis, as shown in Fig. P10.23. Assume that r > R, and expand the expression for the po-
tential, keeping only the second-order terms in Rlr. Graph V(r) and g(r).

ZI

Figure P10.23

10.24. Calculate the gravitational potential and the gravitational field intensity due to a thin circular disk
of mass M and radius R for a point in the plane of the ring. For large distances from the ring, ex-
pand the expression for the potential and find the first-order correction term. Graph V(r) and g(r).
Problems 423

10.25. Calculate the gravitational potential and the gravitational field intensity due to a thin circular disk
of mass M and radius R for point P on the Z-axis at right angles to the plane of the disk. Assume
z> R, and expand the expression for the potential, keeping terms only of the second order in Rlz.
Graph V(z) and g(z).
10.26. Calculate the gravitational potential and the gravitational field intensity due to a thin circular disk
of mass M and radius R at a point P at a distance r from the center of the disk and making an an-
gle 6 with the Z-axis, similar to Fig. P10.23. Assume that r > R, and expand the expression for the
potential, keeping only the second-order terms in Rlr. Graph V(r) and g(r).
10.27. Consider a body that has a cylindrical symmetry with density p(r,ff)for r < R and p = 0 for r >
R. Calculate the gravitational potential at a point (r,ff)far away from the body. (Expand in powers
of Rlr.)
10.28. Consider a system of binary stars, each of mass M and separated by a distance 2r. These stars or-
bit about their common center of mass. A mass m is located at a point P, as shown in Fig. P10.28.
(a) Calculate the gravitational potential and gravitational field intensity and force at point P.
(b) Repeat part (a) if z > r, and z<r.
(c) Suppose the mass m is at point O and then is slightly displaced. Show that it executes simple
harmonic motion.
z,,

Figure P10.28
10.29. A and B are two thin concentric shells of radii R{ and R2 and masses Mi and M2, respectively. A
point mass m is located at a distance r from O. Calculate the gravitational potential and field in-
tensity in the three regions shown in Fig. P10.29. If a mass m is released from infinity, what will
be its speed when it reaches 01

Figure P10.29
424 Gravitational Force and Potential Chap. 10

10.30. By applying Gauss's law, calculate the gravitational field and gravitational potential due to a ho-
mogenous sphere of radius R.
10.31. Calculate the gravitational field intensity and gravitational potential at a distance x from an infinite
sheet of surface density (7 in the XY plane.
10.32. A mass m is placed at a depth h in a tunnel in Earth. Show by using Gauss's law that the force ex-
erted on this mass is due to the mass of the spherical portion below the tunnel. What will be the
force on m at one-half the radius of Earth?
10.33. By using Gauss's law, calculate the gravitational field intensity inside and outside an infinitely long
cylindrical shell of radius R and mass M.
10.34. By using Gauss's law, calculate the gravitational field intensity at a distance x from an infinitely
long cylindrical rod of mass density crper unit length.
10.35. Show that if we consider a mass M outside an enclosed surface, then the net flux through this sur-
face is zero.
10.36. Show that the equations V X g = 0, V - g = -AirGp, and V2V = AirGp are all satisfied by the
gravitational field intensity and gravitational potential in Problems 10.30 and 10.31.
10.37. Show that curl g = 0.

SUGGESTIONS FOR FURTHER READING

BARGER, V., and OLSSON, M., Classical Mechanics, Chapter 7. New York: McGraw-Hill Book Co., 1973.
FRENCH, A. P., Newtonian Mechanics, Chapters 8 and 11. New York: W. W. Norton and Co., Inc., 1971.
KITTEL, C , KNIGHT, W. D., and RUDERMAN, M. A., Mechanics, Berkeley Physics Course, Volume 1, Chap-
ter 9. New York: McGraw-Hill Book Co., 1965.
MARION, J. B., Classical Dynamics, 2nd ed., Chapter 2. New York: Academic Press, Inc., 1970.
STEPHENSON, R. J., Mechanics and Properties of Matter, Chapter 4. New York: John Wiley & Sons, Inc..
1962.
SYMON, K. R., Mechanics, 3rd ed., Chapter 6. Reading, Mass.: Addison-Wesley Publishing Co., 1971.