This thesis explores how an Asteroid Sample Return Mission might make use of solar electric propulsion to send a spacecraft on a journey to the asteroid 1989ML and back. It examines different trajectories that can be used to get an asteroid sample return or similar spacecraft to an interplanetary destination and back in the most fuel-efficient manner. While current plans call for keeping such a spacecraft on the asteroid performing science experiments for approximately 90 days, it is prudent to inquire how lengthening or shortening this time period may affect mission fuel requirements. Using optimal control methods, various mission scenarios have been modeled and simulated. The results suggest that the amount of time that the spacecraft may spend on the asteroid surface can be approximated as a linear function of the available fuel mass. Furthermore, It can be shown that as maximum available thrust is decreased, the radial component of the optimal thrust vector becomes more pronounced.

The Autonomous Multi-Agent Physically Interacting Spacecraft (AMPHIS) test bed examines the problem of multiple spacecraft interacting at close proximity. This thesis contributes to this on-going research by addressing the development of the software architecture for the AMPHIS spacecraft simulator robots and the implementation of a Light Detection and Ranging (LIDAR) unit to be used for state estimation and navigation of the prototype robot. The software modules developed include: user input for simple user tasking; user output for data analysis and animation; external data links for sensors and actuators; and guidance, navigation and control (GNC). The software was developed in the SIMULINK/MATLAB environment as a consistent library to serve as stand alone simulator, actual hardware control on the robot prototype, and any combination of the two. In particular, the software enables hardware-in-the-loop testing to be conducted for any portion of the system with reliable simulation of all other portions of the system. The modularity of this solution facilitates fast proof-of-concept validation for the GNC algorithms. Two sample guidance and control algorithms were developed and are demonstrated here: a Direct Calculus of Variation method, and an artificial potential function guidance method. State estimation methods are discussed, including state estimation from hardware sensors, pose estimation strategies from various vision sensors, and the implementation of a LIDAR unit for state estimation. Finally, the relative motion of the AMPHIS test bed is compared to the relative motion on orbit, including how to simulate the on-orbit behavior using Hill's equations.

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

A lithium vaporizer for a high-power magnetoplasmadynamic (MPD) thruster is modeled using a parallel approach. A one-dimensional, thermal-resistive network is developed and used to calculate the required vaporizer length and power as a function of various parameters. After comparing results predicted by this network model with preheat power data for a 200 kW MPD thruster, we investigate performance over a parameter space of interest for the Advanced Lithium-Fed, Applied-field, Lorentz Force Accelerator (ALFA2) thruster. Heater power sensitivity to cathode tube emissivity, mass flow rate, and vapor superheat are presented. The cold-start heater power for 80 mg/s is found to range from 3.38 to 3.60 kW, corresponding to a vaporizer (axial) length of 18 to 26 cm. The strongest drivers of vaporizer performance are cathode tube emissivity and a conduction heat sink through the mounting flange. Also, for the baseline ALFA2 case, it is shown that increasing the vapor superheat from 100 K to 300 K has the effect of lowering the vaporizer thermal efficiency from 57% to 49%. Also, a finite-volume computational fluid dynamic (CFD) is implemented in FLUENT 6.2 which includes conjugated heat transfer to the solid components of the cathode. This model uses a single-fluid mixture model to simulate the effects of the two-phase vaporizer flow with source terms that model the vaporization. This model provides a solution of higher fidelity by including three-dimensional fluid dynamics such as thermal and momentum boundary layers, as well as calculating a higher resolution temperature distribution throughout the cathode assembly. Results from this model are presented for three mass flow rates of interest (40 mg/s, 80 mg/s, and 120 mg/s). Using a fixed power and length taken from the conceptual ALFA2 design, the dryout point ranges from 12.3-17.6 cm from the base of the cathode assembly for 40 mg/s and 80 mg/s, respectively. For the 120 mg/s case, the two-phase flow never reaches dryout. Finally, results two modeling approaches are compare favorably, with a maximum disagreement of 13.0 percent in prediction of the dryout point and 4.2 percent in predictions of the exit temperature.

