next up previous
Next: Jacobi Ellipsoids Up: Hydrostatics Previous: Equilibrium of a Rotating

Maclaurin Spheroids

One, fairly obvious, way in which the constraint (2.115) can be satisfied is if $ a_2=a_1$ . In other words, if the planet is rotationally symmetric about its axis of rotation. An ellipsoid that is rotationally symmetric about a principal axis--or, equivalently, an ellipsoid with two equal principal radii--is known as a spheroid. In fact, if $ a_2=a_1$ then the cross-section of the planet's outer boundary in any plane passing though the $ x_3$ -axis is an ellipse of major radius $ a_1$ in the direction perpendicular to the $ x_3$ -axis, and minor radius $ a_3$ in the direction parallel to the $ x_3$ -axis. Here, we are assuming that $ a_1>a_3$ : that is, the planet is flattened along its axis of rotation. The degree of flattening is conveniently measured by the eccentricity,

$\displaystyle e_{13} \equiv (1-a_3^{\,2}/a_1^{\,2})^{1/2}.$ (2.118)

Thus, if $ e_{13}=0$ then there is no flattening, and the planet is consequently spherical, whereas if $ e_{13}\rightarrow 1$ then the flattening is complete, and the planet consequently collapses to a disk in the $ x_1$ -$ x_2$ plane.

Let $ u=a_1^{\,2}\,\lambda$ and $ \lambda=e_{13}^{\,2}/z^{\,2}-1$ . Setting $ a_2=a_1$ in Equation (2.114), we obtain

$\displaystyle \frac{\omega^{\,2}}{2\pi\,G\,\rho}$ $\displaystyle =(1-e_{13}^{\,2})^{1/2}\,e_{13}^{\,2}\int_0^\infty \frac{\lambda\,d\lambda}{(1+\lambda)^2\,(1+\lambda -e_{13}^{\,2})^{3/2}}$    
  $\displaystyle = \frac{2\,(1-e_{13}^{\,2})^{1/2}}{e_{13}^{\,3}}\left[\int_0^{e_{...
...}}-(1-e_{13}^{\,2})\int_0^{e_{13}}\frac{z^{\,2}\,dz}{(1-z^{\,2})^{3/2}}\right].$ (2.119)

Performing the integrals, which are standard (Speigel, Liu, and Lipschutz 1999), we find that

