next up previous
Next: Perihelion precession rate Up: Perihelion precession of Mercury Previous: Introduction

Time evolution of eccentricity vector

As we saw in Exercise 3, Mercury's eccentricity vector

$\displaystyle {\bf e} = \frac{\skew{2}\dot{\bf r}\times {\bf h}}{\mu} - \frac{{\bf r}}{r},$ (B.3)

where $ {\bf h} = {\bf r}\times\skew{2}\dot{\bf r}$ is its orbital angular momentum per unit mass, is of magnitude $ e$ , and is always directed from the Sun to Mercury's perihelion. Here, $ e$ is the Mercury's orbital eccentricity. It follows that

$\displaystyle \skew{1}\dot{\bf e} = \frac{\skew{2}\ddot{\bf r}\times {\bf h}}{\...
...2}\dot{\bf r}}{r} + \frac{({\bf r}\cdot\skew{2}\dot{\bf r})\,{\bf r}}{r^{\,3}}.$ (B.4)

Now, $ \dot{\bf h} = {\bf r}\times \skew{2}\ddot{\bf r}= {\bf r}\times \nabla{\cal R}$ , where use has been made of Equation (B.1). Hence, making further use of Equation (B.1), the previous expression reduces to

$\displaystyle \skew{1}\dot{\bf e} = \frac{2\,(\skew{2}\dot{\bf r}\cdot\nabla{\c...
...f r})\,{\nabla\cal R} -({\bf r}\cdot\nabla{\cal R})\,\skew{2}\dot{\bf r}}{\mu}.$ (B.5)

We conclude that the disturbing function causes Mercury's eccentricity vector, $ {\bf e}$ , to evolve in time. Any change in length of this vector gives rise to a change in Mercury's orbital eccentricity, whereas any change in direction produces a rotation of Mercury's perihelion. In the following, because we are only interested in the motion of Mercury's perihelion, we shall focus solely on the change in direction (i.e., the precession) of the $ {\bf e}$ vector.


next up previous
Next: Perihelion precession rate Up: Perihelion precession of Mercury Previous: Introduction
Richard Fitzpatrick 2016-03-31