Appendix/Ramblings/T6CoordinatesPt2

From jetwiki
Revision as of 22:13, 23 July 2021 by Jet53man (talk | contribs) (→Speculation7)
(diff) ← Older revision | Latest revision (diff) | Newer revision → (diff)
Jump to navigation Jump to search

Concentric Ellipsoidal (T6) Coordinates (Part 2)

Orthogonal Coordinates

Speculation5

Spherical Coordinates

rcos⁡θ

=

z,

rsin⁡θ

=

(x2+y2)1/2,

tan⁡φ

=

yx.

Use λ1 Instead of r

Here, as above, we define,

λ12

≡

x2+q2y2+p2z2

Using this expression to eliminate "x" (in favor of λ1) in each of the three spherical-coordinate definitions, we obtain,

r2≡x2+y2+z2

=

λ12−y2(q2−1)−z2(p2−1);

tan2θ≡x2+y2z2

=

1z2[λ12−y2(q2−1)−p2z2];

1tan2φ≡x2y2

=

λ12−p2z2y2−q2.

After a bit of additional algebraic manipulation, we find that,

z2λ12

=

(1+tan2φ)𝒟2,

y2λ12

=

[𝒟2tan2φ−p2tan2φ(1+tan2φ)(1+q2tan2φ)𝒟2],

x2λ12

=

1−q2(y2λ12)−p2(z2λ12),

where,

𝒟2

≡

[(1+q2tan2φ)(p2+tan2θ)−p2(q2−1)tan2φ].

As a check, let's set q2=p2=1, which should reduce to the normal spherical coordinate system.

λ12

→

r2,

      and,      

𝒟2

→

[(1+tan2φ)(1+tan2θ)].

⇒z2λ12

→

11+tan2θ=cos2θ=z2r2;

y2λ12

→

[(1+tan2φ)(1+tan2θ)tan2φ−tan2φ(1+tan2φ)(1+tan2φ)(1+tan2φ)(1+tan2θ)]

 

=

tan2φ(1+tan2φ)[tan2θ(1+tan2θ)]=sin2θsin2φ=y2r2;

x2λ12

→

1−(y2λ12)−(z2λ12),

 

→

1−sin2θsin2φ−cos2θ=−sin2θsin2φ+sin2θ=sin2θcos2φ=x2r2.

Relationship To T3 Coordinates

If we set, q=1, but continue to assume that p>1, we expect to see a representation that resembles our previously discussed, T3 Coordinates. Note, for example, that the new "radial" coordinate is,

λ12

→

(ϖ2+p2z2),

      and,      

𝒟2

→

[(1+tan2φ)(p2+tan2θ)].

⇒p2z2λ12

→

p2(p2+tan2θ)=1(1+p−2tan2θ),

ϖ2λ12=x2λ12+y2λ12

→

1−p2(z2λ12)=[1−1(1+p−2tan2θ)].

We also see that,

ϖ2p2z2

→

(1+p−2tan2θ)[1−1(1+p−2tan2θ)]=p−2tan2θ.

Again Consider Full 3D Ellipsoid

Let's try to replace everywhere, [ϖ/(pz)]2=p−2tan2θ with λ2. This gives,

𝒟2p2

≡

[(1+q2tan2φ)λ2−(q2−1)tan2φ].

which means that,

p2z2λ12

=

(1+tan2φ)[(1+q2tan2φ)λ2−(q2−1)tan2φ]

 

=

(1+tan2φ)/tan2φ[q2λ2(1+q2tan2φ)/(q2tan2φ)−(q2−1)]=1/sin2φ[q2λ2Q2−(q2−1)],

q2y2λ12

=

q2tan2φ(1+q2tan2φ)−q2tan2φ(1+tan2φ)(1+q2tan2φ)[(1+q2tan2φ)λ2−(q2−1)tan2φ]

 

=

q2tan2φ(1+q2tan2φ){1−(1+tan2φ)[(1+q2tan2φ)λ2−(q2−1)tan2φ]}

 

