Two-Body Problem

Fundamental orbital motion under ideal gravitational two-body assumptions

3.1 Two-Body Problem

The two-body problem forms the fundamental foundation of orbital mechanics. It describes the motion of two objects interacting only through their mutual gravitational attraction. In many practical situations—such as a satellite orbiting Earth or a planet orbiting the Sun—the mass of one body is much larger than the other. In these cases the larger body can be treated as the primary gravitational source while the smaller body moves in an orbit determined by gravitational physics.

The importance of the two-body model lies in its simplicity and analytical clarity. Under this assumption, the equations of motion can be solved exactly, and the resulting trajectories are well-known geometric curves called conic sections. These solutions provide the baseline description of orbital motion before additional forces such as atmospheric drag, gravitational perturbations from other bodies, or radiation pressure are considered.

Although real space missions experience many disturbances, the two-body solution remains the starting point for almost all orbital analysis, mission design, and trajectory prediction.

3.1.1 Two-Body Assumption

The two-body model relies on several simplifying assumptions about the gravitational interaction between bodies.

Newton’s Universal Gravitation

The gravitational force between two masses is governed by Newton’s law of universal gravitation. This law states that the attractive force between two bodies is proportional to the product of their masses and inversely proportional to the square of the distance separating them. The direction of the force always lies along the line connecting the centers of the two bodies.

This gravitational attraction is responsible for keeping planets in orbit around the Sun, moons around planets, and satellites around Earth.

$$F = G\frac{m_1 m_2}{r^2}$$

Point-Mass or Spherically Symmetric Bodies

In the two-body formulation, each body is treated either as a point mass or as a perfectly spherical mass distribution. Under this assumption, the gravitational field outside the body behaves exactly as if all of the mass were concentrated at its center.

This simplification allows the gravitational interaction to be described using only the distance between the centers of the bodies.

Central Inverse-Square Force

The gravitational force acts as a central force, meaning it always points toward the center of the primary body. Its magnitude follows an inverse-square relationship with distance.

Because the force is central, it produces no torque about the center of mass. As a result, angular momentum of the orbiting body remains constant. This property leads directly to one of the most important features of orbital motion: the motion remains confined to a single plane.

3.1.2 Equation of Motion

Once the gravitational assumptions are defined, the motion of the orbiting body can be described using Newton’s second law.

Orbital Equation in an Inertial Frame

The position of the orbiting body is typically expressed in an inertial coordinate system such as the Earth-Centered Inertial (ECI) frame. In this frame the acceleration of the body is determined entirely by gravitational attraction.

The resulting equation of motion is a second-order differential equation that relates the satellite’s acceleration to its position relative to the central body.

$$\ddot{\mathbf r} = -\frac{\mu}{r^3}\mathbf r$$

Conic-Section Solution

A remarkable property of the two-body equation of motion is that its solution corresponds to conic sections. Depending on the total mechanical energy of the system, the trajectory can take one of several geometric forms:

  • circle
  • ellipse
  • parabola
  • hyperbola

This result shows that gravitational motion naturally produces conic-shaped trajectories, a fact originally observed by Kepler in his study of planetary motion.

Energy and Angular Momentum Constants

Two important quantities remain constant during two-body motion:

Mechanical energy

The total energy of the orbiting body is the sum of its kinetic energy and gravitational potential energy. The sign and magnitude of this energy determine the type of orbit.

$$\varepsilon = \frac{v^2}{2} - \frac{\mu}{r}$$

Angular momentum

Because the gravitational force is central, angular momentum about the central body is conserved. This conservation leads to a key geometric result: the orbiting body sweeps out equal areas in equal time intervals.

$$\mathbf h = \mathbf r \times \mathbf v$$
$$\dot{\mathbf h} = 0$$
$$h = r^2\dot{\theta}$$

Together, these constants allow the orbit to be described analytically and provide a direct connection between orbital geometry and the underlying physics.

Vis-Viva Equation

Rearranging the energy expression gives the fundamental velocity relation

$$v^2 = \mu\left(\frac{2}{r} - \frac{1}{a}\right)$$

where $a$ is the semi-major axis. This equation connects instantaneous velocity with orbital geometry.

3.1.3 Orbital Geometry

The shape and orientation of an orbit are described through several geometric parameters.

Semi-Major Axis and Eccentricity

The semi-major axis defines the overall size of the orbit. For elliptical motion it represents half the longest diameter of the ellipse.

