User:Tohline/ThreeDimensionalConfigurations/RiemannStype
From VisTrailsWiki
Contents 
Riemann Stype Ellipsoids
 Tiled Menu  Tables of Content  Banner Video  Tohline Home Page  
General Coefficient Expressions
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 

[ EFE, Chapter 3, §17, Eq. (32) ] 
TEST (part 1) Notation: Use in place of . 

Numerical Recipes  
(deg)  (rad)  (deg)  (rad)  
0.9  0.641  50.13357253  0.874995907  32.53852919  0.567904468  0.909025949  0.843118048  0.521450273  0.595131012  0.883418715 
Equilibrium Conditions for Riemann Stype Ellipsoids
Pulling from Chapter 7 — specifically, §48 — of Chandrasekhar's EFE, we understand that the semiaxis ratios, associated with Riemann Stype ellipsoids are given by the roots of the equation,



[ EFE, §48, Eq. (34) ] 
and the associated value of the square of the equilibrium configuration's angular velocity is,



[ EFE, §48, Eq. (33) ] 
where,
(Notice that if we set , this pair of expressions simplifies to the pair we have provided in a separate discussion of the equilibrium conditions for Jacobi ellipsoids.) Following Chandrasekhar's lead and eliminating between these two expressions, we obtain,



[ EFE, §48, Eq. (35) ] 
For a given , this last expression determines the ratios of the axes of the ellipsoids that are compatible with equilibrium; and the value of , that is to be associated with a particular solution of this last expression, then follows from either one of the first two expressions. For convenience of evaluation and for greater clarity, let's rewrite this last (quadratic) equation in the form,



in which case the pair of solutions is,



and the corresponding values of the angular velocity (in units of [G ρ]^{½}) are provided by the expression,



TEST (part 2)  
Direct  Adjoint  
0.9  0.641  0.387793362  0.207337649  0.247245200  3.797483556   0.268008879  ± 1.131374734  15.09117122  ± 0.150771618 
Now, according to Ou (2006), at any coordinate position inside or on the surface of the ellipsoid, , the three components of the velocity as viewed from a frame of rotation that is spinning at the equilibrium configuration's frequency, , are,



where, is an overall scale factor. But, according to §48 of EFE, we see that,



where,



and, 



and is the scalar magnitude of the vorticity vector, . The transformation from EFE's notation to the one used by Ou is, then,



and, 






which, gratifyingly agrees with Ou's equation (17).
Models Examined by Ou (2006)
His Tabulated Model Parameters
Table 1 (see below) lists a subset of the Riemann Stype ellipsoids that were studied by Ou (2006); properties of various socalled Direct configurations can be found in Ou's Table 1, while properties of various Adjoint configurations can be found in his Table 5. Each row of our Table 1 was constructed as follows:
 The pair of axis ratios associated with one of Ou's (2006) uniformdensity, incompressible ellipsoid models (columns 1 and 2 from Ou's Table 1) has been copied into columns 1 and 2 of our table.
 Properties of Direct Configurations …
 The pair of parameter values that is required in order for this to be an equilibrium configuration — as specified by the above set of analytical expressions from EFE — is copied from, respectively, columns 11 and 13 of Ou's Table 1 into columns 3 and 4 of our table; in our table, the "analytic" subscript has been dropped from the column headings.
 The value of the equilibrium configuration's vorticity, — see column 5 of our table — has been determined from the expression,
 Column 6 of our table lists the value of the frequency ratio, .
 Properties of Adjoint Configurations [in order to distinguish from Direct configuration properties, a superscript † has been attached to each parameter name] …
 As listed in column 7 of our Table, the "spin" angular velocity of the adjoint equilibrium configuration has been determined from the vorticity of the direct configuration via the relation,
 As listed in column 10 of our Table, the ratio of the vorticity to the angular velocity in the adjoint equilibrium configuration has been determined from the same ratio in the direct configuration via the relation,
 As indicated, the value of the vorticity in the adjoint equilibrium configuration (column 9 of our table) has been determined from a product of and .
 As listed in column 8 of our table, the value of the parameter, , has been determined from the vorticity in the adjoint equilibrium configuration via the relation,
Table 1: Example Riemann Stype Ellipsoids 

Properties of 
Properties of 

