Orbital Mechanics

Orbital mechanics is the branch of celestial mechanics that answers: where will this object be at time t, given its position and velocity now? The foundational discoveries were Kepler's three laws (1609, 1619) and Newton's derivation of them from universal gravitation (1687). Modern practice adds perturbation theory for multi-body problems and relativistic corrections for high-precision ephemerides. This skill covers the classical core: the six orbital elements, Kepler's laws, the vis-viva equation, orbit types, transfer orbits, gravity assists, and a handful of strategies for computing positions in practice.

Agent affinity: hubble (catalog cross-reference), payne-gaposchkin (binary star orbits and dynamical masses)

Concept IDs: astro-keplers-laws, astro-earth-moon-sun-geometry, astro-planetary-motion

The Two-Body Problem

Two point masses interacting only via gravity produce a relative motion that is exactly solvable. The relative orbit is a conic section — circle, ellipse, parabola, or hyperbola — with the total mass concentrated at one focus.

Assumptions:

• Only two bodies (all others negligible)
• Point masses (or spherically symmetric)
• No radiation pressure, atmosphere, or non-gravitational forces
• Newtonian gravity (relativistic correction needed for Mercury's perihelion, binary pulsars)

When these assumptions hold, the motion is fully determined by six constants — the orbital elements — and an epoch.

The Six Orbital Elements

An orbit in three-dimensional space has six degrees of freedom. There are many choices of six numbers; the classical set is:

For solar-system work the reference plane is usually the ecliptic. For Earth satellites it is the equator. For exoplanets around a host star, the plane is defined by observation geometry.

An alternative sixth element. The true anomaly is time-dependent. For cataloging a fixed orbit you can substitute the time of perihelion passage (T_0) or the mean anomaly at epoch (M_0), and derive nu from them at any later time using Kepler's equation.

Kepler's Three Laws

First law (1609)

Planets move in elliptical orbits with the Sun at one focus.

Consequence: Circular orbits are a special case (e = 0). Most orbits are slightly elliptical — Earth's eccentricity is 0.017. Mercury is 0.206. Pluto is 0.249. Halley's comet is 0.967.

Second law (1609) — Equal areas in equal times

The line from the Sun to a planet sweeps out equal areas in equal time intervals.

Consequence: A planet moves fastest at perihelion (closest to the Sun) and slowest at aphelion. This is equivalent to conservation of angular momentum.

Formula: dA/dt = L / (2m) where L is the orbital angular momentum. For Earth the ratio of perihelion to aphelion speeds is about 1.034, matching the eccentricity.

Third law (1619) — Period and size

The square of the orbital period is proportional to the cube of the semi-major axis:

T^2 = (4 pi^2 / G * M_total) * a^3

For planets around the Sun, using years and AU:

T^2 = a^3

where T is in years and a is in AU. Jupiter has a = 5.2 AU, so T = 5.2^1.5 ~ 11.86 years. Checks.

Binary stars. The same law applies with M_total = M_1 + M_2. Measuring a and T gives the total mass. Measuring the two stars' individual orbits around the barycenter gives the mass ratio, separating the two masses.

The Vis-Viva Equation

The workhorse of orbital mechanics:

v^2 = G * M * (2/r - 1/a)

where v is orbital speed, r is current distance from the focus, and a is the semi-major axis. Vis-viva is Latin for "living force" — an old name for kinetic energy.

Corollaries:

• Circular orbit (r = a): v_circ = sqrt(G M / r). For Earth at 1 AU, this gives 29.78 km/s.
• Escape velocity (a = infinity): v_esc = sqrt(2 G M / r) = v_circ * sqrt(2). For Earth surface, 11.2 km/s.
• Parabolic trajectory (e = 1): v_esc marginally exceeded, approaches zero at infinity.
• Hyperbolic trajectory (a negative in the vis-viva convention): v at infinity is nonzero — you leave and keep going.

Use. Given a spacecraft's position and desired speed, vis-viva tells you what semi-major axis you are on, which fixes period and future positions.

Orbit Types

Circular (e = 0)

Constant speed, constant distance. Low Earth orbit, geostationary orbit, most spacecraft parking orbits. Not strictly realized in nature but a useful idealization.

Elliptical (0 < e < 1)

Closed orbit. All planets, most asteroids, most comets, most moons, most binary stars. The two special points are perihelion (closest to focus) and aphelion (farthest).

For orbits around bodies other than the Sun, the terminology changes: perigee/apogee (Earth), perijove/apojove (Jupiter), periastron/apastron (star), perihelion/aphelion (Sun). The root is the body name.

Parabolic (e = 1)

Marginally unbound. Never repeats. Speed at infinity equals zero. A theoretical boundary case — real "parabolic" comets are typically slightly elliptical or slightly hyperbolic.

Hyperbolic (e > 1)

Unbound. The object visits once and leaves forever. Interstellar objects like 1I/Oumuamua (2017) and 2I/Borisov (2019) follow hyperbolic orbits through the Solar System.

Kepler's Equation

To get position as a function of time in an elliptical orbit, you need to solve Kepler's transcendental equation:

M = E - e * sin(E)

where M is the mean anomaly (linear in time, M = 2 pi * t / T), E is the eccentric anomaly (an intermediate angular parameter), and e is eccentricity.

Solution. Kepler's equation has no closed-form solution. You iterate. Newton's method converges rapidly:

E_{n+1} = E_n - (E_n - e * sin(E_n) - M) / (1 - e * cos(E_n))

Start with E_0 = M for small eccentricities or E_0 = pi for high eccentricities. Typically 3-6 iterations bring it to machine precision.

Once you have E, convert to true anomaly:

tan(nu/2) = sqrt((1+e)/(1-e)) * tan(E/2)

and finally position in the orbital plane is r = a * (1 - e * cos(E)).

Transfer Orbits — Hohmann

The Hohmann transfer is the minimum-energy trajectory between two coplanar circular orbits. Proposed by Walter Hohmann in 1925 — twenty years before anyone could use it.

Procedure. To go from a lower circular orbit of radius r_1 to a higher one of radius r_2:

1. Apply a prograde burn at r_1 that raises apoapsis to r_2. New semi-major axis: a_t = (r_1 + r_2) / 2.
2. Coast on the transfer ellipse half an orbit to apoapsis.
3. Apply a second prograde burn at r_2 that raises periapsis from r_1 to r_2, circularizing.

Delta-v total:

dv_1 = sqrt(G M / r_1) * (sqrt(2 r_2 / (r_1 + r_2)) - 1)
dv_2 = sqrt(G M / r_2) * (1 - sqrt(2 r_1 / (r_1 + r_2)))
dv_total = dv_1 + dv_2

Cost. Earth to Mars Hohmann transfer: about 5.59 km/s total delta-v from low Earth orbit. Transit time: about 259 days.

Trade-offs. Hohmann is energy-optimal but slow. Faster transfers (bi-elliptic, bi-parabolic, or direct high-thrust) spend more propellant for shorter flight time. The Mars launch windows every 26 months arise from the need for Earth and Mars to align for the next Hohmann transfer.

Gravity Assist (Slingshot)

A spacecraft that flies past a moving planet gains or loses heliocentric speed. In the planet's frame, speed in and speed out are equal (elastic encounter); in the Sun's frame, the spacecraft picks up a component of the planet's velocity.

Mathematics. Th