Editing ThreeDimensionalConfigurations/RiemannTypeI

__FORCETOC__ <!-- will force the creation of a Table of Contents -->
<!-- __NOTOC__ will force TOC off -->

=Riemann Type 1 Ellipsoids=
{| class="RiemannType1" style="float:left; margin-right: 20px; border-style: solid; border-width: 3px border-color: black"
|- 
! style="height: 125px; width: 125px; background-color:#ffeeee;" |[[H_BookTiledMenu#Three-Dimensional_Configurations|<b>Type I<br />Riemann<br />Ellipsoids</b>]]
|}
&nbsp;<br />
&nbsp;<br />
&nbsp;<br />
&nbsp;<br />
&nbsp;<br />
&nbsp;<br />
&nbsp;<br />

==Example b1.25c0.470 ==

<div align="center">
<table border="1" align="center" cellpadding="8">
<tr>
  <th align="center">Figure 2a</th>
  <th align="center">Figure 2b</th>
</tr>
<tr>
  <td align="left" bgcolor="lightgrey">
[[File:TypeI b1.25c470A.png|400px|EFE Model b1.25c470A]]
  </td>
  <td align="left" bgcolor="lightgrey">
[[File:TypeI b1.25c470B.png|400px|EFE Model b1.25c470B]]
  </td>
</tr>
</table>
</div>

<span id="Case_I">
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 &beta;<sub>+</sub> and &gamma;<sub>+</sub> have been obtained from eqs. (16) and (17), respectively, as found on p. 131 of Chapter 7, &sect;47.  This model sits along the locus O<sub>2</sub>X<sub>2</sub><sup>(I)</sup> shown in Figure 17 (p. 171) of EFE.</span>

<table border="0" cellpadding="5" align="center">
<tr><td align="left" rowspan="4">
<table border="1" cellpadding="8" align="center">
<tr><th align="center" colspan="2">Example Type I<br />Ellipsoid</th></tr>
<tr>
  <td align="center"><math>~\frac{b}{a} = \frac{a_2}{a_1}</math></td>
  <td align="center">1.25</td>
</tr>
<tr>
  <td align="center"><math>~\frac{c}{a} = \frac{a_3}{a_1}</math></td>
  <td align="center">0.4703</td>
</tr>
<tr>
  <td align="center"><math>~\Omega_2</math></td>
  <td align="center">0.3639</td>
</tr>
<tr>
  <td align="center"><math>~\Omega_3</math></td>
  <td align="center">0.6633</td>
</tr>
<tr>
  <td align="center"><math>~\tan^{-1} \biggl[ \frac{\Omega_3}{\Omega_2} \biggr]</math></td>
  <td align="center">61.25&deg;</td>
</tr>
<tr>
  <td align="center"><math>~\zeta_2</math></td>
  <td align="center">-2.2794</td>
</tr>
<tr>
  <td align="center"><math>~\zeta_3</math></td>
  <td align="center">-1.9637</td>
</tr>
<tr>
  <td align="center"><math>~\tan^{-1} \biggl[ \frac{\zeta_3}{\zeta_2} \biggr]</math></td>
  <td align="center">40.74&deg;</td>
</tr>
<tr>
  <td align="center"><math>~\beta_+</math></td>
  <td align="center">1.13449 (1.13332)</td>
</tr>
<tr>
  <td align="center"><math>~\gamma_+</math></td>
  <td align="center">1.8052</td>
</tr>
</table>
</td>
<th align="center">EFE Table XIII (p. 170)</th>
</tr>
<tr>
<td align="right">
[[File:EFE TableXIIIp170.cropped.png|500px|EFE Table XIII (p. 170)]]
</td>
</tr>
<tr>
<th align="center">EFE Figure 17 (p. 171)</th>
</tr>
<tr>
<td align="right">[[File:EFE Fig17p171.cropped.png|500px|EFE Figure 17 (p. 171)]]</td>
</tr>
</table>

<span id="explanation">In explaining the elements</span> depicted in Figure 2 we will reference two separate cartesian frames:  The (inertial) "lab" frame and a "body" frame affixed to the ellipsoid.  The x-y plane of the inertial frame is identified by the square, grey/black frame and its z-axis is identified by a thin, green, pointed vertical stick.  As depicted in both panels of Figure 2, the x-axis of the ellipsoid's body frame aligns with the x-axis of the inertial frame; in panel 2a, we're generally looking in the plus-x direction toward the clock, and in both panels a black arrow "tip" is affixed to the ellipsoid's negative x-axis.  In this "Type I" ellipsoid &#8212; as opposed to the "S-Type" ellipsoids discussed earlier &#8212; the x-axis 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 y-axis 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 z-axis is permanently aligned with the configuration's shortest axis; in this case, c/a = 0.470.

In ''S-Type'' ellipsoids (discussed previously), the ellipsoidal configuration spins about its shortest axis &#8212; this means that the ellipsoidal configuration's angular velocity vector has no component along either the x-axis or the y-axis (&Omega;<sub>2</sub> = &Omega;<sub>1</sub> = 0).  But for the ''Type I'' ellipsoid shown in Figure 2, &Omega;<sub>2</sub> = 0.3639 while &Omega;<sub>3</sub> = 0.6633.  This means that the configuration's spin axis &#8212; associated with the inertial frame's z-axis (the green pointed stick in both panels of Figure 2) &#8212; is ''not'' aligned with the z-axis of the (purple) ellipsoidal configuration.  As is illustrated best in panel 2a, the ellipsoidal configuration is tipped ''about the x-axis'' (&Omega;<sub>1</sub> = 0) by an angle, 61.25&deg; = tan<sup>-1</sup>(&Omega;<sub>3</sub>/&Omega;<sub>2</sub>) to the spin axis.  When this COLLADA-based 3D scene is animated, the purple ellipsoid spins about the z-axis (green pointed stick) of the inertial frame at a frequency, |&Omega;| = [&Omega;<sub>2</sub><sup>2</sup> + &Omega;<sub>3</sub><sup>2</sup>]<sup>&frac12;</sup> = 0.7566, in units of [&pi; G &rho;]<sup>&frac12;</sup>.

<table border="1" width="80%" align="center" cellpadding="10">
<tr>
<td align="left">
In our COLLADA-based animation, we associate TIME=4 with one clock cycle.  Given that <math>~|\Omega| = 0.7566 [\pi G \rho]^{1 / 2}</math> 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 <math>~1.0/0.7566 = 1.3217</math> cycles.  Hence spinning through 360&deg; is equivalent to TIME = 5.2868; or, TIME = 4 is equivalent to a spin angle of 272.38&deg;.
</td>
</tr>
</table>

<span id="EFEvelocities">Now &#8230; 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&deg; = tan<sup>-1</sup>(&zeta;<sub>3</sub>/&zeta;<sub>2</sub>)  to the inertial frame's z-axis.  <b><font color="red">[The preceding sentence is incorrect!  This is a misinterpretation of how the "tipping angle" is determined.  The correct expression for this "tipping angle" is [[#Tipped|presented below]].]</font></b>  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?  </span>

<ul>
<li>First, we acknowledge that EFE provides an expression for the velocity vector of each fluid element, given its  instantaneous ''body''-coordinate position (x, y, z) = (x<sub>1</sub>, x<sub>2</sub>, x<sub>3</sub>). Namely &#8212; see eq. 154 of Chapter 7, &sect;51 (p. 156)  &#8212;
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~u_1</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~- \biggl[ \frac{a^2}{a^2 + b^2} \biggr] \zeta_3 y + \biggl[ \frac{a^2}{a^2 + c^2} \biggr] \zeta_2 z</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~+\biggl[ \frac{a^2}{ b^2} \biggr]\Omega_3 \gamma  y - \biggl[ \frac{a^2}{c^2} \biggr] \Omega_2 \beta  z \, ,</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~u_2</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~+\biggl[ \frac{b^2}{a^2 + b^2} \biggr] \zeta_3 x</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~- \Omega_3 \gamma  x \, ,</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~u_3</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~- \biggl[ \frac{c^2}{a^2 + c^2} \biggr] \zeta_2 x</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~-\Omega_2 \beta  x \, .</math>
  </td>
</tr>
</table>
But we cannot immediately deduce, from these expressions alone, what the time-dependent behavior is of the three Lagrangian coordinates (x, y, z).
</li>
<li>
Second, we need to determine how to mathematically specify, then display, the coordinates of a plane that is tipped about the x-axis at an arbitrary angle to the equatorial plane of the (purple) ellipsoidal figure.
</li>
</ul>

<span id="TippedPlane">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,</span>
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~1</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\biggl( \frac{x}{a} \biggr)^2 + \biggl( \frac{y}{b} \biggr)^2 + \biggl( \frac{z}{c} \biggr)^2 \, .</math>
  </td>
</tr>
</table>
And the "tipped" plane is defined such that, for any/all values of x, 
<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>~my + z_0 \, .</math>
  </td>
</tr>
</table>
This means,
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~1</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\biggl( \frac{x}{a} \biggr)^2 + \frac{1}{b^2}\biggl( \frac{z-z_0}{m} \biggr)^2 + \biggl( \frac{z}{c} \biggr)^2 </math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~\Rightarrow ~~~\biggl( \frac{x}{a} \biggr)^2</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~1 - \biggl(\frac{z_0}{mb}\biggr)^2 \biggl( \frac{z}{z_0} - 1 \biggr)^2 - \biggl( \frac{z}{c} \biggr)^2 \, .</math>
  </td>
</tr>
</table>
If we specifically set m = 1 &#8212; that is, choose a slope of 45&deg; &#8212; and set z<sub>0</sub> = c &#8212; that is, choose a plane that intersects the north pole of the (purple) ellipsoid &#8212; we then have,
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~\frac{x}{a} </math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\pm \biggl[ 1 - \biggl(\frac{c}{b}\biggr)^2 \biggl( \frac{z}{c} - 1 \biggr)^2 - \biggl( \frac{z}{c} \biggr)^2 \biggr]^{1 / 2} \, .</math>
  </td>
</tr>
</table>
<b><font color="red">[In retrospect, we acknowledge that this was a useful way to explore how to determine the shape of each Lagrangian-fluid-element orbit. It is important to emphasize, however, that our choice of a 45&deg; tip angle was only meant to provide an example test calculation.  Eventually (see [[#Tipped|below]]) we determined that, for this particular Type I example, the correct tipping angle is -19.02&deg;.]</font></b>  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&deg; "tipped" plane.  These points are each represented by a small yellow cube in both panels of Figure 2.

<table border="1" cellpadding="8" align="center">
<tr>
  <th align="center" colspan="3">m = 1  and  z_0 = c</th>
</tr>
<tr>
  <th align="center">z</th>
  <th align="center">y</th>
  <th align="center">|x|</th>
</tr>
<tr>  <td align="center">c</td>  <td align="center">0</td>  <td align="center">0</td> </tr>
<tr>  <td align="center">0.8c</td>  <td align="center">-0.2c</td>  <td align="center">0.5953</td> </tr>
<tr>  <td align="center">0.5c</td>  <td align="center">-0.5c</td>  <td align="center">0.8454</td> </tr>
<tr>  <td align="center">0.25c</td>  <td align="center">-0.75c</td>  <td align="center">0.9262</td> </tr>
<tr>  <td align="center">0</td>  <td align="center">-c</td>  <td align="center">0.9265</td> </tr>
<tr>  <td align="center">-0.25c</td>  <td align="center">-1.25c</td>  <td align="center">0.8464</td> </tr>
<tr>  <td align="center">-0.5c</td>  <td align="center">-1.5c</td>  <td align="center">0.6569</td> </tr>
<tr>  <td align="center">-0.75227c</td>  <td align="center">-1.75227c</td>  <td align="center">0</td> </tr>
</table>


Here are a few thoughts/possible lines-of-attack to consider when addressing the first issue stated above.
<ol type="A">
<li>
It is likely that the angle at which the plane of motion is tilted about the x-axis of the "body" frame is directly related to the ratio, <math>~\zeta_3/\zeta_2</math>.
</li>
<li>
It is likely that the angle at which the plane of motion is tilted with respect to the z-axis of the "body" frame can be deduced by considering the ratio of the velocity components, v<sub>y</sub> and v<sub>z</sub>, when the x-component of the velocity vector is known to be zero; that is, when the fluid element is "crossing the x-axis."
</li>
<li>
Does the closed-orbit motion '''in''' the tipped plane always form a (off-center) ellipse?
</li>
<li>
Once you know the exact angle at which the orbital plane is tipped about the x-axis, an Euler-angle projection of every fluid particle's velocity vector should reduce to in-the-plane motion that has only two components.  Once this transformation is made &#8212; especially if the resulting orbit is always a closed ellipse &#8212; we should be able to determine analytic expressions for the time-dependent 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.
</li>
</ol>

===Regarding thought "A"===

For the example model depicted in Figure 2, we have,
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~\theta_\mathrm{tip} \equiv \tan^{-1}\biggl( \frac{\zeta_3}{\zeta_2} \biggr)</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\tan^{-1}\biggl( \frac{1.9637}{2.2794} \biggr) = 40.74^{\circ}</math>
  </td>
</tr>
</table>


===Regarding thought "B"===
When a fluid element crosses the x-axis we should expect z = 0.  <b><font color="red">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 z-axis.</font></b>

===Regarding thought "D"===
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~v_x\biggr|_\mathrm{tip}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~v_x\biggr|_\mathrm{body}</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~v_y\biggr|_\mathrm{tip}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~v_y|_\mathrm{body}\cdot \cos\theta + v_z|_\mathrm{body}\cdot \sin\theta</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~v_z\biggr|_\mathrm{tip}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~-v_y|_\mathrm{body}\cdot \sin\theta + v_z|_\mathrm{body}\cdot \cos\theta</math>
  </td>
</tr>
</table>

Setting <math>~v_z\biggr|_\mathrm{tip} = 0</math> implies,

<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~v_y|_\mathrm{body}\cdot \sin\theta </math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~v_z|_\mathrm{body}\cdot \cos\theta</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~\Rightarrow ~~~ \tan\theta </math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\biggl(\frac{v_z}{v_y}\biggr)_\mathrm{body}</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~- \biggl[ \frac{c^2}{a^2 + c^2} \biggr] \zeta_2 \cdot \biggl[ \frac{b^2}{a^2 + b^2} \biggr]^{-1} \zeta_3^{-1}</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~-  \frac{\zeta_2}{\zeta_3}\biggl[ \frac{(c/a)^2}{1 + (c/a)^2} \biggr] \biggl[ \frac{(b/a)^2}{1 + (b/a)^2} \biggr]^{-1} \, .</math>
  </td>
</tr>
</table>
For our particular example, then,


<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~\tan\theta</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~-  \frac{2.2794}{1.9637}\biggl[ \frac{(0.4703)^2}{1 + (0.4703)^2} \biggr] \biggl[ \frac{(1.25)^2}{1 + (1.25)^2} \biggr]^{-1} </math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~-  \frac{2.2794}{1.9637}\biggl[ \frac{(0.4703)^2}{1 + (0.4703)^2} \biggr] \biggl[ \frac{(1.25)^2}{1 + (1.25)^2} \biggr]^{-1} = - 0.3448</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~\Rightarrow ~~~ \theta</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\tan^{-1}(-0.3448) = - 19.0^{\circ} \, . </math>
  </td>
</tr>
</table>

Now, given that,
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~x</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~x_t \, ,</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~y</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~y_t \cos\theta - z_t \sin\theta \, ,</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~z</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~z_t \cos\theta + y_t \sin\theta \, ,</math>
  </td>
</tr>
</table>
we can generate the following expressions for the fluid velocity components ''in'' the tipped plane:
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~v_z\biggr|_\mathrm{tip}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~0\, ,</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~v_x\biggr|_\mathrm{tip}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~v_x\biggr|_\mathrm{body}</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~- \biggl[ \frac{a^2}{a^2 + b^2} \biggr] \zeta_3 y + \biggl[ \frac{a^2}{a^2 + c^2} \biggr] \zeta_2 z</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~- \biggl[ \frac{a^2}{a^2 + b^2} \biggr] \zeta_3 \biggl[ y_t \cos\theta - z_t \sin\theta \biggr] + \biggl[ \frac{a^2}{a^2 + c^2} \biggr] \zeta_2 \biggl[ z_t \cos\theta + y_t \sin\theta \biggr]</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
y_t \biggl\{ \biggl[ \frac{a^2}{a^2 + c^2} \biggr] \zeta_2 \sin\theta 
- \biggl[ \frac{a^2}{a^2 + b^2} \biggr] \zeta_3 \cos\theta \biggr\}
+ z_t\biggl\{ \biggl[ \frac{a^2}{a^2 + b^2} \biggr] \zeta_3 \sin\theta 
+ \biggl[ \frac{a^2}{a^2 + c^2} \biggr] \zeta_2 \cos\theta  \biggr\} \, ,
</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~v_y\biggr|_\mathrm{tip}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~v_y|_\mathrm{body}\cdot \cos\theta + v_z|_\mathrm{body}\cdot \sin\theta</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\biggl[ \frac{b^2}{a^2 + b^2} \biggr] \zeta_3 x \cdot \cos\theta - \biggl[ \frac{c^2}{a^2 + c^2} \biggr] \zeta_2 x \cdot \sin\theta</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~x_t \biggl\{ \biggl[ \frac{b^2}{a^2 + b^2} \biggr] \zeta_3 \cdot \cos\theta - \biggl[ \frac{c^2}{a^2 + c^2} \biggr] \zeta_2 \cdot \sin\theta \biggr\} \, .</math>
  </td>
</tr>
</table>

Now, following our [[#Initial_Thoughts|earlier analysis for ''S-Type'' ellipsoids]], let's assume that,
<div align="center">
<math>~x_t = x_\mathrm{max}\cos(\varphi t)</math> &nbsp; &nbsp; &nbsp; and,  &nbsp; &nbsp; &nbsp; <math>~y_t = y_0 - y_\mathrm{max}\sin(\varphi t) \, .</math>
</div>

(Insertion of the extra term, y<sub>0</sub>, acknowledges that the center of the elliptical orbit in the "tipped" plane will be shifted off of the x-axis in the y<sub>t</sub>-direction; the size of this shift should correlate with z<sub>t</sub>, that is, with the vertical location of the tipped plane.)  We have, then,
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~v_y\biggr|_\mathrm{tip} = \frac{dy_t}{dt} = -y_\mathrm{max} \varphi \cos(\varphi t)</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~x_\mathrm{max} \cos(\varphi t) \biggl\{ \biggl[ \frac{b^2}{a^2 + b^2} \biggr] \zeta_3 \cdot \cos\theta - \biggl[ \frac{c^2}{a^2 + c^2} \biggr] \zeta_2 \cdot \sin\theta \biggr\} \, ,</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~v_x\biggr|_\mathrm{tip} = \frac{dx_t}{dt} = -x_\mathrm{max}\varphi \sin(\varphi t)</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\biggl[ y_0 - y_\mathrm{max} \sin(\varphi t) \biggr]\biggl\{ \biggl[ \frac{a^2}{a^2 + c^2} \biggr] \zeta_2 \sin\theta 
- \biggl[ \frac{a^2}{a^2 + b^2} \biggr] \zeta_3 \cos\theta \biggr\}
+ z_t\biggl\{ \biggl[ \frac{a^2}{a^2 + b^2} \biggr] \zeta_3 \sin\theta 
+ \biggl[ \frac{a^2}{a^2 + c^2} \biggr] \zeta_2 \cos\theta  \biggr\} \, .
</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
- y_\mathrm{max} \sin(\varphi t) \biggl\{ \biggl[ \frac{a^2}{a^2 + c^2} \biggr] \zeta_2 \sin\theta 
- \biggl[ \frac{a^2}{a^2 + b^2} \biggr] \zeta_3 \cos\theta \biggr\}
</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
&nbsp;
  </td>
  <td align="left">
<math>~
+ y_0  \biggl\{ \biggl[ \frac{a^2}{a^2 + c^2} \biggr] \zeta_2 \sin\theta 
- \biggl[ \frac{a^2}{a^2 + b^2} \biggr] \zeta_3 \cos\theta \biggr\}
+ z_t\biggl\{ \biggl[ \frac{a^2}{a^2 + b^2} \biggr] \zeta_3 \sin\theta 
+ \biggl[ \frac{a^2}{a^2 + c^2} \biggr] \zeta_2 \cos\theta  \biggr\} \, .
</math>
  </td>
</tr>
</table>
So, in terms of the value of z<sub>t</sub>, the offset in y<sub>t</sub> must be,
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~ - y_0  \biggl\{ \biggl[ \frac{a^2}{a^2 + c^2} \biggr] \zeta_2 \sin\theta 
- \biggl[ \frac{a^2}{a^2 + b^2} \biggr] \zeta_3 \cos\theta \biggr\}
</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
+ z_t\biggl\{ \biggl[ \frac{a^2}{a^2 + b^2} \biggr] \zeta_3 \sin\theta 
+ \biggl[ \frac{a^2}{a^2 + c^2} \biggr] \zeta_2 \cos\theta  \biggr\} \, ;
</math>
  </td>
</tr>
</table>
in which case,
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~ -y_\mathrm{max} \varphi \cos(\varphi t)</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~x_\mathrm{max} \cos(\varphi t) \biggl\{ \biggl[ \frac{b^2}{a^2 + b^2} \biggr] \zeta_3 \cdot \cos\theta - \biggl[ \frac{c^2}{a^2 + c^2} \biggr] \zeta_2 \cdot \sin\theta \biggr\} \, ,</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~ \Rightarrow ~~~ \varphi </math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~-\biggl[ \frac{x_\mathrm{max}}{y_\mathrm{max}} \biggr]  \biggl\{ \biggl[ \frac{b^2}{a^2 + b^2} \biggr] \zeta_3 \cdot \cos\theta - \biggl[ \frac{c^2}{a^2 + c^2} \biggr] \zeta_2 \cdot \sin\theta \biggr\} \, ,</math>
  </td>
</tr>
</table>
and,
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~-x_\mathrm{max}\varphi \sin(\varphi t)</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
- y_\mathrm{max} \sin(\varphi t) \biggl\{ \biggl[ \frac{a^2}{a^2 + c^2} \biggr] \zeta_2 \sin\theta 
- \biggl[ \frac{a^2}{a^2 + b^2} \biggr] \zeta_3 \cos\theta \biggr\}
</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~\Rightarrow ~~~ \varphi </math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\biggl[ \frac{y_\mathrm{max}}{x_\mathrm{max}} \biggr]\biggl\{ \biggl[ \frac{a^2}{a^2 + c^2} \biggr] \zeta_2 \sin\theta 
- \biggl[ \frac{a^2}{a^2 + b^2} \biggr] \zeta_3 \cos\theta \biggr\} \, .
</math>
  </td>
</tr>
</table>
In order for both of these expressions to be simultaneously true, we need,

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

<tr>
  <td align="right">
<math>~\biggl[ \frac{x_\mathrm{max}}{y_\mathrm{max}} \biggr]  \biggl\{ \biggl[ \frac{c^2}{a^2 + c^2} \biggr] \zeta_2 \cdot \sin\theta - \biggl[ \frac{b^2}{a^2 + b^2} \biggr] \zeta_3 \cdot \cos\theta \biggr\}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\biggl[ \frac{y_\mathrm{max}}{x_\mathrm{max}} \biggr]\biggl\{ \biggl[ \frac{a^2}{a^2 + c^2} \biggr] \zeta_2 \sin\theta 
- \biggl[ \frac{a^2}{a^2 + b^2} \biggr] \zeta_3 \cos\theta \biggr\} 
</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~\Rightarrow ~~~ 
\biggl[ \frac{x_\mathrm{max}}{y_\mathrm{max}} \biggr]  \biggl[ \frac{c^2}{a^2 + c^2} \biggr] \zeta_2 \cdot \sin\theta  
-
\biggl[ \frac{x_\mathrm{max}}{y_\mathrm{max}} \biggr]  \biggl[ \frac{b^2}{a^2 + b^2} \biggr] \zeta_3 \cdot \cos\theta 
</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\biggl[ \frac{y_\mathrm{max}}{x_\mathrm{max}} \biggr]
\biggl[ \frac{a^2}{a^2 + c^2} \biggr] \zeta_2 \sin\theta  
-
\biggl[ \frac{y_\mathrm{max}}{x_\mathrm{max}} \biggr]
\biggl[ \frac{a^2}{a^2 + b^2} \biggr] \zeta_3 \cos\theta 
</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~\Rightarrow ~~~ 
\biggl[ \frac{x_\mathrm{max}}{y_\mathrm{max}} \biggr]  \biggl[ \frac{c^2}{a^2 + c^2} \biggr] \zeta_2 \cdot \sin\theta  
-
\biggl[ \frac{y_\mathrm{max}}{x_\mathrm{max}} \biggr]
\biggl[ \frac{a^2}{a^2 + c^2} \biggr] \zeta_2 \sin\theta  
</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\biggl[ \frac{x_\mathrm{max}}{y_\mathrm{max}} \biggr]  \biggl[ \frac{b^2}{a^2 + b^2} \biggr] \zeta_3 \cdot \cos\theta 
-
\biggl[ \frac{y_\mathrm{max}}{x_\mathrm{max}} \biggr]
\biggl[ \frac{a^2}{a^2 + b^2} \biggr] \zeta_3 \cos\theta 
</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~\Rightarrow ~~~ 
\biggl\{ \biggl[ \frac{x_\mathrm{max}}{y_\mathrm{max}} \biggr] c^2 
-
\biggl[ \frac{y_\mathrm{max}}{x_\mathrm{max}} \biggr]
a^2 \biggl\} \frac{\zeta_2 \sin\theta }{(a^2 + c^2)} 
</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\biggl\{ \biggl[ \frac{x_\mathrm{max}}{y_\mathrm{max}} \biggr]  b^2  
-
\biggl[ \frac{y_\mathrm{max}}{x_\mathrm{max}} \biggr]
a^2 \biggr\} \frac{\zeta_3 \cos\theta }{(a^2 + b^2)} \, .
</math>
  </td>
</tr>
</table>


====First Trial Value of &theta;====
Now, from above, we determined that the "tip" angle must obey the expression,

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

<tr>
  <td align="right">
<math>~\tan\theta </math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~- \biggl[ \frac{c^2}{a^2 + c^2} \biggr] \zeta_2 \cdot \biggl[ \frac{b^2}{a^2 + b^2} \biggr]^{-1} \zeta_3^{-1} \, .</math>
  </td>
</tr>
</table>
Hence, we require,
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~- \biggl[ \frac{c^2}{a^2 + c^2} \biggr] \zeta_2 \cdot \biggl[ \frac{b^2}{a^2 + b^2} \biggr]^{-1} \zeta_3^{-1}
\biggl\{ \biggl[ \frac{x_\mathrm{max}}{y_\mathrm{max}} \biggr] c^2 
-
\biggl[ \frac{y_\mathrm{max}}{x_\mathrm{max}} \biggr]
a^2 \biggl\} \frac{\zeta_2 }{(a^2 + c^2)} 
</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\biggl\{ \biggl[ \frac{x_\mathrm{max}}{y_\mathrm{max}} \biggr]  b^2  
-
\biggl[ \frac{y_\mathrm{max}}{x_\mathrm{max}} \biggr]
a^2 \biggr\} \frac{\zeta_3 }{(a^2 + b^2)} 
</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~ \Rightarrow ~~~  
\biggl\{
\biggl[ \frac{y_\mathrm{max}}{x_\mathrm{max}} \biggr]
a^2 
-  \biggl[ \frac{x_\mathrm{max}}{y_\mathrm{max}} \biggr] c^2 
\biggl\} \frac{c^2 \zeta_2^2 }{(a^2 + c^2)^2} 
</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\biggl\{ \biggl[ \frac{x_\mathrm{max}}{y_\mathrm{max}} \biggr]  b^2  
-
\biggl[ \frac{y_\mathrm{max}}{x_\mathrm{max}} \biggr]
a^2 \biggr\} \frac{b^2\zeta_3^2 }{(a^2 + b^2)^2} 
</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~ \Rightarrow ~~~  
\biggl[ \frac{y_\mathrm{max}}{x_\mathrm{max}} \biggr]^2 
\frac{a^2 c^2 \zeta_2^2 }{(a^2 + c^2)^2} 
- 
\frac{c^4 \zeta_2^2 }{(a^2 + c^2)^2} 
</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\frac{b^4 \zeta_3^2 }{(a^2 + b^2)^2} 
- \biggl[ \frac{y_\mathrm{max}}{x_\mathrm{max}} \biggr]^2 
\frac{a^2 b^2\zeta_3^2 }{(a^2 + b^2)^2} \, .
</math>
  </td>
</tr>
</table>
That is,

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

<tr>
  <td align="right">
<math>~  
\biggl[ \frac{y_\mathrm{max}}{x_\mathrm{max}} \biggr]^2 
</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\biggl[ \frac{b^4 \zeta_3^2 }{(a^2 + b^2)^2} 
+ 
\frac{c^4 \zeta_2^2 }{(a^2 + c^2)^2} \biggr]
\biggl[ \frac{a^2 c^2 \zeta_2^2 }{(a^2 + c^2)^2} 
+ 
\frac{a^2 b^2\zeta_3^2 }{(a^2 + b^2)^2} \biggr]^{-1} 
</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\biggl[ b^4 \zeta_3^2 (a^2 + c^2)^2 
+ 
c^4 \zeta_2^2 (a^2 + b^2)^2 \biggr]
\biggl[ a^2 c^2 \zeta_2^2 (a^2 + b^2)^2 
+ 
a^2 b^2\zeta_3^2 (a^2 + c^2)^2 \biggr]^{-1} \, .
</math>
  </td>
</tr>
</table>

So for the example parameters provided above, we have,

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

<tr>
  <td align="right">
<math>~  
\biggl[ \frac{y_\mathrm{max}}{x_\mathrm{max}} \biggr]^2 
</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\biggl[ 14.04 
+ 
1.669 \biggr]
\biggl[ 7.546  
+ 
8.985 \biggr]^{-1} = 0.9502 
</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~  \Rightarrow~~~
\frac{y_\mathrm{max}}{x_\mathrm{max}}  
</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
0.9748
</math>
  </td>
</tr>
</table>


====Second Trial Value of &theta;====

Instead, what if we obtain the tip angle straight from the tabulated expression,
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~\tan\theta</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\frac{\zeta_3}{\zeta_2} \, .</math>
  </td>
</tr>
</table>

Then we have,

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

<tr>
  <td align="right">
<math>~
\biggl\{ \biggl[ \frac{x_\mathrm{max}}{y_\mathrm{max}} \biggr] c^2 
-
\biggl[ \frac{y_\mathrm{max}}{x_\mathrm{max}} \biggr] a^2 \biggl\} 
\frac{1 }{(a^2 + c^2)} 
</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\biggl\{ \biggl[ \frac{x_\mathrm{max}}{y_\mathrm{max}} \biggr]  b^2  
-
\biggl[ \frac{y_\mathrm{max}}{x_\mathrm{max}} \biggr] a^2 \biggr\} 
\frac{1 }{(a^2 + b^2)} 
</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~\Rightarrow ~~~
\biggl[ \frac{x_\mathrm{max}}{y_\mathrm{max}} \biggr] c^2  (a^2 + b^2)
-
\biggl[ \frac{y_\mathrm{max}}{x_\mathrm{max}} \biggr] a^2  (a^2 + b^2)  
</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\biggl[ \frac{x_\mathrm{max}}{y_\mathrm{max}} \biggr]  b^2 (a^2 + c^2)
-
\biggl[ \frac{y_\mathrm{max}}{x_\mathrm{max}} \biggr] a^2 (a^2 + c^2) 
</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~\Rightarrow~~~
\biggl[ \frac{y_\mathrm{max}}{x_\mathrm{max}} \biggr]^2 
</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\biggl[ c^2  (a^2 + b^2)- b^2 (a^2 + c^2) \biggr]
\biggl[ a^2  (a^2 + b^2)  - a^2 (a^2 + c^2) \biggr]^{-1}
</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
[ c^2 - b^2 ]
[ b^2  - c^2 ]^{-1} = -1 \, .
</math>
  </td>
</tr>
</table>

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

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

<tr>
  <td align="right">
<math>~\tan\theta</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~- \frac{\zeta_3}{\zeta_2} ~~~\Rightarrow ~~~ \theta = -40.74^\circ \, .</math>
  </td>
</tr>
</table>

Then we have,

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

<tr>
  <td align="right">
<math>~
- \biggl\{ \biggl[ \frac{x_\mathrm{max}}{y_\mathrm{max}} \biggr] c^2 
-
\biggl[ \frac{y_\mathrm{max}}{x_\mathrm{max}} \biggr] a^2 \biggl\} 
\frac{1 }{(a^2 + c^2)} 
</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\biggl\{ \biggl[ \frac{x_\mathrm{max}}{y_\mathrm{max}} \biggr]  b^2  
-
\biggl[ \frac{y_\mathrm{max}}{x_\mathrm{max}} \biggr] a^2 \biggr\} 
\frac{1 }{(a^2 + b^2)} 
</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~\Rightarrow ~~~
- \biggl[ \frac{x_\mathrm{max}}{y_\mathrm{max}} \biggr] c^2  (a^2 + b^2)
+
\biggl[ \frac{y_\mathrm{max}}{x_\mathrm{max}} \biggr] a^2  (a^2 + b^2)  
</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\biggl[ \frac{x_\mathrm{max}}{y_\mathrm{max}} \biggr]  b^2 (a^2 + c^2)
-
\biggl[ \frac{y_\mathrm{max}}{x_\mathrm{max}} \biggr] a^2 (a^2 + c^2) 
</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~\Rightarrow ~~~
\biggl[ \frac{y_\mathrm{max}}{x_\mathrm{max}} \biggr] [ a^2  (a^2 + b^2)  + a^2 (a^2 + c^2) ]
</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\biggl[ \frac{x_\mathrm{max}}{y_\mathrm{max}} \biggr] [ b^2 (a^2 + c^2) +  c^2  (a^2 + b^2) ]
</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~\Rightarrow ~~~
\biggl[ \frac{y_\mathrm{max}}{x_\mathrm{max}} \biggr]^2 
</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
[ b^2 (a^2 + c^2) +  c^2  (a^2 + b^2) ] [ a^2  (a^2 + b^2)  + a^2 (a^2 + c^2) ]^{-1} \, .
</math>
  </td>
</tr>
</table>

In our particular example, this means,

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

<tr>
  <td align="right">
<math>~
\biggl[ \frac{y_\mathrm{max}}{x_\mathrm{max}} \biggr]^2 
</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
[ 2.4749 ] [ 3.7837 ]^{-1} = 0.6541 
</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~\Rightarrow ~~~
\biggl[ \frac{y_\mathrm{max}}{x_\mathrm{max}} \bigg] 
</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
0.8088 \, . 
</math>
  </td>
</tr>
</table>

===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" &#8212; about the x' = x axis &#8212; away from the z-axis by an angle, &theta;.  We assume that the motion of each Lagrangian fluid element will be restricted to an x'-y' equatorial plane &#8212; that is, we assume that <math>~\dot{z}' \equiv dz'/dt = 0</math> &#8212; 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.
<table border="1" align="center" cellpadding="10">
<tr>
<td align="center">
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~x'</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~x \, ,</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~y'</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~y \cos\theta + (z-z_0)\sin\theta \, ,</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~z'</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~(z-z_0)\cos\theta - y\sin\theta \, .</math>
  </td>
</tr>
</table>

</td>
<td align="center">

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

<tr>
  <td align="right">
<math>~x</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~x' \, ,</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~y</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~y' \cos\theta - z'\sin\theta \, ,</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~z-z_0</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~z'\cos\theta + y'\sin\theta \, .</math>
  </td>
</tr>
</table>

</td>
</tr>

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

<tr>
  <td align="right">
<math>~\dot{x}'</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\dot{x} \, ,</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~\dot{y}'</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\dot{y} \cos\theta + \dot{z}\sin\theta \, ,</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~\cancelto{0}{\dot{z}'}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\dot{z} \cos\theta - \dot{y}\sin\theta \, .</math>
  </td>
</tr>
</table>

</td>
<td align="center">

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

<tr>
  <td align="right">
<math>~\dot{x}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\dot{x}' \, ,</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~\dot{y}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\dot{y}' \cos\theta - \cancelto{0}{\dot{z}'}\sin\theta \, ,</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~\dot{z}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\cancelto{0}{\dot{z}'}\cos\theta + \dot{y}'\sin\theta \, .</math>
  </td>
</tr>
</table>

</td>
</tr>

</table>


<b><font color="red">NOTE:</font></b>  The center of each elliptical orbit is (x', y', z') = (0, y<sub>0</sub>, 0).  In "body" coordinates, then, (x, y, z) = (0, y<sub>0</sub> cos&theta;, z<sub>0</sub> + y<sub>0</sub> sin&theta;).

Focusing on the bottom-right quadrant of this equation-table, we note first that <math>~(\dot{x}, \dot{y}, \dot{z}) = (u_1, u_2, u_3)</math>, as [[#EFEvelocities|provided above from EFE]].  Copying from the upper-right quadrant of this equation-table, let's rewrite these three velocity components in terms of the "tipped" plane coordinates, namely,
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~\dot{x} = u_1</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
- \biggl[ \frac{a^2}{a^2 + b^2} \biggr] \zeta_3 \biggl[ y'\cos\theta - z'\sin\theta\biggr] 
+ \biggl[ \frac{a^2}{a^2 + c^2} \biggr] \zeta_2 \biggl[ z_0 + z'\cos\theta + y'\sin\theta \biggr] \, ,</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~\dot{y} = u_2</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~+\biggl[ \frac{b^2}{a^2 + b^2} \biggr] \zeta_3 x' \, ,</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~\dot{z} = u_3</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~- \biggl[ \frac{c^2}{a^2 + c^2} \biggr] \zeta_2 x' \, .</math>
  </td>
</tr>
</table>


Now, if we assume that each Lagrangian particle executes a closed elliptical orbit ''in the plane'' of the tipped coordinate system (''i.e.,'' <math>~z' = \dot{z}' = 0</math> ), but whose orbit-center may be shifted by an amount, <math>~y_0</math>, away from the z'-axis, we expect &hellip;
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~x'</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~x_\mathrm{max}\cos(\varphi t)</math>
  </td>
<td align="center">&nbsp; &nbsp; &nbsp; and, &nbsp; &nbsp; &nbsp;
  <td align="right">
<math>~y' - y_0</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~y_\mathrm{max}\sin(\varphi t) \, ,</math>
  </td>
</tr>
</table>
which implies,
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~\dot{x}'</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~- x_\mathrm{max}~ \varphi \cdot \sin(\varphi t)</math>
  </td>
<td align="center">&nbsp; &nbsp; &nbsp; and, &nbsp; &nbsp; &nbsp;
  <td align="right">
<math>~\dot{y}' </math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~y_\mathrm{max}~\varphi \cdot \cos(\varphi t) \, .</math>
  </td>
</tr>
</table>

The first of these equation pairs can be plugged directly into the expressions for u<sub>1</sub>, u<sub>2</sub>, and u<sub>3</sub> &#8212; further fleshing out the LHS of the equations in the bottom-right quadrant of the above equation-table &#8212; while the second pair can be used to re-express the RHS of the equations in the bottom-right quadrant of this equation table.  We obtain,
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~
- \biggl[ \frac{a^2}{a^2 + b^2} \biggr] \zeta_3 \biggl\{  \biggl[ y_0 + y_\mathrm{max} \sin(\varphi t) \biggr] \cos\theta - \cancelto{0}{z'}\sin\theta \biggr\} 
+ \biggl[ \frac{a^2}{a^2 + c^2} \biggr] \zeta_2 \biggl\{ z_0 + \cancelto{0}{z'} \cos\theta + \biggl[ y_0 + y_\mathrm{max} \sin(\varphi t) \biggr]\sin\theta \biggr\} </math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~-x_\mathrm{xmax} ~\varphi \cdot \sin(\varphi t) \, ,</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~+\biggl[ \frac{b^2}{a^2 + b^2} \biggr] \zeta_3 x_\mathrm{max}\cos(\varphi t) </math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~y_\mathrm{max}~\varphi \cdot \cos(\varphi t) \cdot \cos\theta \, ,</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~- \biggl[ \frac{c^2}{a^2 + c^2} \biggr] \zeta_2 x_\mathrm{max}\cos(\varphi t) </math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~y_\mathrm{max}~\varphi \cdot \cos(\varphi t) \cdot \sin\theta \, .</math>
  </td>
</tr>
</table>

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

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

<tr>
  <td align="right">
<math>~\tan\theta </math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~- \frac{\zeta_2 }{ \zeta_3 } \biggl[ \frac{a^2 + b^2}{a^2 + c^2} \biggr]\frac{c^2}{b^2} \, ,</math>
  </td>
</tr>
</table>

which gives the "tipping" angle in terms of known parameter values; and second that &hellip;  
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~\biggl[ \frac{x_\mathrm{max}}{y_\mathrm{max}} \biggr] \biggl[ \frac{b^2}{a^2 + b^2} \biggr] \zeta_3 </math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\varphi \cdot \cos\theta \, ,</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~\Rightarrow ~~~ 
\varphi \biggl[ \frac{y_\mathrm{max}}{x_\mathrm{max}} \biggr]  
</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\biggl[ \frac{b^2}{a^2 + b^2} \biggr] \frac{\zeta_3}{\cos\theta} 
\, ,</math>
  </td>
</tr>
</table>
which gives the product of the two unknown quantities, <math>~\varphi</math> and the ratio <math>~y_\mathrm{max}/x_\mathrm{max}</math>, in terms of known parameter values.  And from the first condition, we furthermore find that,
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~-x_\mathrm{max} ~\varphi \cdot \sin(\varphi t) </math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
- \biggl[ \frac{a^2}{a^2 + b^2} \biggr] \zeta_3  \biggl[ y_0 \cos\theta \biggr]   
- \biggl[ \frac{a^2}{a^2 + b^2} \biggr] \zeta_3 \biggl[ y_\mathrm{max} \sin(\varphi t) \biggr] \cos\theta 
+ 
\biggl[ \frac{a^2}{a^2 + c^2} \biggr] \zeta_2 \biggl[ z_0 + y_0 \sin\theta \biggr]
+ 
\biggl[ \frac{a^2}{a^2 + c^2} \biggr] \zeta_2 \biggl[ y_\mathrm{max} \sin(\varphi t) \biggr]\sin\theta  
</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~
\Rightarrow ~~~ x_\mathrm{max} ~\varphi \cdot \sin(\varphi t) 
+ \biggl\{ \biggl[ \frac{a^2}{a^2 + c^2} \biggr] \zeta_2 \sin\theta  
- \biggl[ \frac{a^2}{a^2 + b^2} \biggr] \zeta_3  \cos\theta \biggr\} \biggl[ y_\mathrm{max} \sin(\varphi t) \biggr]
</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>
\biggl\{
\biggl[ \frac{a^2}{a^2 + b^2} \biggr] \zeta_3  \biggl[ y_0 \cos\theta \biggr]   
-
\biggl[ \frac{a^2}{a^2 + c^2} \biggr] \zeta_2 z_0 
-
\biggl[ \frac{a^2}{a^2 + c^2} \biggr] \zeta_2 \biggl[ y_0 \sin\theta \biggr]
\biggr\} \, .
</math>
  </td>
</tr>
</table>
The LHS and the RHS must separately sum to zero, which means &hellip;
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~
\biggl\{ \biggl[ \frac{a^2}{a^2 + b^2} \biggr] \zeta_3  \cos\theta    
-
\biggl[ \frac{a^2}{a^2 + c^2} \biggr] \zeta_2 \sin\theta  \biggr\}y_0
</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>
\biggl[ \frac{a^2}{a^2 + c^2} \biggr] \zeta_2 z_0 \, ,
</math>
  </td>
</tr>
</table>

which gives the ratio, <math>~y_0/z_0</math> in terms of known parameter values; and,

<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~
\biggl[ \frac{x_\mathrm{max}}{ y_\mathrm{max} } \biggr] ~\varphi  
</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\biggl[ \frac{a^2}{a^2 + b^2} \biggr] \zeta_3  \cos\theta  
- \biggl[ \frac{a^2}{a^2 + c^2} \biggr] \zeta_2 \sin\theta  \, ,
</math>
  </td>
</tr>
</table>
which in combination with the just-derived similar product relation can give expressions for both terms, namely,

<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~
\biggl[ \frac{x_\mathrm{max}}{ y_\mathrm{max} } \biggr]^2  
</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\biggl\{ \biggl[ \frac{a^2}{a^2 + b^2} \biggr] \zeta_3  \cos\theta  
- \biggl[ \frac{a^2}{a^2 + c^2} \biggr] \zeta_2 \sin\theta  \biggr\}\biggl[ \frac{a^2 + b^2}{b^2} \biggr] \frac{\cos\theta}{\zeta_3} 
\, ,
</math>
  </td>
</tr>
</table>
and,

<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~
\varphi^2
</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\biggl\{ \biggl[ \frac{a^2}{a^2 + b^2} \biggr] \zeta_3   
- \biggl[ \frac{a^2}{a^2 + c^2} \biggr] \zeta_2 \tan\theta  \biggr\}\biggl[ \frac{b^2}{a^2 + b^2} \biggr] \zeta_3 
\, .
</math>
  </td>
</tr>
</table>

====Example Implementation====

<span id="Tipped">Drawing from</span> the set of [[#Case_I|''Case I'' parameters listed above]],  we will set a = 1, b = 1.25, c = 0.4703, &zeta;<sub>2</sub> = -2.2794, and &zeta;<sub>3</sub> = -1.9637.  This means that,

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

<tr>
  <td align="right">
<math>~\tan\theta </math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~- \frac{\zeta_2 }{ \zeta_3 } \biggl[ \frac{a^2 + b^2}{a^2 + c^2} \biggr]\frac{c^2}{b^2} = - 0.3448 ~~~\Rightarrow ~~~ \theta = -19.02^\circ
\, .</math>
  </td>
</tr>
</table>

We deduce as well that,
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~
\frac{y_0}{z_0}
</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>
\biggl[ \frac{a^2}{a^2 + c^2} \biggr] \zeta_2 \biggl\{ \biggl[ \frac{a^2}{a^2 + b^2} \biggr] \zeta_3  \cos\theta  - \biggl[ \frac{a^2}{a^2 + c^2} \biggr] \zeta_2 \sin\theta  \biggr\}^{-1}
= -1.8666 \cdot \{-0.7245 - 0.6084 \}^{-1} = 1.4004 
\, ;
</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~
\frac{x_\mathrm{max}}{ y_\mathrm{max} }   
</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\biggl\{
\biggl[ \biggl( \frac{a^2}{a^2 + b^2} \biggr) \zeta_3  \cos\theta  
- \biggl( \frac{a^2}{a^2 + c^2} \biggr) \zeta_2 \sin\theta  \biggr] \biggl[ \frac{a^2 + b^2}{b^2} \biggr] \frac{\cos\theta}{\zeta_3} \biggr\}^{1 / 2}
= \{ [-1.3329] \cdot (-0.7895) \}^{1 / 2} = \sqrt{1.0524} = 1.0259 
\, ;
</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~
\varphi
</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\biggl\{ 
\biggl[ \biggl( \frac{a^2}{a^2 + b^2} \biggr) \zeta_3   
- \biggl( \frac{a^2}{a^2 + c^2} \biggr) \zeta_2 \tan\theta  \biggr] \biggl[ \frac{b^2}{a^2 + b^2} \biggr] \zeta_3 \biggr\}^{1 / 2} = 1.2993
\, .
</math>
  </td>
</tr>
</table>

<table border="1" cellpadding="10" align="center" width="80%">
<tr><td align="left">
It is important to keep in mind that each Lagrangian fluid element will complete one full orbit in its "tipped" plane when  <math>~\varphi t = 2\pi</math>, where <math>~t</math> is in units of <math>~[\pi G \rho]^{-1 / 2}</math> and <math>~\varphi</math>  is in units of <math>~[\pi G \rho]^{+1 / 2}</math>.  Note, as well, that in our COLLADA animations, we have adopted the convention that <math>~t = 1</math> when <math>~\mathrm{TIME} = 4</math>.  Hence, one complete orbit will be completed when, 
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~\frac{1}{4} \cdot \mathrm{TIME}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\frac{2\pi}{\varphi} \, .</math>
  </td>
</tr>
</table>
That is to say, in the specific example being used here, on complete orbit is concluded when TIME = 19.343.
</td>
</tr></table>

Referring back to our [[#TippedPlane|earlier geometric prescription of a "tipped plane"]] we insert this new slope having m = tan&theta; = -0.3448, and choose a particular plane by setting a value of z<sub>0</sub>.  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:

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

<tr>
  <td align="right">
<math>~y </math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\frac{1}{\tan\theta} \biggl[ z - z_0 \biggr] \, ;</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~\frac{x}{a} </math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\biggl[ 1 - \biggl(\frac{z_0}{b \tan\theta}\biggr)^2 \biggl( \frac{z}{z_0} - 1 \biggr)^2 - \biggl( \frac{z}{c} \biggr)^2 \biggr]^{1 / 2} \, .</math>
  </td>
</tr>
</table>

For each specified value of z<sub>0</sub>, 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&theta; = -0.4310, this pair is given by the roots of the quadratic expression,

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

<tr>
  <td align="right">
<math>~1 </math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\biggl( \frac{z}{z_0} - 1 \biggr)^2 + \biggl( \frac{z}{c} \biggr)^2  </math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~\Rightarrow ~~~ 0 </math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~z^2\biggl[ \frac{1}{z_0^2} + \frac{1}{c^2} \biggr] -\frac{2z}{z_0} </math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~z \biggl\{ z - \biggl[ \frac{2 z_0 c^2}{(c^2 + z_0^2)} \biggr] \biggr\} \biggl[ \frac{(c^2 + z_0^2)}{c^2 z_0^2} \biggr] \, .</math>
  </td>
</tr>
</table>

The roots are, then, <math>~z_\mathrm{max} = 0</math>, and

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

<tr>
  <td align="right">
<math>~z_\mathrm{min}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\frac{2 z_0 c^2}{(c^2 + z_0^2)} = -0.468514 \, .</math>
  </td>
</tr>
</table>

<span id="ExampleTrajectories">The middle three columns</span> (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.
<table border="1" cellpadding="8" align="center">
<tr>
  <th align="center" colspan="11"><math>~m = \tan\theta = -0.3448 ~~~ \Rightarrow ~~~ \theta = -19.02^\circ</math></th>
</tr>
<tr>
  <td align="center" colspan="3">z_0 = -0.2500</td>
  <th align="center" rowspan="18">&nbsp;</th>
  <td align="center" colspan="3">z_0 = b tan&theta; = -0.4310</td>
  <th align="center" rowspan="18">&nbsp;</th>
  <td align="center" colspan="3">z_0 = -0.6000</td>
</tr>
<tr>
  <th align="center">z</th>
  <th align="center">y</th>
  <th align="center">|x|</th>
  <th align="center">z</th>
  <th align="center">y</th>
  <th align="center">|x|</th>
  <th align="center">z</th>
  <th align="center">y</th>
  <th align="center">|x|</th>
</tr>
<tr>
  <td align="center">0.1564</td><td align="center">-1.1788</td><td align="center">0.0000</td>
  <td align="center" bgcolor="lightblue">0</td>  <td align="center" bgcolor="lightblue">-1.2500</td>  <td align="center" bgcolor="lightblue">0</td> 
  <td align="center">-0.2182</td><td align="center">-1.1073</td><td align="center">0.0000</td>
</tr>
<tr>
  <td align="center">0.1147</td><td align="center">-1.0577</td><td align="center">0.4739</td>
  <td align="center" bgcolor="lightblue">-0.02</td>  <td align="center" bgcolor="lightblue">- 1.1920</td>  <td align="center" bgcolor="lightblue">0.2981</td> 
  <td align="center">-0.2336</td><td align="center">-1.0626</td><td align="center">0.1749</td>
</tr>
<tr>
  <td align="center">0.0729</td><td align="center">-0.9366</td><td align="center">0.6439</td>
  <td align="center" bgcolor="lightblue">-0.05</td>  <td align="center" bgcolor="lightblue">-1.1050</td>  <td align="center" bgcolor="lightblue">0.4553</td>
  <td align="center">-0.2490</td><td align="center">-1.0179</td><td align="center">0.2377</td>
</tr>
<tr>
  <td align="center">0.0312</td><td align="center">-0.8155</td><td align="center">0.7550</td>
  <td align="center" bgcolor="lightblue">-0.08</td>  <td align="center" bgcolor="lightblue">-1.0180</td>  <td align="center" bgcolor="lightblue">0.5548</td>
  <td align="center">-0.2645</td><td align="center">-0.9732</td><td align="center">0.2787</td>
</tr>
<tr>
  <td align="center">-0.0106</td><td align="center">-0.6943</td><td align="center">0.8312</td>
  <td align="center" bgcolor="lightblue">-0.1</td>  <td align="center" bgcolor="lightblue">-0.9600</td>  <td align="center" bgcolor="lightblue">0.6042</td>
  <td align="center">-0.2799</td><td align="center">-0.9285</td><td align="center">0.3068</td>
</tr>
<tr>
  <td align="center">-0.0524</td><td align="center">-0.5732</td><td align="center">0.8817</td>
  <td align="center" bgcolor="lightblue">-0.125</td>  <td align="center" bgcolor="lightblue">-0.8875</td>  <td align="center" bgcolor="lightblue">0.6521</td> 
  <td align="center">-0.2953</td><td align="center">-0.8838</td><td align="center">0.3255</td>
</tr>
<tr>
  <td align="center">-0.0941</td><td align="center">-0.4521</td><td align="center">0.9106</td>
  <td align="center" bgcolor="lightblue">-0.15</td>  <td align="center" bgcolor="lightblue">-0.8150</td>  <td align="center" bgcolor="lightblue">0.6879</td> 
  <td align="center">-0.3107</td><td align="center">-0.8390</td><td align="center">0.3361</td>
</tr>
<tr>
  <td align="center">-0.1359</td><td align="center">-0.3310</td><td align="center">0.9200</td>
  <td align="center" bgcolor="lightblue">-0.175</td>  <td align="center" bgcolor="lightblue">-0.7424</td>  <td align="center" bgcolor="lightblue">0.7133</td> 
  <td align="center">-0.3261</td><td align="center">-0.7943</td><td align="center">0.3396</td>
</tr>
<tr>
  <td align="center">-0.1776</td><td align="center">-0.2099</td><td align="center">0.9106</td>
  <td align="center" bgcolor="lightblue">-0.2</td>  <td align="center" bgcolor="lightblue">- 0.6699</td>  <td align="center" bgcolor="lightblue">0.7293</td>
  <td align="center">-0.3415</td><td align="center">-0.7496</td><td align="center">0.3361</td>
</tr>
<tr>
  <td align="center">-0.2194</td><td align="center">-0.0887</td><td align="center">0.8817</td>
  <td align="center" bgcolor="lightblue">-0.25</td>  <td align="center" bgcolor="lightblue">- 0.5249</td>  <td align="center" bgcolor="lightblue">0.7356</td>
  <td align="center">-0.3570</td><td align="center">-0.7049</td><td align="center">0.3255</td>
</tr>
<tr>
  <td align="center">-0.2612</td><td align="center">0.0324</td><td align="center">0.8312</td>
  <td align="center" bgcolor="lightblue">-0.30</td>  <td align="center" bgcolor="lightblue">- 0.3799</td>  <td align="center" bgcolor="lightblue">0.7076</td> 
  <td align="center">-0.3724</td><td align="center">-0.6602</td><td align="center">0.3068</td>
</tr>
<tr>
  <td align="center">-0.3029</td><td align="center">0.1535</td><td align="center">0.7550</td>
  <td align="center" bgcolor="lightblue">-0.35</td>  <td align="center" bgcolor="lightblue">- 0.2349</td>  <td align="center" bgcolor="lightblue">0.6410</td> 
  <td align="center">-0.3878</td><td align="center">-0.6155</td><td align="center">0.2787</td>
</tr>
<tr>
  <td align="center">-0.3447</td><td align="center">0.2746</td><td align="center">0.6439</td>
  <td align="center" bgcolor="lightblue">-0.4</td>  <td align="center" bgcolor="lightblue">- 0.0899</td>  <td align="center" bgcolor="lightblue">0.5210</td>
  <td align="center">-0.4032</td><td align="center">-0.5708</td><td align="center">0.2377</td>
</tr>
<tr>
  <td align="center">-0.3865</td><td align="center">0.3957</td><td align="center">0.4739</td>
  <td align="center" bgcolor="lightblue">-0.43</td>  <td align="center" bgcolor="lightblue">- 0.0029</td>  <td align="center" bgcolor="lightblue">0.4050</td> 
  <td align="center">-0.4186</td><td align="center">-0.5261</td><td align="center">0.1749</td>
</tr>
<tr>
  <td align="center">-0.4282</td><td align="center">0.5169</td><td align="center">0.0</td>
  <td align="center" bgcolor="lightblue">-0.46</td>  <td align="center" bgcolor="lightblue">+0.0841</td>  <td align="center" bgcolor="lightblue">0.1950</td> 
  <td align="center">-0.4340</td><td align="center">-0.4814</td><td align="center">0.0000</td>
</tr>
<tr>
  <td align="center" colspan="3">&nbsp;</td>
  <td align="center" bgcolor="lightblue">-0.4685</td>  <td align="center" bgcolor="lightblue">+0.1088</td>  <td align="center" bgcolor="lightblue">0.0081</td> 
  <td align="center" colspan="3">&nbsp;</td>
</tr>
</table>


More generally, letting <math>~q \equiv z_0/(b\tan\theta)</math>, the pair of limiting values of z are given by the roots of the quadratic expression,

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

<tr>
  <td align="right">
<math>~1 </math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~q^2 \biggl( \frac{z}{z_0} - 1 \biggr)^2 + \biggl( \frac{z}{c} \biggr)^2  </math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~\Rightarrow ~~~0 </math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~ z^2 \biggl[ \frac{ q^2c^2 + z_0^2)}{z_0^2 q^2 c^2} \biggr]  - \frac{2z}{z_0}    + \biggl[ 1  - \frac{1}{q^2} \biggr] \, .</math>
  </td>
</tr>
</table>

That is,

<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>~ \biggl[ \frac{z_0^2 q^2 c^2}{ 2(q^2 c^2 + z_0^2)} \biggr]\biggl\{ \frac{2}{z_0}  \pm   \biggl[\frac{4}{z_0^2} -  \frac{ 4 (1  - q^{-2})(q^2 c^2 + z_0^2)}{z_0^2 q^2 c^2} \biggr]^{1 / 2} \biggr\} </math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~ \biggl[ \frac{z_0 q^2 c^2}{ (q^2 c^2 + z_0^2)} \biggr]\biggl\{ 1  \pm   \biggl[1 -  \frac{ (1  - q^{-2})(q^2 c^2 + z_0^2)}{ q^2 c^2} \biggr]^{1 / 2} \biggr\} \, .</math>
  </td>
</tr>
</table>
And we see that the constraint set on z<sub>0</sub> is given by the condition,

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

<tr>
  <td align="right">
<math>~1</math>
  </td>
  <td align="center">
<math>~\ge</math>
  </td>
  <td align="left">
<math>~\frac{ (1  - q^{-2})(q^2c^2 + z_0^2)}{q^2c^2} </math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~\Rightarrow ~~~ q^2 c^2</math>
  </td>
  <td align="center">
<math>~\ge</math>
  </td>
  <td align="left">
<math>~(1  - q^{-2})(q^2c^2 + z_0^2) </math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~\ge</math>
  </td>
  <td align="left">
<math>~(q^2 - 1)(c^2 + b^2\tan^2\theta) </math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~\Rightarrow ~~~ 0</math>
  </td>
  <td align="center">
<math>~\ge</math>
  </td>
  <td align="left">
<math>~ (q^2 - 1)b^2\tan^2\theta - c^2</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~\ge</math>
  </td>
  <td align="left">
<math>~ z_0^2 - c^2 - b^2\tan^2\theta</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~\Rightarrow ~~~ z_0^2</math>
  </td>
  <td align="center">
<math>~\le</math>
  </td>
  <td align="left">
<math>~ c^2 + b^2\tan^2\theta \, .</math>
  </td>
</tr>
</table>

Therefore, for this example Type I ellipsoidal configuration, z<sub>0</sub> must lie within the range,
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~-0.650165</math>
  </td>
  <td align="center">
<math>~\le z_0 \le</math>
  </td>
  <td align="left">
<math>~+0.650165 \, .</math>
  </td>
</tr>
</table>
In the [[#ExampleTrajectories|above table]], we have detailed the loci of points along the surface trajectories that correspond to the values, z<sub>0</sub> = - 0.2500 (the leftmost three columns of the table) and z<sub>0</sub> = -0.6000 (the rightmost three columns of the table).  The limiting values of z for these to choices of z<sub>0</sub> are, respectively, [z<sub>min</sub> = -0.4283, z<sub>max</sub> = +0.1564] and  [z<sub>min</sub> = -0.4340, z<sub>max</sub> = -0.2182].

<span id="Figure3">Using COLLADA,</span> 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 &#8212; see the [[#explanation| explanation accompanying Figure 2, above]] &#8212; the ellipsoid is tilted about the x-coordinate axis at an angle of 61.25&deg; 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  [[#ExampleTrajectories|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<sub>0</sub> = -0.25, the smallest corresponds to our choice of z<sub>0</sub> = -0.60, and the one of intermediate size correspond to our choice of z<sub>0</sub> = -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 &theta; = -19.02&deg; to the ''untilted'' equatorial, x-y plane of the purple ellipsoid.

<div align="center">
<table border="1" align="center" cellpadding="8">
<tr>
  <th align="center">Figure 3a</th>
  <th align="center">Figure 3b</th>
</tr>
<tr>
  <td align="left" bgcolor="lightgrey">
[[File:B125c470B.cropped.png|400px|EFE Model b1.25c470B]]
  </td>
  <td align="left" bgcolor="lightgrey">
[[File:B125c470A.cropped.png|400px|EFE Model b1.25c470A]]
  </td>
</tr>
<tr>
  <th align="center" colspan="2">Figure 3c</th>
</tr>
<tr>
  <td align="center" bgcolor="white" colspan="2">
[[File:ProjectedOrbitsFlipped2.png|600px|EFE Model b1.25c470]]
  </td>
</tr>
</table>
</div>

As we have created each orbit in the xml-based COLLADA file, we have first manually typed in the (TIME, x, y, z) coordinates of each orbit that has a ''negative'' value of z<sub>0</sub>.  Specifically, using the Figure 3c projection as a guide/reference, we typed in coordinates for the following orbits:
<ul>
  <li><b>FLUID ELEMENT 1:</b>&nbsp; <math>~z_0 = -0.6000</math></li>
</ul>

<table border="1" cellpadding="8" align="center">
<tr>
  <td align="center">nstep</td>
  <td align="center">TIME</td>
  <td align="center">x</td>
  <td align="center">y</td>
  <td align="center">z</td>
</tr>
<tr>
  <td align="center">1</td>
  <td align="center">0.000</td>
  <td align="center">0.000</td>
  <td align="center">-1.107</td>
  <td align="center">-0.218</td>
</tr>
<tr>
  <td align="center">26</td>
  <td align="center">4.836</td>
  <td align="center">-0.340</td>
  <td align="center">-0.794</td>
  <td align="center">-0.326</td>
</tr>
<tr>
  <td align="center">51</td>
  <td align="center">9.672</td>
  <td align="center">0.000</td>
  <td align="center">-0.481</td>
  <td align="center">-0.434</td>
</tr>
<tr>
  <td align="center">76</td>
  <td align="center">14.507</td>
  <td align="center">+0.340</td>
  <td align="center">-0.794</td>
  <td align="center">-0.326</td>
</tr>
<tr>
  <td align="center">101</td>
  <td align="center">19.343</td>
  <td align="center">0.000</td>
  <td align="center">-1.107</td>
  <td align="center">-0.218</td>
</tr>
</table>

<ul>
  <li><b>FLUID ELEMENT 2:</b>&nbsp; <math>~z_0 = -0.4300</math></li>
</ul>

=See Also=

* Discussion of [[ThreeDimensionalConfigurations/RiemannStype|Ou's Riemann-Like Ellipsoids]]

* [[ThreeDimensionalConfigurations/MeetsCOLLADAandOculusRiftS|Riemann Meets COLLADA &amp; Oculus Rift S]]
** [[Appendix/Ramblings/VirtualReality#Virtual_Reality_and_3D_Printing|Virtual Reality and 3D Printing]]
** [[Appendix/Ramblings/OculusRiftS|Success Importing Animated Scene into Oculus Rift S]]
** [[Appendix/Ramblings/RiemannMeetsOculus|Carefully (Re)Build Riemann Type S Ellipsoids Inside Oculus Rift Environment]]: Example <b>(b/a, c/a) = (0.41, 0.385)</b>
** Other Example S-type Riemann Ellipsoids:
*** <b>[[Appendix/Ramblings/RiemannB90C333|(b/a, c/a) = (0.90, 0.333)]]</b>
*** <b>[[Appendix/Ramblings/RiemannB74C692|(b/a, c/a) = (0.74, 0.692)]]</b>
*** <b>[[Appendix/Ramblings/RiemannB28C256|(b/a, c/a) = (0.28, 0.256)]]</b>


* [[ThreeDimensionalConfigurations/DescriptionOfRiemannTypeI|Description of Riemann Type I Ellipsoids]]
* [[ThreeDimensionalConfigurations/RiemannTypeI#Riemann_Type_1_Ellipsoids|Riemann Type 1 Ellipsoids]] (old introduction)
* [[ThreeDimensionalConfigurations/Challenges#Challenges_Constructing_Ellipsoidal-Like_Configurations|Construction Challenges (Pt. 1)]]
* [[ThreeDimensionalConfigurations/ChallengesPt2|Construction Challenges (Pt. 2)]]
* [[ThreeDimensionalConfigurations/ChallengesPt3|Construction Challenges (Pt. 3)]]
* [[ThreeDimensionalConfigurations/ChallengesPt4|Construction Challenges (Pt. 4)]]
* [[ThreeDimensionalConfigurations/ChallengesPt5|Construction Challenges (Pt. 5)]]
* [[ThreeDimensionalConfigurations/ChallengesPt6|Construction Challenges (Pt. 6)]]

==Chandrasekhar's Detailed Analysis==
* [https://books.google.com/books?id=vo9VAAAAYAAJ&pg=PR6&dq=Gesammelte+Mathematische+Werke+B.+Riemann&source=gbs_selected_pages&cad=3#v=onepage&q=Gesammelte%20Mathematische%20Werke%20B.%20Riemann&f=false 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"
* [https://ui.adsabs.harvard.edu/abs/1965ApJ...142..890C/abstract S. Chandrasekhar (1965)], ApJ, 142, 890 - 961.  ''The Equilibrum and the Stability of the Riemann Ellipsoids. I.'' &#8212; This work is referenced as [[Appendix/References#Appendix_of_EFE|Paper XXV in EFE]] and focuses on S-type Riemann ellipsoids.
* [https://ui.adsabs.harvard.edu/abs/1966ApJ...145..842C/abstract S. Chandrasekhar (1966)], ApJ, 145, 842 - 877.  ''The Equilibrum and the Stability of the Riemann Ellipsoids. II.'' &#8212; This work is referenced as [[Appendix/References#Appendix_of_EFE|Paper XXVIII in EFE]] and focuses on Riemann ellipsoids of Types I, II and III.

Potentially relevant citations:

* [https://ui.adsabs.harvard.edu/abs/1997Ap%26SS.254..269C/abstract T. T. Chia &amp; S. Y. Pung (1997)] Astrophysics and Space Science, 254 (issue 2), pp. 269 - 294, "Dynamical Behaviour of Compressible Homogeneous Uniformly Rotating Ellipsoids with Nonparallel Angular Velocity and Vorticity"
* [https://ui.adsabs.harvard.edu/abs/1989PhLB..218..386R/abstract D. J. Rowe, P. Rochford, &amp; J. Carvalho (1989)], Physics Letters B, 218 (issue 4), pp. 386 - 390, "Vorticity as intrinsic spin in the nuclear collective model.

==Finite-Amplitude Oscillations==

* [https://ui.adsabs.harvard.edu/abs/1967ApJ...149..145R/abstract L. F. Rossner (1967)], ApJ, 149, 145. ''The Finite-Amplitude Oscillations of the Maclaurin Spheroids''  &#8212; This work is referenced as [[Appendix/References#Appendix_of_EFE|Paper XXXVIII in EFE]].
* [https://ui.adsabs.harvard.edu/abs/1968ApJ...152..523F/abstract M. Fujimoto (1968)], ApJ, 152, 523. ''Gravitational Collapse of Rotating Gaseous Ellipsoids''
* [https://ui.adsabs.harvard.edu/abs/1995Ap%26SS.229..215C/abstract T. T. Chia &amp; 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''
* [https://ui.adsabs.harvard.edu/abs/1997Ap%26SS.254..269C/abstract T. T. Chia &amp; 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==
* [https://ui.adsabs.harvard.edu/abs/1970ApJ...162...97H/abstract C. Hunter (1970)], ApJ, 162, 97 - 103.  ''The Disklike Riemann Ellipsoids.''

==Other==

* [[Apps/MaclaurinSpheroids#Maclaurin_Spheroids_.28axisymmetric_structure.29|Properties of Maclaurin Spheroids]]
* [[Apps/MaclaurinSpheroids/GoogleBooks#Excerpts_from_A_Treatise_of_Fluxions|Excerpts from Maclaurin's (1742) ''A Treatise of Fluxions'']]
* [https://ui.adsabs.harvard.edu/abs/2020PhRvL.124e2501S/abstract N. Sensharma, ''et al.'' (2020)], Physical Review Letters, Volume 124, Issue 5, article id.052501.  ''Longitudinal Wobbling Motion in <sup>187</sup>Au.''

{{ SGFfooter }}