Appendix/Ramblings/T6CoordinatesPt2: Difference between revisions

From jetwiki
Jump to navigation Jump to search
No edit summary
No edit summary
Line 742: Line 742:


</td></tr></table>
</td></tr></table>
====Partial Derivatives &amp; Scale Factors====
=====First Coordinate=====
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~\frac{\partial \lambda_1}{\partial x}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\frac{x}{\lambda_1} \, ,</math>
  </td>
</tr>
<tr>
  <td align="right">
<math>~\frac{\partial \lambda_1}{\partial y}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\frac{q^2 y}{\lambda_1} \, ,</math>
  </td>
</tr>
<tr>
  <td align="right">
<math>~\frac{\partial \lambda_1}{\partial z}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\frac{p^2 z}{\lambda_1} \, .</math>
  </td>
</tr>
</table>
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~h_1^{-2}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\biggl( \frac{\partial \lambda_1}{\partial x} \biggr)^2
+ \biggl( \frac{\partial \lambda_1}{\partial y} \biggr)^2
+ \biggl( \frac{\partial \lambda_1}{\partial z} \biggr)^2
</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\biggl( \frac{x}{\lambda_1} \biggr)^2
+ \biggl( \frac{q^2 y}{\lambda_1} \biggr)^2
+ \biggl( \frac{p^2 z}{\lambda_1} \biggr)^2
\, .</math>
  </td>
</tr>
<tr>
  <td align="right">
<math>~\Rightarrow ~~~ h_1</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left" colspan="3">
<math>~
\lambda_1 \ell_{3D} \, ,
</math>
  </td>
</tr>
</table>
where,
<div align="center"><math>\ell_{3D} \equiv (x^2 + q^4y^2 + p^4z^2)^{-1 / 2} \, .</math></div>
As a result, the associated unit vector is,
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~\hat{e}_1</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\hat{\imath} h_1 \biggl( \frac{\partial \lambda_1}{\partial x} \biggr)
+ \hat{\jmath} h_1 \biggl( \frac{\partial \lambda_1}{\partial y} \biggr)
+ \hat{k} h_1 \biggl( \frac{\partial \lambda_1}{\partial z} \biggr)
</math>
  </td>
</tr>
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\hat{\imath} x \ell_{3D}
+ \hat{\jmath} q^2 y\ell_{3D}
+ \hat{k} p^2 z \ell_{3D}
\, .
</math>
  </td>
</tr>
</table>
Notice that,
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~\hat{e}_1 \cdot \hat{e}_1</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
(x^2 + q^4 y^2 + p^4 z^2) \ell_{3D}^2 = 1 \, .
</math>
  </td>
</tr>
</table>
=====Second Coordinate (1<sup>st</sup> Try)=====
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~\frac{\partial \lambda_2}{\partial x}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\frac{1}{pz} \biggl[ \frac{x}{(x^2 + y^2)^{1 / 2}} \biggr] \, ,
</math>
  </td>
</tr>
<tr>
  <td align="right">
<math>~\frac{\partial \lambda_2}{\partial y}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\frac{1}{pz} \biggl[ \frac{y}{(x^2 + y^2)^{1 / 2}} \biggr]
\, ,</math>
  </td>
</tr>
<tr>
  <td align="right">
<math>~\frac{\partial \lambda_2}{\partial z}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
- \frac{(x^2 + y^2)^{1 / 2}}{pz^2}
\, .</math>
  </td>
</tr>
</table>
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~h_2^{-2}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\biggl( \frac{\partial \lambda_2}{\partial x} \biggr)^2
+ \biggl( \frac{\partial \lambda_2}{\partial y} \biggr)^2
+ \biggl( \frac{\partial \lambda_2}{\partial z} \biggr)^2
</math>
  </td>