"A study is presented detailing the simulation of a drag-free follow-on mission to NASA’s Gravity Recovery and Climate Experiment (GRACE). This work evaluates controller performance, as well as thrust, power, and propellant mass requirements for drag-free spacecraft operation at orbital altitudes of 160 - 225 kilometers. In addition, sensitivities to thermospheric wind, GPS signal accuracy and availability of ephemeris data are studied. Orbital dynamics were modeled in Matlab and take into account 2 body gravity effects, J2-J6 non-spherical Earth effects, atmospheric drag and control thrust. A drag model is used in which the drag acceleration is a function of the spacecraft’s relative velocity to the atmosphere, and a “drag parameter,” which includes the spacecraft’s drag coefficient and local mass density of the atmosphere. A MSISE-90 atmospheric model is used to provide local mass densities as well as free stream flow conditions for a Direct Simulation Monte Carlo drag analysis used to validate the spacecraft drag coefficient. The controller is designed around an onboard inertial sensor which uses a freely floating reference mass to measure deviations in the spacecraft position, resulting from non-gravitational forces, from a desired target orbit. Thruster (control actuator) models are based on two different Hall thrusters for providing the orbital along-track acceleration, colloid thrusters for the normal acceleration, and a miniature xenon ion thruster (MiXI) for the cross-track acceleration. The most demanding propulsion requirements correspond to the lowest altitude considered, 160 kilometers. At this altitude the maximum along-track thrust component is calculated to be 98 millinewtons with a required dynamic (throttling) response of 41 mN/s. The maximum position error at this altitude was shown to be in the along-track direction with a magnitude of 3314.9 nanometers and a peak spectral content of 1800 nm/sqrt(Hz) at about 0.1 Hz. At 225 kilometers, the maximum along-track thrust component reduces to 10.3 millinewtons. The maximum dynamic response at this altitude is 4.23 mN/s. The maximum along-track position error is reduced to 367.9 nanometers with a spectral content peak of 40 nm/sqrt(Hz) at 0.1 Hz. For all altitudes, the maximum state errors increase as the mission length increases, however, higher altitude missions show less of a maximum displacement error increase over time than those of lower orbits. The ability of a colloid thruster to control the normal drift is found to be dependent on how frequently the spacecraft state data is updated. Reducing the period between updates from 10 seconds to 1 second reduces the maximum normal state error component from 199 nanometers to less than 32 nanometers, suggesting that spacecraft state update frequency could be a major driver in keeping the spacecraft on the target trajectory. Sensitivity of maximum required thrust and accumulated sensor error to measurement uncertainty is found to be less of a driver than state update frequency. A ‘worst case” thermospheric wind gust was modeled to show the increase on propulsion requirements if such an event were to occur. At 200 kilometers, maximum winds have been measured to be in increase of 650 m/s in the westward direction in the southern pole region. Assuming the majority of the 650 m/s gust occurs over a 4 second time span, the maximum required cross-track thrust at 200 kilometers increases from 1.12 to 2.01 millinewtons. This large increase may drive the thruster choice for a drag-free mission at a similar altitude. For the spacecraft point design considered with a propellant mass fraction of 0.18, the mission lifetime for the 160 km case was calculated to be 0.76 years. This increases 2.27 years at an altitude of 225 km."

Thesis. 1979. M.S.--Massachusetts Institute of Technology. Dept. of Aeronautics and Astronautics. / MICROFICHE COPY AVAILABLE IN ARCHIVES AND AERONAUTICS. / Bibliography: leaves 86-89. / M.S.

Thesis. 1976. M.S.--Massachusetts Institute of Technology. Dept. of Aeronautics and Astronautics. / MICROFICHE COPY AVAILABLE IN ARCHIVES AND AERONAUTICS. / Includes bibliographical references. / by David M. Stofko. / M.S.

