Home / Math and Physics Files / Physics / Mechanics / new frames doc / Unresolved Frames Doc Folders / Hard Questions
Question 5 sun earch doc check RESOLVED
DOCX · 23.3 KB
Open DOCX file
Phil's dated note (PhL 2.28.17) written after he corrected a bad frame change in Section 8.8 (moon-earth) of his frames document. He re-reads the sun-earth (astronomy) document, covering the ecliptic tilt rotation about x, the celestial sphere to local-surface Cartesian transformation by the basis vector method, and the conversion to local spherical angles. He compares the equations with the frames document and concludes the transformations are consistent.
AI-written summary; may contain errors.
Extracted text (machine-read; may contain errors)
How did I change reference frames in sun earth doc? PhL 2.28.17
This is a due diligence test that luckily passed OK.
I just corrected a bad frame change in frames doc Section 8.8 regarding moon earth.
I am now concerned that I might have made this same error elsewhere in my docs. In particular, I am worried about sun earth doc (astronomy doc) where I shift from north pole to local location on the earth. How did I do that frame change in that doc ??? (very nervous right now. ...)
Here is one set of equations
The black equatorial circle is associated with Frame S with x,y,z axes, while the ecliptic red circle is associated with Frame S' with x',y',z' axes. Staring at the figure one sees that
' = Rx(θtmax) and ' = => x' = x
' = Rx(θtmax) (6.1)
where θtmax = 23o. According to the rule that vectors rotate backwards from the unit vectors (see Appendix B), we then have in 2D that
r = Rx(θtmax)r' => = =>
Notice the quote "according to the rule that vectors rotate backwards from the unit vectors". So at least I was aware of this fact at the time, so the last equation is then compatible with the unit vector equations above it. (whew!). I make reference to Appendix B which is my sun-earth take on active versus passive. I will read it right now. // It is just fine! It does not need notation like (V)' because in that doc I don't have vectors like v' which have other definitions. All OK at least at this point in sun earth.
At this point in sun-earth I am using Frame S = celestial sphere and Frame S' = tilted ecliptic as shown in Fig 6.1.
The next Big Deal is Section 7 where I have Frame S = celestial sphere and Frame S' = local system on surface of the earth. So how do I accomplish this transformation?
I first do Cartesian to Cartesian in Section 7 (c) based on this equation
= cosφsinθ + cosφcosθ - sinφ
= sinφsinθ + sinφcosθ + cosφ
= cosθ – sinθ . (7.10)
which I got from a sun-earth appendix. I sure hope this agrees with frames doc on the same thing:
= sinθcosφ + cosθcosφ - sinφ
= sinθsinφ + cosθsinφ + cosφ
= cosθ - sinθ . (E.2.7)
OK that is a very good thing! They agree. Look now at sun earth app A where this is derived. Ouch, I quote it from another of my docs, which is Ref [6] of sun earth which is frames doc!!! But I probably have changed eq nums? I refer to this as Appendix A (A.13.c) of frames doc. YES! That is still correct, I have not changed it! I have this in both App A and App E of frames doc.
I am now reading near (7.10) of sun earth doc. The vector b is earth center to local frame origin on surface of the earth in Fig 7.1, same as in frames doc. I am asking about the sun position in Frame S'.
In (7.1) I have
= e'3 = - e'2 = e'1 , (7.1)
so these are how the Frame S basis vectors like are related to the Frame S' ones like e'3 in Fig 7.1.
I then install these into the triple of equations shown above as E.2.7 to get
= cosφb sinθ1 e'3 – cosφb cosθ1 e'2 – sinφb e'1
= sinφb sinθ1 e'3 – sinφb cosθ1 e'2 + cosφb e'1
= cosθ1 e'3 + sinθ1 e'2 . (7.11)
And then I write
r' ≈ r = x + y + z
= x (cosφbsinθ1 e'3 – cosφbcosθ1 e'2 – sinφb e'1)
+ y (sinφbsinθ1 e'3 – sinφbcosθ1 e'2 + cosφb e'1)
+ z (cosθ1 e'3 + sinθ1 e'2)
= [– x sinφb + y cosφb]e'1
+ [– x cosφbcosθ1 – ysinφbcosθ1 + z sinθ1]e'2
+ [ x cosφbsinθ1 + y sinφbsinθ1 + z cosθ1]e'3
(7.12)
= x' e'1 + y' e'2 + z' e'3
so that the Cartesian sun location coordinates in Frame S' are given by
x' = – x sinφb + y cosφb
y' = – x cosφbcosθ1 – ysinφbcosθ1 + z sinθ1
z' = x cosφbsinθ1 + y sinφbsinθ1+ z cosθ1 . (7.13)
This is "basis vector method" and the result looks good to me. I now have local Cartesian coordinates of the location of the sun. Here (x,y,z) is r in Frame S which is the center of earth aligned frame. I know that the sun's path in this Frame S is given by
x = R cosθt cosφ
y = R cosθt sinφ
z = R sinθt (7.6)
and I then get
x' = R cosθt sin(φ –φb)
y' = R [sinθt sinθ1 – cosθt cosθ1 cos(φ–φb) ]
z' = R [sinθt cosθ1 + cosθt sinθ1 cos(φ-φb)] (7.14)
So these are the sun coordinates in local earth Frame S'. So this concludes my Cartesian → Cartesian transformation, it looks OK.
Then in Section 7 (d) I convert to θ' and φ' as Frame S' spherical coordinates!
x' = r'sinθ'cosφ' y' = r'sinθ'sinφ' z' = r'cosθ'
cosθ' = [sinθt cosθ1 + cosθt sinθ1 cos(φ-φb)]
sinθ' = ρ'/r' = (7.27)
cosφ' = x'/ρ' = [cosθt sin(φ –φb)] /
sinφ' = y'/ρ' = [sinθt sinθ1 – cosθt cosθ1 cos(φ –φb) ]/
tanφ' = y'/x' = [sinθt sinθ1 – cosθt cosθ1 cos(φ –φb) ]/ [cosθt sin(φ –φb)] (7.28)
Therefore, we may write
θ' = cos-1[sinθt cosθ1 + cosθt sinθ1 cos(φ-φb)]
φ' = tan-1 ( [sinθt sinθ1 – cosθt cosθ1 cos(φ –φb) ] / [cosθt sin(φ –φb)] ) . (7.29)
OK enough. I think I can rest easy that I did not do a major screw-up in sun earth on this subject. As noted, I was aware of the back-rotation of basis vectors and forward of regulars. Also, my plot results all came out reasonable so I don't think I could have made such a major error!