Stars and Planets

< < Module 2 > >

—Planet Dynamics—

Chris Ormel

Roadmap module 2

Two-body problem

The solutions of the equation of motion between two gravitating bodies.

|

Orbital resonances

Regular, periodic forcing between planetary bodies, when conjuctions take place at similar points in the orbit


Circular restricted Three-body problem

A special case of the three-body problem with two massive bodies on a circular orbit and a third (massless) test particle. The motion is described in the co-moving frame where the gravitating bodies are at rest.


|

Hill's approximation

Applicable when the secondary mass' is small

Tides

Deformation of an extended body by the perturber's gravity, causing energy to be dissipated

  • Integrals of motion
  • Kepler's equation
  • Orbital elements
  • Guiding center approximation
  • Resonance angle
  • Libration
  • Resonance trapping
  • Jacobean energy
  • Lagrange points
  • Roche lobe overflow
  • Hill radius
  • Roche limit
  • Tidal quality factor
  • Tidal orbital decay
  • Tidal locking
  • Tidal heating

Topics

  • Two body problem
  • Three body problem
  • Tides
  • Resonances

Two body problem

read CO 2.3

 

  • It can be cast in terms of the relative motion and total gravitational mass
  • motion is in a plane

The two-body problem for bodies at a mass ratio . In the inertial frame both bodies orbit around the center of mass (cross). relative motion with one of the bodies put at the center. This representation is far more common (and useful) when . However, motion is not in an inertial coordinate frame. The orbit is elliptical, here with an eccentricity of e=0.7.

Two body problem

read CO 2.3

 

  • It can be cast in terms of the relative motion and total gravitational mass
  • motion is in a plane
  • conserved quantities:
    Energy
    Angular momentum
    Eccentricity vector
  • the orbit

    which describes an ellipse with semi-major axis and semi-minor axis .

The two-body problem describes the relative motion of two bodies m1 and m2. Usually, the most massive body is put at the center. True anomaly ν, semi-major axis a, semi-minor axis pericenter P, apocenter A, eccentricity e.

Two body problem

read CO 2.3

 

  • the orbit

    which describes an ellipse with semi-major axis and semi-minor axis .

  • energy, angular momentum in terms of orbital elements:
  • the orbital period

    where n is the mean motion. This is Kepler's law (not to be confused with Kepler's equation [see below], which describes how is related to time ).

The two-body problem describes the relative motion of two bodies m1 and m2. Usually, the most massive body is put at the center. True anomaly ν, semi-major axis a, semi-minor axis pericenter P, apocenter A, eccentricity e.

Kepler's equation

Three anomalies

  • Mean anomaly (time),

    where t is time, n the mean motion, and t0 the time of pericenter passage

  • True anomaly

    indicates the position of the particle in the orbit

  • Eccentric anomaly E

    another indicator for the position of the orbit (see figure) the eccentric and true anomalies are geometrically related:

Kepler's equation relates the eccentric and mean anomalies

Hence, while we can obtain the mean anomaly (or time) as function of , i.e., , no analytic solutions exist for the reverse problem, . This amounts to solving Kepler's equation — a transcendental equation that must be solved numerically.

eccentric anomaly E, true anomaly ν, semi-major axis a, semi-minor axis b, pericenter P, apocenter A, eccentricity e.

Orbital elements

There are 6 orbital elements:

  1. semi-major axis a
  2. eccentricity e
  3. argument of pericenter ω.

    this is the angle with a reference direction. In the planar case the reference is fixed (e.g., the x-axis); in the 3D case, it is the angle with respect to the ascending node Ω.

    Sometimes, this angle is replaced by the longitude of pericenter

  4. true anomaly

    the position of the particle w.r.t. pericenter. Equivalently, the true longitude Another option is to indicate the positon by the mean anomaly M or by the mean longitude which is not a geometric angle.


for 2D problems, these 4 elements suffice.

Orbital elements

There are 6 orbital elements:

  1. semi-major axis a
  2. eccentricity e
  3. argument of pericenter ω.
  4. true anomaly
  5. inclination i: the angle of the orbit w.r.t. the reference plane.
  6. longitude of the ascending node (Ω in the figure →). The angle to the line in the reference plane where the particle intersects the reference planet in the ascending direction

Guiding center approximation

In the guiding center approximation the orbit solution, as function of mean anomaly M reads, to first order in eccentricity:

This is a superposition of:

  1. uniform circular motion (M is linear with time); and
  2. uniform elliptical motion with the ellipse having a 2:1 aspect ratio.

Note. The 2:1 aspect ratio and the epicycle frequency are valid for a Keplerian potential. This can be generalized to arbitrary potentials.

Note. Historically, pre-Copernican Europeans used circular epicyles to describe the motion of planets under the Ptolemaic model (Earth at center). They did not stop with one epicycle!


Guiding center approximation for a particle in an orbit.
motion with respect to the guiding center.

Orbital resonances

The Galilean satellite system is an example of three-body resonance (Laplace resonance)

The period ratios of the planets obey a conmensurability, where the orbital period ratio is a rational number:

  • same for the second and third planet
  • o is called the resonance order
Q: What are j and o for the Galilean system?
"
(c) Wikipedia

resonance or not?

 

"
"
"

Q: Shown are planets at/near a 3:2 mean motion resonance. Inner body is on a circular orbit, outer on an eccentric. Which are in resonance?

resonance or not?

 

close, but not in resonance. The resonance angle is steadily increasing. Conjunctions move away from apocenter
In resonance. Conjunctions are at apocenter. The resonance angle is constant.
In resonance. The orbit precesses, but conjunctions stay at apocenter. The resonance angle is constant