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.
Computer program includes time minimizing function required for propellant constrained minimum time impulsive transfer between coplanar circular orbits or unrestricted orbits. Basic formulation is extended to include inclined circular or inclined elliptic orbits.
Computer generated graphics are developed to assist in the modeling and assessment of pointing and control systems of orbiting spacecraft. Three-dimensional diagrams are constructed of the Earth and of geometrical models which resemble the spacecraft of interest. Orbital positioning of the spacecraft model relative to the Earth and the orbital ground track are then displayed. A star data base is also available which may be used for telescope pointing and star tracker field-of-views to visually assist in spacecraft pointing and control studies. A geometrical model of the Hubble Space Telescope (HST) is constructed and placed in Earth orbit to demonstrate the use of these programs. Simulated star patterns are then displayed corresponding to the primary mirror's FOV and the telescope's star trackers for various telescope orientations with respect to the celestial sphere.
Computational aerothermodynamic simulations of Orbiter windside tile damage in flight were performed in support of the Space Shuttle Return-to-Flight effort. The simulations were performed for both hypervelocity flight and low-enthalpy wind tunnel conditions and contributed to the Return-to-Flight program by providing information to support a variety of damage scenario analyses. Computations at flight conditions were performed at or very near the peak heating trajectory point for multiple damage scenarios involving damage windside acreage reaction cured glass (RCG) coated silica tile(s). The cavities formed by the missing tile examined in this study were relatively short leading to flow features which indicated open cavity behavior. Results of the computations indicated elevated heating bump factor levels predicted for flight over the predictions for wind tunnel conditions. The peak heating bump factors, defined as the local heating to a reference value upstream of the cavity, on the cavity floor for flight simulation were 67% larger than the peak wind tunnel simulation value. On the downstream face of the cavity the flight simulation values were 60% larger than the wind tunnel simulation values. On the outer mold line (OML) downstream of the cavity, the flight values are about 20% larger than the wind tunnel simulation values. The higher heating bump factors observed in the flight simulations were due to the larger driving potential in terms of energy entering the cavity for the flight simulations. This is evidenced by the larger rate of increase in the total enthalpy through the boundary layer prior to the cavity for the flight simulation.
In support of the Columbia Accident Investigation, inviscid computations of the aerodynamic characteristics for various Shuttle Orbiter damage scenarios were performed using the FELISA unstructured CFD solver. Computed delta aerodynamics were compared with the reconstructed delta aerodynamics in order to postulate a progression of damage through the flight trajectory. By performing computations at hypervelocity flight and CF4 tunnel conditions, a bridge was provided between wind tunnel testing in Langley's 20-Inch CF4 facility and the flight environment experienced by Columbia during re-entry. The rapid modeling capability of the unstructured methodology allowed the computational effort to keep pace with the wind tunnel and, at times, guide the wind tunnel efforts. These computations provided a detailed view of the flowfield characteristics and the contribution of orbiter components (such as the vertical tail and wing) to aerodynamic forces and moments that were unavailable from wind tunnel testing. The damage scenarios are grouped into three categories. Initially, single and multiple missing full RCC panels were analyzed to determine the effect of damage location and magnitude on the aerodynamics. Next is a series of cases with progressive damage, increasing in severity, in the region of RCC panel 9. The final group is a set of wing leading edge and windward surface deformations that model possible structural deformation of the wing skin due to internal heating of the wing structure. By matching the aerodynamics from selected damage scenarios to the reconstructed flight aerodynamics, a progression of damage that is consistent with the flight data, debris forensics, and wind tunnel data is postulated.
This report documents the results of a computational study done to compute the inviscid longitudinal aerodynamic characteristics of the Space Shuttle Orbiter for Mach numbers 10 and 15 at angles of attack of 40, 50, 55, and 60 degrees. These computations were done to provide limited aerodynamic data in support of the Orbiter contingency abort task. The Orbiter had all the control surfaces in the undeflected position. The unstructured grid software FELISA was used for these computations with the equilibrium air option. Normal and axial force coefficients and pitching moment coefficients were computed. The hinge moment coefficients of the body flap and the inboard and outboard elevons were also computed. These results were compared with Orbiter Air Data Book (OADB) data and those computed using GASP. The comparison with the GASP results showed very good agreement in Cm and Ca at all the points. The computed axial force coefficients were smaller than those computed by GASP. There were noticeable differences between the present results and those in the OADB at angles of attack greater than 50 degrees.
NASA’s Advanced Composite Solar Sail System (ACS3) mission consists of a spacecraft that plans to be launched no earlier than April 2024. The spacecraft carries an 81 m2 solar sail that can produce effective ΔV to alter the initial 1000 km sun-synchronous orbit. The main objective of the mission is to demonstrate the capabilities of the solar sail to effectively change the semi-major axis of its initial orbit. Various composite materials were used to produce the sail, together with lightweight booms that will deploy from a 12U CubeSat bus. The ACS3 navigation team at NASA Ames Research Center has built a Flight Dynamics System (FDS) to provide mission navigation and to produce regular ephemeris once in orbit. The FDS can compute the orbit transfers that the spacecraft will perform once the sail is deployed. To achieve that, GPS data is obtained from the spacecraft telemetry and then is used with a Kalman filter and a smoother to obtain an orbit determination solution. The outcome of this process reduces the position and velocity uncertainty in a daily cadence. After that, the state vector output is used to feed a propagation model that includes the updated attitude and orbit of the spacecraft at that given moment. The trajectory model considers the updated attitude plan of the spacecraft as well as the environment conditions such as the solar weather to compute the associated drag and solar radiation pressure. This paper explains in detail the implementation of the FDS, as well as the solar sail solar radiation pressure trajectory model. We also present the results of several potential trajectory models under various assumptions of orbit parameters, attitude, environment, and material properties. In addition, we introduce a trajectory model for potential interplanetary use of an equivalent solar sail in the future. The outcome of this process will be critical to achieving the objective of determining effective semi-major axis change produced by the solar sail.
There is a lack of precise meteoroids orbit from video observations as most of the meteor stations use off-the-shelf CCD cameras. Few meteoroids orbit with precise semi-major axis are available using film photographic method. Precise orbits are necessary to compute the dust flux in the Earth s vicinity, and to estimate the ejection time of the meteoroids accurately by comparing them with the theoretical evolution model. We investigate the use of large CCD sensors to observe multi-station meteors and to compute precise orbit of these meteoroids. An ideal spatial and temporal resolution to get an accuracy to those similar of photographic plates are discussed. Various problems faced due to the use of large CCD, such as increasing the spatial and the temporal resolution at the same time and computational problems in finding the meteor position are illustrated.
A computation is undertaken of the rotation rate of a planet that has accreted from small solid planetesimals, as a function of the eccentricities of the planetesimals' orbits at or near the plane of the planet's orbit. The spin angular momentum generated by a given collision can be either positive or negative; due to the tendency to approximately cancel out, this perturbation process must be pursued to the second order in the ratio of the planetary radius to the planetesimal radial excursions' amplitude. It is found that material accreted from the area of the planet feeding zone's outer edges furnishes the planet positive spin angular momentum.
Timing data and real time orbit determination computer program
The relative importance of solar and trapped proton fluxes in the consideration of shielding requirements for geocentric space missions is analyzed. Using models of these particles, their fluences encountered by spacecraft in circular orbits are computed as functions of orbital altitude and inclination, mission duration, threshold energy (10 to 100 MeV), and risk factor (for solar protons only), and ratios of solar-to-trapped fluences are derived. It is shown that solar protons predominate for low-altitude polar and very high-altitude missions, while trapped protons predominate for missions at low and medium altitudes and low inclinations. It is recommended that if the ratio of solar-to-trapped protons falls between 0.1 and 10, both fluences should be considered in planning shielding systems.
The Preliminary Orbit Determination System (PODS) provides early orbit determination capability in the Trajectory Computation and Orbital Products System (TCOPS) for a Tracking and Data Relay Satellite System (TDRSS)-tracked spacecraft. PODS computes a set of orbit states from an a priori estimate and six tracking measurements, consisting of any combination of TDRSS range and Doppler tracking measurements. PODS uses the homotopy continuation method to solve a set of nonlinear equations, and it is particularly effective for the case when the a priori estimate is not well known. Since range and Doppler measurements produce multiple states in PODS, a screening technique selects the desired state. PODS is executed in the TCOPS environment and can directly access all operational data sets. At the completion of the preliminary orbit determination, the PODS-generated state, along with additional tracking measurements, can be directly input to the differential correction (DC) process to generate an improved state. To validate the computational and operational capabilities of PODS, tests were performed using simulated TDRSS tracking measurements for the Cosmic Background Explorer (COBE) satellite and using real TDRSS measurements for the Earth Radiation Budget Satellite (ERBS) and the Solar Mesosphere Explorer (SME) spacecraft. The effects of various measurement combinations, varying arc lengths, and levels of degradation of the a priori state vector on the PODS solutions were considered.
The nonlinear gyrokinetic theory in the tokamak configuration based on the two-step transform is developed; in the first step, we transform the magnetic potential perturbation to the Hamiltonian part, and in the second step, we transform away the gyroangle-dependent part of the perturbed Hamiltonian. Then the I-transform method is used to decoupled the perturbation part of the motion from the unperturbed motion. The application of the I-transform method to the computation of the guiding-center orbit and the guiding-center distribution function in tokamaks is presented. It is demonstrated that the I-transform method of the orbit computation which involves integrating only along the unperturbed orbit agrees with the conventional method which integrates along the full orbit. A numerical code based on the I-transform method is developed and two numerical examples are given to verify the new method.
The United States plans to establish a permanent manned presence at the Space Station Freedom in low earth orbit (LEO) and then carry out exploration of the solar system from this base. These plans may require orbital cryogenic propellant storage depots. The COOLANT program has been developed to analyze the thermodynamic performance of these depots to support design tradeoff studies. It was developed as part of the Long Term Cryogenic Storage Facility Systems Study for NASA/MSFC. This paper discusses the program structure and capabilities of the COOLANT program. In addition, the results of an analysis of a 200,000 lbm hydrogen/oxygen storage depot tankset using COOLANT are presented.
Orbital debris is a growing problem for the space industry and the world in general, and an important component of the problem is what happens when that debris reenters the Earth’s atmosphere. With more spacecraft opting for fiber-reinforced polymer (FRP) components, we need to understand the thermal destruction process of these materials during an atmospheric reentry and how much of the material can survive to impact the ground. The NASA Orbital Debris Program Office (ODPO) has developed a new charring model for FRP components to be integrated into the latest version of the Object Reentry Survival Analysis Tool (ORSAT). To validate this new model, the ODPO performed several test series using the Inductively Coupled Plasma (ICP) Torch facility at UT Austin. To measure the extent of charred material at different conditions, some of the test samples were scanned using x-ray computed tomography (CT) by the Astromaterials Curation Lab at Johnson Space Center. 3D image analysis was then used to calculate the volume and density of the char in each test sample. This paper presents the image analysis methodology, an assessment of the accuracy of the data analysis, and a comparison with calculations using the ORSAT charring model.
Computer program for transfer orbit trajectory optimization based on Lagrange multipliers, Newton-Raphson method, and power series
Computer program for calculating osculating values of Keplerian elements of satellite orbit
Orbit determination for Lunar Orbiter, describing computer program and studying Doppler residuals
Precise ephemerides have been determined for the U.S. Navy Geosat Exact Repeat Mission (ERM) using an improved gravity-field model, PTGF-4A (Shum et al. 1989). The Geosat orbits were computed in a terrestrial reference system which is tied to the reference system defined by satellite laser ranging (SLR) to Lageos through a survey between the Tranet Doppler receiver and the SLR system located at Wettzell, FRG. The remaining Doppler tracking station coordinates were estimated simultaneously with the geopotential in the PTGF-4A solution. In this analysis, three continuous 17-day Geosat orbits, which were computed using the 46-station Tranet data and global altimeter crossover data, have a crossover residual rms of 20 cm, indicating that the Geosat radial orbit error is of the order of 20 cm. The orbits computed based on data collected by a 7-station OPNET tracking network and crossover data have the same level of accuracy.