Aquesta tesi doctoral està emmarcada en el camp de l'astrodinàmica. Presenta solucions a problemes identificats en el disseny de missions que utilitzen òrbites entorn dels punts de libració, fent servir la teoria de sistemes dinàmics.El problema restringit de tres cossos és un model per estudiar el moviment d'un cos de massa infinitessimal sota l'atracció gravitatòria de dos cossos molt massius. Els cinc punts d'equilibri d'aquest model, en especial L1 i L2, han estat motiu de nombrosos estudis per aplicacions pràctiques en les últimes dècades (SOHO, Genesis...). Genèricament, qualsevol missió en òrbita al voltant del punt L2 del sistema Terra-Sol es veu afectat per ocultacions degudes a l'ombra de la Terra. Si l'òrbita és al voltant de L1, els eclipsis són deguts a la forta influència electromagnètica del Sol. D'entre els diferents tipus d'òrbites de libració, les òrbites de Lissajous resulten de la combinació de dues oscil.lacions perpendiculars. El seu principal avantatge és que les amplituds de les oscil.lacions poden ser escollides independentment i això les fa adapatables als requeriments de cada missió. La necessitat d'estratègies per evitar eclipsis en òrbites de Lissajous entorn dels punts L1 i L2 motivaren la primera part de la tesi. En aquesta part es presenta una eina per la planificació de maniobres en òrbites de Lissajous que no només serveix per solucionar el problema d'evitar els eclipsis, sinó també per trobar trajectòries de transferència entre òrbites d'amplituds diferents i planificar rendez-vous. Per altra banda, existeixen canals de baix cost que uneixen els punts L1 i L2 d'un sistema donat i representen una manera natural de transferir d'una regió de libració a l'altra. Gràcies al seu caràcter hiperbòlic, una òrbita de libració té uns objectes invariants associats: les varietats estable i inestable. Si tenim present que la varietat estable està formada per trajectòries que tendeixen cap a l'òrbita a la qual estan associades quan el temps avança, i que la varietat inestable fa el mateix però enrera en el temps, una intersecció entre una varietat estable i una d'inestable proporciona un camí asimptòtic entre les òrbites corresponents. Un mètode per trobar connexions d'aquest tipus entre òrbites planes entorn de L1 i L2 es presenta a la segona part de la tesi, i s'hi inclouen els resultats d'aplicar aquest mètode als casos dels problemes restringits Sol Terra i Terra-Lluna.La idea d'intersecar varietats hiperbòliques es pot aplicar també en la cerca de camins de baix cost entre les regions de libració del sistema Sol-Terra i Terra-Lluna. Si existissin camins naturals de les òrbites de libració solars cap a les lunars, s'obtindria una manera barata d'anar a la Lluna fent servir varietats invariants, cosa que no es pot fer de manera directa. I a l'inversa, un camí de les regions de libració lunars cap a les solars permetria, per exemple, que una estació fos col.locada en òrbita entorn del punt L2 lunar i servís com a base per donar servei a les missions que operen en òrbites de libració del sistema Sol-Terra. A la tercera part de la tesi es presenten mètodes per trobar trajectòries de baix cost que uneixen la regió L2 del sistema Terra-Lluna amb la regió L2 del sistema Sol-Terra, primer per òrbites planes i més endavant per òrbites de Lissajous, fent servir dos problemes de tres cossos acoblats. Un cop trobades les trajectòries en aquest model simplificat, convé refinar-les per fer-les més realistes. Una metodologia per obtenir trajectòries en efemèrides reals JPL a partir de les trobades entre òrbites de Lissajous en el model acoblat es presenta a la part final de la tesi. Aquestes trajectòries necessiten una maniobra en el punt d'acoblament, que és reduïda en el procés de refinat, arribant a obtenir trajectòries de cost zero quan això és possible. / This PhD. thesis lies within the field of astrodynamics. It provides solutions to problems which have been identified in mission design near libration points, by using dynamical systems theory. The restricted three body problem is a well known model to study the motion of an infinitesimal mass under the gravitational attraction of two massive bodies. Its five equilibrium points, specially L1 and L2, have been the object of several studies aimed at practical applications in the last decades (SOHO, Genesis...). In general, any mission in orbit around L2 of the Sun-Earth system is affected by occultations due to the shadow of the Earth. When the orbit is around L1, the eclipses are caused by the strong electromagnetic influence of the Sun. Among all different types of libration orbits, Lissajous type ones are the combination of two perpendicular oscillations. Its main advantage is that the amplitudes of the oscillations can be chosen independently and this fact makes Lissajous orbits more adaptable to the requirements of each particular mission than other kinds of libration motions. The need for eclipse avoidance strategies in Lissajous orbits around L1 and L2 motivated the first part of the thesis. It is in this part where a tool for planning maneuvers in Lissajous orbits is presented, which not only solves the eclipse avoidance problem, but can also be used for transferring between orbits having different amplitudes and for planning rendez-vous strategies.On the other hand, there exist low cost channels joining the L1 and L2 points of a given sistem, which represent a natural way of transferring from one libration region to the other one. Furthermore, there exist hyperbolic invariant objects, called stable and unstable manifolds, which are associated with libration orbits due to their hyperbolic character. If we bear in mind that the stable manifold of a libration orbit consists of trajectories which tend to the orbit as time goes by, and that the unstable manifold does so but backwards in time, any intersection between a stable and an unstable manifold will provide an asymptotic path between the corresponding libration orbits. A methodology for finding such asymptotic connecting paths between planar orbits around L1 and L2 is presented in the second part of the dissertation, including results for the particular cases of the Sun-Earth and Earth-Moon problems. Moreover, the idea of intersecting hyperbolic manifolds can be applied in the search for low cost paths joining the libration regions of different problems, such as the Sun-Earth and the Earth-Moon ones. If natural paths from the solar libration regions to the lunar ones was found, it would provide a cheap way of transferring to the Moon from the vicinity of the Earth, which is not possible in a direct way using invariant manifolds. And the other way round, paths from the lunar libration regions to the solar ones would allow for the placement of a station in orbit around the lunar L2, providing services to solar libration missions, for instance. In the third part of the thesis, a methodology for finding low cost trajectories joining the lunar L2 region and the solar L2 region is presented. This methodology was developed in a first step for planar orbits and in a further step for Lissajous type orbits, using in both cases two coupled restricted three body problems to model the Sun-Earth-Moon spacecraft four body problem. Once trajectories have been found in this simplified model, it is convenient to refine them to more realistic models. A methodology for obtaining JPL real ephemeris trajectories from the initial ones found in the coupled models is presented in the last part of the dissertation. These trajectories need a maneuver at the coupling point, which can be reduced in the refinement process until low cost connecting trajectories in real ephemeris are obtained (even zero cost, when possible).

