User:Tohline/Appendix/Ramblings/T2Integrals

From VistrailsWiki
Jump to navigation Jump to search
Whitworth's (1981) Isothermal Free-Energy Surface
|   Tiled Menu   |   Tables of Content   |  Banner Video   |  Tohline Home Page   |

Integrals of Motion in T2 Coordinates

Motivated by the HNM82 derivation, in an accompanying chapter we have introduced a new T2 Coordinate System and have outlined a few of its properties. Here we explore whether an analytic prescription of the <math>3^\mathrm{rd}</math> integral of motion can be formulated in a potential that conforms to this set of coordinates.

Review

We begin by summarizing properties of the T2 Coordinate System that we have derived earlier. By defining the dimensionless angle,

<math> \Zeta \equiv \sinh^{-1}\biggl( \frac{qz}{\varpi} \biggr) , </math>

the two key "T2" coordinates can be written as,

<math> \chi_1 </math>

<math>\equiv</math>

<math>B \varpi \cosh\Zeta</math>

      and      

<math> \chi_2 </math>

<math>\equiv</math>

<math>\frac{A \sinh\Zeta}{ \varpi^{q^2-1}}</math>

 

<math>=</math>

<math> B ( \varpi^2 + q^2z^2 )^{1/2} = qB \xi_1 </math>

 

 

<math>=</math>

<math> \frac{Aqz}{\varpi^{q^2}} = Aq \biggl[\frac{1}{\tan\xi_2} \biggr]^{q^2} </math>

Here are a variety of relevant partial derivatives:

 

<math> \frac{\partial}{\partial x} </math>

<math> \frac{\partial}{\partial y} </math>

<math> \frac{\partial}{\partial z} </math>

<math>\chi_1</math>

<math> \biggl(\frac{B^2}{\chi_1}\biggr) x </math>

<math> \biggl(\frac{B^2}{\chi_1}\biggr) y </math>

<math> \biggl(\frac{B^2}{\chi_1}\biggr) q^2 z </math>

<math>\chi_2</math>

<math> - \biggl( \frac{q^3 A z}{\varpi^{q^2+2}} \biggr) x </math>

<math> - \biggl( \frac{q^3 A z}{\varpi^{q^2+2}} \biggr) y </math>

<math> \frac{qA}{\varpi^{q^2}} </math>

<math>\chi_3</math>

<math> - \biggl( \frac{1}{\varpi^{2}} \biggr) y </math>

<math> + \biggl( \frac{1}{\varpi^{2}} \biggr) x </math>

<math> 0 </math>

The scale factors are,

<math>h_1^2</math>

<math>=</math>

<math> \biggl[ \biggl( \frac{\partial\chi_1}{\partial x} \biggr)^2 + \biggl( \frac{\partial\chi_1}{\partial y} \biggr)^2 + \biggl( \frac{\partial\chi_1}{\partial z} \biggr)^2 \biggr]^{-1} </math>

<math>=</math>

<math> \frac{\chi_1^2}{B^4 (\varpi^2 + q^4 z^2)} </math>

<math>=</math>

<math> \frac{\chi_1^2 \ell^2}{B^4} </math>

<math>h_2^2</math>

<math>=</math>

<math> \biggl[ \biggl( \frac{\partial\chi_2}{\partial x} \biggr)^2 + \biggl( \frac{\partial\chi_2}{\partial y} \biggr)^2 + \biggl( \frac{\partial\chi_2}{\partial z} \biggr)^2 \biggr]^{-1} </math>

<math>=</math>

<math> \frac{z^2 \varpi^2 }{\chi_2^2 (\varpi^2 + q^4 z^2)} </math>

<math>=</math>

<math> \frac{z^2 \varpi^2 \ell^2}{\chi_2^2} </math>

<math>h_3^2</math>

<math>=</math>

<math> \biggl[ \biggl( \frac{\partial\chi_3}{\partial x} \biggr)^2 + \biggl( \frac{\partial\chi_3}{\partial y} \biggr)^2 + \biggl( \frac{\partial\chi_3}{\partial z} \biggr)^2 \biggr]^{-1} </math>

<math>=</math>

<math> \varpi^2 </math>

 

 

where,        <math>\ell \equiv (\varpi^2 + q^4 z^2)^{-1/2}</math>.


The position vector is,

<math>\vec{x}</math>

<math>=</math>

<math> \hat{i}x + \hat{j}y + \hat{k}z </math>

<math>=</math>

<math> \hat{e}_1 (h_1 \chi_1) + (1 - q^2) \hat{e}_2 (h_2 \chi_2) . </math>

