Editing Apps/MaclaurinSpheroids

__FORCETOC__ 
<!-- __NOTOC__ will force TOC off -->
=Maclaurin Spheroids=
{| class="MaclaurinSpheroids" style="float:left; margin-right: 20px; border-style: solid; border-width: 3px border-color: black"
|- 
! style="height: 125px; width: 125px; background-color:#ffff99;" |[[H_BookTiledMenu#Spheroidal_.26_Spheroidal-Like|<b>Uniform-Density<br />Maclaurin<br />Spheroids</b>]]
|}

There is no particular reason why one should guess ahead of time that the equilibrium properties of ''any'' rotating, self-gravitating 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 <math>\varpi</math> and <math>~z</math>.

If the spheroid is uniformly rotating, this behavior conspires nicely with the behavior of the [[PGE/RotatingFrame#Centrifugal_and_Coriolis_Accelerations|centrifugal acceleration]] &#8212; which also will be a linear function of <math>\varpi</math> &#8212; 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 [[SSC/Structure/UniformDensity|uniform-density spheres]].

As a consequence of this good fortune, the equilibrium structure of a uniformly rotating, uniform-density <math>~(n = 0)</math>, axisymmetric configuration can be shown to be precisely an oblate spheroid whose internal properties are describable in terms of analytic expressions.  As we show [[Apps/MaclaurinSpheroids/GoogleBooks#Interpreting_Maclaurin.27s_Key_Concluding_Theorem|in an accompanying discussion]], the angular velocity that is required to keep a self-gravitating spheroid of a specified eccentricity in equilibrium can be obtained from a theorem &#8212; we will refer to it as [[Apps/MaclaurinSpheroids/GoogleBooks#MaclaurinTheorem|Maclaurin's Theorem]] &#8212; 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 [https://books.google.com/books?id=dCAOAAAAQAAJ&printsec=frontcover&source=gbs_ge_summary_r&cad=0#v=onepage&q&f=false Volume I] and [https://books.google.com/books?id=xfQ7AQAAMAAJ&printsec=frontcover&source=gbs_ge_summary_r&cad=0#v=onepage&q&f=false Volume II] of Maclaurin's (1742) ''Treatise'' can now be accessed online via Google Books.  Selected excerpts from these two volumes are shown [[Apps/MaclaurinSpheroids/GoogleBooks#Excerpts_from_A_Treatise_of_Fluxions|in our accompanying discussion]].

==Properties of Uniform-Density Spheroids==
===Surface Definition===

Let <math>~a_1</math> be the equatorial radius and <math>~a_3</math> the polar radius of a uniform-density 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 <math>~a_3/a_1</math>, or in terms of the object's eccentricity,

<div align="center">
<math>
 e \equiv \biggl[ 1 - \biggl(\frac{a_3}{a_1}\biggr)^2 \biggr]^{1/2} .
</math> 
</div>

(For an ''oblate'' spheroid, <math>~a_3 \leq a_1</math>; hence, the eccentricity is restricted to the range <math>~0 \leq e \leq 1</math>.)  The meridional cross-section 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 <math>~\varpi = ea_1</math> from the minor <math>~(z)</math> axis.

===Mean Radius===
For purposes of normalization, it will be useful to define the mean radius of the spheroid as,

<div align="center">
<math>
~a_\mathrm{mean} \equiv \biggl[a_1^2 a_3 \biggr]^{1/3} = a_1 (1 - e^2)^{1/6} ,
</math> 
</div>
which is equivalent to the radius of a sphere in the limit <math>~a_3 = a_1</math> <math>~(e=0)</math>.

===Mass===
The total mass of such a spheroid is,

<div align="center">
<math>
 M = \frac{4\pi}{3}~a_1^2 a_3 \rho = \frac{4\pi}{3}~a_1^3 \rho (1 - e^2)^{1/2} .
</math> 
</div>

===Gravitational Potential===
In an accompanying discussion entitled, ''[[ThreeDimensionalConfigurations/HomogeneousEllipsoids#Properties_of_Homogeneous_Ellipsoids|Properties of Homogeneous Ellipsoids]]'', an expression is given for the gravitational potential <math>\Phi(\vec{x})</math> at an internal point or on the surface of an homogeneous ellipsoid with semi-axes <math>~(x,y,z) = (a_1,a_2,a_3)</math>.  For an [[ThreeDimensionalConfigurations/HomogeneousEllipsoids#oblate|homogeneous, oblate spheroid]] in which <math>~a_1 = a_2 \geq a_3</math>, this analytic expression defining the potential reduces to the form,

<div align="center">
<p><math>
 \Phi(\varpi,z) = -\pi G \rho \biggl[ I_\mathrm{BT} a_1^2 - \biggl(A_1 \varpi^2 + A_3 z^2 \biggr) \biggr],
</math>
</p>
[<b>[[Appendix/References#ST83|<font color="red">ST83</font>]]</b>], &sect;7.3, p. 169, Eq. (7.3.1)
</div>

where, the coefficients <math>~A_1</math>, <math>~A_3</math>, and <math>~I_\mathrm{BT}</math> are functions only of the spheroid's eccentricity.  Specifically,

<table align="center" border=0 cellpadding="3">
<tr>
  <td align="right">
<math>
~A_1
</math>
  </td>
  <td align="center">
<math>
~=
</math>
  </td>
  <td align="left">
<math>
\frac{1}{e^2} \biggl[\frac{\sin^{-1}e}{e} - (1-e^2)^{1/2}  \biggr](1-e^2)^{1/2} \, ;
</math>
  </td>
</tr>
<tr>
  <td align="right">
<math>
~A_3
</math>
  </td>
  <td align="center">
<math>
~=
</math>
  </td>
  <td align="left">
<math>
\frac{2}{e^2} \biggl[(1-e^2)^{-1/2} -\frac{\sin^{-1}e}{e} \biggr](1-e^2)^{1/2} \, ;
</math>
  </td>
</tr>
<tr>
  <td align="right">
<math>
~I_\mathrm{BT}
</math>
  </td>
  <td align="center">
<math>
~=
</math>
  </td>
  <td align="left">
<math>
~2A_1 + A_3(1-e^2) = 2 (1-e^2)^{1/2} \biggl[ \frac{\sin^{-1}e}{e}\biggr] \, .
</math>
  </td>
</tr>
<tr>
  <td align="center" colspan="3">
[<b>[[Appendix/References#ST83|<font color="red">ST83</font>]]</b>], &sect;7.3, p. 170, Eqs. (7.3.8)
  </td>
</tr>
</table>

<span id="LimitingValues">Note that these three expressions</span> have the following values in the limit of a sphere <math>~(e=0)</math> or in the limit of an infinitesimally thin disk <math>~(e=1)</math>:

<table align="center" border=1 cellpadding="8">
<tr>
  <td colspan="3" align="center">
'''Table 1:''' &nbsp;[[Appendix/Ramblings/PowerSeriesExpressions#Maclaurin_Spheroid_Index_Symbols|Limiting Values]]
  </td>
</tr>
<tr>
  <td align="center">
&nbsp;
  </td>
  <td align="center">
<b><math>e \rightarrow 0</math></b>
  </td>
  <td align="center">
<b><math>\frac{a_3}{a_1} \rightarrow 0</math></b>
  </td>
</tr>
<tr>
  <td align="center">
<b><math>~\frac{\sin^{-1}e}{e}</math></b>
  </td>
  <td align="center">
<math>~1 + \frac{e^2}{6} + \mathcal{O}\biggl(e^4\biggr)</math>
  </td>
  <td align="center">
<math>~\frac{\pi}{2} - \biggl(\frac{a_3}{a_1}\biggr) +\frac{\pi}{4}\biggl(\frac{a_3}{a_1}\biggr)^2
- \mathcal{O}\biggl(\frac{a_3^3}{a_1^3}\biggr)</math>
  </td>
</tr>
<tr>
  <td align="center">
<b><math>A_1</math></b>
  </td>
  <td align="center">
<math>\frac{2}{3}\biggl[1 - \frac{e^2}{5} - \mathcal{O}\biggl(e^4\biggr)\biggr]</math>
  </td>
  <td align="center">
<math>\frac{\pi}{2} \biggl( \frac{a_3}{a_1}\biggr) - 2\biggl(\frac{a_3}{a_1}\biggr)^2+ \mathcal{O}\biggl(\frac{a_3^3}{a_1^3}\biggr)</math>
  </td>
</tr>
<tr>
  <td align="center">
<b><math>A_3</math></b>
  </td>
  <td align="center">
<math>\frac{2}{3}\biggl[1 + \frac{2e^2}{5} + \mathcal{O}\biggl(e^4\biggr)\biggr]</math>
  </td>
  <td align="center">
<math> 2 - \pi \biggl( \frac{a_3}{a_1}\biggr) + 4\biggl(\frac{a_3}{a_1}\biggr)^2
- \mathcal{O}\biggl(\frac{a_3^3}{a_1^3}\biggr)</math>
  </td>
</tr>
<tr>
  <td align="center">
<b><math>I_\mathrm{BT}</math></b>
  </td>
  <td align="center">
<math>2</math>
  </td>
  <td align="center">
<math>0</math>
  </td>
</tr>
</table>

====Example Equi-gravitational-potential Contours====
As an example, let's examine the gravitational potential everywhere inside (and  on the surface) of the oblate spheroid whose properties are presented in the first row of model data in  [[ThreeDimensionalConfigurations/HomogeneousEllipsoids#Table1|Table 1 of our accompanying discussion of the properties of homogeneous ellipsoids]].  That is, let's examine a model with <math>~a_1 = 1.0</math> and &hellip;
<table border="0" align="center" width="80%">
<tr>
  <td align="center"><math>~\frac{a_3}{a_1} = 0.582724 \, ,</math></td>
  <td align="center"><math>~e = 0.81267 \, ,</math></td>
  <td align="center">&nbsp;</td>
</tr>
<tr>
  <td align="center"><math>~A_1 = A_2 = 0.51589042 \, ,</math></td>
  <td align="center"><math>~A_3 = 0.96821916 \, ,</math></td>
  <td align="center"><math>~I_\mathrm{BT} = 1.360556 \, .</math></td>
</tr>
</table>

In the meridional <math>~(\varpi, z)</math> plane, the surface of this oblate-spheroidal configuration &#8212; identified by the thick, solid-black curve below, in Figure 1 &#8212; is defined by the expression,
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~\frac{\varpi^2}{a_1^2} + \frac{z^2}{a_3^2}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left" colspan="2">
<math>~1 </math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~\Rightarrow ~~~ z</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\pm a_3 [1 -  \varpi^2]^{1 / 2}  \, ,</math>
  </td>
  <td align="right">&nbsp; &nbsp; &nbsp; &nbsp; for <math>~0 \le | \varpi | \le 1 \, .</math></td>
</tr>
</table>
Throughout the interior of this configuration, each associated <math>~\Phi_\mathrm{eff}</math> = constant, equipotential surface is defined by the expression,

<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~\phi_\mathrm{choice} \equiv \frac{\Phi_\mathrm{eff}}{\pi G \rho} +  I_\mathrm{BT}a_1^2 </math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left" colspan="1">
<math>~\biggl( A_1 - \frac{\omega_0^2}{2\pi G \rho}\biggr) \varpi^2 + A_3 z^2   </math>
  </td>
</tr>
</table>

(Notice that, written in this manner, <math>~\phi_\mathrm{choice}</math> assumes its minimum value (zero) when <math>~(\varpi, z) = (0, 0)</math>, that is, at the center of the configuration.)  This means that,
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~z </math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\pm \frac{1}{\sqrt{A_3}} \biggl[ \phi_\mathrm{choice} - \biggl( A_1 - \frac{\omega_0^2}{2\pi G \rho}\biggr) \varpi^2\biggr]^{1 / 2} \, .  </math>
  </td>
</tr>
</table>

----

'''No Rotation'''

When we do not consider the effects of rotation and plot, instead, just the equi-gravitational-potential surfaces, then 
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~z </math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\pm \frac{1}{\sqrt{A_3}} \biggl[ \phi_\mathrm{choice} - A_1  \varpi^2\biggr]^{1 / 2} \, .  </math>
  </td>
</tr>
</table>

Because we know that the <math>~\Phi_\mathrm{grav}</math> = constant surfaces are all less flattened than the configuration itself, we should expect that the largest value of the potential that will arise inside &#8212; actually, on the surface of &#8212; the flattened spheroidal configuration will be found at <math>~(\varpi, z) = (1, 0)</math>, that is, when,
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~0 </math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\pm \frac{1}{\sqrt{A_3}} \biggl[ \phi_\mathrm{choice} - A_1  \biggr]^{1 / 2} \, .  </math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~\Rightarrow  ~~~\phi_\mathrm{choice}\biggr|_\mathrm{max}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~ A_1  \, .  </math>
  </td>
</tr>
</table>
So we will plot various equipotential surfaces having, <math>~0 < \phi_\mathrm{choice} < A_1 \, ,</math> recognizing that they will each cut through the equatorial plane <math>~(z = 0)</math> at the radial coordinate given by,
<div align="center">
<math>~\varpi = \sqrt{\phi_\mathrm{choice}/A_1} \, .</math>
</div>
Next, we recognize that the largest equipotential surface that fits entirely within the surface of the oblate spheroidal configuration has the value of the potential that is found on the symmetry axis and at the pole of the spheroid, that is, at <math>~(\varpi, z) = (0, a_3) \, .</math>  For this case we find,
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~\phi_\mathrm{choice}\biggr|_\mathrm{mid}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~a_3^2 A_3 \, .</math>
  </td>
</tr>
</table>
Hence, all equipotential surfaces having <math>~0 < \phi_\mathrm{choice} \le a_3^2 A_3</math> will lie entirely within the spheroid.  But equipotential surfaces having <math>~a_3^2 A_3 < \phi_\mathrm{choice} \le A_1</math> will cut through the surface of the spheroid at the value of <math>~\varpi</math> where "the two values of z<sup>2</sup> match," that is, where,
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~a_3^2(1-\varpi^2)</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\frac{1}{A_3} \biggl[ \phi_\mathrm{choice} - A_1  \varpi^2\biggr] </math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~\Rightarrow~~~\varpi </math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\biggl[ \frac{\phi_\mathrm{choice} - a_3^2 A_3}{ A_1   -  a_3^2 A_3 } \biggr]^{1 / 2} \, .</math>
  </td>
</tr>
</table>

Therefore, for the example model parameters specified  above, our selection of equipotential surfaces to plot should be guided by the following constraints.
<table border="1" align="center" cellpadding="8" width="80%">
<tr>
  <td align="center" colspan="2">Equipotential Contour Lies &hellip;</td>
</tr>
<tr>
  <td align="center" colspan="1" width="50%">Entirely Inside Spheroid's Surface</td>
  <td align="center" colspan="1">Partially Outside Spheroid's Surface</td>
</tr>
<tr>
  <td align="center" colspan="1"><math>~0 < \phi_\mathrm{choice} < 0.32878</math></td>
  <td align="center" colspan="1"><math>~0.32878 < \phi_\mathrm{choice} < 0.51589</math></td>
</tr>
<tr>
  <td align="center" colspan="1"><math>~0 \le \varpi \le (\phi_\mathrm{choice}/0.51589)^{1 / 2}</math></td>
  <td align="center" colspan="1"><math>~\biggl[ \frac{\phi_\mathrm{choice} - 0.32878}{ 0.18711 } \biggr]^{1 / 2} \le \varpi \le (\phi_\mathrm{choice}/0.51589)^{1 / 2}</math></td>
</tr>
<tr><td align="left" colspan="2">
<div align="center">'''Figure 1: &nbsp; Meridional Plane Cross-section'''<br />
[[File:MacAtJacBifurcationJustGravity01.png|550px|center|Maclaurin Spheroid Cross-section at Jacobi Bifurcation]]</div>
''Solid black curve'': Surface of oblate spheroid having a<sub>3</sub>/a<sub>1</sub> = 0.582724.  ''Dashed curves'':  Equi-gravitational-potential contours plotted in increments of <math>~\Delta\phi_\mathrm{choice} = 0.075</math>; specifically, <math>~\phi_\mathrm{choice}</math> = 0.029 (black), 0.104 (dark blue), 0.179 (red), 0.254 (light blue), 0.329 (green), 0.404 (purple), and 0.479 (orange).
</td>
</tr>
</table>


----

'''With Rotation'''
This expression is only applicable to our physical problem under the following conditions:
<ol>
<li>
The argument of the square root must not be negative, that is, <math>~\varpi</math> must be confined to the range,
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~0 \le | \varpi |</math>
  </td>
  <td align="center">
<math>~\le</math>
  </td>
  <td align="left">
<math>~\biggl[ \phi_\mathrm{choice}\biggl( A_1 - \frac{\omega_0^2}{2\pi G \rho}\biggr)^{-1 }\biggr]^{1 / 2} \, .</math>
  </td>
</tr>
</table>
Note that, in turn, in order to ensure that the argument of ''this'' square root is not negative, we should only explore rotation rates for which <math>~\omega_0^2/(2\pi G \rho) \le A_1 \, .</math>
</li>
<li>
In order that our equipotential surface be relevant only to the ''interior'' of our configuration, for every allowed value of <math>~\varpi ,</math> the value of <math>~z</math> corresponding to the potential surface must be less than or equal to the value of  <math>~z</math> at the surface of the configuration.  That is,
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~a_3^2(1-\varpi^2) </math>
  </td>
  <td align="center">
<math>~\le</math>
  </td>
  <td align="left">
<math>~\frac{1}{A_3} \biggl[ \phi_\mathrm{choice} - \biggl( A_1 - \frac{\omega_0^2}{2\pi G \rho}\biggr) \varpi^2\biggr]</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~\Rightarrow ~~~\biggl( A_1 - \frac{\omega_0^2}{2\pi G \rho} - a_3^2 A_3^2 \biggr) \varpi^2 </math>
  </td>
  <td align="center">
<math>~\le</math>
  </td>
  <td align="left">
<math>~\phi_\mathrm{choice} - a_3^2 A_3^2 </math>
  </td>
</tr>
</table>
</li>
</ol>

==Equilibrium Structure==

===Governing Relations===
To obtain the equilibrium structure of Maclaurin spheroids, we will adopt the [[AxisymmetricConfigurations/SolutionStrategies#Technique|technique outlined earlier for determining the structure of axisymmetric configurations]].  Specifically, the algebraic expression,
<div align="center">
<math>~H + \Phi_\mathrm{eff} = C_\mathrm{B}</math> ,
</div>
must be solved in conjunction with the [[AxisymmetricConfigurations/PGE#Governing_Equations|Poisson equation written in cylindrical coordinates for axisymmetric configurations]], namely,
<div align="center">
<math>
~\frac{1}{\varpi} \frac{\partial }{\partial\varpi} \biggl[ \varpi \frac{\partial \Phi}{\partial\varpi} \biggr] + \frac{\partial^2 \Phi}{\partial z^2} = 4\pi G \rho . 
</math>
</div>

===Expression for Effective Potential===
For any value of the eccentricity, <math>~e</math>, the [[Apps/MaclaurinSpheroids#Gravitational_Potential|above expression for the gravitational potential]] satisfies this two-dimensional 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 [[AxisymmetricConfigurations/SolutionStrategies#SRPtable|simple rotation profiles]].  Together, these relations give us the relevant expression for the effective potential, namely,
<div align="center">
<math>
\Phi_\mathrm{eff}(\varpi,z) = \Phi + \Psi =  -\pi G \rho \biggl[ I_\mathrm{BT} a_1^2 - \biggl(A_1 \varpi^2 + A_3 z^2 \biggr) \biggr] - \frac{1}{2}\varpi^2 \omega_0^2 .
</math>
</div> 

Hence, the enthalpy throughout the configuration must be given by the expression,
<div align="center">
<math>
H(\varpi,z) = C_\mathrm{B} + \pi G \rho \biggl[ I_\mathrm{BT} a_1^2 - \biggl(A_1 \varpi^2 + A_3 z^2 \biggr) \biggr] + \frac{1}{2}\varpi^2 \omega_0^2 .
</math>
</div> 
This expression contains two constants, <math>~C_\mathrm{B}</math> and <math>~\omega_0</math>, 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 <math>~H=0</math> 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 &#8212; ''i.e.'', at <math>(\varpi,z) = (a_1,0)</math> &#8212;  and along the symmetry axis at the pole &#8212; ''i.e.'', at <math>(\varpi,z) = (0,a_3)</math>. From the boundary condition at the pole, we derive the Bernoulli constant, specifically,
<div align="center">
<math>
C_\mathrm{B} = - \pi G \rho \biggl[ I_\mathrm{BT} a_1^2 - A_3 a_3^2  \biggr] = 
- \pi G \rho a_1^2 \biggl[ I_\mathrm{BT} - A_3 (1-e^2)  \biggr] ;
</math>
</div> 
and from the boundary condition in the equatorial plane we derive the rotational angular velocity, specifically,
<div align="center" id="EquilibriumFrequency">
<table align="center" border="0" cellpadding="5">
<tr>
  <td align="right">
<math>
\frac{1}{2}a_1^2 \omega_0^2 
</math>
  </td>
  <td align="center">
<math>
~=
</math>
  </td>
  <td align="left">
<math>
- C_\mathrm{B} - \pi G \rho \biggl[ I_\mathrm{BT} a_1^2 - A_1 a_1^2 \biggr] 
</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>
\Rightarrow ~~~~~ \omega_0^2
</math>
  </td>
  <td align="center">
<math>
~=
</math>
  </td>
  <td align="left">
<math>
2\pi G \rho \biggl[ A_1 - A_3 (1-e^2) \biggr] \, .
</math>
  </td>
</tr>
<tr>
  <td align="center" colspan="3">
[<b>[[Appendix/References#T78|<font color="red">T78</font>]]</b>], &sect;4.5, p. 86, Eq. (52)<br />
[<b>[[Appendix/References#ST83|<font color="red">ST83</font>]]</b>], &sect;7.3, p. 172, Eq. (7.3.18)
  </td>
</tr>
</table> 
</div>

Plugging these constants into the expression for the enthalpy results in the desired solution,
<table align="center" border="0" cellpadding="5">
<tr>
  <td align="right">
<math>
H(\varpi,z)
</math>
  </td>
  <td align="center">
<math>
~=
</math>
  </td>
  <td align="left">
<math>
\pi G \rho a_1^2 A_3 (1-e^2)\biggl[1 - \biggl( \frac{\varpi}{a_1} \biggr)^2 - \biggl( \frac{z}{a_3} \biggr)^2  \biggr] .
</math>
  </td>
</tr>
</table>

[[File:CommentButton02.png|right|100px|Comment by J. E. Tohline:  In Tassoul (1978), the leading coefficient in the expression for the pressure &#8212; and, hence, the central pressure &#8212; is too large by a factor of 2.]]We know from our [[SR#Barotropic_Structure|separate discussion of supplemental, barotropic equations of state]] that, for a uniform-density, <math>~n = 0</math> polytropic configuration, the pressure is related to the enthalpy via the expression, <math>~P = H\rho</math>.  Hence, we conclude that,
<table align="center" border="0" cellpadding="5">
<tr>
  <td align="right">
<math>
P(\varpi,z)
</math>
  </td>
  <td align="center">
<math>
~=
</math>
  </td>
  <td align="left">
<math>
\pi G \rho^2 a_1^2 A_3 (1-e^2)\biggl[1 - \biggl( \frac{\varpi}{a_1} \biggr)^2 - \biggl( \frac{z}{a_3} \biggr)^2  
\biggr] 
</math>
  </td>
</tr>
<tr>
  <td align="center" colspan="3">
[<b>[[Appendix/References#T78|<font color="red">T78</font>]]</b>], &sect;4.5, p. 86, Eq. (51)<br />
[<b>[[Appendix/References#ST83|<font color="red">ST83</font>]]</b>], &sect;7.3, p. 172, Eqs. (7.3.16) &amp; (7.3.17)
  </td>
</tr>

<tr>
  <td align="right">
<math>
\Rightarrow ~~~~ \frac{P}{P_0}
</math>
  </td>
  <td align="center">
<math>
~=
</math>
  </td>
  <td align="left">
<math>
\frac{3}{2} A_3 (1-e^2)^{2/3}\biggl[1 - \biggl( \frac{\varpi}{a_1} \biggr)^2 - \biggl( \frac{z}{a_3} \biggr)^2  
\biggr] ,
</math>
  </td>
</tr>
</table>
where,
<div align="center">
<math>~
P_0 \equiv \frac{2}{3}\pi G \rho^2 a_\mathrm{mean}^2 = \frac{2}{3}\pi G \rho^2 a_1^2 (1-e^2)^{1/3} .
</math>
</div>

=See Also=
* [http://en.wikipedia.org/wiki/Lane-Emden_equation Lane-Emden equation]
* [http://en.wikipedia.org/wiki/Polytrope Polytrope]


{{ SGFfooter }}