=

q2tan2φ(1+q2tan2φ){1−(1+tan2φ)/tan2φ[q2λ2(1+q2tan2φ)/(q2tan2φ)−(q2−1)]}

 

=

1Q2{1−1/sin2φ[q2λ2Q2−(q2−1)]}=1Q2[][[]−1sin2φ],

x2λ12

=

1−q2tan2φ(1+q2tan2φ){1−(1+tan2φ)[(1+q2tan2φ)λ2−(q2−1)tan2φ]}

 

 

−(1+tan2φ)[(1+q2tan2φ)λ2−(q2−1)tan2φ]

 

=

1−q2tan2φ(1+q2tan2φ)−{(1+tan2φ)[(1+q2tan2φ)λ2−(q2−1)tan2φ]}

 

=

1−1Q2−{1/sin2φ[q2λ2Q2−(q2−1)]}=1Q2[]{Q2[]−[]−Q2sin2φ},

x2+q2y2λ12

=

1−[1+q2tan2φ(1+q2tan2φ)]{(1+tan2φ)[(1+q2tan2φ)λ2−(q2−1)tan2φ]}.

Now, notice that,

q2y2Q2λ12

=

1−1[]sin2φ,

and,

x2λ12+1Q2

=

1−1[]sin2φ.

Hence,

x2λ12+1Q2

=

q2y2Q2λ12

⇒0

=

Q4(q2y2λ12)−Q2(x2λ12)−1

 

=

Q4−Q2(x2q2y2)−(λ12q2y2),

where,

Q2

≡

1+q2tan2φq2tan2φ.

Solving the quadratic equation, we have,

Q2

=

12{(x2q2y2)±[(x2q2y2)2+4(λ12q2y2)]1/2}

 

=

(x22q2y2){1±[1+4(λ12q2y2x4)]1/2}.

Tentative Summary

λ1

≡

(x2+q2y2+p2z2)1/2,

λ2

≡

(x2+y2)1/2pz,

λ3=Q2

≡

(x22q2y2){1±[1+4(λ12q2y2x4)]1/2}.

Partial Derivatives & Scale Factors

First Coordinate

∂λ1∂x

=

xλ1,

∂λ1∂y

=

q2yλ1,

∂λ1∂z

=

p2zλ1.

h1−2

=

(∂λ1∂x)2+(∂λ1∂y)2+(∂λ1∂z)2

=

(xλ1)2+(q2yλ1)2+(p2zλ1)2.

⇒h1

=

λ1ℓ3D,

where,

ℓ3D≡(x2+q4y2+p4z2)−1/2.

As a result, the associated unit vector is,

e^1

=

ı^h1(∂λ1∂x)+ȷ^h1(∂λ1∂y)+k^h1(∂λ1∂z)

 

=

ı^xℓ3D+ȷ^q2yℓ3D+k^p2zℓ3D.

Notice that,

e^1⋅e^1

=

(x2+q4y2+p4z2)ℓ3D2=1.


Second Coordinate (1st Try)

∂λ2∂x

=

1pz[x(x2+y2)1/2],

∂λ2∂y

=

1pz[y(x2+y2)1/2],

∂λ2∂z

=

−(x2+y2)1/2pz2.

h2−2

=

(∂λ2∂x)2+(∂λ2∂y)2+(∂λ2∂z)2

 

=

{1pz[x(x2+y2)1/2]}2+{1pz[y(x2+y2)1/2]}2+{(x2+y2)1/2pz2}2

 

=

{[x2(x2+y2)p2z2]}+{[y2(x2+y2)p2z2]}+{(x2+y2)p2z4}

 

=

1p2z2+(x2+y2)p2z4=(x2+y2+z2)p2z4

⇒h2

=

pz2r

As a result, the associated unit vector is,

e^2

=

ı^h2(∂λ2∂x)+ȷ^h2(∂λ2∂y)+k^h2(∂λ2∂z)

 

