next up previous
Next: Roche Ellipsoids Up: Hydrostatics Previous: Maclaurin Spheroids

Jacobi Ellipsoids

If $ a_2\neq a_1$ (i.e., if the outer boundary of the rotating body is ellipsoidal, rather than spheroidal) then the constraint (2.115) can only be satisfied when

$\displaystyle \int_0^\infty\left[\frac{a_1^{\,2}\,a_2^{\,2}}{(a_1^{\,2}+u)\,(a_2^{\,2}+u)}- \frac{a_3^{\,2}}{(a_3^{\,2}+u)}\right]\frac{du}{\mit\Delta} = 0.$ (2.133)

Without loss of generality, we can assume that $ a_1\geq a_2\geq a_3$ . Let

$\displaystyle a_2$ $\displaystyle =a_1\,\cos\beta,$ (2.134)
$\displaystyle a_3$ $\displaystyle =a_1\,\cos\gamma,$ (2.135)

where $ \gamma\geq \beta$ . It follows that the cross-sections of the planet's outer boundary in the $ x_1$ -$ x_2$ and $ x_1$ -$ x_3$ planes are ellipses of eccentricities

$\displaystyle e_{12}$ $\displaystyle =(1-a_3^{\,2}/a_2^{\,2})^{1/2}= \sin\beta,$ (2.136)
$\displaystyle e_{13}$ $\displaystyle = (1-a_3^{\,2}/a_1^{\,2})^{1/2} = \sin\gamma,$ (2.137)

respectively. It is also helpful to define

$\displaystyle \alpha = \sin^{-1}(\sin\beta/\sin\gamma).$ (2.138)

Let $ \sin^2\theta = [a_1^{\,2}/(a_1^{\,2}+u)]\,\sin^2\gamma$ . Here, $ u=0$ corresponds to $ \theta=\gamma$ , and $ u=\infty$ to $ \theta=0$ . Equations (2.115) and (2.114) transform to (Darwin 1886)

$\displaystyle E(\gamma,\alpha)- 2\,F(\gamma,\alpha) +\left[\frac{1+(\sin\alpha\,\tan\beta\,\cos\gamma)^2}{\cos^2\alpha}\right] E(\gamma,\alpha)$      
$\displaystyle -\frac{\sin\alpha\,\tan\beta\,\cos\gamma\,(1+\sin^2\beta)}{\cos^2\alpha}$ $\displaystyle = 0,$ (2.139)

and

$\displaystyle \widehat{\omega}^{\,2} =2\left[\frac{F(\gamma,\alpha)-E(\gamma,\a...
...n^3\gamma\,\cos^2\alpha}-\frac{\cos^2\beta}{\tan^2\gamma\,\cos^2\alpha}\right],$ (2.140)

respectively, where

$\displaystyle E(\gamma,\alpha)$ $\displaystyle = \int_0^\gamma (1-\sin^2\alpha \,\sin^2\theta)^{1/2}\,d\theta,$ (2.141)
$\displaystyle F(\gamma,\alpha)$ $\displaystyle = \int_0^\gamma \frac{d\theta}{(1-\sin^2\alpha \,\sin^2\theta)^{1/2}},$ (2.142)

are special functions known as incomplete elliptic integrals (Abramowitz and Stegun 1965). The integral $ \alpha _0$ , defined in Equation (2.107), transforms to

$\displaystyle \alpha_0 = \frac{2\,(\cos\beta\,\cos\gamma)^{1/3}}{a_0\,\sin\gamma}\,F(\gamma,\alpha).$ (2.143)

Finally, making use of some of the analysis in the previous two sections, the normalized angular momentum, and normalized mechanical energy, of the planet can be written

$\displaystyle \widehat{L}$ $\displaystyle = \frac{\sqrt{6}}{10}\,\frac{1+\cos^2\beta}{(\cos\beta\,\cos\gamma)^{2/3}}\,\widehat{\omega},$ (2.144)
$\displaystyle \widehat{E}$ $\displaystyle =-\frac{3}{5}\,\frac{(\cos\beta\,\cos\gamma)^{1/3}}{\sin\gamma}\,...
...}\,\frac{1+\cos^2\beta}{(\cos\beta\,\cos\gamma)^{2/3}}\,\widehat{\omega}^{\,2},$ (2.145)

respectively.

