Loading [MathJax]/jax/output/HTML-CSS/jax.js

Monday, March 17, 2014

Improving the estimate of the Sun's rotational speed

To estimate the Sun's rotational speed in the AU lab, we took two points on the Sun's edges and used them to estimate the Sun's rotation speed.  See the picture!



In particular, we took the two points with the greatest difference in speed and assumed they were on a line perpendicular to the Sun's rotation axis, and hence gave a good estimate of the rotation speed.  This is because, for a solid body rotating at some angular frequency ω, the rotation speed of a point on the surface is
v=rω,
with r as shown on the picture below.



Since we are not accounting for any angle dependence in our computation (i.e., we used r=RSun), we have taken it that γ=90.  But in reality it may not be.  Using two pairs of points we can determine its actual value, and thence correct the rotation velocity estimate.  For the pair of points c and d, we have the velocity difference Δv, which we define as Δvij=vivj, to be

Δvcd=2rω=2RSunωsinγ,
while for the pair of points e and f we have
Δvef=2xω=2RSunωcosγ.

Note that x=RSunsin(90γ)=RSuncosγ; you should convince yourself that the angles are really as I've drawn them.  In so doing, you may want to remember that vertical angles are equal.

With our two equations for two different Δvs, we can solve for γ as
γ=tan1[ΔvcdΔvef].

The true rotation velocity difference is then just
Δvtrue=Δvcd1sinγ=Δvcd[Δv2cd+Δv2efΔvcd]=Δv2cd+Δv2ef.

Note the true velocity difference is larger than or equal to our uncorrected result, since sinγ1; the last result also clearly shows this!  To get to the last result, I used the identity that
sin[tan1(XY)]=XX2+Y2.


No comments:

Post a Comment