=

ı^[xzr(x2+y2)1/2]+ȷ^[yzr(x2+y2)1/2]−k^[(x2+y2)1/2r].

Notice that,

e^2⋅e^2

=

[x2z2r2(x2+y2)]+[y2z2r2(x2+y2)]+[(x2+y2)r2]=1.

Let's check to see if this "second" unit vector is orthogonal to the "first."

e^1⋅e^2

=

xℓ3D[xzr(x2+y2)1/2]+q2yℓ3D[yzr(x2+y2)1/2]−p2zℓ3D[(x2+y2)1/2r]

 

=

ℓ3D{[x2zr(x2+y2)1/2]+[q2y2zr(x2+y2)1/2]−[p2z(x2+y2)1/2r]}

 

=

zℓ3Dr(x2+y2)1/2{[x2]+[q2y2]−[p2(x2+y2)]}

 

≠

0 .


Second Coordinate (2nd Try)

Let's try,

λ2

=

[(x2+q2y2+𝔣⋅p2z2)1/2pz],

⇒∂λ2∂x

=

xpz(x2+q2y2+𝔣⋅p2z2)1/2=xp2z2λ2,

∂λ2∂y

=

q2ypz(x2+q2y2+𝔣⋅p2z2)1/2=q2yp2z2λ2,

∂λ2∂z

=

𝔣⋅p2zpz(x2+q2y2+𝔣⋅p2z2)1/2−(x2+q2y2+𝔣⋅p2z2)1/2pz2=1p2z2λ2(𝔣⋅p2z)−λ2z=1p2z2λ2(𝔣⋅p2z−p2zλ22).

Hence,

h2−2

=

(∂λ2∂x)2+(∂λ2∂y)2+(∂λ2∂z)2

 

=

[xp2z2λ2]2+[q2yp2z2λ2]2+[𝔣zλ2−λ2z]2

 

=

[x2+q4y2p4z4λ22]+[1zλ2(𝔣−λ22)]2

 

=

1p4z4λ22[x2+q4y2+p4z2(𝔣−λ22)2]

⇒h2

=

p2z2λ2[x2+q4y2+p4z2(𝔣−λ22)2]1/2.

So, the associated unit vector is,

e^2

=

ı^h2(∂λ2∂x)+ȷ^h2(∂λ2∂y)+k^h2(∂λ2∂z)

 

=

ı^{x[x2+q4y2+p4z2(𝔣−λ22)2]1/2}+ȷ^{q2y[x2+q4y2+p4z2(𝔣−λ22)2]1/2}+k^{p2z(𝔣−λ22)[x2+q4y2+p4z2(𝔣−λ22)2]1/2}.

Checking orthogonality …

e^1⋅e^2

=

xℓ3D{x[x2+q4y2+p4z2(𝔣−λ22)2]1/2}+q2yℓ3D{q2y[x2+q4y2+p4z2(𝔣−λ22)2]1/2}+p2zℓ3D{p2z(𝔣−λ22)[x2+q4y2+p4z2(𝔣−λ22)2]1/2}

 

=

ℓ3D[x2+q4y2+p4z2(𝔣−λ22)2]1/2{x2+q4y2+p4z2(𝔣−λ22)}.

 

=

ℓ3D[x2+q4y2+p4z2(𝔣−λ22)2]1/2{x2+q4y2+p4z2(𝔣−λ22)}.

If 𝔣=0, we have …

p2z(𝔣−λ22)

    →   

[−p2zλ22]𝔣=0=−p2z[(x2+q2y2+𝔣⋅0p2z2)1/2pz]2=−(x2+q2y2)z,

which, in turn, means …

[x2+q4y2+p4z2(𝔣0−λ22)2]1/2

=

[x2+q4y2+p4z2λ24]1/2

 

=

{x2+q4y2+p4z2[(x2+q2y2+𝔣⋅0p2z2)1/2pz]4}1/2

 

=