(1)  (2)  (3)  (4)  (5)  (6)  (7)  (8)  (9)  (10) 
0.90  0.795  1.14704  0.43181  0.86842  0.75709  0.43181  1.14704  +2.30682  5.3422 
0.641  1.13137  0.15077   0.30322   0.26801   0.15077  1.13137  2.27531   15.0913  
0.590  1.10661  0.06406  0.12883  0.11642  0.06406  1.10661  +2.22552  34.7411  
0.564  1.09034  0.02033  0.04089  0.03750  0.02033  1.09034  +2.19279  107.86  
0.538  1.07148   0.02324  +0.04674  +0.04362  +0.02324   1.07148  +2.15487  +92.722  
0.487  1.02639   0.10880  +0.21881  +0.21318  +0.10880  1.02639  +2.06418  +18.972  
0.333  0.79257   0.39224  +0.78884  +0.99529  +0.39224  0.79257  +1.59395  +4.06370  
0.28  0.256  0.80944  0.03668  0.14127  0.17453  0.03668  0.80944  +3.11750  84.992 
0.245083  0.796512^{a}  0.0  0.0  0.0  0.0  …  …  
0.231  0.77651   0.04714  +0.18156  +0.23381  +0.04714  0.77651  +2.99067  +63.442  
0.205  0.72853   0.13511  +0.52037  +0.71427  +0.13511  0.72853  +2.80588  +20.7674  
^{a}According to Table IV (p. 103) of EFE, the square of the angular velocity of this Jacobi ellipsoid is, ; from this value, we find that, . 
Our Parameter Determinations
The parameter values that have been posted above in our Table 1 are typically given with five digits of precision. This is because, as explained, the values were determined from the analytically determined values, and , that were provided by Ou (2006) with only five digit accuracy. Our Table 2 (shown immediately below) provides values of this same set of model parameters to better than eleven digits accuracy. We calculated these parameter values by following the steps detailed in earlier subsections of this chapter and, as a foundation, using doubleprecision versions of Numerical Recipes algorithms to evaluate the special functions, and . As an example, the above pair of brief tables titled, TEST (part 1) and TEST (part 2) detail all of the intermediate steps that were used in order to determine the highprecision parameter values specifically for the model having the axisratio pair . This table of higher precision parameter values was primarily generated in order to convince ourselves that we understood from first principles how to accurately determine the properties of Riemann Stype ellipsoids; the lowerprecision parameter values that we derived from Ou's work provided a handy means of crosschecking these "first principles" determinations.
In generating our Table 2, we wondered what the approriate signs were of the various model parameters — especially when part of our objective is to distinguish between direct and adjunct configurations. We took the following approach: First we decided that the spin frequency of every direct configuration should be positive. (Evidently, Ou made this same choice.)
Table 2: Example Riemann Stype Ellipsoids (doubleprecision evaluation) 

Properties of 
Properties of 

(1)  (2)  (3)  (4)  (5)  (6)  (7)  (8)  (9)  (10) 
0.90  0.795  +1.147036091720  +0.431809451699  0.868416786194  0.757096320116  0.431809460593  1.147036104571  +2.306817054749  5.342210487323 
0.641  +1.131374738327  +0.150771621841  0.303218483925  0.268008886644  0.150771621877  1.131374730590  +2.275320291519  15.091170863305  
0.590  +1.106612583610  +0.064060198174  0.128832176328  0.116420305902  0.064060197762  1.106612576964  +2.225520849228  34.741086358509  
0.564  +1.090339840378  +0.020334563779  0.040895067155  0.037506716440  0.020334563809  1.090339837153  +2.192794561386  107.8358300897  
0.538  +1.071485625744  0.023236834336  +0.046731855720  +0.043614077664  +0.023236835120  1.071485656401  +2.154876708984  +92.735376233270  
0.487  +1.026387311947  0.108799837242  +0.218808561563  +0.213183225210  +0.108799835209  1.026387320039  +2.064178943634  +18.972261524065  
0.333  +0.792566980901  0.392440787995  +0.789242029190  +0.995804846843  +0.392440793882  0.792566979129  +1.593940258026  +4.061606964516  
0.74  0.692  +1.132149  +0.385992  0.385992  1.132149  
0.41  0.385  +0.971082  +0.141594  0.141594  0.971082  
0.333  +0.929630138695  +0.003311666790  0.009435019456  0.010149218281  0.003311666699  0.929630099681  +2.648538827896  799.7601146950  
0.28  0.256  +0.809436834686  +0.036676037913  0.141255140305  0.174510396110  0.036676038521  0.809436833116  +3.117488145828  85.000678306244 
0.245083  0.796512^{a}  0.0  0.0  0.0  0.0  …  …  
0.231  +0.776514825339  0.047142035397  +0.181564182043  +0.233819345828  +0.047142037070  0.776514835457  +2.990691423416  +63.440011724689  
0.205  +0.728526018042  0.135108121071  +0.520359277725  +0.714263156392  +0.135108125079  0.728526039364  +2.805866003036  +20.767558718483  
^{a}According to Table IV (p. 103) of EFE, the square of the angular velocity of this Jacobi ellipsoid is, ; from this value, we find that, . 
Feeding a 3D Animation
Initial Thoughts
Let's examine the elliptical trajectory of a Lagrangian particle that is moving in the equatorial plane of a Riemann SType ellipsoid. As viewed in a frame that is spinning about the Zaxis at angular frequency, , the trajectory is defined by,
r^{2} 