</tr>
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\biggl\{ \frac{1}{pz} \biggl[ \frac{x}{(x^2 + y^2)^{1 / 2}} \biggr] \biggr\}^2
+ \biggl\{ \frac{1}{pz} \biggl[ \frac{y}{(x^2 + y^2)^{1 / 2}} \biggr] \biggr\}^2
+ \biggl\{ \frac{(x^2 + y^2)^{1 / 2}}{pz^2} \biggr\}^2
</math>
  </td>
</tr>
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\biggl\{ \biggl[ \frac{x^2}{(x^2 + y^2)p^2 z^2} \biggr] \biggr\}
+ \biggl\{ \biggl[ \frac{y^2}{(x^2 + y^2)p^2 z^2} \biggr] \biggr\}
+ \biggl\{ \frac{(x^2 + y^2)}{p^2 z^4} \biggr\}
</math>
  </td>
</tr>
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\frac{1}{p^2 z^2}
+ \frac{(x^2 + y^2)}{p^2 z^4}
=
\frac{(x^2 + y^2 + z^2)}{p^2 z^4}
</math>
  </td>
</tr>
<tr>
  <td align="right">
<math>~\Rightarrow~~~h_2</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\frac{p z^2}{r }
</math>
  </td>
</tr>
</table>
As a result, the associated unit vector is,
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~\hat{e}_2</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\hat{\imath} h_2 \biggl( \frac{\partial \lambda_2}{\partial x} \biggr)
+ \hat{\jmath} h_2 \biggl( \frac{\partial \lambda_2}{\partial y} \biggr)
+ \hat{k} h_2 \biggl( \frac{\partial \lambda_2}{\partial z} \biggr)
</math>
  </td>
</tr>
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\hat{\imath}  \biggl[ \frac{xz}{r(x^2 + y^2)^{1 / 2}} \biggr]
+ \hat{\jmath}  \biggl[ \frac{yz}{r(x^2 + y^2)^{1 / 2}} \biggr]
- \hat{k} \biggl[ \frac{(x^2 + y^2)^{1 / 2}}{r} \biggr] \, .
</math>
  </td>
</table>
Notice that,
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~\hat{e}_2 \cdot \hat{e}_2</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\biggl[ \frac{x^2 z^2}{r^2(x^2 + y^2)} \biggr]
+ \biggl[ \frac{y^2 z^2}{r^2(x^2 + y^2)} \biggr]
+ \biggl[ \frac{(x^2 + y^2)}{r^2} \biggr]
= 1 \, .
</math>
  </td>
</tr>
</table>
Let's check to see if this "second" unit vector is orthogonal to the "first."
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~\hat{e}_1 \cdot \hat{e}_2</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
x\ell_{3D} \biggl[ \frac{xz}{r(x^2 + y^2)^{1 / 2}} \biggr]
+ q^2 y\ell_{3D} \biggl[ \frac{yz}{r(x^2 + y^2)^{1 / 2}} \biggr]
- p^2 z \ell_{3D} \biggl[ \frac{(x^2 + y^2)^{1 / 2}}{r} \biggr]
</math>
  </td>
</tr>
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\ell_{3D} \biggl\{ \biggl[ \frac{x^2z}{r(x^2 + y^2)^{1 / 2}} \biggr]
+ \biggl[ \frac{q^2 y^2 z}{r(x^2 + y^2)^{1 / 2}} \biggr]
-  \biggl[ \frac{p^2 z(x^2 + y^2)^{1 / 2}}{r} \biggr]
\biggr\}
</math>
  </td>
</tr>
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\frac{z\ell_{3D}}{r (x^2 + y^2)^{1 / 2}} \biggl\{ \biggl[ x^2\biggr]
+ \biggl[ q^2 y^2 \biggr]
-  \biggl[ p^2 (x^2 + y^2) \biggr]
\biggr\}
</math>
  </td>
</tr>
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~\ne</math>
  </td>
  <td align="left">
<math>~
0 \ .
</math>
  </td>
</tr>
</table>
=====Second Coordinate (2<sup>nd</sup> Try)=====
Let's try,
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~\lambda_2</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\biggl[\frac{(x^2 + q^2y^2 + \mathfrak{f}\cdot p^2 z^2)^{1 / 2}}{pz}  \biggr]
\, ,
</math>
  </td>
