## Tuesday, September 21, 2021

### Revisiting The Delaunay Variables of Celestial Mechanics

Perturbation diagram to for Delaunay variable analysis

We already saw the Delaunay variables play an important role in celestial mechanics, especially in perturbation theory.  This is concerned with how the presence of a third body (say m 3  in the accompanying diagram) perturbs or disturbs the motion of m 2.  Given the complexity, it is useful to revisit these concepts.  These Delaunay variables are defined by the equations:

ℓ  =   n(t  - T)

L =  Ö(m a)

G  =   L    Ö (1 -   e 2)

H =   G cos i

g   =  w

h  =   Ω

The symbols used are standard, e.g.  w   is the argument of the perihelion, W is the longitude of the ascending node, e is the eccentricity, a is the semi-major axis of the orbit,   and i is the inclination of the orbit.   Also, n is the mean motion (i.e. of the planet or other body per day),  T is the date of perihelion, t is any future date assigned,  and m  is the reduced mass, i.e.

=   1/ (1/m 2   + 1/ m 3)

As with all similar problems in Hamiltonian mechanics, the approach begins with the applicable Hamiltonian for the system (see diagram).  In its most general form we write:

=   ½ å 3 i=1  P i 2    -  m / r  +  R

After modification and as a function of the Delaunay variables (L, ):

H  (L,  ℓ  )  =   -  m 2 / 2 L 2  -  2  3    {1/  D (L,  ℓ  )  +   r ·  r 3   /   r  }

A particular set of differential equations will always be associated with the new Hamiltonian, e.g.

L’ =    -  H  /  ℓ        ℓ    =     H  /  L

g ‘  =   H  /  G          h’ =    H   H

G’ =   -    H  /  g     H’  =  -    H  /  h

A huge simplification arrives  if R = 0., which then reduces to the 2-body problem, which we've examined in previous blog posts.   Note that R in the second Hamiltonian is expressed:

R   =   2  3    {1/  D (L,  ℓ  )  +   r ·  r 3   /   r  }

Now, for the case at hand (consult diagram), we may write:

D 2     =  r  2   +   r  2    -   2 ·  r 3  cos S

And:  r ·  r 3  =   r   r 3  cos S

Assume now that 3  is greater than r (perturbing an inner planet or body by an outer one):

1 /D      =  1/r [ 1 + (r / r 3)  2  -  2 (r / r 3)  cos S] 1/2

Then the above can be rewritten as the sum:

½ å ¥ r= 0  (r / r 3i  P i

=   1/  r 3  [ o   +   (r / r 3)    + (r / r 3)  2     +  ...]

The P i   are functions of the angle S and are called Legendre polynomials.  In this case, the first three may be written:

o   =   1,      P 1     =   cos S,    P 2     =   ½ ( cos 2 S   -   1)