Difference between revisions of "User:Tohline/ThreeDimensionalConfigurations/ChallengesPt4"

From VistrailsWiki
Jump to navigation Jump to search
Line 868: Line 868:
</tr>
</tr>
</table>
</table>
In order for the <math>~\boldsymbol{k}'</math> component to be zero in the tipped plane, we must choose the tipping angle such that,
<span id="ThetaDef">In order</span> for the <math>~\boldsymbol{k}'</math> component to be zero in the tipped plane, we must choose the tipping angle such that,
<table border="0" cellpadding="5" align="center">
<table border="0" cellpadding="5" align="center">



Revision as of 22:18, 25 May 2021

Challenges Constructing Ellipsoidal-Like Configurations (Pt. 4)

This chapter extends the accompanying chapters titled, Construction Challenges (Pt. 1), (Pt. 2), and (Pt. 3). The focus here is on firming up our understanding of the relationships between various "tilted" Cartesian coordinate frames.

Whitworth's (1981) Isothermal Free-Energy Surface
|   Tiled Menu   |   Tables of Content   |  Banner Video   |  Tohline Home Page   |

The Plan

Intersection Expression

STEP #1

First, we present the mathematical expression that describes the intersection between the surface of an ellipsoid and a plane having the following properties:

  • The plane cuts through the ellipsoid's z-axis at a distance, <math>~z_0</math>, from the center of the ellipsoid;
  • The line of intersection is parallel to the x-axis of the ellipsoid; and,
  • The line that is perpendicular to the plane and passes through the z-axis at <math>~z_0</math> is tipped at an angle, <math>~\theta</math>, to the z-axis.

As is illustrated in Figure 1, we will use the line referenced in this third property description to serve as the z'-axis of a Cartesian grid that is tipped at the angle, <math>~\theta</math>, with respect to the body frame; and we will align the x' axis with the x-axis, so it should be clear that the z'-axis lies in the y-z plane of the ellipsoid.

Figure 1

<math>~x</math>

<math>~=</math>

<math>~x' \, ,</math>

<math>~y</math>

<math>~=</math>

<math>~ y' \cos\theta - z'\sin\theta \, ,</math>

<math>~(z - z_0)</math>

<math>~=</math>

<math>~ z' \cos\theta + y'\sin\theta \, .</math>

Primed Coordinates

<math>~x'</math>

<math>~=</math>

<math>~x \, ,</math>

<math>~y'</math>

<math>~=</math>

<math>~ y \cos\theta + (z - z_0) \sin\theta \, ,</math>

<math>~z'</math>

<math>~=</math>

<math>~ (z-z_0) \cos\theta - y \sin\theta \, .</math>


As has been shown in our accompanying discussion, we obtain the following,

Intersection Expression

<math>~1 - \frac{x^2}{a^2} </math>

<math>~=</math>

<math>~y^2 \biggl[\frac{c^2 + b^2\tan^2\theta}{b^2c^2} \biggr] + y \biggl[ \frac{2z_0 \tan\theta}{c^2} \biggr] + \frac{z_0^2}{c^2} \, , </math>

as long as z0 lies within the range,

<math>~z_0^2</math>

<math>~\le</math>

<math>~c^2 + b^2\tan^2\theta \, .</math>

Rewriting this "intersection expression" in terms of the tipped (primed) coordinate frame gives us,

