Appendix/Ramblings/T3Integrals

From jetwiki
Jump to navigation Jump to search


Integrals of Motion in T3 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 offer a modest redefinition of the second radial coordinate in an effort to bring even more symmetry to the definition of the position vector, x.

Definition

By defining the dimensionless angle,

Zsinh1(qzϖ),

the two key "T3" coordinates will be written as,

λ1

ϖcoshZ=(ϖ2+q2z2)1/2

      and      

λ2

ϖ[sinhZ]1/(1q2)=[ϖq2qz]1/(q21)

Relevant Partial Derivatives

Here are some relevant partial derivatives:

 

x

y

z

λ1

xλ1

yλ1

q2zλ1

λ2

1(q21)[ϖq21sinhZ]q2/(q21)(q3zϖq2+2)x
=q2(q21)[ϖq2qz]1/(q21)(xϖ2)
=q2(q21)[ϖqz]1/(q21)(xϖ)

1(q21)[ϖq21sinhZ]q2/(q21)(q3zϖq2+2)y
=q2(q21)[ϖq2qz]1/(q21)(yϖ2)
=q2(q21)[ϖqz]1/(q21)(yϖ)

1(q21)[ϖq21sinhζ]q2/(q21)qϖq2
=1(q21)[ϖq2qz]1/(q21)1z
=1(q21)[ϖq2qzq2]1/(q21)

λ3

yϖ2

+xϖ2

0

Alternatively, partials can be taken with respect to the cylindrical coordinates, ϖ, z and ϕ. (Incidentally, I have reversed the traditional order of the ϕ and z coordinates in an attempt to parallelize structure between cylindrical and T3 coordinates since λ3ϕ.)

 

ϖ

z

ϕ

λ1

ϖλ1

q2zλ1

0

λ2

q2q21(ϖqz)1/(q21)

1q21(ϖq2qzq2)1/(q21)

0

λ3

0

0

1

Furthermore, the inverted partials are

 

λ1

λ2

λ3

ϖ

ϖ2λ1

(q21)q2ϖz22/λ2

0

z

q2z2λ1

(q21)ϖ2z2/λ2

0

ϕ

0

0

1

Scale Factors

The scale factors are,

h12

=

[(λ1x)2+(λ1y)2+(λ1z)2]1

=

λ122

 

 

h22

=

[(λ2x)2+(λ2y)2+(λ2z)2]1

=

(q21)2(ϖzλ2)2

 

 

h32

=

[(λ3x)2+(λ3y)2+(λ3z)2]1

=

ϖ2

 

 

where,        (ϖ2+q4z2)1/2.

Direction Cosines

The following table contains expressions for the nine direction cosines.

Direction Cosines for T3 Coordinates
γni=hnλnxi

 

i

n

x

y

q2z

q2xzϖ

q2yzϖ

ϖ

yϖ

+xϖ

0

where: (ϖ2+q4z2)1/2

Orthogonality Condition

Next, let's use the example orthogonality condition derived elsewhere in connection with our overview of direction cosines. Specifically, let's see if Equation DC.02 is satisfied.

λ1ϖλ2ϖ+λ1zλ2z

=

ϖλ1[q2q21(ϖqz)1/(q21)]q2zλ1[1q21(ϖq2qzq2)1/(q21)]

 

=

q2(q21)λ1{ϖ[(ϖqz)1/(q21)]z[(ϖq2qzq2)1/(q21)]}

 

=

q2(q21)λ1[q1/(q21)ϖ1+1/(q21)z1/(q21)q1/(q21)ϖq2/(q21)z1q2/(q21)]

 

=

0.

Hence, the key orthogonality condition defined by Equation DC.02 is satisfied. MF53 also gives us relationships that should apply between the various direction cosines if the coordinate system is orthogonal. Let's check a few cases to see whether γmn=Mmn, where "Mmn is the minor of γmn in the determinant |γmn|":

M11=γ22γ33γ23γ32=(ϖ)xϖ=+x.

M12=γ23γ31γ21γ33=ϖ(yϖ)=+y.

M13=γ21γ32γ22γ31=(q2xzϖ)xϖ(q2yzϖ)(yϖ)=q2z.

M31=γ12γ23γ13γ22=y(ϖ)q2z(q2yzϖ)=y2ϖ(ϖ2+q4z2)=yϖ.

M33=γ11γ22γ12γ21=x(q2yzϖ)y(q2xzϖ)=0.

All of these beautifully obey the relationship, γmn=Mmn.

Position Vector

The position vector is,

x

=

