🛰️ Orbital Mechanics & Kepler's Laws

Virtual Laboratory | Satellite Communication Engineering

🎯 Laboratory Objectives

Upon completion of this virtual laboratory, students will be able to:

  • Demonstrate and verify Kepler's three laws of planetary motion in satellite orbital contexts, understanding how they govern all artificial satellite trajectories around Earth.
  • Calculate orbital parameters including semi-major axis, eccentricity, orbital period, and velocity for various satellite missions using fundamental orbital mechanics equations.
  • Visualize and compare GEO, MEO, LEO, and HEO orbital characteristics to understand the engineering trade-offs involved in selecting orbits for different communication applications.
  • Determine satellite visibility windows and ground track patterns for different orbit types, enabling Earth station antenna pointing and communication scheduling.
  • Analyze the relationship between orbital altitude and coverage area using spherical geometry to design optimal satellite constellations for global or regional coverage.

📚 Theory

Fundamental principles governing satellite orbital motion and their application to communication satellite systems.

🌍 Introduction to Orbital Mechanics

Orbital mechanics (or astrodynamics) is the application of ballistics and celestial mechanics to the practical problems concerning the motion of rockets and other spacecraft. For satellite communication engineers, understanding orbital mechanics is essential because the position, velocity, and visibility of a communication satellite directly determine link quality, coverage, and system design parameters.

A satellite in orbit around Earth is in a state of continuous free fall toward Earth, but its tangential velocity is sufficient that it continuously "misses" the Earth, resulting in a closed elliptical path. The motion is governed primarily by Earth's gravitational field, with perturbations from the Moon, Sun, atmospheric drag, and Earth's oblateness.

1. Kepler's First Law (Law of Ellipses)

"The orbit of every satellite is an ellipse with the Earth at one of the two foci."

For artificial satellites, this means that unless the orbit is perfectly circular (a special case of an ellipse where both foci coincide), the satellite's distance from Earth's center varies between a minimum (perigee) and maximum (apogee).

Major Axis (2a) Minor Axis (2b) Earth F₁ F₂ Perigee (rₚ) Apogee (rₐ) Center (C) a = semi-major axis

Figure 1: Elliptical orbit geometry showing Earth at one focus, with perigee and apogee distances.

r = a(1 - e²) / (1 + e·cos θ)
where r = radial distance, a = semi-major axis, e = eccentricity, θ = true anomaly
rₚ = a(1 - e)    and    rₐ = a(1 + e)
Perigee and apogee distances from Earth's center

2. Kepler's Second Law (Law of Equal Areas)

"A line joining a satellite and the Earth sweeps out equal areas during equal intervals of time."

This law is a consequence of the conservation of angular momentum. As a satellite moves from perigee toward apogee, it slows down; as it moves from apogee toward perigee, it speeds up. The velocity is maximum at perigee and minimum at apogee.

Earth F₂ Area A₁ Δt = 1 hour Area A₂ Δt = 1 hour A₁ = A₂ vₘₐₓ vₘᵢₙ Perigee Apogee

Figure 2: Kepler's Second Law — Equal areas swept in equal times. Satellite moves faster near perigee, slower near apogee.

dA/dt = L / (2m) = constant
where L = angular momentum, m = satellite mass. Areal velocity is constant.

3. Kepler's Third Law (Law of Harmonies)

"The square of the orbital period of a satellite is directly proportional to the cube of the semi-major axis of its orbit."

This law allows engineers to determine the orbital period from the orbit size, or conversely, to design an orbit with a specific period. For geostationary satellites, this law precisely determines the required altitude for a 24-hour period.

