Joel E. Tohline
tohline@rouge.phys.lsu.edu
[ Download PDF file dated ]
As soon as rotation is introduced into the problem, two general complications arise:
Furthermore, even in the case of spherically symmetric objects, it is rare that an equilibrium solution can be determined analytically (cf., our discussion herein of spherical polytropes). Hence, as a rule, the equilibrium properties of rotating, selfgravitating configurations are not describable in terms of analytical functions. (As is discussed in detail in an accompanying chapter, Maclaurin Spheroids are a notable exception to this rule and thereby play an extremely important role in our understanding of the structure and stability of rotating, selfgravitating objects.)
As a foundation for the discussion of the structure of axisymmetric objects sporting a variety of degrees of compressibility and an assortment of angular momentum distributions, we present here the set of equations that govern the equilibrium of rotating, selfgravitating fluid systems.
We seek to construct equilibrium, selfgravitating configurations that are axisymmetric and have simple rotation. Such structures are most easily described mathematically in cylindrical coordinates. (See the accompanying appendix for a variety of relevant expressions.) For example, to implement the "axisymmetric" constraint in cylindrical coordinates, we set ¶_{j} ® 0 in all spatial operators. By doing this, we derive the following
Ñ  =  e_{v}¶_{v} + e_{z}¶_{z}  [Equation III.F.1] 
Ñ×F  =  (1/v)¶_{v}(v F_{v}) + ¶_{z}F_{z}  [Equation III.F.2] 
Ñ^{2}G  =  (1/v)¶_{v}(v ¶_{v}G) + ¶_{z}(¶_{z}G)  [Equation III.F.3] 
(v×Ñ)F  =  e_{v} [ v_{v}¶_{v}F_{v} + v_{z}¶_{z}F_{v}  v_{j}F_{j}/v ] +  
e_{j} [ v_{v}¶_{v}F_{j} + v_{z}¶_{z}F_{j} + v_{j}F_{v}/v ] +  
e_{z} [ v_{v}¶_{v}F_{z} + v_{z}¶_{z}F_{z} ]  [Equation III.F.4] 
The vector velocity in cylindrical coordinates quite generally takes the form:
[Appendix Equation ?]
[Equation III.F.6]
and, furthermore, that the angular velocity w is only a function of the cylindrical radial coordinate, v. In summary, then, for all equilibrium axisymmetric objects considered throughout this H_Book, we will adopt the following
where,
It may very well be possible to construct equilibrium, axisymmetric models with steadystate flow fields different from the "simple rotation" prescribed here. On physical grounds, it has been argued that w should be uniform on cylinders  i.e., w should not be a function of z  in steadystate systems. But steadystate flows with v_{v} ¹ 0 and v_{z} ¹ 0  such as meridional circulation  may exist for many different types of configurations. In practice, though, systems with more complicated flows have not been studied extensively. (Almost certainly because they are significantly more difficult to construct than models with simple rotation.) Without apology, we will restrict our discussions to configurations with "simple rotation," as just defined. 
With simple rotation imposed, the convective operator [III.F.4] further reduces to the form: 
All equilibrium axisymmetric figures must be in steadystate as viewed by an observer in the inertial reference frame. Hence, in each of the principal governing equations that serve, as a group, to define the configuration's equilibrium structure, we must set ¶_{t} ® 0.
Specifically, in Euler's equation,
and the imposition of simple rotation means, furthermore,
Hence, from Euler's equation we derive the following condition for
e_{v} [ v w^{2} ]  =  ( 1/r )ÑP + ÑF 
=  e_{v} [ ( 1/r ) ¶_{v}P + ¶_{v}F ] + e_{z} [ ( 1/r ) ¶_{z}P + ¶_{z}F ] . 
[Equation III.F.11]
By setting ¶_{t} ® 0 and imposing simple rotation, the equation of continuity and the first law of thermodynamics both are automatically satisfied for any axisymmetric distribution of the mass density and specific entropy. Hence, these two principal governing equations need not be dealt with directly when constructing equilibrium, axisymmetric configurations.
where q is a scalar physical variable. What does the variable q represent in the case of the equation of continuity? What does it represent in the case of the first law of thermodynamics? Show that by setting ¶_{t} ® 0 and imposing simple rotation, the lefthandside of this equation automatically goes to zero for any axisymmetric functional form of the variable q.

Finally, from the mathematical expression listed above [III.F.3] for the twodimensional, axisymmetric Laplacian in cylindrical coordinates, we deduce directly that the gravitational potential in the expression for hydrostatic equilibrium can be determined from the 
(1/v)¶_{v}(v ¶_{v}F) + ¶_{z}(¶_{z}F) = 4pGr .
[Equation III.F.13]
The equilibrium structure of selfgravitating, axisymmetric objects in simple rotation can be determined by solving simultaneously the following set of coupled partial differential equations in conjunction with a chosen barotropic equation of state:
[Equation III.F.14]
[Equation III.F.15]
In terms of the enthalpy, it becomes possible to rewrite the condition for hydrostatic equilibrium [III.F.11] in the following form: 
[Equation III.F.19]
Furthermore, for configurations in simple rotation [III.F.7], it usually is possible to construct a scalar centrifugal potential y(v) such that, 
where h_{o} is a normalization constant. Specifically, for any reasonable function w(v),
[Equation III.F.21]
satisfies this requirement.
Show that the centrifugal potential y(v) for any object that is rotating uniformly with angular velocity W_{o} is,
y =  v^{2} / 2,
[Equation III.F.22] if h_{o} = W_{o}.

With the centrifugal potential defined in this way [III.F.21], the statement of hydrostatic equilibrium for axisymmetric configurations in simple rotation [III.F.19] may be written, 
[Equation III.F.24]
where C_{o} is a constant. It often is easier to construct rotating equilibrium configurations by demanding that this algebraic equation [III.F.24] be satisfied in conjunction with the axisymmetric Poisson equation than it is to construct rotating equilibria by solving the outlined set of coupled partial differential equations [III.F.13], [III.F.14], and [III.F.15] by a bruteforce technique. 
Throughout much of this H_Book, we will rely upon a numerical "selfconsistentfield" (SCF) technique to construct rotating equilibrium configurations having a variety of different degrees of compressibility and a variety of different geometries. As was mentioned briefly in the context of spherical polytropic configurations, through an SCF technique, one iterates back and forth between solutions to the two equations [III.F.13 & III.F.24] that govern the equilibrium structure of selfgravitating, axisymmetric objects in simple rotation until a satisfactory, "selfconsistent" solution to both equations has been determined. 
In particular, we will rely heavily upon the specific SCF technique that was developed by Hachisu in the mid1980's. The accompanying pages describe in detail how an HSCF (Hachisu selfconsistentfield) numerical algorithm is constructed, and provide a guide to the various wwwaccessible HSCF utilities that have been developed at LSU.
Home Page  Preface  Context  Applications  Appendices  Search Index 