The eccentricity describes the shape of the orbit and indicates how much the trajectory deviates from a perfect circle. Small eccentricity values correspond to nearly circular orbits, while larger values produce increasingly elongated paths.

Periapsis and Apoapsis

The distance between the orbiting body and the central body varies throughout the orbit.

  • Periapsis is the point of closest approach.
  • Apoapsis is the point of greatest distance.

These points define the minimum and maximum orbital radius.

Orbit Classification

The eccentricity parameter determines the class of orbit:

Orbit Type Eccentricity
Circular orbit $e = 0$
Elliptical orbit $0 < e < 1$
Parabolic trajectory $e = 1$
Hyperbolic trajectory $e > 1$

Elliptical orbits correspond to bound motion, while parabolic and hyperbolic trajectories represent escape paths.

General Orbit Equation

The conic-section orbit can be written in polar form as

$$r = \frac{p}{1 + e\cos f}$$

where

  • $p$ is the semi-latus rectum,
  • $e$ is the eccentricity,
  • $f$ is the true anomaly.

Flight-Path Angle

The flight-path angle describes the direction of motion relative to the local radial direction from the central body. It indicates whether the satellite is moving outward, inward, or tangentially along the orbit.

This angle plays an important role in determining the velocity components of the spacecraft and in understanding orbital transfers.

$$\tan\gamma = \frac{e\sin f}{1 + e\cos f}$$

3.1.4 Anomalies and Time

To describe the position of an orbiting body as a function of time, orbital mechanics uses several angular parameters called anomalies.

True Anomaly ($\nu$)

The true anomaly measures the instantaneous angular position of the spacecraft relative to periapsis. It is defined as the angle between the direction of periapsis and the current position vector of the orbiting body.

This angle provides a direct geometric description of the satellite’s location along the orbit.

Eccentric Anomaly ($E$)

Although the true anomaly has clear geometric meaning, it is often more convenient mathematically to use the eccentric anomaly. This angle is defined through a reference circle associated with the ellipse and simplifies many orbital equations.

The eccentric anomaly allows the orbital radius and Cartesian coordinates of the orbiting body to be expressed in compact mathematical form.

$$r = a(1 - e\cos E)$$

Mean Anomaly and Kepler’s Equation

The mean anomaly provides a linear measure of orbital time. It increases uniformly with time and represents the fraction of the orbital period that has elapsed since periapsis passage.

$$M = n(t - t_0)$$
$$n = \sqrt{\frac{\mu}{a^3}}$$

The relationship between mean anomaly and eccentric anomaly is given by Kepler’s equation, which must typically be solved numerically.

$$M = E - e\sin E$$

Together, these anomaly relationships allow the position of a satellite to be computed at any time along its orbit.

3.1.5 Position and Velocity Solutions

Once the orbital geometry and anomaly parameters are known, the satellite’s state can be determined.

State from Anomalies

Using the orbital anomalies and orbital elements, the position of the spacecraft can be calculated relative to the central body. These relationships convert angular position along the orbit into Cartesian coordinates.

Orbital Radius Equation

The distance from the central body varies along the orbit according to the conic-section equation. This equation relates the orbital radius to the true anomaly and the orbit’s geometric parameters.

$$r = \frac{p}{1 + e\cos f}$$

It provides a direct description of how the satellite moves closer to or farther from the central body during its orbit.

Velocity Magnitude and Direction

The spacecraft velocity consists of two components:

  • a radial component directed toward or away from the central body,
  • a transverse component perpendicular to the radial direction.

The magnitude of the velocity depends on both orbital energy and instantaneous position. The direction of motion is determined by the geometry of the orbit and the flight-path angle.

$$v_r = \dot r$$
$$v_\theta = r\dot\theta$$
$$v_r = \frac{\mu}{h}e\sin f$$
$$v_\theta = \frac{\mu}{h}(1 + e\cos f)$$

State Vector in the Orbital Frame

These relationships form the basis for computing the complete orbital state vector, which includes both position and velocity.

$$\mathbf r = r(\cos f\,\hat{\mathbf i}_p + \sin f\,\hat{\mathbf j}_p)$$
$$\mathbf v = \dot r\,\hat{\mathbf i}_p + \frac{h}{r}\,\hat{\mathbf j}_p$$

These expressions define the spacecraft state vector $(\mathbf r,\mathbf v)$, which completely describes the orbit at a given time.