This thesis evaluates decommissioning strategies for the International Space Station ISS. A permanent solution is attempted by employing energy efficient invariant manifolds that arise in the circular restricted three body problem CRTBP to transport the ISS from its low Earth orbit LEO to a lunar orbit. Although the invariant manifolds provide efficient transport, getting the the ISS onto the manifolds proves quite expensive, and the trajectories take too long to complete. Therefore a more practical, although temporary, solution consisting of an optimal re-boost maneuver with the European Space Agency's automated transfer vehicle ATV is proposed. The optimal re-boost trajectory is found using control parameterization and the sequential quadratic programming SQP algorithm. The model used for optimization takes into account the affects of atmospheric drag and gravity perturbations. The optimal re-boost maneuver produces a satellite lifetime of approximately ninety-five years using a two ATV strategy.

The general station-keeping problem is a focal topic when considering any spacecraft mission application. Recent missions are increasingly requiring complex trajectories to satisfy mission requirements, necessitating the need for accurate station-keeping controllers. An ideal controller reliably corrects for spacecraft state error, minimizes the required propellant, and is computationally efficient. To that end, this investigation assesses the effectiveness of several controller formulations in the circular restricted three-body model. Particularly, a spacecraft is positioned in a L₁ southern halo orbit within the Sun-Earth Moon Barycenter system. To prevent the spacecraft from departing the vicinity of this reference halo orbit, the Floquet mode station-keeping approach is introduced and evaluated. While this control strategy generally succeeds in the station-keeping objective, a breakdown in performance is observed proportional to increases in state error. Therefore, a new hybrid controller is developed which leverages Floquet mode and reinforcement learning. The hybrid controller is observed to efficiently determine corrective maneuvers that consistently recover the reference orbit for all evaluated scenarios. A comparative analysis of the performance metrics of both control strategies is conducted, highlighting differences in the rates of success and the expected propellant costs. The performance comparison demonstrates a relative improvement in the ability of the hybrid controller to meet the mission objectives, and suggests the applicability of reinforcement learning to the station-keeping problem.

With the growth of human interest in the Lunar region, methods of enabling Lunar access including surface and Low Lunar Orbit (LLO) from periodic orbit in the Lunar region is becoming more important. The current investigation explores the Lunar access capabilities of these periodic orbits. Impact trajectories originating from the 9:2 Lunar Synodic Resonant (LSR) Near Rectilinear Halo Orbit (NRHO) are determined through explicit propagation and mapping of initial conditions formed by applying small maneuvers at locations across the orbit. These trajectories yielding desirable Lunar impact final conditions are then used to converge impacting transfers from the NRHO to Shackleton crater near the Lunar south pole. The stability of periodic orbits in the Lunar region is analyzed through application of a stability index and time constant. The Lunar access capabilities of the Lunar region periodic orbits found to be sufficiently unstable are then analyzed through impact and periapse maps. Using the impact data, candidate periodic orbits are incorporated in the the NRHO to Shackleton crater mission design to control mission geometry. Finally, the periapse map data is used to determine periodic orbits with desirable apse conditions that are then used to design NRHO to LLO transfer trajectories.

