Appendix/Ramblings/T6CoordinatesPt3
Concentric Ellipsoidal (T6) Coordinates (Part 3)[edit]
Best Thus Far[edit]
Part A[edit]
| Direction Cosine Components for T6 Coordinates | ||||||||||||||
| --- | --- | --- | --- | --- | ||||||||||
|
||||||||||||||
|
Try …
|
|
|
|
|
|
|
|
|
|
|
|
In this case we find,
|
|
|
|
|
|
|
|
|
|
|
|
The scale factor is, then,
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Part B (25 February 2021)[edit]
Now, from above, we know that,
|
|
|
|
|
Example: |
|||
| 2.14037 | 1.39187 | 0.04623 | 3.57847 |
As an aside, note that,
|
|
|
|
|
|
|
|
We realize that this ratio of lengths may also be written in the form,
|
|
|
|
|
Same Example, but Different Expression: |
||||
| 0.67620 | 0.94359 | 1.87054 | 0.08813 | 3.57847 |
Let's try …
|
|
|
|
|
|
|
|
|
|
|
|
This means that the relevant scale factor is,
|
|
|
|
|
|
|
|
and the three associated direction cosines are,
|
|
|
|
|
|
|
|
|
|
|
|
These direction cosines exactly match what is required in order to ensure that the coordinate, , is everywhere orthogonal to both and . GREAT! The resulting summary table is, therefore:
| Direction Cosine Components for T10 Coordinates | ||||||||||||||
| --- | ||||||||||||||
|
||||||||||||||
Try …
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
This gives,
|
|
|
|
|
|
|
|
Or, given that,
|
|
|
|
we can also write,
|
|
|
|
Similarly,
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Understanding the Volume Element[edit]
Let's see if the expression for the volume element makes sense; that is, does
|
|
|
|
First, let's make sure that we understand how to relate the components of the Cartesian line element with the components of our T10 coordinates.
Line Element[edit]
MF53 claim that the following relation gives the various expressions for the scale factors; we will go ahead and incorporate the expectation that, since our coordinate system is orthogonal, the off-diagonal elements are zero.
|
|
|
|
Let's see. The first term on the RHS is,
|
|
|
|
|
|
|
|
|
|
|
|
the other two terms assume easily deduced, similar forms. When put together and after regrouping terms, we can write,
|
|
|
|
|
|
|
|
|
|
|
|
Given that this summation should also equal the square of the Cartesian line element, , we conclude that the three terms enclosed inside each of the pair of brackets must sum to unity. Specifically, from the coefficient of , we can write,
|
|
|
|
Using this relation to replace in each of the other two bracketed expressions, we find for the coefficients of and , respectively,
|
|
|
|
|
|
|
|
We can use the first of these two expressions to solve for in terms of , namely,
|
|
|
|
|
|
|
|
Analogously, the second of these two expressions gives,
|
|
|
|
Eliminating between the two gives the desired overall expression for , namely,
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
… Not sure this is headed anywhere useful!
Volume Element[edit]
|
|
|
|
|
|
|
|
COLLADA[edit]
Here we try to use the 3D-visualization capabilities of COLLADA to test whether or not the three coordinates associated with the T6 Coordinate system are indeed orthogonal to one another. We begin by making a copy of the Inertial17.dae text file, which we obtain from an accompanying discussion. When viewed with the Mac's Preview application, this group of COLLADA-based instructions displays a purple ellipsoid with axis ratios, (b/a, c/a) = (0.41, 0.385). This means that we are dealing with an ellipsoid for which,
|
|
|
|
and, |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
First Trial[edit]
| First Trial (specified variable values have bgcolor="pink") |
|||||
| x | y | z | |||
| 0.5 | 0.35493 | 0.00000 | 1 | 0.46052 | 2.11310 |
Unit Vectors[edit]
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Tangent Plane[edit]
From our above derivation, the plane that is tangent to the ellipsoid's surface at is given by the expression,
|
|
|
|
For this First Trial, we have (for all values of , given that ) …
|
|
|
|
|
|
|
|
So let's plot a segment of the tangent plane whose four corners are given by the coordinates,
| Corner | x | y | z |
| A | x_0 - 0.25 = +0.25 | 0.41408 | -0.25 |
| B | x_0 + 0.25 = +0.75 | 0.29577 | -0.25 |
| C | x_0 - 0.25 = +0.25 | 0.41408 | +0.25 |
| D | x_0 + 0.25 = +0.75 | 0.29577 | +0.25 |
Now, in order to give some thickness to this tangent-plane, let's adjust the four corner locations by a distance of in the direction.
Eight Corners of Tangent Plane[edit]
Corner 1: Shift surface-point location by in the direction, by in the direction, and by by in the direction. This gives …
|
|
|
|
Second Trial[edit]
| Second Trial … [specified variable values have bgcolor="pink"] |
|||||
| x_0 | y_0 | z_0 | |||
| 0.5 | 0.35493 | 0.00000 | 1 | 0.46052 | 2.11310 |
Generic Unit Vector Expressions[edit]
Let's adopt the notation,
|
|
|
|
for, |
Then, for the T6 Coordinate system, we have,
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
| Second Trial | |||
| x | y | z | |
| 0.23026 | 0.97313 | 0.0 | |
| 0.0 | 0.0 | -1.0 | |
| - 0.97313 | 0.23026 | 0.0 | |
What are the coordinates of the eight corners of a thin tangent-plane? Let's say that we want the plane to extend …
- From to in the direction … here we set ;
- From to in the direction … here we set ;
- From to in the direction … here we set .
|
|
|
|
|
|
|
|
|
|
|
|
| Tangent Plane Schematic |
|
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
†In the figure on the left, vertex 0 is hidden from view behind the (yellow) solid rectangle. |
|||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Third Trial[edit]
GoodPlane01[edit]
| Third Trial … [specified variable values have bgcolor="pink"] |
|||||
| x_0 | y_0 | z_0 | |||
| 0.8 | 0.24600 | 0.00000 | 1 | 0.59959 | 2.34146 |
Again, for the T6 Coordinate system, we have,
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
| Third Trial | ||||
| x | y | z | ||
| 0.47967 | 0.87745 | 0.0 | 0.02 | |
| 0.0 | 0.0 | -1.0 | 0.25 | |
| - 0.87753 | 0.47952 | 0.0 | 0.25 | |
In constructing the Tangent-Plane (TP) for a 3D COLLADA display, we first move from the point that is on the surface of the ellipsoid, , to
| vertex "m" |
Components | |||
| 0 |
|
|
|
|
|
0.8 - 0.02 (0.47952) - 0.25 (-0.87752) = 1.00979 |
0.24590 - 0.02 (0.87753) - 0.25 (0.47952) = 0.10847 |
- 0.25 (-1.0) = + 0.25 |
||
| 1 |
|
|
|
|
|
0.8 - 0.02 (0.47952) - 0.25 (-0.87752) = 1.00979 |
0.24590 - 0.02 (0.87753) - 0.25 (0.47952) = 0.10847 |
+ 0.25 (-1.0) = - 0.25 |
||
| 2 |
|
|
|
|
|
0.8 - 0.02 (0.47952) + 0.25 (-0.87752) = 0.56307 |
0.24590 - 0.02 (0.87753) + 0.25 (0.47952) = 0.34823 |
- 0.25 (-1.0) = + 0.25 |
||
| 3 |
|
|
|
|
|
0.8 - 0.02 (0.47952) + 0.25 (-0.87752) = 0.57103 |
0.24590 - 0.02 (0.87753) + 0.25 (0.47952) = 0.34823 |
+ 0.25 (-1.0) = - 0.25 |
||
| 4 |
0.8 + 0.02 (0.47952) - 0.25 (-0.87752) = 1.0290 |
0.24590 + 0.02 (0.87753) - 0.25 (0.47952) = 0.1436 |
- 0.25 (-1.0) = + 0.25 |
|
| 5 |
0.8 + 0.02 (0.47952) - 0.25 (-0.87752) = 1.0290 |
0.24590 + 0.02 (0.87753) - 0.25 (0.47952) = 0.1436 |
+ 0.25 (-1.0) = - 0.25 |
|
| 6 |
0.8 + 0.02 (0.47952) + 0.25 (-0.87752) = 0.59021 |
0.24590 + 0.02 (0.87753) + 0.25 (0.47952) = 0.38333 |
- 0.25 (-1.0) = + 0.25 |
|
| 7 |
0.8 + 0.02 (0.47952) + 0.25 (-0.87752) = 0.59021 |
0.24590 + 0.02 (0.87753) + 0.25 (0.47952) = 0.38333 |
+ 0.25 (-1.0) = - 0.25 |
|
| Tangent Plane Schematic | Vertex Locations via Excel |
| Tangent Plane Schematic | |
GoodPlane02[edit]
| Tangent Plane Schematic | |
GoodPlane03[edit]
| Tangent Plane Schematic | Tangent Plane Schematic |
|
CAPTION: The image on the right differs from the image on the left in only one way — = 0.1 instead of 0.25. It illustrates more clearly that the (longest) coordinate axis is not parallel to the z-axis when |
|
GoodPlane04[edit]
| Tangent Plane Schematic |
Further Exploration[edit]
Let's set:
|
|
|
|
and, |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Next, let's examine the curve that results from varying while and are held fixed. From the expression for we appreciate that,
|
|
|
|
and from the expression for we have,
|
|
|
|
Hence, the relationship between and is,
|
|
|
|
|
|
|
|
|
Alternatively,
Hence, the relationship between and is,
|
Here are some example values …
| and, | ||||||
| 1st Solution | 2nd Solution | lambda_3 coordinate | ||||
| 0.01 | 0.407825695 | 0.0995168 | - | - | ||
| 0.03 | 0.40481851 | 0.138 | - | - | ||
| 0.04 | 0.40309223 | 0.1503934 | - | - | ||
| 0.08 | 0.393779065 | 0.1854283 | - | - | ||
| 0.12 | 0.37990705 | 0.2103761 | - | - | ||
| 0.16 | 0.36067787 | 0.23111 | 1.04123×10-4 | 0.9095546 | ||
| 0.2 | 0.33500747 | 0.2500033 | 2.23778×10-4 | 0.85448 | ||
| 0.22 | 0.31923525 | 0.2592611 | 3.36653 ×10-4 | 0.82065 | ||
| 0.24 | 0.30106924 | 0.2686685 | 5.2327 ×10-4 | 0.78192 | ||
| 0.26 | 0.2799962 | 0.2784963 | 8.53243 ×10-4 | 0.73752 | ||
| 0.28 | 0.25521147 | 0.2891526 | 1.491545 ×10-3 | 0.68634 | ||
| 0.3 | 0.22530908 | 0.3013752 | 2.89262 ×10-3 | 0.62671 | ||
| 0.32 | 0.1873233 | 0.3168808 | 6.6223 ×10-3 | 0.55579 | ||
| 0.34 | 0.13149897 | 0.3423994 | 2.09221 ×10-2 | 0.46637 | ||
| 0.343 | 0.1191543 | 0.3490285 | 0.026458 | 0.4496 | ||
| 0.344 | 0.1145 | 0.3517 | 0.02880 | 0.4435 | ||
| 0.345 | 0.1093972 | 0.354688 | 0.03155965 | 0.4371186 | ||
| 0.3485 | 0.0847372 | 0.3713588 | 0.0480478 | 0.4085204 | ||
See Also[edit]
|
Appendices: | VisTrailsEquations | VisTrailsVariables | References | Ramblings | VisTrailsImages | myphys.lsu | ADS | |