Appendix/Ramblings/T9Coordinates

From JETohlineWiki
Jump to navigation Jump to search


T9 Coordinates[edit]

Establish 1st and 3rd Coordinates[edit]

λ1

≡

(x2+q2y2+p2z2)1/2,

λ3

≡

[q2yxq2]1/(q2−1)=q2/(q2−1)y1/(q2−1)xq2/(1−q2).

∂λ3∂x

=

(q21−q2)λ3x,

∂λ3∂y

=

(1q2−1)λ3y,

h3−2

=

[(q21−q2)λ3x]2+[(1q2−1)λ3y]2=λ32(q2−1)2[x2+q4y2x2y2]

⇒h3

=

(q2−1)λ3[xy(x2+q4y2)1/2].

Hence, the three direction-cosines are,

γ31=h3(∂λ3∂x)

=

(q2−1)λ3[xy(x2+q4y2)1/2](q21−q2)λ3x=−[q2y(x2+q4y2)1/2],

γ32=h3(∂λ3∂y)

=

(q2−1)λ3[xy(x2+q4y2)1/2](1q2−1)λ3y=[x(x2+q4y2)1/2].

And the position vector is given by the expression,

x→

=

e^1(γ11x+γ12y+γ13z)+e^2(γ21x+γ22y+γ23z)+e^3(γ31x+γ32y+γ33z)

 

=

e^1[x2ℓ3D+q2y2ℓ3D+p2z2ℓ3D]+e^2(γ21x+γ22y+γ23z)+e^3[−q2xy(x2+q4y2)1/2+xy(x2+q4y2)1/2]

 

=

e^1[λ12ℓ3D]+e^2(γ21x+γ22y+γ23z)−e^3[(q2−1)xy(x2+q4y2)1/2].

Guess 2nd Coordinate[edit]

Unspecified Coefficients[edit]

λ2

≡

(ax2+by2+ez2)−1/2(ax2+by2)1/2

 

=

[(ax2+by2)(ax2+by2+ez2)]1/2

 

=

[1+ez2(ax2+by2)]−1/2.

∂λ2∂x

=

−12[1+ez2(ax2+by2)]−3/2∂∂x[ez2(ax2+by2)−1]

 

=

−ez22[1+ez2(ax2+by2)]−3/2(−1)[(ax2+by2)−2]2ax

 

=

axez2(ax2+by2)2[1+ez2(ax2+by2)]−3/2

 

=

axez2λ23(ax2+by2)2,

∂λ2∂y

=

−12[1+ez2(ax2+by2)]−3/2∂∂y[ez2(ax2+by2)−1]

 

=

−ez22[1+ez2(ax2+by2)]−3/2(−1)[(ax2+by2)−2]2by

 

=

byez2(ax2+by2)2[1+ez2(ax2+by2)]−3/2

 

=

byez2λ23(ax2+by2)2,

∂λ2∂z

=

−12[1+ez2(ax2+by2)]−3/2∂∂z[ez2(ax2+by2)−1]

 

=

−ezλ23(ax2+by2).

Hence,

h2−2

=

[axez2λ23(ax2+by2)2]2+[byez2λ23(ax2+by2)2]2+[ezλ23(ax2+by2)]2

 

=

[ezλ23(ax2+by2)2]2[(axz)2+(byz)2+(ax2+by2)2]

⇒h2

=

[(ax2+by2)2ezλ23][(axz)2+(byz)2+(ax2+by2)2]−1/2.

The three direction-cosines are, then,

γ21=h2(∂λ2∂x)

=

axez2λ23(ax2+by2)2[(ax2+by2)2ezλ23][(axz)2+(byz)2+(ax2+by2)2]−1/2

 

=

axz[(axz)2+(byz)2+(ax2+by2)2]−1/2

γ22=h2(∂λ2∂y)

=

byez2λ23(ax2+by2)2[(ax2+by2)2ezλ23][(axz)2+(byz)2+(ax2+by2)2]−1/2

 

=

byz[(axz)2+(byz)2+(ax2+by2)2]−1/2

γ23=h2(∂λ2∂z)

=

−ezλ23(ax2+by2)[(ax2+by2)2ezλ23][(axz)2+(byz)2+(ax2+by2)2]−1/2

 

=

−(ax2+by2)[(axz)2+(byz)2+(ax2+by2)2]−1/2.