Jia Peng. / Thesis (M.Phil.)--Chinese University of Hong Kong, 2005. / Includes bibliographical references (leaves 85-87). / Abstracts in English and Chinese. / Abstract --- p.iv / Chapter 1 --- Introduction --- p.1 / Chapter 1.1 --- Overview of Output Regulation Problem --- p.2 / Chapter 1.2 --- Attitude Tracking Control of Rigid Spacecraft --- p.3 / Chapter 1.3 --- Overview of Continuous-time Nonlinear H∞ Control --- p.4 / Chapter 1.4 --- Overview of Discrete-time Nonlinear Hq∞ Control --- p.6 / Chapter 1.5 --- Flight Control in Windshears --- p.8 / Chapter 1.6 --- Nonlinear Benchmark System --- p.9 / Chapter 1.7 --- Outline of the Work --- p.11 / Chapter 2 --- Attitude Control and Asymptotic Disturbance Rejection of Rigid Spacecraft --- p.12 / Chapter 2.1 --- Model Description --- p.12 / Chapter 2.2 --- Problem Formulation --- p.16 / Chapter 2.3 --- Preliminaries of General Framework for Global Robust Output Regulation --- p.17 / Chapter 2.4 --- Application of Global Robust Output Regulation --- p.21 / Chapter 2.4.1 --- Case I: without unknown parameters --- p.21 / Chapter 2.4.2 --- Case II: with unknown parameters --- p.26 / Chapter 2.5 --- Simulation --- p.34 / Chapter 2.5.1 --- Case I: without unknown parameters --- p.34 / Chapter 2.5.2 --- Case II: with unknown parameters --- p.36 / Chapter 2.6 --- Conclusions --- p.38 / Chapter 3 --- Application of Approximation Continuous-time Nonlinear H∞ Control Law --- p.45 / Chapter 3.1 --- Preliminaries of Approximation Continuous-time Nonlinear Hq∞ Control Law --- p.45 / Chapter 3.2 --- Disturbance Attenuation of Flight Control System in Windshears --- p.50 / Chapter 3.2.1 --- Design of Control Law --- p.51 / Chapter 3.2.2 --- Computer Simulation --- p.56 / Chapter 3.3 --- Conclusions --- p.57 / Chapter 4 --- Application of Approximation Discrete-time Nonlinear H∞ Control Law --- p.65 / Chapter 4.1 --- Preliminaries of Approximation Discrete-time Nonlinear H∞ Control Law --- p.66 / Chapter 4.2 --- Explicit Expression of u --- p.69 / Chapter 4.3 --- Disturbance Attenuation of RTAC System --- p.73 / Chapter 4.4 --- Computer Simulation --- p.78 / Chapter 4.5 --- Conclusions --- p.80 / Chapter 5 --- Conclusions --- p.83 / Bibliography --- p.85 / A Programs --- p.88 / Vita --- p.112

Modeling the atmospheric entry of spacecraft is challenging because of the large number of physical phenomena that occur during the process. In order to study thermal protection systems, engineers rely on high fidelity solvers to provide accurate predictions of both the thermochemical environment surrounding the heat shield, and its material response. Therefore, it is necessary to guarantee that the numerical models are correctly implemented and thoroughly validated. In recent years, a high-fidelity modeling tool has been developed at the University of Kentucky for the purpose of studying atmospheric entry. The objective of this work is to verify and validate this code. The verification consists of the development of an automated regression testing utility. It is intended to both aid code developers in the debugging process, as well as verify the correct implementation of the numerical models as these are developed. The validation process will be performed through comparison to relevant ablation experiments, namely arc-jet tests. Two modules of the code are used: fluid dynamics, and material response. First the fluid dynamics module is verified against both computational and experimental data on two distinct arc-jet tests. The material response module is then validated against arc-jet test data using PICA.

Testing is an integral part of a satellite’s development, requirements verification and risk mitigation efforts. A robust test program serves to verify construction, integration and assembly workmanship, ensures component, subsystem and system level functionality and reduces risk of mission or capability loss on orbit. The objective of this thesis was to develop a detailed test program for pico-satellites with a focus on the Cal Poly CubeSat architecture. The test program established a testing baseline from which other programs or users could tailor to meet their needs. Inclusive of the test program was a detailed decomposition of discrete and derived test requirements compiled from the CubeSat and Launch Vehicle communities, military guidelines, and industry standards. The test requirements were integrated into a methodical, efficient and risk adverse test flow for verification.