i^x+j^y+k^z

=

e^1(h1λ1)+e^2(h2λ2).

Vector Derivatives

For orthogonal coordinate systems, the time-rate-of-change of the three unit vectors are given by the expressions,

ddte^1

=

e^2A+e^3B

ddte^2

=

e^1A+e^3C

ddte^3

=

e^1Be^2C

where,

A

λ˙2h1h2λ1λ˙1h2h1λ2

B

λ˙3h1h3λ1λ˙1h3h1λ3

C

λ˙3h2h3λ2λ˙2h3h2λ3

Another way of expressing this involves Christoffel symbols and is most easily written using index notation.

ddte^c=hahc Γbcaλ˙b e^a  (ac)

Here, we have been admittedly slopping with the placement and notation of indices in order to best accommodate the notation we have been using up to here. The b index is summed over all the coordinates. The a index is summed over all coordinates EXCEPT the c coordinate. The c index is NOT summed over because it is a free index, meaning that it can equal any of the coordinates depending on which unit vector you want to differentiate.

Writing this out for each of the individual unit vectors, and striking through terms that are automatically zero when dealing with an orthogonal coordinate system, produces

ddte^1

=

h2h1(Γ112λ˙1+Γ212λ˙2+Γ312λ˙3)Ae^2+h3h1(Γ113λ˙1+Γ213λ˙2+Γ313λ˙3)Be^3

ddte^2

=

h1h2(Γ121λ˙1+Γ221λ˙2+Γ321λ˙3)Ae^1+h3h2(Γ123λ˙1+Γ223λ˙2+Γ323λ˙3)Ce^3

ddte^3

=

h1h3(Γ131λ˙1+Γ231λ˙2+Γ331λ˙3)Be^1+h2h3(Γ132λ˙1+Γ232λ˙2+Γ332λ˙3)Ce^2

where, quite generally, the 27 Christoffel symbols are,

Γ111

=

1h1h1

Γ121=Γ211

=

2h1h1

Γ131=Γ311

=

0

Γ221

=

h2h11h2h1

Γ231=Γ321

=

0

Γ331

=

h3h11h3h1

Γ112

=

h1h22h1h2

Γ122=Γ212

=

1h2h2

Γ132=Γ312

=

0

Γ222

=

2h2h2

Γ232=Γ322

=

0

Γ332

=

h3h22h3h2

Γ113

=

0

Γ123=Γ213

=

0

Γ133=Γ313

=

1h3h3

Γ223

=

0

Γ233=Γ323

=

2h3h3

Γ333

=

0

Notice that it is the Christoffel symbols labeled with all three of the coordinate indices that are automatically zero for orthogonal coordinate systems.

For additional details surrounding the Christoffel symbols (the significant role they play in the field equations, and how they can be calculated) visit this page coming soon.

Time-Derivative of Position and Velocity Vectors

In general for an orthogonal coordinate system, the velocity vector can be written as,

v=e^1(h1λ˙1)+e^2(h2λ˙2)+e^3(h3λ˙3).

So, in general, the time-rate-of-change of the velocity vector is,

dvdt

=

e^1[d(h1λ˙1)dt]+(h1λ˙1)de^1dt+e^2[d(h2λ˙2)dt]+(h2λ˙2)de^2dt+e^3[d(h3λ˙3)dt]+(h3λ˙3)de^3dt

 

=

e^1[d(h1λ˙1)dt]+(h1λ˙1)[e^2A+e^3B]+e^2[d(h2λ˙2)dt]+(h2λ˙2)[e^1A+e^3C]+e^3[d(h3λ˙3)dt]+(h3λ˙3)[e^1Be^2C]

 

=

e^1[d(h1λ˙1)dtA(h2λ˙2)B(h3λ˙3)]+e^2[d(h2λ˙2)dt+A(h1λ˙1)C(h3λ˙3)]+e^3[d(h3λ˙3)dt+B(h1λ˙1)+C(h2λ˙2)]

Now, for the T3 coordinate system the position vector has a similar form, specifically,

x=e^1(h1λ1)+e^2(h2λ2).

By analogy, then, the time-rate-of-change of the position vector is,

dxdt

=

e^1[d(h1λ1)dt]+(h1λ1)de^1dt+e^2[d(h2λ2)dt]+(h2λ2)de^2dt

 

=

e^1[d(h1λ1)dt]+(h1λ1)[e^2A+e^3B]+e^2[d(h2λ2)dt]+(h2λ2)[e^1A+e^3C]

 

=

