Difference between revisions of "User:Tohline/Appendix/Ramblings/Hybrid Scheme Implications"

From VistrailsWiki
Jump to navigation Jump to search
Line 3,530: Line 3,530:
   </td>
   </td>
</tr>
</tr>
<tr><td align="center" colspan="3">[ [[User:Tohline/Appendix/References#EFE|EFE]], <font color="#00CC00">Chapter 7, &sect;48, Eq. (33)</font> ]</td></tr>
</table>
</table>




{{LSU_HBook_footer}}
{{LSU_HBook_footer}}

Revision as of 23:12, 10 September 2020

Implications of Hybrid Scheme

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

Background

Key H_Book Chapters

[Ref01]   Inertial-Frame Euler Equation

[Ref02]   Traditional Description of Rotating Reference Frame

[Ref03]   Hybrid Advection Scheme

[Ref04]   Riemann S-type Ellipsoids

[Ref05]   Korycansky and Papaloizou (1996)

Principal Governing Equations

Quoting from [Ref01] … Among the principal governing equations we have included the inertial-frame,

Lagrangian Representation
of the Euler Equation,

LSU Key.png

<math>\frac{d\vec{v}}{dt} = - \frac{1}{\rho} \nabla P - \nabla \Phi</math>

[EFE], Chap. 2, §11, p. 20, Eq. (38)
[BLRY07], p. 13, Eq. (1.55)

Shifting into a rotating frame characterized by the angular velocity vector,

<math>~\vec{\Omega}_f \equiv \hat\mathbf{k} \Omega_f \, ,</math>

and applying the operations that are specified in the first few subsections of [Ref02], we recognize the following relationships …

<math>~\vec{v}_\mathrm{inertial}</math>

<math>~=</math>

<math>~\vec{v}_\mathrm{rot} + {\vec\Omega}_f \times \vec{x} \, ,</math>

<math>~\biggl[ \frac{d \vec{v}}{dt} \biggr]_\mathrm{inertial}</math>

<math>~=</math>

<math>~ \biggl[ \frac{d \vec{v}}{dt} \biggr]_\mathrm{rot} + 2{\vec\Omega}_f \times {\vec{v}}_\mathrm{rot} + {\vec\Omega}_f \times ({\vec\Omega}_f \times \vec{x}) </math>

 

<math>~=</math>

<math>~ \biggl[ \frac{d \vec{v}}{dt} \biggr]_\mathrm{rot} + 2{\vec\Omega}_f \times {\vec{v}}_\mathrm{rot} - \frac{1}{2} \nabla | {\vec\Omega}_f \times \vec{x}|^2 </math>

 

<math>~=</math>

<math>~ \biggl[ \frac{\partial \vec{v}}{\partial t} \biggr]_\mathrm{rot} + ({\vec{v}}_\mathrm{rot} \cdot \nabla){\vec{v}}_\mathrm{rot} + 2{\vec\Omega}_f \times {\vec{v}}_\mathrm{rot} - \frac{1}{2} \nabla | {\vec\Omega}_f \times \vec{x}|^2 \, .</math>

Making this substitution on the left-hand-side of the above-specified "Lagrangian Representation of the Euler Equation," we obtain what we have referred to also in [Ref02] as the,

Eulerian Representation
of the Euler Equation
as viewed from a Rotating Reference Frame

<math>\biggl[\frac{\partial\vec{v}}{\partial t}\biggr]_\mathrm{rot} + ({\vec{v}}_\mathrm{rot}\cdot \nabla) {\vec{v}}_\mathrm{rot}= - \frac{1}{\rho} \nabla P - \nabla \biggl[\Phi - \frac{1}{2}|{\vec{\Omega}}_f \times \vec{x}|^2 \biggr] - 2{\vec{\Omega}}_f \times {\vec{v}}_\mathrm{rot} \, .</math>

This form of the Euler equation also appears early in [Ref05], where we set up a discussion of the paper by Korycansky & Papaloizou (1996, ApJS, 105, 181; hereafter KP96). But, for now, let's back up a couple of steps and retain the total time derivative on the left-hand-side. That is, let's select as the foundation expression the,

Lagrangian Representation
of the Euler Equation
as viewed from a Rotating Reference Frame

<math>~\biggl[ \frac{d \vec{v}}{dt} \biggr]_\mathrm{rot} </math>

<math>~=</math>

<math>~- \frac{1}{\rho} \nabla P - \nabla \Phi - 2{\vec\Omega}_f \times {\vec{v}}_\mathrm{rot} - {\vec\Omega}_f \times ({\vec\Omega}_f \times \vec{x}) \, ,</math>

[EFE], Chap. 2, §12, p. 25, Eq. (62)

which also serves as the foundation of most of our [Ref03] discussions.

Exercising the Hybrid Scheme

Focus on Tracking Angular Momentum

Let's begin by using <math>~\mathbf{u'}</math>, instead of <math>~{\vec{v}}_\mathrm{rot}</math>, to represent the fluid velocity vector as viewed from the rotating frame of reference. Our foundation expression becomes,

<math>~\frac{d \bold{u'}}{dt} </math>

<math>~=</math>

<math>~- \frac{1}{\rho} \nabla P - \nabla \Phi - 2{\vec\Omega}_f \times \bold{u}' - {\vec\Omega}_f \times ({\vec\Omega}_f \times \vec{x}) \, ,</math>

where we appreciate that we can move from the Lagrangian to an Eulerian representation by employing the operator substitution,

<math>~\frac{d}{dt}</math>

<math>~\rightarrow</math>

<math>~\frac{\partial}{\partial t} + \mathbf{u'} \cdot \nabla </math>

Next, using [Ref03] as a guide, let's focus on tracking angular momentum. We need to break the vector momentum equation, as well as the velocity vectors, into their <math>~(\bold{\hat{e}}_\varpi, \bold{\hat{e}}_\varphi, \bold{\hat{k}})</math> components.

NOTE: For the time being, we will write the velocity vector in terms of generic components, namely,

<math>~\bold{u}' = \bold{\hat{e}}_\varpi u'_\varpi + \bold{\hat{e}}_\varphi u'_\varphi + \bold{\hat{k}}u'_z \, .</math>

But, eventually, we want to explicitly insert the rotating-frame velocity that underpins the equilibrium properties of Riemann S-type ellipsoids. In Chap. 7, §47, Eq. 1 (p. 130) of [EFE], this is given in Cartesian coordinates, so we will need to convert his expressions to the equivalent cylindrical-coordinate components.

The time-derivative on the left-hand-side of our foundation expression becomes,

<math> \frac{d\mathbf{u'}}{dt} </math>

<math>~=~</math>

<math> \frac{d}{dt} [ \mathbf{\hat{e}}_\varpi u'_\varpi + \mathbf{\hat{e}}_\varphi u'_\varphi + \mathbf{\hat{k}} u'_z ] </math>

 

<math>~=~</math>

<math> \mathbf{\hat{e}}_\varpi \frac{d u'_\varpi}{dt} + \mathbf{\hat{e}}_\varphi \frac{d u'_\varphi}{dt} + \mathbf{\hat{k}} \frac{d u'_z}{dt} + ( u'_\varpi) \frac{d}{dt}\mathbf{\hat{e}}_\varpi + ( u'_\varphi) \frac{d}{dt}\mathbf{\hat{e}}_\varphi </math>

 

<math>~=~</math>

<math> \mathbf{\hat{e}}_\varpi \frac{d u'_\varpi}{dt} + \mathbf{\hat{e}}_\varphi \frac{d u'_\varphi}{dt} + \mathbf{\hat{k}} \frac{d u'_z}{dt} + \mathbf{\hat{e}}_\varphi(u'_\varpi) \frac{u'_\varphi}{\varpi} - \mathbf{\hat{e}}_\varpi(u'_\varphi) \frac{u'_\varphi}{\varpi} \, . </math>

We also recognize that, when expressed in cylindrical coordinates,

<math> ~{\vec{\Omega}}_f \times \vec{x} </math>

<math>~=~</math>

<math> {\hat\mathbf{k}} \Omega_f\times (\mathbf{\hat{e}}_\varpi \varpi + \mathbf{\hat{k}}z) = \mathbf{\hat{e}}_\varphi \Omega_f \varpi \, , </math>

<math> {\vec{\Omega}}_f \times ({\vec{\Omega}}_f \times \vec{x}) </math>

<math>~=~</math>

<math> \hat{\mathbf{k}} \Omega_f \times ( \mathbf{\hat{e}}_\varphi \Omega_f \varpi ) = - \mathbf{\hat{e}}_\varpi \Omega_f^2 \varpi \, , </math>

<math> {\vec{\Omega}}_f \times {\mathbf{u'}} </math>

<math>~=~</math>

<math> {\hat\mathbf{k}} \Omega_f\times (\mathbf{\hat{e}}_\varpi u'_\varpi + \mathbf{\hat{e}}_\varphi u'_\varphi + \mathbf{\hat{k}}u'_z) = \mathbf{\hat{e}}_\varphi \Omega_f u'_\varpi - \mathbf{\hat{e}}_\varpi \Omega_f u'_\varphi \, , </math>

<math> {\vec{v}}_\mathrm{inertial} </math>

<math>~=~</math>

<math> \mathbf{u'} + \mathbf{\hat{e}}_\varphi \Omega_f \varpi \, . </math>

The set of scalar momentum-component equations is obtained by "dotting" each unit vector into the vector equation.

<math>\mathbf{\hat{e}}_\varpi:</math>

<math>~\frac{d u'_\varpi}{dt} - \frac{(u'_\varphi)^2}{\varpi} </math>

<math>~=</math>

<math>~- \mathbf{\hat{e}}_\varpi \cdot \frac{\nabla P}{\rho} - \mathbf{\hat{e}}_\varpi \cdot \nabla \Phi + 2 \biggl[ \Omega_f u'_\varphi \biggr] + \Omega_f^2 \varpi </math>

<math>~\Rightarrow ~~~ \frac{d u'_\varpi}{dt} </math>

<math>~=</math>

<math>~- \mathbf{\hat{e}}_\varpi \cdot \frac{\nabla P}{\rho} - \mathbf{\hat{e}}_\varpi \cdot \nabla \Phi + \frac{1}{\varpi} \biggl[ (u'_\varphi)^2 + 2 \Omega_f u'_\varphi \varpi + \Omega_f^2 \varpi^2 \biggr]</math>

 

<math>~=</math>

<math>~ - \mathbf{\hat{e}}_\varpi \cdot \frac{\nabla P}{\rho} - \mathbf{\hat{e}}_\varpi \cdot \nabla \Phi + \frac{1}{\varpi} (u'_\varphi + \Omega_f \varpi)^2 \, ; </math>

<math>\mathbf{\hat{e}}_\varphi:</math>

<math>~\frac{d u'_\varphi}{dt} + \frac{u'_\varpi u'_\varphi}{\varpi} </math>

<math>~=</math>

<math>~- \mathbf{\hat{e}}_\varphi \cdot \frac{\nabla P}{\rho} - \mathbf{\hat{e}}_\varphi \cdot \nabla \Phi - 2\biggl[ \Omega_f u'_\varpi \biggr] </math>

(mult. thru by ϖ)   <math>~\Rightarrow ~~~\frac{d (\varpi u'_\varphi )}{dt} </math>

<math>~=</math>

<math>~- \mathbf{\hat{e}}_\varphi \cdot \frac{\varpi \nabla P}{\rho} - \mathbf{\hat{e}}_\varphi \cdot \varpi \nabla \Phi - 2 \Omega_f \varpi u'_\varpi \, ; </math>

<math>\mathbf{\hat{k}}:</math>

<math>~\frac{d u'_z}{dt} </math>

<math>~=</math>

<math>~- \mathbf{\hat{k}} \cdot \frac{\nabla P }{\rho} - \mathbf{\hat{k}} \cdot \nabla \Phi \, . </math>

Now, recalling that <math>~\mathbf{u'} = (\mathbf{v} - \mathbf{\hat{e}}_\varphi \varpi \Omega_f)</math>, let's make the substitutions …

<math>~u'_\varpi \rightarrow v_\varpi \, ,</math>     

<math>~u'_\varphi \rightarrow (v_\varphi - \varpi\Omega_f) \, ,</math>      and,      

<math>~u'_z \rightarrow v_z \, .</math>

This mapping gives,

<math>\mathbf{\hat{e}}_\varphi:</math>

<math>~\frac{d [\varpi v_\varphi - \varpi^2 \Omega_f]}{dt} </math>

<math>~=</math>

<math>~- \mathbf{\hat{e}}_\varphi \cdot \frac{\varpi \nabla P}{\rho} - \mathbf{\hat{e}}_\varphi \cdot \varpi \nabla \Phi - 2 \Omega_f \varpi v_\varpi \, ; </math>

<math>~\Rightarrow ~~~ \frac{d (\varpi v_\varphi )}{dt} </math>

<math>~=</math>

<math>~- \mathbf{\hat{e}}_\varphi \cdot \frac{\varpi \nabla P}{\rho} - \mathbf{\hat{e}}_\varphi \cdot \varpi \nabla \Phi \, ; </math>

<math>~\Rightarrow ~~~ \frac{1}{\varpi} ~\frac{d (\varpi v_\varphi )}{dt} </math>

<math>~=</math>

<math>~- \mathbf{\hat{e}}_\varphi \cdot \biggl[ \frac{\nabla P}{\rho} + \nabla \Phi \biggr] \, ; </math>

<math>\mathbf{\hat{k}}:</math>

<math>~\frac{d v_z}{dt} </math>

<math>~=</math>

<math>~- \mathbf{\hat{k}} \cdot \biggl[ \frac{\nabla P }{\rho} + \nabla \Phi \biggr] \, . </math>

<math>\mathbf{\hat{e}}_\varpi:</math>

<math>~\frac{d v_\varpi}{dt} </math>

<math>~=</math>

<math>~ - \mathbf{\hat{e}}_\varpi \cdot \biggl[ \frac{\nabla P}{\rho} + \nabla \Phi \biggr] + \frac{v_\varphi^2}{\varpi} \, ; </math>

Steady-State Velocity Field for Jacobi Ellipsoids

In steady-state, the (Lagrangian time-derivative) operator on the left-hand-side of all three component equations maps to the following operator:

<math>~\mathbf{u'} \cdot \nabla</math>

<math>~=</math>

<math>~\sum_{i=1}^3 u'_i \frac{\partial}{\partial x_i} \, ,</math>

        (in Cartesian coordinates);

<math>~\mathbf{u'} \cdot \nabla</math>

<math>~=</math>

<math>~ u'_\varpi \frac{\partial}{\partial \varpi} + \frac{u'_\varphi}{\varpi} \frac{\partial}{\partial \varphi} + u'_z \frac{\partial}{\partial z} \, ,</math>

        (in cylindrical coordinates);

We know, as well, that,

<math>~u'_\varpi = u'_x \cos\varphi + u'_y \sin\varphi \, ,</math>

      and,      

<math>~u'_\varphi = u'_y \cos\varphi - u'_x \sin\varphi \, .</math>

Hence, the cylindrical-coordinate-based operator may be rewritten as,

<math>~\mathbf{u'} \cdot \nabla</math>

<math>~=</math>

<math>~ ( u'_x \cos\varphi + u'_y \sin\varphi ) \frac{\partial}{\partial \varpi} + ( u'_y \cos\varphi - u'_x \sin\varphi )\frac{1}{\varpi} \frac{\partial}{\partial \varphi} + u'_z \frac{\partial}{\partial z} \, .</math>

Drawing from [ Ref04 ] … As Ou(2006) has pointed out, the velocity field of a Riemann S-type ellipsoid as viewed from a frame rotating with angular velocity <math>~{\vec{\Omega}}_f = \boldsymbol{\hat{k}} \Omega_f</math> takes the following form:

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

<math>~=</math>

<math>~\lambda \biggl[ \boldsymbol{\hat{\imath}} \biggl(\frac{a}{b}\biggr)y - \boldsymbol{\hat{\jmath}} \biggl(\frac{b}{a}\biggr)x \biggr] \, ,</math>

Ou(2006), p. 550, §2, Eq. (3)

where <math>~\lambda</math> is a constant that determines the magnitude of the internal motion of the fluid, and the origin of the x-y coordinate system is at the center of the ellipsoid. This velocity field, <math>~\mathbf{u'}</math>, is designed so that velocity vectors everywhere are always aligned with elliptical stream lines by demanding that they be tangent to the equi-effective-potential contours, which are concentric ellipses. Hence, for Riemann S-type ellipsoids, we have,

<math>~u'_x = \lambda\biggl(\frac{a}{b}\biggr)y = \lambda\biggl(\frac{a}{b}\biggr)\varpi \sin\varphi \, ;</math>

       

<math>~u'_y = -\lambda\biggl(\frac{b}{a}\biggr)x = -\lambda\biggl(\frac{b}{a}\biggr)\varpi \cos\varphi \, ;</math>

       

<math>~u'_z = 0 \, .</math>

So, for the velocity flow that underpins Riemann S-type ellipsoids, the cylindrical-coordinate-based operator is

<math>~\mathbf{u'} \cdot \nabla</math>

<math>~=</math>

<math>~ \biggl[ \lambda\biggl(\frac{a}{b}\biggr)\varpi \sin\varphi \cos\varphi -\lambda\biggl(\frac{b}{a}\biggr)\varpi \cos\varphi \sin\varphi \biggr] \frac{\partial}{\partial \varpi} + \biggl[ -\lambda\biggl(\frac{b}{a}\biggr)\varpi \cos\varphi \cos\varphi - \lambda\biggl(\frac{a}{b}\biggr)\varpi \sin\varphi \sin\varphi \biggr] \frac{1}{\varpi} \frac{\partial}{\partial \varphi} </math>

 

<math>~=</math>

<math>~ \biggl[ \biggl(\frac{a}{b}\biggr) - \biggl(\frac{b}{a}\biggr) \biggr]\lambda \varpi \sin\varphi \cos\varphi \frac{\partial}{\partial \varpi} - \biggl[ \biggl(\frac{b}{a}\biggr) \cos^2\varphi + \biggl(\frac{a}{b}\biggr) \sin^2\varphi \biggr]\lambda \frac{\partial}{\partial \varphi} \, .</math>

And, given that,

<math>~\mathbf{\hat{e}}_\varphi \Omega_f \varpi</math>

<math>~=</math>

<math>~ \Omega_f \varpi \biggl[ \boldsymbol{\hat{\jmath}} \cos\varphi - \boldsymbol{\hat{\imath}} \sin\varphi \biggr] \, , </math>

the inertial-frame velocity components are,

<math>~ v_x = \lambda\biggl(\frac{a}{b}\biggr)\varpi \sin\varphi - \Omega_f \varpi \sin\varphi = \biggl[ \lambda\biggl(\frac{a}{b}\biggr) - \Omega_f \biggr]\varpi\sin\varphi \, ;</math>

       

<math>~ v_y = -\lambda\biggl(\frac{b}{a}\biggr) \varpi\cos\varphi + \Omega_f\varpi \cos\varphi = \biggl[\Omega_f -\lambda\biggl(\frac{b}{a}\biggr) \biggr]\varpi\cos\varphi \, ;</math>

       

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

That is,

<math>~v_\varpi = v_x \cos\varphi + v_y \sin\varphi</math>

<math>~=</math>

<math>~ \biggl[ \lambda\biggl(\frac{a}{b}\biggr) - \Omega_f \biggr]\varpi\sin\varphi \cos\varphi + \biggl[\Omega_f -\lambda\biggl(\frac{b}{a}\biggr) \biggr]\varpi\cos\varphi \sin\varphi </math>

 

<math>~=</math>

<math>~ \biggl[ \biggl(\frac{a}{b}\biggr) - \biggl(\frac{b}{a}\biggr) \biggr]\lambda \varpi\sin\varphi \cos\varphi \, ; </math>

<math>~v_\varphi = v_y \cos\varphi - v_x \sin\varphi</math>

<math>~=</math>

<math>~ \biggl[\Omega_f -\lambda\biggl(\frac{b}{a}\biggr) \biggr]\varpi\cos^2\varphi - \biggl[ \lambda\biggl(\frac{a}{b}\biggr) - \Omega_f \biggr]\varpi\sin^2\varphi </math>

 

<math>~=</math>

<math>~ \Omega_f \varpi -\lambda \varpi \biggl[ \biggl(\frac{b}{a}\biggr) \cos^2\varphi +\biggl(\frac{a}{b}\biggr) \sin^2\varphi \biggr] \, . </math>

Note, as well, that,

<math>~\frac{v_\varphi^2}{\varpi}</math>

<math>~=</math>

<math>~ \frac{1}{\varpi} \biggl\{ \Omega_f \varpi -\lambda \varpi \biggl[ \biggl(\frac{b}{a}\biggr) \cos^2\varphi +\biggl(\frac{a}{b}\biggr) \sin^2\varphi \biggr] \biggr\}^2 </math>

 

<math>~=</math>

<math>~ \varpi \biggl\{ \Omega_f^2 - 2\lambda \Omega_f \biggl[ \biggl(\frac{b}{a}\biggr) \cos^2\varphi +\biggl(\frac{a}{b}\biggr) \sin^2\varphi \biggr] + \lambda^2 \biggl[ \biggl(\frac{b}{a}\biggr) \cos^2\varphi +\biggl(\frac{a}{b}\biggr) \sin^2\varphi \biggr]^2 \biggr\} \, . </math>


Finally, then, we find that the left-hand-side of the momentum-component expressions are,

<math>\mathbf{\hat{k}}:</math>

<math>~\frac{d v_z}{dt} </math>

<math>~=</math>

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

<math>\mathbf{\hat{e}}_\varpi:</math>

<math>~\frac{d v_\varpi}{dt} </math>

<math>~=</math>

<math>~ \biggl\{ \biggl[ \biggl(\frac{a}{b}\biggr) - \biggl(\frac{b}{a}\biggr) \biggr]\lambda \varpi \sin\varphi \cos\varphi \frac{\partial}{\partial \varpi} - \biggl[ \biggl(\frac{b}{a}\biggr) \cos^2\varphi + \biggl(\frac{a}{b}\biggr) \sin^2\varphi \biggr]\lambda \frac{\partial}{\partial \varphi} \biggr\}\biggl[ \biggl(\frac{a}{b}\biggr) - \biggl(\frac{b}{a}\biggr) \biggr]\lambda \varpi\sin\varphi \cos\varphi </math>

 

 

<math>~=</math>

<math>~ \lambda \biggl[ \biggl(\frac{a}{b}\biggr) - \biggl(\frac{b}{a}\biggr) \biggr]\biggl\{ \biggl[ \biggl(\frac{a}{b}\biggr) - \biggl(\frac{b}{a}\biggr) \biggr]\lambda \varpi \sin\varphi \cos\varphi \frac{\partial}{\partial \varpi} \biggl[ \varpi\sin\varphi \cos\varphi \biggr] - \biggl[ \biggl(\frac{b}{a}\biggr) \cos^2\varphi + \biggl(\frac{a}{b}\biggr) \sin^2\varphi \biggr]\lambda \frac{\partial}{\partial \varphi}\biggl[ \varpi\sin\varphi \cos\varphi \biggr] \biggr\} </math>

 

 

<math>~=</math>

<math>~ \lambda \biggl[ \biggl(\frac{a}{b}\biggr) - \biggl(\frac{b}{a}\biggr) \biggr]\biggl\{ \biggl[ \biggl(\frac{a}{b}\biggr) - \biggl(\frac{b}{a}\biggr) \biggr]\lambda \varpi \sin^2\varphi \cos^2\varphi + \biggl[ \biggl(\frac{b}{a}\biggr) \cos^2\varphi + \biggl(\frac{a}{b}\biggr) \sin^2\varphi \biggr] \lambda\varpi \biggl[ \sin^2\varphi - \cos^2\varphi \biggr] \biggr\} </math>

 

 

<math>~=</math>

<math>~ \lambda^2 \varpi \biggl( \frac{a}{b} - \frac{b}{a} \biggr) \biggl[ \biggl( \frac{a}{b}\biggr) \sin^4\varphi - \biggl(\frac{b}{a}\biggr) \cos^4\varphi \biggr] \, ; </math>

<math>\mathbf{\hat{e}}_\varphi:</math>

<math>~\frac{1}{\varpi} ~\frac{d (\varpi v_\varphi )}{dt} </math>

<math>~=</math>

<math>~ \frac{1}{\varpi}\biggl\{ \biggl[ \biggl(\frac{a}{b}\biggr) - \biggl(\frac{b}{a}\biggr) \biggr]\lambda \varpi \sin\varphi \cos\varphi \frac{\partial}{\partial \varpi} - \biggl[ \biggl(\frac{b}{a}\biggr) \cos^2\varphi + \biggl(\frac{a}{b}\biggr) \sin^2\varphi \biggr]\lambda \frac{\partial}{\partial \varphi} \biggr\} \biggl\{ \biggl[\Omega_f -\lambda\biggl(\frac{b}{a}\biggr) \biggr]\varpi^2\cos^2\varphi - \biggl[ \lambda\biggl(\frac{a}{b}\biggr) - \Omega_f \biggr]\varpi^2\sin^2\varphi \biggr\} </math>

 

<math>~=</math>

<math>~ \biggl[ \biggl(\frac{a}{b}\biggr) - \biggl(\frac{b}{a}\biggr) \biggr] 2\varpi \lambda \sin\varphi \cos\varphi \biggl\{ \biggl[\Omega_f -\lambda\biggl(\frac{b}{a}\biggr) \biggr]\cos^2\varphi - \biggl[ \lambda\biggl(\frac{a}{b}\biggr) - \Omega_f \biggr]\sin^2\varphi \biggr\} </math>

 

 

<math>~+ \biggl[ \biggl(\frac{b}{a}\biggr) \cos^2\varphi + \biggl(\frac{a}{b}\biggr) \sin^2\varphi \biggr] \biggl[ \frac{a}{b} - \frac{b}{a} \biggr] 2\varpi \lambda^2 \sin\varphi \cos\varphi </math>

 

<math>~=</math>

<math>~ \biggl( \frac{a}{b} - \frac{b}{a} \biggr) 2\varpi \lambda \sin\varphi \cos\varphi \biggl\{ \biggl[\Omega_f -\lambda\biggl(\frac{b}{a}\biggr) \biggr]\cos^2\varphi - \biggl[ \lambda\biggl(\frac{a}{b}\biggr) - \Omega_f \biggr]\sin^2\varphi +2\lambda\biggl[ \biggl(\frac{b}{a}\biggr) \cos^2\varphi + \biggl(\frac{a}{b}\biggr) \sin^2\varphi \biggr] \biggr\} </math>

 

<math>~=</math>

<math>~ \biggl( \frac{a}{b} - \frac{b}{a} \biggr) 2\varpi \lambda \sin\varphi \cos\varphi \biggl\{ \biggl[\Omega_f +\lambda\biggl(\frac{b}{a}\biggr)\biggr]\cos^2\varphi + \biggl[ \Omega_f + \lambda\biggl( \frac{a}{b}\biggr) \biggr]\sin^2\varphi \biggr\} \, . </math>

Try Again

An example equation of motion is,

<math>~\bold{a} = - \frac{1}{\rho}\bold\nabla P - \bold\nabla \Phi_\mathrm{grav} \, .</math>

Here we will focus only on the left-hand-side of this equation, examining various ways the vector acceleration may be mathematically expressed. We will consider, in particular, building a model in a curvilinear (cylindrical), rather than a Cartesian, coordinate base; and viewing the model's evolution in a rotating, rather than inertial, frame of reference.

Inertial Frame

As viewed from a cylindrical-coordinate-based <math>~(\varpi, \varphi, z)</math> inertial reference frame, we are interested in specifying the location,

<math>~\bold{x} = \mathbf{\hat{e}}_\varpi \varpi + \bold{\hat{k}} z \, ,</math>

[BT87], p. 646, Appendix §1.B.2, Eq. (1B-18)

of a Lagrangian fluid element at time <math>~t = 0</math> — hereafter denoted by the subscript, <math>~0</math> — as well as at later times. Although the position vector, <math>~\bold{x}</math>, does not explicitly display a dependence on the azimuthal coordinate angle, <math>~\varphi</math>, it is important to realize that the orientation in space of the unit vector, <math>~\bold{\hat{e}}_\varpi</math>, does depend on the value of this coordinate angle.

At any point in time, the instantaneous velocity of this Lagrangian fluid element will correspond precisely with the (total) time-derivative of its instantaneous position vector, that is,

<math>~\bold{v}</math>

<math>~\equiv</math>

<math>~\frac{d\bold{x}}{dt}</math>

<math>~=</math>

<math>~\bold{\hat{e}}_\varpi \frac{d\varpi}{dt} + \bold{\hat{k}} \frac{dz}{dt} + \varpi \frac{d \bold{\hat{e}}_\varpi}{dt}</math>

 

 

 

<math>~=</math>

<math>~\bold{\hat{e}}_\varpi \frac{d\varpi}{dt} + \bold{\hat{k}} \frac{dz}{dt} + \varpi \biggl[ \bold{\hat{e}}_\varphi \frac{d\varphi}{dt} \bigg] \, .</math>

[BT87], p. 647, Appendix §1.B.2, Eq. (1B-23)

In carrying out this time differentiation, the last term on the right-hand-side accounts for the aforementioned dependence of <math>~\bold{\hat{e}}_\varpi</math> on <math>~\varphi</math>. Similarly, the following component breakdown of the Lagrangian fluid element's acceleration takes into account the dependence of <math>~\bold{\hat{e}}_\varphi</math> on <math>~\varphi</math>:

<math>~\bold{a}</math>

<math>~\equiv</math>

<math>~\frac{d\bold{v}}{dt}</math>

<math>~=</math>

<math>~ \bold{\hat{e}}_\varpi \frac{d^2\varpi}{dt^2} + \bold{\hat{k}} \frac{d^2z}{dt^2} + \bold{\hat{e}}_\varphi \biggl[\frac{d\varpi}{dt} \cdot \frac{d\varphi}{dt} + \varpi \frac{d^2\varphi}{dt^2}\biggr] + \varpi \frac{d\varphi}{dt} \biggl[ \frac{d\bold{\hat{e}}_\varphi}{dt} \biggr] + \frac{d\varpi}{dt} \biggl[ \frac{d\bold{\hat{e}}_\varpi}{dt} \biggr] </math>

 

<math>~=</math>

<math>~ \bold{\hat{e}}_\varpi \frac{d^2\varpi}{dt^2} + \bold{\hat{k}} \frac{d^2z}{dt^2} + \bold{\hat{e}}_\varphi \biggl[\frac{d\varpi}{dt} \cdot \frac{d\varphi}{dt} + \varpi \frac{d^2\varphi}{dt^2}\biggr] + \varpi \frac{d\varphi}{dt} \biggl[- \bold{\hat{e}}_\varpi \frac{d\varphi}{dt} \biggr] + \frac{d\varpi}{dt} \biggl[ \bold{\hat{e}}_\varphi \frac{d\varphi}{dt} \biggr] </math>

 

<math>~=</math>

<math>~ \bold{\hat{e}}_\varpi \biggl[\frac{d^2\varpi}{dt^2} - \varpi \biggl(\frac{d\varphi}{dt}\biggr)^2 \biggr] + \bold{\hat{e}}_\varphi \biggl[ 2 \biggl( \frac{d\varpi}{dt} \cdot \frac{d\varphi}{dt} \biggr) + \varpi \frac{d^2\varphi}{dt^2}\biggr] + \bold{\hat{k}} \frac{d^2z}{dt^2} \, . </math>

[BT87], p. 647, Appendix §1.B.2, Eq. (1B-24)

Let's rewrite the velocity vector as,

<math>~\bold{v}</math>

<math>~=</math>

<math>~\bold{\hat{e}}_\varpi \dot\varpi + \bold{\hat{e}}_\varphi \varpi \dot\varphi + \bold{\hat{k}} \dot{z} \, ,</math>

and (the second line of) this acceleration expression as,

<math>~~\bold{a} \equiv \frac{d\bold{v}}{dt} = \bold{\hat{e}}_\varpi \frac{d \dot\varpi}{dt} + \bold{\hat{e}}_\varphi \frac{d}{dt}\biggl[\varpi \dot\varphi \biggr] + \bold{\hat{k}} \frac{d \dot{z}}{dt} + \underbrace{ \dot\varpi \biggl[ \bold{\hat{e}}_\varphi \frac{d\varphi}{dt} \biggr] - \varpi \dot\varphi \biggl[\bold{\hat{e}}_\varpi \frac{d\varphi}{dt} \biggr] }_\text{curvature terms}\, . </math>

Now, if <math>~\bold{B}</math> is a vector quantity that characterizes some property of a fluid element — such as momentum density, velocity, or vorticity — the difference between the Lagrangian and Eulerian time-derivatives of that vector quantity is given by the expression,

<math>~\frac{d\bold{B}}{dt} - \frac{\partial \bold{B}}{\partial t}</math>

<math>~=</math>

<math>~(\bold{v} \cdot \bold\nabla)\bold{B} \, ,</math>

where the various elements of this right-hand-side mathematical operator can be obtained by replacing <math>~\bold{A}</math> with <math>~\bold{v}</math> in the so-called convective operator.

Convective Operator in Cylindrical Coordinates

<math>~(\bold{A} \cdot \bold\nabla) \bold{B}</math>

<math>~=</math>

<math>~ \bold{\hat{e}}_\varpi \biggl[ A_\varpi \frac{\partial B_\varpi}{\partial \varpi} + \frac{A_\varphi }{\varpi}\frac{\partial B_\varpi}{\partial \varphi} + A_z \frac{\partial B_\varpi}{\partial z} - \frac{A_\varphi B_\varphi}{\varpi} \biggr] </math>

 

 

<math>~ + \bold{\hat{e}}_\varphi \biggl[ A_\varpi \frac{\partial B_\varphi}{\partial \varpi} + \frac{A_\varphi }{\varpi}\frac{\partial B_\varphi}{\partial \varphi} + A_z \frac{\partial B_\varphi}{\partial z} + \frac{A_\varphi B_\varpi}{\varpi} \biggr] </math>

 

 

<math>~ + \bold{\hat{e}}_z \biggl[ A_\varpi \frac{\partial B_z}{\partial \varpi} + \frac{A_\varphi }{\varpi}\frac{\partial B_z}{\partial \varphi} + A_z \frac{\partial B_z}{\partial z} \biggr] \, . </math>

[BT87], p. 651, Appendix §1.B.3, Eq. (1B-54)

We will adopt the following, more compact notation:

<math>~(\bold{A} \cdot \bold\nabla) \bold{B}</math>

<math>~=</math>

<math>~ \bold{\hat{e}}_\varpi \biggl[ \mathcal{L}_A B_\varpi- \frac{A_\varphi B_\varphi}{\varpi} \biggr] + \bold{\hat{e}}_\varphi \biggl[ \mathcal{L}_A B_\varphi + \frac{A_\varphi B_\varpi}{\varpi} \biggr] + \bold{\hat{e}}_z \biggl[ \mathcal{L}_A B_z \biggr] \, , </math>

where the operator,

<math>~\mathcal{L}_A</math>

<math>~\equiv</math>

<math>~ \biggl[ A_\varpi \frac{\partial }{\partial \varpi} + \frac{A_\varphi }{\varpi}\frac{\partial }{\partial \varphi} + A_z \frac{\partial }{\partial z} \biggr] \, . </math>

In particular, if we are examining the behavior of the fluid velocity <math>~(\bold{B} \rightarrow \bold{v} )</math>, we find that,

<math>~\frac{d\bold{v}}{dt} - \frac{\partial \bold{v}}{\partial t}</math>

<math>~=~(\bold{v} \cdot \bold\nabla)\bold{v} </math>

 

<math>~=~ \bold{\hat{e}}_\varpi \biggl[ \mathcal{L}_v \dot\varpi \biggr] + \bold{\hat{e}}_\varphi \biggl[ \mathcal{L}_v (\varpi \dot\varphi) \biggr] + \bold{\hat{e}}_z \biggl[ \mathcal{L}_v \dot{z} \biggr] ~+ \underbrace{\bold{\hat{e}}_\varphi (\dot\varphi \dot\varpi ) -~ \bold{\hat{e}}_\varpi ( \varpi {\dot\varphi}^2 ) }_\text{curvature terms} \, , </math>

where the operator,

<math>~\mathcal{L}_v</math>

<math>~\equiv</math>

<math>~ \biggl[ v_\varpi \frac{\partial }{\partial \varpi} + \frac{v_\varphi }{\varpi}\frac{\partial }{\partial \varphi} + v_z \frac{\partial }{\partial z} \biggr] \, . </math>

Notice that the pair of "curvature terms" that appear in this expression are identical to the pair of curvature terms that appear in the acceleration expression, above. We conclude, therefore, that for each of the three separate (cylindrical-coordinate-based) components of the vector acceleration, the relationship between the Lagrangian (total) and Eulerian (partial) time derivative is, respectively,

<math>~\bold{\hat{e}}_\varpi</math>:     

<math>~\frac{d\dot\varpi}{dt} - \varpi {\dot\varphi}^2</math>

<math>~=</math>

<math>~ \frac{\partial \dot\varpi}{\partial t} + \biggl[\mathcal{L}_v \dot\varpi \biggr] - \varpi {\dot\varphi}^2 \, ; </math>

<math>~\bold{\hat{e}}_\varphi</math>:     

<math>~\frac{d (\varpi \dot\varphi ) }{dt} + \dot\varpi \dot\varphi</math>

<math>~=</math>

<math>~ \frac{\partial (\varpi \dot\varphi ) }{\partial t} + \biggl[ \mathcal{L}_v (\varpi \dot\varphi ) \biggr] + \dot\varpi \dot\varphi \, ; </math>

<math>~\bold{\hat{k}}</math>:     

<math>~\frac{d \dot{z} }{dt}</math>

<math>~=</math>

<math>~ \frac{\partial \dot{z} }{\partial t} + \biggl[ \mathcal{L}_v \dot{z} \biggr] \, . </math>

Rotating Frame

Drawing from an accompanying discussion of rotating reference frames, let's build our model in a cylindrical coordinate system that is spinning about its <math>~\bold{\hat{k}}</math>-axis with a time-independent angular velocity, <math>~\bold\Omega_f = \bold{\hat{k}} \Omega_f</math>. Furthermore, let's use <math>~\bold{u}</math> — instead of <math>~\bold{v}</math> — to represent the velocity as viewed in the rotating frame. We know that,

<math>~\bold{v}</math>

<math>~=</math>

<math>~ \bold{u} + \bold\Omega_f \times \bold{x}_\mathrm{rot} </math>

<math>~=</math>

<math>~ \bold{u} + \bold{\hat{e}}_\varphi \varpi \Omega_f \, , </math>

and,

<math>~\bold{a} = \frac{d\bold{v}}{dt} = \frac{d\bold{u}}{dt} - \overbrace{ \biggl[ \underbrace{(- 2\bold\Omega_f \times \bold{u}) }_\text{Coriolis} ~+~ \underbrace{(- \bold\Omega_f \times (\bold\Omega_f \times \bold{x}_\mathrm{rot} )) }_\text{Centrifugal} \biggr] }^\text{Fictitious accelerations} \, . </math>

[BT87], p. 664, Appendix §1.D.3, Eq. (1D-43)

Note that, in the particular case being considered here,

<math>~\bold{a}_\mathrm{fict}</math>

<math>~=</math>

<math>~ - 2\bold\Omega_f \times \bold{u} - \bold\Omega_f \times (\bold\Omega_f \times \bold{x}_\mathrm{rot} ) </math>

 

<math>~=</math>

<math>~ - 2\bold\Omega_f \times [ \bold{v} - \bold\Omega_f \times \bold{x}_\mathrm{rot} ] - \bold\Omega_f \times (\bold\Omega_f \times \bold{x}_\mathrm{rot} ) </math>

 

<math>~=</math>

<math>~ - 2\bold\Omega_f \times \bold{v} + \bold\Omega_f \times (\bold\Omega_f \times \bold{x}_\mathrm{rot} ) </math>

 

<math>~=</math>

<math>~ - 2\bold\Omega_f \times \bold{v} + \bold{\hat{k}}\Omega_f \times (\bold{\hat{e}}_\varphi \varpi \Omega_f ) </math>

 

<math>~=</math>

<math>~ - 2\Omega_f [\bold{\hat{e}}_\varphi v_\varpi - \bold{\hat{e}}_\varpi v_\varphi] - (\bold{\hat{e}}_\varpi \varpi \Omega_f^2 ) \, . </math>

 

<math>~=</math>

<math>~ + \bold{\hat{e}}_\varpi [2\Omega_f v_\varphi - \varpi \Omega_f^2 ] - \bold{\hat{e}}_\varphi 2\Omega_f v_\varpi \, . </math>

We may therefore also write,

<math>~\bold{a} + \bold{a}_\mathrm{fict} = \frac{d\bold{u}}{dt}</math>

<math>~= \frac{\partial \bold{u}}{\partial t} + (\bold{u} \cdot \bold\nabla)\bold{u} </math>

 

<math>~=~ \frac{\partial \bold{u}}{\partial t} + \bold{\hat{e}}_\varpi \biggl[ \mathcal{L}_u u_\varpi \biggr] + \bold{\hat{e}}_\varphi \biggl[ \mathcal{L}_u u_\varphi \biggr] + \bold{\hat{e}}_z \biggl[ \mathcal{L}_u u_z \biggr] ~+ \underbrace{\bold{\hat{e}}_\varphi \biggl( \frac{ u_\varpi u_\varphi}{\varpi} \biggr) -~ \bold{\hat{e}}_\varpi \biggl( \frac{u_\varphi^2}{\varpi} \biggr) }_\text{curvature terms} \, , </math>

where the operator,

<math>~\mathcal{L}_u</math>

<math>~\equiv</math>

<math>~ \biggl[ u_\varpi \frac{\partial }{\partial \varpi} + \frac{u_\varphi }{\varpi}\frac{\partial }{\partial \varphi} + u_z \frac{\partial }{\partial z} \biggr] \, . </math>

In numerical simulations that are carried out on a cylindrical grid and in a rotating reference frame, it is customary to group the "curvature terms" with the fictitious acceleration to obtain,

<math>~\frac{\partial \bold{u}}{\partial t} + \bold{\hat{e}}_\varpi \biggl[ \mathcal{L}_u u_\varpi \biggr] + \bold{\hat{e}}_\varphi \biggl[ \mathcal{L}_u u_\varphi \biggr] + \bold{\hat{e}}_z \biggl[ \mathcal{L}_u u_z \biggr] </math>

<math>~=</math>

<math>~ \bold{a} + \bold{a}_\mathrm{fict} ~- \bold{\hat{e}}_\varphi \biggl( \frac{ u_\varpi u_\varphi}{\varpi} \biggr) +~ \bold{\hat{e}}_\varpi \biggl( \frac{u_\varphi^2}{\varpi} \biggr) </math>

 

<math>~=</math>

<math>~ \bold{a} + \bold{\hat{e}}_\varpi \biggl[2\Omega_f v_\varphi - \varpi \Omega_f^2 + \biggl( \frac{u_\varphi^2}{\varpi} \biggr) \biggr] - \bold{\hat{e}}_\varphi \biggl[ 2\Omega_f v_\varpi ~+ \biggl( \frac{ u_\varpi u_\varphi}{\varpi} \biggr)\biggr] </math>

 

<math>~=</math>

<math>~ \bold{a} + \bold{\hat{e}}_\varpi \biggl[2\Omega_f (u_\varphi + \varpi \Omega_f) - \varpi \Omega_f^2 + \biggl( \frac{u_\varphi^2}{\varpi} \biggr) \biggr] - \bold{\hat{e}}_\varphi \biggl[ 2\varpi \Omega_f + u_\varphi \biggr] \frac{u_\varpi}{\varpi} </math>

 

<math>~=</math>

<math>~ \bold{a} + \bold{\hat{e}}_\varpi \biggl[ \varpi\Omega_f + u_\varphi \biggr]^2 \frac{1}{\varpi} - \bold{\hat{e}}_\varphi \biggl[ 2\varpi \Omega_f + u_\varphi \biggr] \frac{u_\varpi}{\varpi} \, , </math>

treating the ensemble as an additional "source" of acceleration.

Example from the Literature
(see an accompanying related discussion)

We begin with the version of the Euler equation that has just been derived, namely,

<math>~\frac{\partial \bold{u}}{\partial t} + \bold{\hat{e}}_\varpi \biggl[ \mathcal{L}_u u_\varpi \biggr] + \bold{\hat{e}}_\varphi \biggl[ \mathcal{L}_u u_\varphi \biggr] + \bold{\hat{e}}_z \biggl[ \mathcal{L}_u u_z \biggr] </math>

<math>~=</math>

<math>~ \bold{a} + \bold{\hat{e}}_\varpi \biggl[ \varpi\Omega_f + u_\varphi \biggr]^2 \frac{1}{\varpi} - \bold{\hat{e}}_\varphi \biggl[ 2\varpi \Omega_f + u_\varphi \biggr] \frac{u_\varpi}{\varpi} \, . </math>



In examining and rearranging terms in each of the three components of this Euler equation, we will recognize that,

<math>~\mathcal{L}_u u_i</math>

<math>~=</math>

<math>~(\bold{u} \cdot \bold\nabla)u_i \, ;</math>

and that, after multiplying the standard Lagrangian representation of the continuity equation through by <math>~u_i</math>, we have,

<math>~0</math>

<math>~=</math>

<math>~ u_i \biggl[\frac{d \rho}{dt} + \rho \bold\nabla \cdot \bold{v} \biggr] = u_i \biggl[\frac{d \rho}{dt} + \rho \bold\nabla \cdot \bold{u} + \rho \cancelto{0}{\bold\nabla \cdot (\bold{\hat{e}}_\varphi \varpi \Omega_f)} \biggr] </math>

 

<math>~=</math>

<math>~u_i \biggl[\frac{\partial \rho}{\partial t} + \bold\nabla \cdot (\rho \bold{u}) \biggr]</math>

 

<math>~=</math>

<math>~ \frac{\partial \rho u_i}{\partial t} - \rho \frac{\partial u_i}{\partial t} + \bold\nabla \cdot (\rho u_i \bold{u}) - \rho (\bold{u}\cdot \bold\nabla) u_i </math>

<math>~\Rightarrow~~~ \rho \frac{\partial u_i}{\partial t} + \rho \biggl[ \mathcal{L}_u u_i \biggr] </math>

<math>~=</math>

<math>~ \frac{\partial \rho u_i}{\partial t} + \bold\nabla \cdot (\rho u_i \bold{u}) \, . </math>



Vertical Component:  Multiplying the <math>~\bold{\hat{k}}</math> component of this Euler equation through by <math>~\rho</math>, gives,

<math>~\rho \frac{\partial u_z}{\partial t} + \rho \biggl[ \mathcal{L}_u u_z \biggr] </math>

<math>~=</math>

<math>~\rho a_z </math>

<math>~\Rightarrow~~~ \frac{\partial \rho u_z}{\partial t} + \bold\nabla \cdot (\rho u_z \bold{u}) </math>

<math>~=</math>

<math>~\rho a_z \, .</math>

Norman & Wilson (1978), ApJ, 224, pp. 497 - 511, §III.b, Eq. (5)
New & Tohline (1997), ApJ, 490, pp. 311 - 237, §2, Eq. (3)

Radial Component:  Multiplying the <math>~\bold{\hat{e}}_\varpi</math> component of this Euler equation through by <math>~\rho</math>, gives,

<math>~\rho\frac{\partial u_\varpi}{\partial t} + \rho\biggl[ \mathcal{L}_u u_\varpi \biggr] </math>

<math>~=</math>

<math>~ \rho a_\varpi + \biggl[u_\varphi + \Omega_f \varpi \biggr]^2 \frac{\rho}{\varpi} </math>

<math>~\Rightarrow ~~~ \frac{\partial \rho u_\varpi}{\partial t} + \bold\nabla \cdot (\rho u_\varpi \bold{u})</math>

<math>~=</math>

<math>~ \rho a_\varpi + \biggl[\frac{(\rho \varpi u_\varphi)^2}{\rho \varpi^3} + \rho\Omega_f^2 \varpi + \frac{ 2\Omega_f (\rho \varpi u_\varphi )}{\varpi} \biggr] \, . </math>

Norman & Wilson (1978), ApJ, 224, pp. 497 - 511, §III.b, Eq. (6)
New & Tohline (1997), ApJ, 490, pp. 311 - 237, §2.2, Eq. (11)

Azimuthal Component:  Multiplying the <math>~\bold{\hat{e}}_\varphi</math> component of this Euler equation through by <math>~\rho</math>, gives,

<math>~\rho\frac{\partial u_\varphi}{\partial t} + \rho\biggl[ \mathcal{L}_u u_\varphi \biggr] </math>

<math>~=</math>

<math>~ \rho a_\varphi - \biggl[ 2\varpi \Omega_f + u_\varphi \biggr] \frac{\rho u_\varpi}{\varpi} </math>

<math>~\Rightarrow ~~~ \frac{\partial \rho u_\varphi}{\partial t} + \bold\nabla \cdot (\rho u_\varphi \bold{u})</math>

<math>~=</math>

<math>~ \rho a_\varphi - \biggl[ 2\varpi \Omega_f + u_\varphi \biggr] \frac{\rho u_\varpi}{\varpi} \, . </math>

Then, multiplying through by <math>~\varpi</math>, gives,

<math>~ \varpi \frac{\partial \rho u_\varphi}{\partial t} + \varpi \bold\nabla \cdot (\rho u_\varphi \bold{u})</math>

<math>~=</math>

<math>~ \rho \varpi a_\varphi - \biggl[ 2\varpi \Omega_f + u_\varphi \biggr] \rho u_\varpi </math>

<math>~\Rightarrow~~~ \frac{\partial \rho \varpi u_\varphi}{\partial t} + \bold\nabla \cdot (\rho \varpi u_\varphi \bold{u}) - \rho u_\varphi \cancelto{u_\varpi}{\biggl[ \mathcal{L}_u \varpi \biggr]} </math>

<math>~=</math>

<math>~ \rho \varpi a_\varphi - 2\Omega_f \varpi \rho u_\varpi - \rho u_\varphi u_\varpi </math>

<math>~\Rightarrow~~~ \frac{\partial \rho \varpi u_\varphi}{\partial t} + \bold\nabla \cdot (\rho \varpi u_\varphi \bold{u}) </math>

<math>~=</math>

<math>~ \rho \varpi a_\varphi - 2\Omega_f \varpi \rho u_\varpi \, . </math>

Norman & Wilson (1978), ApJ, 224, pp. 497 - 511, §III.b, Eq. (7)
New & Tohline (1997), ApJ, 490, pp. 311 - 237, §2.2, Eq. (12)

Hybrid Scheme

In our hybrid scheme, we will continue to use <math>~\mathcal{L}_u</math> — that is, an advection operator that incorporates the rotating-frame velocity, <math>~\bold{u}</math> — but we will switch all other velocity references to the inertial-frame velocity, <math>~\bold{v}</math>, and its components. This will be done via the above-declared mapping,

<math>~\bold{u}</math>

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

<math>~[\bold{v} - \bold{\hat{e}}_\varphi \varpi \Omega_f] \, ,</math>

that is, <math>~u_\varpi ~\rightarrow~ v_\varpi</math>, <math>~u_z ~\rightarrow~ v_z</math>, and <math>~u_\varphi ~\rightarrow~ (v_\varphi - \varpi \Omega_f)</math>. The Euler equation becomes,

<math>~\frac{\partial [\bold{v} - \bold{\hat{e}}_\varphi \varpi \Omega_f]}{\partial t} + \bold{\hat{e}}_\varpi \biggl[ \mathcal{L}_u v_\varpi \biggr] + \bold{\hat{e}}_\varphi \biggl[ \mathcal{L}_u (v_\varphi - \varpi \Omega_f)\biggr] + \bold{\hat{e}}_z \biggl[ \mathcal{L}_u v_z \biggr] </math>

<math>~=</math>

<math>~ \bold{a} + \bold{\hat{e}}_\varpi \biggl[ \frac{v_\varphi^2}{\varpi} \biggr] - \bold{\hat{e}}_\varphi \biggl[ \varpi \Omega_f + v_\varphi \biggr] \frac{v_\varpi}{\varpi} \, , </math>

where we recognize that,

<math>~\mathcal{L}_u v_i</math>

<math>~=</math>

<math>~(\bold{u}\cdot \bold\nabla) v_i \, .</math>

Now, given that,

<math>~\frac{\partial}{\partial t} \biggl[ \bold{\hat{e}}_\varphi \varpi \Omega_f \biggr]</math>

<math>~=</math>

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

and,

<math>~\mathcal{L}_u (\varpi \Omega_f)</math>

<math>~=</math>

<math>~\biggl[ u_\varpi \frac{\partial }{\partial \varpi} + \frac{u_\varphi }{\varpi}\frac{\partial }{\partial \varphi} + u_z \frac{\partial }{\partial z} \biggr] (\varpi \Omega_f) = u_\varpi \Omega_f = v_\varpi \Omega_f \, , </math>

the Euler equation becomes,

<math>~\frac{\partial \bold{v} }{\partial t} + \bold{\hat{e}}_\varpi \biggl[ \mathcal{L}_u v_\varpi \biggr] + \bold{\hat{e}}_\varphi \biggl[ \mathcal{L}_u v_\varphi - v_\varpi \Omega_f \biggr] + \bold{\hat{e}}_z \biggl[ \mathcal{L}_u v_z \biggr] </math>

<math>~=</math>

<math>~ \bold{a} + \bold{\hat{e}}_\varpi \biggl[ \frac{v_\varphi^2}{\varpi} \biggr] - \bold{\hat{e}}_\varphi \biggl[ \varpi \Omega_f + v_\varphi \biggr] \frac{v_\varpi}{\varpi} </math>

<math>~\Rightarrow ~~~ \frac{\partial \bold{v} }{\partial t} + \bold{\hat{e}}_\varpi \biggl[ \mathcal{L}_u v_\varpi \biggr] + \bold{\hat{e}}_\varphi \biggl[ \mathcal{L}_u v_\varphi \biggr] + \bold{\hat{e}}_z \biggl[ \mathcal{L}_u v_z \biggr] </math>

<math>~=</math>

<math>~ \bold{a} + \bold{\hat{e}}_\varpi \biggl[ \frac{v_\varphi^2}{\varpi} \biggr] - \bold{\hat{e}}_\varphi \biggl[ \frac{v_\varpi v_\varphi}{\varpi} \biggr] \, . </math>

Also, if we multiply the standard Lagrangian representation of the continuity equation through by <math>~v_i</math>, we have,

<math>~0</math>

<math>~=</math>

<math>~ v_i \biggl[\frac{d \rho}{dt} + \rho \bold\nabla \cdot \bold{v} \biggr] = v_i \biggl[\frac{d \rho}{dt} + \rho \bold\nabla \cdot \bold{u} + \rho \cancelto{0}{\bold\nabla \cdot (\bold{\hat{e}}_\varphi \varpi \Omega_f)} \biggr] </math>

 

<math>~=</math>

<math>~v_i \biggl[\frac{\partial \rho}{\partial t} + \bold\nabla \cdot (\rho \bold{u}) \biggr]</math>

 

<math>~=</math>

<math>~ \frac{\partial \rho v_i}{\partial t} - \rho \frac{\partial v_i}{\partial t} + \bold\nabla \cdot (\rho v_i \bold{u}) - \rho (\bold{u}\cdot \bold\nabla) v_i </math>

<math>~\Rightarrow~~~ \rho \frac{\partial v_i}{\partial t} + \rho \biggl[ \mathcal{L}_u v_i \biggr] </math>

<math>~=</math>

<math>~ \frac{\partial \rho v_i}{\partial t} + \bold\nabla \cdot (\rho v_i \bold{u}) \, . </math>



Vertical Component:  Multiplying the <math>~\bold{\hat{k}}</math> component of our modified Euler equation through by <math>~\rho</math>, then incorporating this version of the continuity equation, gives,

<math>~\rho \frac{\partial v_z}{\partial t} + \rho \biggl[ \mathcal{L}_u v_z \biggr] </math>

<math>~=</math>

<math>~\rho a_z </math>

<math>~\Rightarrow~~~ \frac{\partial \rho v_z}{\partial t} + \bold\nabla \cdot (\rho v_z \bold{u}) </math>

<math>~=</math>

<math>~\rho a_z \, .</math>

Radial Component:  Multiplying the <math>~\bold{\hat{e}}_\varpi</math> component of our modified Euler equation through by <math>~\rho</math>, then incorporating the continuity equation, gives,

<math>~\rho \frac{\partial v_\varpi}{\partial t} + \rho \biggl[ \mathcal{L}_u v_\varpi \biggr] </math>

<math>~=</math>

<math>~\rho a_\varpi + \frac{v_\varphi^2}{\varpi} </math>

<math>~\Rightarrow~~~ \frac{\partial \rho v_\varpi}{\partial t} + \bold\nabla \cdot (\rho v_\varpi \bold{u}) </math>

<math>~=</math>

<math>~\rho a_\varpi + \frac{v_\varphi^2}{\varpi} \, .</math>

Azimuthal Component:  Multiplying the <math>~\bold{\hat{e}}_\varphi</math> component of our modified Euler equation through by <math>~\rho</math>, then incorporating the continuity equation, gives,

<math>~\rho \frac{\partial v_\varphi}{\partial t} + \rho \biggl[ \mathcal{L}_u v_\varphi \biggr] </math>

<math>~=</math>

<math>~\rho a_\varphi - \frac{\rho v_\varphi v_\varpi}{\varpi} </math>

<math>~\Rightarrow~~~ \frac{\partial \rho v_\varphi}{\partial t} + \bold\nabla \cdot (\rho v_\varphi \bold{u}) </math>

<math>~=</math>

<math>~\rho a_\varphi - \frac{\rho v_\varphi v_\varpi}{\varpi} \, .</math>

Then, multiplying through by <math>~\varpi</math>, we have,

<math>~\rho \varpi a_\varphi - \rho v_\varphi v_\varpi</math>

<math>~=</math>

<math>~ \varpi \frac{\partial \rho v_\varphi}{\partial t} + \varpi \bold\nabla \cdot (\rho v_\varphi \bold{u}) </math>

 

<math>~=</math>

<math>~ \frac{\partial(\rho \varpi v_\varphi )}{\partial t} - (\rho v_\varphi )\cancelto{0}{\frac{\partial \varpi}{\partial t}} + \bold\nabla \cdot (\rho \varpi v_\varphi \bold{u}) - \rho v_\varphi (\bold{u}\cdot \bold\nabla) \varpi </math>

 

<math>~=</math>

<math>~ \frac{\partial(\rho \varpi v_\varphi )}{\partial t} + \bold\nabla \cdot (\rho \varpi v_\varphi \bold{u}) - \rho v_\varphi \biggl[ \cancelto{v_\varpi}{\mathcal{L}_u \varpi} \biggr] </math>

<math>~ \Rightarrow~~~ \frac{\partial(\rho \varpi v_\varphi )}{\partial t} + \bold\nabla \cdot (\rho \varpi v_\varphi \bold{u}) </math>

<math>~=</math>

<math>~\rho \varpi a_\varphi \, .</math>

Compare

Here we consider which formalism is best suited for modeling a fully three-dimensional, nonaxisymmetric configuration that is spinning about (usually) its shortest axis with a uniform and time-independent frequency and which, when viewed from a frame that is rotating with that frequency, exhibits a nontrivial but nevertheless steady-state internal flow. Examples are Riemann S-type ellipsoids, and binary stars in circular orbits.

It is most desirable to choose a formalism that recognizes the steady-state nature of the flow. In the vast majority of cases being considered here, this rules out using any scheme that is designed around an inertial-frame coordinate base. (As a counterexample, Dedekind ellipsoids can be constructed in the inertial frame because <math>~\Omega_f = 0</math> for all models along the Dedekind equilibrium sequence.) It is quite reasonable, however, to adopt a rotating, cylindrical-coordinate base as has been described above and as is summarized immediately below.

Traditional Rotating, Cylindrical-Coordinate Summary
<math>~\bold{\hat{k}:}</math>

<math>~ \frac{\partial (\rho u_z)}{\partial t} + \bold\nabla \cdot (\rho u_z \bold{u}) </math>

<math>~=</math>

<math>~\bold{\hat{k}}\cdot (\rho \bold{a}) \, ;</math>

<math>~\bold{\hat{e}_\varpi:}</math>

<math>~ \frac{\partial ( \rho u_\varpi )}{\partial t} + \bold\nabla \cdot (\rho u_\varpi \bold{u})</math>

<math>~=</math>

<math>~\bold{\hat{e}}_\varpi \cdot (\rho \bold{a}) + \biggl[\frac{(\rho \varpi u_\varphi)^2}{\rho \varpi^3} + \rho\Omega_f^2 \varpi + \frac{ 2\Omega_f (\rho \varpi u_\varphi )}{\varpi} \biggr] \, ; </math>

<math>~\bold{\hat{e}_\varphi:}</math>

<math>~ \frac{\partial (\rho \varpi u_\varphi )}{\partial t} + \bold\nabla \cdot (\rho \varpi u_\varphi \bold{u}) </math>

<math>~=</math>

<math>~\bold{\hat{e}}_\varphi \cdot (\rho \varpi \bold{a}) - 2\Omega_f \varpi \rho u_\varpi \, . </math>

In this scheme, all of the velocities and associated momentum densities in all three components of the Euler equation are expressed in terms of the rotating-frame velocity vector, <math>~\bold{u}</math>, or its cylindrical-coordinate-based components, <math>~(u_\varpi, v_\varphi, v_z)</math>. When the configuration's distorted (nonaxisymmetric) shape is largely supported by rapid rotation, this scheme provides an advantage over other — for example, inertial-frame-based — schemes because the fraction of the fluid's total momentum that is being advected across the grid is often quite small. There is a penalty to be paid, however. Additional "source" terms appear on the right-hand-side of the radial- and azimuthal-component expressions; they are nonlinear in the velocity and introduce cross-talk between the component expressions.

Hybrid Scheme Summary
<math>~\bold{\hat{k}:}</math>

<math>~ \frac{\partial (\rho v_z)}{\partial t} + \bold\nabla \cdot (\rho v_z \bold{u}) </math>

<math>~=</math>

<math>~\bold{\hat{k}} \cdot (\rho \bold{a}) \, ;</math>

<math>~\bold{\hat{e}_\varpi:}</math>

<math>~ \frac{\partial (\rho v_\varpi)}{\partial t} + \bold\nabla \cdot (\rho v_\varpi \bold{u}) </math>

<math>~=</math>

<math>~\bold{\hat{e}}_\varpi \cdot (\rho \bold{a}) + \frac{v_\varphi^2}{\varpi} \, ;</math>

<math>~\bold{\hat{e}_\varphi:}</math>

<math>~ \frac{\partial(\rho \varpi v_\varphi )}{\partial t} + \bold\nabla \cdot (\rho \varpi v_\varphi \bold{u}) </math>

<math>~=</math>

<math>~\bold{\hat{e}}_\varphi \cdot (\rho \varpi \bold{a}) \, .</math>

For Riemann S-type ellipsoids,

<math>~u_x</math>

<math>~=</math>

<math>~ \lambda \biggl(\frac{a}{b}\biggr) y \, , </math>

<math>~u_y</math>

<math>~=</math>

<math>~ - \lambda \biggl(\frac{b}{a}\biggr) x \, , </math>

<math>~v_\varpi</math>

<math>~=</math>

<math>~ \lambda \biggl[\frac{a}{b} - \frac{b}{a} \biggr] x y (x^2 + y^2)^{-1 / 2} \, , </math>

<math>~\varpi v_\varphi</math>

<math>~=</math>

<math>~ - \biggl[ \lambda \biggl(\frac{b}{a}\biggr) - \Omega_f\biggr]x^2 - \biggl[ \lambda \biggl(\frac{a}{b}\biggr) - \Omega_f\biggr]y^2 \, . </math>

Hence,

<math>~\bold\nabla \cdot (v_\varpi \bold{u})</math>

<math>~=</math>

<math>~ \frac{\partial }{\partial x}\biggl\{ v_\varpi u_x \biggr\} + \frac{\partial }{\partial y}\biggl\{ v_\varpi u_y \biggr\} </math>

 

<math>~=</math>

<math>~ \lambda^2 \biggl(\frac{a}{b}\biggr) \biggl[\frac{a}{b} - \frac{b}{a} \biggr] \frac{\partial }{\partial x}\biggl\{ x y^2 (x^2 + y^2)^{-1 / 2} \biggr\} - \lambda^2 \biggl(\frac{b}{a}\biggr) \biggl[\frac{a}{b} - \frac{b}{a} \biggr] \frac{\partial }{\partial y}\biggl\{ x^2 y (x^2 + y^2)^{-1 / 2} \biggr\} </math>

 

<math>~=</math>

<math>~ y^2 \lambda^2 \biggl(\frac{a}{b}\biggr) \biggl[\frac{a}{b} - \frac{b}{a} \biggr] \biggl\{ (x^2 + y^2)^{-1 / 2} - x^2 (x^2 + y^2)^{-3 / 2}\biggr\} - x^2 \lambda^2 \biggl(\frac{b}{a}\biggr) \biggl[\frac{a}{b} - \frac{b}{a} \biggr] \biggl\{ (x^2 + y^2)^{-1 / 2} - y^2 (x^2 + y^2)^{-3 / 2} \biggr\} </math>

 

<math>~=</math>

<math>~\frac{\lambda^2}{\varpi^3} \biggl[\frac{a}{b} - \frac{b}{a} \biggr] \biggl\{ y^4 \biggl(\frac{a}{b}\biggr) - x^4 \biggl(\frac{b}{a}\biggr) \biggr\} </math>

And,

<math>~\bold\nabla \cdot (\varpi v_\varphi \bold{u})</math>

<math>~=</math>

<math>~ \frac{\partial}{\partial x} \biggl\{- \lambda\biggl(\frac{a}{b}\biggr)y\biggl[ \lambda \biggl(\frac{b}{a}\biggr) - \Omega_f\biggr]x^2 \biggr\} + \frac{\partial}{\partial y} \biggl\{ \lambda\biggl(\frac{b}{a}\biggr)x \biggl[ \lambda \biggl(\frac{a}{b}\biggr) - \Omega_f\biggr]y^2 \biggr\} </math>

 

<math>~=</math>

<math>~ \biggl\{- 2 \lambda\biggl(\frac{a}{b}\biggr)\biggl[ \lambda \biggl(\frac{b}{a}\biggr) - \Omega_f\biggr]x y \biggr\} + \biggl\{ 2 \lambda\biggl(\frac{b}{a}\biggr)\biggl[ \lambda \biggl(\frac{a}{b}\biggr) - \Omega_f\biggr]x y \biggr\} </math>

 

<math>~=</math>

<math>~ \biggl\{ -\biggl(\frac{a}{b}\biggr)\biggl[ \lambda \biggl(\frac{b}{a}\biggr) - \Omega_f\biggr] + \biggl(\frac{b}{a}\biggr)\biggl[ \lambda \biggl(\frac{a}{b}\biggr) - \Omega_f\biggr] \biggr\} 2 \lambda x y </math>

 

<math>~=</math>

<math>~ \biggl\{-\biggl[ \lambda - \biggl(\frac{a}{b}\biggr)\Omega_f\biggr] + \biggl[ \lambda - \biggl(\frac{b}{a}\biggr)\Omega_f\biggr] \biggr\} 2 \lambda x y </math>

 

<math>~=</math>

<math>~ 2 \lambda xy \Omega_f \biggl[\frac{a}{b} - \frac{b}{a} \biggr] \, . </math>

To within an additive constant — see, for example, our associated discussion of Maclaurin spheroids — the gravitational potential and the enthalpy are, respectively,

<math>~\Phi_\mathrm{grav}</math>

<math>~=</math>

<math>~\pi G \rho \biggl[ A_1 x^2 + A_2 y^2 + A_3 z^2 \biggr] = \pi G \rho \biggl[ A_1 \varpi^2 \cos^2\varphi + A_2 \varpi^2 \sin^2\varphi + A_3 z^2 \biggr] \, , </math>

<math>~H</math>

<math>~=</math>

<math>~ H_0 \biggl[1 - \frac{x^2}{a^2} - \frac{y^2}{b^2} - \frac{z^2}{c^2} \biggr] = H_0 \biggl[1 - \frac{\varpi^2\cos^2\varphi}{a^2} - \frac{\varpi^2 \sin^2\varphi}{b^2} - \frac{z^2}{c^2} \biggr] \, .</math>

Hence,

<math>~\bold{a} = - \bold\nabla (H + \Phi_\mathrm{grav} )</math>

<math>~=</math>

<math>~-\biggl[ \bold{\hat{e}}_\varpi \frac{\partial}{\partial \varpi} \biggl( H + \Phi_\mathrm{grav}\biggr) + \frac{\bold{\hat{e}}_\varphi}{\varpi} \frac{\partial}{\partial \varphi} \biggl( H + \Phi_\mathrm{grav}\biggr) + \bold{\hat{k}} \frac{\partial}{\partial z} \biggl( H + \Phi_\mathrm{grav}\biggr) \biggr] </math>

 

<math>~=</math>

<math>~ - \bold{\hat{e}}_\varpi \biggl[ -2H_0 \varpi \biggl( \frac{\cos^2\varphi}{a^2} + \frac{\sin^2\varphi}{b^2}\biggr) + 2\pi G \rho \varpi \biggl(A_1\cos^2\varphi + A_2\sin^2\varphi \biggr) \biggr] </math>

 

 

<math>~ - \frac{\bold{\hat{e}}_\varphi}{\varpi} \biggl[ -2H_0\varpi^2 \sin\varphi \cos\varphi\biggl( \frac{1}{b^2} - \frac{1}{a^2}\biggr) + 2\pi G \rho \varpi^2 \sin\varphi \cos\varphi \biggl( A_2 - A_1 \biggr) \biggr] </math>

 

 

<math>~ - \bold{\hat{k}} \biggl[ - \frac{2H_0 z}{c^2} + 2\pi G \rho A_3 z \biggr] \, . </math>



Vertical Component:   Because <math>~v_z = 0</math>, it must be true that <math>~\bold{\hat{k}}\cdot \bold{a} = 0</math>. This, in turn means,

<math>~H_0 </math>

<math>~=</math>

<math>~\pi G \rho c^2 A_3 \, . </math>

Azimuthal Component:   In steady-state, the partial time-derivative must be zero, so we require,

<math>~ \bold\nabla \cdot (\varpi v_\varphi \bold{u}) </math>

<math>~=</math>

<math>~\varpi \bold{\hat{e}}_\varphi \cdot \bold{a}</math>

<math>~ \Rightarrow~~~ 2 \lambda xy \Omega_f \biggl[\frac{a}{b} - \frac{b}{a} \biggr] </math>

<math>~=</math>

<math>~

\biggl[ 2H_0\varpi^2 \sin\varphi \cos\varphi\biggl( \frac{1}{b^2} - \frac{1}{a^2}\biggr) 

- 2\pi G \rho \varpi^2 \sin\varphi \cos\varphi \biggl( A_2 - A_1 \biggr) \biggr] </math>

 

<math>~=</math>

<math>~ 2xy \biggl[ H_0\biggl( \frac{1}{b^2} - \frac{1}{a^2}\biggr) - \pi G \rho \biggl( A_2 - A_1 \biggr) \biggr] </math>

 

<math>~=</math>

<math>~ 2\pi G \rho xy \biggl[ c^2 A_3\biggl( \frac{1}{b^2} - \frac{1}{a^2}\biggr) - \biggl( A_2 - A_1 \biggr) \biggr] </math>

<math>~ \Rightarrow~~~ \lambda \Omega_f \biggl[\frac{a^2 - b^2}{ab} \biggr] </math>

<math>~=</math>

<math>~ \pi G \rho \biggl[ c^2 A_3\biggl( \frac{a^2 - b^2}{a^2 b^2} \biggr) + \biggl( A_1 - A_2 \biggr) \biggr] </math>

 

<math>~=</math>

<math>~ \pi G \rho \biggl[ ( A_1 - A_2 ) - c^2 A_3\biggl( \frac{b^2 - a^2}{a^2 b^2} \biggr) \biggr] </math>

<math>~ \Rightarrow~~~ - a b \lambda \Omega_f </math>

<math>~=</math>

<math>~ \pi G \rho \biggl[ \frac{( A_1 - A_2 )a^2b^2}{b^2 - a^2} - c^2 A_3 \biggr] \, . </math>

This expression can be used either (a) to give <math>~\Omega_f</math> in terms of a set of known quantities and the unknown parameter, <math>~\lambda</math>; or (b) to give <math>~\lambda</math> in terms of a set of known quantities and the unknown parameter, <math>~\Omega_f</math>.

Recognizing from here that, <math>~f \equiv \zeta/\Omega_f</math> and

<math>~\lambda \biggl( \frac{a}{b} + \frac{b}{a} \biggr)</math>

<math>~=</math>

<math>~ - \zeta = - f\Omega_f </math>

<math>~\Rightarrow ~~~ \lambda </math>

<math>~=</math>

<math>~ - \biggl( \frac{ab}{a^2 + b^2} \biggr)f \Omega_f \, , </math>

this last expression can be rewritten as,

<math>~ \biggl( \frac{a^2 b^2}{a^2 + b^2} \biggr)f \Omega_f^2 </math>

<math>~=</math>

<math>~ \pi G \rho \biggl[ \frac{( A_1 - A_2 )a^2b^2}{b^2 - a^2} - c^2 A_3 \biggr] \, . </math>

[ EFE, Chapter 7, §48, Eq. (34) ]

Radial Component:   In steady-state, this partial time-derivative also must be zero, so we require,

<math>~ \bold\nabla \cdot (v_\varpi \bold{u}) </math>

<math>~=</math>

<math>~\bold{\hat{e}}_\varpi \cdot \bold{a} + \frac{v_\varphi^2}{\varpi} </math>

<math>~\Rightarrow ~~~ \frac{\lambda^2}{\varpi^3} \biggl[\frac{a}{b} - \frac{b}{a} \biggr] \biggl\{ y^4 \biggl(\frac{a}{b}\biggr) - x^4 \biggl(\frac{b}{a}\biggr) \biggr\} </math>

<math>~=</math>

<math>~\frac{2\pi G \rho }{\varpi}\biggl[ c^2 A_3 \biggl( \frac{x^2}{a^2} + \frac{y^2}{b^2}\biggr) - \biggl(A_1 x^2 + A_2y^2 \biggr) \biggr] + \frac{1}{\varpi^3} \biggl\{ -\biggl[\lambda \biggl(\frac{b}{a}\biggr) - \Omega_f\biggr]x^2 - \biggl[ \lambda \biggl(\frac{a}{b}\biggr) - \Omega_f\biggr]y^2 \biggr\}^2 </math>

<math>~\Rightarrow ~~~ \lambda^2 \biggl[\frac{a}{b} - \frac{b}{a} \biggr] \biggl[ y^4 \biggl(\frac{a}{b}\biggr) - x^4 \biggl(\frac{b}{a}\biggr) \biggr] </math>

<math>~=</math>

<math>~2\pi G \rho (x^2 + y^2) \biggl[ x^2 \biggl( \frac{c^2 A_3}{a^2} - A_1 \biggr) + y^2 \biggl( \frac{c^2 A_3}{b^2} - A_2 \biggr) \biggr] + \biggl\{ \biggl[ \lambda \biggl(\frac{b}{a}\biggr) - \Omega_f\biggr]x^2 + \biggl[ \lambda \biggl(\frac{a}{b}\biggr) - \Omega_f\biggr]y^2 \biggr\}^2 </math>

 

<math>~=</math>

<math>~ 2\pi G \rho \biggl\{ x^4 \biggl( \frac{c^2 A_3}{a^2} - A_1 \biggr) +x^2 y^2 \biggl( \frac{c^2 A_3}{b^2} - A_2 + \frac{c^2 A_3}{a^2} - A_1 \biggr) +y^4 \biggl( \frac{c^2 A_3}{b^2} - A_2 \biggr) \biggr\} </math>

 

 

<math>~ + \biggl[ \lambda \biggl(\frac{b}{a}\biggr) - \Omega_f\biggr]^2x^4 + 2x^2 y^2 \biggl[ \lambda \biggl(\frac{b}{a}\biggr) - \Omega_f\biggr] \biggl[ \lambda \biggl(\frac{a}{b}\biggr) - \Omega_f\biggr] + \biggl[ \lambda \biggl(\frac{a}{b}\biggr) - \Omega_f\biggr]^2y^4 </math>

<math>~\Rightarrow ~~~ 2\pi G \rho \biggl\{ x^4 \biggl( \frac{c^2 A_3}{a^2} - A_1 \biggr) +y^4 \biggl( \frac{c^2 A_3}{b^2} - A_2 \biggr) \biggr\} </math>

<math>~=</math>

<math>~ \lambda^2 \biggl[\frac{a^2}{b^2} - 1 \biggr] y^4 +\biggl[ \lambda \biggl(\frac{a}{b}\biggr) - \Omega_f\biggr]^2y^4 - \lambda^2 \biggl[1 - \frac{b^2}{a^2} \biggr] x^4 - \biggl[ \lambda \biggl(\frac{b}{a}\biggr) - \Omega_f\biggr]^2x^4 </math>

 

 

<math>~ + 2x^2 y^2 \biggl\{ \biggl[ \lambda \biggl(\frac{b}{a}\biggr) - \Omega_f\biggr] \biggl[ \lambda \biggl(\frac{a}{b}\biggr) - \Omega_f\biggr] + \pi G \rho \biggl[ \frac{c^2 A_3 (a^2 + b^2)}{a^2b^2} - A_1 - A_2 \biggr] \biggr\} </math>

<math>~\Rightarrow ~~~ x^4 \biggl\{ 2\pi G \rho \biggl( \frac{c^2 A_3}{a^2} - A_1 \biggr) + \lambda^2 \biggl[1 - \frac{b^2}{a^2} \biggr] + \biggl[ \lambda^2 \biggl(\frac{b^2}{a^2}\biggr) - 2\lambda\biggl(\frac{b}{a}\biggr) \Omega_f + \Omega_f^2 \biggr] \biggr\} </math>

<math>~=</math>

<math>~y^4 \biggl\{ \lambda^2 \biggl[\frac{a^2}{b^2} - 1 \biggr] + \biggl[\lambda^2\biggl( \frac{a^2}{b^2}\biggr) - 2\lambda \Omega_f \biggl(\frac{a}{b}\biggr) + \Omega_f^2 \biggr] -2\pi G \rho \biggl( \frac{c^2 A_3}{b^2} - A_2 \biggr) \biggr\} </math>

 

 

<math>~ + 2x^2 y^2 \biggl\{ \biggl[ \lambda \biggl(\frac{b}{a}\biggr) - \Omega_f\biggr] \biggl[ \lambda \biggl(\frac{a}{b}\biggr) - \Omega_f\biggr] + \pi G \rho \biggl[ \frac{c^2 A_3 (a^2 + b^2)}{a^2b^2} - A_1 - A_2 \biggr] \biggr\} </math>

Note that,

<math>~ - a b \lambda \Omega_f </math>

<math>~=</math>

<math>~ \pi G \rho \biggl[ \frac{( A_1 - A_2 )a^2b^2}{b^2 - a^2} - c^2 A_3 \biggr] </math>

<math>~ \Rightarrow~~~ c^2 A_3 </math>

<math>~=</math>

<math>~ ab \biggl( \frac{\lambda \Omega_f }{\pi G \rho} \biggr) + \frac{( A_1 - A_2 )a^2b^2}{b^2 - a^2} \, . </math>

Temporary

Radial Component:   In steady-state, this partial time-derivative also must be zero, so we require,

<math>~ \bold\nabla \cdot (v_\varpi \bold{u}) </math>

<math>~=</math>

<math>~\bold{\hat{e}}_\varpi \cdot \bold{a} + \frac{v_\varphi^2}{\varpi} </math>

<math>~\Rightarrow ~~~ \frac{\lambda^2}{\varpi^3} \biggl[\frac{a}{b} - \frac{b}{a} \biggr] \biggl\{ y^4 \biggl(\frac{a}{b}\biggr) - x^4 \biggl(\frac{b}{a}\biggr) \biggr\} </math>

<math>~=</math>

<math>~ \frac{2\pi G \rho }{\varpi}\biggl[ c^2 A_3 \biggl( \frac{x^2}{a^2} + \frac{y^2}{b^2}\biggr) - \biggl(A_1 x^2 + A_2y^2 \biggr) \biggr] + \frac{1}{\varpi^3} \biggl\{ -\biggl[\lambda \biggl(\frac{b}{a}\biggr) - \Omega_f\biggr]x^2 - \biggl[ \lambda \biggl(\frac{a}{b}\biggr) - \Omega_f\biggr]y^2 \biggr\}^2 </math>

<math>~\Rightarrow ~~~ \lambda^2 \biggl\{ a^2 y^4 - b^2 x^4 \biggr\} </math>

<math>~=</math>

<math>~ \frac{2\pi G \rho b^2 \varpi^2}{(a^2 - b^2)} \biggl[ c^2 A_3 x^2 - A_1 a^2 x^2 \biggr] + \frac{2\pi G \rho a^2 \varpi^2}{(a^2 - b^2)} \biggl[ c^2 A_3 y^2 - A_2 b^2 y^2 \biggr] +\biggl[\frac{a^2 b^2}{a^2 - b^2} \biggr] \biggl\{ \biggl[\lambda \biggl(\frac{b}{a}\biggr) - \Omega_f\biggr]x^2 + \biggl[ \lambda \biggl(\frac{a}{b}\biggr) - \Omega_f\biggr]y^2 \biggr\}^2 </math>

<math>~\Rightarrow ~~~ \lambda^2 (a^2 - b^2) \biggl\{ a^2 y^4 - b^2 x^4 \biggr\} </math>

<math>~=</math>

<math>~ 2\pi G \rho b^2 \varpi^2 x^2 \biggl[ c^2 A_3 - A_1 a^2 \biggr] + 2\pi G \rho a^2 \varpi^2 y^2 \biggl[ c^2 A_3 - A_2 b^2 \biggr] + \biggl\{ \biggl[\lambda b^2 - a b \Omega_f \biggr]x^2 + \biggl[ \lambda a^2 - a b \Omega_f\biggr]y^2 \biggr\}^2 \, . </math>

Also,

<math>~ - a b \lambda \Omega_f </math>

<math>~=</math>

<math>~ \pi G \rho \biggl[ \frac{( A_1 - A_2 )a^2b^2}{b^2 - a^2} - c^2 A_3 \biggr] </math>

<math>~\Rightarrow~~~ c^2 A_3 </math>

<math>~=</math>

<math>~ \frac{( A_1 - A_2 )a^2b^2}{b^2 - a^2} + \frac{a b \lambda \Omega_f }{\pi G \rho} \, . </math>

Hence,

<math>~ \lambda^2 (a^2 - b^2) \biggl\{ a^2 y^4 - b^2 x^4 \biggr\} - \biggl\{ \biggl[\lambda b^2 - a b \Omega_f \biggr]x^2 + \biggl[ \lambda a^2 - a b \Omega_f\biggr]y^2 \biggr\}^2 </math>

<math>~=</math>

<math>~ 2\pi G \rho b^2 \varpi^2 x^2 \biggl[ \frac{( A_1 - A_2 )a^2b^2}{b^2 - a^2} + \frac{a b \lambda \Omega_f }{\pi G \rho} - A_1 a^2 \biggr] + 2\pi G \rho a^2 \varpi^2 y^2 \biggl[ \frac{( A_1 - A_2 )a^2b^2}{b^2 - a^2} + \frac{a b \lambda \Omega_f }{\pi G \rho} - A_2 b^2 \biggr] </math>

 

<math>~=</math>

<math>~ \frac{2\pi G \rho b^2 \varpi^2 x^2 }{(b^2 - a^2) } \biggl[ ( A_1 - A_2 )a^2b^2 - A_1 a^2 (b^2 - a^2) \biggr] + \frac{ 2\pi G \rho a^2 \varpi^2 y^2 }{ (b^2 - a^2)} \biggl[ ( A_1 - A_2 )a^2b^2 - A_2 b^2(b^2 - a^2) \biggr] + 2 a b \lambda \Omega_f \varpi^2 (a^2 y^2 + b^2 x^2) </math>

 

<math>~=</math>

<math>~ \frac{2\pi G \rho a^2 b^2 \varpi^4 }{(b^2 - a^2) } \biggl[ A_1 a^2 - A_2 b^2 \biggr] + 2 a b \lambda \Omega_f \varpi^2 (a^2 y^2 + b^2 x^2) </math>

<math>~\Rightarrow ~~~ \frac{2\pi G \rho }{(a^2 - b^2) } \biggl[ A_1 a^2 - A_2 b^2 \biggr]a^2 b^2 \varpi^4 </math>

<math>~=</math>

<math>~ 2 a b \lambda \Omega_f (x^2 + y^2) (a^2 y^2 + b^2 x^2) - \lambda^2 (a^2 - b^2) [ a^2 y^4 - b^2 x^4 ] + \biggl\{ \biggl[\lambda b^2 - a b \Omega_f \biggr]^2 x^4 + 2\biggl[\lambda b^2 - a b \Omega_f \biggr]\biggl[ \lambda a^2 - a b \Omega_f\biggr]y^2x^2 + \biggl[ \lambda a^2 - a b \Omega_f\biggr]^2 y^4 \biggr\} </math>

 

<math>~=</math>

<math>~ 2 a b \lambda \Omega_f ( a^2 x^2 y^2 + b^2 x^4 + a^2 y^4 + b^2x^2y^2 ) - \lambda^2 (a^2 - b^2) [ a^2 y^4 - b^2 x^4 ] </math>

 

 

<math>~ + \biggl[\lambda^2 b^4 - 2a b \Omega_f \lambda b^2 + a^2 b^2 \Omega_f^2 \biggr] x^4 + 2\biggl[\lambda^2 a^2 b^2 - \lambda a b^3 \Omega_f - \lambda a^3 b\Omega_f + a^2 b^2 \Omega_f^2 \biggr]y^2x^2 + \biggl[ \lambda^2 a^4 - 2a b \Omega_f \lambda a^2 + a^2 b^2 \Omega_f^2 \biggr] y^4 </math>

 

<math>~=</math>

<math>~ x^2 y^2 \biggl\{ 2 \biggl[ \lambda^2 a^2 b^2 - \lambda a b^3 \Omega_f - \lambda a^3 b\Omega_f + a^2 b^2 \Omega_f^2 \biggr] + 2ab\lambda \Omega_f (a^2 + b^2) \biggr\} </math>

 

 

<math>~ + \biggl[\lambda^2 b^4 - 2a b \Omega_f \lambda b^2 + a^2 b^2 \Omega_f^2 + 2ab^3 \lambda \Omega_f + \lambda^2(a^2-b^2)b^2\biggr] x^4 + \biggl[ \lambda^2 a^4 - 2a b \Omega_f \lambda a^2 + a^2 b^2 \Omega_f^2 + 2ab\lambda \Omega_f a^2 - \lambda^2(a^2-b^2)a^2\biggr] y^4 </math>

 

<math>~=</math>

<math>~ 2a^2 b^2 x^2 y^2 \biggl[ \lambda^2 + \Omega_f^2 \biggr] + \biggl[ \Omega_f^2 + \lambda^2 \biggr] a^2 b^2 x^4 + \biggl[ \Omega_f^2 + \lambda^2 \biggr] a^2 b^2y^4 </math>

 

<math>~=</math>

<math>~ \biggl[ \lambda^2 + \Omega_f^2\biggr]a^2 b^2 \biggl[ x^4 +2x^2y^2 + y^4 \biggr] </math>

<math>~\Rightarrow ~~~ \frac{2\pi G \rho }{(a^2 - b^2) } \biggl[ A_1 a^2 - A_2 b^2 \biggr] </math>

<math>~=</math>

<math>~ \biggl[ \lambda^2 + \Omega_f^2\biggr] \, . </math>

[ EFE, Chapter 7, §48, Eq. (33) ]


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