Direction Cosine Components for T9 Coordinates
n λn hn ∂λn∂x ∂λn∂y ∂λn∂z γn1 γn2 γn3
(1) (2) (3) (4) (5) (6) (7) (8) (9)
1 (x2+q2y2+p2z2)1/2 λ1ℓ3D xλ1 q2yλ1 p2zλ1 (x)ℓ3D (q2y)ℓ3D (p2z)ℓ3D
2A --- ℓ3D(x2+q4y2)1/2 xp2z(x2+q4y2) q2yp2z(x2+q4y2) −1 xp2zℓ3D(x2+q4y2)1/2 q2yp2zℓ3D(x2+q4y2)1/2 −(x2+q4y2)ℓ3D(x2+q4y2)1/2
2B [1+ez2(ax2+by2)]−1/2 [(ax2+by2)2ezλ23]𝔍2B axez2λ23(ax2+by2)2 byez2λ23(ax2+by2)2 −ezλ23(ax2+by2) axz𝔍2B byz𝔍2B −(ax2+by2)𝔍2B
3 [q2yxq2]1/(q2−1) (q2−1)xyλ3(x2+q4y2)1/2 −(q2q2−1)λ3x (q2q2−1)λ3q2y 0 −q2y(x2+q4y2)1/2 x(x2+q4y2)1/2 0

ℓ3D

≡

(x2+q4y2+p4z2)−1/2,

𝔍0−2≡(x2+q4y2)ℓ3D2

=

[(xp2z)2+(q2yp2z)2+(x2+q4y2)2],

𝔍2B

≡

[(axz)2+(byz)2+(ax2+by2)2]−1/2.


This table titled, "Direction Cosine Components for T9 Coordinates," contains the following information:

  1. The (first) row labeled, n=1, correctly details the scale-factor, h1, and the unit vector expression, e^1=(ı^γ11+ȷ^γ12+k^γ13), that result from the given specification of the λ1 coordinate. By design, the unit vector, e^1, is everywhere normal to the "surface" of the ellipsoid.
  2. The (fourth) row labeled, n=3, correctly details the scale-factor, h3, and the unit vector expression, e^3=(ı^γ31+ȷ^γ32+k^γ33), that result from the given specification of the λ2 coordinate. By design, this unit vector, e^3, has no vertical component — that is, γ33=0 — and, by design, it is everywhere perpendicular to the "surface-normal" unit vector, e^1.
  3. We desire a unit vector, e^2, that is mutually orthogonal to the other two unit vectors; this has been accomplished by examining their cross-product, namely, we have set e^2A=e^3×e^1. Determined in this manner, the expressions for the three direction-cosine components of e^2A have been written in the last three columns of the (second) row labeled, n=2A. While we are confident that the correct specification of e^2 is,

    e^2A=(ı^γ21+ȷ^γ22+k^γ23)

    =

    ı^(xp2z)𝔍0+ȷ^(q2yp2z)𝔍0−k^(x2+q4y2)𝔍0,

    as yet (18 February 2021), we have been unable to determine an expression for the coordinate, λ2A(x,y,z), from which all three of these direction-cosine expressions can be simultaneously derived — hence, the dashes in the second column of row 2A. The expression for h2A that has been presented in the third column of row 2A (and framed in pink) is a guess which, when divided into each of the three direction cosines, gives respectively the three guessed partial-derivative expressions shown in columns 4, 5, and 6 of row 2A.

  4. The second column of row 2B contains a guess (framed in yellow) for the coordinate expression, λ2B; this expression contains three unspecified scalar coefficients, a, b and e. The remaining columns of this row contain the three partial derivatives, the associated scale factor, and the three direction cosines that result from this guessed coordinate expression. If we can find values of the three scalar coefficients that give (row 2B) expressions for the three direction cosines that perfectly match the direction cosines written in row 2A, then we will be able to state that λ2B is — at least one form of — our sought-after third coordinate expression.

Yellow-Framed Guess 2B[edit]

Referencing the 2B table row, above, we are looking for coefficient values that map,

axz→xp2z,

     

byz→q2yp2z,

     

(ax2+by2)→(x2+q4y2),

and that map,

𝔍2B≡[(axz)2+(byz)2+(ax2+by2)2]−1/2

=