{x2+q4y2+[(x2+q2y2)2z2]}1/2

 

=

(x2+q4y2)1/2[1+(x2+q2y2)z2]1/2

 

=

(x2+q4y2)1/2z[z2+(x2+q2y2)]1/2,

and,

e^1⋅e^2

=

ℓ3D[x2+q4y2+p4z2(𝔣−λ22)2]1/2{x2+q4y2+p4z2(𝔣−λ22)}.

Speculation6

Determine λ2

This is very similar to the above, Speculation2. Try,

λ2

=

xy1/q2z2/p2,

in which case,

∂λ2∂x

=

λ2x,

∂λ2∂y

=

xz2/p2(1q2)y1/q2−1=λ2q2y,

∂λ2∂z

=

−2λ2p2z.

The associated scale factor is, then,

h2

=

[(∂λ2∂x)2+(∂λ2∂y)2+(∂λ2∂z)2]−1/2

 

=

[(λ2x)2+(λ2q2y)2+(−2λ2p2z)2]−1/2

 

=

1λ2[1x2+1q4y2+4p4z2]−1/2

 

=

1λ2[(q4y2p4z2+x2p4z2+4x2q4y2x2q4y2p4z2]−1/2

 

=

1λ2[xq2yp2z𝒟].

where,

𝒟

≡

(q4y2p4z2+x2p4z2+4x2q4y2)1/2.

The associated unit vector is, then,

e^2

=

ı^h2(∂λ2∂x)+ȷ^h2(∂λ2∂y)+k^h2(∂λ2∂z)

 

=

xq2yp2z𝒟{ı^(1x)+ȷ^(1q2y)+k^(−2p2z)} .

Recalling that the unit vector associated with the "first" coordinate is,

e^1

≡

ı^(xℓ3D)+ȷ^(q2yℓ3D)+k^(p2zℓ3D),

where,

ℓ3D

=

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

let's check to see whether the "second" unit vector is orthogonal to the "first."

e^1⋅e^2

=

(xq2yp2z)ℓ3D𝒟[1+1−2]=0.

Hooray!


Direction Cosines for Third Unit Vector

Now, what is the unit vector, e^3, that is simultaneously orthogonal to both these "first" and the "second" unit vectors?

e^3≡e^1×e^2

=

ı^[(e1y)(e2z)−(e2y)(e1z))]+ȷ^[(e1z)(e2x)−(e2z)(e1x))]+k^[(e1x)(e2y)−(e2x)(e1y))]

 

=

(xq2yp2z)ℓ3D𝒟{ı^[(−2q2yp2z)−(p2zq2y)]+ȷ^[(p2zx)−(−2xp2z)]+k^[(xq2y)−(q2yx)]}

 

=

(xq2yp2z)ℓ3D𝒟{−ı^[2q4y2+p4z2q2yp2z]+ȷ^[p4z2+2x2xp2z]+k^[x2−q4y2xq2y]}

 

=

ℓ3D𝒟{−ı^[x(2q4y2+p4z2)]+ȷ^[q2y(p4z2+2x2)]+k^[p2z(x2−q4y2)]}.

Is this a valid unit vector? First, note that …

(ℓ3D𝒟)−2

=

(q4y2p4z2+x2p4z2+4x2q4y2)(x2+q4y2+p4z2)

 

=

(x2q4y2p4z2+x4p4z2+4x4q4y2)+(q8y4p4z2+x2q4y2p4z2+4x2q8y4)+(q4y2p8z4+x2p8z4+4x2q4y2p4z2)

 

=

6x2q4y2p4z2+x4(p4z2+4q4y2)+q8y4(p4z2+4x2)+p8z4(x2+q4y2).

Then we have,

(ℓ3D𝒟)−2e^3⋅e^3

=

[x(2q4y2+p4z2)]2+[q2y(p4z2+2x2)]2+[p2z(x2−q4y2)]2

 

=

