The development of a self-contained, onboard, fully autonomous trajectory guidance tool for spacecraft is presented. To be considered completely autonomous requires the capability to both identify an appropriate startup solution, and then use that solution to target a set of user-defined path and endpoint constraints. To minimize the cost of flight software development and validation, both the generation of the startup solution and the targeting algorithm are designed to be as computationally efficient as possible. This study addresses both the determination of a startup arc and the subsequent targeting process. The first part of the investigation considers the targeting algorithm. Linear targeting through differential corrections is a well-known approach for identifying feasible solutions that meet specified mission and trajectory constraints. However, to date, these methods relied on the assumption that the associated control inputs were impulsive in nature. This research focuses on the theoretical development and numerical validation of a generalized linear targeting algorithm capable of accommodating finite periods of continuous control action for a wide range of applications. Examples are presented to illustrate the general concept and to contrast the performance of this new targeting process against more classical impulsive targeting methods. The second section of the study introduces a novel approach utilizing artificial potential function methods to identify suitable startup solutions. Although common in other types of path planning, these methods have not yet been used for orbital or interplanetary trajectory design, primarily due to their inherent suboptimality. However, results show that this issue can be addressed with relative ease by the targeting algorithm. / text

In recent times, as the National Aeronautics and Space Administration (NASA) focuses on establishing a sustained presence in cislunar space, there has been an increase in planned missions to the cislunar vicinity for lunar exploration. Due to this increase in planned missions, the use of cislunar structures available in the Circular Restricted Three-Body Problem (CRTBP) has become of greater interest. Traditionally, transfers that leverage CRTBP structures in the cislunar vicinity have been generated as point designs. As a consequence of the non-linearity of this model, transitioning these point designs to other epochs or mission scenarios is non-trivial. Hence, a trade space of transfer solutions, that leverage the underlying dynamics, is of interest for rapid mission design. In this study, numerical methods and dynamical systems theory are leveraged to extract available dynamical structures in the model, which are subsequently exploited for transfer design. A trade space of relatively low time of flight, two-maneuver transfers, from a 500 km altitude Low Earth Orbit (LEO) to the Earth-Moon L1 Lyapunov orbit family is generated and analyzed.

Since the discovery of Ceres, asteroids have been of immense scientific interest and intrigue. They hold answers to many of the fundamental questionsabout the formation and evolution of the Solar System. Therefore, a missionsurveying the asteroid belt with close encounter of carefully chosen asteroidswould be of immense scientific benefit. The trajectory of such an asteroidtour mission needs to be designed such that asteroids of a wide range ofcompositions and sizes are encountered; all with an extremely limited ∆Vbudget.This thesis presents a novel heuristic algorithm to optimize trajectoriesfor an asteroid tour mission with close range flybys (≤ 1000 km). The coresearch algorithm efficiently decouples combinatorial (i.e. choosing the asteroids to flyby)and continuous optimization (i.e. optimizing critical maneuversand events) of what is essentially a mixed integer programming problem.Additionally, different methods to generate a healthy initial population forthe combinatorial optimization are presented.The algorithm is used to generate a set of 1800 feasible trajectories withina 2029+ launch frame. A statistical analysis of these set of trajectories isperformed and important metrics for the search are set based on the statistics.Trajectories allowing flybys to prominent families of asteroids like Flora andNysa with ∆V as low as 4.99 km/s are obtained.Two modified implementations of the algorithm are presented. In a firstiteration, a large sample of trajectories is generated with a limited numberof encounters to the most scientifically interesting targets. While, a posteriori, trajectories are filled in with as many small targets as possible. Thisis achieved in two different ways, namely single step extension and multiplestep extension. The former fills in the trajectories with small targets in onestep, while the latter optimizes the trajectory by filling in with one asteroid per step. The thesis also presents detection of asteroids for successfullyperforming flybys. A photometric filter is developed which prunes out badlyilluminated asteroids. The best trajectory is found to perform well againstthis filter such that nine out of the ten planned flybys are feasible.

