Motion of Orbiting Satellites - College Physics

SECTION 1

Historical Context & Motivation

The motion of celestial bodies has captivated human inquiry for millennia, but a rigorous mechanical understanding of orbits emerged only through centuries of painstaking observation and theoretical synthesis. The ancient Greeks proposed geocentric models with complex epicycles to account for the apparent retrograde motion of planets, yet these constructions remained largely kinematic descriptions rather than dynamical explanations. The decisive shift began in the sixteenth and seventeenth centuries, when astronomers and natural philosophers sought to uncover the force law that produces orbital trajectories. Today, satellite mechanics sits at the intersection of gravitational physics, energy conservation, and angular momentum — the very core of rotational mechanics.

1609 • Kepler's First Two Laws

Johannes Kepler published Astronomia Nova, demonstrating that Mars follows an elliptical orbit with the Sun at one focus and that a radius vector sweeps equal areas in equal times — the first quantitative orbital laws.

1687 • Newton's Principia

Isaac Newton's Philosophiæ Naturalis Principia Mathematica derived Kepler's laws from a universal inverse-square gravitational force, unifying terrestrial and celestial mechanics under a single framework.

1957 • Sputnik — First Artificial Satellite

The Soviet Union launched Sputnik 1 into low Earth orbit, validating centuries of orbital theory and inaugurating the space age. Its 96-minute period matched Newtonian predictions to remarkable accuracy.

1963 • Geosynchronous Orbit Realized

Syncom 2 achieved geosynchronous orbit at approximately 35,786 km altitude, confirming the specific orbital radius predicted by equating gravitational and centripetal forces for a 24-hour period — a concept first proposed by Arthur C. Clarke in 1945.

The central question that orbital mechanics addresses is deceptively simple: given a gravitational field, what determines the shape, speed, and stability of an orbit? The answer weaves together Newton's second law, the conservation of mechanical energy, and the conservation of angular momentum. Understanding these principles not only explains the motion of planets and moons but also underpins the design of every communications satellite, GPS constellation, and interplanetary probe in operation today.

SECTION 2

Core Principles & Definitions

Satellite orbital mechanics rests on a small number of powerful ideas that, taken together, fully characterize the motion of a body under a central gravitational force. These principles apply equally to a spacecraft in low Earth orbit and to the Moon's path around the Earth, because the underlying physics is scale-invariant. The key is recognizing that gravity provides the centripetal acceleration required for circular (or, more generally, elliptical) motion, and that two conservation laws — energy and angular momentum — constrain the orbit's geometry and the satellite's speed at every point along its trajectory.

Gravitational Force as Centripetal Force

For a satellite of mass m in a circular orbit of radius r, the gravitational attraction of the central body (mass M) supplies exactly the centripetal force: GMm/r² = mv²/r. This balance uniquely determines the orbital speed for a given radius.

Conservation of Angular Momentum

Because gravity is a central force (directed along the line connecting the two bodies), it exerts zero torque about the center. Hence L = mvr sin θ remains constant throughout the orbit. In an ellipse, the satellite speeds up at periapsis and slows at apoapsis to conserve L.

Conservation of Mechanical Energy

The total mechanical energy E = ½mv² − GMm/r is conserved because gravity is conservative. For bound (closed) orbits, E < 0. The total energy uniquely determines the semi-major axis of the orbit via E = −GMm/(2a).

Kepler's Third Law

The square of the orbital period is proportional to the cube of the semi-major axis: T² = (4π²/GM)a³. This emerges directly from combining the centripetal condition with the definition of period and generalizes naturally from circular to elliptical orbits.

Weightlessness in Orbit

An orbiting satellite and its occupants are in continuous free fall. Gravity has not vanished — indeed it is nearly as strong in low Earth orbit as on the surface — but both the satellite and everything inside it share the same gravitational acceleration, producing the sensation of weightlessness.

✦

KEY TAKEAWAY

Think of an orbiting satellite as a ball rolling around the inside of a frictionless bowl: at every point the bowl's curved surface pushes the ball inward just enough to keep it on its circular path. In orbit, gravity replaces the bowl, providing exactly the inward (centripetal) acceleration needed. If the ball moves faster, it rides higher on the bowl; if slower, it drops lower. The same trade-off between speed and radius governs real orbits, constrained by the conservation of energy and angular momentum.

