Interplanetary Trajectory Astrodynamics | Patched Conic Approximation

Porkchop plots are diagrams showing the delta-v required for missions as a function of launch date and arrival date. They are named for their resemblance to pork chops. The horizontal axis shows launch dates, the vertical axis shows arrival dates, and color intensity shows required delta-v. Blue regions indicate favorable launch windows with low delta-v requirements. Launch windows for Mars missions occur approximately every 26 months (about 2 years) when Earth and Mars are favorably positioned. The windows are shaped like butterfly wings due to the orbital inclination difference between Earth and Mars. The time difference between these windows can be up to 200 days.

To determine the interval between consecutive launch windows, solve for the smallest t > 0 such that the angular relationship between Mars and Earth repeats. Set α_Mars(t) - α_Earth(t) = α₀ + N×2π, where α₀ is the initial angular separation (44°). The solution yields t = 1/(1/T_Earth - 1/T_Mars) = 779 days. This represents the synodic period—the time between successive opportunities for interplanetary transfers between the same planets.

Porkchop plots are contour plots generated through brute-force optimization that display Earth departure dates on the x-axis and Mars arrival dates on the y-axis, with contour lines representing constant travel time, delta-v, departure C3, and arrival V-infinity; these plots help identify optimal launch windows by revealing regions of minimum delta-v requirements, where the gap in the plot separates trajectories with true anomaly changes less than 180 degrees (short-way, shorter travel times) from those greater than 180 degrees (long-way, longer travel times).

This segment explains how to interpret the Trajectory Planner's output, known as a 'pork chop plot.' The colored areas represent valid launch dates and times, while black areas indicate invalid solutions. The horizontal axis represents time, and the vertical axis shows required maneuver parameters (prograde, outward, plane change). The presenter demonstrates clicking within colored areas to find valid solutions, noting that the ideal launch window is typically in the center of the colored region. Different dates within the valid window require different Delta-V amounts, with edge dates requiring significantly more energy than ideal dates.

Launch windows are specific time periods when rockets must launch to successfully reach their destinations, determined by orbital mechanics and mission requirements; for low Earth orbit missions, they depend on matching orbital planes with targets like the ISS, while interplanetary missions rely on pork chop plots and synodic periods (e.g., Mars every 780 days) to find optimal trajectories, and lunar missions require precise timing based on the Moon's position and landing site illumination needs.

The restricted three-body problem considers two massive bodies (like Earth and Moon) and a third body of negligible mass that does not affect the motion of the massive bodies. This problem was crucial for designing the Apollo lunar missions, where spacecraft trajectories were calculated using the gravitational fields of Earth and Moon. The restricted three-body problem has special points (Lagrange points) where a small body can remain stationary relative to the two massive bodies. Points L4 and L5 form equilateral triangles with the two massive bodies and are stable equilibrium points. Points L1, L2, and L3 are unstable equilibrium points.

Lagrange points are equilibrium positions that emerge as solutions to the restricted three-body problem. For a system with the Sun, Earth, and a small body (like a probe or asteroid), there are five Lagrange points (L1-L5). Points L1, L2, and L3 lie on the line connecting the two main bodies and are unstable equilibrium points—small displacements cause the body to drift toward either the Sun or Earth. Points L4 and L5 form equilateral triangles with the two main bodies and are stable equilibrium points where a body can orbit without drifting toward either main body.

Lagrange points (L1-L5) are five special locations in the restricted three-body problem where a small object experiences zero net gravitational force relative to the two larger bodies. L1, L2, and L3 lie along the line connecting the two primary bodies, while L4 and L5 form equilateral triangles with them. Objects at these points can orbit in synchrony with the two primaries without propulsion, though they are generally unstable equilibrium points.

The restricted three-body problem reveals five special equilibrium points where gravitational forces balance to allow stable or quasi-stable positioning. L1 lies inside Earth's orbit where reduced gravity allows slower orbital speed matching Earth. L2 extends 1.5 million km beyond Earth where combined solar and Earth gravity provides enhanced pull for matching orbital velocity. L3 sits on the far side of the Sun, always hidden from Earth. L4 and L5 form triangular points 60 degrees ahead and behind Earth, providing stable equilibrium where objects remain trapped indefinitely. These points enable missions like Herschel and James Webb to maintain position relative to Earth while observing space without atmospheric interference.

Lagrange points are five equilibrium positions in the circular restricted three-body problem where a small third body can remain stationary relative to two larger bodies orbiting each other. These points were mathematically derived by Euler (linear solutions L1, L2, L3) and Lagrange (equilateral triangle solutions L4, L5). The collinear points (L1, L2, L3) are unstable, while the triangular points (L4, L5) can be stable depending on the mass ratio of the two primary bodies. These points have practical applications in space exploration, such as positioning solar observatories at L1 and L2, and the James Webb Space Telescope at L2, as well as the discovery of Trojan asteroids at Jupiter's L4 and L5 points.

Low-thrust propulsion systems, characterized by accelerations of one million G's or less and typically using electric propulsion like Hall thrusters, offer significant advantages over chemical propulsion including reduced satellite mass (achieving four-fifths of the mass required for chemical transfers), higher specific impulses (up to 3,000-5,000 seconds versus hundreds for chemical systems), and continuous controllability for trajectory corrections; however, they require longer transfer times and present complex mathematical optimization challenges that specialized software tools like Astrogator address through collocation methods and homotopy-based approaches to find optimal trajectories for applications such as geostationary orbit transfers and interplanetary missions.