</tr>
<tr>
  <td align="right">
<math>~\Rightarrow~~~\frac{\partial \lambda_2}{\partial x}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\frac{x}{pz(x^2 + q^2y^2 + \mathfrak{f}\cdot p^2 z^2)^{1 / 2} } = \frac{x}{p^2 z^2 \lambda_2}
\, ,
</math>
  </td>
</tr>
<tr>
  <td align="right">
<math>~\frac{\partial \lambda_2}{\partial y}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\frac{q^2 y}{pz(x^2 + q^2y^2 + \mathfrak{f}\cdot p^2 z^2)^{1 / 2} } = \frac{q^2y}{p^2 z^2 \lambda_2}
\, ,
</math>
  </td>
</tr>
<tr>
  <td align="right">
<math>~\frac{\partial \lambda_2}{\partial z}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\frac{\mathfrak{f}\cdot p^2z}{pz(x^2 + q^2y^2 + \mathfrak{f}\cdot p^2 z^2)^{1 / 2} }
-
\frac{(x^2 + q^2y^2 + \mathfrak{f}\cdot p^2 z^2)^{1 / 2}}{pz^2} 
= \frac{1}{p^2z^2 \lambda_2 } \biggl( \mathfrak{f}\cdot p^2z \biggr) - \frac{\lambda_2 }{z}
= \frac{1}{p^2z^2 \lambda_2 } \biggl( \mathfrak{f}\cdot p^2z  - p^2z \lambda_2^2 \biggr)
\, .
</math>
  </td>
</tr>
</table>
Hence,
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~h_2^{-2}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\biggl( \frac{\partial \lambda_2}{\partial x} \biggr)^2
+ \biggl( \frac{\partial \lambda_2}{\partial y} \biggr)^2
+ \biggl( \frac{\partial \lambda_2}{\partial z} \biggr)^2
</math>
  </td>
</tr>
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\biggl[ \frac{x}{p^2 z^2 \lambda_2} \biggr]^2
+ \biggl[ \frac{q^2y}{p^2 z^2 \lambda_2} \biggr]^2
+ \biggl[ \frac{ \mathfrak{f} }{z \lambda_2 }  - \frac{\lambda_2 }{z}\biggr]^2
</math>
  </td>
</tr>
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\biggl[ \frac{x^2 + q^4 y^2}{p^4 z^4 \lambda_2^2} \biggr]
+ \biggl[ \frac{1}{z\lambda_2}\biggl( \mathfrak{f} - \lambda_2^2 \biggr) \biggr]^2
</math>
  </td>
</tr>
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~ \frac{1}{p^4 z^4 \lambda_2^2}
\biggl[ x^2 + q^4 y^2 + p^4 z^2 (\mathfrak{f} - \lambda_2^2)^2 \biggr]
</math>
  </td>
</tr>
<tr>
  <td align="right">
<math>~\Rightarrow ~~~ h_2</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\frac{p^2 z^2 \lambda_2}{[ x^2 + q^4 y^2 + p^4 z^2 (\mathfrak{f} - \lambda_2^2)^2 ]^{1 / 2} } \, .
</math>
  </td>
</tr>
</table>
So, the associated unit vector is,
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~\hat{e}_2</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\hat{\imath} h_2 \biggl( \frac{\partial \lambda_2}{\partial x} \biggr)
+ \hat{\jmath} h_2 \biggl( \frac{\partial \lambda_2}{\partial y} \biggr)
+ \hat{k} h_2 \biggl( \frac{\partial \lambda_2}{\partial z} \biggr)
</math>
  </td>