[z2(a2x2+b2y2)+(ax2+by2)2]−1/2,

into the expression,

𝔍0≡ℓ3D(x2+q4y2)1/2

=

[(x2+q4y2+p4z2)(x2+q4y2)]−1/2

 

=

[p4z2(x2+q4y2)+(x2+q4y2)2]−1/2

 

=

[x2p4z2+q4y2p4z2+(x2+q4y2)2]−1/2.

A portion of these mappings are accomplished by setting a=p2 and b=q2p2, but this pair of specified coefficient values does not satisfy other mappings. Alternatively, a separate subset of mappings — but, again, not all mappings — is satisfied by setting a=1 and b=q4. So the yellow-framed guess 2B does not provide a correct second-coordinate expression.

Blue-Framed Guess 2C[edit]

Let's try,

λ2B

=

(cz)2e(ax2+by2)−f.


Direction Cosine Components for Additional T9 Coordinate Guesses
n
 
(1)
λn
 
(2)
hn
 
(3)
∂λn/∂x
 
(4)
∂λn/∂y
 
(5)
∂λn/∂z
 
(6)
γn1
 
(7)
γn2
 
(8)
γn3
 
(9)
2A --- --- --- --- --- xp2z𝔍0 q2yp2z𝔍0 −(x2+q4y2)𝔍0
2B [1+ez2(ax2+by2)]−1/2 [(ax2+by2)2ezλ23]𝔍2B axez2λ23(ax2+by2)2 byez2λ23(ax2+by2)2 −ezλ23(ax2+by2) axz𝔍2B byz𝔍2B −(ax2+by2)𝔍2B
2C (cz)2e(ax2+by2)−f [z(ax2+by2)2eλ2C]𝔍2C −2afxλ2C(ax2+by2) −2bfyλ2C(ax2+by2) 2eλ2Cz −(afxze)𝔍2C −(bfyze)𝔍2C (ax2+by2)𝔍2C

ℓ3D

≡

(x2+q4y2+p4z2)−1/2,

𝔍0−2≡(x2+q4y2)ℓ3D2

=

[(xp2z)2+(q2yp2z)2+(x2+q4y2)2],

𝔍2B

≡

[(axz)2+(byz)2+(ax2+by2)2]−1/2,

𝔍2C

≡

[(xz)2(fae)2+(yz)2(fbe)2+(ax2+by2)2]−1/2.


Examining just the expression for 𝔍2C, we see that we definitely need: a=1 and b=q4. Also, we need,

afe

=

p2

    ⇒   

fe

=

p2;

and, separately we need,

fbe

=

q2p2

    ⇒   

fe

=

p2q2.

These cannot simultaneously be satisfied. So the blue-framed guess 2C does not provide a correct second-coordinate expression. But we are very close! We need one additional scalar coefficient degree of freedom.

Next Thought[edit]

More generally, this leads to an expression for the scale-factor of the form,

h2−2

∼

𝔊(x,y,z)[x2p4z2+q4y2p4z2+(x2+q4y2)2]

⇒∂λ2∂x

∼

[𝔊(x,y,z)]1/2[xp2z];

     

∂λ2∂y

∼

[𝔊(x,y,z)]1/2[q2yp2z];

     

∂λ2∂z

∼

[𝔊(x,y,z)]1/2[x2+q4y2].



Now, if we set, a=p2 and b=q2p2, we have,

(axz)2+(byz)2+(ax2+by2)2

=

(xp2z)2+(q2p2yz)2+(p2x2+q2p2y2)2

 

=

p4[q4y2z2+x2z2+(x2+q2y2)2],

in which case,

Complete Orthogonality Check[edit]

If the set of unit vectors is indeed orthogonal, then we must find that,

γmn

=

Mmn,

where the quantity Mmn is the minor of γmn in the determinant, |γmn|. (Note: This last expression is true only for right-handed coordinate systems. If the coordinate system is left-handed, we should find, γmn=−Mmn.) More specifically, for any right-handed, orthogonal curvilinear coordinate system we should find:

γ11

=

γ22γ33−γ23γ32,

γ12

=

γ23γ31−γ21γ33,

γ13

=

γ21γ32−γ22γ31,

   

γ21

=

γ32γ13−γ33γ12,

γ22

=

γ33γ11−γ31γ13,

γ23

=

γ31γ12−γ32γ11,

   

