[Physics] Satellite in Elliptical orbit

celestial-mechanicsnewtonian-gravitynewtonian-mechanicsorbital-motion

When finding the period of a satellite orbiting the earth we equate the centripetal force to the gravitational force

$$\frac{mv^2}{r} = \frac{-GMm}{r^2}$$
If I understood well the $r$ cancels into the $r^2$ because the distance from the earth at any point on the orbit equals the radius of the orbit.
Now what if we have an ellipse instead of a circle?

And by the way why does Kepler's $1^{st}$ law say In elliptical orbits

Best Answer

Your equation is an attempt to apply Newton's 2nd Law $F=ma$. The force is central so the acceleration is purely radial. In plane polar co-ordinates the radial component of acceleration has two terms: $\ddot r$ and $-r\dot\theta^2$. The 1st term is zero if the particle is constrained to move in a circle. The 2nd term is the centripetal acceleration.

In your equation, the $\ddot r$ term is missing. You clearly expect elliptical orbits, so $\ddot r \ne 0$. You also have the wrong sign for the $r\dot\theta^2$ term. The centripetal acceleration $r\dot\theta^2$ and the gravitational force $G\frac{Mm}{r^2}$ are both directed inwards, so they should have the same sign.

So the equation of motion which allows for elliptical motion should be
$m(\ddot r - r\dot\theta^2) = -\frac{GMm}{r^2}$.

The variable $\dot\theta$ can be eliminated using the fact that angular momentum is conserved $(mr^2\dot\theta=constant)$ since there is no tangential component of force. You then have a differential equation involving only $r$ and $t$. It does not have a simple solution. Alternatively you can remove the time-dependence and solve for $r$ as a function of $\theta$ using the substitution $u=\frac{1}{r}$. This does have a simple solution, which is the equation of an ellipse. You can find the details of the calculation on many websites, eg Farside Physics.

If you only want to find out how the speed $v$ varies with radius $r$, the solution is much simpler. You apply Conservation of Energy. Kinetic energy is $\frac12mv^2$ and gravitational potential energy is $-\frac{GMm}{r}$, so
$\frac12mv_0^2-\frac{GMm}{r_0}=\frac12mv^2-\frac{GMm}{r}$
Given the speed at any particular radius, you can calculate the speed at any other radius. However, this equation does not tell you the shape of the orbit.

Related Question