<div>After completion of a resupply mission to NASA’s proposed Lunar Orbital Platform - Gateway, safe disposal of the Logistics Module is required. One potential option is disposal to heliocentric space. This investigation includes an exploration of the trajectory escape dynamics from an Earth-Moon L2 Near Rectilinear Halo Orbit (NRHO). The effects of the solar gravitational perturbations are assessed in the Bicircular Restricted 4-Body Problem (BCR4BP), as defined in the Earth-Moon rotating frame and in the Sun-B1 rotating frame, where B1 is the Earth-Moon barycenter. Disposal trajectories candidates are classified in three outcomes: direct escape, in direct escapes and captures.</div><div>Characteristics of each outcome is defined in terms of three parameters: the location of the apoapses within to the Sun-B1 rotating frame, a characteristic Hamiltonian value, and the osculating eccentricity with respect to the Earth-Moon barycenter. Sample trajectories are presented for each outcome. Low-cost disposal options are introduced.</div>

The purpose of this thesis is to design a machine learning algorithm platform that provides expanded knowledge of mission availability through a launch season by improving trajectory resolution and introducing launch mission forecasting. The specific scenario addressed in this paper is one in which data is provided for four deterministic translational maneuvers through a mission to a Near Rectilinear Halo Orbit (NRHO) with a 9:2 synodic frequency. Current launch availability knowledge under NASA’s Orion Orbit Performance Team is established by altering optimization variables associated to given reference launch epochs. This current method can be an abstract task and relies on an orbit analyst to structure a mission based off an established mission design methodology associated to the performance of Orion and NASA's Space Launch System. Introducing a machine learning algorithm trained to construct mission scenarios within the feasible range of known trajectories reduces the required interaction of the orbit analyst by removing the needed step of optimizing the orbit to fit an expected translational response required of the spacecraft. In this study, k-Nearest Neighbor and Bayesian Linear Regression successfully predicted classical orbital elements for the launch windows observed. However both algorithms had limitations due to their approaches to model fitting. Training machine learning algorithms off of classical orbital elements introduced a repetitive approach to reconstructing mission segments for different arrival opportunities through the launch window and can prove to be a viable method of launch window scan generation for future missions.

<p>A dynamical understanding of orbits in the Earth-Moon neighborhood that can sustain long-term activities and trajectories that link locations of interest forms a critical foundation for the creation of infrastructure to support a lasting presence in this region of space. In response, this investigation aims to identify and exploit fundamental dynamical motion in the vicinity of a candidate ‘hub’ orbit, the L2 southern 9:2 lunar synodic resonant near rectilinear halo orbit (NRHO), while incorporating realistic mission constraints. The strategies developed in this investigation are, however, not restricted to this particular orbit but are, in fact, applicable to a wide variety of stable and nearly-stable cislunar orbits. Since stable and nearly-stable orbits that may lack useful manifold structures are of interest for long-term activities in cislunar space due to low orbit maintenance costs, strategies to alternatively initiate transfer design into and out of these orbits are necessary. Additionally, it is crucial to understand the complex behaviors in the neighborhood of any candidate hub orbit. In this investigation, a bifurcation analysis is used to identify periodic orbit families in close proximity to the hub orbit that may possess members with favorable stability properties, i.e., unstable orbits. Stability properties are quantified using a metric defined as the stability index. Broucke stability diagrams, a tool in which the eigenvalues of the monodromy matrix are recast into two simple parameters, are used to identify bifurcations along orbit families. Continuation algorithms, in combination with a differential corrections scheme, are used to compute new families of periodic orbits originating at bifurcations. These families possess unstable members with associated invariant manifolds that are indeed useful for trajectory design. Members of the families nearby the L2 NRHOs are demonstrated to persist in a higher-fidelity ephemeris model. </p><p><br></p> <p>Transfers based on the identified nearby dynamical structures and their associated manifolds are designed. To formulate initial guesses for transfer trajectories, a Poincaré mapping technique is used. Various sample trajectory designs are produced in this investigation to demonstrate the wide applicability of the design methodology. Initially, designs are based in the circular restricted three-body problem, however, geometries are demonstrated to persist in a higher-fidelity ephemeris model, as well. A strategy to avoid Earth and Moon eclipse conditions along many-revolution quasi-periodic ephemeris orbits and transfer trajectories is proposed in response to upcoming mission needs. Lunar synodic resonance, in combination with careful epoch selection, produces a simple eclipse-avoidance technique. Additionally, an integral-type eclipse avoidance path constraint is derived and incorporated into a differential corrections scheme as well. Finally, transfer trajectories in the circular restricted three-body problem and higher-fidelity ephemeris model are optimized and the geometry is shown to persist.</p>

