next up previous
Next: McCullough's Formula Up: Gravitational Potential Theory Previous: Potential Outside a Uniform


Rotational Flattening

Let us consider the equilibrium configuration of a self-gravitating spheroid, composed of uniform density incompressible fluid, which is rotating steadily about some fixed axis. Let $M$ be the total mass, $a$ the mean radius, $\epsilon$ the ellipticity, and ${\mit\Omega}$ the angular velocity of rotation. Furthermore, let the axis of rotation coincide with the axis of symmetry, which is assumed to run along the $z$-axis.

Let us transform to a non-inertial frame of reference which co-rotates with the spheroid about the $z$-axis, and in which the spheroid consequently appears to be stationary. From Chapter 7, the problem is now analogous to that of a non-rotating spheroid, except that the surface acceleration is written ${\bf g} = {\bf g}_g + {\bf g}_c$, where ${\bf g}_g=-\nabla\Phi(r,\theta)$ is the gravitational acceleration, and ${\bf g}_c$ the centrifugal acceleration. The latter acceleration is of magnitude $r\,\sin\theta\,\Omega^2$, and is everywhere directed away from the axis of rotation (see Figure 40 and Chapter 7). The acceleration thus has components

$\displaystyle {\bf g}_c = r\,\Omega^2\,\sin\theta\,(\sin\theta,\,\cos\theta,\,0)$     (912)

in spherical polar coordinates. It follows that ${\bf g}_c = - \nabla\chi$, where
\begin{displaymath}
\chi(r,\theta) = - \frac{\Omega^2\,r^2}{2}\,\sin^2\theta= -\frac{\Omega^2\,r^2}{3}\,\left[1-P_2(\cos\theta)\right]
\end{displaymath} (913)

can be thought of as a sort of centrifugal potential. Hence, the total surface acceleration is
\begin{displaymath}
{\bf g} = - \nabla(\Phi+\chi).
\end{displaymath} (914)

As before, the criterion for an equilibrium state is that the surface lie at a constant total potential, so as to eliminate tangential surface forces which cannot be balanced by internal pressure. Hence, assuming that the surface satisfies Equation (901), the equilibrium configuration is specified by

\begin{displaymath}
\Phi(a_\theta,\theta)+\chi(a_\theta,\theta) = c,
\end{displaymath} (915)

where $c$ is a constant. It follows from Equations (911) and (913) that, to first-order in $\epsilon$ and ${\mit\Omega}^2$,
\begin{displaymath}
- \frac{G\,M}{a} \left[1+ \frac{4}{15}\,\epsilon\,P_2(\cos\t...
...c{\mit\Omega^2\,a^2}{3}\left[1-P_2(\cos\theta)\right]\simeq c,
\end{displaymath} (916)

which yields
\begin{displaymath}
\epsilon = \frac{5}{4}\,\frac{{\mit\Omega}^2\,a^3}{G\,M}.
\end{displaymath} (917)

We conclude, from the above expression, that the equilibrium configuration of a (relatively slowly) rotating self-gravitating mass distribution is an oblate spheroid: i.e., a sphere which is slightly flattened along its axis of rotation. The degree of flattening is proportional to the square of the rotation rate. Now, from (901), the mean radius of the spheroid is $a$, the radius at the poles (i.e., along the axis of rotation) is $a_p=a\,(1-2\,\epsilon/3)$, and the radius at the equator (i.e., perpendicular to the axis of rotation) is $a_e = a\,(1+\epsilon/3)$--see Figure 40. Hence, the degree of rotational flattening can be written

\begin{displaymath}
\frac{a_e-a_p}{a} = \epsilon = \frac{5}{4}\,\frac{{\mit\Omega}^2\,a^3}{G\,M}.
\end{displaymath} (918)

Figure 40: Rotational flattening.
\begin{figure}
\epsfysize =2.25in
\centerline{\epsffile{Chapter12/fig12.02.eps}}
\end{figure}

Now, for the Earth, $a=6.37\times 10^6\,{\rm m}$, $\Omega = 7.27\times 10^{-5}\,{\rm rad./s}$, and $M = 5.97\times 10^{24}\,{\rm kg}$. Thus, we predict that

\begin{displaymath}
\epsilon = 0.00429,
\end{displaymath} (919)

corresponding to a difference between equatorial and polar radii of
\begin{displaymath}
\Delta a = a_e-a_p = \epsilon\,a = 27.3\,{\rm km}.
\end{displaymath} (920)

In fact, the observed degree of flattening of the Earth is $\epsilon=0.00335$, corresponding to a difference between equatorial and polar radii of $21.4\,{\rm km}$. The main reason that our analysis has overestimated the degree of rotational flattening of the Earth is that it models the terrestrial interior as a uniform density incompressible fluid. In reality, the Earth's core is much denser than its crust (see Exercise 12.1).

For Jupiter, $a=6.92\times 10^7\,{\rm m}$, $\Omega = 1.76\times 10^{-4}\,{\rm rad./s}$, and $M=1.90\times 10^{27}\,{\rm kg}$. Hence, we predict that

\begin{displaymath}
\epsilon= 0.101.
\end{displaymath} (921)

Note that this degree of flattening is much larger than that of the Earth, due to Jupiter's relatively large radius (about 10 times that of Earth), combined with its relatively short rotation period (about 0.4 days). In fact, the polar flattening of Jupiter is clearly apparent from images of this planet. The observed degree of polar flattening of Jupiter is actually $\epsilon=0.065$. Our estimate of $\epsilon$ is probably slightly too large because Jupiter, which is mostly gaseous, has a mass distribution which is strongly concentrated at its core (see Exercise 12.1).


next up previous
Next: McCullough's Formula Up: Gravitational Potential Theory Previous: Potential Outside a Uniform
Richard Fitzpatrick 2011-03-31