Other Potentially Useful Differential Relations

In examining the equation of motion and searching for analytic representations of the <math>3^\mathrm{rd}</math> integral of motion, we will need to know how each of the scale factors varies along each of the coordinate directions. Since our T2 coordinate system is an orthogonal system of coordinates, in general we can write,

<math> \frac{\partial h_k}{\partial\chi_j} = \sum_{i=1}^3 h_j^2 \biggl( \frac{\partial\chi_j}{\partial x_i} \biggr) \frac{\partial h_k}{\partial x_i} , </math>

where <math>x_i</math> are the three Cartesian coordinates. Analytic expressions for the first partial derivative in each term on the RHS, <math>\partial\chi_j/\partial x_i</math>, can be obtained from the table shown above. To derive expressions for the second partial derivative in each RHS term, the following differential relations also will be useful:

 

<math> \frac{\partial}{\partial x} </math>

<math> \frac{\partial}{\partial y} </math>

<math> \frac{\partial}{\partial z} </math>

<math>\varpi</math>

<math> \frac{x}{\varpi} </math>

<math> \frac{y}{\varpi} </math>

<math> 0 </math>

<math>\ell</math>

<math> - x\ell^3 </math>

<math> - y \ell^3 </math>

<math> - q^4 z \ell^3 </math>

Putting all of these expressions together, we derive the following:


 

<math> \frac{\partial}{\partial \chi_1} </math>

<math> \frac{\partial}{\partial \chi_2} </math>

<math> \frac{\partial}{\partial \chi_3} </math>

<math>h_1</math>

<math> \frac{\ell}{B^2} - h_1^2 \ell^3 (\varpi^2 + q^6 z^2) </math>

<math> (1-q^2)q^2 h_1 h_2^2 \ell^2 \chi_2 </math>

<math> 0 </math>

<math>h_2</math>

<math> q^2 h_1^2 h_2 \ell^2 \chi_1 </math>

<math> - \frac{q^4 h_2 \ell^4 z^2}{\chi_2} \biggl[ ( \varpi^2 + q^2 z^2 ) + ( \varpi^2 + q^4 z^2 ) \biggr] </math>

<math> 0 </math>

Time Derivatives

Assuming no mistakes have been made in our derivation of the above expressions, the time-derivative of the two key scale factors become,

<math>\frac{dh_1}{dt}</math>

<math>=</math>

<math>\frac{\partial h_1}{\partial\chi_1} \dot{\chi}_1 + \frac{\partial h_1}{\partial\chi_2} \dot{\chi}_2</math>

 

<math>=</math>

<math> \biggl[ \frac{\ell}{B^2} - h_1^2 \ell^3 (\varpi^2 + q^6 z^2) \biggr]\dot{\chi}_1 + \biggl[ (1-q^2)q^2 h_1 h_2^2 \ell^2 \chi_2 \biggr]\dot{\chi}_2  ; </math>

<math>\frac{dh_2}{dt}</math>

<math>=</math>

<math>\frac{\partial h_2}{\partial\chi_1} \dot{\chi}_1 + \frac{\partial h_2}{\partial\chi_2} \dot{\chi}_2</math>

 

<math>=</math>

<math> \biggl[ q^2 h_1^2 h_2 \ell^2 \chi_1 \biggr]\dot{\chi}_1 - \frac{q^4 h_2 \ell^4 z^2}{\chi_2} \biggl[ ( \varpi^2 + q^2 z^2 ) + ( \varpi^2 + q^4 z^2 ) \biggr]\dot{\chi}_2 . </math>

If the potential is only a function of the first T2 coordinate, <math>\chi_1</math>, then the <math>2^\mathrm{nd}</math> component of the equation of motion states,

<math> \frac{d}{dt}( h_2\dot{\chi}_2 ) + h_1 \dot{\chi}_1 \biggl[ \dot{\chi}_2 \frac{1}{h_1}\frac{\partial h_2}{\partial \chi_1} - \dot{\chi}_2 \frac{1}{h_2}\frac{\partial h_1}{\partial \chi_2} \biggr] = 0 . </math>

In other words,

<math> \frac{d}{dt}( h_2\dot{\chi}_2 ) = q^2 h_1^2 h_2 \ell^2 \dot{\chi}_1 \chi_2^2\frac{d}{dt}\biggl(\frac{\chi_1}{\chi_2}\biggr) . </math>

And the time derivative of a quantity that resembles the traditional angular momentum is,

