Browsing by Author "Fitzgerald, Riley McCrea"
Now showing 1 - 12 of 12
Results Per Page
Sort Options
- Analysis of Low-Energy Lunar Transfers in a High-Fidelity Dynamics ModelTorchia, Patrick Jason (Virginia Tech, 2023-07-03)Renewed interest in returning to the Moon, emboldened by recent directives and missions by NASA, has necessitated the establishment of lunar infrastructure to support continuous human presence. With that, the objective of making this return more cost effective has gained significant importance. Low energy lunar transfers are more efficient ways to reach the Moon than the traditional Hohmann-type transfer. These trajectories leverage the multi-body gravitational effects to reduce overall delta-v requirements, in some cases removing the capture delta-v completely. While the time of flight for these transfers can be much longer than a Hohmann-type transfer, the chaotic design space of these transfers can enable large changes in arrival conditions at the Moon for small changes in initial conditions. Many investigations of these transfers take place in simplified dynamical models, such as the Planar Circular Restricted Three Body Problem, with very few higher-fidelity models being implemented. This approach is good to understand the dynamics of these trajectories as well as provide initial guesses for higher-fidelity models; but approximating the dynamics heavily make these models less applicable to mission design. This thesis aims to investigate the application of a higher-order model to simulate these trajectories. STK Astrogator was used to recreate the NASA GRAIL trajectory; and from the recreated trajectory, a nominal trajectory absent of mid-course corrections was established. This nominal trajectory was used to perform parametric and variational studies of departure and arrival conditions as well as compare to a nominal trajectory in a reduced-fidelity model. An investigation into the post launch correction burn requirements following launch vehicle under-performance was completed. Utilizing low energy transfers proved beneficial to adjusting arrival conditions for low delta-v requirements. All arrival inclinations are reasonably achievable for around 255 m/s. Using 255 m/s as a baseline, right ascension of the ascending node could be reached in a 40 degree range and argument of periapsis in a 50 degree range. Lunar insertion arrival can be varied by 7 hours on either side for less than 80 m/s. Trans-lunar injection epoch can be varied by 7 hours on either side of nominal departure for less than 4 m/s. Orbit radius and initial velocity are the most expensive errors to correct. These trajectories can be tuned to reduce the overall mid-course correction delta-v requirement for differing arrival inclinations if other orbital elements are relaxed. A relationship between placement of post-launch correction maneuver for velocity or radius errors was found. Comparing the trajectory in STK to the Inclined Bi-Elliptic Restricted Four Body problem, revealed that timing of the trajectory is variable while keeping the same arrival and departure conditions. However, solar radiation pressure cannot be ignored for more accurate simulation of these trajectories. This investigation has shown that low energy lunar transfers are a viable method to reach the Moon and their chaotic nature can be leveraged to relax restrictions in the design space.
- Analysis of Transfer Trajectories Utilizing Sequential Saturn-Titan AerocapturesPayne, Isaac Lee (Virginia Tech, 2023-07-03)This thesis aims to investigate the potential of a transfer orbit using successive aerocaptures at Saturn and Titan to establish a science orbit around Titan. Titan is an Earth-like moon with a dense atmosphere and organic compounds present. It has many similarities with Earth that are useful to study such as superrotation. Superrotation is when the atmosphere rotates faster than the body it surrounds. In order to study Titan, we need to establish an orbit around it. The Saturn system is distant from Earth, 8.5 Astronomical Units (AU) which makes it difficult to reach from a time and velocity point of view. We propose to use an aerocapture at Saturn to intercept Titan with lower relative velocity in order to perform an aerocapture at Titan. The analysis was performed in primarily MATLAB to simulate the orbits. The results of this showed that we can aerocapture a spacecraft at Saturn and arrive at Titan within roughly 4 to 8 km/s relative velocity regardless of the incoming hyperbolic excess velocity at the Saturn system. This can be improve upon by using intermediate transfer orbits, such as bi-elliptics, to arrive with even lower relative velocities to Titan of as low as 1 km/s. The drag acceleration experienced during the Saturn aerocapture had peak values of between 0.2 and 1.4 g's and acceleration over 50% of the peak is experienced between 6.8 and 8 minutes. This capture method has the potential to make Titan more easily accessible and allow for scientific study of a clear target for improving our understanding of Earth-like processes on other bodies in our solar system.
- Boundary Resilience: A New Approach to Analyzing Behavior in Complex SystemsWilhelm, Julia Claire Wolf (Virginia Tech, 2024-04-30)Systems engineering has many subdisciplines which would be useful to study in terms of complex system behavior. However, it is the interactions between a complex system and its operating environment which drive the motivation for this analysis. Specifically, this work introduces a new approach to assessing these interactions called "boundary resilience." While classical resilience theory measures a system's internal reaction to adverse event, boundary resilience evaluates the impacts such an event may have on the surrounding environment. As the scope of this analysis is quite large, it was deemed appropriate to conduct a case study to determine the fundamental tenants of boundary resilience. SpaceX's satellite Internet mega-constellation (StarLink) was chosen due to its large potential to impact the space environment as well as its size and complexity. This study produced two boundary resilience measures, one for local boundary resilience of a single component and one for the global boundary behavior of the entire system. The local metric measures the likelihood of an adverse event occurring at that boundary location as well as its potential to impact the surrounding environment. The global boundary resilience metric reflects a nonlinear relationship among the system components.
- Evaluating the Capability of ICON-MIGHTI to Detect Plasma Bubbles in the IonosphereLech, Brenden (Virginia Tech, 2024-12-09)The MIGHTI airglow imager onboard the ICON spacecraft in LEO was built to make remote thermospheric windspeed measurements at low latitudes. The MIGHTI team, when reviewing the data, observed variations in day-to-day brightness potentially indicative of plasma bubbles: regions of low-density E-region plasma which rise through the F-region and cause radio scintillation that interferes with communications and GPS performance. Here, we explore the possibility of MIGHTI observing plasma bubbles by using its red-line airglow measurements to attempt to detect this phenomenon. Small-scale structuring indicative of plasma bubbles is searched for by comparing measurements between MIGHTI's two identical imagers, which make remote airglow measurements at the same region from perpendicular directions. The usability of the two imagers for this purpose is assessed, given they are not calibrated to measure absolute airglow brightness, and it is determined that the level of disagreement between them does not prevent these comparisons. The evolution of the ionosphere in the time between the two instruments' measurements is accounted for using seasonal medians of expected behavior. Co-located measurements where the two MIGHTI imagers disagreed significantly were found, filtering out disagreements in measurement not likely to have a significant underlying ionospheric cause, although none were indicative of plasma bubble observations. These significantly differing measurements were most common shortly after dusk and in regions near the equator, especially between -30 to 70 degrees longitude. Simulations show the lack of definitive plasma bubble detections is likely due to MIGHTI's long image exposure time averaging out the effect of plasma bubbles as ICON orbits. More is now known about the potential for making comparative red-line airglow measurements between MIGHTI's imagers, and this information could be used in future work to explore larger-scale ionospheric structuring within the MIGHTI dataset.
- Lunar Laser Ranging for Autonomous Cislunar Spacecraft NavigationZaffram, Matthew (Virginia Tech, 2023-08-15)The number of objects occupying orbital regimes beyond Geosynchronous Earth Orbit and cislunar space are expected to grow in the coming years; Especially with the Moon reemerging as latest frontier in the race for space exploration and technological superiority. In order to support this growth, new methods of autonomously navigating in cislunar space are necessary to reduce demand and reliance on ground based tracking infrastructure. Periodic orbits about the first libration point offer favorable vantage points for scientific or military spacecraft missions involving the Earth or Moon. This thesis develops a new autonomous spacecraft navigation method for cislunar space and analyzes its performance applied to Lyapunov and halo orbits around $L_1$. This method uses existing lunar ranging retroreflectors (LRRR) installed on the Moon's surface in the 1960s and 1970s. A spacecraft can make laser ranging measurements to the LRRR to estimate its orbit states. A simulation platform was created to test this concept in the circular restricted three body problem and evaluate its performance. This navigation method was found to be successful for a subset of Lyapunov and halo orbits when cycling the five measurement targets. Simulation data showed that sub-kilometer position estimation and sub 2 centimeter per second velocity accuracies are achievable without receiving any state updates from external sources.
- Modeling and Analysis of a Thermospheric Density Measurement System Based on Torque EstimationAceto, Christopher James (Virginia Tech, 2023-07-12)This thesis models and analyzes an in-situ method for measuring the density of the thermosphere at low Earth orbit (LEO) altitudes in real time. As satellites orbit in the thermosphere, the sparse yet present air perturbs their orbits via the drag force. The drag force is poorly characterized and has a significant effect at LEO altitudes relative to other forces, making this perturbation force one of the greatest uncertainties in LEO orbit propagation. A steadily increasing number of satellites orbit at LEO altitudes, so for safety, it is critical to accurately track these satellites to avoid collisions. Therefore, better knowledge of the drag force is required. The drag force depends directly on the air mass density in the thermosphere, and current knowledge of the thermospheric density is limited. Models exist to describe the variations in density over time, but due to the many unpredictable factors which affect the thermosphere, the best of these models are only accurate to within 10%. Also, currently available techniques to measure the thermospheric density can only return time-averaged measurements, which causes inaccuracies in orbit propagation due to local density variations. Some planned in-situ density measurement missions rely on measuring acceleration caused by the drag force, but this requires a highly accurate accelerometer to be able to separate the drag force from other stronger forces acting on a satellite. The Satellite Producing Aerodynamic Torque to Understand LEO Atmosphere (SPATULA) concept was introduced as an alternative method, which infers density based on measurements of the drag torque. In the rotational regime, drag produces the strongest torque at LEO altitudes by far, making it possible to acquire accurate density measurements with inexpensive, commercially available sensors and actuators on a SPATULA spacecraft. This thesis expands upon a preliminary study of the SPATULA concept. A SPATULA spacecraft's dynamics are modeled in three dimensions, and a novel method is introduced for modeling the dependence of external torques on the geometry and attitude of the spacecraft. In addition to the dynamics model, discrete-time algorithms for guidance, system state filtering, attitude control, and density estimation are developed for the six degrees of freedom case. The MathWorks tools MATLAB and Simulink are used to simulate the physics and system models. The simulations are used to evaluate the performance of the SPATULA system's density measurements and compare them to conventional methods. It is found that the accuracy and bandwidth of the SPATULA system have a significant dependence on the assumed accuracy of the torque models in the system's filter. When the bandwidth is set to avoid significant phase shift errors, the SPATULA system can produce real-time measurements of density accurate over a minimum time scale of about 60 seconds, and the density error has a standard deviation of about 2 x 10^-14 kg/m^3. This accuracy is about 6 times better than the best thermospheric models, and it is also better than reported accuracies of most other density measurement methods. If bandwidth is sacrificed, the density error standard deviation can be decreased by a factor of 4. This introduces additional error due to phase shift delays, but these can be corrected with signal processing techniques. With the higher accuracy, the SPATULA system loses its real-time ability, but the data it produces would still provide excellent insight for improving thermospheric models. With high accuracy and low cost, the SPATULA concept is a promising path to pursue toward improving thermospheric density knowledge.
- Multiple Gravity Assists for Low Energy Transport in the Planar Circular Restricted 3-Body ProblemWerner, Matthew Allan (Virginia Tech, 2022-06-23)Much effort in recent times has been devoted to the study of low energy transport in multibody gravitational systems. Despite continuing advancements in computational abilities, such studies can often be demanding or time consuming in the three-body and four-body settings. In this work, the Hamiltonian describing the planar circular restricted three-body problem is rewritten for systems having small mass parameters, resulting in a 2D symplectic twist map describing the evolution of a particle's Keplerian motion following successive close approaches with the secondary. This map, like the true dynamics, admits resonances and other invariant structures in its phase space to be analyzed. Particularly, the map contains rotational invariant circles reminiscent of McGehee's invariant tori blocking transport in the true phase space, adding a new quantitative description to existing chaotic zone estimates about the secondary. Used in a patched three-body setting, the map also serves as a tool for investigating transfer trajectories connecting loose captures about one secondary to the other without any propulsion systems. Any identified initial conditions resulting in such a transfer could then serve as initial guesses to be iterated upon in the continuous system. In this work, the projection of the McGehee torus within the interior realm is identified and quantified, and a transfer from Earth to Venus is exemplified.
- Preliminary Design of a Titan-Orbiting Stellar Occultation MissionWagner, Nathan John (Virginia Tech, 2022-06-09)This thesis serves to provide a conceptual mission design for a Titan-orbiting stellar occultation mission. Titan has a significant atmosphere much like Earth's. An improved understanding of Titan's atmosphere could provide valuable information about the evolution of Earth's climate. Titan's atmosphere is known to be in a state of superrotation, wherein the atmosphere rotates significantly faster than the surface beneath. The details of the creation and sustainment of this extreme state on Titan in terms of angular momentum exchange remain unknown despite current theories and models. These unknowns, alongside inconsistencies between current models with observations from the Cassini mission, call for an urgent need for Titan atmospheric observations able to resolve atmospheric waves. The science objectives driving the mission design include maximizing the number of measurements, the latitude versus longitude coverage, the latitude versus local solar time coverage, and the mission duration. These measurement needs can be met by a Titan orbiter utilizing a refractive stellar occultation technique. Refractive stellar occultation observes starlight bending through an atmosphere as stars set behind a body. The observed bending profile can be inverted to infer density, temperature, and pressure profiles. This research uses Systems Tool Kit (STK) as a simulation tool to predict measurement coverage for various orbits. The orbital radius was determined to be the driving independent variable which set all other design variables, including the orbital plane which was uniquely selected for a given orbital radius to maximize the number of occultations. The results of this study show that a lower orbital radius is desired as this produces the best combination of measurement number and distribution. This orbital plane should be closely aligned with the Milky Way galactic plane to see the most stars occult. For the lowest sustainable orbital altitude, Low Titan Orbit (LTO) at 1200 km, the orbital plane should be nearly polar to maximize the number of occultations and latitude coverage. The optimal orbit selection (defined by orbital elements a = 3775 km, e = 0, i = 85 degrees, Ω = 87 degrees, ω = 0 degrees, and ν = 0 degrees) for a single satellite can produce nearly 400 stellar occultation opportunities per orbit and provide full latitude versus longitude coverage. A single satellite shows gaps in latitude versus local solar time coverage at mid-latitudes normal to the satellite ground track which may inhibit the diagnosis of the angular momentum flux associated with thermal tides. If necessary, a second satellite in an orbit orthogonal to the first is suggested to close coverage gaps to provide full local time coverage over a Titan day. The optimal orbit selection of this second satellite (defined by orbital elements a = 3775 km, e = 0, i = 5.3 degrees, Ω = 5.9 degrees, ω = 0 degrees, and ν = 0 degrees) provides an additional 343 occultation opportunities per orbit and increases latitude versus local solar time coverage by a factor of 1.5. The understanding of Titan's Earth-like atmosphere could provide insight into climate evolution here on Earth. This concept proposes a novel approach to improving this understanding.
- Reactive, Autonomous, Markovian Sensor Tasking in Communication Starved EnvironmentsKadan, Jonathan Evan (Virginia Tech, 2024-01-02)The current Space Traffic Management (STM) community was not prepared for the exponential increase in the resident space object (RSO) population that has taken place over the last several years. The combination of poor communication infrastructure and long scheduling lead times of the Space Surveillance Network (SSN) prevent any type of reactive sensor tasking, which is required in event of anomaly detection. This dissertation was designed to survey extensions to the classical notions of covariance based sensor tasking strategies and develop a methodology for evaluating these techniques. A suboptimal partially observable Markov decision process (POMDP) was used as the simulation framework to test various reward functions and decision algorithms while enabling autonomous, reactive sensor tasking. The goal of this work was used the developed evaluation methodology to perform statistical analyses to determine which metrics were most reliable and efficient for Space Traffic Management (STM) of the geosynchronous Earth orbit (GEO) resident space object (RSO) catalog. Hypotheses were tested against simulations of 873 resident space object (RSO) in geosynchronous Earth orbit (GEO) being tracked by 18 heterogeneous, geographically disperse ground-based electro-optical (EO) sensors. This dissertation evaluates the ability of various sensor tasking metrics to produce rewards that maximize geosynchronous Earth orbit (GEO) catalog coverage capability of a sensor network under realistic communication restrictions.
- Using a Curvilinear Coordinate System for Satellite Relative MotionMidas, Alex Matthew (Virginia Tech, 2024-02-23)The number of dynamics needed to model the motion between a Chief and Deputy satellites has grown greatly since the introduction of the Hill, Clohessy-Wilshire (HCW) equations of motion were introduced. The models have grown to include various things like perturbations, specifically drag, J2, and solar radiation pressure. Dynamics models have also been developed that use True Anomaly as the independent variable instead of time. A lot of work has been put forth to also include cases where the Chief is in an eccentric orbit. While these models have increased the fidelity of relative dynamics these models become very complicated to implement. That is why the HCW equations remain extremely popular after all these developments. However, their simplicity causes issues when there is In-Track separation between the Chief and Deputy satellites. The error in the dynamics increases as this separation increases which leads to a typical constraint that the separation between the Chief and Deputy needs to be much smaller than the radius of the Chief's orbit. That is where this works starts, by examining into ways to increase the accuracy in the HCW equations as the In-Track separation between the Chief and Deputy grows. In which, this will be done by using a curvilinear coordinate system. Furthermore, a technique of using a Virtual Chief satellite will by employed to allow for the HCW equations to be valid for cases where the Chief is in an eccentric orbit.
- Using the Circular Restricted Three-Body Problem to Design an Earth-Moon Orbit Architecture for Asteroid MiningMunson Jr., Mark Allan (Virginia Tech, 2024-06-05)Engineering and technical challenges exist with the material transport of natural resources in space. One aspect of this transport problem is the design of an orbit architecture in the Earth-Moon system (EMS) that facilitates these resources through the mining cycle. In this thesis, it is proposed to use the Circular Restricted 3-Body Problem (CR3BP) to design an orbit architecture composed of L3 Lyapunov orbits, hyperbolic invariant stable and unstable manifolds, and geosynchronous (GEO) orbits. A single shooting method (SSM) and natural parameter continuation (NPC) numerical algorithm is used to compute a family of L3 Lyapunov orbits. Invariant Manifold Theory (IMT) is leveraged to find the set of feasible hyperbolic invariant stable and unstable manifolds associated with a L3 Lyapunov orbit. Ideal L3 Lyapunov orbits are chosen to construct an orbit architecture based off favorable metrics like orbital period, Jacobi Constant, and stability index. Manifolds that enter the GEO and xGEO (beyond GEO) volumes are identified. Finally, a ∆V analysis for GEO to manifold transfer is conducted. An achievement of this study is the computation of stable L3 Lyapunov orbits. The primary contribution of this paper lies in its modeling of a L3 Lyapunov orbit architecture using the CR3BP.
- A Viable Orbital Debris Mitigation Mission using Active Debris RemovalSmeltzer, Stanley Logan (Virginia Tech, 2023-06-28)Currently, the Low Earth Orbit (LEO) space environment contains a growing number of orbital debris objects. This growing orbital debris population increases collision probabilities between both orbital debris and functioning satellites. A phenomenon known as Kessler Syndrome can be induced if these collisions occur. Kessler Syndrome states that these collisions can lead to an exponential increase in the orbital debris population, which could dangerously impede future space missions. Current literature outlines the necessity of stabilizing the near-Earth environment debris population and introduces the concept of active debris removal (ADR). The use of ADR on five orbital debris objects per year was found to be a requirement to achieve stability within the orbital debris population. A viable mission architecture is henceforth explored to utilize ADR for near-future execution to further develop research for orbital debris mitigation missions. The larger orbital debris objects are found in many different orbital regimes and are primarily composed of spent rocket bodies and retired satellites. Different orbital debris ranking schemes have been developed based on the population density in these different regimes, which are linked to higher collision probabilities. Using these ranking schemes, a set of target objects are selected to be investigated for this mission design that was composed of target objects with similar orbital characteristics that were not launched by the Commonwealth of Independent States (CIS) to minimize legal barriers. Different ADR capture and removal methods are inspected to find the optimal methods for this mission. An Analytical Hierarchy Process (AHP) has been used to assess these different methods, which utilizes comparisons of the different methods among a set of weighted criteria. A net capture method with a low thrust chemical engine for removal is identified as the optimal ADR method. The use of a laser detumbling system is also selected to stabilize target objects with a high rotation rate. A rendezvous and deorbit orbital analysis are conducted using both a low fidelity tool (for preliminary results) and a high fidelity tool (for more precise results). The rendezvous analysis is used to select a mission architecture that was composed of two different chaser satellites which rendezvous with the five different target objects by taking advantage of nodal precession. The deorbit analysis investigates different decay timelines and found the delta-v estimates that would be required to deorbit the target objects within the same year that they were captured in. These two orbital analyses provide valuable insight to the mission timeline, delta-v estimates, and approximate mass requirement for the chaser satellite and deorbit kits. The results of the target selection process, ADR selection process, and the rendezvous and deorbit analyses are meant to provide an initial concept and analysis for a near-future ADR mission. These approximate results provide insight and information to further develop orbital debris mitigation research to help solve the orbital debris population growth challenge for future space missions.