<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">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>