Orbital mechanics - Trajectory Computation and Maneuver Design

Calculating Trajectories and Orbit Determination Introduction Determining a spacecraft's trajectory and calculating orbital transfers are fundamental problems in astrodynamics. Whether you're planning an interplanetary mission or maneuvering a satellite, you need reliable methods to predict where a spacecraft will be and when it will arrive. This section covers the essential tools: Kepler's equation for linking time to orbital position, universal variable methods that work across all orbit types, and the patched conic approximation that simplifies complex multi-body problems into manageable pieces. Kepler's Equation Kepler's equation is the bridge between time and position in an orbit. It relates three quantities: mean anomaly, eccentric anomaly, and eccentricity. What is anomaly? Anomaly is a measure of position along an orbit. The true anomaly $\theta$ is the angle from periapsis to the spacecraft's current location as measured from the central body. However, $\theta$ is inconvenient for calculating time. Instead, we use eccentric anomaly $E$, which is defined geometrically through an auxiliary circle and relates more directly to time. The mean anomaly $M$ is the position a spacecraft would have if it orbited at constant angular velocity. Kepler's equation connects these: $$M = E - e \sin E$$ where $e$ is the eccentricity. This equation is profound because $M$ increases linearly with time, so we can easily find how much time has elapsed. However, finding $E$ from $M$ requires solving this transcendental equation, which cannot be done algebraically. How to use it for time-of-flight calculations: You know the true anomaly $\theta$ where the spacecraft currently is or will be Convert $\theta$ to eccentric anomaly $E$ using: $\tan(E/2) = \sqrt{\frac{1-e}{1+e}} \tan(\theta/2)$ Calculate mean anomaly $M = E - e \sin E$ Multiply by orbital period to get elapsed time: $t = \frac{M}{2\pi} \cdot T{orbit}$ Solving Kepler's equation numerically: Since the equation is transcendental (involving both polynomial and trigonometric terms), we must solve it iteratively. Newton's method is the standard approach: $$E{new} = E{old} - \frac{M - E{old} + e \sin E{old}}{-1 + e \cos E{old}}$$ Starting with an initial guess (often $E0 = M$), you iterate until the solution converges. This typically takes only a few iterations. Universal Variable Formulation While Kepler's equation works beautifully for elliptical orbits, it breaks down for parabolic (escape) trajectories and becomes awkward for hyperbolic (unbound) orbits. The universal variable formulation elegantly solves this problem by using a single set of equations for all orbit types. The universal variable $\chi$ (chi) replaces the eccentric anomaly and is defined such that: For elliptical orbits: $\chi = \sqrt{a} E$ where $a$ is the semi-major axis For parabolic and hyperbolic orbits: it's defined through related geometric quantities The beauty of this approach is that the time-of-flight equation becomes: $$\Delta t = \sqrt{\frac{p^3}{\mu}} \left[ \chi - e \sin \chi + f(\chi) \right]$$ where $f(\chi)$ is a specific function and $\mu$ is the gravitational parameter. This single equation works whether your orbit is circular, elliptical, parabolic, or hyperbolic. This makes the universal variable method indispensable for trajectory design, where you might transition between different orbit types or work with escape trajectories. Patched Conic Approximation Calculating exact trajectories for a spacecraft moving under the gravity of multiple bodies (Earth, Sun, and Mars, for example) is computationally complex. The patched conic approximation provides an elegant simplification: assume that in each region of space, only one gravitating body dominates, and solve the two-body problem within each region. How it works: Consider an Earth-to-Mars transfer mission: Near Earth (sphere of influence of Earth): The spacecraft leaves Earth under Earth's gravity alone. Calculate the escape velocity or orbital transfer needed. Between planets (under the Sun's gravity): Once far enough from Earth that the Sun's gravity dominates, treat the spacecraft as orbiting the Sun and calculate the Hohmann transfer or other interplanetary trajectory. Near Mars (sphere of influence of Mars): When the spacecraft approaches Mars and Mars's gravity becomes dominant, calculate the approach trajectory and entry into Mars orbit. The boundaries between these regions are defined by the spheres of influence—regions where each body's gravitational effect is strongest. Why this matters: The patched conic approximation provides quick estimates of delta-v requirements and transfer times, which are essential for mission planning. However, because it ignores gravitational interactions between bodies in overlapping regions, the results are approximations. Precise navigation requires numerical integration of the full multi-body problem using high-fidelity computational tools. Orbital Maneuvers and Transfers What Is an Orbital Maneuver? An orbital maneuver is any change to a spacecraft's orbit achieved through propulsion. This includes raising or lowering orbits, changing inclination, and escaping to interplanetary space. Understanding efficient transfer methods is crucial because every kilogram of propellant costs thousands of dollars to launch. The measure of efficiency is delta-v ($\Delta v$), the total change in velocity required. A deep-space maneuver refers to a propulsive burn conducted far from Earth, such as a mid-course correction on the way to Mars or a trajectory adjustment near the Sun. These differ from near-Earth maneuvers by operating outside Earth's significant gravitational influence. Hohmann Transfer The Hohmann transfer is the workhorse of orbital mechanics. It transfers a spacecraft between two coplanar circular orbits using exactly two impulsive burns (instantaneous velocity changes) and requires the minimum delta-v of any two-impulse transfer method. How it works: Imagine a spacecraft in a circular orbit of radius $r1$ around Earth, and you want to reach a circular orbit of radius $r2$ where $r2 > r1$ (climbing to a higher orbit). First burn at periapsis: Fire the engines at the current orbit (which becomes periapsis of the transfer orbit). This raises the apoapsis to $r2$. Coast: Follow an elliptical transfer orbit halfway around. This takes time $\Delta t = \pi \sqrt{\frac{a{transfer}^3}{\mu}}$ where $a{transfer} = \frac{r1 + r2}{2}$. Second burn at apoapsis: Once you reach $r2$, fire again to circularize the final orbit. The total delta-v required is: $$\Delta v{total} = \Delta v1 + \Delta v2$$ where each burn can be calculated from the velocities needed: $$\Delta v1 = \sqrt{\frac{\mu}{r1}} \left( \sqrt{\frac{2r2}{r1 + r2}} - 1 \right)$$ $$\Delta v2 = \sqrt{\frac{\mu}{r2}} \left( 1 - \sqrt{\frac{2r1}{r1 + r2}} \right)$$ Why is this optimal? The Hohmann transfer is the minimum-energy two-impulse solution because the transfer ellipse is tangent to both circular orbits. Any faster transfer requires larger velocity changes, and any path that isn't tangent wastes energy. Bi-elliptic Transfer Sometimes, using more than two burns can actually require less total delta-v. The bi-elliptic transfer demonstrates this counterintuitive result. In a bi-elliptic transfer, you perform three burns: First burn: Raise apoapsis to some intermediate altitude $rb$ where $rb > r2$. Coast: Follow the first ellipse from $r1$ to $rb$. Second burn at $rb$: Raise the periapsis of a new ellipse to $r2$. Coast: Follow the second ellipse from $rb$ to $r2$. Third burn: Circularize at $r2$. The advantage emerges when the orbital radius ratio is large—roughly when $\frac{r2}{r1} > 12$. At these extreme altitude changes, the bi-elliptic transfer requires less total delta-v because the intermediate apoapsis allows more efficient velocity changes at both ends. The two ellipses give more flexibility in distributing the energy cost. However, bi-elliptic transfers take much longer (you make two half-revolutions around the intermediate orbit). For most practical missions, the time penalty outweighs the delta-v savings, but in special cases—particularly deep-space missions—bi-elliptic transfers are valuable. High-Thrust vs. Low-Thrust Transfers Spacecraft can be propelled by different engine types, and this dramatically affects transfer strategy. High-thrust transfers use conventional chemical rockets that produce a large acceleration over a short time. Hohmann and bi-elliptic transfers assume impulsive (instantaneous) burns, which are a good approximation for high-thrust systems. The advantages are fast transfers: a Hohmann transfer to Mars takes only a few months. The disadvantage is high delta-v requirements for large orbital changes. Low-thrust transfers use electric propulsion engines that produce continuous, gentle acceleration. Ion drives, for example, can operate continuously throughout the journey. Because the acceleration is small and spread over days or weeks, the entire trajectory curves continuously rather than following two sharp impulses. Low-thrust transfers are more complex to optimize but offer major advantages: Much lower propellant consumption (lower delta-v) Longer transfer times, but acceptable for cargo missions Can achieve orbits that high-thrust methods cannot The strategy often involves slowly spiraling outward, with continuous thrust applied near apogee. This exploits the Oberth effect: propulsive burns are most effective (produce the largest velocity change) when applied at high speed, such as at periapsis. A low-thrust engine continuously applying force near apogee will gradually raise the apoapsis, then as the orbit expands and velocity decreases everywhere, apply more burns to continue spiraling. Plane-Change Maneuvers So far we've discussed changing orbital altitude while staying in the same orbital plane. Plane-change maneuvers alter the orbital inclination or shift the orbital plane—essential for reaching different destinations like polar orbits. Finding the optimal location: The most efficient place to change inclination is at the ascending node (where the orbit crosses the equator moving northward) or descending node (moving southward)—the points where the two orbital planes intersect. At these nodes, you don't waste delta-v rotating your trajectory out of the original plane unnecessarily. Why apoapsis is preferable: Among the nodes, apoapsis is optimal because velocity is lowest there. Since delta-v cost for a plane change is proportional to the velocity, performing the burn when velocity is minimal minimizes the delta-v required. For a plane change of inclination $\Delta i$, the delta-v is approximately: $$\Delta v = 2 v \sin(\Delta i / 2)$$ where $v$ is the orbital speed. Lower $v$ means lower delta-v. The Oberth Effect twist: There's a subtle optimization here. A small inclination change can sometimes be incorporated into the periapsis raise burn due to the Oberth effect: propulsive burns produce larger velocity changes when applied at higher speeds (at periapsis). Combining a small plane change with a major periapsis burn at high speed is sometimes more efficient than performing a separate plane-change burn later at lower speed. This requires careful trajectory analysis to optimize. <extrainfo> For very large plane changes, multiple small burns distributed around the orbit may be more efficient than a single burn at a node, though this is rarely used in practice due to complexity. </extrainfo>