SECTION 3

Visual Explanation — Circular & Elliptical Orbits

The dashed cyan circle represents a circular orbit at constant radius, with the gravitational force vector (cyan arrow) always directed radially inward and the velocity vector (amber arrow) always tangent. The solid pink ellipse shows an elliptical orbit with the central body at one focus. At periapsis the satellite is closest and fastest; at apoapsis it is farthest and slowest — a direct consequence of angular momentum conservation.

The diagram above illustrates the two canonical orbit types encountered in introductory mechanics. For the circular orbit (dashed cyan), the satellite maintains a fixed distance from the central body and moves with a constant speed determined entirely by the orbital radius. The gravitational force always points radially inward, perpendicular to the velocity, so it changes the direction of motion without altering the speed — the hallmark of uniform circular motion. For the elliptical orbit (pink), the radial distance varies continuously. Kepler's second law (equal areas in equal times) follows directly from the constancy of angular momentum: when r decreases, v must increase so that the product mvr (for the perpendicular component) remains unchanged. Energy conservation simultaneously accounts for the speed variation: kinetic energy increases as the satellite falls deeper into the gravitational potential well, and vice versa.

SECTION 4

Mathematical Framework

We now formalize the principles introduced qualitatively. Throughout this section, M denotes the mass of the central body, m the satellite mass, G the universal gravitational constant, r the orbital radius (center-to-center distance), and v the orbital speed. We assume m ≪ M so that the central body can be treated as stationary.

Orbital Speed for a Circular Orbit

CENTRIPETAL CONDITION

GMm / r² = mv² / r ⟹ v = √(GM / r)

Setting the gravitational force equal to the centripetal force and solving for v. Note that the satellite mass m cancels — the orbital speed depends only on the central body's mass and the orbital radius.

Gravitational Potential Energy & Total Energy

TOTAL MECHANICAL ENERGY (CIRCULAR ORBIT)

E = K + U = ½mv² − GMm/r = −GMm / (2r)

Substituting v² = GM/r into the kinetic energy gives K = GMm/(2r). Since U = −GMm/r, the total energy is negative for all bound orbits, signifying that the satellite is gravitationally trapped. The kinetic energy equals half the magnitude of the potential energy — a result known as the virial theorem for a 1/r² force.

Kepler's Third Law

PERIOD–RADIUS RELATION

T = 2πr / v = 2πr / √(GM/r) = 2π√(r³ / GM) ⟹ T² = (4π² / GM) r³

For circular orbits the period scales as r^3/2. The same relation holds for elliptical orbits with r replaced by the semi-major axis a.

Angular Momentum

ORBITAL ANGULAR MOMENTUM

L = mvr (circular) or L = mvr sinθ (general)

For a circular orbit, v is always perpendicular to r, so sinθ = 1. Because gravity exerts zero torque about the center of force, L is conserved. Substituting v = √(GM/r) gives L = m√(GMr), showing that higher orbits carry more angular momentum per unit mass.

SECTION 5

Energy Diagram & Orbit Classification

A powerful way to classify orbital trajectories is through the effective potential energy diagram. By introducing an effective potential U_eff(r) = −GMm/r + L²/(2mr²), we reduce the two-dimensional orbital problem to an equivalent one-dimensional radial problem. The term L²/(2mr²) acts as a centrifugal barrier that prevents the satellite from collapsing to r = 0 (provided L ≠ 0). The shape of U_eff determines whether an orbit is bound or unbound, circular or elliptical.

The effective potential U_eff (solid purple curve) results from superposing the attractive gravitational potential (dashed cyan) with the repulsive centrifugal term. The minimum of U_eff (amber dot) corresponds to a stable circular orbit. A total energy between U_eff,min and zero (pink dashed line) gives a bound elliptical orbit that oscillates between r_min and r_max.

Orbit Type	Total Energy E	Eccentricity e	Trajectory
Circular	E = U_eff,min < 0	e = 0	Closed, constant radius
Elliptical	U_eff,min < E < 0	0 < e < 1	Closed, varying radius
Parabolic	E = 0	e = 1	Open, just barely escapes
Hyperbolic	E > 0	e > 1	Open, unbound flyby

Classification of Keplerian orbits by total energy and eccentricity

