Appendix/Ramblings/T6CoordinatesPt2: Difference between revisions
No edit summary |
No edit summary |
||
| Line 1,747: | Line 1,747: | ||
<font color="red">'''Hooray!'''</font> | <font color="red">'''Hooray!'''</font> | ||
====Direction Cosines for <i>Third</i> Unit Vector==== | |||
Now, what is the unit vector, <math>~\hat{e}_3</math>, that is simultaneously orthogonal to both these "first" and the "second" unit vectors? | |||
<table border="0" cellpadding="5" align="center"> | |||
<tr> | |||
<td align="right"> | |||
<math>~\hat{e}_3 \equiv \hat{e}_1 \times \hat{e}_2</math> | |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
\hat\imath \biggl[ ( e_{1y} )( e_{2z}) - ( e_{2y} )( e_{1z}) ) \biggl] | |||
+ \hat\jmath \biggl[ ( e_{1z} )( e_{2x}) - ( e_{2z} )( e_{1x}) ) \biggl] | |||
+ \hat{k} \biggl[ ( e_{1x} )( e_{2y}) - ( e_{2x} )( e_{1y}) ) \biggl] | |||
</math> | |||
</td> | |||
</tr> | |||
<tr> | |||
<td align="right"> | |||
| |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~\frac{(x q^2 y p^2 z) \ell_{3D}}{\mathcal{D}} | |||
\biggl\{ | |||
\hat\imath \biggl[ \biggl( -\frac{2 q^2y}{p^2 z} \biggr) - \biggl( \frac{p^2z}{q^2y} \biggr) \biggl] | |||
+ \hat\jmath \biggl[ \biggl( \frac{p^2z}{x} \biggr) - \biggl(-\frac{2x}{p^2z} \biggr) \biggl] | |||
+ \hat{k} \biggl[ \biggl( \frac{x}{q^2y} \biggr) - \biggl( \frac{q^2y}{x} \biggr) \biggl] | |||
\biggr\} | |||
</math> | |||
</td> | |||
</tr> | |||
<tr> | |||
<td align="right"> | |||
| |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~\frac{(x q^2 y p^2 z) \ell_{3D}}{\mathcal{D}} | |||
\biggl\{ | |||
-\hat\imath \biggl[ \frac{2 q^4y^2 + p^4z^2}{q^2 y p^2 z} \biggl] | |||
+ \hat\jmath \biggl[ \frac{p^4z^2 + 2x^2}{xp^2 z} \biggl] | |||
+ \hat{k} \biggl[ \frac{x^2 - q^4y^2}{x q^2y} \biggl] | |||
\biggr\} | |||
</math> | |||
</td> | |||
</tr> | |||
<tr> | |||
<td align="right"> | |||
| |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~\frac{\ell_{3D}}{\mathcal{D}} | |||
\biggl\{ | |||
-\hat\imath \biggl[ x(2 q^4y^2 + p^4z^2 ) \biggl] | |||
+ \hat\jmath \biggl[ q^2 y(p^4z^2 + 2x^2 ) \biggl] | |||
+ \hat{k} \biggl[ p^2z( x^2 - q^4y^2 ) \biggl] | |||
\biggr\} \, . | |||
</math> | |||
</td> | |||
</tr> | |||
</table> | |||
Is this a valid unit vector? First, note that … | |||
<table border="0" cellpadding="5" align="center"> | |||
<tr> | |||
<td align="right"> | |||
<math>~\biggl( \frac{\ell_{3D}}{\mathcal{D}} \biggr)^{-2}</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 ) | |||
</math> | |||
</td> | |||
</tr> | |||
<tr> | |||
<td align="right"> | |||
| |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
(x^2 q^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 + x^2 q^4y^2 p^4 z^2 + 4x^2q^8y^4) | |||
+(q^4 y^2 p^8 z^4 + x^2 p^8 z^4 + 4x^2q^4y^2 p^4 z^2) | |||
</math> | |||
</td> | |||
</tr> | |||
<tr> | |||
<td align="right"> | |||
| |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
6x^2 q^4 y^2 p^4 z^2 + x^4(p^4 z^2 + 4 q^4y^2) | |||
+ q^8 y^4(p^4 z^2 + 4x^2) | |||
+p^8z^4(x^2 + q^4 y^2 )\, . | |||
</math> | |||
</td> | |||
</tr> | |||
</table> | |||
<span id="Eureka">Then we have,</span> | |||
<table border="0" cellpadding="5" align="center"> | |||
<tr> | |||
<td align="right"> | |||
<math>~\biggl( \frac{\ell_{3D}}{\mathcal{D}} \biggr)^{-2}\hat{e}_3 \cdot \hat{e}_3</math> | |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
\biggl[ x(2 q^4y^2 + p^4z^2 ) \biggl]^2 | |||
+ | |||
\biggl[ q^2 y(p^4z^2 + 2x^2 ) \biggl]^2 | |||
+ | |||
\biggl[ p^2z( x^2 - q^4y^2 ) \biggl]^2 | |||
</math> | |||
</td> | |||
</tr> | |||
<tr> | |||
<td align="right"> | |||
| |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
x^2(4 q^8y^4 + 4q^4y^2p^4z^2 + p^8z^4 ) | |||
+ | |||
q^4 y^2(p^8z^4 + 4x^2p^4z^2 + 4x^4 ) | |||
+ | |||
p^4z^2( x^4 - 2x^2q^4 y^2 + q^8y^4 ) | |||
</math> | |||
</td> | |||
</tr> | |||
<tr> | |||
<td align="right"> | |||
| |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
4 x^2 q^8y^4 + 4x^2 q^4y^2p^4z^2 + x^2 p^8z^4 | |||
+ | |||
q^4 y^2p^8z^4 + 4x^2q^4 y^2p^4z^2 + 4x^4q^4 y^2 | |||
+ | |||
x^4p^4z^2 - 2x^2q^4 y^2p^4z^2 + q^8y^4p^4z^2 | |||
</math> | |||
</td> | |||
</tr> | |||
<tr> | |||
<td align="right"> | |||
| |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
6x^2 q^4y^2p^4z^2 | |||
+ p^8z^4 (x^2 +q^4 y^2) | |||
+ x^4(4q^4 y^2 + p^4z^2) | |||
+ q^8 y^4(4 x^2 + p^4z^2 ) | |||
</math> | |||
</td> | |||
</tr> | |||
<tr> | |||
<td align="right"> | |||
| |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
\biggl( \frac{\ell_{3D}}{\mathcal{D}} \biggr)^{-2} \, , | |||
</math> | |||
</td> | |||
</tr> | |||
</table> | |||
which means that, <math>~\hat{e}_3\cdot \hat{e}_3 = 1</math>. <font color="red">'''Hooray! Again (11/11/2020)!'''</font> | |||
<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> | |||
Let's double-check whether this "third" unit vector is orthogonal to both the "first" and the "second" unit vectors. | |||
<table border="0" cellpadding="5" align="center"> | |||
<tr> | |||
<td align="right"> | |||
<math>~\hat{e}_1 \cdot \hat{e}_3</math> | |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~\frac{\ell_{3D}^2}{\mathcal{D}} | |||
\biggl\{ | |||
-x \biggl[ x(2 q^4y^2 + p^4z^2 ) \biggl] | |||
+ q^2 y \biggl[ q^2 y(p^4z^2 + 2x^2 ) \biggl] | |||
+ p^2 z \biggl[ p^2z( x^2 - q^4y^2 ) \biggl] | |||
\biggr\} | |||
</math> | |||
</td> | |||
</tr> | |||
<tr> | |||
<td align="right"> | |||
| |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~\frac{\ell_{3D}^2}{\mathcal{D}} | |||
\biggl\{ | |||
- (2 x^2q^4y^2 + x^2p^4z^2 ) | |||
+ (q^4 y^2 p^4z^2 + 2x^2 q^4 y^2) | |||
+ ( x^2p^4z^2 - q^4y^2 p^4z^2 ) | |||
\biggr\} | |||
</math> | |||
</td> | |||
</tr> | |||
<tr> | |||
<td align="right"> | |||
| |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
0 \, , | |||
</math> | |||
</td> | |||
</tr> | |||
</table> | |||
and, | |||
<table border="0" cellpadding="5" align="center"> | |||
<tr> | |||
<td align="right"> | |||
<math>~\hat{e}_2 \cdot \hat{e}_3</math> | |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~\frac{\ell_{3D}}{\mathcal{D}} \cdot \frac{x q^2 y p^2 z}{\mathcal{D}} | |||
\biggl\{ | |||
- \biggl[ (2 q^4y^2 + p^4z^2 ) \biggl] | |||
+ \biggl[ (p^4z^2 + 2x^2 ) \biggl] | |||
- \biggl[ 2( x^2 - q^4y^2 ) \biggl] | |||
\biggr\} | |||
</math> | |||
</td> | |||
</tr> | |||
<tr> | |||
<td align="right"> | |||
| |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
0 \, . | |||
</math> | |||
</td> | |||
</tr> | |||
</table> | |||
<font color="red">'''Q. E. D.'''</font> | |||
{{ SGFfooter }} | {{ SGFfooter }} | ||
Revision as of 16:44, 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.
|
Appendices: | VisTrailsEquations | VisTrailsVariables | References | Ramblings | VisTrailsImages | myphys.lsu | ADS | |