Aerospace & Aeronautical Engineering / Keplerian Orbit Mechanics
FREE TO EXPLORE
LABORATORY 41 · AEROSPACE & AERONAUTICAL ENGINEERING

Keplerian Orbit Mechanics theory

Understand keplerian orbit mechanics: governing aerodynamic equations, shock relations, stability derivatives, and propulsion models.

Aerodynamic and physical foundations of Keplerian Orbit Mechanics.

Orbital mechanics and celestial trajectories are founded upon Johannes Kepler's empirical planetary laws and Sir Isaac Newton's universal law of gravitation. For a spacecraft of mass m orbiting a massive primary body (such as Earth, mass M >> m), the motion is modeled by the fundamental two-body equation: d^2 r / dt^2 = - (mu / r^3) * r, where mu = G * M is the standard gravitational parameter (mu_Earth = 3.986004 x 10^14 m^3/s^2). Integrating this equation yields conservation of specific angular momentum h = r x v and conservation of specific mechanical orbital energy: epsilon = 0.5 * v^2 - mu / r = - mu / (2 * a), where a is the semi-major axis of the conic section orbit. This energy conservation principle gives rise to the celebrated Vis-Viva equation, which explicitly dictates orbital velocity v at any radial distance r from the primary: v^2 = mu * (2 / r - 1 / a) = 2 * mu * (1 / r) - mu / a. For an eccentric Earth orbit of semi-major axis a = 12,000 km, plotting the squared orbital velocity v^2 against the inverse geocentric radius 1 / r demonstrates a strictly linear relationship whose gradient directly extracts twice the gravitational parameter 2 * mu = 7.972 x 10^14 m^3/s^2. In this experiment, simulated deep-space telemetry tracking of an orbital satellite across radial stations verifies the Vis-Viva energy formulation.

v^2 = 2\mu \left(\frac{1}{r}\right) - \frac{\mu}{a} \implies \left(v^2 \times 10^{-7}\right) = m \cdot \left(\frac{1}{r} \times 10^7\right) + b

Viva voce questions & analytical theory

What are Johannes Kepler's Three Laws of Planetary Motion?

1. Law of Ellipses: Every planet moves in an elliptical orbit with the Sun at one of the two foci. 2. Law of Equal Areas: A line segment joining a planet and the Sun sweeps out equal areas during equal intervals of time (conservation of angular momentum). 3. Law of Harmonies: The square of the orbital period T is directly proportional to the cube of the semi-major axis a: T^2 = (4*pi^2 / mu) * a^3.

What is the physical meaning of the Vis-Viva equation in orbital mechanics?

The Vis-Viva equation (Latin for 'living force') v^2 = mu * (2/r - 1/a) is the direct mathematical manifestation of the conservation of mechanical energy in a gravitational field. It equates kinetic energy (0.5 * v^2) and gravitational potential energy (-mu / r) to constant total orbital energy (-mu / (2*a)).

What are the six classical Keplerian orbital elements required to uniquely define an orbit in 3D space?

1. Semi-major axis (a) — defines orbit size; 2. Eccentricity (e) — defines orbit shape; 3. Inclination (i) — defines orbital plane tilt relative to equator; 4. Right Ascension of the Ascending Node (Omega) — defines orbital plane orientation in space; 5. Argument of Periapsis (omega) — defines orientation of ellipse in its plane; and 6. True Anomaly (nu) — defines satellite position along the orbit at time t.

What distinguishes circular, elliptical, parabolic, and hyperbolic orbits in terms of energy and eccentricity?

Circular orbit: e = 0, epsilon < 0 (bound). Elliptical orbit: 0 < e < 1, epsilon < 0 (bound). Parabolic orbit: e = 1, epsilon = 0 (escape trajectory with zero excess velocity). Hyperbolic orbit: e > 1, epsilon > 0 (unbound flyby trajectory with positive hyperbolic excess velocity v_inf).

What is escape velocity and how is it derived from the Vis-Viva equation?

Escape velocity v_esc is the minimum speed required to escape the gravitational pull of a primary body to infinity with zero remaining velocity. Setting semi-major axis a approaches infinity (parabolic trajectory) in the Vis-Viva equation yields v_esc = sqrt(2 * mu / r) = sqrt(2) * v_circular. At Earth's surface, v_esc = 11.186 km/s.

What is a geostationary orbit (GEO) and what orbital radius is required?

A geostationary orbit is a circular, equatorial orbit (i = 0 deg) whose orbital period matches Earth's sidereal rotation period (23 hours, 56 minutes, 4 seconds = 86,164 seconds). Applying Kepler's Third Law yields an orbital radius r_GEO = (mu * T^2 / (4*pi^2))^(1/3) = 42,164 km (altitude 35,786 km). A satellite in GEO appears motionless in the sky.

Experimental precautions and wind tunnel safety

  • Use an Earth-Centered Inertial (ECI J2000) frame; rotating Earth-Centered Earth-Fixed (ECEF) frames introduce fictitious Coriolis accelerations.
  • Account for Earth oblateness (J2 geopotential perturbation) when performing multi-orbit ground track integration.
  • Verify that perigee radius r_p exceeds Earth atmospheric entry radius (r_p > 6500 km) to prevent orbital decay.
  • Maintain double-precision (64-bit IEEE 754) floating-point calculations to avoid numerical orbital energy drift.

Theoretical assumptions & flow limits

  • Pure Keplerian two-body formulation neglects non-spherical Earth gravity harmonics (J2, J3, J4).
  • Third-body gravitational perturbations from the Moon and Sun become noticeable over multi-week timescales.
  • Residual atmospheric drag leeches orbital energy at perigee altitudes below 500 km.
  • Solar radiation pressure exerts minor non-gravitational acceleration on high area-to-mass satellites.
↳
Connect aerodynamic theory to empirical telemetry.

Linear regression of orbital velocity squared against inverse radius yields a slope 2*mu = 7.9720 x 10^14 m^3/s^2 corresponding to Earth mu = 3.9860 x 10^14 m^3/s^2 with R^2 > 0.999999.