where . (The surface of the relevant ellipsoid is associated with the value, .)
Let's choose a pair of axis ratios — for example, and — then, from Table 1 of our above discussion, draw the associated value of either or that corresponds to the Jacobilike equilibrium configuration — in this example, and . Then, for any point inside of the ellipsoid, the fluid's velocity components (as viewed from the rotating frame of reference) are,

and, 

Alternatively, we have,

and, 

Now, each Lagrangian fluid element's motion is oscillatory in both the and coordinate directions. So let's see how this plays out. Suppose,

and, 

Then,

and, 

Hence our functional representation of the timedependent behavior of both and works perfectly if, for each orbit inside of or on the surface of the configuration, we set and if the ratio . Hooray!
Preferred Normalizations
Let's do this again, assuming that and both have units of length and that has the unit of time. Then, let's use to normalize lengths and use to normalize time. We therefore have,

and, 

NOTE: When implementing in an xmlbased COLLADA (3D animation) file, we associate with . Hence we can everywhere replace with (in radians) or (in degrees) .

Next, let's normalize the velocities such that and the total mass, , are both assumed to be the same in every examined Riemann ellipsoid. In particular, we will normalize to,





in which case we have,

and, 

Finally, setting, means,

and, 

SType Ellipsoid Example b41c385
The EFE model that we chose to use in our first successful construction of a COLLADAbased, 3D and interactive animation had the following properties (model selected from the above table):
Figure 1 displays a pair of stillframe images of this (purple) ellipsoidal configuration after the ellipsoid has completed precisely five (counterclockwise) spin cycles. (The snapshots have been taken at the same point in time, but from two different "camera" viewing angles.) The cycle of the "wall mounted" clock is based on the fundamental, EFEadopted frequency of [π G ρ]^{½}. In the lefthand image (labeled Figure 1a), the time on the clock appears to be about 9:08. This means that as the ellipsoid has completed five spin cycles, the clock has completed approximately [9 + 8/60] ≈ 9.13 cycles. In other words, the ratio (ellipsoidtoclock) of these two frequencies is,



This matches the tabulated value of presented above. Now, let's examine the motion of an example Lagrangian fluid element, which has been marked in the 3D scene by a red "arrow" riding in the equatorial plane and along the surface of the (purple) ellipsoidal figure. At time zero, the fluid marker was placed at the end of the longest axis of the ellipsoid that was nearest the "wall clock"; then, as time progressed and the ellipsoidal figure turned counterclockwise, the fluid marker moved clockwise and completed less than one full "orbit" in the same time that the ellipsoidal figure completed five full spin cycles. In the righthand image (labeled Figure 1b), we can see that relative to the ellipsoidal figure, the fluid marker has moved through approximately threequarters of its assigned elliptical "orbit"; let's say, 73% of one full cycle. This means that the ratio of the Lagrangian fluid element's orbital frequency to the frequency of the wallclock is,



This matches the tabulated value of presented above.
Type I Ellipsoid Example b1.25c0.470
We have pulled the parameters for this example Type I Ellipsoid from Table XIII(a) on p. 170 of EFE; in addition, the superior values of β_{+} and γ_{+} have been obtained from eqs. (16) and (17), respectively, as found on p. 131 of Chapter 7, §47. This model sits along the locus O_{2}X_{2}^{(I)} shown in Figure 17 (p. 171) of EFE.

