Finding Friends
spherical trig
2024-01-08 · 30m

I wrote the last question in the 2023 NAC, called Finding Friends. The premise of the question is that there is an astronomer named Andrew trying to chase a star χ\chi, and the task is to compute some interesting metrics related to whether or not he will ever be able to catch χ\chi.

The original problem was very long, and the final version that ended up being in the exam was condensed by a few parts in order to not make the exam too long or impossible. I figure there is no better place to dump the entire question than my blog, so here is the entire question and solution that I originally submitted.

Question: Finding Friends (solution)

In the future, Earth is perfectly spherical and completely covered with water. Sailor Andrew is on a boat initially positioned at the equator in such that a way that γ\gamma, the vernal equinox, is directly overhead. He observes a star χ\chi with right ascension and declination (α,δ)=(15,10)(\alpha, \delta)=(15^{\circ}, 10^{\circ}). He wants to chase χ\chi, and he needs your help.

You may assume the following numerical values for each constant:

  • RR_{\oplus} is Earth’s radius (same as the current day value).
  • ω=2π/(24 hrs)\omega = 2\pi/(24\text{ hrs}) is Earth’s angular frequency.
  • His boat is very fast; its maximum speed is v=ωR/4\vert v\vert = \omega R_{\oplus}/4.
  • The vernal equinox γ\gamma has celestial coordinates (0,0)(0^{\circ}, 0^{\circ}).
  • (α,δ)=(15,10)(\alpha, \delta) = (15^{\circ}, 10^{\circ}).

Section 1

(a): 2 points (sol)

Calculate the angular distance between γ\gamma and χ\chi. Express your answer in degrees.

In his starting position, χ\chi is somewhere above his horizon. As he sails to different parts of the planet and time passes, χ\chi moves around his night sky. He wants sail his boat in order to get χ\chi as close to his zenith as possible. To make sure he is as efficient as possible, he wants to make a few calculations first.

(b): 2 points (sol)

If he does not move his boat from his starting position, calculate the minimum zenith distance that χ\chi reaches. Express your answer in degrees.

Now, he considers what would happen if he sails straight north at all times.

(c): 3 points (sol)

Starting from his intial position when γ\gamma is overhead, he directs his boat straight north and starts sailing with speed v\vert v\vert. Does χ\chi first cross his north-south meridian or east-west meridian? Briefly explain why this is the case.

(d): 3 points (sol)

Calculate the zenith distance of χ\chi when it first crosses his north-south meridian. Express your answer in degrees.

(e): 15 points (sol)

Calculate the zenith distance of χ\chi when it first crosses his east-west meridian. Express your answer in degrees.

Let XX be the point in his trajectory (in the celestial sphere) at which χ\chi reaches its minimum zenith distance, and NN be the north celestial pole. Note that XX is a point in the celestial sphere, and not Earth.

(f): 6 points (sol)

Argue that NγX+χXNπ/2\angle N\gamma X+\angle \chi XN \approx \pi/2, where both angles are spherical angles. Then, calculate χXN\angle \chi XN. Express your answer in degrees.

Note: For this part only, you may assume that the declination of XX is small, hence the approximation.

(g): 15 points (sol)

Calculate the minimum zenith distance that χ\chi reaches. Express your answer in degrees.

For this problem:

  • You may no longer assume that the declination of XX is small. However, you may use the value of χXN\angle{\chi XN} derived in part (f).
  • You are given that the right ascension of XX is within (α±2)(\alpha\pm 2^{\circ}).
  • His trajectory in the celestial sphere is not necessarily a great circle. Be careful when computing angles and lengths!

Section 2

His next consideration is to see what happens if he steers his ship directly towards χ\chi at all times. As before, assume he starts from his initial position with γ\gamma overhead, and his boat always moves with speed v\vert v\vert. If χ\chi passes below his horizon, he stops his ship.

