Appendix/Ramblings/T6CoordinatesPt2: Difference between revisions
No edit summary |
|||
| Line 2,934: | Line 2,934: | ||
</tr> | </tr> | ||
</table> | </table> | ||
====Fiddle Around==== | |||
Let … | |||
<table border="0" cellpadding="5" align="center"> | |||
<tr> | |||
<td align="right"> | |||
<math>~\mathcal{L}_x</math> | |||
</td> | |||
<td align="center"> | |||
<math>~\equiv</math> | |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
- \biggl[ | |||
B \cdot \frac{\partial A}{\partial x} | |||
- | |||
A \cdot \frac{\partial B}{\partial x} | |||
\biggr] | |||
</math> | |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~8x \biggl(q^4y^2 + \frac{p^4z^2}{2} \biggr)^2 | |||
</math> | |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~\frac{2}{x} \biggl[ 2x \biggl(q^4y^2 + \frac{p^4z^2}{2} \biggr)\biggr]^2 | |||
</math> | |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~8x~\mathfrak{F}_x(y,z) | |||
</math> | |||
</td> | |||
</tr> | |||
<tr> | |||
<td align="right"> | |||
<math>~\mathcal{L}_y</math> | |||
</td> | |||
<td align="center"> | |||
<math>~\equiv</math> | |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
- \biggl[ | |||
B \cdot \frac{\partial A}{\partial y} | |||
- | |||
A \cdot \frac{\partial B}{\partial y} | |||
\biggr] | |||
</math> | |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
8q^4y\biggl( x^2 + \frac{p^4z^2}{2} \biggr)^2 | |||
</math> | |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
\frac{2}{y} \biggl[ 2q^2y\biggl( x^2 + \frac{p^4z^2}{2} \biggr)\biggr]^2 | |||
</math> | |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~8y~\mathfrak{F}_y(x,z) | |||
</math> | |||
</td> | |||
</tr> | |||
<tr> | |||
<td align="right"> | |||
<math>~\mathcal{L}_z</math> | |||
</td> | |||
<td align="center"> | |||
<math>~\equiv</math> | |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
-\biggl[ | |||
B \cdot \frac{\partial A}{\partial z} | |||
- | |||
A \cdot \frac{\partial B}{\partial z} | |||
\biggr] | |||
</math> | |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
2p^4 z \biggl(x^2 - q^4y^2 \biggr)^2 | |||
</math> | |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
\frac{2}{z} \biggl[p^2 z \biggl(x^2 - q^4y^2 \biggr)\biggr]^2 | |||
</math> | |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~8z~\mathfrak{F}_z(x,y) | |||
</math> | |||
</td> | |||
</tr> | |||
</table> | |||
With this shorthand in place, we can write, | |||
<table border="0" cellpadding="5" align="center"> | |||
<tr> | |||
<td align="right"> | |||
<math>~\hat{e}_3</math> | |||
</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> | |||
<tr> | |||
<td align="right"> | |||
| |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
\frac{1}{(AB)^{1 / 2}} | |||
\biggl\{ | |||
-\hat\imath \biggl[ \frac{x \mathcal{L}_x}{2} \biggl]^{1 / 2} | |||
+ \hat\jmath \biggl[ \frac{y \mathcal{L}_y}{2} \biggl]^{1 / 2} | |||
+ \hat{k} \biggl[ \frac{z \mathcal{L}_z}{2} \biggl]^{1 / 2} | |||
\biggr\} | |||
\, . | |||
</math> | |||
</td> | |||
</tr> | |||
</table> | |||
We therefore also recognize that, | |||
<table border="0" cellpadding="5" align="center"> | |||
<tr> | |||
<td align="right"> | |||
<math>~h_3 \biggl(\frac{\partial \lambda_3}{\partial x}\biggr)</math> | |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
-\frac{1}{(AB)^{1 / 2}} \biggl[ \frac{x \mathcal{L}_x}{2} \biggl]^{1 / 2} | |||
</math> | |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
-\frac{1}{(AB)^{1 / 2}} \biggl[ 4x^2 ~\mathfrak{F}_x \biggl]^{1 / 2} \, , | |||
</math> | |||
</td> | |||
</tr> | |||
<tr> | |||
<td align="right"> | |||
<math>~h_3 \biggl(\frac{\partial \lambda_3}{\partial y}\biggr)</math> | |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
\frac{1}{(AB)^{1 / 2}} \biggl[ \frac{y \mathcal{L}_y}{2} \biggl]^{1 / 2} | |||
</math> | |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
\frac{1}{(AB)^{1 / 2}} \biggl[ 4y^2~\mathfrak{F}_y \biggl]^{1 / 2} \, , | |||
</math> | |||
</td> | |||
</tr> | |||
<tr> | |||
<td align="right"> | |||
<math>~h_3 \biggl(\frac{\partial \lambda_3}{\partial z}\biggr)</math> | |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
\frac{1}{(AB)^{1 / 2}} \biggl[ \frac{z \mathcal{L}_z}{2} \biggl]^{1 / 2} | |||
</math> | |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
\frac{1}{(AB)^{1 / 2}} \biggl[ 4z^2~\mathfrak{F}_z \biggl]^{1 / 2} \, . | |||
</math> | |||
</td> | |||
</tr> | |||
</table> | |||
Now, if — and it is a BIG "if" — <math>~h_3 = h_0(AB)^{-1 / 2}</math>, then we have, | |||
<table border="0" cellpadding="5" align="center"> | |||
<tr> | |||
<td align="right"> | |||
<math>~h_0 \biggl(\frac{\partial \lambda_3}{\partial x}\biggr)</math> | |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
-\biggl[ 4x^2 ~\mathfrak{F}_x \biggl]^{1 / 2} | |||
</math> | |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
-2x \biggl[ \mathfrak{F}_x \biggl]^{1 / 2} \, , | |||
</math> | |||
</td> | |||
</tr> | |||
<tr> | |||
<td align="right"> | |||
<math>~h_0 \biggl(\frac{\partial \lambda_3}{\partial y}\biggr)</math> | |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
\biggl[ 4y^2~\mathfrak{F}_y \biggl]^{1 / 2} | |||
</math> | |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
2y \biggl[ \mathfrak{F}_y \biggl]^{1 / 2} \, , | |||
</math> | |||
</td> | |||
</tr> | |||
<tr> | |||
<td align="right"> | |||
<math>~h_0 \biggl(\frac{\partial \lambda_3}{\partial z}\biggr)</math> | |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
\biggl[ 4z^2~\mathfrak{F}_z \biggl]^{1 / 2} | |||
</math> | |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
2z\biggl[ \mathfrak{F}_z \biggl]^{1 / 2} \, , | |||
</math> | |||
</td> | |||
</tr> | |||
</table> | |||
<table border="0" cellpadding="5" align="center"> | |||
<tr> | |||
<td align="right"> | |||
<math>~\Rightarrow~~~ h_0 \lambda_3</math> | |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
-x^2 \biggl[ \mathfrak{F}_x \biggl]^{1 / 2} | |||
+ y^2 \biggl[ \mathfrak{F}_y \biggl]^{1 / 2} | |||
+ z^2 \biggl[ \mathfrak{F}_z \biggl]^{1 / 2} \, . | |||
</math> | |||
</td> | |||
</tr> | |||
</table> | |||
But if this is the correct expression for <math>~\lambda_3</math> and its three partial derivatives, then it must be true that, | |||
<table border="0" cellpadding="5" align="center"> | |||
<tr> | |||
<td align="right"> | |||
<math>~h_3^{-2}</math> | |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
\biggl(\frac{\partial \lambda_3}{\partial x}\biggr)^2 | |||
+ | |||
\biggl(\frac{\partial \lambda_3}{\partial y}\biggr)^2 | |||
+ | |||
\biggl(\frac{\partial \lambda_3}{\partial z}\biggr)^2 | |||
</math> | |||
</td> | |||
</tr> | |||
<tr> | |||
<td align="right"> | |||
<math>~\Rightarrow ~~~ \biggl( \frac{h_3}{h_0}\biggr)^{-2}</math> | |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
4x^2 \biggl[ \mathfrak{F}_x \biggl] | |||
+ | |||
4y^2 \biggl[ \mathfrak{F}_y \biggl] | |||
+ | |||
4z^2\biggl[ \mathfrak{F}_z \biggl] | |||
</math> | |||
</td> | |||
</tr> | |||
<tr> | |||
<td align="right"> | |||
| |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
4x^2 \biggl[ \biggl(q^4y^2 + \frac{p^4z^2}{2} \biggr)^2 \biggl] | |||
+ | |||
4y^2 \biggl[ q^4\biggl( x^2 + \frac{p^4z^2}{2} \biggr)^2\biggl] | |||
+ | |||
4z^2\biggl[ \frac{p^4}{4}\biggl(x^2 - q^4y^2 \biggr)^2 \biggl] | |||
</math> | |||
</td> | |||
</tr> | |||
<tr> | |||
<td align="right"> | |||
| |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
x^2 (2q^4y^2 + p^4z^2 )^2 | |||
+ | |||
q^4 y^2( 2x^2 + p^4z^2 )^2 | |||
+ | |||
p^4 z^2 (x^2 - q^4y^2 )^2 | |||
</math> | |||
</td> | |||
</tr> | |||
</table> | |||
Well … the right-hand side of this expression is identical to the right-hand side of the [[#Eureka|above expression]], where we showed that it equals <math>~(\ell_{3D}/\mathcal{D})^{-2}</math>. That is to say, we are now showing that, | |||
<table border="0" cellpadding="5" align="center"> | |||
<tr> | |||
<td align="right"> | |||
<math>~\biggl( \frac{h_3}{h_0}\biggr)^{-2}</math> | |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~\biggl( \frac{\ell_{3D}}{\mathcal{D}} \biggr)^{-2} = [AB]</math> | |||
</td> | |||
</tr> | |||
<tr> | |||
<td align="right"> | |||
<math>~\Rightarrow ~~~ \frac{h_3}{h_0}</math> | |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~(AB)^{-1 / 2} \, .</math> | |||
</td> | |||
</tr> | |||
</table> | |||
And this is ''precisely'' what, just a few lines above, we hypothesized the functional expression for <math>~h_3</math> ought to be. <font color="red">'''EUREKA!'''</font> | |||
{{ SGFfooter }} | {{ SGFfooter }} | ||
Revision as of 17:53, 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!
|
Appendices: | VisTrailsEquations | VisTrailsVariables | References | Ramblings | VisTrailsImages | myphys.lsu | ADS | |