e^1[d(h1λ1)dtA(h2λ2)]+e^2[d(h2λ2)dt+A(h1λ1)]+e^3[B(h1λ1)+C(h2λ2)]

Derived Identity for T3 Coordinates

Looking at the "e^2" component of this last expression, we have,

e^2dxdt=d(h2λ2)dt+A(h1λ1).

But we also know that,

e^2v=h2λ˙2.

Hence, it must be true that,

For T3 Coordinates

λ2dh2dt=A(h1λ1)

or, equivalently,

A=λ2h1λ1dh2dt

At some point, this identity needs to be checked by taking various partial derivatives of the scale factors and plugging them into the generic definition of A, given above. (Actually, this shouldn't be necessary because in January, 2009, we derived the same equation-of-motion result shown below while using the uglier expression for A. So this must be a correct identity in the context of T3 coordinates.)

Implications of Equation of Motion

Looking now at the "e^2" component of the acceleration (which we will set equal to zero), and assuming no motion in the 3rd component direction, we have,

e^2dvdt=d(h2λ˙2)dt+A(h1λ˙1)=0

d(h2λ˙2)dt=A(h1λ˙1)

or, inserting the relation derived above for A in terms of dh2/dt for T3 coordinates,

EOM.01

d(h2λ˙2)dt=(λ2λ˙1λ1)dh2dt

d(h2λ1λ˙2)dt=λ˙1d(h2λ2)dt.

Or, equivalently (but perhaps more perversely),

EOM.02

λ˙2λ2[1h2λ˙2d(h2λ˙2)dt]=λ˙1λ1[1h2dh2dt]

dlnλ2dt[dln(h2λ˙2)dt]=dlnλ1dt[dlnh2dt]


Note that Equation EOM.01 is equivalent to the simplest form of the conservation that we derived — actually, Jay Call derived it — back in January, 2009. Specifically,

λ1d(h2λ˙2)dt=λ2λ˙1dh2dt

The 64-thousand dollar question is, "Can we turn any of these expressions into a form which states that the total time-derivative of some function equals zero?"


Let's look at the 1st component of the equation of motion.

e^1dvdt=d(h1λ˙1)dtA(h2λ˙2)=1h1Φλ1

d(h1λ˙1)dt=A(h2λ˙2)1h1Φλ1=[h2λ2h1λ1]λ˙2dh2dt1h1Φλ1

Logarithmic Derivatives of Scale Factors

Given the collection of expressions detailed in our derivation up to this point, the logarithmic derivatives of the scale factors are:

lnh1lnλ1

=

(qh1h2λ2λ1)2

lnh1lnλ2

=

+(qh1h2λ2λ1)2

lnh2lnλ1

=

+(qh12)2

lnh2lnλ2

=

(qh12)2

This means, for example, that,

lnh1lnλ2=lnh1lnλ1,

and,

lnh2lnλ2=lnh2lnλ1.

It also means that,

dlnh2dt=(lnh2lnλ1)λ˙1λ1+(lnh2lnλ2)λ˙2λ2=(qh12)2[λ˙1λ1λ˙2λ2]

1h14dlnh2dt=dln(λ1/λ2)q2dt;

and,

dlnh1dt=(lnh1lnλ1)λ˙1λ1+(lnh1lnλ2)λ˙2λ2=(qh1h2λ2λ1)2[λ˙1λ1λ˙2λ2]

(λ1h1h2λ2)2dlnh1dt=dln(λ1/λ2)q2dt.

Hence, the logarithmic time-derivatives of the two key scale factors can be related to one another via the expression,

T3 Coordinates

(h1λ1)2dlnh1dt+(h2λ2)2dlnh2dt=0.

Also, plugging the relevant logarithmic derivatives of the key scale factors directly into the definition of A given in our discussion of vector derivatives, above, we can confirm the relationship derived earlier between A and dh2/dt. We also can relate A directly to dh1/dt. Specifically,

T3 Coordinates

A=(h2λ2h1λ1)dlnh2dt=+(h1λ1h2λ2)dlnh1dt.

It must therefore also be true that,

T3 Coordinates

(h1λ1)2lnh1λj=(h2λ2)2lnh2λj

Special Case (quadratic)

On a separate page we examine the structure of the second component of the equation of motion in the special "quadratic" case of q2=2. Jay Call is independently investigating this same special case on yet another wiki page so that we can ultimately compare and check each others' results.

See Also

Tiled Menu

Appendices: | VisTrailsEquations | VisTrailsVariables | References | Ramblings | VisTrailsImages | myphys.lsu | ADS |