An integral part of NASA's vision for the coming years is a sustained infrastructure in cislunar space. The current baseline trajectory for this facility is a Near Rectilinear Halo Orbit (NRHO), a periodic orbit in the Circular Restricted Three-Body Problem. One of the goals of the facility is to serve as a proving ground for human spaceflight operations in deep space. Thus, this investigation focuses on transfers between the baseline NRHO and a family of periodic orbits that originate from a period-doubling bifurcation along the halo family. This new family of orbits has been termed the ``butterfly" family. This investigation also provides an overview of the evolution for a large subset of the butterfly family. Transfers to multiple subsets of the family are found by leveraging different design strategies and techniques from dynamical systems theory. The different design strategies are discussed in detail, and the transfers to each of these regions are compared in terms of propellant costs and times of flight.

<p dir="ltr">In the pursuit of establishing a sustainable space economy within the cislunar region, it is vital to formulate transfer design strategies that uncover economically viable highways between different regions of the space domain. The inherent complexity of spacecraft dynamics in the cislunar space poses challenges in determining feasible transfer options. However, the motion characterized by known dynamical structures modeled through the circular restricted three-body problem (CR3BP) aids in the identification of pathways with reasonable maneuver costs and flight times. A framework is proposed that incorporates a quasi-periodic orbit (QPOs) as an option to design transfer scenarios. This investigation focuses on the construction of transfers between periodic orbits. The framework is exemplified by the construction of pathways between an L2 9:2 synodic resonant Near-Rectilinear Halo Orbit (NRHO) and a planar Moon-centered Distant Retrograde Orbit (DRO). The innate difference in the geometries of the departure and arrival orbits of the sample case, along with the lack of natural flows towards and away from them, imply that links between these orbits may necessitate costly maneuvers. A strategy is formulated that leverages the stable and unstable manifolds associated with intermediate periodic orbits and quasi-periodic orbits to construct end-toend trajectories. As part of this strategy, a systematic methodology is outlined to streamline the determination of transfer options provided by the 5-dimensional manifolds associated with a QPO family. This approach reveals multiple local basins of solutions, both interior and exterior-types, characterized by selected intermediate orbits. The construction of transfers informed by the manifolds associated with QPOs is more intricate than those based on periodic orbits. However, QPO-derived solutions allow for the recognition of alternative local basins of solutions and often offer more cost-effective transfer options when compared to trajectories designed using periodic orbits that underlie the QPOs.</p>

