This preview shows pages 1–3. Sign up to view the full content.
This preview has intentionally blurred sections. Sign up to view the full version.
View Full Document
Unformatted text preview: Physics 341: Problem Set #3 Solutions 1. For this problem you can use my orbitps03 spreadsheet, or you can write your own program. You can download an Excel version of the spreadsheet from http://www.physics.rutgers.edu/ugrad/341/orbitps03.xls To use the spreadsheet, you will need to edit the fields that are shaded yellow. (a) For this problem we will use natural units for the solar system, measuring lengths in AU and time in years. Use Newtons precise version of Keplers Third Law to show that GM = 4 2 AU 3 yr 2 . This value is fixed in the spreadsheet (cell B3). Newton showed that Keplers Third Law is P 2 = 4 2 GM a 3 = GM = 4 2 a 3 P 2 By definition, the Earth goes around the Sun with a semimajor axis a = 1 AU in a period of P = 1 yr. And so GM = 4 2 (1 AU) 3 (1 yr) 2 = 4 2 AU 3 yr 2 39 . 4784 AU 3 yr 2 (b) We will compute the orbit of Eris, the infamous tenth planet that caused Pluto to be demoted to dwarf planet status. Eris (also called 2003 UB 313 in the textbook) has a semimajor axis a = 68 . 048 AU and eccentricity e = 0 . 4336 based on the best current observations. Calculate the orbital period P of Eris (in years). To determine the period we can use Keplers Third Law: P = 2 a 3 / 2 ( GM ) 1 / 2 Recall that because of our (good) choice of units GM = 4 2 AU 3 yr 2 , so P = 2 a 3 / 2 (4 2 AU 3 yr 2 ) 1 / 2 = (68 . 048 AU) 3 / 2 1 AU 3 / 2 yr 1 = 561 . 33 yr (c) We can choose the initial conditions to have time t = 0 and angle = 0 , with coordinates centered on the Sun. Lets start at perihelion (closest approach to the Sun), so that v r = 0 at t = 0 . Determine the remaining initial conditions ( ` and r ) you need to reproduce the orbit. From class, the specific angular momentum is ` = p GMa (1 e 2 ) = (39 . 4784 AU 3 yr 2 ) (68 . 048 AU) (1 . 4336 2 ) 1 / 2 = 46 . 705 AU 2 yr 1 1 Also from the class notes (see the ellipse diagram in Figure 2 of the notes from Lecture 4), the perihelion distance is r = r p = a ae = a (1 e ) = 68 . 048 AU (1 . 4336) = 38 . 542 AU (d) Plug these initial conditions into the spreadsheet (or your own program) and plot Eriss orbit. You will need to adjust the time step t , to cover one full period. My spreadsheet computes about 1000 time steps, so to cover one full period we should adjust the time step to t = P/ 1000 = 0 . 56 yr. I will use t = 0 . 57 yr just to make sure we cover slightly more than one orbit. Then plugging in our results from part (c) yields the following orbit: !"#$% ( )"$* !"#$ !%$ !&$ $ &$ %$ "#$ !"#$ !%$ !&$ $ &$ %$ "#$ + ,./ 0,./ (e) Write down formulas for the perihelion distance, aphelion distance, and semimi nor axis in terms of just a and e . Calculate these quantities for Eris and compare them to the values in your spreadsheet orbit. How well do they agree?...
View
Full
Document
This note was uploaded on 10/18/2011 for the course PHYSICS 341 taught by Professor Keeton during the Fall '08 term at Rutgers.
 Fall '08
 Keeton
 Physics

Click to edit the document details