User:Tohline/VE/RiemannEllipsoids
From VisTrailsWiki
Contents

SteadyState 2^{nd}Order Tensor Virial Equations
 Tiled Menu  Tables of Content  Banner Video  Tohline Home Page  
Summary
Drawing from our accompanying discussion of virial equations as viewed from a rotating frame of reference, here we employ the 2^{nd}order tensor virial equation (TVE),



to determine the equilibrium conditions of uniformdensity ellipsoids that have semiaxes, and an internal velocity field, (as prescribed below), that preserves this specified ellipsoidal shape, as viewed from a frame of reference that is rotating with angular velocity, . Because each of the indices, and , run from 1 to 3, inclusive, this TVE appears to provide nine equilibrium constraints; and once the values of the density and the three semiaxes are specified, there appear to be seven unknowns: and the three pairs of velocityfield components , , In practice, however, only five constraints are relevant/independent because, as is encapsulated in …
Riemann's Fundamental Theorem
… nontrivial solutions are obtained only if no more than two of the three pairs of velocityfield components are different from zero. 
Following EFE, we will set , in which case the only applicable TVE constraint relations are the five identified in the following table of equations.
Indices  Each Associated 2^{nd}Order TVE Expression  









General Coefficient Expressions
In the context of our discussion of configurations that are triaxial ellipsoids, we begin by adopting the subscript notation to identify which semiaxis length is the (largest, mediumlength, smallest). As has been detailed in an accompanying chapter, the gravitational potential anywhere inside or on the surface of an homogeneous ellipsoid may be given analytically in terms of the following three coefficient expressions:









where, and are incomplete elliptic integrals of the first and second kind, respectively, with arguments,

and 

Specific Case of a_{1} > a_{2} > a_{3}
When we discuss configurations in which — such as Jacobi, Dedekind, or most Riemann SType ellipsoids — we must adopt the associations, , , and . This means that the coefficients, , , and are defined by the expressions,









where, the arguments of the incomplete elliptic integrals are,

and 

[ EFE, Chapter 3, §17, Eq. (32) ] 
Specific Case of a_{1} > a_{3} > a_{2}
When we discuss configurations in which — these are usually referred to in EFE as prolate SType Riemann ellipsoids — we must instead adopt the associations, , , and . This means that the coefficients, , , and are defined by the expressions,









where, the arguments of the incomplete elliptic integrals of the first and second kind are,

and 

[ EFE, Chapter 7, §48d, footnote to Table VII (p. 143) ] 
NOTE: All irrotational ellipsoids belong to this category of configurations.
Specific Case of a_{2} > a_{1} > a_{3}
When we discuss configurations in which — for example, most Riemann ellipsoids of Types I, II, & III — we must instead adopt the associations, , , and . This means that the coefficients, , , and are defined by the expressions,









where, the arguments of the incomplete elliptic integrals are,

and 

Oblate Spheroids [a_{2} = a_{1} > a_{3}]
Starting with the case of and setting , we recognize, first, that . Hence, we have,



Adopted (Internal) Velocity Field
EFE (p. 130) states that the … kinematical requirement, that the motion , associated with , preserves the ellipsoidal boundary, leads to the following expressions for its components:









[ EFE, Chapter 7, §47, Eq. (1) ] 
Equilibrium Expressions
[EFE §11(b), p. 22] Under conditions of a stationary state, [the tensor virial equation] gives,



[This] provides six integral relations which must obtain whenever the conditions are stationary.
When viewing the (generally ellipsoidal) configuration from a rotating frame of reference, the 2^{nd}order TVE takes on the more general form:



[ EFE, Chapter 2, §12, Eq. (64) ] 
EFE (p. 57) also shows that … The potential energy tensor … for a homogeneous ellipsoid is given by



[ EFE, Chapter 3, §22, Eq. (128) ] 
where