γ31

=

γ12γ23−γ13γ22,

γ32

=

γ13γ21−γ11γ23,

γ33

=

γ11γ22−γ12γ21.

and the position vector is,

x→

=

ı^x+ȷ^y+k^z

 

=

e^1(γ11x+γ12y+γ13z)+e^2(γ21x+γ22y+γ23z)+e^3(γ31x+γ32y+γ33z).

For (κ1, κ2, κ3)[edit]

γ22γ33−γ23γ32

=

[q2yp2zℓ3D(x2+q4y2)1/2][0]+[(x2+q4y2)ℓ3D(x2+q4y2)1/2][x(x2+q4y2)1/2]=xℓ3D=γ11.       Yes!

γ23γ31−γ21γ33

=

[(x2+q4y2)1/2ℓ3D][q2y(x2+q4y2)1/2]−[xp2zℓ3D(x2+q4y2)1/2][0]=q2yℓ3D=γ12.       Yes!

γ21γ32−γ22γ31

=

[xp2zℓ3D(x2+q4y2)1/2][x(x2+q4y2)1/2]+[q2yp2zℓ3D(x2+q4y2)1/2][q2y(x2+q4y2)1/2]

 

=

[x2p2zℓ3D(x2+q4y2)]+[q4y2p2zℓ3D(x2+q4y2)]=p2zℓ3D=γ13.       Yes!

γ32γ13−γ33γ12

=

[x(x2+q4y2)1/2][p2zℓ3D]−[0][q2yℓ3D]=[xp2zℓ3D(x2+q4y2)1/2]=γ21.       Yes!

γ33γ11−γ31γ13

=

[0][xℓ3D]+[q2y(x2+q4y2)1/2][p2zℓ3D]=[q2yp2zℓ3D(x2+q4y2)1/2]=γ22.       Yes!

γ31γ12−γ32γ11

=

−[q2y(x2+q4y2)1/2][q2yℓ3D]−[x(x2+q4y2)1/2][xℓ3D]=−(x2+q4y2)1/2ℓ3D=γ23.       Yes!

γ12γ23−γ13γ22

=

−[q2yℓ3D][(x2+q4y2)1/2ℓ3D]−[p2zℓ3D][q2yp2zℓ3D(x2+q4y2)1/2]

 

=

−q2yℓ3D2(x2+q4y2)1/2[x2+q4y2+p4z2]=−q2y(x2+q4y2)1/2=γ31.       Yes!

γ13γ21−γ11γ23

=

[p2zℓ3D][xp2zℓ3D(x2+q4y2)1/2]+[xℓ3D][(x2+q4y2)1/2ℓ3D]

 

=

xℓ3D2(x2+q4y2)1/2[x2+q4y2+p4z2]=x(x2+q4y2)1/2=γ32.       Yes!

γ11γ22−γ12γ21

=

[xℓ3D][q2yp2zℓ3D(x2+q4y2)1/2]−[q2yℓ3D][xp2zℓ3D(x2+q4y2)1/2]=0=γ33.       Yes!

Given that the prescribed interrelationships between all nine direction cosines are satisfied, we conclude that the (κ1,κ2,κ3) coordinate system is an orthogonal one. Accordingly, the position vector is,

x→

=

e^1{x2ℓ3D+q2y2ℓ3D+p2z2ℓ3D}

 

 

+e^2{x2+q2y2−1p2(x2+q4y2)}p2zℓ3D(x2+q4y2)1/2

 

 

+e^3{−xq2y(x2+q4y2)1/2+xy(x2+q4y2)1/2}

 

=

e^1[κ12ℓ3D]+e^2[x2+q2y2−1p2(x2+q4y2)]p2zℓ3D(x2+q4y2)1/2+e^3[(1−q2)xy(x2+q4y2)1/2]

 

=

e^1[κ12ℓ3D]+e^2[x2(1−1p2)−q2y2(q2−1)]p2zℓ3D(x2+q4y2)1/2−e^3[(q2−1)xy(x2+q4y2)1/2].

For (κ1, κ4, κ5)[edit]

γ42γ53−γ43γ52

=

[xq2yp2z𝒟(1q2y)][p2zℓ3D𝒟(x2−q4y2)]+[xq2yp2z𝒟(2p2z)][q2yℓ3D𝒟(p4z2+2x2)]

 