SECTION 6

Worked Example — Geostationary Orbit

A geostationary satellite orbits Earth in the equatorial plane with a period equal to Earth's rotational period (T = 24.0 h = 86 400 s). Determine the orbital radius, speed, angular momentum per unit mass, and total mechanical energy per unit mass. Use M_E = 5.972 × 10²⁴ kg and G = 6.674 × 10⁻¹¹ N·m²/kg².

Geostationary Orbit Calculation

1
Step 1 — Find the Orbital Radius from Kepler's Third Law

Starting from T² = (4π²/GM)r³, solve for r: r = (GMT²/(4π²))^1/3. Compute GM = (6.674 × 10⁻¹¹)(5.972 × 10²⁴) = 3.986 × 10¹⁴ m³/s². Then r³ = (3.986 × 10¹⁴)(86 400)² / (4π²) = (3.986 × 10¹⁴)(7.4649 × 10⁹) / 39.478 = 7.530 × 10²² m³. Taking the cube root yields the orbital radius.

r ≈ 4.224 × 10⁷ m ≈ 42 200 km (about 35 800 km above Earth's surface)
2
Step 2 — Compute the Orbital Speed

Using v = 2πr/T = 2π(4.224 × 10⁷ m) / (86 400 s):

v ≈ 3.07 × 10³ m/s ≈ 3.07 km/s
3
Step 3 — Angular Momentum per Unit Mass

For a circular orbit, L/m = vr = (3.07 × 10³)(4.224 × 10⁷):

L/m ≈ 1.30 × 10¹¹ m²/s
4
Step 4 — Total Mechanical Energy per Unit Mass

E/m = −GM/(2r) = −(3.986 × 10¹⁴) / (2 × 4.224 × 10⁷):

E/m ≈ −4.72 × 10⁶ J/kg
5
Step 5 — Interpret Results

The negative total energy confirms the satellite is gravitationally bound. The orbital speed of about 3.07 km/s is far below Earth's escape speed at that altitude (≈ 4.34 km/s), consistent with a bound circular orbit. The orbital radius of ≈ 6.6 Earth radii is the unique radius at which the period matches Earth's rotation, enabling the satellite to hover over a fixed point on the equator.

SECTION 7

Strengths & Limitations of the Keplerian Model

The Newtonian two-body framework we have developed is extraordinarily powerful for a wide range of orbital problems, but it rests on idealizations that may break down under certain conditions. Understanding these limitations is essential for appreciating when more sophisticated models are required, and when the simple framework is entirely sufficient.

Aspect	Strengths	Limitations
Two-body assumption	Yields exact analytic solutions (conic sections). Sufficient for most introductory satellite problems.	Breaks down for multi-body systems (e.g., Sun–Earth–Moon). Perturbation theory or numerical integration required.
Point-mass approximation	Valid when r ≫ radii of both bodies. Simplifies mathematics enormously.	Earth's oblateness (J₂ effect) causes orbital precession for LEO satellites; cannot be ignored in precision applications.
No drag	Good approximation for high-altitude orbits (GEO, MEO) where atmospheric density is negligible.	LEO satellites (< 600 km) experience measurable atmospheric drag, causing orbital decay and eventual re-entry.
Newtonian gravity	Accurate to many decimal places for Earth-orbiting satellites and most solar-system applications.	General relativistic corrections needed for Mercury's perihelion precession, GPS clock synchronization, and strong-field regimes.

Applicability of the Keplerian two-body model

✦

KEY TAKEAWAY

The Keplerian model is like a perfectly tuned guitar playing in a quiet room — it captures the dominant harmonics beautifully. In reality, perturbations from other bodies, atmospheric drag, and relativistic effects add small dissonances, much like background noise and string imperfections in a real performance. For most undergraduate-level problems and many engineering applications, the idealized model provides remarkably accurate predictions, but precision satellite operations demand the 'noise-cancellation' of perturbation theory.

SECTION 8

Connection to Advanced Theory

The orbit theory developed in this lesson serves as a launching pad for several more advanced treatments in both physics and engineering. The most immediate generalization involves orbital transfers — maneuvers that shift a satellite from one orbit to another by changing its energy and angular momentum through impulsive thrust. The Hohmann transfer, the simplest and most fuel-efficient two-impulse transfer between coplanar circular orbits, relies directly on the energy–radius and angular-momentum relations derived above. Beyond classical mechanics, general relativity introduces corrections that are small but measurable: GPS satellites, for instance, require both special and general relativistic time corrections amounting to about 38 microseconds per day, without which position errors would accumulate at roughly 10 km/day.

Feature	Newtonian Orbital Mechanics	Advanced / Relativistic Treatment
Gravity model	Inverse-square force, instantaneous action at a distance	Spacetime curvature (GR); gravitational waves propagate at c
Orbit precession	Closed ellipses (no precession in pure two-body)	Perihelion advance (Mercury: 43″ per century from GR)
Time & clocks	Absolute time, identical clock rates everywhere	Gravitational time dilation; clocks run faster at higher altitude
Multi-body problems	Perturbation theory or numerical N-body simulations	Lagrangian points, chaos in 3-body problem; restricted 3-body problem

Classical vs. advanced orbital mechanics

If you continue to courses in astrodynamics or general relativity, you will see how the energy and angular momentum framework generalizes seamlessly. The Schwarzschild metric introduces an effective potential that closely mirrors the Newtonian version but includes an additional −L²GM/(mc²r³) term, responsible for orbital precession. The conceptual machinery — effective potentials, conserved quantities, and orbit classification by energy — carries over almost verbatim. Mastering the Newtonian treatment therefore gives you a deep and lasting foundation for the advanced theory.

SECTION 9

Practice Problems

Question 1: CONCEPTUAL

A satellite is in a stable circular orbit around Earth. An astronaut inside the satellite releases a wrench from rest (relative to the satellite). Explain, using Newton's laws and the concept of free fall, why the wrench appears to float rather than falling to the floor of the satellite. Does gravity act on the wrench?

Show answer

Yes, gravity absolutely acts on the wrench — it experiences a gravitational acceleration of approximately g = GM/r² directed toward Earth's center. However, the satellite itself undergoes the same gravitational acceleration, so the satellite and the wrench share the same free-fall trajectory. In the non-inertial reference frame of the satellite, there is no relative acceleration between the wrench and the walls; hence the wrench appears to float. This is precisely the condition of weightlessness — not the absence of gravity, but the absence of a normal force. The distinction is important: gravity has not vanished; rather, every object in the satellite follows the same geodesic (or, in Newtonian language, the same parabolic/circular free-fall path).

Question 2: BASIC CALCULATION

The International Space Station (ISS) orbits at an altitude of approximately 408 km above Earth's surface. Given Earth's radius R_E = 6.371 × 10⁶ m and GM = 3.986 × 10¹⁴ m³/s², calculate the ISS's orbital speed and orbital period.

Show answer

The orbital radius is r = R_E + h = 6.371 × 10⁶ + 4.08 × 10⁵ = 6.779 × 10⁶ m. The orbital speed is v = √(GM/r) = √(3.986 × 10¹⁴ / 6.779 × 10⁶) = √(5.879 × 10⁷) ≈ 7.67 × 10³ m/s ≈ 7.67 km/s. The period is T = 2πr/v = 2π(6.779 × 10⁶)/(7.67 × 10³) ≈ 5 550 s ≈ 92.5 minutes. This matches the well-known ISS period, confirming the validity of the simple Keplerian model for LEO calculations.

Question 3: INTERMEDIATE

A satellite is in a circular orbit of radius r₁ = 7 000 km around Earth. Mission control commands a single tangential thruster burn to place the satellite into an elliptical transfer orbit with apoapsis at r₂ = 14 000 km. Using energy conservation, find the speed of the satellite immediately after the burn (i.e., at the periapsis of the transfer ellipse). Compare it to the original circular orbital speed.

Show answer

The original circular speed is v_circ = √(GM/r₁) = √(3.986 × 10¹⁴ / 7.0 × 10⁶) ≈ 7.546 km/s. The semi-major axis of the transfer ellipse is a = (r₁ + r₂)/2 = (7 000 + 14 000)/2 = 10 500 km = 1.05 × 10⁷ m. The total energy per unit mass of the transfer orbit is E/m = −GM/(2a) = −(3.986 × 10¹⁴)/(2 × 1.05 × 10⁷) = −1.898 × 10⁷ J/kg. At periapsis (r = r₁), we apply E/m = ½v² − GM/r₁, so ½v² = E/m + GM/r₁ = −1.898 × 10⁷ + 5.694 × 10⁷ = 3.796 × 10⁷ J/kg, giving v = √(2 × 3.796 × 10⁷) ≈ 8.714 km/s. This is about 1.17 km/s faster than the original circular speed — this Δv is the impulse the thruster must provide.

Question 4: APPLIED

A GPS satellite orbits Earth with a period of exactly 11 hours 58 minutes (half a sidereal day, T = 43 080 s). (a) Determine the orbital radius. (b) Calculate the orbital angular momentum per unit mass. (c) If the satellite's total mass is 2 000 kg, compute its total mechanical energy and compare it to the energy required to launch it from Earth's surface (approximate the surface launch energy as ½mv_esc² where v_esc = 11.2 km/s).

Show answer

(a) From r = (GMT²/(4π²))^(1/3): r³ = (3.986 × 10¹⁴)(43 080)² / (4π²) = (3.986 × 10¹⁴)(1.856 × 10⁹)/39.478 = 1.875 × 10²² m³, so r ≈ 2.656 × 10⁷ m ≈ 26 560 km. This corresponds to an altitude of about 20 190 km, consistent with the known GPS constellation altitude. (b) L/m = vr = (2πr/T)r = 2πr²/T = 2π(2.656 × 10⁷)²/43 080 ≈ 1.030 × 10¹¹ m²/s. (c) E = −GMm/(2r) = −(3.986 × 10¹⁴)(2000)/(2 × 2.656 × 10⁷) ≈ −1.50 × 10¹⁰ J. The launch energy is roughly ½mv²_esc = ½(2000)(11 200)² ≈ 1.25 × 10¹¹ J. The orbital energy magnitude (1.50 × 10¹⁰ J) is about 12% of the launch energy, highlighting that most energy goes into escaping the deep surface potential well, while 'parking' the satellite in its orbit costs comparatively less.

Question 5: CRITICAL THINKING

Consider a satellite in an elliptical orbit with periapsis distance r_p and apoapsis distance r_a around a planet of mass M. Using only conservation of energy and conservation of angular momentum (without invoking Kepler's laws directly), derive expressions for the speeds at periapsis (v_p) and apoapsis (v_a) in terms of G, M, r_p, and r_a. Then show that the ratio v_p/v_a = r_a/r_p.

Show answer

At both periapsis and apoapsis the velocity is purely tangential (perpendicular to the radius), so angular momentum conservation gives L = mv_pr_p = mv_ar_a, immediately yielding v_p/v_a = r_a/r_p. For the individual speeds, use energy conservation: ½v_p² − GM/r_p = ½v_a² − GM/r_a. Substituting v_a = v_pr_p/r_a and solving gives v_p = √[2GMr_a / (r_p(r_p + r_a))]. Symmetrically, v_a = √[2GMr_p / (r_a(r_p + r_a))]. These elegant results show that the speeds at the orbital extremes are completely determined by the two turning-point radii and the central mass, with no explicit reference to the semi-major axis or eccentricity needed.

SUMMARY

Lesson Summary

Summary

The motion of orbiting satellites is governed by the interplay of three fundamental ideas. Newton's law of gravitation provides the centripetal force that curves the satellite's path, uniquely determining the orbital speed v = √(GM/r) for a circular orbit. Conservation of mechanical energy (E = ½mv² − GMm/r = −GMm/(2a)) classifies orbits as bound (E < 0) or unbound (E ≥ 0) and connects the total energy to the semi-major axis. Conservation of angular momentum (L = mvr sinθ = const) explains why satellites speed up at periapsis and slow at apoapsis, and underpins Kepler's second law of equal areas.

Together, these principles yield Kepler's third law (T² ∝ r³), predict the geostationary orbit radius of ≈ 42 200 km, and enable the design of orbital transfer maneuvers. The effective potential diagram unifies orbit classification — circular, elliptical, parabolic, and hyperbolic — into a single visual framework. While the Newtonian two-body model has limitations (it neglects perturbations, drag, and relativistic effects), it remains the essential foundation for all of orbital mechanics and a showcase for the power of conservation laws in physics.