<math>~1 - \frac{(x')^2}{a^2} </math>

<math>~=</math>

<math>~(y' \cos\theta - z' \sin\theta)^2 \biggl[\frac{c^2 + b^2\tan^2\theta}{b^2c^2} \biggr] + (y' \cos\theta - z' \sin\theta) \biggl[ \frac{2z_0 \tan\theta}{c^2} \biggr] + \frac{z_0^2}{c^2} \, . </math>

STEP #2

As viewed from the tipped coordinated frame, the curve that is identified by this intersection should be an

Off-Center Ellipse

<math>~1</math>

<math>~=</math>

<math>~\biggl[\frac{x'}{x_\mathrm{max}} \biggr]^2 + \biggl[\frac{y' - y_c}{y_\mathrm{max}} \biggr]^2 </math>

 

<math>~=</math>

<math>~\biggl[\frac{x'}{x_\mathrm{max}} \biggr]^2 + \biggl[\frac{(y')^2 - 2y' y_c + y_c^2}{y^2_\mathrm{max}} \biggr] \, ,</math>

that lies in the x'-y' plane — that is, <math>~z' = 0</math>. Let's see if the intersection expression can be molded into this form.

<math>~1 - \frac{z_0^2}{c^2} - \frac{(x')^2}{a^2} </math>

<math>~=</math>

<math>~(y')^2 \biggl[\frac{c^2 + b^2\tan^2\theta}{b^2c^2} \biggr]\cos^2\theta + 2y' \biggl[ \frac{z_0 \sin\theta}{c^2} \biggr] </math>

 

<math>~=</math>

<math>~\biggl[\frac{c^2 + b^2\tan^2\theta}{b^2c^2} \biggr]\cos^2\theta \biggl\{ (y')^2 - 2y' \biggl[ \frac{-z_0 \sin\theta}{c^2 \cos^2\theta} \biggr]\biggl[\frac{b^2c^2}{c^2 + b^2\tan^2\theta} \biggr] \biggr\}</math>

 

<math>~=</math>

<math>~\kappa^2 \biggl[ (y')^2 - 2y' \underbrace{\biggl( \frac{-z_0 \sin\theta}{c^2 \kappa^2} \biggr)}_{y_c} \biggr] \, ,</math>

RESULT 3
(same as Result 1, but different from Result 2, below)

<math>~\frac{y_c}{z_0}</math>

<math>~=</math>

<math>~ -\frac{\sin\theta}{c^2\kappa^2} </math>

where,

<math>~\kappa^2</math>

<math>~\equiv</math>

<math>~ \frac{c^2 \cos^2\theta + b^2 \sin^2\theta}{b^2c^2} \, . </math>

Dividing through by <math>~\kappa^2</math>, then adding <math>~y_c^2</math> to both sides gives,

<math>~(y')^2 - 2y' y_c + y_c^2</math>

<math>~=</math>

<math>~\underbrace{\biggl[ \frac{1}{\kappa^2} - \frac{z_0^2}{c^2 \kappa^2} + y_c^2 \biggr]}_{y^2_\mathrm{max}} - \frac{(x')^2}{a^2\kappa^2} \, .</math>

Finally, we have,

<math>\frac{1}{y^2_\mathrm{max}} \biggl[ (y')^2 - 2y' y_c + y_c^2 \biggr]</math>

<math>~=</math>

<math>~1 - (x')^2 \underbrace{\biggl[ \frac{1}{a^2\kappa^2 y_\mathrm{max}^2} \biggr]}_{ 1/x^2_\mathrm{max} } \, .</math>

So … the intersection expression can be molded into the form of an off-center ellipse if we make the following associations:

<math>~\frac{y_c}{z_0}</math>

<math>~=</math>

<math>~- \frac{\sin\theta}{c^2 \kappa^2} \, ,</math>

<math>~y_\mathrm{max}^2</math>

<math>~=</math>

<math>~\frac{1}{\kappa^2}\biggl[ 1 - \frac{z_0^2}{c^2 } - \frac{z_0 \sin\theta}{c^2} \biggr] \, ,</math>

<math>~x_\mathrm{max}^2</math>

<math>~=</math>

<math>~a^2 \biggl[ 1 - \frac{z_0^2}{c^2 } - \frac{z_0 \sin\theta}{c^2} \biggr] \, .</math>

Note as well that,

<math>~\biggl( \frac{x_\mathrm{max}}{y_\mathrm{max}}\biggr)^2</math>

<math>~=</math>

<math>~a^2\kappa^2 = \frac{a^2}{b^2 c^2} \biggl[ c^2 \cos^2\theta + b^2 \sin^2\theta \biggr] \, .</math>

Lagrangian Trajectory and Velocities

We presume that the off-center ellipse that is defined by the intersection expression identifies the trajectory of a Lagrangian fluid element. If this is the case, there are a couple of ways that the velocity — both the amplitude and its vector orientation — can be derived.

STEP #3

If the intersection expression identifies a Lagrangian trajectory, then the velocity vector must be tangent to the off-center ellipse at every location. At each <math>~(x', y')</math> coordinate location, the slope of the above-defined off-center ellipse is,

<math>~\frac{dy'}{dx'}</math>

<math>~=</math>

<math>~ \biggl( \frac{y_\mathrm{max}}{x_\mathrm{max}} \biggr)^2 \frac{x'}{(y_c - y')} \, . </math>

From this expression we deduce that the x'- and y'- components of the velocity vector are, respectively,

<math>~\frac{\boldsymbol{\hat\imath'} \cdot \boldsymbol{u'} }{ [\boldsymbol{u'}\cdot \boldsymbol{u'}]^{1 / 2} }</math>

<math>~=</math>

<math>~ \frac{1}{u'_0} \biggl( \frac{x_\mathrm{max}}{y_\mathrm{max}}\biggr) (y_c - y') \, , </math>

      and,      

<math>~\frac{\boldsymbol{\hat\jmath'} \cdot \boldsymbol{u'} }{ [\boldsymbol{u'}\cdot \boldsymbol{u'}]^{1 / 2} }</math>

<math>~=</math>

<math>~ \frac{1}{u'_0} \biggl( \frac{y_\mathrm{max}}{x_\mathrm{max}}\biggr) x' \, , </math>

where the position-dependent — and, hence also, the time-dependent — length scale,

<math>~u'_0</math>

<math>~\equiv</math>

<math>~ \biggl\{ \biggl[ \biggl( \frac{x_\mathrm{max}}{y_\mathrm{max}}\biggr) (y_c - y') \biggr]^2 + \biggl[ \frac{1}{|u'|} \biggl( \frac{y_\mathrm{max}}{x_\mathrm{max}}\biggr) x' \biggr]^2 \biggr\}^{1 / 2} </math>

 

<math>~=</math>

<math>~ \frac{1}{x_\mathrm{max} y_\mathrm{max}} \biggl[ x_\mathrm{max}^4 ( y_c - y')^2 + y_\mathrm{max}^4 (x')^2 \biggr]^{1 / 2} \, . </math>

STEP #4

As a function of time, the x'-y' coordinates and associated velocity components of each Lagrangian fluid element are given by the expressions,

<math>~x'</math>

<math>~=</math>

<math>~x_\mathrm{max}\cos(\dot\varphi t)</math>

      and,      

<math>~y' - y_c</math>

<math>~=</math>

<math>~y_\mathrm{max}\sin(\dot\varphi t) \, ,</math>

<math>~\dot{x}'</math>

<math>~=</math>

<math>~- x_\mathrm{max}~ \dot\varphi \cdot \sin(\dot\varphi t) = (y_c - y') \biggl[ \frac{x_\mathrm{max}}{y_\mathrm{max}} \biggr] \dot\varphi </math>

      and,      

<math>~\dot{y}' </math>

<math>~=</math>

<math>~y_\mathrm{max}~\dot\varphi \cdot \cos(\dot\varphi t) = x' \biggl[ \frac{y_\mathrm{max}}{x_\mathrm{max}}\biggr] \dot\varphi \, .</math>

This means that the (dimensional) velocity vector is,

<math>~\boldsymbol{u'}</math>

<math>~=</math>

<math>~ \boldsymbol{\hat\imath'} \dot{x}' + \boldsymbol{\hat\jmath'} \dot{y}' </math>

 

<math>~=</math>

<math>~ \boldsymbol{\hat\imath'} \biggl[ (y_c - y') \biggl( \frac{x_\mathrm{max}}{y_\mathrm{max}} \biggr) \dot\varphi \biggr] + \boldsymbol{\hat\jmath'} \biggl[ x' \biggl( \frac{y_\mathrm{max}}{x_\mathrm{max}}\biggr) \dot\varphi \biggr] </math>

<math>~\Rightarrow~~~\boldsymbol{u'} \cdot \boldsymbol{u'}</math>

<math>~=</math>

<math>~ \biggl[ (y_c - y') \biggl( \frac{x_\mathrm{max}}{y_\mathrm{max}} \biggr) \dot\varphi \biggr]^2 + \biggl[ x' \biggl( \frac{y_\mathrm{max}}{x_\mathrm{max}}\biggr) \dot\varphi \biggr]^2 </math>

 

<math>~=</math>

<math>~\frac{\dot\varphi^2}{x_\mathrm{max}^2 y_\mathrm{max}^2} \biggl[ (y_c - y')^2 x_\mathrm{max}^4 + (x')^2 y_\mathrm{max}^4 \biggr] </math>

 

<math>~=</math>

<math>~ (u_0')^2 \dot\varphi^2 \, . </math>

Riemann Flow

STEP #5

Example Type I
Ellipsoid
(see also)
<math>~\frac{b}{a} = \frac{a_2}{a_1}</math> 1.25
<math>~\frac{c}{a} = \frac{a_3}{a_1}</math> 0.4703
<math>~\Omega_2</math> 0.3639
<math>~\Omega_3</math> 0.6633
<math>~\tan^{-1} \biggl[ \frac{\Omega_3}{\Omega_2} \biggr]</math> 61.25°
<math>~\zeta_2</math> -2.2794
<math>~\zeta_3</math> -1.9637
<math>~\tan^{-1} \biggl[ \frac{\zeta_3}{\zeta_2} \biggr]</math> 40.74°
<math>~\beta_+</math> 1.13449 (1.13332)
<math>~\gamma_+</math> 1.8052

As we have summarized in an accompanying discussion of Riemann Type 1 ellipsoids — see also our separate discussion[EFE] provides an expression for the velocity vector of each fluid element, given its instantaneous body-coordinate position (x, y, z) = (x1, x2, x3) — see his Eq. (154), Chapter 7, §51 (p. 156). As viewed from the rotating body coordinate frame, the three component expressions are,

<math>~\dot{x} = u_1 = \boldsymbol{\hat\imath} \cdot \boldsymbol{u}</math>

<math>~=</math>

<math>~\biggl(\frac{a}{b}\biggr)^2 \gamma \Omega_3 y - \biggl(\frac{a}{c}\biggr)^2 \beta \Omega_2 z</math>

<math>~=</math>

<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>

<math>~\dot{y} = u_2 = \boldsymbol{\hat\jmath} \cdot \boldsymbol{u}</math>

<math>~=</math>

<math>~- \gamma \Omega_3 x</math>

<math>~=</math>

<math>~+\biggl[ \frac{b^2}{a^2 + b^2} \biggr] \zeta_3 x \, ,</math>

<math>~\dot{z} = u_3 = \boldsymbol{\hat{k}} \cdot \boldsymbol{u}</math>

<math>~=</math>

<math>~+ \beta \Omega_2 x</math>

<math>~=</math>

<math>~- \biggl[ \frac{c^2}{a^2 + c^2} \biggr] \zeta_2 x \, ,</math>

where,

<math>~\beta</math>

<math>~=</math>

<math>~ - \biggl[ \frac{c^2}{a^2 + c^2} \biggr] \frac{\zeta_2}{\Omega_2} </math>

      and,      

<math>~\gamma</math>

<math>~=</math>

<math>~ - \biggl[ \frac{b^2}{a^2 + b^2} \biggr] \frac{\zeta_3}{\Omega_3} \, . </math>


In order to transform Riemann's velocity vector from the body frame (unprimed) to the "tipped orbit" frame (primed coordinates), we use the following mappings of the three unit vectors:

<math>~\boldsymbol{\hat\imath}</math>

<math>~\rightarrow</math>

<math>~\boldsymbol{\hat\imath'} \, ,</math>

<math>~\boldsymbol{\hat\jmath}</math>

<math>~\rightarrow</math>

<math>~\boldsymbol{\hat\jmath'}\cos\theta - \boldsymbol{\hat{k}'}\sin\theta \, ,</math>

<math>~\boldsymbol{\hat{k}}</math>

<math>~\rightarrow</math>

<math>~\boldsymbol{\hat\jmath'}\sin\theta + \boldsymbol{\hat{k}'}\cos\theta \, .</math>

<math>~x</math>

<math>~\rightarrow</math>

<math>~x' \, ,</math>

<math>~y</math>

<math>~\rightarrow</math>

<math>~y' \cos\theta - z' \sin\theta \, ,</math>

<math>~z - z_0</math>

<math>~\rightarrow</math>

<math>~y' \sin\theta + z'\cos\theta \, .</math>

In the tipped frame, we find,

<math>~\boldsymbol{u'}_\mathrm{EFE}</math>

<math>~=</math>

<math>~ \boldsymbol{\hat\imath'} \biggl[ \biggl(\frac{a}{b}\biggr)^2 \gamma \Omega_3 ( y'\cos\theta - z'\sin\theta ) - \biggl(\frac{a}{c}\biggr)^2 \beta \Omega_2 ( z_0 + y'\sin\theta + z'\cos\theta ) \biggr] + [\boldsymbol{\hat{k}'}\sin\theta -\boldsymbol{\hat\jmath'}\cos\theta ] \biggl[ \gamma \Omega_3 x' \biggr] + [ \boldsymbol{\hat\jmath'}\sin\theta + \boldsymbol{\hat{k}'}\cos\theta] \biggl[ \beta \Omega_2 x' \biggr] </math>

 

<math>~=</math>

<math>~ \boldsymbol{\hat\imath'} \biggl[ \biggl(\frac{a}{b}\biggr)^2 \gamma \Omega_3 ( y'\cos\theta - z'\sin\theta ) - \biggl(\frac{a}{c}\biggr)^2 \beta \Omega_2 ( z_0 + y'\sin\theta + z'\cos\theta ) \biggr] + \boldsymbol{\hat\jmath'} \biggl[\beta \Omega_2 x' \cdot \sin\theta - \gamma \Omega_3 x' \cdot \cos\theta \biggr] + \boldsymbol{\hat{k}'}\biggl[ \beta \Omega_2 x' \cdot \cos\theta + \gamma \Omega_3 x' \cdot \sin\theta \biggr] \, . </math>

In order for the <math>~\boldsymbol{k}'</math> component to be zero in the tipped plane, we must choose the tipping angle such that,

<math>~\tan\theta</math>

<math>~=</math>

<math>~ -\frac{\beta\Omega_2}{\gamma \Omega_3} = -0.344793 ~~~\Rightarrow~~~ \theta = -19.0238^\circ \, . </math>

And if we examine the flow only in the tipped x'-y' plane, then we should set <math>~z' = -z_0/\cos\theta</math>. These two constraints lead to the velocity expression,

<math>~\boldsymbol{u'}_\mathrm{EFE}</math>

<math>~=</math>

<math>~ \boldsymbol{\hat\imath'} \biggl[ \biggl(\frac{a}{b}\biggr)^2 \gamma \Omega_3 ( y'\cos\theta + z_0\tan\theta ) - \biggl(\frac{a}{c}\biggr)^2 \beta \Omega_2 ( y'\sin\theta ) \biggr] + \boldsymbol{\hat\jmath'} \underbrace{\biggl[\beta \Omega_2 \sin\theta - \gamma \Omega_3 \cos\theta \biggr]}_{\dot\varphi y_\mathrm{max}/x_\mathrm{max}}x' \, . </math>

As we have indicated, this <math>~\boldsymbol{\hat\jmath'}</math> component will match our Step #4 velocity expression if,

<math>~\biggl(\frac{y_\mathrm{max}}{x_\mathrm{max}} \biggr)\dot\varphi</math>

<math>~=</math>

<math>~ \beta \Omega_2 \sin\theta - \gamma \Omega_3 \cos\theta </math>

 

<math>~=</math>

<math>~ \biggl[\frac{1}{1 + \tan^2\theta}\biggr]^{1 / 2} \biggl[\beta\Omega_2 \tan\theta - \gamma\Omega_3\biggr] </math>

 

<math>~=</math>

<math>~ -\biggl[\frac{\gamma^2 \Omega_3^2}{\gamma^2 \Omega_3^2 + \beta^2\Omega_2^2}\biggr]^{1 / 2} \biggl[ \frac{\beta^2\Omega_2^2}{\gamma \Omega_3} + \gamma\Omega_3\biggr] </math>

 

<math>~=</math>

<math>~ -~\biggl[ \beta^2\Omega_2^2 + \gamma^2\Omega_3^2\biggr]^{1 / 2} \, . </math>

Rewriting the <math>~\boldsymbol{\hat\imath'}</math> component, we have,

<math>~\biggl(\frac{a}{b}\biggr)^2 \gamma \Omega_3 ( y'\cos\theta + z_0\tan\theta ) - \biggl(\frac{a}{c}\biggr)^2 \beta \Omega_2 ( y'\sin\theta )</math>

<math>~=</math>

<math>~ y' \biggl[ \biggl(\frac{a}{b}\biggr)^2 \gamma \Omega_3 \cos\theta - \biggl(\frac{a}{c}\biggr)^2 \beta \Omega_2 \sin\theta \biggr] + \biggl(\frac{a}{b}\biggr)^2 \gamma \Omega_3 ( z_0\tan\theta ) </math>

 

<math>~=</math>

<math>~ y' \biggl[ \biggl(\frac{a}{b}\biggr)^2 \gamma \Omega_3 - \biggl(\frac{a}{c}\biggr)^2 \beta \Omega_2 \tan\theta \biggr]\biggl[\frac{1}{1 + \tan^2\theta}\biggr]^{1 / 2} - \biggl(\frac{a}{b}\biggr)^2 \beta \Omega_2 z_0 </math>

 

<math>~=</math>

<math>~ y' \biggl[ \biggl(\frac{a}{b}\biggr)^2 \gamma^2 \Omega_3^2 + \biggl(\frac{a}{c}\biggr)^2 \beta^2 \Omega_2^2 \biggr] \biggl[\frac{1}{\gamma^2 \Omega_3^2 + \beta^2\Omega_2^2}\biggr]^{1 / 2} - \biggl(\frac{a}{b}\biggr)^2 \beta \Omega_2 z_0 \, . </math>

Hence,

<math>~\boldsymbol{u'}_\mathrm{EFE}</math>

<math>~=</math>

<math>~ \boldsymbol{\hat\imath'} \biggl\{ y' \underbrace{\biggl[ \biggl(\frac{a}{b}\biggr)^2 \gamma^2 \Omega_3^2 + \biggl(\frac{a}{c}\biggr)^2 \beta^2 \Omega_2^2 \biggr] \biggl[\frac{1}{\gamma^2 \Omega_3^2 + \beta^2\Omega_2^2}\biggr]^{1 / 2}}_{-\dot\varphi x_\mathrm{max}/y_\mathrm{max}} - \biggl(\frac{a}{b}\biggr)^2 \beta \Omega_2 z_0\biggr\} + \boldsymbol{\hat\jmath'} \biggl[ \dot\varphi \biggl(\frac{y_\mathrm{max}}{x_\mathrm{max}} \biggr) \biggr]x' \, . </math>

That is to say, we need to set,

<math>~\dot\varphi \biggl( \frac{x_\mathrm{max}}{y_\mathrm{max}}\biggr)</math>

<math>~=</math>

<math>~-~ \biggl[ c^2 \gamma^2 \Omega_3^2 + b^2 \beta^2 \Omega_2^2 \biggr] \biggl[\frac{1}{\gamma^2 \Omega_3^2 + \beta^2\Omega_2^2}\biggr]^{1 / 2} \biggl( \frac{a^2}{b^2 c^2}\biggr) \, . </math>

When this is combined with the constraint set by the <math>~\boldsymbol{\hat\jmath'}</math> component, we find,

<math>~\dot\varphi^2 = \dot\varphi \biggl( \frac{x_\mathrm{max}}{y_\mathrm{max}}\biggr) \cdot \dot\varphi \biggl( \frac{y_\mathrm{max}}{x_\mathrm{max}}\biggr)</math>

<math>~=</math>

<math>~+ \biggl[ c^2 \gamma^2 \Omega_3^2 + b^2 \beta^2 \Omega_2^2 \biggr] \biggl( \frac{a^2}{b^2 c^2}\biggr) = (1.29930)^2 \, , </math>

and, hence,

<math>~\biggl( \frac{x_\mathrm{max}}{y_\mathrm{max}}\biggr)^2</math>

<math>~=</math>

<math>~ \biggl[ c^2 \gamma^2 \Omega_3^2 + b^2 \beta^2 \Omega_2^2 \biggr] \biggl[\frac{1}{\gamma^2 \Omega_3^2 + \beta^2\Omega_2^2}\biggr]^{1 / 2} \biggl( \frac{a^2}{b^2 c^2}\biggr) \biggl[ \beta^2\Omega_2^2 + \gamma^2\Omega_3^2\biggr]^{-1 / 2} </math>

 

<math>~=</math>

<math>~ \biggl[\frac{c^2 \gamma^2 \Omega_3^2 + b^2 \beta^2 \Omega_2^2 }{\gamma^2 \Omega_3^2 + \beta^2\Omega_2^2}\biggr] \biggl( \frac{a^2}{b^2 c^2}\biggr) = ( 1.02585 )^2 \, . </math>

Finally, then,

<math>~y_c \biggl( \frac{x_\mathrm{max}}{y_\mathrm{max}}\biggr) \dot\varphi</math>

<math>~=</math>

<math>~ -\biggl(\frac{a}{b}\biggr)^2 \beta \Omega_2 z_0 </math>

<math>~\Rightarrow ~~~ \frac{y_c}{z_0} </math>

<math>~=</math>

<math>~ -\biggl(\frac{a}{b}\biggr)^2 \beta \Omega_2 \biggl\{ \biggl( \frac{x_\mathrm{max}}{y_\mathrm{max}}\biggr) \dot\varphi \biggr\}^{-1} </math>

 

<math>~=</math>

<math>~ c^2 \beta \Omega_2 \biggl[ c^2 \gamma^2 \Omega_3^2 + b^2 \beta^2 \Omega_2^2 \biggr]^{-1} \biggl[\gamma^2 \Omega_3^2 + \beta^2\Omega_2^2\biggr]^{1 / 2} = 0.19823 \, . </math>

RESULT 2
(different from Result 1 and Result 3, above)

<math>~\frac{y_0}{z_0}</math>

<math>~=</math>

<math>~ -\frac{\sin\theta}{b^2\kappa^2} </math>

Reconcile Results Differences

Before calling upon any of Riemann's model parameters, from geometric considerations alone we should be able to determine exactly what the expression is for any off-center ellipse that results from slicing — at a tipped angle — the chosen ellipsoid.

Figure 1

<math>~x</math>

<math>~=</math>

<math>~x' \, ,</math>

<math>~y</math>

<math>~=</math>

<math>~ y' \cos\theta - z'\sin\theta \, ,</math>

<math>~(z - z_0)</math>

<math>~=</math>

<math>~ z' \cos\theta + y'\sin\theta \, .</math>

Primed Coordinates

<math>~x'</math>

<math>~=</math>

<math>~x \, ,</math>

<math>~y'</math>

<math>~=</math>

<math>~ y \cos\theta + (z - z_0) \sin\theta \, ,</math>

<math>~z'</math>

<math>~=</math>

<math>~ (z-z_0) \cos\theta - y \sin\theta \, .</math>

Body Frame

As has been shown in our accompanying discussion, we obtain the following,

Intersection Expression

<math>~1 - \frac{x^2}{a^2} </math>

<math>~=</math>

<math>~y^2 \biggl[\frac{c^2 + b^2\tan^2\theta}{b^2c^2} \biggr] + y \biggl[ \frac{2z_0 \tan\theta}{c^2} \biggr] + \frac{z_0^2}{c^2} \, , </math>

as long as z0 lies within the range,

<math>~z_0^2</math>

<math>~\le</math>

<math>~c^2 + b^2\tan^2\theta \, .</math>

Given the values of <math>~a, b, c, \theta, z_0</math>, we can immediately map this to the tipped plane to obtain the surface-intersection function, <math>~x'(y')</math>. Specifically we find,

<math>~x'</math>

<math>~=</math>

<math>~a \biggl\{ 1 - \biggl[ (y'\cos\theta - \cancelto{0}{z'} \sin\theta)^2 \biggl( \frac{c^2 + b^2\tan^2\theta}{b^2c^2} \biggr) + (y'\cos\theta - \cancelto{0}{z'} \sin\theta) \biggl( \frac{2z_0 \tan\theta}{c^2} \biggr) + \frac{z_0^2}{c^2} \biggr] \biggr\}^{1 / 2}</math>

 

<math>~=</math>

<math>~a \biggl\{ 1 - \biggl[ (y'\cos\theta)^2 \biggl( \frac{c^2 + b^2\tan^2\theta}{b^2c^2} \biggr) + y'\cos\theta \biggl( \frac{2z_0 \tan\theta}{c^2} \biggr) + \frac{z_0^2}{c^2} \biggr] \biggr\}^{1 / 2}</math>

Along the tipped <math>~y'</math> axis, two points will mark the ends of the x'-y' orbit; they are identified by the roots of this Intersection Expression when x = 0. That is, by the roots of,

<math>~0</math>

<math>~=</math>

<math>~y^2 \biggl[\frac{\kappa^2}{\cos^2\theta} \biggr] + y \biggl[ \frac{2z_0 \tan\theta}{c^2} \biggr] + \biggl[\frac{z_0^2}{c^2} -1 \biggr] \, ,</math>

where,

<math>~\kappa^2</math>

<math>~\equiv</math>

<math>~ \frac{c^2 \cos^2\theta + b^2 \sin^2\theta}{b^2c^2} \, . </math>

The roots are …

Scratch notes:

<math>~y_\mathrm{edge}</math>

<math>~=</math>

<math>~\frac{B}{2A} \biggl\{\pm \biggl[1 - \frac{4AC}{B^2} \biggr]^{1 / 2} -1 \biggr\}</math>

where,

<math>~ \frac{4AC}{B^2} = \frac{4\kappa^2}{\cos^2\theta} \biggl[\frac{z_0^2}{c^2} -1 \biggr] \biggl[ \frac{2z_0 \tan\theta}{c^2} \biggr]^{-2} = \frac{\kappa^2}{\cos^2\theta} \biggl[ \frac{c^2(z_0^2 - c^2)}{z_0^2 \tan^2\theta} \biggr] = \biggl[ \frac{c^2 \kappa^2 (z_0^2 - c^2)}{z_0^2 \sin^2\theta} \biggr] </math>


<math>~y_\mathrm{edge}</math>

<math>~=</math>

<math>~\biggl[ \frac{z_0 \sin\theta \cos\theta}{c^2 \kappa^2} \biggr]\biggl\{\pm \biggl[1 - \frac{c^2 \kappa^2 (z_0^2 - c^2)}{z_0^2 \sin^2\theta} \biggr]^{1 / 2} -1 \biggr\}</math>

 

<math>~=</math>

<math>~\biggl[ \frac{\cos\theta}{c^2 \kappa^2} \biggr]\biggl\{\pm \biggl[z_0^2 \sin^2\theta + c^2 \kappa^2 (c^2 - z_0^2 ) \biggr]^{1 / 2} - z_0 \sin\theta \biggr\} \, .</math>

Hereafter we will use <math>~y_p</math> to denote the "plus" root, and <math>~y_m</math> to denote the "minus" root; that is,

<math>~y_p</math>

<math>~\equiv</math>

<math>~\biggl[ \frac{\cos\theta}{c^2 \kappa^2} \biggr]\biggl\{\biggl[z_0^2 \sin^2\theta + c^2 \kappa^2 (c^2 - z_0^2 ) \biggr]^{1 / 2} - z_0 \sin\theta \biggr\} \, ,</math>

<math>~y_m</math>

<math>~\equiv</math>

<math>~- \biggl[ \frac{\cos\theta}{c^2 \kappa^2} \biggr]\biggl\{\biggl[z_0^2 \sin^2\theta + c^2 \kappa^2 (c^2 - z_0^2 ) \biggr]^{1 / 2} + z_0 \sin\theta \biggr\} \, .</math>

  • If tanθ is positive …
    • Quadrant 1 ⇒ both sinθ and cosθ are positive.
    • Quadrant 3 ⇒ both sinθ and cosθ are negative.
  • If tanθ is negative …
    • Quadrant 2 ⇒ sinθ is positive but cosθ is negative.
    • Quadrant 4 ⇒ sinθ is negative while cosθ are positive.

Taking numerical values from our chosen Example Model — and using z0 = ± 0.25 as our 1st test cases — we have:

<math>~\tan\theta</math>

<math>~=</math>

<math>~ -\frac{\beta\Omega_2}{\gamma \Omega_3} = -0.344793 ~~~\Rightarrow~~~ \theta = -19.0238^\circ \, . </math>

<math>~\kappa^2</math>

<math>~\equiv</math>

<math>~ \frac{c^2 \cos^2\theta + b^2 \sin^2\theta}{b^2c^2} = 1.05238\, . </math>

<math>~\mathrm{TERM1}</math>

<math>~\equiv</math>

<math>~z_0^2 \sin^2\theta + c^2\kappa^2(c^2-z_0^2) =0.043577 </math>    (for z0 = ± 0.25).

z0 = +0.25 and θ in QUADRANT #4 (sinθ is negative while cosθ is positive):

<math>~y_p</math>

<math>~\equiv</math>

<math>~\biggl[ \frac{\cos\theta}{c^2 \kappa^2} \biggr]\biggl\{\biggl[ \mathrm{TERM1} \biggr]^{1 / 2} - z_0 \sin\theta \biggr\} = 1.1788 \, ,</math>

<math>~y_m</math>

<math>~\equiv</math>

<math>~- \biggl[ \frac{\cos\theta}{c^2 \kappa^2} \biggr]\biggl\{\biggl[ \mathrm{TERM1} \biggr]^{1 / 2} + z_0 \sin\theta \biggr\} = - 0.5169 \, .</math>

<math>~z_p</math>

<math>~\equiv</math>

<math>~-c \biggl[ 1 - \biggl( \frac{y_p}{b} \biggr)^2 \biggr]^{1 / 2} = -0.15645 \, ,</math>

<math>~z_m</math>

<math>~\equiv</math>

<math>~c \biggl[ 1 - \biggl( \frac{y_m}{b} \biggr)^2 \biggr]^{1 / 2} = 0.42821 \, ,</math>

<math>~|z_0|_\mathrm{max} = 0.63792</math>
<math>~z_0</math> <math>~y_p</math> <math>~z_p</math> <math>~y_m</math> <math>~z_m</math>
<math>~+ 0.6379</math> <math>~+ 0.8511</math> <math>~+ 0.3444</math> <math>~+ 0.8379</math> <math>~+ 0.3490</math>
<math>~+ 0.60</math> <math>~+ 1.1073</math> <math>~+ 0.2182</math> <math>~+ 0.4814</math> <math>~+ 0.4340</math>
<math>~+ 0.55</math> <math>~+ 1.1950</math> <math>~+ 0.1380</math> <math>~+ 0.2613</math> <math>~+ 0.4599</math>
<math>~+ 0.50</math> <math>~+ 1.2342</math> <math>~+ 0.0744</math> <math>~+ 0.0897</math> <math>~+ 0.4691</math>
<math>~+ 0.45</math> <math>~+ 1.2489</math> <math>~- 0.0194</math> <math>~- 0.0574</math> <math>~+ 0.4698</math>
<math>~+ 0.40</math> <math>~+ 1.2474</math> <math>~- 0.0301</math> <math>~- 0.1883</math> <math>~+ 0.4649</math>
<math>~+ 0.35</math> <math>~+ 1.2338</math> <math>~- 0.0754</math> <math>~- 0.3071</math> <math>~+ 0.4559</math>
<math>~+ 0.30</math> <math>~+ 1.2105</math> <math>~- 0.1174</math> <math>~- 0.4161</math> <math>~+ 0.4435</math>
+ 0.25 + 1.1788 - 0.1564 - 0.5169 + 0.4282
<math>~+ 0.20</math> <math>~+ 1.1399</math> <math>~- 0.1930</math> <math>~- 0.6103</math> <math>~+ 0.4104</math>
<math>~+ 0.15</math> <math>~+ 1.0943</math> <math>~- 0.2273</math> <math>~- 0.6971</math> <math>~+ 0.3904</math>
<math>~+ 0.10</math> <math>~+ 1.0426</math> <math>~- 0.2595</math> <math>~- 0.7778</math> <math>~+ 0.3682</math>
<math>~+ 0.05</math> <math>~+ 0.9849</math> <math>~- 0.2896</math> <math>~- 0.8525</math> <math>~+ 0.3439</math>
<math>~+ 0.00</math> <math>~+ 0.9216</math> <math>~- 0.3177</math> <math>~- 0.9216</math> <math>~+ 0.3177</math>
<math>~- 0.05</math> <math>~+ 0.8525</math> <math>~- 0.3439</math> <math>~- 0.9849</math> <math>~+ 0.2895</math>
<math>~- 0.10</math> <math>~+ 0.7778</math> <math>~- 0.3682</math> <math>~- 1.0426</math> <math>~+ 0.2595</math>
etc.
<math>~- 0.6379</math> <math>~- 0.8379</math> <math>~- 0.3490</math> <math>~- 0.8511</math> <math>~- 0.3444</math>


The y-coordinate of the center of the orbit will lie halfway between these two "edges". That is (numerical example given for the case of z0 = + 0.25),

<math>~y_\mathrm{center}</math>

<math>~=</math>

<math>~\frac{1}{2}\biggl[ y_p + y_m \biggr] = +0.3310 \, .</math>

Correspondingly, the vertical (z) location of the orbit center in the body frame is given by the expression,

<math>~z_\mathrm{center}</math>

<math>~=</math>

<math>~y_\mathrm{center}\tan\theta + z_0 = +0.1359 \, .</math>

The line that runs parallel to the x-axis — lies perpendicular to the y-z plane — and that passes through this "center" location intersects the surface of the ellipsoid at a value of x that is obtained by plugging <math>~y_\mathrm{center}</math> into the Intersection Expression. That is,

<math>~\frac{x_\mathrm{surf}}{a}</math>

<math>~=</math>

<math>~\biggl\{ 1 - y^2_\mathrm{center} \biggl[\frac{\kappa^2}{\cos^2\theta} \biggr] - y_\mathrm{center} \biggl[ \frac{2z_0 \tan\theta}{c^2} \biggr] - \frac{z_0^2}{c^2} \biggr\}^{1 / 2} = 0.92001 \, ;</math>

along the relevant orbit, this is also the maximum value of x.

Tipped Plane

We use the expressions in the right panel of Figure 1, above to transform (x, y, z) body coordinates to (x', y', z') tipped-frame coordiantes. Given that <math>~x' = x</math>, we appreciate that the x'-coordinate of the orbit center is zero, and the maximum value <math>~(x'_\mathrm{surf})</math> that is encountered along the relevant orbit is the same as the value of <math>~x_\mathrm{surf}</math> as determined in the body (unprimed) frame.

The z'-coordinate of the orbit center is,

<math>~z'_\mathrm{center}</math>

<math>~=</math>

<math>~ (z_\mathrm{center} - z_0)\cos\theta - y_\mathrm{center}\sin\theta </math>

 

<math>~=</math>

<math>~ (y_\mathrm{center}\tan\theta) \cos\theta - y_\mathrm{center}\sin\theta </math>

 

<math>~=</math>

<math>~ 0 \, . </math>

Finally, in the tipped plane, the maximum (subscript_p) and minimum (subscript_m) values of y' are,

<math>~y'_p</math>

<math>~=</math>

<math>~ (z_p - z_0)\sin\theta + y_p\cos\theta = 1.2469 \, , </math>

<math>~y'_m</math>

<math>~=</math>

<math>~ (z_m - z_0)\sin\theta + y_m\cos\theta = -0.5467 \, , </math>

<math>~\Rightarrow ~~~ y'_\mathrm{center}</math>

<math>~=</math>

<math>~ \frac{1}{2}\biggl[ y_p' + y_m'\biggr] = 0.3501\, . </math>

<math>~\Rightarrow ~~~ y'_\mathrm{surf}</math>

<math>~=</math>

<math>~ y'_p - y'_\mathrm{center} = 0.8968\, . </math>

So, in the tipped, x'-y' plane, the elliptical orbit that corresponds to z0 = + 0.25 should be described by the expression,

<math>~1</math>

<math>~=</math>

<math>~ \biggl[\frac{x'}{x_\mathrm{surf}} \biggr]^2 + \biggl[ \frac{(y' - y'_\mathrm{center} ) }{y'_\mathrm{surf}} \biggr]^2 \, . </math>


NOTE:

  • <math>~y'_\mathrm{center}/z_0 = 0.350096/0.25 = 1.40038 \, .</math>
  • <math>~x_\mathrm{surf}/y'_\mathrm{surf} = 0.92001/0.896819 = 1.0259 \, .</math>

At a minimum these numerical values should be compared to the Riemann model parameters computed in our earlier (Pt. 2) examination.

Various Coordinate Frames

Tipped Orbit Planes

Summary

In a separate discussion, we have shown that, as viewed from a frame that "tumbles" with the (purple) body of a Type 1 Riemann ellipsoid, each Lagrangian fluid element moves along an elliptical path in a plane that is tipped by an angle <math>~\theta</math> about the x-axis of the body. As viewed from the (primed) coordinates associated with this tipped plane, by definition, z' = constant and dz'/dt = 0, and the planar orbit is defined by the expression for an,

Off-Center Ellipse

<math>~1</math>

<math>~=</math>

<math>~\biggl[\frac{x'}{x_\mathrm{max}} \biggr]^2 + \biggl[\frac{y' - y_c(z')}{y_\mathrm{max}} \biggr]^2 \, .</math>

Tipped Orbit Frame (yellow, primed)

Tipped Orbital Planes

Given that b/a = 1.25 and c/a = 0.4703 for our chosen Example Type I Ellipsoid, we find that, <math>~\theta = - 1.18122 ~\mathrm{rad} = -67.68^\circ</math>.

Notice that the offset, <math>~y_c</math>, is a function of the tipped plane's vertical coordinate, <math>~z'</math>. As a function of time, the x'-y' coordinates and associated velocity components of each Lagrangian fluid element are given by the expressions,

<math>~x'</math>

<math>~=</math>

<math>~x_\mathrm{max}\cos(\dot\varphi t)</math>

      and,      

<math>~y' - y_c</math>

<math>~=</math>

<math>~y_\mathrm{max}\sin(\dot\varphi t) \, ,</math>

<math>~\dot{x}'</math>

<math>~=</math>

<math>~- x_\mathrm{max}~ \dot\varphi \cdot \sin(\dot\varphi t) = (y_c - y') \biggl[ \frac{x_\mathrm{max}}{y_\mathrm{max}} \biggr] \dot\varphi </math>

      and,      

<math>~\dot{y}' </math>

<math>~=</math>

<math>~y_\mathrm{max}~\dot\varphi \cdot \cos(\dot\varphi t) = x' \biggl[ \frac{y_\mathrm{max}}{x_\mathrm{max}}\biggr] \dot\varphi \, .</math>

We have determined that (numerical value given for our chosen example Type I ellipsoid),

<math>~\tan\theta</math>

<math>~=</math>

<math>~ - \frac{b^2 \beta \Omega_2}{c^2 \gamma \Omega_3} = -2.43573\, , </math>

This definition of tan(θ) is inconsistent with all others!

where, as has also been specified in an accompanying discussion,

<math>~\beta</math>

<math>~=</math>

<math>~ - \biggl[ \frac{c^2}{a^2 + c^2} \biggr] \frac{\zeta_2}{\Omega_2} = +1.13451 </math>

      and,      

<math>~\gamma</math>

<math>~=</math>

<math>~ - \biggl[ \frac{b^2}{a^2 + b^2} \biggr] \frac{\zeta_3}{\Omega_3} = +1.80518\, . </math>

We also have determined that,

<math>~\biggl[\frac{x_\mathrm{max}}{y_\mathrm{max}} \biggr]^4</math>

<math>~=</math>

<math>~ \frac{a^4 (c^4 \gamma^2 \Omega_3^2 + b^4 \beta^2 \Omega_2^2)}{b^4 c^4(\gamma^2\Omega_3^2 + \beta^2\Omega_2^2)} ~~~\Rightarrow ~~~\frac{x_\mathrm{max}}{y_\mathrm{max}} = 1.26218 \, , </math>

<math>~{\dot\varphi}^4 </math>

<math>~=</math>

<math>~ \frac{a^4}{b^4 c^4} \biggl(\gamma^2\Omega_3^2 + \beta^2\Omega_2^2 \biggr) (c^4 \gamma^2 \Omega_3^2 + b^4 \beta^2 \Omega_2^2) ~~~\Rightarrow ~~~ \dot\varphi = 1.59862\, , </math>

<math>~\frac{y_c}{z_0}</math>

<math>~=</math>

<math>~ -\sin\theta ~~~\Rightarrow~~~~ \frac{y_c}{z_0} = -0.92507 \, .</math>

Demonstration

In order to transform a vector from the "tipped orbit" frame (primed coordinates) to the "body" frame (unprimed), we use the following mappings of the three unit vectors:

<math>~\boldsymbol{\hat\imath'}</math>

<math>~\rightarrow</math>

<math>~\boldsymbol{\hat\imath} \, ,</math>

<math>~\boldsymbol{\hat\jmath'}</math>

<math>~\rightarrow</math>

<math>~\boldsymbol{\hat\jmath}\cos\theta + \boldsymbol{\hat{k}}\sin\theta \, ,</math>

<math>~\boldsymbol{\hat{k}'}</math>

<math>~\rightarrow</math>

<math>~-\boldsymbol{\hat\jmath}\sin\theta + \boldsymbol{\hat{k}}\cos\theta \, .</math>

Given that, by design in our "tipped orbit" frame, there is no vertical motion — that is, <math>~\dot{z}' = 0</math> — mapping the (primed coordinate) velocity to the body (unprimed) coordinate is particularly straightforward. Specifically,

<math>~\boldsymbol{u'}</math>

<math>~=</math>

<math>~ \boldsymbol{\hat\imath'} \dot{x}' + \boldsymbol{\hat\jmath'} \dot{y}' + \boldsymbol{\hat{k}'} \cancelto{0}{\dot{z}'} </math>

 

<math>~~~\rightarrow~~</math>

<math>~ \boldsymbol{\hat\imath} \dot{x}' + [\boldsymbol{\hat\jmath}\cos\theta + \boldsymbol{\hat{k}}\sin\theta] \dot{y}' </math>

 

<math>~=</math>

<math>~ \boldsymbol{\hat\imath} \biggl\{ (y_c - y') \biggl( \frac{x_\mathrm{max}}{y_\mathrm{max}}\biggr) \dot\varphi \biggr\} + [\boldsymbol{\hat\jmath}\cos\theta + \boldsymbol{\hat{k}}\sin\theta] \biggl\{ x' \biggl( \frac{y_\mathrm{max}}{x_\mathrm{max}}\biggr) \dot\varphi \biggr\} \, . </math>

Recognizing, as before, that the relevant coordinate mapping is,

<math>~x</math>

<math>~=</math>

<math>~x' \, ,</math>

<math>~y</math>

<math>~=</math>

<math>~ y' \cos\theta - z'\sin\theta \, ,</math>

<math>~(z - z_0)</math>

<math>~=</math>

<math>~ z' \cos\theta + y'\sin\theta \, .</math>

Primed Coordinates

<math>~x'</math>

<math>~=</math>

<math>~x \, ,</math>

<math>~y'</math>

<math>~=</math>

<math>~ y \cos\theta + (z - z_0) \sin\theta \, ,</math>

<math>~z'</math>

<math>~=</math>

<math>~ (z-z_0) \cos\theta - y \sin\theta \, .</math>

we have,

<math>~\boldsymbol{u'}</math>

<math>~~~\rightarrow~~~</math>

<math>~ \boldsymbol{\hat\imath} \dot\varphi \biggl( \frac{x_\mathrm{max}}{y_\mathrm{max}}\biggr)\biggl\{y_c - y\cos\theta - (z - z_0)\sin\theta\biggr\} + \boldsymbol{\hat\jmath} \dot\varphi \biggl( \frac{y_\mathrm{max}}{x_\mathrm{max}}\biggr) \biggr\{ x\cos\theta \biggr\} + \boldsymbol{\hat{k}} \dot\varphi \biggl( \frac{y_\mathrm{max}}{x_\mathrm{max}}\biggr) \biggr\{ x\sin\theta \biggr\} \, . </math>

Therefore, as viewed from the body frame, the individual velocity components are,

<math>~\boldsymbol{\hat\jmath} \cdot \boldsymbol{u'}</math>

<math>~=</math>

<math>~ \dot\varphi \biggl( \frac{y_\mathrm{max}}{x_\mathrm{max}}\biggr) x\cos\theta </math>

 

<math>~=</math>

<math>~x \cdot \biggl\{ \frac{b^4 c^4(\gamma^2\Omega_3^2 + \beta^2\Omega_2^2)}{a^4 (c^4 \gamma^2 \Omega_3^2 + b^4 \beta^2 \Omega_2^2)} \frac{a^4}{b^4 c^4} \biggl(\gamma^2\Omega_3^2 + \beta^2\Omega_2^2 \biggr) (c^4 \gamma^2 \Omega_3^2 + b^4 \beta^2 \Omega_2^2) \biggr\}^{1 / 4} \frac{{c^2 \gamma \Omega_3} }{[c^4 \gamma^2 \Omega_3^2 + b^4\beta^2 \Omega_2^2]^{1 / 2}} </math>

Test Orthogonality

Let's see if the Riemann fluid velocity vector is everywhere tangent to our identified off-center elliptical (Lagrangian-particle) orbital trajectories. We can determine this by first deriving an expression for the vector normal to the trajectories, then see if the dot product of this vector and the velocity vector is everywhere zero.

The unit vector that is normal to a trajectory is obtained from the (appropriately normalized) gradient of the function,

<math>~f(x', y')</math>

<math>~\equiv</math>

<math>~\biggl[\frac{x'}{x_\mathrm{max}} \biggr]^2 + \biggl[\frac{y' - y_c}{y_\mathrm{max}} \biggr]^2 - 1\, .</math>

We find, first,

<math>~\nabla f</math>

<math>~=</math>

<math>~ \boldsymbol{\hat\imath'} \biggl[\frac{2x'}{x^2_\mathrm{max}} \biggr] + \boldsymbol{\hat\jmath'} \biggl[\frac{2(y' - y_c)}{y^2_\mathrm{max}} \biggr] \, , </math>

where, drawing from Step #2, above,

<math>~\frac{y_c}{z_0}</math>

<math>~=</math>

<math>~- \frac{\sin\theta}{c^2 \kappa^2} \, ,</math>

<math>~z_0^2</math>

<math>~\le</math>

<math>~c^2 + b^2\tan^2\theta \, ,</math>

<math>~y_\mathrm{max}^2</math>

<math>~=</math>

<math>~\frac{1}{\kappa^2}\biggl[ 1 - \frac{z_0^2}{c^2 } - \frac{z_0 \sin\theta}{c^2} \biggr] \, ,</math>

<math>~x_\mathrm{max}^2</math>

<math>~=</math>

<math>~a^2 \biggl[ 1 - \frac{z_0^2}{c^2 } - \frac{z_0 \sin\theta}{c^2} \biggr] \, .</math>

Note as well that,

<math>~\biggl( \frac{x_\mathrm{max}}{y_\mathrm{max}}\biggr)^2</math>

<math>~=</math>

<math>~a^2\kappa^2 = \frac{a^2}{b^2 c^2} \biggl[ c^2 \cos^2\theta + b^2 \sin^2\theta \biggr] \, .</math>

Next, we pull from Step #3, above,

<math>~\boldsymbol{u'}_\mathrm{EFE}</math>

<math>~=</math>

<math>~ \boldsymbol{\hat\imath'} \biggl[ \biggl(\frac{a}{b}\biggr)^2 \gamma \Omega_3 ( y'\cos\theta + z_0\tan\theta ) - \biggl(\frac{a}{c}\biggr)^2 \beta \Omega_2 ( y'\sin\theta ) \biggr] + \boldsymbol{\hat\jmath'} \biggl[\beta \Omega_2 \sin\theta - \gamma \Omega_3 \cos\theta \biggr]x' \, . </math>

The dot product of these two vectors is,

<math>~\nabla f \cdot \boldsymbol{u'}_\mathrm{EFE}</math>

<math>~=</math>

<math>~ \boldsymbol{\hat\imath'} \biggl[ \biggl(\frac{a}{b}\biggr)^2 \gamma \Omega_3 ( y'\cos\theta + z_0\tan\theta ) - \biggl(\frac{a}{c}\biggr)^2 \beta \Omega_2 ( y'\sin\theta ) \biggr] + \boldsymbol{\hat\jmath'} \biggl[\beta \Omega_2 \sin\theta - \gamma \Omega_3 \cos\theta \biggr]x' \, . </math>

See Also


Whitworth's (1981) Isothermal Free-Energy Surface

© 2014 - 2021 by Joel E. Tohline
|   H_Book Home   |   YouTube   |
Appendices: | Equations | Variables | References | Ramblings | Images | myphys.lsu | ADS |
Recommended citation:   Tohline, Joel E. (2021), The Structure, Stability, & Dynamics of Self-Gravitating Fluids, a (MediaWiki-based) Vistrails.org publication, https://www.vistrails.org/index.php/User:Tohline/citation