<div> <div> <div> <p>As cislunar and outer space exploration regains worldwide popularity, the low-thrust spacecraft technology, whether in the form of solar sails, electric propulsion or nuclear propulsion, has seen a major increase in the last two decades, as new technologies arise that not only seek for a reduction of the size of the spacecraft —and/or the payloads— but also to minimize the cost of spaceflights, while trying to approach further destinations in our solar system. Mission designers are being challenged with the need to develop new strategies to generate rapid and informed initial guesses for low-thrust spacecraft trajectory design, that are easily converged into fully continuous solutions in position, velocity and mass states, in a high-fidelity dynamical model that incorporates the true ephemerides and perturbations of the gravitational attracting bodies acting on the spacecraft as it navigates through space. </p> <p>In an effort to explore further mission options for spacecraft traveling in the lunar vicinity, new interest arises into the problem of constructing a general framework for the initial guess generation of low-thrust trajectories in cislunar space, that is independent of the force models in which the orbits of interest are de ned. Given the efficiency of the low-thrust engines, most vehicles are equipped to perform further exploration of the cislunar space after completion of their primary science and technology demonstrations in orbits around the Moon. In this investigation, a generalized strategy for constructing initial guesses for low-thrust spacecraft traveling between lunar orbits that exist within the context of multiple dynamical models is presented. These trajectories are converged as mass-optimal solutions in lower fidelity model, that are easily transitioned and validated in the higher-fidelity ephemeris model, and, achieve large orbital plane changes while evolving entirely within the cislunar region. </p> <p>The robustness of the initial guess generation of the spacecraft’s path, depends highly on the fidelity of the dynamical model utilized to construct such trajectories, as well as on the numerical techniques employed to converge and propagate them into continuous solutions. Other researchers have extensively investigated novel techniques for the generation of initial guesses for the low-thrust spacecraft trajectory design problem including, but not limited to, patched conics strategies, methodologies for the transformation of impulsive burns into nite burns, the orbit chaining framework and, more recently, artificial intelligence schemes. This investigation develops an adaptive orbit chaining type approach that relies on the energy parametrization of periodic orbits that exist within the context of the circular restricted three-body problem, to construct informed initial guess for the low-thrust spacecraft trajectory.</p> <div> <div> <div> <p>A variety of multiple transfer applications for vehicles traveling between orbits in the cislunar region is explored for a wide range of low-thrust spacecraft with varying thrust acceleration magnitude. The examples presented in this investigation are consistent with the low-thrust parameters of previously own missions that utilized the same propulsion capabilities, such as, the DAWN mission and the Japanese Hayabusa missions 1 and 2. The trajectories presented in this work are optimized for either propellant consumption or time- of-flight in the lower-fidelity model, and later transitioned into a higher-fidelity ephemeris model that includes the gravitational attraction of the Sun, the Earth and the Moon. </p> <p>Two strategies are explored for the transition of trajectories from a lower-fidelity model to the higher-fidelity ephemeris model, both of which are successful in retaining the transfer geometry. The framework presented in this investigation is further applied to the upcoming NASA Lunar IceCube (LIC) mission to explore possible extended mission options once its primary science and technology demonstration objectives are achieved. It is demonstrated in this investigation that the strategies developed and presented in this work are not only applicable to the specific low-thrust vehicles explored, but it is applicable to any spacecraft with any type of propulsion technology. Furthermore, the energy-informed adaptive algorithm is easily transition to generate trajectories in a range of varying dynamical models. </p> </div> </div> </div> </div> </div> </div>

<div>In recent decades the revolutionary possibilities of low-thrust electric propulsion have been demonstrated by the success of missions such as Dawn and Hayabusa 1 and 2. The efficiency of low-thrust engines reduces the propellant mass required to achieve mission objectives and this benefit is frequently worth the additional time of flight incurred, particularly for robotic spacecraft. However, low-thrust trajectory design poses a challenging optimal control problem. At each instant in time, spacecraft control parameters that minimize an objective, typically propellant consumption or time of flight, must be determined. The characteristics of low-thrust optimal solutions are often unintuitive, making it difficult to develop an <i>a priori</i> estimate for the state and control history of a spacecraft that can be used to initialize an optimization algorithm. This investigation seeks to develop a low-thrust trajectory design framework to address this challenge by combining the existing techniques of orbit chaining and direct collocation. Together, these two methods offer a novel approach for low-thrust trajectory design that is intuitive, flexible, and robust.</div><div><br></div><div>This investigation presents a framework for the construction of orbit chains and the convergence of these initial guesses to optimal low-thrust solutions via direct collocation. The general procedure is first demonstrated with simple trajectory design problems which show how dynamical structures, such as periodic orbits and invariant manifolds, are employed to assemble orbits chains. Following this, two practical mission design problems demonstrate the applicability of this framework to real world scenarios. An orbit chain and direct collocation approach is utilized to develop low-thrust transfers for the planned Gateway spacecraft between a variety of lunar and libration point orbits (LPOs). Additionally, the proposed framework is applied to create a systematic method for the construction of transfers for the Lunar IceCube spacecraft from deployment to insertion upon its destination orbit near the Moon. Three and four-body dynamical models are leveraged for preliminary trajectory design in the first and second mission design applications, respectively, before transfers are transitioned to an ephemeris model for validation. Together, these realistic sample applications, along with the early examples, demonstrate that orbit chaining and direct collocation constitute an intuitive, flexible, and robust framework for low-thrust trajectory design. </div>

Analytical solutions to complex trajectory design problems are scarce, since only a few specific cases allow for closed-form solutions. The main purpose of this dissertation is to design simple algorithms for trajectory design using continuous thrust, with a focus on low-thrust applications. By “simple” here we seek to achieve algorithms that either admit an analytical solution, or require minimal input by the user and minimal computation time. The three main contributions of this dissertation are: designing Lyapunov-based closed-loop guidance laws for orbit transfers, finding semi-analytical solutions using a constant magnitude thrust, and perturbation theory for approximate solutions to low-thrust problems. The technical aspect that these problems share in common is that they all use, fully or partially, a thrusting model in which the energy of the system is kept constant. Many orbit transfer problems are shown to be solved with this thrusting protocol. / text