x2(4q8y4+4q4y2p4z2+p8z4)+q4y2(p8z4+4x2p4z2+4x4)+p4z2(x4−2x2q4y2+q8y4)

 

=

4x2q8y4+4x2q4y2p4z2+x2p8z4+q4y2p8z4+4x2q4y2p4z2+4x4q4y2+x4p4z2−2x2q4y2p4z2+q8y4p4z2

 

=

6x2q4y2p4z2+p8z4(x2+q4y2)+x4(4q4y2+p4z2)+q8y4(4x2+p4z2)

 

=

(ℓ3D𝒟)−2,

which means that, e^3⋅e^3=1.    Hooray! Again (11/11/2020)!

Direction Cosine Components for T6 Coordinates
n λn hn ∂λn∂x ∂λn∂y ∂λn∂z γn1 γn2 γn3
1 (x2+q2y2+p2z2)1/2 λ1ℓ3D xλ1 q2yλ1 p2zλ1 (x)ℓ3D (q2y)ℓ3D (p2z)ℓ3D
2 xy1/q2z2/p2 1λ2[xq2yp2z𝒟] λ2x λ2q2y −2λ2p2z xq2yp2z𝒟(1x) xq2yp2z𝒟(1q2y) xq2yp2z𝒟(−2p2z)
3 --- --- --- --- --- −ℓ3D𝒟[x(2q4y2+p4z2)] ℓ3D𝒟[q2y(p4z2+2x2)] ℓ3D𝒟[p2z(x2−q4y2)]

ℓ3D

≡

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

𝒟

≡

(q4y2p4z2+x2p4z2+4x2q4y2)1/2.


Let's double-check whether this "third" unit vector is orthogonal to both the "first" and the "second" unit vectors.

e^1⋅e^3

=

ℓ3D2𝒟{−x[x(2q4y2+p4z2)]+q2y[q2y(p4z2+2x2)]+p2z[p2z(x2−q4y2)]}

 

=

ℓ3D2𝒟{−(2x2q4y2+x2p4z2)+(q4y2p4z2+2x2q4y2)+(x2p4z2−q4y2p4z2)}

 

=

0,

and,

e^2⋅e^3

=

ℓ3D𝒟⋅xq2yp2z𝒟{−[(2q4y2+p4z2)]+[(p4z2+2x2)]−[2(x2−q4y2)]}

 

=

0.

Q. E. D.

Search for Third Coordinate Expression

Let's try …

λ3

=

𝒟nℓ3Dm

 

=

(q4y2p4z2+x2p4z2+4x2q4y2)n/2(x2+q4y2+p4z2)−m/2

⇒∂λ3∂xi

=

ℓ3Dm[n2(q4y2p4z2+x2p4z2+4x2q4y2)n/2−1]∂∂xi[(q4y2p4z2+x2p4z2+4x2q4y2)]

 

 

+𝒟n[−m2(x2+q4y2+p4z2)−m/2−1]∂∂xi[(x2+q4y2+p4z2)]

 

=

𝒟nℓ3Dm[n2(q4y2p4z2+x2p4z2+4x2q4y2)−1]∂∂xi[(q4y2p4z2+x2p4z2+4x2q4y2)]

 

 

+𝒟nℓ3Dm[−m2(x2+q4y2+p4z2)−1]∂∂xi[(x2+q4y2+p4z2)].


Hence,

1𝒟nℓ3Dm⋅∂λ3∂x

=

n2𝒟2∂∂x[q4y2p4z2+x2p4z2+4x2q4y2]−mℓ3D22∂∂x[x2+q4y2+p4z2]

 

=

x{n𝒟2[p4z2+4q4y2]−mℓ3D2}

 

=

x{n(p4z2+4q4y2)(q4y2p4z2+x2p4z2+4x2q4y2)−m(x2+q4y2+p4z2)}

This is overly cluttered! Let's try, instead …

A≡ℓ3D−2

=

(x2+q4y2+p4z2),

      and,      

B≡𝒟2

=

