Appendix/Ramblings/T6CoordinatesPt2: Difference between revisions
| (3 intermediate revisions by the same user not shown) | |||
| Line 3,455: | Line 3,455: | ||
===Speculation7=== | ===Speculation7=== | ||
<table border="1" cellpadding="8" align="center"> | |||
<tr> | |||
<td align="center" colspan="9">'''Direction Cosine Components for T6 Coordinates'''</td> | |||
</tr> | |||
<tr> | |||
<td align="center"><math>~n</math></td> | |||
<td align="center"><math>~\lambda_n</math></td> | |||
<td align="center"><math>~h_n</math></td> | |||
<td align="center"><math>~\frac{\partial \lambda_n}{\partial x}</math></td> | |||
<td align="center"><math>~\frac{\partial \lambda_n}{\partial y}</math></td> | |||
<td align="center"><math>~\frac{\partial \lambda_n}{\partial z}</math></td> | |||
<td align="center"><math>~\gamma_{n1}</math></td> | |||
<td align="center"><math>~\gamma_{n2}</math></td> | |||
<td align="center"><math>~\gamma_{n3}</math></td> | |||
</tr> | |||
<tr> | |||
<td align="center"><math>~1</math></td> | |||
<td align="center"><math>~(x^2 + q^2 y^2 + p^2 z^2)^{1 / 2} </math></td> | |||
<td align="center"><math>~\lambda_1 \ell_{3D}</math></td> | |||
<td align="center"><math>~\frac{x}{\lambda_1}</math></td> | |||
<td align="center"><math>~\frac{q^2 y}{\lambda_1}</math></td> | |||
<td align="center"><math>~\frac{p^2 z}{\lambda_1}</math></td> | |||
<td align="center"><math>~(x) \ell_{3D}</math></td> | |||
<td align="center"><math>~(q^2 y)\ell_{3D}</math></td> | |||
<td align="center"><math>~(p^2z) \ell_{3D}</math></td> | |||
</tr> | |||
<tr> | |||
<td align="center"><math>~2</math></td> | |||
<td align="center"><math>~\frac{x y^{1/q^2}}{ z^{2/p^2}}</math></td> | |||
<td align="center"><math>~\frac{1}{\lambda_2}\biggl[\frac{x q^2 y p^2 z}{ \mathcal{D}}\biggr] </math></td> | |||
<td align="center"><math>~\frac{\lambda_2}{x}</math></td> | |||
<td align="center"><math>~\frac{\lambda_2}{q^2 y}</math></td> | |||
<td align="center"><math>~-\frac{2\lambda_2}{p^2 z}</math></td> | |||
<td align="center"><math>~\frac{x q^2 y p^2 z}{\mathcal{D}} \biggl(\frac{1}{x}\biggr)</math></td> | |||
<td align="center"><math>~ \frac{x q^2 y p^2 z}{\mathcal{D}} \biggl(\frac{1}{q^2y}\biggr)</math></td> | |||
<td align="center"><math>~\frac{x q^2 y p^2 z}{\mathcal{D}} \biggl(-\frac{2}{p^2z}\biggr)</math></td> | |||
</tr> | |||
<tr> | |||
<td align="center"><math>~3</math></td> | |||
<td align="center">---</td> | |||
<td align="center">---</td> | |||
<td align="center">---</td> | |||
<td align="center">---</td> | |||
<td align="center">---</td> | |||
<td align="center"><math>~-\frac{\ell_{3D}}{\mathcal{D}}\biggl[ x(2 q^4y^2 + p^4z^2 ) \biggl]</math></td> | |||
<td align="center"><math>~\frac{\ell_{3D}}{\mathcal{D}}\biggl[ q^2 y(p^4z^2 + 2x^2 ) \biggl]</math></td> | |||
<td align="center"><math>~\frac{\ell_{3D}}{\mathcal{D}}\biggl[ p^2z( x^2 - q^4y^2 ) \biggl]</math></td> | |||
</tr> | |||
<tr> | |||
<td align="left" colspan="9"> | |||
<table border="0" cellpadding="8" align="center"> | |||
<tr> | |||
<td align="right"> | |||
<math>~\ell_{3D}</math> | |||
</td> | |||
<td align="center"> | |||
<math>~\equiv</math> | |||
</td> | |||
<td align="left"> | |||
<math>~(x^2 + q^4y^2 + p^4 z^2 )^{- 1 / 2} \, ,</math> | |||
</td> | |||
</tr> | |||
<tr> | |||
<td align="right"> | |||
<math>~\mathcal{D}</math> | |||
</td> | |||
<td align="center"> | |||
<math>~\equiv</math> | |||
</td> | |||
<td align="left"> | |||
<math>~(q^4 y^2 p^4 z^2 + x^2 p^4 z^2 + 4x^2q^4y^2)^{1 / 2} \, .</math> | |||
</td> | |||
</tr> | |||
</table> | |||
</td> | |||
</tr> | |||
</table> | |||
On my white-board I have shown that, if | |||
<div align="center"> | |||
<math>~\lambda_3 \equiv \ell_{3D} \mathcal{D} \, ,</math> | |||
</div> | |||
then everything will work out as long as, | |||
<table border="0" cellpadding="5" align="center"> | |||
<tr> | |||
<td align="right"> | |||
<math>~\mathcal{L}</math> | |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
\biggl(\frac{\mathcal{D}}{\ell_{3D}} \biggr)^2 \frac{1}{\ell_{3D}^4} \, , | |||
</math> | |||
</td> | |||
</tr> | |||
</table> | |||
where, | |||
<table border="0" cellpadding="5" align="center"> | |||
<tr> | |||
<td align="right"> | |||
<math>~\mathcal{L}</math> | |||
</td> | |||
<td align="center"> | |||
<math>~\equiv</math> | |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
x^2 (2q^4 y^2 + p^4z^2 )^4 | |||
+ | |||
q^8 y^2 (2x^2 + p^4 z^2)^4 | |||
+ | |||
p^8z^2( x^2 - q^4y^2)^4 | |||
</math> | |||
</td> | |||
</tr> | |||
<tr> | |||
<td align="right"> | |||
| |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
x^2 (4q^8 y^4 + 4q^4y^2p^4z^2 + p^8z^4 )^2 | |||
+ | |||
q^8 y^2 (4x^4 + 4x^2p^4z^2 + p^8 z^4)^2 | |||
+ | |||
p^8z^2( x^4 - 2x^2q^4y^2 + q^8y^4)^2 | |||
</math> | |||
</td> | |||
</tr> | |||
<tr> | |||
<td align="right"> | |||
| |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
x^2 | |||
[16q^{16}y^8 + 16q^{12}y^6p^4z^2 + 4q^8y^4p^8z^4 + 16q^{12}y^6p^4z^2 + 16q^8y^4p^8z^4 + 4q^4y^2p^{12}z^6 + 4q^8y^4p^8z^4 + 4q^4y^2 p^{12}z^6 + p^{16}z^8] | |||
</math> | |||
</td> | |||
</tr> | |||
<tr> | |||
<td align="right"> | |||
| |||
</td> | |||
<td align="center"> | |||
| |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
+ | |||
q^8 y^2 | |||
[16x^8 + 16x^6p^4z^2 + 4x^4p^8z^4 + 16x^6p^4z^2 + 16x^4p^8z^4 + 4x^2p^{12}z^6 + 4x^4p^8z^4 + 4x^2p^{12}z^6 + p^{16}z^8] | |||
</math> | |||
</td> | |||
</tr> | |||
<tr> | |||
<td align="right"> | |||
| |||
</td> | |||
<td align="center"> | |||
| |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
+ | |||
p^8z^2 | |||
[x^8 - 2x^6q^4y^2 + x^4q^8y^4 - 2x^6q^4y^2 + 4x^4q^8y^4 - 2x^2q^{12}y^6 + x^4q^8y^4 - 2x^2q^{12}y^6 + q^{16}y^8] | |||
</math> | |||
</td> | |||
</tr> | |||
<tr> | |||
<td align="right"> | |||
| |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
x^2 | |||
[16q^{16}y^8 + 32q^{12}y^6p^4z^2 + 24q^8y^4p^8z^4 + 8q^4y^2p^{12}z^6 + p^{16}z^8] | |||
</math> | |||
</td> | |||
</tr> | |||
<tr> | |||
<td align="right"> | |||
| |||
</td> | |||
<td align="center"> | |||
| |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
+ | |||
q^8 y^2 | |||
[16x^8 + 32x^6p^4z^2 + 24x^4p^8z^4 + 8x^2p^{12}z^6 + p^{16}z^8] | |||
</math> | |||
</td> | |||
</tr> | |||
<tr> | |||
<td align="right"> | |||
| |||
</td> | |||
<td align="center"> | |||
| |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
+ | |||
p^8z^2 | |||
[x^8 - 4x^6q^4y^2 + 6x^4q^8y^4 - 4x^2q^{12}y^6 + q^{16}y^8] | |||
</math> | |||
</td> | |||
</tr> | |||
</table> | |||
Let's check this out. | |||
<table border="0" cellpadding="5" align="center"> | |||
<tr> | |||
<td align="right"> | |||
<math>~\mathrm{RHS}~\equiv \biggl(\frac{\mathcal{D}}{\ell_{3D}} \biggr)^2 \frac{1}{\ell_{3D}^4}</math> | |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~(q^4 y^2 p^4 z^2 + x^2 p^4 z^2 + 4x^2q^4y^2)(x^2 + q^4y^2 + p^4 z^2 )^3</math> | |||
</td> | |||
</tr> | |||
<tr> | |||
<td align="right"> | |||
| |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~(x^2 + q^4y^2 + p^4 z^2 )(q^4 y^2 p^4 z^2 + x^2 p^4 z^2 + 4x^2q^4y^2)[x^4 + 2x^2q^4y^2 + 2x^2p^4z^2 +q^8y^4 + 2q^4y^2p^4z^2 + p^8z^4]</math> | |||
</td> | |||
</tr> | |||
<tr> | |||
<td align="right"> | |||
| |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
[(6x^2q^4 y^2 p^4 z^2 + x^4 p^4 z^2 + 4x^4q^4y^2) | |||
+ | |||
(q^8 y^4 p^4 z^2 + 4x^2q^8 y^4) | |||
+ | |||
(q^4 y^2 p^8 z^4 + x^2 p^8 z^4 )] | |||
</math> | |||
</td> | |||
</tr> | |||
<tr> | |||
<td align="right"> | |||
| |||
</td> | |||
<td align="center"> | |||
| |||
</td> | |||
<td align="left"> | |||
<math>~\times [x^4 + 2x^2q^4y^2 + 2x^2p^4z^2 +q^8y^4 + 2q^4y^2p^4z^2 + p^8z^4]</math> | |||
</td> | |||
</tr> | |||
<tr> | |||
<td align="right"> | |||
| |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
[6x^2q^4 y^2 p^4 z^2 + x^4(p^4 z^2 + 4q^4y^2) | |||
+ | |||
q^8 y^4(p^4 z^2 + 4x^2) | |||
+ | |||
p^8 z^4 (q^4 y^2 + x^2 )] | |||
</math> | |||
</td> | |||
</tr> | |||
<tr> | |||
<td align="right"> | |||
| |||
</td> | |||
<td align="center"> | |||
| |||
</td> | |||
<td align="left"> | |||
<math>~\times [x^4 + 2x^2q^4y^2 + 2x^2p^4z^2 +q^8y^4 + 2q^4y^2p^4z^2 + p^8z^4]</math> | |||
</td> | |||
</tr> | |||
</table> | |||
{{ SGFfooter }} | {{ SGFfooter }} | ||
Latest revision as of 18:13, 23 July 2021
Concentric Ellipsoidal (T6) Coordinates (Part 2)
Orthogonal Coordinates
Speculation5
Spherical Coordinates
|
|
|
|
|
|
|
|
|
|
|
|
Use λ1 Instead of r
Here, as above, we define,
|
|
|
|
Using this expression to eliminate "x" (in favor of λ1) in each of the three spherical-coordinate definitions, we obtain,
|
|
|
|
|
|
|
|
|
|
|
|
After a bit of additional algebraic manipulation, we find that,
|
|
|
|
|
|
|
|
|
|
|
|
where,
|
|
|
|
|
As a check, let's set , which should reduce to the normal spherical coordinate system.
|
Relationship To T3 Coordinates
If we set, , but continue to assume that , we expect to see a representation that resembles our previously discussed, T3 Coordinates. Note, for example, that the new "radial" coordinate is,
|
|
|
|
and, |
|
|
|
|
|
|
|
|
|
|
|
We also see that,
|
|
|
|
Again Consider Full 3D Ellipsoid
Let's try to replace everywhere, with . This gives,
|
|
|
|
which means that,
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Now, notice that,
|
|
|
|
and,
|
|
|
|
Hence,
|
|
|
|
|
|
|
|
|
|
|
|
where,
|
|
|
|
Solving the quadratic equation, we have,
|
|
|
|
|
|
|
|
|
Tentative Summary
|
Partial Derivatives & Scale Factors
First Coordinate
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
where,
As a result, the associated unit vector is,
|
|
|
|
|
|
|
|
Notice that,
|
|
|
|
Second Coordinate (1st Try)
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
As a result, the associated unit vector is,
|
|
|
|
|
|
|
|
Notice that,
|
|
|
|
Let's check to see if this "second" unit vector is orthogonal to the "first."
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Second Coordinate (2nd Try)
Let's try,
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Hence,
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
So, the associated unit vector is,
|
|
|
|
|
|
|
|
Checking orthogonality …
|
|
|
|
|
|
|
|
|
|
|
|
If , we have …
|
|
|
|
which, in turn, means …
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
and,
|
|
|
|
Speculation6
Determine λ2
This is very similar to the above, Speculation2. Try,
|
|
|
|
in which case,
|
|
|
|
|
|
|
|
|
|
|
|
The associated scale factor is, then,
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
where,
|
|
|
|
The associated unit vector is, then,
|
|
|
|
|
|
|
|
Recalling that the unit vector associated with the "first" coordinate is,
|
|
|
|
where,
|
|
|
|
let's check to see whether the "second" unit vector is orthogonal to the "first."
|
|
|
|
Hooray!
Direction Cosines for Third Unit Vector
Now, what is the unit vector, , that is simultaneously orthogonal to both these "first" and the "second" unit vectors?
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Is this a valid unit vector? First, note that …
|
|
|
|
|
|
|
|
|
|
|
|
Then we have,
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
which means that, . Hooray! Again (11/11/2020)!
| Direction Cosine Components for T6 Coordinates | ||||||||||||||
| --- | --- | --- | --- | --- | ||||||||||
|
||||||||||||||
Let's double-check whether this "third" unit vector is orthogonal to both the "first" and the "second" unit vectors.
|
|
|
|
|
|
|
|
|
|
|
|
and,
|
|
|
|
|
|
|
|
Q. E. D.
Search for Third Coordinate Expression
Let's try …
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Hence,
|
|
|
|
|
|
|
|
|
|
|
|
This is overly cluttered! Let's try, instead …
|
Now, let's assume that,
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Looking ahead …
Then, for example,
|
As a result, we have,
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
and,
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
and,
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Wow! Really close! (13 November 2020)
Just for fun, let's see what we get for . It is given by the expression,
|
|
|
|
|
|
|
|
|
|
|
|
Fiddle Around
Let …
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
With this shorthand in place, we can write,
|
|
|
|
|
|
|
|
We therefore also recognize that,
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Now, if — and it is a BIG "if" — , then we have,
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
But if this is the correct expression for and its three partial derivatives, then it must be true that,
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
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 . That is to say, we are now showing that,
|
|
|
|
|
|
|
|
And this is precisely what, just a few lines above, we hypothesized the functional expression for ought to be. EUREKA!
Summary
In summary, then …
|
|
|
|
|
|
|
|
|
|
|
|
and,
|
|
|
|
|
|
|
|
No! Once again this does not work. The direction cosines — and, hence, the components of the unit vector — are not correct!
Speculation7
| Direction Cosine Components for T6 Coordinates | ||||||||||||||
| --- | --- | --- | --- | --- | ||||||||||
|
||||||||||||||
On my white-board I have shown that, if
then everything will work out as long as,
|
|
|
|
where,
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Let's check this out.
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Appendices: | VisTrailsEquations | VisTrailsVariables | References | Ramblings | VisTrailsImages | myphys.lsu | ADS | |