<p dir="ltr">The Global Navigation Satellite System (GNSS) is heavily relied on for the navigation of Earth satellites. For satellites in cislunar space and beyond, GNSS is not readily available. As a result, other sources such as NASA's Deep Space Network (DSN) must be relied on for navigation. However, DSN is overburdened and can only support a small number of satellites at a time. Furthermore, communication with external sources can become interrupted or deprived in these environments. Given NASA's current efforts towards cislunar space operations and the expected increase in cislunar satellite traffic, there will be a need for more autonomous navigation options in cislunar space and beyond.</p><p dir="ltr">In this thesis, a navigation system capable of accurate and computationally efficient orbit determination in these communication-deprived environments is proposed and investigated. The emphasis on computational efficiency is in support of cubesats which are constrained in size, cost, and mass; this makes navigation even more challenging when resources such as GNSS signals or ground station tracking become unavailable.</p><p dir="ltr">The proposed navigation system, which is called GRAVNAV in this thesis, involves a two-satellite formation orbiting a planet. The primary satellite hosts an Extended Kalman Filter (EKF) and is capable of measuring the relative position of the secondary satellite; accurate attitude estimates are also available to the primary satellite. The relative position measurements allow the EKF to estimate the absolute position and velocity of both satellites. In this thesis, the proposed navigation system is investigated in the two-body and three-body problems.</p><p dir="ltr">The two-body analysis illuminates the effect of the gravity model error on orbit determination performance. High-fidelity gravity models can be computationally expensive for cubesats; however, celestial bodies such as the Earth and Moon have non-uniform and highly-irregular gravity fields that require complex models to describe the motion of satellites orbiting in their gravity field. Initial results show that when a second-order zonal harmonic gravity model is used, the orbit determination accuracy is poor at low altitudes due to large gravity model errors while high-altitude orbits yield good accuracy due to small gravity model errors. To remedy the poor performance for low-altitude orbits, a Gravity Model Error Compensation (GMEC) technique is proposed and investigated. Along with a special tuning model developed specifically for GRAVNAV, this technique is demonstrated to work well for various geocentric and lunar orbits.</p><p><br></p><p dir="ltr">In addition to the gravity model error, other variables affecting the state estimation accuracy are also explored in the two-body analysis. These variables include the six Keplerian orbital elements, measurement accuracy, intersatellite range, and satellite formation shape. The GRAVNAV analysis shows that a smaller intersatellite range results in increased state estimation error. Despite the intersatellite range bounds, semimajor axis, measurement model, and measurement errors being identical for both orbits, the satellite formation shape also has a strong influence on orbit determination accuracy. Formations that place both satellites in different orbits significantly outperform those that place both satellites in the same orbit.</p><p dir="ltr">The three-body analysis primarily focuses on characterizing the unique behavior of GRAVNAV in Near Rectilinear Halo Orbits (NRHOs). Like the two-body analysis, the effect of the satellite formation shape is also characterized and shown to have a similar impact on the orbit determination performance. Unlike the two-body problem, however, different orbits possess different stability properties which are shown to significantly affect orbit determination performance. The more stable NRHOs yield better GRAVNAV performance and are also less sensitive to factors that negatively impact performance such as measurement error, process noise, and decreased intersatellite range.</p><p dir="ltr">Overall, the analyses in this thesis show that GRAVNAV yields accurate and computationally efficient orbit determination when GMEC is used. This, along with the independence of GRAVNAV from GNSS signals and ground-station tracking, shows that GRAVNAV has good potential for navigation in cislunar space and beyond.</p>

<div>As destinations of missions in both human and robotic spaceflight become more exotic, a foundational understanding the dynamical structures in the gravitational environments enable more informed mission trajectory designs. One particular type of structure, quasi-periodic orbits, are examined in this investigation. Specifically, efficient computation of quasi-periodic orbits and leveraging quasi-periodic orbits as trajectory design alternatives in the Earth-Moon and Sun-Earth systems. First, periodic orbits and their associated center manifold are discussed to provide the background for the existence of quasi-periodic motion on n-dimensional invariant tori, where n corresponds to the number of fundamental frequencies that define the motion. Single and multiple shooting differential corrections strategies are summarized to compute families 2-dimensional tori in the Circular Restricted Three-Body Problem (CR3BP) using a stroboscopic mapping technique, originally developed by Howell and Olikara. Three types of quasi-periodic orbit families are presented: constant energy, constant frequency ratio, and constant mapping time families. Stability of quasi-periodic orbits is summarized and characterized with a single stability index quantity. For unstable quasi-periodic orbits, hyperbolic manifolds are computed from the differential of a discretized invariant curve. The use of quasi-periodic orbits is also demonstrated for destination orbits and transfer trajectories. Quasi-DROs are examined in the CR3BP and the Sun-Earth-Moon ephemeris model to achieve constant line of sight with Earth and avoid lunar eclipsing by exploiting orbital resonance. Arcs from quasi-periodic orbits are leveraged to provide an initial guess for transfer trajectory design between a planar Lyapunov orbit and an unstable halo orbit in the Earth-Moon system. Additionally, quasi-periodic trajectory arcs are exploited for transfer trajectory initial guesses between nearly stable periodic orbits in the Earth-Moon system. Lastly, stable hyperbolic manifolds from a Sun-Earth L₁ quasi-vertical orbit are employed to design maneuver-free transfer from the LEO vicinity to a quasi-vertical orbit.</div>