Difference between revisions of "User:Jaycall/KillingVectorApproach"

From VistrailsWiki
Jump to navigation Jump to search
(Move subsection "application to T3 coordinates" to entirely new page)
(→‎Strategic Choice of Characteristic Vector: Add a couple of right-justified equation numbers)
Line 108: Line 108:
So we're looking for a vector <math>\vec{C}</math> such that
So we're looking for a vector <math>\vec{C}</math> such that


<span id="CV.01"><table align="right" border="1" cellpadding="10" width="10%">
<tr><th><font color="darkblue">CV.01</font></th></tr>
</table></span>
<div align="center">
<div align="center">
<math>
<math>
Line 246: Line 249:
And limiting our interest to motion within the meridional plane (setting <math>\dot{\lambda}_3 = 0</math>) and simplifying
And limiting our interest to motion within the meridional plane (setting <math>\dot{\lambda}_3 = 0</math>) and simplifying


<span id="CV.02"><table align="right" border="1" cellpadding="10" width="10%">
<tr><th><font color="darkblue">CV.02</font></th></tr>
</table></span>
<table align="center" border="1" cellpadding="5">
<table align="center" border="1" cellpadding="5">
<tr>
<tr>
   <td align="right">
   <td align="right">
<math>
<math>
\frac{d\Xi}{dt}  
\frac{d\Xi}{dt} \equiv \frac{d\ln C_2}{dt}
</math>
</math>
   </td>
   </td>

Revision as of 21:25, 6 June 2010

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

Killing Vector Approach to Integrals of Motion Problem

Motivated by our recent work in relativistic fluid dynamics, one insightful approach to the integrals of motion problem involves rewriting the equations of motion in terms of an arbitrary vector field, which we're currently calling a characteristic vector. The benefit here is that we can produce a mathematical condition that tells us what combintations of variables should be grouped together in order to produce conserved quantities. The details follow.

Newtonian Equations of Motion

Using index summation notation (wherein repeated indices are summed over), it is possible to write the Newtonian equations of motion for a test particle in a fixed, static potential in an arbitrary orthogonal coordinate system. The equation describing motion in the <math>i</math>-dimension is:

<math> \frac{d}{dt} \left( m \ {h_i}^2 \dot{\lambda}_i \right) = m \ {h_k}^2 \Gamma^k_{ij} \dot{\lambda}_j \dot{\lambda}_k - m \ \partial_i \Phi </math>

where <math>m</math> is the mass of the test particle, the <math>\lambda</math>s are the chosen coordinates, and the <math>\dot{\lambda}</math>s are the coordinate velocities (or the total time-derivatives of the coordinates), the <math>h</math>s are the metric scale factors associated with each of the coordinates, <math>\Phi</math> is the gravitational potential, and the <math>\Gamma</math>s are the Christoffel symbols.

Christoffel Symbols

For an orthogonal coordinate system, the Christoffel symbols can be written (depending on which, if any, indices are repeated) as

Christoffel Symbols
(for orthogonal coordinate systems)

<math> \Gamma^k_{ij} = 0 </math>

for

<math> i \ne j \ne k </math>

<math> \Gamma^k_{ii} = - \frac{h_i}{h_k} \frac{\partial_k h_i}{h_k} </math>

for

<math> i \ne k </math>

<math> \Gamma^i_{ij} = \Gamma^i_{ji} = \frac{\partial_i h_j}{h_j} </math>

Shouldn't this be:   <math>\frac{\partial_j h_i}{h_i}</math> ?

<math> \Gamma^i_{ii} = \frac{\partial_i h_i}{h_i} </math>

otherwise

Identification of Conserved Quantity

If, through some miracle, the right-hand side (RHS) of the equation of motion should equal zero, then the combination of variables inside the parentheses must by definition be a conserved quantity. But good luck guessing a coordinate system in which this miracle will occur!

Introduce Characteristic Vector

The trick is to introduce an arbitrary vector field, which will be used to build a location-dependent weighted linear combination of the equations of motion. Since the equations of motion associated with each of the coordinates collectively form a single vector equation, we will just form an inner product of our characteristic vector with each side of the vector equation of motion. The result is

<math> C_i \frac{d}{dt} \left( m \ {h_i}^2 \dot{\lambda}_i \right) = m \ {h_k}^2 \Gamma^k_{ij} \dot{\lambda}_j \dot{\lambda}_k C_i - m \ C_i \ \partial_i \Phi </math>

