User:Tohline/Apps/MaclaurinSpheroids
From VisTrailsWiki
Contents 
Maclaurin Spheroids (axisymmetric structure)
 Tiled Menu  Tables of Content  Banner Video  Tohline Home Page  
There is no particular reason why one should guess ahead of time that the equilibrium properties of any rotating, selfgravitating configuration should be describable in terms of analytic functions. As luck would have it, however, the gravitational potential at the surface of and inside an homogeneous spheroid is expressible analytically. (The potential is constant on concentric spheroidal surfaces that generally have a different axis ratio from the spheroidal mass distribution.) Furthermore, the gradient of the gravitational potential is separable in cylindrical coordinates, proving to be a simple linear function of both and .
If the spheroid is uniformly rotating, this behavior conspires nicely with the behavior of the centrifugal acceleration — which also will be a linear function of — to permit an analytic (and integrable) prescription of the pressure gradient. Not surprisingly, it resembles the functional form of the pressure gradient that is required to balance the gravitational force in uniformdensity spheres.
As a consequence of this good fortune, the equilibrium structure of a uniformly rotating, uniformdensity , axisymmetric configuration can be shown to be precisely an oblate spheroid whose internal properties are describable in terms of analytic expressions. As we show in an accompanying discussion, the angular velocity that is required to keep a selfgravitating spheroid of a specified eccentricity in equilibrium can be obtained from a theorem — we will refer to it as Maclaurin's Theorem — that was derived from purely geometric arguments over 270 years ago by Colin Maclaurin (1742) in A Treatise of Fluxions. This result has been enumerated in many subsequent publications (e.g., Tassoul 1978; Chandrasekhar 1987). It should be appreciated that both Volume I and Volume II of Maclaurin's (1742) Treatise can now be accessed online via Google Books. Selected excerpts from these two volumes are shown in our accompanying discussion.
Properties of UniformDensity Spheroids
Surface Definition
Let be the equatorial radius and the polar radius of a uniformdensity object whose surface is defined precisely by an oblate spheroid. The degree of flattening of the object may be parameterized in terms of the axis ratio , or in terms of the object's eccentricity,
(For an oblate spheroid, ; hence, the eccentricity is restricted to the range .) The meridional crosssection of such a spheroid is an ellipse with the same eccentricity. The foci of this ellipse lie in the equatorial plane of the spheroid at a distance from the minor axis.
Mean Radius
For purposes of normalization, it will be useful to define the mean radius of the spheroid as,
which is equivalent to the radius of a sphere in the limit .
Mass
The total mass of such a spheroid is,
Gravitational Potential
In an accompanying discussion entitled, Properties of Homogeneous Ellipsoids, an expression is given for the gravitational potential at an internal point or on the surface of an homogeneous ellipsoid with semiaxes . For an homogeneous, oblate spheroid in which , this analytic expression defining the potential reduces to the form,
where, the coefficients , , and are functions only of the spheroid's eccentricity. Specifically,









Note that these three expressions have the following values in the limit of a sphere or in the limit of an infinitesimally thin disk :
Limiting Values 













Equilibrium Structure
Governing Relations
To obtain the equilibrium structure of Maclaurin spheroids, we will adopt the technique outlined earlier for determining the structure of axisymmetric configurations. Specifically, the algebraic expression,
,
must be solved in conjunction with the Poisson equation written in cylindrical coordinates for axisymmetric configurations, namely,
Expression for Effective Potential
For any value of the eccentricity, , the above expression for the gravitational potential satisfies this twodimensional Poisson equation. Furthermore, an algebraic expression defining the centrifugal potential inside a uniformly rotating configuration can be drawn from our accompanying table that summarizes the properties of various simple rotation profiles. Together, these relations give us the relevant expression for the effective potential, namely,
Hence, the enthalpy throughout the configuration must be given by the expression,
This expression contains two constants, and , that can be determined from relevant boundary conditions.
Apply Boundary Conditions
The enthalpy should go to zero everywhere on the surface of the spheroid. By pinning the surface down at two points and setting at both of these locations, we can determine the two unknown constants in the above expression. We choose to pin down the edge of the configuration in the equatorial plane — i.e., at — and along the symmetry axis at the pole — i.e., at . From the boundary condition at the pole, we derive the Bernoulli constant, specifically,
and from the boundary condition in the equatorial plane we derive the rotational angular velocity, specifically,






Plugging these constants into the expression for the enthalpy results in the desired solution,



We know from our separate discussion of supplemental, barotropic equations of state that, for a uniformdensity, polytropic configuration, the pressure is related to the enthalpy via the expression, . Hence, we conclude that,






where,
[ T78, §4.5, Eq. (51)^{1} ]
Footnotes
 From our derivation, it is clear that the leading coefficient in the pressure function can be written as . In T78, this coefficient is given as, , which is too large by a factor of .
Related Wikipedia Discussions
© 2014  2019 by Joel E. Tohline 