Chapter 2: Problem 6
Planetary orbits The orbit in space of one body around another, such as a planet around the Sun, need not be circular. In general it takes the form of an ellipse, with the body sometimes closer in and sometimes further out. If you are given the distance \(\ell_{1}\) of closest approach that a planet makes to the Sun, also called its perihelion, and its linear velocity \(v_{1}\) at perihelion, then any other property of the orbit can be calculated from these two as follows. a) Kepler's second law tells us that the distance \(\ell_{2}\) and velocity \(v_{2}\) of the planet at its most distant point, or aphelion, satisfy \(\ell_{2} v_{2}=\ell_{1} v_{1}\). At the same time the total energy, kinetic plus gravitational, of a planet with velocity \(v\) and distance \(r\) from the Sun is given by $$ E=\frac{1}{2} m v^{2}-G \frac{m M}{r}, $$ where \(m\) is the planet's mass, \(M=1.9891 \times 10^{30} \mathrm{~kg}\) is the mass of the Sun, and \(G=6.6738 \times 10^{-11} \mathrm{~m}^{3} \mathrm{~kg}^{-1} \mathrm{~s}^{-2}\) is Newton's gravitational constant. Given that energy must be conserved, show that \(v_{2}\) is the smaller root of the quadratic equation $$ v_{2}^{2}-\frac{2 G M}{v_{1} \ell_{1}} v_{2}-\left[v_{1}^{2}-\frac{2 G M}{\ell_{1}}\right]=0 . $$ Once we have \(v_{2}\) we can calculate \(\ell_{2}\) using the relation \(\ell_{2}=\ell_{1} v_{1} / v_{2}\). b) Given the values of \(v_{1}, \ell_{1}\), and \(\ell_{2}\), other parameters of the orbit are given by simple formulas can that be derived from Kepler's laws and the fact that the orbit is an ellipse: $$ \begin{aligned} \text { Semi-major axis: } & a=\frac{1}{2}\left(\ell_{1}+\ell_{2}\right), \\ \text { Semi-minor axis: } & b=\sqrt{\ell_{1} \ell_{2}}, \\ \text { Orbital period: } & T=\frac{2 \pi a b}{\ell_{1} v_{1}}, \\ \text { Orbital eccentricity: } & e=\frac{\ell_{2}-\ell_{1}}{\ell_{2}+\ell_{1}} . \end{aligned} $$ Write a program that asks the user to enter the distance to the Sun and velocity at perihelion, then calculates and prints the quantities \(\ell_{2}, v_{2}, T\), and \(e\). c) Test your program by having it calculate the properties of the orbits of the Earth (for which \(\ell_{1}=1.4710 \times 10^{11} \mathrm{~m}\) and \(v_{1}=3.0287 \times 10^{4} \mathrm{~ms}^{-1}\) ) and Halley's comet \(\left(\ell_{1}=8.7830 \times 10^{10} \mathrm{~m}\right.\) and \(\left.v_{1}=5.4529 \times 10^{4} \mathrm{~m} \mathrm{~s}^{-1}\right)\). Among other things, you should find that the orbital period of the Earth is one year and that of Halley's comet is about 76 years.
Short Answer
Step by step solution
Key Concepts
These are the key concepts you need to understand to accurately answer the question.