As before, assume he starts from his initial position when γ\gamma is overhead. Provide a brief argument for why this steering strategy will take him at least one day (i.e., 2π/ω2\pi/\omega) until he can get χ\chi directly overhead using this method.

Relative to the celestial sphere, he has two velocities: one due to the rotation of earth, vrv_r, and one due to his boat’s velocity, vv. Let ζ\zeta be the counterclockwise angle from vrv_{r} to vv. Let θ\theta be the counterclockwise angle from vrv_{r} to vtv_t, where vtv_t is his net velocity in the celestial sphere.

Now we will attempt to characterize how the zenith distance changes over his trajectory. At any particular point in his path, let ϕ\phi be his latitude.

(h): 2 points (sol)

Find an expression for tanθ\tan \theta in terms of ζ\zeta, vv, ω\omega, RR_{\oplus}, and ϕ\phi. Some of these constants have known numerical values, but you do not have to substitute them in.

(Hint: Consider the components of vtv_t. Your final expression should involve all the listed constants.)

(i): 7 points (sol)

In his current trajectory, prove that a small change in the zenith distance of χ\chi, Δz\Delta z, corresponds to a horizontal displacement in the celestial sphere by an angle of

Δλ=Δzcosθcos(ζθ)cosϕ.\Delta \lambda = \frac{\vert \Delta z\vert \cos{\theta}}{\cos{(\zeta - \theta)}\cos \phi}.

Note that Δz\Delta z and Δλ\Delta \lambda are small angular changes.

(j): 6 points (sol)

As before, assume that he starts with γ\gamma overhead and steers his boat towards χ\chi at all times. At some point before χ\chi dips below the horizon for the first time, it reaches a minimum zenith distance. At this moment in time, compute cosϕsinθ\cos{\phi}\sin{\theta}. Your answer should be numerical.

(Hint: you may find your expression from part (h) useful.)

Section 3

He is not satisfied with these results. He feels there should be an optimal way to get χ\chi directly overhead, as fast as possible. Help him figure out the fastest way to get χ\chi to his zenith.

(k): 2 points (sol)

Given that he is allowed to maneuver his boat however he wants throughout the duration of his journey, and that he starts from his initial position with γ\gamma overhead, show that χ\chi cannot be at his zenith the first time it culminates, i.e., crosses his north-south meridian.

(l): 7 points (sol)

Let tt be the fastest time he can get χ\chi to his zenith. Compute the quantity (t20 hrs)(t-\text{20 hrs}), in seconds.

Solution: Finding Friends (question)

Section 1

(a) (ques)

Let NN be the north celestial pole and dd be our desired distance. Then, the spherical law of cosines in Nγχ\triangle{N\gamma\chi} gives

cosd=cos90cos(90δ)+sin90sin(90δ)cosα,\cos{d} = \cos{90^{\circ}}\cos{(90^{\circ}-\delta)} + \sin{90^{\circ}}\sin{(90^{\circ}-\delta)}\cos{\alpha},

which implies d=arccos(cosδcosα)17.96d = \arccos{(\cos{\delta}\cos{\alpha})}\approx 17.96^{\circ}.

(b) (ques)

Since χ\chi is to his east, it has not yet culminated. Also, since he is not moving, χ\chi achieves its minimum zenith distance as it culminates.

One way to figure out this minimum zenith distance is to consider the view from his horizon. χ\chi moves through the small circle representing the declination line at δ\delta. Since he is located on the equator, the celestial equator from his perspective passes through the zenith. Therefore, the minimum zenith distance is the distance between the declination line at δ\delta, and the celestial equator, which is δ=10\delta = 10^{\circ}.

An alternative approach is to consider his movement through the celestial sphere. Culmination occurs when his location and χ\chi have the same right ascension. At this point, their distance is δ\delta, so this is the minimum zenith distance.

(c) (ques)

χ\chi starts in the northeastern quadrant of the sky. When χ\chi crosses the north-south meridian, it changes from the east to west hemisphere. When χ\chi crosses the the east-west meridian, it changes from the north to south hemisphere.