Approved for public release, distribution is unlimited / This thesis explores using a direct pseudospectral method for the solution of optimal control problems with mixed dynamics. An easy to use MATLAB optimization package known as DIDO is used to obtain the solutions. The modeling of both low thrust interplanetary trajectories as well as aerocapture trajectories is detailed and the solutions for low thrust minimum time and minimum fuel trajectories are explored with particular emphasis on verification of the optimality of the obtained solution. Optimal aerocpature trajectories are solved for rotating atmospheres over a range of arrival Vinfinities. Solutions are obtained using various performance indexes including minimum fuel, minimum heat load, and minimum total aerocapture mass. Finally, the problem formulation and solutions for the mixed dynamic problem of low thrust trajectories with a terminal aerocapture maneuver is addressed yielding new trajectories maximizing the total scientific mass at arrival. / Lieutenant, United States Navy

An automatic algorithm for accurate numerical gradient calculations has been developed. The algorithm is based on both finite differences and Chebyshev interpolation approximations. The novelty of the method is an automated tuning of the step size perturbation required for both methods. This automation guaranties the best possible solution using these approaches without the requirement of user inputs. The algorithm treats the functions as a black box, which makes it extremely useful when general and complex problems are considered. This is the case of spacecraft trajectory design problems and complex optimization systems. An efficient procedure for the automatic implementation is presented. Several examples based on an Earth-Moon free return trajectory are presented to validate and demonstrate the accuracy of the method. A state transition matrix (STM) procedure is developed as a reference for the validation of the method. / text

While the interest in the Martian moons increases, the low-thrust propulsion technology is expected to enable novel mission scenarios but is associated with unique trajectory design challenges. Accordingly, the current investigation introduces a multi-phase low-thrust design framework. The trajectory of a potential spacecraft that departs from the Earth vicinity to reach both of the Martian moons, is divided into four phases. To describe the motion of the spacecraft under the influence of gravitational bodies, the two-body problem (2BP) and the Circular-Restricted Three Body Problem (CR3BP) are employed as lower-fidelity models, from which the results are validated in a higher-fidelity ephemeris model. For the computation and optimization of low-thrust trajectories, direct collocation algorithm is introduced. Utilizing the dynamical models and the numerical scheme, the low-thrust trajectory design challenge associated each phase is located and tackled separately. For the heliocentric leg, multiple optimal control problems are formulated between the planets in heliocentric space over different departure and arrival epochs. A contour plot is then generated to illustrate the trade-off between the propellant consumption and the time of flight. For the tour of the Martian moons, the science orbits for both moons are defined. Then, a new algorithm that interfaces the Q-law guidance scheme and direct collocation algorithm is introduced to generate low-thrust transfer trajectories between the science orbits. Finally, an end-to-end trajectory is produced by merging the piece-wise solutions from each phase. The validity of the introduced multi-phase formulation is confirmed by converging the trajectories in a higher-fidelity ephemeris model.<br>

<p>Many opportunities for frequent transit between the lunar vicinity and the heliocentric region will arise in the near future, including servicing missions to space telescopes and proposed missions to various asteroids and other destinations in the solar system. The overarching goal of this investigation is the development a framework for periodic and transit options in the Earth-Moon-Sun system. Rather than overlapping different dynamical models to capture the dynamics of the cislunar and heliocentric region, this analysis leverages a four-body dynamical model, the Bicircular Restricted Four-Body Problem (BCR4BP), that includes the dynamical structures that exist due to the combined influences of the Earth, the Moon, and the Sun. The BCR4BP is an intermediate step in fidelity between the CR3BP and the higher-fidelity ephemeris model. The results demonstrate that dynamical structures from the Earth-Moon-Sun BCR4BP provide valuable information on the flow between cislunar and heliocentric spaces. </p> <p><br></p> <p>Dynamical structures associated with periodic and bounded motion within the BCR4BP are successfully employed to construct transfers between the 9:2 NRHO and locations of interest in heliocentric space. The framework developed in this analysis is effective for transit between any cislunar orbit and the Sun-Earth libration point regions; a current important use case for this capability involves departures from the NRHOs, orbits that possess complex dynamics and near-stable properties. Leveraging this methodology, one-way trajectories from the lunar vicinity to a destination orbit in heliocentric space are constructed, as well as round-trip trajectories that returns to the NRHO after completion of the objectives in heliocentric space. The challenges of such trajectory design include the phasing of the trajectory with respect to the Earth, the Moon, the Sun, on both the outbound and inbound legs of the trajectory. Applications for this trajectory include servicing missions to a space telescope in heliocentric space, where the initial and final locations of the mission is the Gateway near the Moon. Lastly, the results of this analysis demonstrate that the properties and geometry of the periodic orbits, bounded motion, and transfers that are delivered from the BCR4BP are maintained when the trajectories are transitioned to the higher-fidelity ephemeris model. </p>