EFE Table XIII (p. 170)  

EFE Figure 17 (p. 171)  
In explaining the elements depicted in Figure 2 we will reference two separate cartesian frames: The (inertial) "lab" frame and a "body" frame affixed to the ellipsoid. The xy plane of the inertial frame is identified by the square, grey/black frame and its zaxis is identified by a thin, green, pointed vertical stick. As depicted in both panels of Figure 2, the xaxis of the ellipsoid's body frame aligns with the xaxis of the inertial frame; in panel 2a, we're generally looking in the plusx direction toward the clock, and in both panels a black arrow "tip" is affixed to the ellipsoid's negative xaxis. In this "Type I" ellipsoid — as opposed to the "SType" ellipsoids discussed earlier — the xaxis of the body frame aligns permanently with the intermediate axis of the (purple) ellipsoidal configuration; following EFE, we assign the value, a = 1, to this intermediate length. The body frame's yaxis is permanently aligned with the configurations's longest axis; in the particular model displayed in Figure 2, b/a = 1.25. And the body frame's zaxis is permanently aligned with the configuration's shortest axis; in this case, c/a = 0.470.
In SType ellipsoids (discussed previously), the ellipsoidal configuration spins about its shortest axis — this means that the ellipsoidal configuration's angular velocity vector has no component along either the xaxis or the yaxis (Ω_{2} = Ω_{1} = 0). But for the Type I ellipsoid shown in Figure 2, Ω_{2} = 0.3639 while Ω_{3} = 0.6633. This means that the configuration's spin axis — associated with the inertial frame's zaxis (the green pointed stick in both panels of Figure 2) — is not aligned with the zaxis of the (purple) ellipsoidal configuration. As is illustrated best in panel 2a, the ellipsoidal configuration is tipped about the xaxis (Ω_{1} = 0) by an angle, 61.25° = tan^{1}(Ω_{3}/Ω_{2}) to the spin axis. When this COLLADAbased 3D scene is animated, the purple ellipsoid spins about the zaxis (green pointed stick) of the inertial frame at a frequency, Ω = [Ω_{2}^{2} + Ω_{3}^{2}]^{½} = 0.7566, in units of [π G ρ]^{½}.
In our COLLADAbased animation, we associate TIME=4 with one clock cycle. Given that for this particular example Type I Riemann ellipsoid, we recognize that the number of clock cycles required to complete one spin about the green axis is cycles. Hence spinning through 360° is equivalent to TIME = 5.2868; or, TIME = 4 is equivalent to a spin angle of 272.38°. 
Now … how do we determine the motion of Lagrangian fluid elements due to the nonzero internal vorticity of the ellipsoidal structure? Evidently the vorticity vector is tipped at an angle of 40.74° = tan^{1}(ζ_{3}/ζ_{2}) to the inertial frame's zaxis. [The preceding sentence is incorrect! This is a misinterpretation of how the "tipping angle" is determined. The correct expression for this "tipping angle" is presented below.] This presumably means that, as viewed from the frame of reference of the spinning ellipsoidal figure, each Lagrangian fluid element moves along a closed (approximately elliptical?) path in a plane that is perpendicular to this vorticity vector. How we determine these various trajectories precisely?
 First, we acknowledge that EFE provides an expression for the velocity vector of each fluid element, given its instantaneous bodycoordinate position (x, y, z) = (x_{1}, x_{2}, x_{3}). Namely — see eq. 154 of Chapter 7, §51 (p. 156) —
But we cannot immediately deduce, from these expressions alone, what the timedependent behavior is of the three Lagrangian coordinates (x, y, z).
 Second, we need to determine how to mathematically specify, then display, the coordinates of a plane that is tipped about the xaxis at an arbitrary angle to the equatorial plane of the (purple) ellipsoidal figure.
Regarding the second issue, let's identify the coordinates of a closed curve (an ellipse?) that identifies the the intersection between a "tipped" plane and the surface of the underlying ellipsoid. The surface of the ellipsoid is defined by the equation,



And the "tipped" plane is defined such that, for any/all values of x,



This means,