The constraint (2.139) is obviously satisfied in the limit $ \beta\rightarrow 0$ , because this implies that $ \alpha\rightarrow0$ and $ E(\gamma,\alpha)$ , $ F(\gamma,\alpha)\rightarrow \gamma$ . Of course, this limit corresponds to the axisymmetric Maclaurin spheroids discussed in the previous section. Carl Jacobi (1804-1851), in 1834, was the first researcher to obtain the very surprising result that Equation (2.139) also has non-axisymmetric ellipsoidal solutions characterized by $ \beta>0$ . These solutions are known as the Jacobi ellipsoids in his honor. The properties of the Jacobi ellipsoids, as determined from Equations (2.139), (2.140), (2.144), and (2.145), are set out in Table 2.2, and illustrated in Figures 2.5 and 2.6. It can be seen that the sequence of Jacobi ellipsoids bifurcates from the sequence of Maclaurin spheroids when $ e_{13}=0.81267$ . Moreover, there are no Jacobi ellipsoids with $ e_{13}< 0.81267$ . However, as $ e_{13}$ increases above this critical value, the eccentricity, $ e_{12}$ , of the Jacobi ellipsoids in the $ x_1$ -$ x_2$ plane grows rapidly, approaching unity as $ e_{13}$ approaches unity. Thus, in the limit $ e_{13}\rightarrow 1$ , in which a Maclaurin spheroid collapses to a disk in the $ x_1$ -$ x_2$ plane, a Jacobi ellipsoid collapses to a line running along the $ x_1$ -axis. Note, from Figures 2.5 and 2.6, that, at fixed $ e_{13}$ , the Jacobi ellipsoids have lower angular velocity and angular momentum than Maclaurin spheroids (with the same mass and volume). Furthermore, as is the case for a Maclaurin spheroid, there is a maximum angular velocity that a Jacobi ellipsoid can have (i.e., $ \widehat{\omega}=0.43257$ ), but no maximum angular momentum.


Table 2.2: Properties of the Jacobi ellipsoids.
$ e_{12}$ $ e_{13}$ $ \widehat{\omega}$ $ \widehat{L}$ $ -\widehat{E}$ $ e_{12}$ $ e_{13}$ $ \widehat{\omega}$ $ \widehat{L}$ $ -\widehat{E}$
0.00 0.81267 0.43257 0.30375 0.50452 0.60 0.85585 0.42827 0.30984 0.50138
0.05 0.81293 0.43257 0.30375 0.50459 0.65 0.86480 0.42609 0.31296 0.49975
0.10 0.81372 0.43257 0.30375 0.50459 0.70 0.87510 0.42288 0.31760 0.49734
0.15 0.81504 0.43256 0.30377 0.50458 0.75 0.88705 0.41807 0.32462 0.49372
0.20 0.81691 0.43253 0.30380 0.50457 0.80 0.90102 0.41069 0.33562 0.48814
0.25 0.81934 0.43248 0.30388 0.50453 0.85 0.91761 0.39879 0.35390 0.47908
0.30 0.82237 0.43237 0.30402 0.50445 0.90 0.93778 0.37787 0.38783 0.46295
0.35 0.82603 0.43220 0.30427 0.50432 0.95 0.96340 0.33353 0.46860 0.42782
0.40 0.83037 0.43191 0.30468 0.50410 0.96 0.96950 0.31776 0.50078 0.41499
0.45 0.83544 0.43146 0.30532 0.50376 0.97 0.97605 0.29691 0.54672 0.39771
0.50 0.84131 0.43078 0.30628 0.50326 0.98 0.98317 0.26722 0.62003 0.37241
0.55 0.84808 0.42976 0.30772 0.50250 0.99 0.99101 0.21809 0.76872 0.32842


Figure 2.7: Normalized mechanical energy of a Maclaurin spheroid (solid) and a Jacobi ellipsoid (dashed) versus the normalized angular momentum.
\begin{figure}
\epsfysize =3in
\centerline{\epsffile{Chapter02/fig5.03.eps}}
\end{figure}

Figure 2.7 shows the mechanical energy of the Maclaurin spheroids and Jacobi ellipsoids plotted as a function of their angular momentum. It can be seen that the Jacobi ellipsoid with a given angular momentum has a lower energy that the corresponding Maclaurin spheroid (i.e., the spheroid with the same angular momentum, mass, and volume). This is significant because, in the presence of a small amount of dissipation (i.e., viscosity), we would generally expect an isolated fluid system to slowly evolve toward the equilibrium state with the lowest energy, subject to any global constraints on the system. For the case of a weakly viscous, isolated, rotating, liquid planet, the relevant constraints are that the mass, volume, and net angular momentum of the system cannot spontaneously change. Thus, we expect such a planet to evolve toward the equilibrium state with the lowest energy for a given mass, volume, and angular momentum. This suggests, from Figure 2.7, that at relatively high angular momentum (i.e., $ \widehat{L}>0.30375$ , $ e_{13}>0.81267$ ), when the Jacobi ellipsoid solutions exist, they are stable equilibrium states (because there is no lower energy state to which the system can evolve), whereas the Maclaurin spheroids are unstable. On the other hand, at relatively low angular momentum (i.e., $ \widehat{L}<0.30375$ , $ e_{13}< 0.81267$ ), when there are no Jacobi ellipsoid solutions, the Maclaurin spheroids are stable equilibrium states (again, because there is no lower energy state to which they can evolve). These predictions are borne out by the results of direct stability analysis performed on the Maclaurin spheroids and Jacobi ellipsoids (Chandrasekhar 1969). In fact, such stability studies demonstrate that the Maclaurin spheroids are unstable in the presence of weak dissipation for $ e_{13}>0.81267$ , and unconditionally unstable for $ e_{13}>0.95289$ . The Jacobi ellipsoids, on the other hand, are unconditionally stable for $ e_{13}<0.93858$ , but are unconditionally unstable for $ e_{13}>0.93858$ , evolving toward lower energy ``pear shaped'' equilibria (which are, themselves, unstable in the presence of weak dissipation).


next up previous
Next: Roche Ellipsoids Up: Hydrostatics Previous: Maclaurin Spheroids
Richard Fitzpatrick 2016-03-31