At time tt, his location in the celestial sphere is given by (ωt,ωt/4)(\omega t, \omega t/4), so (α,α/4)(\alpha, \alpha/4) is a point in his trajectory. Since α/4<δ\alpha/4 < \delta, χ\chi culminates in his northern sky. Thus, χ\chi passes the north-south meridian before passing the east-west meridian.

(d) (ques)

χ\chi culminates when it has the same right ascension as Andrew, so ωt=α\omega t = \alpha. So, the zenith distance is given by

δα46.25.\delta - \frac{\alpha}{4} \approx 6.25^{\circ}.

(e) (ques)

Let DD be his position in the celestial sphere. When χ\chi first crosses the east-west meridian, DD must satisfy NDχ=90\angle ND\chi = 90^{\circ}, as shown in the diagram below. The red line in the diagram represents his trajectory in the celestial sphere (this trajectory is not a great circle, so we cannot use it in our calculations).

sphere

Now, chasing angles:

  • The perpendicular condition implies that DD is the apex of another great circle (call Γ\Gamma) that passes through both DD and χ\chi. Γ\Gamma intersects the equator at a point AA, such that DA=90DA = 90^{\circ}.
  • As labelled in the diagram, let BB and CC be the feet of χ\chi and DD from NN, respectively (so that χBA=DCA=90\angle \chi BA = \angle DCA = 90^{\circ}). Since DD is the apex of Γ\Gamma, DA=CA=90DA = CA = 90^{\circ}.
  • We know γB=α\gamma B = \alpha and γC=ωt\gamma C = \omega t, so AB=90+αωtAB = 90^{\circ}+\alpha-\omega t.
  • We also know that χB=δ\chi B = \delta, and DC=ωt/4DC = \omega t/4.

Our goal is to find the zenith distance, which is χD\chi D. To find this arc, we’ll need ωt\omega t, which will let us solve for the right ascension of DD. Let ε=χAB=DAC\varepsilon = \angle{\chi AB} = \angle{DAC}. The four-parts formula in both χAB\triangle{\chi AB} and DAC\triangle{DAC} gives

tanε=tanχBsinAB=tanDCsinAC.\tan{\varepsilon} = \frac{\tan \chi B}{\sin AB} = \frac{\tan DC}{\sin AC}.

Substituting known values, this reduces to the equation

tanδ=cos(ωtα)tan(ωt4)    ωt=4arctan(tanδcos(ωtα)).\tan \delta = \cos{(\omega t - \alpha)}\tan{\left(\frac{\omega t}{4}\right)}\implies \omega t = 4\arctan \left(\frac{\tan \delta}{\cos(\omega t - \alpha)}\right).

We can use this to numerically iterate for ωt\omega t. Our initial guess does not need to be good; any value from around 1010^{\circ} to 7070^{\circ} converges in a few iterations, giving us ωt46.97\omega t\approx 46.97^{\circ}.

Finally, a few more angles:

  • χND=BNC=ωtα\angle \chi ND = \angle BNC = \omega t - \alpha.
  • Nχ=90δN\chi = 90^{\circ}-\delta.
  • ND=90ωt/4ND = 90^{\circ} - \omega t/4.

These are all known values. The spherical law of cosines in NχD\triangle{N\chi D} now gives

χD=arccos(sinδsin(ωt4)+cos(ωtα)cosδcos(ωt4))31.42.\chi D = \arccos{\left(\sin \delta \sin \left(\frac{\omega t}{4}\right) + \cos{(\omega t-\alpha)}\cos \delta \cos{\left(\frac{\omega t}{4}\right)}\right)}\approx 31.42^{\circ}.

(f) (ques)