[ EFE, Chapter 3, §22, Eq. (129) ] 
is the moment of inertia tensor. Expressions for all nine components of the kinetic energy tensor, are derived in Appendix E, below; and expressions for each of the six Coriolis components can be found in Appendices B, C, & D.
The Three Diagonal Elements
For , we have,


















Once we choose the values of the (semi) axis lengths of an ellipsoid — from which the value of can be immediately determined — along with a specification of , this equation has the following five unknowns: . Similarly, for ,















This gives us a second equation, but an additional pair of (for a total of seven) unknowns: . For the third diagonal element — that is, for — we have,















This gives us three equations vs. seven unknowns.
OffDiagonal Elements
Notice that the offdiagonal components of both and are zero. Hence, the equilibrium expression that is dictated by each offdiagonal component of the 2^{nd}order TVE is,



For example — as is explicitly illustrated on p. 130 of EFE — for and ,




[ EFE, Chapter 7, §47, Eq. (3) ] 
whereas for and ,




[ EFE, Chapter 7, §47, Eq. (4) ] 
Given our adoption of a uniformdensity configuration whose surface has a precisely ellipsoidal shape and, along with it, our adoption of the above specific prescription for the internal velocity field, , we recognize that,
This has allowed us to set to zero one of the integrals in each of these last two expressions. In what follows, we will benefit from recognizing, as well, that,

Our first offdiagonal element is, then,












The second is,












How Solution is Obtained
Adding this pair of governing expressions we obtain,






[ EFE, Chapter 7, §47, Eq. (6) ] 
and subtracting the pair gives,






[ EFE, Chapter 7, §47, Eq. (7) ] 
Various Degrees of Simplification
Riemann Ellipsoids of Types I, II, & III
In this, most general, case, the two vectors and are not parallel to any of the principal axes of the ellipsoid, and they are not aligned with each other, but they both lie in the plane — that is to say, . For a given specified density and choice of the three semiaxes , all five of the expressions displayed in our above Summary Table must be used in order to determine the equilibrium configuration's associated values of the five unknowns: . Here we show how these five unknowns can be derived from the five constraint equations, closely following the analysis that is presented in §47 (pp. 129  132) of [ EFE ].
Constraints Due to OffDiagonal Elements
We begin by subtracting the constraint equation provided by the first offdiagonal element from the constraint equation provided by the second offdiagonal element . This gives,












[ EFE, Chapter 7, §47, Eq. (11) ] 
Adding the two instead gives,









[ EFE, Chapter 7, §47, Eq. (10) ] 
The first of these relations cleanly gives an expression for the frequency ratio, , in terms of the other frequency ratio, . This allows us to rewrite the second relation in terms of the ratio, , alone. We obtain,









ASIDE: Alternatively, given that,
the quadratic equation that governs the value of the frequency ratio, is …
Now, in our discussion of Riemann SType ellipsoids, there is also a quadratic equation that governs the equilibrium frequency ratio, . It is, specifically,
Notice that the first and third terms of this quadratic equation exactly match the first and third terms of the quadratic equation, which we have just derived, that governs the same frequency ratio in Riemann ellipsoids of Types I, II & III. Does the second term match? That is, is the coefficient of the linear term the same in both quadratic relations? Well, …
Even appreciating that we can make the substitution, , I don't see any way that this coefficient expression can be manipulated to match the associated coefficient in the other expression, namely, . 
This is a quadratic equation whose solution gives,



For the other frequency ratio we therefore find,









SUMMARY: Riemann Ellipsoids of Types I, II, & III
As is emphasized in EFE (Chapter 7, §47, p. 131) "… the signs in front of the radicals, in the two expressions, go together. Furthermore, "the two roots … correspond to the fact that, consistent with Dedekind's theorem, two states of internal motions are compatible with the same external figure." As has also been pointed out in EFE (Chapter 7, §51, p. 158), from the steps that have led to the development and solution of the above pair of quadratic equations we can demonstrate that the following relations also hold:

Constraints Due to Diagonal Elements
Next, to simplify manipulations, let's replace the frequency ratios by these newly defined — and known — parameters, and , in the three diagonalelement expressions that are written out in our above Summary Table.
Indices  Rewritten DiagonalElement Expressions  





Using the element to preplace in the other two expressions, we obtain,



and,



Inserting the various relations highlighted above, these two expressions may be rewritten as,






and,









Together, then,















[ EFE, Chapter 7, §51, Eq. (170) ] 
where,



[ EFE, Chapter 3, §21, Eqs. (105) & (107) ] 
Similarly, given that (see just above),






we have,









[ EFE, Chapter 7, §51, Eq. (171) ] 
Finally, looking back at the constraint and recognizing that,



we find,












Riemann SType Ellipsoids
In this case, we assume that and are aligned with each other and, as well, are aligned with the axis; that is to say, in addition to setting we also set . So, there are only three unknowns — — and they can be determined by ignoring offaxis expressions and simultaneously solving the diagonal element expressions displayed in our above Summary Table. Furthermore, two of the three diagonalelement expressions can be simplified because we are setting . The three relevant equilibrium constraints are:
Indices  2^{nd}Order TVE Expressions that are Relevant to Riemann SType Ellipsoids  





The component expression immediately identifies the value of one of the unknowns, namely,



From the remaining pair of diagonalelement expressions, we therefore have,



and,



Multiplying the first of these two expressions through by and the second through by , then subtracting the second from the first gives,















[ EFE, Chapter 7, §48, Eq. (30) ] 
Note that — as EFE has done and as we have recorded in a related discussion — the first term on the righthandside of this last expression can be expressed more compactly in terms of the coefficient, .
Alternatively, dividing the first expression through by and the second by , then adding the pair of expressions gives,









If we divide through by 2, then replace the product, , in this expression by the relation derived immediately above, we have,












[ EFE, Chapter 7, §48, Eq. (29) ] 
It has become customary to characterize each Riemann SType ellipsoid by the value of its equilibrium frequency ratio,



in which case the relevant pair of constraint equations becomes,
and,

These two equations can straightforwardly be combined to generate a quadratic equation for the frequency ratio, . Then, once the value of has been determined, either expression can be used to determine the corresponding equilibrium value for in the unit of . The fact that the value of is determined from the solution of a quadratic equation underscores the realization that, for a given specification of the ellipsoidal geometry , if an equilibrium exists — i.e., if the solution for is real rather than imaginary — then two equally valid, and usually different (i.e., nondegenerate), values of will be realized. This means that two different underlying flows — one direct and the other adjoint — will sustain the shape of the ellipsoidal configuration, as viewed from a frame that is rotating about the axis with frequency, .
Jacobi and Dedekind Ellipsoids
Describe …
Maclaurin Spheroids
Describe …
Appendices: Various Integrals Over Ellipsoid Volume
Throughout this set of appendices, we work with a uniformdensity ellipsoid whose surface is defined by the expression,



Appendix A: Volume
Here we seek to find the volume of the ellipsoid via the Cartesian integral expression,



Preliminaries
First, we will integrate over and specify the integration limits via the expression,



second, we will integrate over and specify the integration limits via the expression,



third, we will integrate over and set the limits of integration as .
Carry Out the Integration
Following thestrategy that has just been outlined, we have,


















Appendix B: Coriolis Component u_{1}x_{2}






























[ EFE, Chapter 7, §47, p. 130, Eq. (9a) ] 
Appendix C: Coriolis Component u_{1}x_{3}
Here we will additionally make use of the integration limits,



Integration over the relevant Coriolis component gives,



























[ EFE, Chapter 7, §47, p. 130, Eq. (9b) ] 
Appendix D: The Other Four Coriolis Components
It follows that,
























Appendix E: Kinetic Energy Components
Diagonal Elements















Similarly,


















OffDiagonal Elements


















[ EFE, Chapter 7, §47, p. 130, Eq. (8) ] 
Similarly,
























And, finally,






and, 



See Also
© 2014  2020 by Joel E. Tohline 