<p dir="ltr">This research considers orbit maintenance strategies for multi-body orbits in the context of the Earth-Moon Circular Restricted Three Body Problem (CR3BP). Dynamical requirements for successful long-term orbit maintenance strategies are highlighted.</p>

<div> <div> <div> <p>Resonant orbits in a multi-body environment have been investigated in the past to aid the understanding of perceived chaotic behavior in the solar system. The invariant manifolds associated with resonant orbits have also been recently incorporated into the design of trajectories requiring reduced maneuver costs. Poincaré sections are now also extensively utilized in the search for novel, maneuver-free trajectories in various systems. This investigation employs dynamical systems techniques in the computation and characterization of resonant orbits in the higher-fidelity Circular Restricted Three-Body model. Differential corrections and numerical methods are widely leveraged in this analysis in the determination of orbits corresponding to different resonance ratios. The versatility of resonant orbits in the design of low cost trajectories to support exploration for several planet-moon systems is demonstrated. The efficacy of the resonant orbits is illustrated via transfer trajectory design in the Earth-Moon, Saturn-Titan, and the Mars-Deimos systems. Lastly, Poincaré sections associated with different resonance ratios are incorporated into the search for natural, maneuver-free trajectories in the Saturn-Titan system. To that end, homoclinic and heteroclinic trajectories are constructed. Additionally, chains of periodic orbits that mimic the geometries for two different resonant ratios are examined, i.e., periodic orbits that cycle between different resonances are determined. The tools and techniques demonstrated in this investigation are useful for the design of trajectories in several different systems within the CR3BP. </p> </div> </div> </div>

<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>An increasing interest in lunar exploration calls for low-cost techniques of reaching the Moon. Ballistic lunar transfers are long duration trajectories that leverage solar perturbations to reduce the multi-body energy of a spacecraft upon arrival into cislunar space. An investigation is conducted to explore methods of constructing ballistic lunar transfers. The techniques employ dynamical systems theory to leverage the underlying dynamical flow of the multi-body regime. Ballistic lunar transfers are governed by the gravitational influence of the Earth-Moon-Sun system; thus, multi-body gravity models are employed, i.e., the circular restricted three-body problem (CR3BP) and the bicircular restricted four-body problem (BCR4BP). The Sun-Earth CR3BP provides insight into the Sun’s effect on transfers near the Earth. The BCR4BP offers a coherent model for constructing end-to-end ballistic lunar transfers. Multiple techniques are employed to uncover ballistic transfers to conic and multi-body orbits in cislunar space. Initial conditions to deliver the spacecraft into various orbits emerge from Periapse Poincaré maps. From a chosen geometry, families of transfers from the Earth to conic orbits about the Moon are developed. Instantaneous equilibrium solutions in the BCR4BP provide an approximate for the theoretical minimum lunar orbit insertion costs, and are leveraged to create low-cost solutions. Trajectories to the <i>L</i>2 2:1 synodic resonant Lyapunov orbit, <i>L</i>2 2:1 synodic resonant Halo orbit, and the 3:1 synodic resonant Distant Retrograde Orbit (DRO) are investigated.</p>

<p>The 9:2 lunar synodic resonant near rectilinear halo orbit (NRHO) is the chosen orbit for the Gateway, a future lunar space station constructed by the National Aeronautics and Space Administration (NASA) as well as several commercial and international partners. Designing trajectories in this sensitive lunar region combined with the absence of a singular systematic methodology to approach mission design poses challenges as researchers attempt to design transfers to and from this nearly stable orbit. This investigation builds on previous research in Poincar\'e mapping strategies to design transfers from the 9:2 NRHO using higher-dimensional maps and maps with non-state variables. First, Poincar\'e maps are applied to planar transfers to demonstrate the utility of hyperplanes and establish that maps with only two or three dimensions are required in the planar problem. However, with the addition of two state variables, the spatial problem presents challenges in visualizing the full state. Higher-dimensional maps utilizing glyphs and color are employed for spatial transfer design involving the 9:2 NRHO. The visualization of all required dimensions on one plot accurately reveals low cost transfers into both a 3:2 planar resonant orbit and an L2 vertical orbit. Next, the application of higher-dimensional maps is extended beyond state variables. Visualizing time-of-flight on a map axis enables the selection of faster transfers. Additionally, glyphs and color depicting angular momentum rather than velocity lead to transfers with nearly tangential maneuvers. Theoretical minimum maneuvers occur at tangential intersections, so these transfers are low cost. Finally, a map displaying several initial and final orbit options, discerned through the inclusion of Jacobi constant on an axis, eliminates the need to recompute a map for each initial and final orbit pair. Thus, computation time is greatly reduced in addition to visualizing more of the design space in one plot. The higher-dimensional mapping strategies investigated are relevant for transfer design or other applications requiring the visualization of several dimensions simultaneously. Overall, this investigation outlines Poincar\'e mapping strategies for transfer scenarios of different design space dimensions and represents initial research into non-state variable mapping methods.</p>