Low-thrust engines require more delta-V than single-burn trajectories, creating uncertainty in mission planning. The 28,000 m/s velocity capability must account for extra delta-V spent during extended burns. Launch timing must occur well before optimal windows (35 days for Mars) to allow Earth SOI exit and solar orbit insertion. The strategy involves cycling up before plotting to other planets, exiting Earth's SOI in the correct direction, and plotting in solar orbit where extended burns can be performed. Continuous prograde thrust naturally circularizes orbits, simplifying trajectory management.

Trajectory optimization is essential because current launch vehicles cannot directly deliver large spacecraft to outer planets. Three trajectory categories exist by complexity: impulsive transfers (simplest, parameter optimization), low-thrust electric propulsion (continuous optimization), and impulsive transfers with gravity assists (most difficult). The Hohmann transfer provides efficient transfers between circular orbits using two impulses, while Lambert's problem solves non-Hohmann transfers analytically. Optimization objectives typically minimize propellant use, enabling smaller launch vehicles and maximizing scientific payload or mission lifetime.

The Kustaanheimo-Stiefel (KS) transformation linearizes orbital dynamics by expressing spacecraft position as a quaternion scaled from the unit x-axis, converting two-body motion into a simple harmonic oscillator in four dimensions; this linearization enables efficient trajectory optimization for low-thrust transfers like GTO-to-GEO, where the augmented Lagrangian solver can converge to optimal solutions without requiring initial guesses, as demonstrated by successful 30-day, 60-day, and 100-day transfer simulations achieving geostationary orbit with zero initial guess warm-starts.

Ion thrusters produce extremely low thrust—equivalent to holding a single sheet of paper—compared to chemical rockets that can lift a 30-story building. Yet this apparent weakness is actually their strength. While chemical rockets burn fuel for only about 8 minutes before coasting, ion thrusters can operate continuously for months or years. In the vacuum of space with no atmospheric drag, this sustained low thrust accumulates over time, gradually building velocity. After one year of continuous operation, an ion-powered spacecraft reaches speeds chemical rockets can only dream of. This approach represents humanity's transition from fighting against environmental constraints to harmonizing with the fundamental physical laws of the cosmos.

This comprehensive section covers the fundamental concepts of orbital perturbations. Keplerian orbits assume two-body problems with spherical symmetry, but real-world effects cause deviations. Major perturbation sources include Earth's oblateness (J2 effect), lunar/solar gravity, radiation pressure, and atmospheric drag. The relative magnitudes are: Earth's oblateness > lunar/solar > solar radiation > atmospheric drag. Three methods address perturbations: updating Keplerian elements from observations, using expanded elements, or numerical integration. GPS uses daily updates; GLONASS uses numerical integration. The mathematical framework uses zonal harmonics J_k and Legendre polynomials P_k, with J2 dominating. The mean velocity correction factor is [1 + (3/2)(J₂r²/a²)(sin²I)]. For 100 km altitude at 55° inclination, this yields ~0.06°/day precession. Effects diminish with altitude as Earth appears more point-like.

The J2 perturbation is the dominant gravitational effect caused by Earth's oblateness (its equatorial bulge), which causes the gravitational acceleration vector to point toward the equatorial plane rather than directly toward Earth's center; this perturbation significantly affects orbital elements such as the argument of perigee and right ascension of the ascending node, causing them to change over time, and must be included in accurate orbit propagation models alongside the two-body gravitational force.

The J2 perturbation is the dominant gravitational effect caused by Earth's oblateness (its equatorial bulge, where the polar radius is about 40 km less than the equatorial radius), which creates a perturbing acceleration on spacecraft in low Earth orbit. This perturbation can be derived by expanding the geopotential function using Legendre polynomials up to degree 2 and order 0, yielding the expression for the perturbing acceleration in the orbital reference frame as fr = -3μJ2RE²/(2r⁴)[1 - 3sin²i sin²u], fθ = -3μJ2RE²/(2r⁴)sin²i sin(2u), and fz = -3μJ2RE²/(2r⁴)sin(2i)sinu, where RE is Earth's mean equatorial radius (6378 km), μ is Earth's gravitational parameter, i is the orbital inclination, and u is the argument of latitude. This perturbation causes secular variations in the orbital elements, particularly affecting the right ascension of the ascending node and the argument of perigee, while having no effect on the semi-major axis or eccentricity.

For a sphere of uniform mass, gravitational potential is V = -μ/r with acceleration a = -∇V = -μr̂/r². Earth is an oblate spheroid, flattened at poles and bulging at equator, requiring modification of the spherical potential. The gravitational potential uses spherical coordinates and is expressed as a series: V = -μ/r [1 + Σ J_k (R/r)^k P_k(cos θ)], where J_k are dimensionless constants from satellite observations, R is Earth's equatorial radius, and P_k are Legendre polynomials. J2 is the most important term, representing primary flattening. Computing the J2 gradient in Cartesian coordinates requires the chain rule: ∂P/∂x = (∂P/∂r)(∂r/∂x) + (∂P/∂θ)(∂θ/∂x), with ∂r/∂x = x/r and ∂θ/∂x = -xz/r³. The final J2 perturbation acceleration components are: a_x = -3J2μx/r² × (5z²/r² - 1)/r², a_y = -3J2μy/r² × (5z²/r² - 1)/r², a_z = -3J2μz/r² × (5z²/r² - 3)/r². The perturbation is zero at the equator and changes sign at approximately ±35.26° latitude.

This section explains Earth's J2 effect (equatorial bulge of 44 km) and its impact on orbital elements. J2 creates torque affecting Right Ascension of Ascending Node (RAAN) and argument of perigee. For Sun-synchronous orbits (95-105° inclination), J2 causes RAAN precession matching Earth's solar orbit, enabling consistent equatorial passes. For Molniya orbits (63.4°/116.6° inclination), J2 fixes perigee position, allowing satellites to remain over northern hemisphere for most of their 12-hour period.