(q4y2p4z2+x2p4z2+4x2q4y2)

⇒∂A∂x

=

2x,

∂A∂y

=

2q4y,

∂A∂z

=

2p4z;

∂B∂x

=

2x(4q4y2+p4z2),

∂B∂y

=

2q4y(p4z2+4x2),

∂B∂z

=

2p4z(q4y2+x2).


Now, let's assume that,

λ3

≡

(AB)1/2,

⇒∂λ3∂xi

=

12(AB)1/2⋅∂A∂xi−A1/22B3/2⋅∂B∂xi

 

=

λ32AB[B⋅∂A∂xi−A⋅∂B∂xi].

⇒[2ABλ3]∂λ3∂xi

=

(q4y2p4z2+x2p4z2+4x2q4y2)⋅∂A∂xi−(x2+q4y2+p4z2)⋅∂B∂xi.

Looking ahead …

h3−2

=

{λ32AB[B⋅∂A∂x−A⋅∂B∂x]}2+{λ32AB[B⋅∂A∂y−A⋅∂B∂y]}2+{λ32AB[B⋅∂A∂z−A⋅∂B∂z]}2

⇒[2ABλ3]2h3−2

=

[B⋅∂A∂x−A⋅∂B∂x]2+[B⋅∂A∂y−A⋅∂B∂y]2+[B⋅∂A∂z−A⋅∂B∂z]2

⇒[λ32AB]h3

=

{[B⋅∂A∂x−A⋅∂B∂x]2+[B⋅∂A∂y−A⋅∂B∂y]2+[B⋅∂A∂z−A⋅∂B∂z]2}−1/2

Then, for example,

γ31≡h3(∂λ3∂x)

=

[B⋅∂A∂x−A⋅∂B∂x]{[B⋅∂A∂x−A⋅∂B∂x]2+[B⋅∂A∂y−A⋅∂B∂y]2+[B⋅∂A∂z−A⋅∂B∂z]2}−1/2

As a result, we have,

⇒[2ABλ3]∂λ3∂x

=

2x[(q4y2p4z2+x2p4z2+4x2q4y2)−(x2+q4y2+p4z2)(4q4y2+p4z2)]

 

=

2x[(q4y2p4z2+x2p4z2+4x2q4y2)−(4x2q4y2+x2p4z2+4q8y4+q4y2p4z2+4q4y2p4z2+p8z4)]

 

=

2x[−(4q8y4+4q4y2p4z2+p8z4)]

 

=

−2x(2q4y2+p4z2)2

 

=

−8x(q4y2+p4z22)2

⇒[AB]∂ln⁡λ3∂ln⁡x

=

−[2x(q4y2+p4z22)]2;

and,

⇒[2ABλ3]∂λ3∂y

=

2q4y[(q4y2p4z2+x2p4z2+4x2q4y2)−(x2+q4y2+p4z2)(p4z2+4x2)]

 

=

2q4y[(q4y2p4z2+x2p4z2+4x2q4y2)−(x2p4z2+4x4+q4y2p4z2+4x2q4y2+p8z4+4x2p4z2)]

 

=

−2q4y(4x4+p8z4+4x2p4z2)

 

=

−2q4y(2x2+p4z2)2

⇒[AB]∂ln⁡λ3∂ln⁡y

=

−[2q2y(x2+p4z22)]2;

and,

⇒[2ABλ3]∂λ3∂z

=

2p4z[(q4y2p4z2+x2p4z2+4x2q4y2)−(x2+q4y2+p4z2)(q4y2+x2)]

 

=

2p4z[(q4y2p4z2+x2p4z2+4x2q4y2)−(x2q4y2+x4+q8y4+x2q4y2+q4y2p4z2+x2p4z2)]

 

=

2p4z[(2x2q4y2)−(x4+q8y4)]

 

=

−2p4z[x4+q8y4−2x2q4y2]

 

=

−2p4z(x2−q4y2)2

⇒[AB]∂ln⁡λ3∂ln⁡z