$\displaystyle \frac{\omega^{\,2}}{2\pi\,G\,\rho} = \left(\frac{3-2\,e_{13}^{\,2...
...^{1/2}\,\sin^{-1} e_{13} - \left(\frac{3}{e_{13}^{\,2}}\right)(1-e_{13}^{\,2}).$ (2.120)

This famous result was first obtained by Colin Maclaurin (1698-1746) in 1742. Finally, in order to calculate the potential energy, (2.116), we need to evaluate

$\displaystyle \alpha_0$ $\displaystyle = \frac{1}{a_1}\int_0^\infty\frac{d\lambda}{(1+\lambda)\,(1+\lamb...
...{\,2})^{1/2}} =\frac{2}{a_1\,e_{13}}\int_0^{e_{13}}\frac{dz}{(1-z^{\,2})^{1/2}}$    
  $\displaystyle = \frac{2}{a_1}\,\frac{\sin^{-1} e_{13}}{e_{13}}.$ (2.121)


Table 2.1: Properties of the Maclaurin spheroids.
               
$ e_{13}$ $ \widehat{\omega}$ $ \widehat{L}$ $ -\widehat{E}$ $ e_{13}$ $ \widehat{\omega}$ $ \widehat{L}$ $ -\widehat{E}$
               
0.00 0.00000 0.00000 0.60000 0.60 0.31729 0.18037 0.56233
0.05 0.02582 0.01266 0.59980 0.65 0.34484 0.20286 0.55320
0.10 0.05168 0.02540 0.59919 0.70 0.37239 0.22834 0.54200
0.15 0.07758 0.03830 0.59817 0.75 0.39967 0.25792 0.52800
0.20 0.10357 0.05144 0.59672 0.80 0.42612 0.29345 0.51001
0.25 0.12967 0.06491 0.59479 0.85 0.45046 0.33833 0.48587
0.30 0.15591 0.07882 0.59236 0.90 0.46932 0.39994 0.45107
0.35 0.18231 0.09329 0.58936 0.95 0.47045 0.50074 0.39272
0.40 0.20889 0.10846 0.58572 0.96 0.46472 0.53194 0.37485
0.45 0.23567 0.12450 0.58135 0.97 0.45418 0.57123 0.35273
0.50 0.26267 0.14163 0.57612 0.98 0.43475 0.62486 0.32351
0.55 0.28989 0.16013 0.56986 0.99 0.39389 0.71209 0.27916


Let $ e_{13}=\sin\gamma$ . Thus, $ \gamma =0$ corresponds to no rotational flattening, and $ \gamma=\pi/2$ to complete flattening. Moreover, $ a_1=a_0\,(\cos\gamma)^{-1/3}$ and $ a_3=a_0\,(\cos\gamma)^{2/3}$ , where $ a_0=(a_1\,a_2\,a_3)^{1/3}= (3\,V/4\pi)^{1/3}$ is the mean radius. It is also helpful to define $ \widehat{\omega}=\omega/(2\pi\,G\,\rho)^{1/2}$ , $ \widehat{L}=L/(G\,M^{\,3}\,a_0)^{1/2}$ , and $ \widehat{E}= E/(G\,M^{\,2}/a_0)$ . The previous analysis leads to the following set of equations that specify the properties of the so-called Maclaurin spheroids:

$\displaystyle \widehat{\omega}^{\,2}$ $\displaystyle =\frac{\cos\gamma}{\sin^2\gamma}\left[(1+2\,\cos^2\gamma)\,\frac{\gamma}{\sin\gamma}-3\,\cos\gamma\right],$ (2.122)
$\displaystyle \widehat{L}$ $\displaystyle = \frac{\sqrt{6}}{5}\,\frac{(\cos\gamma)^{-1/6}}{\sin\gamma}\left[(1+2\,\cos^2\gamma)\,\frac{\gamma}{\sin\gamma}-3\,\cos\gamma\right]^{1/2},$ (2.123)
$\displaystyle \widehat{E}$ $\displaystyle = -\frac{3}{10}\,\frac{(\cos\gamma)^{1/3}}{\sin^2\gamma}\left[(1-4\,\cos^2\gamma)\,\frac{\gamma}{\sin\gamma}+3\,\cos\gamma\right].$ (2.124)

These properties are set out in Table 2.1.

Figure: Normalized angular velocity squared of a Maclaurin spheroid (solid) and a Jacobi ellipsoid (dashed) versus the eccentricity $ e_{13}$ in the $ x_1$ -$ x_3$ plane.
\begin{figure}
\epsfysize =3in
\centerline{\epsffile{Chapter02/fig5.01.eps}}
\end{figure}

In the limit, $ \gamma\rightarrow 0$ , in which the planet is relatively slowly rotating (i.e., $ \widehat{\omega}\ll 1$ ), and its degree of flattening consequently slight, Equations (2.122)-(2.124) reduce to

$\displaystyle e_{13}$ $\displaystyle \simeq \frac{\sqrt{15}}{2}\,\widehat{\omega},$ (2.125)
$\displaystyle \widehat{L}$ $\displaystyle \simeq \frac{\sqrt{6}}{5}\,\widehat{\omega},$ (2.126)
$\displaystyle \widehat{E}$ $\displaystyle \simeq -\frac{3}{5}.$ (2.127)

In other words, in the limit of relatively slow rotation, when the planet is almost spherical, its eccentricity becomes directly proportional to its angular velocity. In this case, it is more conventional to parameterize angular velocity in terms of

$\displaystyle m = \frac{\omega^{\,2}\,a_0}{g_0} =\frac{3}{2}\,\widehat{\omega}^{\,2},$ (2.128)

where $ g_0=G\,M^{\,2}/a_0^{\,2}$ is the mean surface gravitational acceleration. Furthermore, the degree of rotational flattening is more conveniently expressed in terms of the ellipticity,

$\displaystyle \epsilon = \frac{a_1-a_3}{a_0}\simeq \frac{e_{13}^{\,2}}{2}.$ (2.129)

Thus, it follows from (2.125) that

$\displaystyle \epsilon \simeq \frac{5}{4}\,m.$ (2.130)

For the case of the Earth ( $ \widehat{\omega}=7.27\times 10^{-5}\,{\rm rad.\,s^{-1}}$ , $ a_0=6.37\times 10^6\,{\rm m}$ , $ g_0=9.81\,{\rm m\,s^{-1}}$ --Yoder 1995), we obtain

$\displaystyle m \simeq \frac{1}{291}.$ (2.131)

Thus, it follows that, were the Earth homogeneous, its figure would be a spheroid, flattened at the poles, of ellipticity

$\displaystyle \epsilon \simeq \frac{5}{4}\,\frac{1}{291}\simeq \frac{1}{233}.$ (2.132)

This result was first obtained by Newton. The actual ellipticity of the Earth is about $ 1/294$ (Yoder 1995), which is substantially smaller than Newton's prediction. The discrepancy is due to the fact that the Earth is strongly inhomogeneous, being much denser at its core than in its outer regions.

Figure: Normalized angular momentum of a Maclaurin spheroid (solid) and a Jacobi ellipsoid (dashed) versus the eccentricity $ e_{13}$ in the $ x_1$ -$ x_3$ plane.
\begin{figure}
\epsfysize =3in
\centerline{\epsffile{Chapter02/fig5.02.eps}}
\end{figure}

Figures 2.5 and 2.6 illustrate the variation of the normalized angular velocity, $ \widehat{\omega}$ , and angular momentum, $ \widehat{L}$ , of a Maclaurin spheroid with its eccentricity, $ e_{13}$ , as predicted by Equations (2.122)-(2.124). It can be seen, from Figure 2.5, that there is a limit to how large the normalized angular velocity of such a spheroid can become. The limiting value corresponds to $ \widehat{\omega} = 0.47399$ , and occurs when $ e_{13}=0.92995$ . For values of $ \widehat{\omega}$ lying below $ 0.47399$ there are two possible Maclaurin spheroids, one with an eccentricity less than $ 0.92995$ , and one with an eccentricity greater than $ 0.92995$ . Note, however, from Figure 2.6, that despite the fact that the angular velocity, $ \widehat{\omega}$ , of a Maclaurin spheroid varies in a non-monotonic manner with the eccentricity, $ e_{13}$ , the angular momentum, $ \widehat{L}$ , increases monotonically with $ e_{13}$ , becoming infinite in the limit $ e_{13}\rightarrow 1$ . It follows that there is no upper limit to the angular momentum of a Maclaurin spheroid.


next up previous
Next: Jacobi Ellipsoids Up: Hydrostatics Previous: Equilibrium of a Rotating
Richard Fitzpatrick 2016-01-22