Fact-checked by Grok 2 weeks ago

Patched conic approximation

The patched conic approximation is a simplified analytical method in astrodynamics for modeling spacecraft trajectories through the gravitational influences of multiple celestial bodies, such as the Sun and planets, by dividing the path into discrete segments treated as two-body problems using conic sections (ellipses, hyperbolas, or parabolas). This approach approximates the complex n-body problem by assuming that within each segment, the gravity of one dominant body (e.g., a planet or the Sun) overwhelmingly controls the spacecraft's motion, while perturbations from other bodies are negligible. The method relies on the concept of a for each , defined as the radial distance from the body where its gravitational pull equals that of (or the ), marking the boundaries where segments are "patched" together. For an interplanetary transfer, the is typically subdivided into three main phases: a from the departure 's SOI (with the at the focus), a heliocentric elliptic phase (with at the focus) connecting the departure and arrival SOIs, and a capture or flyby into the target 's SOI (with the target at the focus). This patching ensures continuity of position and velocity at the SOI boundaries, allowing for efficient calculation of key parameters like launch energy (), arrival velocity (V∞), and transfer times. Widely applied in preliminary mission design for interplanetary probes, the approximation underpins efficient transfers such as the Hohmann ellipse, which minimizes propellant use by tangentally linking circular planetary orbits. It has been instrumental in planning historic missions, including those to Mars and outer planets, by enabling rapid trade studies on launch windows and trajectory options without full numerical integration. However, its accuracy diminishes for closely spaced bodies like Earth and the Moon, where the departure planet's SOI encompasses the target (e.g., Earth's SOI extends beyond lunar distance), necessitating more advanced models like the circular restricted three-body problem. Despite these limitations, the method remains a foundational tool due to its computational simplicity and sufficient precision for early-phase astrodynamics.

Fundamentals

Definition and Purpose

The in astrodynamics involves predicting the motion of a under the simultaneous gravitational influences of multiple celestial bodies, which is computationally intensive and lacks closed-form analytical solutions for more than two bodies. The patched conic approximation addresses this complexity by segmenting the trajectory into a series of two-body problems, where the spacecraft's motion relative to one dominant body is modeled using conic sections—such as ellipses for bound orbits, parabolas for marginal escapes, or hyperbolas for hyperbolic encounters—while neglecting perturbations from distant bodies until transitioning to the next segment. This hierarchical approach divides space into regions of gravitational dominance, often delineated by spheres of influence, allowing the overall path to be "patched" together at boundary points. The method originated in the late 1950s during the early planning of space missions, with foundational work by Y. V. Egorov in 1958 introducing it as an approximate solution to the restricted . It gained prominence in the through applications to lunar and interplanetary trajectories, as detailed in analyses like those by Lagerstrom and Kevorkian on matched-conic approximations for two fixed force centers. By 1971, the technique was comprehensively described in standard astrodynamics texts, reflecting its role in enabling the Apollo program's trajectory designs. The primary purpose of the patched conic approximation is to facilitate rapid preliminary design of trajectories by providing an efficient, analytical framework for estimating key parameters like velocity changes and transfer times. For instance, it models heliocentric motion dominated by , interrupted by brief planetary perturbations during close encounters, allowing mission planners to iterate designs without full numerical n-body integrations. This approximation builds on prerequisite concepts from classical , including Keplerian orbits, which describe ideal two-body motion under inverse-square gravitational attraction, and conic sections as the geometric solutions to these dynamics.

Spheres of Influence

The sphere of influence (SOI) is defined as the volume surrounding a secondary celestial body, such as a , within which the gravitational from that body on a third body (e.g., a ) dominates over the from the primary body (e.g., the Sun). This concept, introduced by in the late , delineates regions where two-body approximations are valid for trajectory modeling in multi-body systems. The radius of the SOI is given by Laplace's formula: r_{\text{SOI}} = a \left( \frac{m}{M} \right)^{2/5} where a is the semi-major axis of the secondary body's around the primary, m is the mass of the secondary body, and M is the mass of the primary body. This formula arises from in the restricted , where the SOI radius is derived by equating the ratios of the central to perturbing s for both the primary and secondary bodies: the (primary) central over its perturbing effect equals the planetary (secondary) central over its perturbing effect, under the approximation that the distance from the secondary is much smaller than its orbital distance from the primary. For orbiting , with a \approx 1.496 \times 10^8 km, m \approx 5.972 \times 10^{24} kg, and M \approx 1.989 \times 10^{30} kg, the SOI radius is approximately 925,000 km (about 145 radii). The choice of SOI size balances computational simplicity with reasonable accuracy in trajectory approximations, as the SOI radii are small relative to interplanetary distances—for instance, Earth's SOI is only about 0.006 compared to the Sun-Earth distance of 1 —allowing the neglect of distant perturbations outside the SOI without significant error in most mission planning scenarios. An alternative to Laplace's SOI is the Hill sphere, which approximates the region where a secondary body's can hold objects in stable orbits against the primary's forces, with radius r_H \approx a \left( \frac{m}{3M} \right)^{1/3}; however, Laplace's formulation is more prevalent in patched conic methods due to its alignment with perturbation-based switching criteria for interplanetary transfers.

Method

Trajectory Segmentation

In the patched conic approximation, the spacecraft's overall trajectory is divided into discrete segments at the boundaries of spheres of influence (SOIs), where the gravitational dominance shifts from one celestial body to another. Entry and exit points are identified precisely at these SOI radii, marking transitions in the dominant . Within a planet's SOI, the motion is modeled in geocentric (or planetocentric) coordinates as or elliptic orbits relative to that body, while outside the SOI—typically in interplanetary space—the trajectory switches to heliocentric coordinates dominated by the Sun's gravity. This division simplifies the multi-body problem by treating each segment as an isolated two-body interaction, assuming the influence of other bodies is negligible within the defined region. Each segment is approximated as a conic orbit solution to the , grounded in Kepler's laws, which describe the elliptical (or ) paths under a central inverse-square gravitational , conserving and . segments ( e > 1) represent departure from a , where the escapes the SOI with excess , or arrival at a target, involving a hyperbolic encounter flyby. Elliptic segments ($0 \leq e < 1) model cruise phases under solar gravity, such as Hohmann-like transfers between planetary SOIs, where the follows a bound orbit relative to the Sun. Orbit determination within segments relies on the vis-viva equation for energy-based calculations: v = \sqrt{ \mu \left( \frac{2}{r} - \frac{1}{a} \right) } Here, v is the speed at radial distance r, \mu is the gravitational parameter of the dominant body, and a is the semi-major axis, which distinguishes elliptic (a > 0) from hyperbolic (a < 0) regimes by relating total mechanical energy to orbital shape. To handle multiple gravitational bodies, the approximation employs a hierarchical structure, prioritizing the Sun as the primary influence for interplanetary legs, with planets or moons as secondary perturbers within their respective SOIs. Minor perturbations from distant bodies are ignored in each segment, reducing the n-body dynamics to a sequence of two-body solutions, such as Sun-spacecraft for heliocentric cruise or planet-spacecraft for periapsis passages. This nested approach extends naturally to systems like Earth-Moon or Jupiter-moons, where inner SOIs are treated relative to the outer body's frame.

Patching Procedure

The patching procedure in the patched conic approximation involves joining individual conic trajectory segments at the boundaries of spheres of influence (SOI) to form a continuous overall path, ensuring continuity in position and velocity while transitioning between dominant gravitational bodies. This process relies on coordinate transformations from planetocentric to heliocentric reference frames (or vice versa) and assumes an instantaneous switch at the SOI edge, neglecting higher-order perturbations during the transition. Position continuity is achieved by matching the spacecraft's location relative to the central body at the patch point, while velocity continuity is enforced through vector addition that accounts for the relative motion between the spacecraft and the planet. The procedure typically follows a structured sequence of steps to propagate and match segments. First, the trajectory segment within a planet's SOI—such as a or capture —is propagated analytically to the SOI boundary using two-body centered on that , yielding the 's position and planetocentric at the edge. Second, the heliocentric vector of the is computed by vectorially adding the 's heliocentric to the 's asymptotic relative to the (the hyperbolic excess , \vec{v}_\infty), ensuring the outgoing segment aligns with the heliocentric transfer . For arrival trajectories, the process is reversed: the heliocentric incoming is subtracted from the 's to obtain \vec{v}_\infty, which defines a approach tangent to the SOI, with the parameters solved to match the required periapsis conditions. In non-coplanar scenarios, where the transfer orbit plane differs from the 's , rotation matrices are applied to transform components between frames, aligning the incoming/outgoing asymptotes while preserving the magnitude of \vec{v}_\infty. The key matching equation for the asymptotic is \vec{v}_\infty = \vec{v}_\text{helio} - \vec{v}_\text{planet}, where \vec{v}_\text{helio} is the spacecraft's in the heliocentric frame from the transfer conic, and \vec{v}_\text{planet} is the 's heliocentric at the patch point; this vector difference determines the orientation and energy of the local . Across patches, energy is conserved within each conic segment via the , and is maintained relative to the dominant body, though the instantaneous transition implies no explicit conservation enforcement between segments beyond the match. Impulsive adjustments, if required for maneuvers like plane changes, are incorporated at the patch point but treated separately from the geometric matching. Sources of error in the patching procedure arise primarily from the assumption of an instantaneous transition at the SOI boundary, which overlooks gradual gravitational perturbations near the edge and can introduce small discontinuities in higher-fidelity n-body simulations. These approximations are typically accurate to within 1-2% for interplanetary transfers but degrade for close planetary encounters or highly perturbed environments.

Applications

Interplanetary Transfers

The patched conic approximation plays a central role in designing interplanetary transfers by simplifying the computation of heliocentric trajectories between planetary departure and arrival points. At its core, this involves solving , which determines the connecting two position vectors in the solar gravitational field over a specified time-of-flight, assuming two-body motion dominated by . This heliocentric leg is then patched with departure and arrival hyperbolas relative to the planets' spheres of influence, enabling efficient delta-v budgeting for the overall mission. The approach assumes instantaneous velocity changes at sphere-of-influence boundaries, providing a first-order estimate suitable for preliminary design. A baseline for such transfers is the Hohmann transfer, which represents the minimum-energy elliptic arc between two coplanar circular orbits, such as those of and Mars. The initial delta-v impulse, \Delta v_1 = \sqrt{\frac{\mu}{r_1}} \left( \sqrt{\frac{2 r_2}{r_1 + r_2}} - 1 \right), raises the perigee velocity to enter the transfer ellipse, where \mu is the solar gravitational parameter, r_1 the departure radius, and r_2 the arrival radius; the arrival delta-v is symmetric but subtractive for capture. To identify optimal launch windows, porkchop plots are generated, contouring C_3 (twice the squared hyperbolic excess velocity) against departure and arrival dates, revealing low-energy opportunities like Type I (short arc) or Type II (long arc) transfers recurring every synodic period. For Earth-Mars missions, these plots typically highlight windows with C_3 minima around 8-12 km²/s² every 26 months. For more distant targets, such as outer planets, bi-elliptic transfers extend the Hohmann concept by introducing an intermediate apogee beyond the target orbit, potentially reducing total delta-v by up to 10-15% when the radius ratio exceeds about 12, though they require longer flight times. This approximation is integrated into mission planning software like NASA's General Mission Analysis Tool (GMAT), where patched conics provide initial guesses for before refining with higher-fidelity models. In modern applications, the method supports cislunar transfers in NASA's , approximating delta-v for and lunar orbit insertion, and preliminary analyses for Mars sample return concepts, where it estimates round-trip budgets exceeding 10 km/s while accounting for ascent from the Martian surface. Patching also briefly incorporates planetary escape and insertion maneuvers to match heliocentric velocities.

Mission Example

A representative application of the patched conic approximation is a Hohmann-like transfer from () to Mars orbit during the 2026 launch window, which aligns with favorable planetary alignments for minimum-energy trajectories. This scenario assumes a spacecraft departing from a 300 km circular parking orbit around , targeting insertion into a low Mars orbit after the heliocentric cruise phase. The method divides the trajectory into three segments: an departure hyperbola, a heliocentric elliptical transfer , and a Mars arrival hyperbola, with velocity matching at the boundaries of the spheres of influence (SOIs). The trajectory begins with the departure , where the receives a trans-Mars injection (TMI) to achieve from 's . This has a C_3 \approx 12 km²/s², corresponding to a excess (v_\infty) of approximately 3.46 km/s relative to at the SOI boundary, located at about 925,000 km from 's center. The total delta-v for this segment, including the initial from , is roughly 3.2 km/s, leveraging the for efficiency. Upon exiting Earth's SOI, the velocity is patched to the heliocentric ellipse, an elliptical with a semi-major axis of approximately 1.26 , tangent to at perihelion (1 ) and Mars' orbit at aphelion (1.524 ). The transfer has an eccentricity of about 0.208 and requires no further burns during the coast phase, with a (TOF) of 259 days. The launch occurs when Mars is positioned approximately 44° ahead of in heliocentric longitude to ensure tangential arrival. The overall interplanetary delta-v contribution is about 0.4 km/s, bringing the cumulative total to ~3.6 km/s. At Mars' SOI entry, roughly 578,000 from Mars' center, the incoming heliocentric is patched to the arrival , with a v_\infty of approximately 2.7 /s relative to Mars. This requires a subsequent Mars insertion (MOI) burn or to capture into , but the patching assumes instantaneous adjustment at the SOI. To illustrate the patching, consider the vectors: at Earth's SOI exit, the geocentric v_\infty (magnitude 3.46 /s, directed outward along the departure asymptote) is added vectorially to Earth's heliocentric (~29.8 /s tangential) to yield the initial heliocentric (~32.6 /s). Similarly, at Mars' SOI entry, the heliocentric (~21.5 /s) is subtracted from Mars' heliocentric (~24.1 /s tangential) to obtain the Mars-relative v_\infty (magnitude 2.7 /s, inclined for hyperbolic approach). This addition/subtraction at SOI boundaries simplifies the into solvable two-body segments without numerical propagation. Assuming a November 2026 launch, the spacecraft arrives at Mars in August 2027, demonstrating the method's utility for timely missions. This approach mirrors trajectories in ongoing Mars exploration efforts, such as follow-on missions to the Perseverance rover, including sample return concepts that build on similar patched conic designs for efficient resource use.

Limitations and Extensions

Accuracy Considerations

The patched conic approximation introduces primary errors stemming from its neglect of third-body gravitational perturbations within spheres of influence (SOIs), such as the Moon's influence on during Earth escape trajectories, where the method assumes purely two-body between the and while ignoring lunar perturbations that can alter the hyperbolic exit by several percent in close approaches. Additionally, boundary mismatches arise due to the finite size of SOIs, as the abrupt switch from heliocentric to planetocentric at the SOI edge overlooks the gradual transition of gravitational dominance, leading to discontinuities in and that propagate along the trajectory. Quantitative assessments reveal that the approximation typically yields errors below 1% in () and departure excess velocity for interplanetary transfers to Mars, but arrival excess velocity (V∞) discrepancies can reach 8% when compared to full integrations, with errors generally below 5% in delta-v budgets for outer solar system missions like those to due to larger SOI scales relative to heliocentric distances. However, errors escalate to 5-10% or more near orbital resonances or low-energy swing-bys, where neglected perturbations cause semi-major axis deviations exceeding 10% of the Sun-planet distance in about 5% of cases, and studies demonstrate sensitivity to SOI choice, with variations in the Laplace-defined altering predicted miss distances by tens of thousands of kilometers. The method falters in problematic scenarios involving multi-body dominance, such as trajectories near Lagrange points or Trojan asteroids, where stable equilibria like L4 and L5 points require restricted dynamics that the two-body patching cannot capture, resulting in unphysical predictions or overshoots of hundreds of thousands of kilometers. For instance, missions like the Space Telescope's at Earth-Sun demand circular restricted modeling for accuracy, as patched conics fail to account for the balanced perturbations essential to such Lissajous-family orbits. These inaccuracies are commonly mitigated through mid-course corrections using onboard propulsion, which adjust the trajectory after launch to compensate for unmodeled perturbations, as demonstrated in simulated Apollo-era transfers where such maneuvers reduced final errors to within targeting tolerances despite initial patched conic basing.

Alternatives

One prominent alternative to the patched conic approximation is full of the n-body , exemplified by Cowell's method, which directly solves the equations governing dynamics under gravitational influences from multiple bodies using high-order integrators such as Runge-Kutta schemes. This approach delivers high-fidelity trajectories by accounting for all perturbations without simplifying assumptions, but it is computationally intensive, often requiring significant processing time for long-duration simulations. In comparison, the patched conic method's analytical solutions enable rapid preliminary assessments, making numerical integration suitable for scenarios demanding precision over speed. For missions involving close gravitational encounters, such as planetary flybys or point transfers, the restricted (RTBP)—particularly its circular variant (CR3BP)—provides a targeted alternative by modeling the motion of a massless under the influence of two massive primaries orbiting each other circularly. This framework reveals dynamical structures like invariant manifolds and equilibrium points, with the Jacobi constant serving as a that delineates accessible energy surfaces and zero-velocity curves, restricting motion to feasible regions in . Unlike the two-body focus of patched conics, RTBP captures resonant effects and chaotic behaviors near the primaries, enhancing accuracy for resonant or low-energy transfers. Hybrid methods bridge these approaches by leveraging patched conics for initial segmentation and heliocentric/planetocentric phasing, followed by numerical to incorporate perturbations and refine the , as employed in mission design tools for interplanetary navigation. For long-duration missions with continuous accelerations, such as operations, extensions incorporate averaged perturbed models that secularly average non-gravitational forces over orbital periods, reducing computational load while maintaining fidelity for multi-year heliocentric arcs. These , often iterated with optimization algorithms, balance the speed of analytical approximations with the realism of numerical methods. Selection between methods depends on mission phase and requirements: patched conics excel in early feasibility studies and sweeps for their low computational cost and ease of implementation, whereas or RTBP-based models are chosen for operational validation in high-precision contexts, such as deep space probes or GPS constellation maintenance, where errors from simplifications could impact success. Hybrid strategies are increasingly standard for complex missions, starting with patched conics to generate candidate trajectories before numerical refinement ensures compliance with constraints like fuel budgets or arrival windows.

References

  1. [1]
    [PDF] aas 07-160 comparison of a simple patched conic trajectory code to ...
    The patched conic approximation subdivides the planetary mission into three distinct trajectories and patches them together to create a single trajectory path.
  2. [2]
    Patched Conic - an overview | ScienceDirect Topics
    During the patched-conic approximation, we “turn off” the sun's gravitational pull when the spacecraft is inside the SOI of the smaller body (earth or ...
  3. [3]
    None
    Below is a merged summary of the Patched-Conic Approximation from *Fundamentals of Astrodynamics* by Bate, Mueller, and White, consolidating all information from the provided segments into a single, comprehensive response. To maximize detail and clarity, I will use a table in CSV format to organize the key aspects (Definition, Purpose, Simplification, Prerequisite Concepts, Historical Context, and Page References) across the different segments, followed by a narrative summary and a list of useful URLs. This approach ensures all details are retained and easily accessible.
  4. [4]
    [PDF] Patched Conic Approximations
    The patched conic approximation for interplanetary transfers assumes that the sphere of influence of a planet has an infinite radius when observed from the ...Missing: astrodynamics | Show results with:astrodynamics
  5. [5]
    [PDF] X "'1 - NASA Technical Reports Server (NTRS)
    Since the development of the patched conic technique, several methods of approximating the solution to the restricted problem of three bodies have been ...
  6. [6]
    Matched-conic approximation to the two fixed force-center problem
    Matched-conic approximation to the two fixed force-center problem. Lagerstrom, P. A.; ;; Kevorkian, J. Abstract. Planar motion of a particle of negligible mass ...Missing: patched | Show results with:patched
  7. [7]
  8. [8]
    [PDF] Effect of Coordinate Switching on Simulation Accuracy of Translunar ...
    Aug 9, 2008 · Figure 3.1: Variation in the size of the Sphere of Influence. When Laplace first derived the radius of this sphere, he used a value of φ = 0◦.
  9. [9]
    None
    Below is a merged summary of the Patched Conic Approximation based on all provided segments. To retain all information in a dense and organized manner, I will use a combination of narrative text and a table in CSV format to capture detailed specifics (e.g., equations, examples, and numerical values). The narrative provides an overview, while the table consolidates key details for clarity and completeness.
  10. [10]
    [PDF] Dynamical Systems, the Three-Body Problem and Space Mission ...
    Apr 25, 2011 · the patched conic approximation (or patched two-body approximation), dis- cussed in Bate, Mueller, and White [1971]. The strategy of the ...
  11. [11]
    [PDF] Interplanetary Mission Design Handbook:
    MIDAS uses Lambert's Theorem to calculate Earth to Mars ballistic transfer trajectories in which two-body conic motion and a central force field are assumed.
  12. [12]
    [PDF] 19690023905.pdf
    The case studies included here are variations of the patched conic mode ... of the path (H) with the sphere of influence boundary at a value of (t) Sm d .
  13. [13]
    [PDF] 19680010566.pdf - NASA Technical Reports Server (NTRS)
    Hohmann transfer is always the absolute optimum among two-impulse transfers. However, if the ratio between orbit radii is sufficiently large, then there exists ...
  14. [14]
    [PDF] Using the General Mission - Analysis Tool (GMAT)
    Feb 7, 2017 · ◇ Start with patched-conic initial guess for each segment. ◇ GMAT targeting sequence used to find smooth solution from segmented initial ...
  15. [15]
    Orthogonal Approximation of Invariant Manifolds in the Circular ...
    Apr 6, 2023 · The NASA-led Artemis program constitutes a series ... The earliest method of trajectory design between two bodies was the patched conic approach, ...
  16. [16]
    [PDF] A Search for Viable Venus and Jupiter Sample Return Mission ...
    The patched-conic approximation is an industry- accepted method when making a first-ordered analysis for interplanetary trajectories. The method has the.
  17. [17]
    [PDF] Interplanetary Mission Design Handbook:
    Oct 17, 2020 · Two classes of missions are normally used to described Earth-Mars transfers. ... Earth Departure Mars Arrival Earth Arrival Earth Arrival. Excess ...
  18. [18]
    Missions - Atomic Rockets
    Mar 29, 2023 · One of the cheapest maneuvers is called a "Hohmann transfer." Using one to go from Terra to Mars takes about 5,700 meters per second of delta-V ...
  19. [19]
    Lunar gravity assists using patched-conics approximation, three and ...
    Jul 1, 2019 · The goal of the present paper is to study the errors given by the patched-conics approximation in a lunar gravity assist maneuver.
  20. [20]
    A comparison of the “patched-conics approach” and the restricted ...
    The main goal of the present paper is to compare the results obtained with the well-known dynamics given by the restricted problem and the “patched-conics” ...
  21. [21]
    [PDF] Extending the Patched-Conic Approximation to the Restricted Four ...
    Using the vis-viva equation once again, the Earth-relative speed ν0 that the vehicle must have in order to initiate the hyperbolic transfer orbit at t0 ...
  22. [22]
    [PDF] nasa cr-66819 final report simulated trajectories error analysis ...
    based on patched conic trajectories. problem. - + A. - % A. SI. The ... the third midcourse correction, the actual trajectory still misses the nominal ...
  23. [23]
    [PDF] Evaluation of a new method of integrating the orbital equations of ...
    Two well-known methods of setting up the equations to be integrated are Cowell's method in which the total accelerations acting on the spacecraft are integrated ...
  24. [24]
    [PDF] 1 NASA CR-2255
    All these techniques are claimed to be much faster than numerical integration and considerably more accurate than the well known patched-conic approximation.
  25. [25]
    [PDF] AAS 22-624 OPTIMAL CONTROL IN THE CIRCULAR RESTRICTED ...
    Orbit design in the Circular Restricted Three-Body Problem (CR3BP) is a promising avenue, due to the beneficial properties of halo orbits. Modern space missions ...
  26. [26]
    [PDF] trajectory design in the spatial circular restricted - Purdue Engineering
    This dissertation is about trajectory design in the spatial circular restricted three-body problem, using higher-dimensional Poincaré maps.
  27. [27]
    [PDF] Extending the hybrid methodology for orbit propagation by fitting ...
    Apr 18, 2019 · The hybrid methodology for orbit propagation is a technique that allows improving the accuracy of any propagator for predicting the future ...Missing: patched conics
  28. [28]
    A Single-Averaged Model for the Solar Radiation Pressure Applied ...
    May 18, 2023 · In this work, we investigate the process of space debris mitigation from the GEO region using a solar sail.
  29. [29]
    The hybrid patched conic technique as applied to the lunar trajectory ...
    Hybrid patched conic technique as iterative procedure for generating translunar and transearth trajectories emphasizing computing time saving.Missing: propagation | Show results with:propagation
  30. [30]
    [PDF] ac round oordinat !# "%$ st ! & s - Purdue Engineering
    In an initial patched conic analysis, a timing condition is employed to construct an arbitrary number of trajectory arcs, which are then connected at lunar ...
  31. [31]
    Direct transfer trajectory design options for interplanetary orbiter ...
    Apr 1, 2017 · The proposed iterative patched conic method identifies the four distinct transfer trajectory design options available for a direct orbiter mission.