=

−4[(p4z24)(x2−q4y2)2]

 

=

−[2(p2z2)(x2−q4y2)]2.

Wow!   Really close! (13 November 2020)

Just for fun, let's see what we get for h3. It is given by the expression,

h3−2

=

(∂λ3∂x)2+(∂λ3∂y)2+(∂λ3∂z)2

 

=

{λ3ABx[2x(q4y2+p4z22)]2}2+{λ3ABy[2q2y(x2+p4z22)]2}2+{λ3ABz[2(p2z2)(x2−q4y2)]2}2

⇒[ABλ3]2h3−2

=

{1x2[2x(q4y2+p4z22)]4}+{1y2[2q2y(x2+p4z22)]4}+{1z2[2(p2z2)(x2−q4y2)]4}

Fiddle Around

Let …

ℒx

≡

−[B⋅∂A∂x−A⋅∂B∂x]

=

8x(q4y2+p4z22)2

=

2x[2x(q4y2+p4z22)]2

=

8x𝔉x(y,z)

ℒy

≡

−[B⋅∂A∂y−A⋅∂B∂y]

=

8q4y(x2+p4z22)2

=

2y[2q2y(x2+p4z22)]2

=

8y𝔉y(x,z)

ℒz

≡

−[B⋅∂A∂z−A⋅∂B∂z]

=

2p4z(x2−q4y2)2

=

2z[p2z(x2−q4y2)]2

=

8z𝔉z(x,y)

With this shorthand in place, we can write,

e^3

=

ℓ3D𝒟{−ı^[x(2q4y2+p4z2)]+ȷ^[q2y(p4z2+2x2)]+k^[p2z(x2−q4y2)]}

 

=

1(AB)1/2{−ı^[xℒx2]1/2+ȷ^[yℒy2]1/2+k^[zℒz2]1/2}.

We therefore also recognize that,

h3(∂λ3∂x)

=

−1(AB)1/2[xℒx2]1/2

=

−1(AB)1/2[4x2𝔉x]1/2,

h3(∂λ3∂y)

=

1(AB)1/2[yℒy2]1/2

=

1(AB)1/2[4y2𝔉y]1/2,

h3(∂λ3∂z)

=

1(AB)1/2[zℒz2]1/2

=

1(AB)1/2[4z2𝔉z]1/2.

Now, if — and it is a BIG "if" — h3=h0(AB)−1/2, then we have,

h0(∂λ3∂x)

=

−[4x2𝔉x]1/2

=

−2x[𝔉x]1/2,

h0(∂λ3∂y)

=

[4y2𝔉y]1/2

=

2y[𝔉y]1/2,

h0(∂λ3∂z)

=

[4z2𝔉z]1/2

=

2z[𝔉z]1/2,

⇒h0λ3

=

−x2[𝔉x]1/2+y2[𝔉y]1/2+z2[𝔉z]1/2.

But if this is the correct expression for λ3 and its three partial derivatives, then it must be true that,

h3−2

=

(∂λ3∂x)2+(∂λ3∂y)2+(∂λ3∂z)2

⇒(h3h0)−2

=

4x2[𝔉x]+4y2[𝔉y]+4z2[𝔉z]

 

=

4x2[(q4y2+p4z22)2]+4y2[q4(x2+p4z22)2]+4z2[p44(x2−q4y2)2]

 

=

x2(2q4y2+p4z2)2+q4y2(2x2+p4z2)2+p4z2(x2−q4y2)2

Well … the right-hand side of this expression is identical to the right-hand side of the above expression, where we showed that it equals (ℓ3D/𝒟)−2. That is to say, we are now showing that,

(h3h0)−2

=

(ℓ3D𝒟)−2=[AB]

⇒h3h0

=

(AB)−1/2.

And this is precisely what, just a few lines above, we hypothesized the functional expression for h3 ought to be. EUREKA!

Summary

In summary, then …

λ3

≡