The zenith distance χX\chi X is minimized when χXγ=90\angle \chi X\gamma = 90^{\circ}. Let XX' be the foot of XX from NN, so that NXγ=90\angle NX'\gamma = 90^{\circ}. Then, if we let β=NXχ\beta = \angle NX\chi:

  • N,X,XN, X, X' lie on the same great circle, NXχ+χXγ+γXX=180\angle NX\chi + \angle \chi X\gamma + \angle \gamma X X' = 180^{\circ}, which implies γXX=90β\angle \gamma XX' = 90-\beta.
  • Since the latitude of XX is small, and XXγ\triangle{XX'\gamma} is right, we can say that γXX+XγX90\angle \gamma XX' + \angle X\gamma X'\approx 90^{\circ}, so XγXβ\angle X\gamma X' \approx \beta. This approximation is valid since γX=ωt=4XX\gamma X' = \omega t = 4X'X, so XXγ\triangle{XX'\gamma} can be approximated by a flat right triangle.
  • Since NγX=90\angle N\gamma X' = 90^{\circ}, β\beta and NγX\angle N\gamma X add to 9090^{\circ}, as desired.

To compute β\beta, we can use our approximation of XXγ\triangle{XX'\gamma} as a flat right triangle to say that tanβ1/4\tan{\beta}\approx 1/4, so β14\beta\approx 14^{\circ}.

(g) (ques)

Consider NXχ\triangle{NX\chi}. We know the lengths of two sides, Nχ=90δN\chi = 90^{\circ}-\delta and NX=90ωt/4NX = 90^{\circ}-\omega t/4, and the values of a nested and adjacent angle, N=ωtα\angle N = \omega t-\alpha and X=β\angle X = \beta. We want to find the length of the zenith angle, which is XχX\chi. Unfortunately, we can’t directly solve for XχX\chi, but we may first try solving for ωt\omega t using the information we already have, which will give us enough information about the triangle to solve for XχX\chi.

Applying the four parts formula:

sin(ωt4)cos(ωtα)=cos(ωt4)cot(90δ)sin(ωtα)cotβ.\sin{\left(\frac{\omega t}{4}\right)}\cos{(\omega t - \alpha)} = \cos{\left(\frac{\omega t}{4}\right)}\cot{(90^{\circ}-\delta)} - \sin{(\omega t - \alpha)}\cot{\beta}.

As in part (e), we can isolate ωt\omega t and try to solve this equation with numerical iteration:

ωt=4arctan(tanδcos(ωtα)sin(ωtα)cotβcos(ωt/4)cos(ωtα)).\omega t = 4\arctan{\left(\frac{\tan \delta}{\cos{(\omega t-\alpha)}} - \frac{\sin{(\omega t- \alpha)}\cot \beta}{\cos{(\omega t/4)}\cos{(\omega t - \alpha)}}\right)}.

Unfortunately, this diverges for most values. Luckily, we’re given that ωt\omega t is within (α±2)(\alpha\pm 2^{\circ}), which is a good enough range for us to still find a solution at ωt16.45\omega t\approx 16.45^{\circ}.

Finally, the law of sines in XχN\triangle{X\chi N} gives

χX=arcsin(cosδsin(ωtα)sinβ)6.06.\chi X = \arcsin{\left(\frac{\cos{\delta}\sin{(\omega t - \alpha)}}{\sin{\beta}}\right)}\approx 6.06^{\circ}.

Section 2

(h) (ques)

The vertical component of vtv_t is given by vsinζv\sin{\zeta}, while the horizontal component is given by vcosζ+ωRcosϕv\cos{\zeta} + \omega R_{\oplus}\cos{\phi}. Therefore,

tanθ=vsinζvcosζ+ωRcosϕ.\tan{\theta} = \frac{v\sin{\zeta}}{v\cos{\zeta} + \omega R_{\oplus}\cos{\phi}}.

(i) (ques)

The key for this part is to draw good diagrams. Below is a larger picture of the trajectory:

sphere

Below we decompose the trajectory:

sphere

In the diagram, xtx_t and xt+1x_{t+1} represent two points in his trajectory that are close together. Note that vv always points at χ\chi, while vrv_r always points in the horizontal direction. From the diagram, we can see how Δz\Delta z and Δλ\Delta \lambda are related:

Δzcos(ζθ)=Δλcosϕcosθ    Δλ=Δzcosθcos(ζθ)cosϕ.\frac{\vert \Delta z\vert}{\cos{(\zeta - \theta)}} = \frac{\Delta \lambda \cos{\phi}}{\cos{\theta}}\implies \Delta \lambda = \frac{\vert \Delta z\vert\cos{\theta}}{\cos{(\zeta - \theta)}\cos \phi}.

(j) (ques)

The zenith distance approaches its minimum value when the vectors vv and vtv_t are perpendicular, since this indicates that it will no longer be possible for him to move closer to χ\chi due to Earth’s rotation.

In other words, this happens when ζθ=90\zeta - \theta = 90^{\circ}. Substituting this into our expression from part (h),

sinθcosθ=vcosθvsinθ+ωRcosϕ,\frac{\sin \theta}{\cos \theta} = \frac{v\cos \theta}{-v\sin \theta + \omega R_{\oplus}\cos{\phi}},

which implies

cosϕsinθ=vωR=14.\cos{\phi}\sin{\theta} = \frac{v}{\omega R_{\oplus}} = \frac{1}{4}.

Remark. Note that the condition ζθ=90\zeta - \theta = 90^{\circ} agrees with the equation derived in part (i). As we approach the minimum zenith angle, the amount of horizontal distance needed to travel in order to change the zenith angle by a fixed Δz\Delta z effectively approaches infinity.

Section 3

(k) (ques)

In order for χ\chi to be directly overhead, he needs to have the same declination and right ascension as χ\chi.

The amount of time he needs to increase his latitude to δ\delta is bounded below by

δ/(v/R)=4δω,\delta/(v/R_{\oplus}) = \frac{4\delta}{\omega},

which is when he travels straight north at all times.

On the other hand, the amount of time it takes him to have the same right ascension is bounded above by

α/(ωv/R)=4α3ω,\alpha/(\omega - v/R_{\oplus}) = \frac{4\alpha}{3\omega},

which is when he travels straight west at all times. Since 4α/3<4δ4\alpha/3 < 4\delta, he cannot get χ\chi to his zenith when it first culminates.

(l) (ques)

From part (h), he cannot get χ\chi directly overhead the first time it culminates. Therefore, his best strategy is to try to get χ\chi to his zenith the next time it appears above his horizon. After time tt, the change in his right ascension due to the rotation of the earth is ωt\omega t. Therefore, to optimize his distance travelled, he should travel along a great circle (defined on the surface of the Earth) that increases his longitude by (2π+αωt)(2\pi + \alpha - \omega t), and increases his latitude by δ\delta.

Let NN be the north celestial pole, and YY be his stopping point, which has coordinates (2π+αωt,δ)(2\pi + \alpha-\omega t, \delta). Note that the coordinates of YY are his longitude and latitude on Earth, and not in the celestial sphere. Also, Yγ=ωt/4Y\gamma = \omega t/4, since the maximum speed his boat can travel is ωR/4\omega R_{\oplus}/4. Now, the spherical law of cosines in NYγ\triangle{NY\gamma} gives

cos(ωt4)=cosδcos(αωt),\cos{\left(\frac{\omega t}{4}\right)} = \cos{\delta}\cos{(\alpha - \omega t)},

so

ωt=4arccos(cosδcos(αωt)).\omega t = 4\arccos(\cos \delta \cos{(\alpha - \omega t)}).

Now we can numerically iterate. Using t=20 hrst = \text{20 hrs} as a first guess, we get ωt=300.2\omega t = 300.2^{\circ}, which gives (t20 hrs)48 s(t-\text{20 hrs})\approx 48\text{ s}. To account for approximation errors, a large range of values in this general vicinity are acceptable.