next up previous
Next: The Laplace vector Up: Perturbations and more on Previous: General formulae for perturbations

Perturbation of the energy and the angular momentum

The energy for the Kepler motion is

\begin{displaymath}
E = \frac{1}{2}({\bf v}\cdot {\bf v}) -\frac{\mu}{r}
\end{displaymath}

Then,

\begin{displaymath}
\dot{E} =
\frac{1}{2}(\dot{{\bf v}}\cdot{\bf v}+{\bf v}\cdot\dot{{\bf v}})+\frac{\mu}{r^2}\dot{r}
\end{displaymath}

or

\begin{displaymath}
\dot{E} = \frac{1}{2}(\dot{{\bf v}}\cdot{\bf v})+\frac{\mu}{r^2}\dot{r}
\end{displaymath}

After substituting the equations of motion for $\dot{{\bf v}}$ and rearranging terms, we get

\begin{displaymath}
\dot{E} = \left( -\frac{\mu}{r^3}{\bf r} + {\bf F}\right)\cdot{\bf v}
+ \frac{\mu}{r^2}\dot{r}
\end{displaymath}

We need now $r\dot{r} = {\bf r}\cdot {\bf v}$, and we get

\begin{displaymath}
\dot{E} = \left( -\frac{\mu}{r^r}{\bf r} + {\bf F}\right)\cdot{\bf v}
+ \frac{\mu}{r^3}({\bf r}\cdot {\bf v})
\end{displaymath}

which simplifies to

\begin{displaymath}
\dot{E} = {\bf F}\cdot {\bf v}
\end{displaymath}

This is what we expected, i.e. the change of the energy is proportional to the perturbation. In particular, since ${\bf v}$ is tangential to the orbit, only the component of ${\bf F}$ tangential to the orbit changes the energy of the orbit. Also note that when ${\bf F}= 0$, $E$ is conserved. Conserved quantities are often called first inetrals of the motion.

Next, recall that the angular momentum is defined as

\begin{displaymath}
{\bf A}= {\bf r}\times {\bf v},
\end{displaymath}

and was also a conserved quantity (or first integral) of the unperturbed Kepler motion. To investigate the effect of the perturbation we compute its time derivative

\begin{displaymath}
\dot{{\bf A}} = \dot{{\bf r}} \times {\bf v}+ {{\bf r}} \times \dot{{\bf v}}
\end{displaymath}

which after substituting $\dot{\bf v}$ again yields

\begin{displaymath}
\dot{{\bf A}} = {\bf v}\times {\bf v}+ {{\bf r}} \times \left(
-\frac{\mu}{r^3}{\bf r} + {\bf F}\right).
\end{displaymath}

Here $ {\bf v}\times {\bf v}= 0$ and $ {\bf r}\times {\bf r}= 0$, leading to

\begin{displaymath}
\dot{{\bf A}} = {\bf r}\times {\bf F}
\end{displaymath}

Note that only components of ${\bf F}$ which are perpendicular to ${\bf r}$ contribute to the change of the angular momentum. The quantity ${\bf r}\times{\bf F}$ is known as the torque in mechanics.


next up previous
Next: The Laplace vector Up: Perturbations and more on Previous: General formulae for perturbations
Werner Horn 2006-06-06