If we specifically set m = 1 — that is, choose a slope of 45° — and set z_{0} = c — that is, choose a plane that intersects the north pole of the (purple) ellipsoid — we then have,



[In retrospect, we acknowledge that this was a useful way to explore how to determine the shape of each Lagrangianfluidelement orbit. It is important to emphasize, however, that our choice of a 45° tip angle was only meant to provide an example test calculation. Eventually (see below) we determined that, for this particular Type I example, the correct tipping angle is 19.02°.] Remembering that for the chosen ellipsoidal model (a, b, c) = (1.000, 1.250, 0.4703), the following table identifies the coordinates of a collection of points that simultaneously lie on the surface of the ellipsoid and within the chosen 45° "tipped" plane. These points are each represented by a small yellow cube in both panels of Figure 2.
m = 1 and z_0 = c  

z  y  x 
c  0  0 
0.8c  0.2c  0.5953 
0.5c  0.5c  0.8454 
0.25c  0.75c  0.9262 
0  c  0.9265 
0.25c  1.25c  0.8464 
0.5c  1.5c  0.6569 
0.75227c  1.75227c  0 
Here are a few thoughts/possible linesofattack to consider when addressing the first issue stated above.
 It is likely that the angle at which the plane of motion is tilted about the xaxis of the "body" frame is directly related to the ratio, .
 It is likely that the angle at which the plane of motion is tilted with respect to the zaxis of the "body" frame can be deduced by considering the ratio of the velocity components, v_{y} and v_{z}, when the xcomponent of the velocity vector is known to be zero; that is, when the fluid element is "crossing the xaxis."
 Does the closedorbit motion in the tipped plane always form a (offcenter) ellipse?
 Once you know the exact angle at which the orbital plane is tipped about the xaxis, an Eulerangle projection of every fluid particle's velocity vector should reduce to intheplane motion that has only two components. Once this transformation is made — especially if the resulting orbit is always a closed ellipse — we should be able to determine analytic expressions for the timedependent behavior of each fluid element's x and y coordinate positions in the orbit plane and, subsequently, the (x, y, z) coordinate positions in the "body" frame.
Regarding thought "A"
For the example model depicted in Figure 2, we have,



Regarding thought "B"
When a fluid element crosses the xaxis we should expect z = 0. This needs to be rethought because the center of the closed elliptical trajectory of most fluid elements will be offset from the body's inherent zaxis.
Regarding thought "D"









Setting implies,












For our particular example, then,









Now, given that,









we can generate the following expressions for the fluid velocity components in the tipped plane:
























Now, following our earlier analysis for SType ellipsoids, let's assume that,
and,
(Insertion of the extra term, y_{0}, acknowledges that the center of the elliptical orbit in the "tipped" plane will be shifted off of the xaxis in the y_{t}direction; the size of this shift should correlate with z_{t}, that is, with the vertical location of the tipped plane.) We have, then,












So, in terms of the value of z_{t}, the offset in y_{t} must be,



in which case,






and,






In order for both of these expressions to be simultaneously true, we need,












First Trial Value of θ
Now, from above, we determined that the "tip" angle must obey the expression,



Hence, we require,









That is,






So for the example parameters provided above, we have,






Second Trial Value of θ
Instead, what if we obtain the tip angle straight from the tabulated expression,



Then we have,












And this is not physically acceptable. What if, instead, we set,



Then we have,












In our particular example, this means,