Next, we bring <math>C_i</math> inside the total time-derivative on the left-hand side (LHS). This produces an additional term, which we promptly move over to the RHS and include as part of what is commonly referred to as the source (since it's the source of any change in the quantity in parentheses).

<math> \frac{d}{dt} \left( m \ {h_i}^2 \dot{\lambda}_i C_i \right) = m \ {h_i}^2 \dot{\lambda}_i \dot{C}_i + m \ {h_k}^2 \Gamma^k_{ij} \dot{\lambda}_j \dot{\lambda}_k C_i - m \ C_i \ \partial_i \Phi </math>

By doing this, we have formed a new conservative quantity; that is, a new quantity in the parentheses which will be conserved if and only if the source is zero. The utility is in the fact that we should be able to force the source associated with this new conservative quantity to go to zero by choosing the right characteristic vector. Once we find the right characteristic vector, we'll be able to use it directly to build the conserved quantity <math>m \ {h_i}^2 \dot{\lambda}_i C_i</math>.

Strategic Choice of Characteristic Vector

So we're looking for a vector <math>\vec{C}</math> such that

CV.01

<math> {h_i}^2 \dot{\lambda}_i \dot{C}_i + {h_k}^2 \Gamma^k_{ij} \dot{\lambda}_j \dot{\lambda}_k C_i - C_i \ \partial_i \Phi = 0 . </math>

The good news is that this is a single scalar equation constraining the three independent components of <math>\vec{C}</math>, so several families of solutions should exist. Consequently, we can afford to be choosey!

In the problem we're exploring using T3 coordinates, we already know a conserved quantity associated with the azimuthal coordinate — angular momentum. We're looking for an additional (independent) conserved quantity associated with the two radial coordinates. So it makes sense to look for one of the solutions that has <math>C_3 = 0</math>. Furthermore, we'd like to avoid dealing with the unknown potential function (which varies only with <math>\lambda_1</math>) as much as possible, so as long as we're being choosey, let's look for a solution that has <math>C_1 = 0</math>. We now have three conditions on three components. The condition constraining <math>C_2</math> is

<math> {h_2}^2 \dot{\lambda}_2 \dot{C}_2 + {h_k}^2 \Gamma^k_{2j} \dot{\lambda}_j \dot{\lambda}_k C_2 = 0 . </math>

This can be rewritten more simply as

<math> \dot{C}_2 = - \left( \frac{{h_k}^2}{{h_2}^2} \Gamma^k_{2j} \frac{\dot{\lambda}_j \dot{\lambda}_k}{\dot{\lambda}_2} \right) C_2 . </math>

Perhaps the best approach to solving this condition is separation of variables.

 

<math> \int \frac{dC_2}{C_2} </math>

<math> = </math>

<math> - \int \left( \frac{{h_k}^2}{{h_2}^2} \Gamma^k_{2j} \frac{\dot{\lambda}_j \dot{\lambda}_k}{\dot{\lambda}_2} \right) dt </math>

<math> \Rightarrow </math>

<math> \ln C_2 </math>

<math> = </math>

<math> - \int \left( \frac{{h_k}^2}{{h_2}^2} \Gamma^k_{2j} \frac{\dot{\lambda}_j \dot{\lambda}_k}{\dot{\lambda}_2} \right) dt </math>

<math> \Rightarrow </math>

<math> C_2 </math>

<math> = </math>

<math> \exp \left\{ - \int \left( \frac{{h_k}^2}{{h_2}^2} \Gamma^k_{2j} \frac{\dot{\lambda}_j \dot{\lambda}_k}{\dot{\lambda}_2} \right) dt \right\} . </math>

The final necessary step will be to write the quantity in parentheses as the exact derivative of some other quantity, call it <math>\Xi</math>. The long-sought conserved quantity will then be <math> m {h_2}^2 \dot{\lambda}_2 \exp \left( - \Xi \right) </math> , where

<math> \Xi \equiv \int \left( \frac{{h_k}^2}{{h_2}^2} \Gamma^k_{2j} \frac{\dot{\lambda}_j \dot{\lambda}_k}{\dot{\lambda}_2} \right) dt . </math>


Question from Joel

Is the following logic correct?

Suppose we examine the term in which <math>j=2</math> and <math>k=1</math>. The integral becomes,

<math> \Xi\biggr|_{j=2,k=1} = \int \left( \frac{{h_1}^2}{{h_2}^2} \Gamma^1_{22} \frac{\dot{\lambda}_2 \dot{\lambda}_1}{\dot{\lambda}_2} \right) dt = \int \left( \frac{{h_1}^2}{{h_2}^2} \Gamma^1_{22} \right) d\lambda_1 </math>

Given that the scale factors and the Christoffel symbols are expressible entirely in terms of <math>\lambda_1</math> and <math>\lambda_2</math>, one could imagine a situation — for example in the quadratic case of <math>q^2=2</math> — in which this integral could be completed analytically. (I presume that wherever <math>\lambda_2</math> appears inside this integral, it can be treated as a constant because the two coordinates are independent of one another.)

Response from Jay

Yes, I believe you've worked this term out correctly, but as you can see in what follows, it cancels out with the other cross term. Writing out each of the terms in <math>\Xi</math> gives,

<math> \frac{d\Xi}{dt} = \frac{{h_1}^2}{{h_2}^2} \Gamma^1_{21} \frac{\dot{\lambda}_1 \dot{\lambda}_1}{\dot{\lambda}_2} + \frac{{h_1}^2}{{h_2}^2} \Gamma^1_{22} \frac{\dot{\lambda}_2 \dot{\lambda}_1}{\dot{\lambda}_2} + \frac{{h_1}^2}{{h_2}^2} \Gamma^1_{23} \frac{\dot{\lambda}_3 \dot{\lambda}_1}{\dot{\lambda}_2} + \frac{{h_2}^2}{{h_2}^2} \Gamma^2_{21} \frac{\dot{\lambda}_1 \dot{\lambda}_2}{\dot{\lambda}_2} + \frac{{h_2}^2}{{h_2}^2} \Gamma^2_{22} \frac{\dot{\lambda}_2 \dot{\lambda}_2}{\dot{\lambda}_2} + \frac{{h_2}^2}{{h_2}^2} \Gamma^2_{23} \frac{\dot{\lambda}_3 \dot{\lambda}_2}{\dot{\lambda}_2} + \frac{{h_3}^2}{{h_2}^2} \Gamma^3_{21} \frac{\dot{\lambda}_1 \dot{\lambda}_3}{\dot{\lambda}_2} + \frac{{h_3}^2}{{h_2}^2} \Gamma^3_{22} \frac{\dot{\lambda}_2 \dot{\lambda}_3}{\dot{\lambda}_2} + \frac{{h_3}^2}{{h_2}^2} \Gamma^3_{23} \frac{\dot{\lambda}_3 \dot{\lambda}_3}{\dot{\lambda}_2} . </math>

Plugging in values for the Christoffel symbols leads to the expression

<math> \frac{d\Xi}{dt} = \frac{{h_1}^2}{{h_2}^2} \frac{\partial_2 h_1}{h_1} \frac{\dot{\lambda}_1 \dot{\lambda}_1}{\dot{\lambda}_2} - \frac{{h_1}^2}{{h_2}^2} \frac{h_2}{h_1} \frac{\partial_1 h_2}{h_1} \frac{\dot{\lambda}_2 \dot{\lambda}_1}{\dot{\lambda}_2} + \frac{{h_2}^2}{{h_2}^2} \frac{\partial_1 h_2}{h_2} \frac{\dot{\lambda}_1 \dot{\lambda}_2}{\dot{\lambda}_2} + \frac{{h_2}^2}{{h_2}^2} \frac{\partial_2 h_2}{h_2} \frac{\dot{\lambda}_2 \dot{\lambda}_2}{\dot{\lambda}_2} + \frac{{h_3}^2}{{h_2}^2} \frac{\partial_2 h_3}{h_3} \frac{\dot{\lambda}_3 \dot{\lambda}_3}{\dot{\lambda}_2} . </math>

And limiting our interest to motion within the meridional plane (setting <math>\dot{\lambda}_3 = 0</math>) and simplifying

CV.02

<math> \frac{d\Xi}{dt} \equiv \frac{d\ln C_2}{dt} </math>

<math> = </math>

<math> \frac{h_1}{h_2} \frac{\partial_2 h_1}{h_2} \frac{{\dot{\lambda}_1}^2}{\dot{\lambda}_2} - \cancel{\frac{\partial_1 h_2}{h_2} \dot{\lambda}_1} + \cancel{\frac{\partial_1 h_2}{h_2} \dot{\lambda}_1} + \frac{\partial_2 h_2}{h_2} \dot{\lambda}_2 </math>

 

<math> = </math>

<math> \biggl[ \biggl(\frac{h_1 \dot{\lambda}_1}{h_2 \dot{\lambda}_2}\biggr)^2 \frac{\partial \ln h_1}{\partial\ln\lambda_2} + \frac{\partial \ln h_2}{\partial \ln\lambda_2}\biggr] \frac{d \ln{\lambda}_2}{dt} </math>

Application to T3 Coordinates

On a separate page, we have applied this "characteristic vector" approach specifically to T3 coordinates.

 

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