SEARCH · Engineering Papers
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
The details of tha assessment of autonomous Mars ascent and orbit rendezvous for earth return missions are presented. Analyses addressing navigation system assessments, trajectory planning, targeting approaches, flight control guidance strategies, and performance sensitivities are included. Tradeoffs in the analysis and design process are discussed.
When charged particles spiral along a large constant magnetic field, their trajectories are scattered by random components that are superposed on the guiding field. In the simplest analysis of this situation, scattering causes the particles to diffuse parallel to the guiding field. At the next level of approximation, moving pulses that correspond to a coherent mode of propagation are present, but they are represented by delta-functions whose infinitely narrow width makes no sense physically and is inconsistent with the finite duration of coherent pulses observed in solar energetic particle events. To derive a more realistic description, the transport problem is formulated in terms of 4 x 4 matrices, which derive from a representation of the particle distribution function in terms of eigenfunctions of the scattering operator, and which lead to useful approximations that give explicit predictions of the detailed evolution not only of the coherent pulses, but also of the diffusive wake. More specifically, the new description embodies a simple convolution of a narrow Gaussian with the solutions above that involve delta-functions, but with a slightly reduced coherent velocity. The validity of these approximations, which can easily be calculated on a desktop computer, has been exhaustively confirmed by comparison with results of Monte Carlo simulations which kept track of 50 million particles and which were carried out on the Maspar computer at Goddard Space Flight Center.
An integrated methodology has been assembled to compute the engine performance, takeoff and landing trajectories, and community noise levels for a subsonic commercial aircraft. Where feasible, physics-based noise analysis methods have been used to make the results more applicable to newer, revolutionary designs and to allow for a more direct evaluation of new technologies. The methodology is intended to be used with approximation methods and risk analysis techniques to allow for the analysis of a greater number of variable combinations while retaining the advantages of physics-based analysis. Details of the methodology are described and limited results are presented for a representative subsonic commercial aircraft.
The purpose of this paper is to study the dynamic behavior of a three-wave journal bearing using a transient approach. The transient analysis permits the determination of the rotor behavior after the fractional frequency whirl appears. The journal trajectory is determined by solving a set of nonlinear equations of motion using the Runge-Katta method. The fluid film forces are computed by integrating the transient Reynolds equation at each time step location of the shaft with respect to the bearing. Because of the large values of the rotational speeds, turbulent effects were included in the computations. The influence of the temperature on the viscosity was also considered. Numerical results were compared to experimenta1 results obtained at the NASA Glenn Research Center. Comparisons of the theoretical results with experimental data were found to be in good agreement. The numerical and experimental results showed that the fluid film of a three-wave journal bearing having a diameter of 30 mm, a length of 27 mm, and a wave amplitude ratio greater than 0.15 is stable even at rotational speeds of 60,000 RPM. For lower wave amplitude ratios, the threshold speed at which the fluid film becomes unstable depends on the wave amplitude and on the supply pocket pressure. Even if the fluid film is unstable, the wave bearing maintains the whirl orbit inside the bearing clearance.
Dragonfly is a proposed spacecraft and mission that would send a mobile robotic rotorcraft lander to Titan, the largest moon of Saturn, in order to study prebiotic chemistry and extraterrestrial habitability at various locations. Titan is unique in having an abundant, complex, and diverse carbon-rich chemistry on the surface of a water-ice-dominated world with an interior water ocean, making it a high priority target for astrobiology and origin of life studies. The mission was initially proposed in April 2017 to NASA's New Frontiers program by the Johns Hopkins Applied Physics Laboratory. In December 2017, it was selected as one of two finalists (out of twelve proposals) to further refine the mission's concept. NASA Ames Research Center and NASA Langley Research Center are partnering as the leads for the Dragonfly's entry system to provide the completed EDL Assembly. The aerothermal analysis for Dragonfly utilizes four simulation tools from NASA Ames Research Center. Traj for calculating the trajectory, DPLR 4.04.0 for calculating the flowfield around the vehicle and convective heating, NEQAIR V15.0 for calculating the radiative heating, and FIAT for calculating the material response and thermal protection system (TPS) sizing for the heatshield. The entry conditions are relatively benign and can readily be accommodated with a tiled PICA heatshield similar to MSL and a number of flight proven materials for the backshell. This work will demonstrate that the aerothermal entry environments can be readily solved using heritage materials and techniques.
Development of Detect-and-avoid (DAA) systems’ maneuver guidance requirements at RTCA has centered on tactical maneuvers away from the intruder. Recapturing the flight-plan route the initial maneuver is not directly taken into account by most DAA maneuver guidance algorithms. This work demonstrates potential challenges and inefficiencies that can arise in recapturing the route after resolving a conflict. Horizontal resolution trajectories for a test matrix of encounters with varying aircraft speeds and geometries are computed by minimizing either flight time or deviation to serve as a baseline. Turn directions computed from a reference DAA algorithm coupled with a pilot selection model and a second algorithm called Autoresolver (AR) are compared to these baseline resolution trajectories. Baseline results show that turning into the intruder yields favorable cost for most encounters. Additional analysis of pilot response data shows that only three-fourths of pilots’ horizontal maneuvers turn into intruders, a percentage much lower than the baseline results. Improvement to a DAA guidance algorithm based on findings in this work is discussed.
Development of Detect-and-avoid (DAA) systems’ maneuver guidance requirements at RTCA has centered on tactical maneuvers away from the intruder. Recapturing the flight-plan route the initial maneuver is not directly taken into account by most DAA maneuver guidance algorithms. This work demonstrates potential challenges and inefficiencies that can arise in recapturing the route after resolving a conflict. Horizontal resolution trajectories for a test matrix of encounters with varying aircraft speeds and geometries are computed by minimizing either flight time or deviation to serve as a baseline. Turn directions computed from a reference DAA algorithm coupled with a pilot selection model and a second algorithm called Autoresolver (AR) are compared to these baseline resolution trajectories. Baseline results show that turning into the intruder yields favorable cost for most encounters. Additional analysis of pilot response data shows that only three-fourths of pilots’ horizontal maneuvers turn into intruders, a percentage much lower than the baseline results. Improvement to a DAA guidance algorithm based on findings in this work is discussed.
This paper presents analysis of when and how well, on the approach to previously unvisited small bodies, the trajectory-design relevant physical parameters are determined.
The capabilities of a SmallSat-class spacecraft targeting the outer solar system and using a combined chemical and electric propulsion system are explored. The development of compact hybrid rockets has enabled high-thrust engines to be packaged tightly enough to fit on cubesat and SmallSat spacecraft. These hybrid rockets provide 10’s-100 N of thrust depending on the propellant load & >300 s of specific impulse and have been demonstrated in both ambient and vacuum environments. Advancements in low-power long-life Hall thruster technologies have provided the potential for significantly greater propellant throughputs, enabling their use as a primary propulsion element on interplanetary spacecraft. In a recent characterization test campaign, the MaSMi-DM Hall thruster demonstrated power throttling from 150 – 1000 W with >1500 s of specific impulse available at >500 W and >40% total thrust efficiency available at >300 W; peak values of 1940 s and 53% were observed. A notional low-mass spacecraft employing a combined hybrid rocket and low-power electric propulsion system was designed and used for mission concept analysis targeting the outer solar system. Using an imposed wet mass limit of 400 kg, mission trajectories to Saturn and Uranus were generated. Orbit capture with >40% of the launch mass was shown to be possible at either target, with mission transfer times of 7.5 years and 13.5 years for Saturn and Uranus, respectively. Significant follow-on mission activities near Saturn (e.g. to Titan & Enceladus) were also possible by carrying extra propellant mass while remaining under the total wet mass limit.
The Origins, Spectral Interpretation, Resource Identification, and Security–Regolith Explorer (OSIRIS-REx) mission collected a sample from the surface of the near-Earth asteroid (101955) Bennu in late 2020. Bennu challenged the team with a surface that was much rockier than expected, resulting in modifications to the prelaunch design of the Touch And Go (TAG) sequence. Following enhancements in onboard trajectory correction, ground-based navigation, and maneuver execution error modeling, the spacecraft was delivered to the chosen TAG site within 1 m of the target, and a sample was successfully collected on the first attempt. This paper provides a comprehensive description of all flight dynamics aspects of TAG trajectory planning and execution. It also describes hazard map generation and how that combined with error analysis results to predict the probability of safe contact before TAG and the onboard wave-off determination during TAG.
Sinai-Ruelle-Bowen (SRB) measures are limiting stationary distributions describing the statistical behavior of chaotic dynamical systems. Directional derivatives of SRB measure densities conditioned on unstable manifolds are critical in the sensitivity analysis of hyperbolic chaos. These derivatives, known as the SRB density gradients, are by-products of the regularization of Lebesgue integrals appearing in the original linear response expression. In this paper, we propose a novel trajectory- driven algorithm for computing the SRB density gradient defined for systems with high-dimensional unstable manifolds. We apply the concept of measure preservation together with the chain rule on smooth manifolds. Due to the recursive one-step nature of our derivations, the proposed procedure is memory-efficient and can be naturally integrated with existing Monte Carlo schemes widely used in computational chaotic dynamics. Here, we numerically show the exponential convergence of our scheme, analyze the computational cost, and present its use in the context of Monte Carlo integration.
The fuel tank isolation as a result of leak or rupture can leave an Orion Crew Module with only half of the loaded propellant for ISS return atmospheric entry. To assess the feasibility of returning under this condition, an analysis of various entry control options with deliberate degradation of control performance was performed. The study determined that a ballistic entry without a raise burn, a steeper flight path trajectory, relaxed atmospheric pitch/yaw rate damping, and degraded touchdown control could achieve 2-sigma requirements compliance with a 2-sigma fuel usage that is less than half of the liftoff propellant loading. The results of this analysis indicate that an entry with only half the nominal propellant load is feasible following a tank isolation fault.
Trajectory control, landing, surface navigation, and path selection studies for Mars exploration
On January 15th, 2006, 0556 UTC, the Stardust spacecraft released a 42-kg Sample Return Canister (SRC) along a trajectory intended to impact a target on the Air Force Utah Test and Training Range (UTTR), near Dugway, Utah. Assurances of a successful SRC delivery to UTTR depended on identifying (and mitigating, if possible) a myriad of error sources. These sources included atmospheric effects, maneuver execution, Orbit Determination (OD) uncertainties and (Delta)V induced by the firing of the unbalanced Reaction Control System (RCS) thrusters needed for deadband attitude control. Every mm/s in prediction error at the TCM-19 epoch would amount to missing the target by approximately one kilometer. This paper will describe the work performed in analyzing and predicting the levels of (Delta)V caused by the attitude deadbanding, as well as prediction performance.
A pair of human-in-the-loop simulation evaluations of a distributed air/ground separation assurance system have been conducted to investigate the function allocation between humans and automation systems as well as ground-based and airborne agents in the Next Generation Air Transportation System and beyond. This paper focuses on an analysis of certain critical conflicts observed between self-separating aircraft and ground-managed traffic in the same airspace. The principal cause of each conflict is identified and potential mitigations are discussed, such as: the sharing of trajectory intent information between the ground and the air; more cautious trajectory planning by the self-separating aircraft; and more equitable rules-of-the-road between the self-separating aircraft and ground-managed aircraft. This analysis will inform the ongoing design of an airborne separation assurance automation tool.
The development of guidance approximations for the atmospheric (aeroglide) portion of the minimum fuel, orbital plane change, trajectory optimization problem is described. Asymptotic methods are used to reduce the two point, boundary value, optimization problem to a turning point problem from the bank angle control. The turning point problem solution, which yields an approximate optimal control policy, is given in terms of parabolic cylinder functions, which are tabulated, and integral expressions, which must be numerically computed. Comparisons of the former, over their region of validity, with optimal control solutions show good qualitative agreement. Additional work and analysis is needed to compute the guidance approximation work.
The NASA Conjunction Assessment Risk Analysis (CARA) team has recently implemented updated software to calculate the probability of collision (P (sub c)) for Earth-orbiting satellites. The algorithm can employ complex dynamical models for orbital motion, and account for the effects of non-linear trajectories as well as both position and velocity uncertainties. This “3D P (sub c)” method entails computing a 3-dimensional numerical integral for each estimated probability. Our analysis indicates that the 3D method provides several new insights over the traditional “2D P (sub c)” method, even when approximating the orbital motion using the relatively simple Keplerian two-body dynamical model. First, the formulation provides the means to estimate variations in the time derivative of the collision probability, or the probability rate, R (sub c). For close-proximity satellites, such as those orbiting in formations or clusters, R (sub c) variations can show multiple peaks that repeat or blend with one another, providing insight into the ongoing temporal distribution of risk. For single, isolated conjunctions, R (sub c) analysis provides the means to identify and bound the times of peak collision risk. Additionally, analysis of multiple actual archived conjunctions demonstrates that the commonly used “2D P (sub c)” approximation can occasionally provide inaccurate estimates. These include cases in which the 2D method yields negligibly small probabilities (e.g., P (sub c)) is greater than 10 (sup -10)), but the 3D estimates are sufficiently large to prompt increased monitoring or collision mitigation (e.g., P (sub c) is greater than or equal to 10 (sup -5)). Finally, the archive analysis indicates that a relatively efficient calculation can be used to identify which conjunctions will have negligibly small probabilities. This small-P (sub c) screening test can significantly speed the overall risk analysis computation for large numbers of conjunctions.
Icarus, an in-house material response solver developed at NASA Ames, is applied to perform 2-D bondline-temperature analysis on the heatshield surface of the latest Dragonfly entry system design for a mission to Titan. Two axisymmetric domains are used: the near-shoulder region and the entire heatshield. The heatshield consists of a PICA TPS bonded to a layered system of Aluminum honeycomb wrapped by an M55J carbon-fiber facesheet. The backshell, which is also included for more accurate analysis, consists of similar materials layered. Icarus simulations using both domains indicate that the in-plane thermal conductivity property of the M55J facesheet is observed to play a dominant role in bondline temperature. Similarly, the maximum bondline temperature is found around the PICA-tile-interface region instead of the near-shoulder region or the stagnation point based on the current trajectory using the orthotropic properties of M55J. Therefore, these findings indicate the significance of 2-D or higher-dimensional material modeling analysis in fully understanding how bondline temperature behaves across the heatshield and obtaining a basis for TPS design locations. Verification at different radial locations ahead of the shoulder, where thermal conduction is close to 1-D, between Icarus (2-D), Icarus (1-D), and FIAT is performed and good agreement is observed. With the entire heatshield domain, the total mass loss of the aeroshell materials is approximately 3% of their initial mass potentially due to thermal decomposition of phenolic resins inside PICA under the current high-heating environment.