=

xℓ3D𝒟2[(x2−q4y2)p4z2+(p4z2+2x2)2q4y2]

 

=

xℓ3D𝒟2[x2p4z2+4x2q4y2+q4y2p4z2]=xℓ3D=γ11.       Yes!

γ43γ51−γ41γ53

=

[xq2yp2z𝒟(2p2z)][xℓ3D𝒟(2q4y2+p4z2)]−[xq2yp2z𝒟(1x)][p2zℓ3D𝒟(x2−q4y2)]

 

=

q2yℓ3D𝒟2[2x2(2q4y2+p4z2)−p4z2(x2−q4y2)]

 

=

q2yℓ3D𝒟2[4x2q4y2+x2p4z2+q4y2p4z2]=q2yℓ3D=γ12.       Yes!

γ41γ52−γ42γ51

=

[xq2yp2z𝒟(1x)][q2yℓ3D𝒟(p4z2+2x2)]+[xq2yp2z𝒟(1q2y)][xℓ3D𝒟(2q4y2+p4z2)]

 

=

p2zℓ3D𝒟2[q4y2(p4z2+2x2)+x2(2q4y2+p4z2)]

 

=

p2zℓ3D𝒟2[q4y2p4z2+4x2q4y2+x2p4z2]=p2zℓ3D=γ13.       Yes!

γ52γ13−γ53γ12

=

[q2yp2zℓ3D2𝒟(p4z2+2x2)]−[q2yp2zℓ3D2𝒟(x2−q4y2)]

 

=

q2yp2zℓ3D2𝒟[x2+q4y2+p4z2]=q2yp2z𝒟=γ41.       Yes!

γ53γ11−γ51γ13

=

[xp2zℓ3D2𝒟(x2−q4y2)]+[xp2zℓ3D2𝒟(2q4y2+p4z2)]

 

=

xp2zℓ3D2𝒟[x2+q4y2+p4z2]=xp2z𝒟=γ42.       Yes!

γ51γ12−γ52γ11

=

−[xq2yℓ3D2𝒟(2q4y2+p4z2)]−[xq2yℓ3D2𝒟(p4z2+2x2)]

 

=

−2xq2yℓ3D2𝒟[x2+q4y2+p4z2]=−2xq2y𝒟=γ43.       Yes!

γ12γ43−γ13γ42

=

−[2xq4y2ℓ3D𝒟]−[xp4z2ℓ3D𝒟]=−[xℓ3D𝒟](2q4y2+p4z2)=γ51.       Yes!

γ13γ41−γ11γ43

=

[q2yp4z2ℓ3D𝒟]+[2x2q2yℓ3D𝒟]=[q2yℓ3D𝒟](p4z2+2x2)=γ52.       Yes!

γ11γ42−γ12γ41

=

[x2p2zℓ3D𝒟]−[q4y2p2zℓ3D𝒟]=[p2zℓ3D𝒟](x2−q4y2)=γ53.       Yes!

Given that the prescribed interrelationships between all nine direction cosines are satisfied, we conclude that the (κ1,κ4,κ5) coordinate system is an orthogonal one. Accordingly, the position vector is,

x→

=

e^1(γ11x+γ12y+γ13z)+e^2(γ41x+γ42y+γ43z)+e^3(γ51x+γ52y+γ53z)

 

=

e^1{x2ℓ3D+q2y2ℓ3D+p2z2ℓ3D}

 

 

+e^2{xq2yp2z𝒟+xyp2z𝒟−2xq2yz𝒟}

 

 

+e^3{−ℓ3Dx2𝒟(2q4y2+p4z2)+ℓ3Dq2y2𝒟(p4z2+2x2)+ℓ3Dp2z2𝒟(x2−q4y2)}

 

=

e^1{κ12ℓ3D}+e^2{q2p2+p2−2q2}xyz𝒟

 

 

+e^3{−x2(2q4y2+p4z2)+q2y2(p4z2+2x2)+p2z2(x2−q4y2)}ℓ3D𝒟

 

=

e^1{κ12ℓ3D}+e^2{q2p2+p2−2q2}xyz𝒟−e^3{2x2q2y2(q2−1)+x2p2z2(p2−1)+q2y2p2z2(q2−p2)}ℓ3D𝒟.

See Also[edit]


Tiled Menu

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