<math> \frac{d}{dt}( h_2^2 \dot{\chi}_2 ) = \frac{q^2 \ell^2 h_2^2}{\chi_2} \biggl\{ (h_1 \chi_2 \dot{\chi}_1)^2 - q^2 \ell^2 z^2 \biggl[(\varpi^2 + q^2 z^2) + (\varpi^2 + q^4 z^2) \biggr]\dot{\chi}_2^2 \biggr\} . </math>


First Special Case (quadratic)

As has been discussed in an accompanying chapter, when <math>q^2=2</math> the product of <math>\chi_1</math> and <math>\chi_2</math> shows up as a key quantity when inverting the coordinate definitions. In particular, by defining,

<math>\Chi_q \equiv \frac{2\chi_1 \chi_2}{AB} = 2\cosh \Zeta \sinh\Zeta = \sinh(2\Zeta) ,</math>

and its companion,

<math>\Upsilon_q \equiv \sqrt{1 + \Chi_q^2} = \cosh^2\Zeta + \sinh^2\Zeta = \cosh(2\Zeta) ,</math>

we can write,

<math> \varpi^2 </math>

<math> = </math>

<math> \frac{A^2}{2\chi_2^2} \biggl[ \sqrt{1 + \Chi_q^2} - 1 \biggr] </math>

<math> = </math>

<math> \biggl( \frac{A}{\chi_2} \biggr)^2 \biggl[\frac{1}{2} (\Upsilon - 1) \biggr] </math>

<math> = </math>

<math> \biggl( \frac{A}{\chi_2} \biggr)^2 \sinh^2\Zeta </math>

<math> = </math>

<math> \biggl( \frac{\chi_1}{B} \biggr)^2 \frac{1}{\cosh^2\Zeta} </math>

<math> z^2 </math>

<math> = </math>

<math> \frac{A^2}{4 \chi_2^2}\biggl[ \frac{1}{2} \Chi_q^2 +1 - \sqrt{1 + \Chi_q^2}\biggr] </math>

<math> = </math>

<math> \frac{1}{2} \biggl( \frac{A}{\chi_2} \biggr)^2 \biggl[\frac{1}{2}(\Upsilon - 1)\biggr]^2 </math>

<math> = </math>

<math> \frac{1}{2} \biggl( \frac{A}{\chi_2} \biggr)^2 \sinh^4\Zeta </math>

<math> = </math>

<math> \frac{1}{2} \biggl( \frac{\chi_1}{B} \biggr)^2 \tanh^2\Zeta </math>

<math> \ell^{-2} </math>

<math> = </math>

<math> \frac{A^2}{2\chi_2^2} \biggl[ \Chi_q^2 + 1 - \sqrt{ 1 + \Chi_q^2} \biggr] </math>

<math> = </math>

<math> \biggl( \frac{A}{\chi_2} \biggr)^2 \biggl[ \frac{1}{2} (\Upsilon - 1)\biggr] \Upsilon </math>

<math> = </math>

<math> \biggl( \frac{A}{\chi_2} \biggr)^2 \sinh^2\Zeta ( \cosh^2\Zeta + \sinh^2\Zeta ) </math>

<math> = </math>

<math> \biggl( \frac{\chi_1}{B} \biggr)^2( 1 + \tanh^2\Zeta ) </math>

Hence, potentially useful expressions for the scale factors are,

<math> 2B^2 h_2^2 </math>

<math> = </math>

<math> 2 \biggl[\frac{B z \varpi \ell}{\chi_2} \biggr]^2 </math>

<math> = </math>

<math> \biggl( \frac{\chi_1}{\chi_2} \biggr)^2 \biggl[ \frac{\tanh^2\Zeta}{1+\tanh^2\Zeta} \biggr](1-\tanh^2\Zeta) ; </math>

<math> 2B^2 h_1^2 </math>

<math> = </math>

<math> 2 \biggl[\frac{\chi_1 \ell}{B} \biggr]^2 </math>

<math> = </math>

<math> \frac{2}{1+\tanh^2\Zeta} . </math>

Notice that, because it is expressible entirely in terms of <math>\Chi_q</math>, the variable <math>\Zeta</math> is a function only of the product of the two key coordinates. Hence, the scale factor <math>h_1</math> is only a function of the product of the coordinates while <math>h_2</math> depends on the ratio, as well as the product of the two coordinates. In this context, note that <math>\Zeta</math> can be derived from the product of the two key coordinates via one of the following two relations:


<math> \frac{\chi_1\chi_2}{BA} </math>

<math> = </math>

<math> \sinh\Zeta ~(1 + \sinh^2\Zeta)^{1/2} </math>

 

<math> = </math>

<math> \cosh\Zeta ~(\cosh^2\Zeta - 1)^{1/2} . </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