Appendix/Ramblings/T6CoordinatesPt3
Concentric Ellipsoidal (T6) Coordinates (Part 3)
Best Thus Far
Part A
| Direction Cosine Components for T6 Coordinates | ||||||||||||||
| --- | --- | --- | --- | --- | ||||||||||
|
||||||||||||||
|
Try …
|
|
|
|
|
|
|
|
|
|
|
|
In this case we find,
|
|
|
|
|
|
|
|
|
|
|
|
The scale factor is, then,
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Part B (25 February 2021)
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
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
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
|
|
|
|
|
|
|
|
COLLADA
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
| First Trial (specified variable values have bgcolor="pink") |
|||||
| x | y | z | |||
| 0.5 | 0.35493 | 0.00000 | 1 | 0.46052 | 2.11310 |
Unit Vectors
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Tangent Plane
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
Corner 1: Shift surface-point location by in the direction, by in the direction, and by by in the direction. This gives …
|
|
|
|
Second Trial
| 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
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
GoodPlane01
| 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
| Tangent Plane Schematic | |
GoodPlane03
| 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
| Tangent Plane Schematic |
Further Exploration
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
|
Appendices: | VisTrailsEquations | VisTrailsVariables | References | Ramblings | VisTrailsImages | myphys.lsu | ADS | |