−x2[𝔉x]1/2+y2[𝔉y]1/2+z2[𝔉z]1/2

 

=

−x2[(q4y2+p4z22)2]1/2+y2[q4(x2+p4z22)2]1/2+z2[p44(x2−q4y2)2]1/2

 

=

−x2(q4y2+p4z22)+q2y2(x2+p4z22)+p2z22(x2−q4y2),

and,


h3

=

(AB)−1/2

 

=

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

No! Once again this does not work. The direction cosines — and, hence, the components of the e^3 unit vector — are not correct!

Speculation7

Direction Cosine Components for T6 Coordinates
n λn hn ∂λn∂x ∂λn∂y ∂λn∂z γn1 γn2 γn3
1 (x2+q2y2+p2z2)1/2 λ1ℓ3D xλ1 q2yλ1 p2zλ1 (x)ℓ3D (q2y)ℓ3D (p2z)ℓ3D
2 xy1/q2z2/p2 1λ2[xq2yp2z𝒟] λ2x λ2q2y −2λ2p2z xq2yp2z𝒟(1x) xq2yp2z𝒟(1q2y) xq2yp2z𝒟(−2p2z)
3 --- --- --- --- --- −ℓ3D𝒟[x(2q4y2+p4z2)] ℓ3D𝒟[q2y(p4z2+2x2)] ℓ3D𝒟[p2z(x2−q4y2)]

ℓ3D

≡

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

𝒟

≡

(q4y2p4z2+x2p4z2+4x2q4y2)1/2.


On my white-board I have shown that, if

λ3≡ℓ3D𝒟,

then everything will work out as long as,

ℒ

=

(𝒟ℓ3D)21ℓ3D4,

where,

ℒ

≡

x2(2q4y2+p4z2)4+q8y2(2x2+p4z2)4+p8z2(x2−q4y2)4

 

=

x2(4q8y4+4q4y2p4z2+p8z4)2+q8y2(4x4+4x2p4z2+p8z4)2+p8z2(x4−2x2q4y2+q8y4)2

 

=

x2[16q16y8+16q12y6p4z2+4q8y4p8z4+16q12y6p4z2+16q8y4p8z4+4q4y2p12z6+4q8y4p8z4+4q4y2p12z6+p16z8]

 

 

+q8y2[16x8+16x6p4z2+4x4p8z4+16x6p4z2+16x4p8z4+4x2p12z6+4x4p8z4+4x2p12z6+p16z8]

 

 

+p8z2[x8−2x6q4y2+x4q8y4−2x6q4y2+4x4q8y4−2x2q12y6+x4q8y4−2x2q12y6+q16y8]

 

=

x2[16q16y8+32q12y6p4z2+24q8y4p8z4+8q4y2p12z6+p16z8]

 

 

+q8y2[16x8+32x6p4z2+24x4p8z4+8x2p12z6+p16z8]

 

 

+p8z2[x8−4x6q4y2+6x4q8y4−4x2q12y6+q16y8]

Let's check this out.

RHS≡(𝒟ℓ3D)21ℓ3D4

=

(q4y2p4z2+x2p4z2+4x2q4y2)(x2+q4y2+p4z2)3

 

=

(x2+q4y2+p4z2)(q4y2p4z2+x2p4z2+4x2q4y2)[x4+2x2q4y2+2x2p4z2+q8y4+2q4y2p4z2+p8z4]

 

=

[(6x2q4y2p4z2+x4p4z2+4x4q4y2)+(q8y4p4z2+4x2q8y4)+(q4y2p8z4+x2p8z4)]

 

 

×[x4+2x2q4y2+2x2p4z2+q8y4+2q4y2p4z2+p8z4]

 

=

[6x2q4y2p4z2+x4(p4z2+4q4y2)+q8y4(p4z2+4x2)+p8z4(q4y2+x2)]

 

 

×[x4+2x2q4y2+2x2p4z2+q8y4+2q4y2p4z2+p8z4]

Tiled Menu

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