</tr>
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\hat{\imath}  \biggl\{ \frac{x}{[ x^2 + q^4 y^2 + p^4 z^2 (\mathfrak{f} - \lambda_2^2)^2 ]^{1 / 2} } \biggr\}
+ \hat{\jmath}  \biggl\{ \frac{q^2 y}{[ x^2 + q^4 y^2 + p^4 z^2 (\mathfrak{f} - \lambda_2^2)^2 ]^{1 / 2} } \biggr\}
+ \hat{k} \biggl\{ \frac{p^2z(\mathfrak{f}-\lambda_2^2)}{[ x^2 + q^4 y^2 + p^4 z^2 (\mathfrak{f} - \lambda_2^2)^2 ]^{1 / 2} } \biggr\} \, .
</math>
  </td>
</table>
Checking orthogonality &hellip;
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~\hat{e}_1 \cdot \hat{e}_2</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
x\ell_{3D} \biggl\{ \frac{x}{[ x^2 + q^4 y^2 + p^4 z^2 (\mathfrak{f} - \lambda_2^2)^2 ]^{1 / 2} } \biggr\}
+ q^2 y\ell_{3D} \biggl\{ \frac{q^2 y}{[ x^2 + q^4 y^2 + p^4 z^2 (\mathfrak{f} - \lambda_2^2)^2 ]^{1 / 2} } \biggr\}
+ p^2 z \ell_{3D} \biggl\{ \frac{p^2z(\mathfrak{f}-\lambda_2^2)}{[ x^2 + q^4 y^2 + p^4 z^2 (\mathfrak{f} - \lambda_2^2)^2 ]^{1 / 2} } \biggr\}
</math>
  </td>
</tr>
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\frac{\ell_{3D}}{ [ x^2 + q^4 y^2 + p^4 z^2 (\mathfrak{f} - \lambda_2^2)^2 ]^{1 / 2} }
\biggl\{ x^2 + q^4y^2 + p^4 z^2 (\mathfrak{f} - \lambda_2^2)\biggr\} \, .
</math>
  </td>
</tr>
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\frac{\ell_{3D}}{ [ x^2 + q^4 y^2 + p^4 z^2 (\mathfrak{f} - \lambda_2^2)^2 ]^{1 / 2} }
\biggl\{ x^2 + q^4y^2 + p^4 z^2 (\mathfrak{f} - \lambda_2^2)\biggr\} \, .
</math>
  </td>
</tr>
</table>
If <math>~\mathfrak{f} = 0</math>, we have &hellip;
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~p^2 z (\mathfrak{f} - \lambda_2^2) </math>
  </td>
  <td align="center">
&nbsp; &nbsp; <math>~~~\rightarrow ~~~ </math>&nbsp; &nbsp;
  </td>
  <td align="left">
<math>~
\biggl[- p^2 z \lambda_2^2\biggr]_{\mathfrak{f} = 0}
=
- p^2 z\biggl[\frac{(x^2 + q^2y^2 + \cancelto{0}{\mathfrak{f} \cdot }p^2 z^2 )^{1 / 2}}{pz}  \biggr]^2
=
- \frac{(x^2 + q^2y^2  )}{z} \, ,</math>
  </td>
</tr>
</table>
which, in turn, means &hellip;
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~[ x^2 + q^4 y^2 + p^4 z^2 (\cancelto{0}{\mathfrak{f}} - \lambda_2^2)^2 ]^{1 / 2}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
[ x^2 + q^4 y^2 + p^4 z^2 \lambda_2^4 ]^{1 / 2}
</math>
  </td>
</tr>
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\biggl\{ x^2 + q^4 y^2 + p^4 z^2 \biggl[\frac{(x^2 + q^2y^2 + \cancelto{0}{\mathfrak{f}\cdot} p^2 z^2)^{1 / 2}}{pz}  \biggr]^4 \biggr\}^{1 / 2}
</math>
  </td>
</tr>
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\biggl\{ x^2 + q^4 y^2 + \biggl[\frac{(x^2 + q^2y^2 )^{2}}{z^2}  \biggr] \biggr\}^{1 / 2}
</math>
  </td>
</tr>
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
(x^2 + q^4 y^2)^{1 / 2} \biggl[ 1 + \frac{(x^2 + q^2y^2 )}{z^2} \biggr]^{1 / 2}
</math>
  </td>