Try Again
Methodical Derivation of Orbital Parameters
Let the unprimed coordinates (x, y, z) represent the body frame of the ellipsoid, and the primed coordinates (x', y', z') represent a coordinate system in which the z'axis is "tipped" — about the x' = x axis — away from the zaxis by an angle, θ. We assume that the motion of each Lagrangian fluid element will be restricted to an x'y' equatorial plane — that is, we assume that — but in general its velocity vector in the unprimed "body" frame will have the three nonzero components as specified in EFE and above. Here are some relevant transformations between these two coordinate systems.





NOTE: The center of each elliptical orbit is (x', y', z') = (0, y_{0}, 0). In "body" coordinates, then, (x, y, z) = (0, y_{0} cosθ, z_{0} + y_{0} sinθ).
Focusing on the bottomright quadrant of this equationtable, we note first that , as provided above from EFE. Copying from the upperright quadrant of this equationtable, let's rewrite these three velocity components in terms of the "tipped" plane coordinates, namely,









Now, if we assume that each Lagrangian particle executes a closed elliptical orbit in the plane of the tipped coordinate system (i.e., ), but whose orbitcenter may be shifted by an amount, , away from the z'axis, we expect …



and, 



which implies,



and, 



The first of these equation pairs can be plugged directly into the expressions for u_{1}, u_{2}, and u_{3} — further fleshing out the LHS of the equations in the bottomright quadrant of the above equationtable — while the second pair can be used to reexpress the RHS of the equations in the bottomright quadrant of this equation table. We obtain,









Combining the second and third of these conditions, we find first that,



which gives the "tipping" angle in terms of known parameter values; and second that …






which gives the product of the two unknown quantities, and the ratio , in terms of known parameter values. And from the first condition, we furthermore find that,






The LHS and the RHS must separately sum to zero, which means …



which gives the ratio, in terms of known parameter values; and,



which in combination with the justderived similar product relation can give expressions for both terms, namely,



and,



Example Implementation
Drawing from the set of Case I parameters listed above, we will set a = 1, b = 1.25, c = 0.4703, ζ_{2} = 2.2794, and ζ_{3} = 1.9637. This means that,



We deduce as well that,









It is important to keep in mind that each Lagrangian fluid element will complete one full orbit in its "tipped" plane when , where is in units of and is in units of . Note, as well, that in our COLLADA animations, we have adopted the convention that when . Hence, one complete orbit will be completed when,
That is to say, in the specific example being used here, on complete orbit is concluded when TIME = 19.343. 
Referring back to our earlier geometric prescription of a "tipped plane" we insert this new slope having m = tanθ = 0.3448, and choose a particular plane by setting a value of z_{0}. Then we can map out a locus of points that show the intersection of the chosen plane with the surface of the ellipsoid by specifying various values of z, then calculating the corresponding values of y and x via the respective equations:






For each specified value of z_{0}, the relevant range of z values is given by the pair of values for which x/a = 0. For example, if we set z_0 = b tanθ = 0.4310, this pair is given by the roots of the quadratic expression,









The roots are, then, , and



The middle three columns (table cells having bgcolor="lightblue") in the following table list values of x and y that correspond to various values of z that lie between this pair of limiting values.
z_0 = 0.2500  z_0 = b tanθ = 0.4310  z_0 = 0.6000  
z  y  x  z  y  x  z  y  x  
0.1564  1.1788  0.0000  0  1.2500  0  0.2182  1.1073  0.0000  
0.1147  1.0577  0.4739  0.02   1.1920  0.2981  0.2336  1.0626  0.1749  
0.0729  0.9366  0.6439  0.05  1.1050  0.4553  0.2490  1.0179  0.2377  
0.0312  0.8155  0.7550  0.08  1.0180  0.5548  0.2645  0.9732  0.2787  
0.0106  0.6943  0.8312  0.1  0.9600  0.6042  0.2799  0.9285  0.3068  
0.0524  0.5732  0.8817  0.125  0.8875  0.6521  0.2953  0.8838  0.3255  
0.0941  0.4521  0.9106  0.15  0.8150  0.6879  0.3107  0.8390  0.3361  
0.1359  0.3310  0.9200  0.175  0.7424  0.7133  0.3261  0.7943  0.3396  
0.1776  0.2099  0.9106  0.2   0.6699  0.7293  0.3415  0.7496  0.3361  
0.2194  0.0887  0.8817  0.25   0.5249  0.7356  0.3570  0.7049  0.3255  
0.2612  0.0324  0.8312  0.30   0.3799  0.7076  0.3724  0.6602  0.3068  
0.3029  0.1535  0.7550  0.35   0.2349  0.6410  0.3878  0.6155  0.2787  
0.3447  0.2746  0.6439  0.4   0.0899  0.5210  0.4032  0.5708  0.2377  
0.3865  0.3957  0.4739  0.43   0.0029  0.4050  0.4186  0.5261  0.1749  
0.4282  0.5169  0.0  0.46  +0.0841  0.1950  0.4340  0.4814  0.0000  
0.4685  +0.1088  0.0081 
More generally, letting , the pair of limiting values of z are given by the roots of the quadratic expression,






That is,






And we see that the constraint set on z_{0} is given by the condition,


















Therefore, for this example Type I ellipsoidal configuration, z_{0} must lie within the range,



In the above table, we have detailed the loci of points along the surface trajectories that correspond to the values, z_{0} =  0.2500 (the leftmost three columns of the table) and z_{0} = 0.6000 (the rightmost three columns of the table). The limiting values of z for these to choices of z_{0} are, respectively, [z_{min} = 0.4283, z_{max} = +0.1564] and [z_{min} = 0.4340, z_{max} = 0.2182].
Using COLLADA, we have constructed an animated and interactive 3D scene that displays in purple the surface of our example Type I ellipsoid; panels a and b of Figure 3 show what this ellipsoid looks like when viewed from two different perspectives. (As a reminder — see the explanation accompanying Figure 2, above — the ellipsoid is tilted about the xcoordinate axis at an angle of 61.25° to the equilibrium spin axis, which is shown in green.) Yellow markers also have been placed in this 3D scene at each of the coordinate locations specified in the above table. From the perspective presented in Figure 3b, we can immediately identify three separate, nearly circular trajectories; the largest one corresponds to our choice of z_{0} = 0.25, the smallest corresponds to our choice of z_{0} = 0.60, and the one of intermediate size correspond to our choice of z_{0} = 0.4310. When viewed from the perspective presented in Figure 3a, we see that these three trajectories define three separate planes; each plane is tipped at an angle of θ = 19.02° to the untilted equatorial, xy plane of the purple ellipsoid.
As we have created each orbit in the xmlbased COLLADA file, we have first manually typed in the (TIME, x, y, z) coordinates of each orbit that has a negative value of z_{0}. Specifically, using the Figure 3c projection as a guide/reference, we typed in coordinates for the following orbits:
 FLUID ELEMENT 1:
nstep  TIME  x  y  z 
1  0.000  0.000  1.107  0.218 
26  4.836  0.340  0.794  0.326 
51  9.672  0.000  0.481  0.434 
76  14.507  +0.340  0.794  0.326 
101  19.343  0.000  1.107  0.218 
 FLUID ELEMENT 2:
See Also
 Discussion of Ou's RiemannLike Ellipsoids
 Riemann Meets COLLADA & Oculus Rift S
 Virtual Reality and 3D Printing
 Success Importing Animated Scene into Oculus Rift S
 Carefully (Re)Build Riemann Type S Ellipsoids Inside Oculus Rift Environment: Example (b/a, c/a) = (0.41, 0.385)
 Other Example Stype Riemann Ellipsoids:
Chandrasekhar's Detailed Analysis
 Bernhard Riemann (1876) Gesammelte Mathematische Werke und Wissenschaftlicher, especially Chapter X (p. 168) titled (something along the following line), "A Contribution to Research on Rotating Ellipsoidal Fluids"
 S. Chandrasekhar (1965), ApJ, 142, 890  961. The Equilibrum and the Stability of the Riemann Ellipsoids. I. — This work is referenced as Paper XXV in EFE and focuses on Stype Riemann ellipsoids.
 S. Chandrasekhar (1966), ApJ, 145, 842  877. The Equilibrum and the Stability of the Riemann Ellipsoids. II. — This work is referenced as Paper XXVIII in EFE and focuses on Riemann ellipsoids of Types I, II and III.
FiniteAmplitude Oscillations
 L. F. Rossner (1967), ApJ, 149, 145. The FiniteAmplitude Oscillations of the Maclaurin Spheroids — This work is referenced as Paper XXXVIII in EFE.
 M. Fujimoto (1968), ApJ, 152, 523. Gravitational Collapse of Rotating Gaseous Ellipsoids
 T. T. Chia & S. Y. Pung (1995), Astrophysics and Space Science, 229, issue 2, 215  233. Effects of Variations of Parallel Angular Velocity and Vorticity on the Oscillations of Compressible Homogeneous Rotating Ellipsoids
 T. T. Chia & S. Y. Pung (1997), Astrophysics and Space Science, 254, 269  294. Dynamical Behaviour of Compressible Homogeneous Uniformly Rotating Ellipsoids with Nonparallel Angular Velocity and Vorticity
In the Context of Galaxy Disks
 C. Hunter (1970), ApJ, 162, 97  103. The Disklike Riemann Ellipsoids.
Other
© 2014  2020 by Joel E. Tohline 