Wednesday, November 6, 2013

Math Drives Astronomy (Pt. 2)



















By the time the aspiring Astronomy major reaches his sophomore year, he will be facing even more math in courses such as Spherical Astronomy.  The emphasis here is on astronomical time (sidereal time and position - the latter based on different coordinate systems used (equatorial, horizontal, ecliptic).

Each such coordinate system is defined by a different set of poles and equator. For example, the coordinate system depicted above - the equatorial system - is based on the projection of the Earth's N. and S. poles into the sky, as well as its equator. The poles then become the North and South Celestial poles, and the equator becomes the celestial equator. If these poles are defined respectively at +90 degrees (NCP) and -90 degrees (SCP) and the celestial equator at 0 degrees, then a consistent system of celestial latitude and longitude  can be defined in a consistent system, to locate any celestial object. We call the longitude coordinate (Θ) the Right Ascension (R.A.) while the celestial latitude coordinate is called declination and is measured in degrees north or south of the celestial equator. (In the diagram, the complementary angle of the declination is shown, φ, which we call the zenith distance.)

One of the first things the astronomy sophomore learns is how to find directions around the celestial sphere, including how to relate the R.A. to time and time keeping.  He must also do exercises showing he can find the hour angle - and distinguish it from the Right Ascension. This starts with identifying the R.A. of his local meridian (the imaginary celestial longitude passing through his zenith or highest point.)  The Right Ascension of the star is clearly equal to the local sidereal time (L.S.T.) plus the hour angle. Thus, we can write:  HA = RA of observer meridian - RA of object
Of course, before he gets very far, he will be adept at sketching any number of diagrams to solve time and position problems, such as the diagram below - with perspective looking down onto the North Celestial pole:

For example, in the diagram shown (Fig. 2) if the star's R.A. is 7h 00m and the observer is at a local sidereal time of 6h 00m, then the hour angle becomes: HA = 6h - 7h = -1h or -15 degrees. (Since every hour of longitude corresponds to 15 degrees angle, i.e. Earth turns through 15 degrees every hour, 360 degrees in 24 hours.)

The budding astronomy student will also have to know how to sketch a three dimensional diagram to enable conversion between coordinates, say from the horizontal (observer -based) to the celestial sphere. This will always include what is called the fundamental "astronomical triangle" (shaded region of Fig. 3) from which spherical trig relationships can be obtained and conversions can be made to different coordinate systems.

Using spherical trig, the student can write the law of sines and law of cosines for spherical triangles (such as shown in Fig. 3)which are analogs of the law of sines and cosines for triangles in plane trig.

We have for the law of sines:

Sin A/ sin a = sin B/ sin b = sin C/ sin c

where A, B, C denote ANGLES and a,b,c denote measured arcs. (Note: we could also have written these by flipping the numerators and denominators).

We have for the law of cosines:

cos a = cos b cos c + sin b sin c cos A

Where a, b, c have the same meanings, and of course, we could write the same relationship out for any included angle.

Now, we use Fig. 3, for a celestial sphere application, in which we use the spherical trig relations to obtain an astronomical measurement.

Using the angles shown in Fig. 3 each of the angles for the law of cosines (given above) can be found. They are as follows:

cos a = cos (90 deg - decl.)

where decl. = declination

cos b = cos (90 deg - Lat)

where 'Lat' denotes the latitude. (Recall from Fig. 1 if φ is polar distance (which can also be zenith distance) then φ = (90 - Lat))

cos c = cos z

where z here is the zenith distance.

sin b = sin (90 deg - Lat)

sin c = sin z

and finally,

cos A = cos A

where A is azimuth.

Let's say we want to find the declination of the star if the observer's latitude is 45 degrees N, and the azimuth of the star is measured to be 60 degrees, with its zenith distance z = 30 degrees. Then one would solve for cos a:

cos a = cos (90 deg - decl.)=

cos (90 deg - Lat) cos z + sin (90 deg - Lat) sin z cos (A)

cos (90 deg - decl.) =
 
cos (90 - 45) cos 30 + sin (90 - 45) sin 30 cos 60
 
And:

cos (90 deg - decl.)= cos (45) cos 30 + sin (45) sin 30 cos 60

We know, or can use tables or calculator to find:

cos 45 = [2]^½/ 2

cos 30 = [3]^½/ 2

sin 45 = [2]^½/ 2

sin 30 = ½

cos 60 = ½

Then:
cos (90 deg - decl.)= {([2]^½/ 2 )([3]^½/ 2)} + {[2]^½/ 2} (½) (½)

cos (90 deg - decl.)= [6]^½/ 4 + [2]^½/ 8
 
= {2[6]^½ + [2]^½}/ 8

cos (90 deg - decl.)= 0.789

arc cos (90 deg - decl.)= 37.9 deg

Then:

decl. = 90 deg - 37.9 deg = 52.1 deg

Or,  decl. (star) = + 52.1 degrees
 
The student will also become familiar with matrix methods for converting between the different coordinate systems.
 
The basic principle involves relating the Cartesian coordinates (rectilinear) of a point on the celestial sphere (diagram) to the curvilinear coordinates measured in the primary and secondary reference planes. One has then, for example:

(x)
(y)
(z) u,v =


(cos v .....cos u)
(cos v .....sin u)
(sin v..............)

After conversion the curvilinear coordinates may be calculated according to:

u = arctan (y/x) and v = arcsin (z)

Now, consider conventional orthogonal matrices of 3 x 3 dimensions, given as functions: R1(Θ), R2(Θ) and R3(Θ), to rotate the general system by the angle Θ about axes x, y and z, respectively. Thus we obtain:

R1(Θ) =

(1..........0................0)
(0.....cos(Θ)..... sin(Θ))
(0......-sin(Θ)....cos(Θ))


R2(Θ) =

(cos (Θ)......0........- sin(Θ))
(0................1...............0.. )
(sin(Θ)........0.......cos(Θ) )


R3(Θ) =


(cos(Θ)..........sin(Θ)..........0)
(-sin (Θ)......cos(Θ)...........0)
(0 ..................0..................1)


To fix ideas, say you wish to obtain the horizontal coordinates (A, a) for some object in the sky and you know its R.A. and decl. from an almanac.  Then the procedure is fairly straightforward, and entails writing:
 
R3(Θ) = R3(-180 deg)

R2(Θ) = R2(90 - lat.)

so that:

(x)
(y)
(z) A,a = R3(-180 deg) R2(90 - lat.) (XYZ(h, decl.))
 
where : (XYZ(h, decl.)) =

(x)
(y)
(z) h,decl.
 
Bear in mind: R3(-180 deg) =

(cos 180........-sin 180................0)
(sin 180........cos 180................0)
(0......................0......................1)


And:R3(Θ) =


(cos(Θ)..........sin(Θ)..........0)
(-sin (Θ)......cos(Θ)...........0)
(0 ..................0..................1)

Therefore: R3(-180) =

(-1.....0.......0)
(0.......-1.....0)
(0.......0.......1)

Obtaining R2(Θ) = R2(90 - lat.) is just as easy, if one recalls the basic trig identity:

cos (90 - φ) = sin (φ)
 
Not surprisingly, it is the sophomore year at which point most wannabe Astronomy majors change their minds and either drop out entirely or change their major - say to a less mathematically demanding subject. At the Univ. of South Florida, nearly 60% had dropped out by the end of the quarter which featured spherical astronomy.
 

No comments:

Post a Comment