</tr>
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\frac{(x^2 + q^4 y^2)^{1 / 2}}{z} \biggl[ z^2 + (x^2 + q^2y^2 ) \biggr]^{1 / 2} \, ,
</math>
  </td>
</tr>
</table>
and,
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~\hat{e}_1 \cdot \hat{e}_2</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\frac{\ell_{3D}}{ [ x^2 + q^4 y^2 + p^4 z^2 (\mathfrak{f} - \lambda_2^2)^2 ]^{1 / 2} }
\biggl\{ x^2 + q^4y^2 + p^4 z^2 (\mathfrak{f} - \lambda_2^2)\biggr\} \, .
</math>
  </td>
</tr>
</table>


{{ SGFfooter }}
{{ SGFfooter }}

Revision as of 16:37, 23 July 2021

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,

r2x2+y2+z2

=

λ12y2(q21)z2(p21);

tan2θx2+y2z2

=

1z2[λ12y2(q21)p2z2];

1tan2φx2y2

=

λ12p2z2y2q2.

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

=

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

where,

𝒟2

[(1+q2tan2φ)(p2+tan2θ)p2(q21)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),

 

1sin2θ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+p2tan2θ),

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

1p2(z2λ12)=[11(1+p2tan2θ)].

We also see that,

ϖ2p2z2

(1+p2tan2θ)[11(1+p2tan2θ)]=p2tan2θ.

Again Consider Full 3D Ellipsoid

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

𝒟2p2

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

which means that,

p2z2λ12

=

(1+tan2φ)[(1+q2tan2φ)λ2(q21)tan2φ]

 

=

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

q2y2λ12

=

q2tan2φ(1+q2tan2φ)q2tan2φ(1+tan2φ)(1+q2tan2φ)[(1+q2tan2φ)λ2(q21)tan2φ]

 

=

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

 

=

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

 

=

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

x2λ12

=

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

 

 

(1+tan2φ)[(1+q2tan2φ)λ2(q21)tan2φ]

 

=

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

 

=

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

x2+q2y2λ12

=

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

Now, notice that,

q2y2Q2λ12

=

11[]sin2φ,

and,

x2λ12+1Q2

=

11[]sin2φ.

Hence,

x2λ12+1Q2

=

q2y2Q2λ12

0

=

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

 

=

Q4Q2(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

λ1x

=

xλ1,

λ1y

=

q2yλ1,

λ1z

=

p2zλ1.

h12

=

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

=

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

h1

=

λ13D,

where,

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

As a result, the associated unit vector is,

e^1

=

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

 

=

ı^x3D+ȷ^q2y3D+k^p2z3D.

Notice that,

e^1e^1

=

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


Second Coordinate (1st Try)

λ2x

=

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

λ2y

=

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

λ2z

=

(x2+y2)1/2pz2.

h22

=

(λ2x)2+(λ2y)2+(λ2z)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(λ2x)+ȷ^h2(λ2y)+k^h2(λ2z)

 

=

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

Notice that,

e^2e^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^1e^2

=

x3D[xzr(x2+y2)1/2]+q2y3D[yzr(x2+y2)1/2]p2z3D[(x2+y2)1/2r]

 

=

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

 

=

z3Dr(x2+y2)1/2{[x2]+[q2y2][p2(x2+y2)]}

 

0 .


Second Coordinate (2nd Try)

Let's try,

λ2

=

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

λ2x

=

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

λ2y

=

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

λ2z

=

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

Hence,

h22

=

(λ2x)2+(λ2y)2+(λ2z)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(λ2x)+ȷ^h2(λ2y)+k^h2(λ2z)

 

=

ı^{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^1e^2

=

x3D{x[x2+q4y2+p4z2(𝔣λ22)2]1/2}+q2y3D{q2y[x2+q4y2+p4z2(𝔣λ22)2]1/2}+p2z3D{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^1e^2

=

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

Tiled Menu

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