T² = (4π² / μ) · a³
where T = orbital period (s), a = semi-major axis (m), μ = GMₑ = 3.986 × 10¹⁴ m³/s² (Earth's gravitational parameter)
T = 2π √(a³/μ) = 5063.5 × √(a³)   [when a is in Earth radii]
Simplified form for quick calculations. For circular orbits, a = r = Rₑ + h.

4. Orbital Velocity

The velocity of a satellite at any point in its orbit can be determined from the vis-viva equation:

v = √[μ(2/r - 1/a)]
Vis-viva equation. For circular orbits (r = a), this reduces to v = √(μ/r).

5. Orbital Classifications for Communication Satellites

Orbit TypeAltitude RangePeriodKey Applications
LEO (Low Earth Orbit)200 – 2,000 km~90 min – 2 hrsIridium, Starlink, ISS, remote sensing
MEO (Medium Earth Orbit)2,000 – 35,786 km2 – 12 hrsGPS, Galileo, GLONASS, O3b
GEO (Geostationary)35,786 km23h 56m 4sIntelsat, DTH broadcast, weather
HEO (Highly Elliptical)Perigee ~500 km, Apogee ~40,000 km~12 – 24 hrsMolniya, Tundra (high-latitude coverage)

6. Visibility and Coverage Geometry

The coverage area (footprint) of a satellite depends on its orbital altitude and the minimum elevation angle required for communication. Using spherical geometry:

Earth Rₑ = 6371 km Satellite h = altitude r = Rₑ + h SSP Coverage θ ε Key Relations: SSP = Subsatellite Point cos(θ+ε) = Rₑ·cos(ε)/r Coverage angle: θ Slant range: d d² = Rₑ² + r² - 2Rₑr·cos(θ)

Figure 3: Satellite coverage geometry — coverage area is centered around the subsatellite point (SSP), the point on Earth's surface directly below the satellite.

θ = arccos[Rₑ/(Rₑ+h) · cos(ε)] - ε
Coverage half-angle θ from Earth's center, given altitude h and minimum elevation angle ε.
Coverage Area = 2πRₑ²(1 - cos θ)
Surface area on Earth visible to the satellite (spherical cap area).

7. Ground Tracks

A ground track is the path traced on Earth's surface directly below a satellite. For non-geostationary orbits, the ground track is determined by:

  • Orbital inclination (i): The angle between the orbital plane and Earth's equatorial plane. Determines the maximum latitude reached.
  • Orbital period: Determines how quickly the satellite completes one orbit.
  • Earth's rotation: Causes the ground track to shift westward by approximately 360° × (T / 86164) per orbit, where T is the orbital period in seconds and 86164s is the sidereal day.

For a circular orbit, the ground track forms a sinusoidal pattern when projected on a 2D map. The maximum latitude equals the orbital inclination for direct orbits, or 180° - i for retrograde orbits.

🧪 Laboratory Procedure

Follow these steps to complete the virtual laboratory experiments.

1

Study the Theory Section

Carefully review the theory section covering Kepler's three laws, orbital parameters, and coverage geometry. Ensure you understand the key equations: vis-viva equation, Kepler's Third Law, and coverage angle formula. Take note of the standard values: Earth's radius Rₑ = 6,371 km, gravitational parameter μ = 3.986 × 10¹⁴ m³/s².

2

Launch Simulation 1: Kepler's Laws Verification

Navigate to Simulation 1. Set the semi-major axis to 8,000 km and eccentricity to 0.3. Observe the elliptical orbit and verify that Earth is at one focus. Run the animation and confirm that the satellite moves faster near perigee and slower near apogee. Record the velocity at perigee and apogee. Verify that equal areas are swept in equal times by measuring the area swept in two different 1-hour intervals.

3

Launch Simulation 2: Orbital Parameter Calculator

Use the calculator to determine orbital parameters for the following scenarios: (a) LEO satellite at 800 km altitude with e = 0.001, (b) MEO satellite at 20,200 km altitude (GPS) with e = 0.02, (c) GEO satellite at 35,786 km with e = 0.0001. Record the orbital period, velocity at perigee/apogee, and specific mechanical energy for each case.

4

Launch Simulation 3: Orbit Type Comparison

Compare GEO, MEO, LEO, and HEO orbits side-by-side. Observe the differences in orbital period, ground track patterns, and Earth coverage. For the HEO orbit, set perigee = 500 km and apogee = 40,000 km. Note how the satellite "hovers" near apogee, making it useful for high-latitude communications. Record your observations on coverage characteristics.

5

Launch Simulation 4: Visibility & Ground Track Analyzer

Set an Earth station location (latitude, longitude) and satellite orbital parameters. Run the simulation to determine visibility windows — the time periods when the satellite is above the minimum elevation angle. Vary the orbital inclination and observe how ground track patterns change. For a LEO satellite at 45° inclination, calculate the maximum contact duration per pass and the number of passes per day.

6

Launch Simulation 5: Coverage Area Calculator

Input various orbital altitudes (from 200 km to 45,000 km) and minimum elevation angles (5°, 10°, 20°, 30°). Calculate and plot the coverage half-angle θ and coverage area for each combination. Determine the minimum number of satellites required for continuous global coverage at each altitude. Compare your results with known constellation designs (Iridium: 66 satellites at 780 km; GPS: 24 satellites at 20,200 km).

7

Data Collection & Analysis

Compile all simulation results into organized tables. Plot graphs of: (a) Orbital period vs. altitude, (b) Orbital velocity vs. altitude, (c) Coverage area vs. altitude for different elevation angles, (d) Number of satellites for global coverage vs. altitude. Perform error analysis by comparing simulation results with theoretical calculations.

8

Prepare Laboratory Report

Follow the report writing guidelines provided in the next section. Ensure all objectives are addressed, all data is presented with proper units and significant figures, and conclusions are supported by evidence from the simulations.

🔬 Interactive Simulations

Perform the following virtual experiments using the interactive simulation tools below.

Simulation 1: Kepler's Laws Verification

Visualize an elliptical satellite orbit and verify Kepler's Laws. Observe velocity changes and equal area sweeping.

8,000 km
0.30
3x
Perigee Distance
--
km
Apogee Distance
--
km
Velocity at Perigee
--
km/s
Velocity at Apogee
--
km/s
Orbital Period
--
hours
Current Velocity
--
km/s

Simulation 2: Orbital Parameter Calculator

Calculate complete orbital parameters for any satellite orbit. Enter altitude and eccentricity to compute all key values.

800 km
800 km

Simulation 3: Orbit Type Comparison

Compare GEO, MEO, LEO, and HEO orbits side-by-side. Visualize relative sizes, periods, and ground tracks.

ParameterLEOMEOGEOHEO
Altitude400 km20,200 km35,786 km500 × 40,000 km
Period1.5 hrs12 hrs23.9 hrs12 hrs
Velocity7.67 km/s3.87 km/s3.07 km/s1.5 – 10 km/s
Latency (RTT)~2.7 ms~135 ms~239 msVariable
Ground TrackMoving fastFigure-8StationaryLong dwell at apogee

Simulation 4: Visibility Window & Ground Track Analyzer

Determine when a satellite is visible from a ground station and visualize its ground track over Earth.

800 km
45°
10°

Simulation 5: Coverage Area Calculator

Calculate satellite coverage area based on altitude and minimum elevation angle. Determine constellation size for global coverage.

35,786 km

📝 Guidelines for Report Writing

Follow these guidelines to prepare a comprehensive and well-structured laboratory report.

1. Report Structure

  • Title Page: Include experiment title, student name, registration number, course code, date, and department.
  • Abstract/Executive Summary: A brief paragraph (150-200 words) summarizing the objectives, methods, key results, and conclusions.
  • Introduction: State the purpose of the experiment and its relevance to satellite communication engineering.
  • Theory: Present the fundamental equations and principles used, with proper citations. Include diagrams where appropriate.
  • Procedure: Describe the step-by-step methodology followed during the virtual experiments.
  • Results & Discussion: Present all simulation data in tables and graphs. Discuss the physical significance of your findings.
  • Conclusion: Summarize key findings and state whether objectives were achieved.
  • References: List all textbooks, papers, and online resources consulted (IEEE format preferred).

2. Data Presentation Requirements

  • All numerical values must include proper units and appropriate significant figures (typically 3-4 significant figures).
  • Tabulate comparison data for different orbit types (LEO, MEO, GEO, HEO) using a consistent format.
  • Include graphs with properly labeled axes, titles, and legends. Use different colors or line styles for distinct data series.
  • For Kepler's Third Law verification, plot T² vs. a³ and confirm linearity through the origin.
  • For coverage analysis, plot coverage area vs. altitude for multiple elevation angles on the same graph.
  • Include screenshots of simulation outputs where they support your analysis.

3. Analysis & Discussion Expectations

  • Kepler's Laws: Discuss how your simulation results verify each law. Explain any discrepancies between theoretical and simulated values.
  • Orbital Parameters: Compare calculated values with published data for real satellites (e.g., ISS, GPS, Intelsat). Calculate percentage errors.
  • Orbit Comparison: Analyze the engineering trade-offs between orbit types for communication applications. Consider latency, coverage, launch cost, and atmospheric drag.
  • Visibility Analysis: Discuss how orbital inclination affects visibility from different latitudes. Explain why polar orbits provide global coverage.
  • Coverage Geometry: Derive the minimum number of satellites theoretically required for continuous global coverage and compare with actual constellation designs.

4. Sample Questions for Discussion

  • Why does a geostationary orbit have zero inclination and zero eccentricity? What happens if these conditions are not met?
  • Calculate the velocity reduction required to move a satellite from a circular LEO to a transfer ellipse (Hohmann transfer) and then to GEO.
  • Explain why atmospheric drag is a significant concern for LEO satellites but negligible for GEO satellites.
  • A communication satellite at 800 km altitude must maintain a minimum elevation angle of 10°. What is the maximum slant range, and how does this affect link budget design?
  • Design a constellation of LEO satellites at 1,200 km altitude with 60° inclination that provides continuous global coverage. How many satellites are needed, and what is the orbital spacing?

5. Marking Rubric (Suggested)

ComponentWeightCriteria
Theory & Background15%Completeness, accuracy, proper referencing
Procedure Description10%Clarity, completeness, reproducibility
Data Presentation20%Tables, graphs, units, significant figures
Results & Analysis30%Correct calculations, physical insight, error analysis
Discussion & Conclusions15%Critical thinking, connections to real systems
Report Quality10%Organization, grammar, formatting, professionalism

6. Recommended References

  • Roddy, D. (2006). Satellite Communications (4th ed.). McGraw-Hill.
  • Maral, G., & Bousquet, M. (2009). Satellite Communications Systems (5th ed.). Wiley.
  • Pratt, T., Bostian, C., & Allnutt, J. (2003). Satellite Communications (2nd ed.). Wiley.
  • Wertz, J. R., & Larson, W. J. (1999). Space Mission Analysis and Design (3rd ed.). Microcosm Press.
  • ITU-R Recommendations: P.618, P.676, P.531 (for propagation models).
  • NASA Orbital Mechanics resources: https://orbitaldebris.jsc.nasa.gov/