SEARCH · Search NASA
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.
An analytic solution for the J2 perturbed equatorial orbit is obtained in terms of elliptic functions and integrals. The necessary equations for computing the position and elocity vectors, and the time are given in terms of known functions. The perturbed periapsis and apoapsis distances are determined from the roots of a characteristic cubic.
A crawler system is designed to move a low-g/variable-g laboratory module along a tether between the Space Station and an attached space platform. An analysis is made of the effects of control law parameter change on the displacement, velocity, and acceleration of the crawler system. The control law is then modified by the addition of a constant-velocity section and the values of distance traveled, velocity, and acceleration are analyzed as a function of time. The power and torque equations are derived for a crawler system moving along a tether in orbit and numerical values of power and torque required for each prescribed movement are calculated versus time for four different cases using the control laws. A two-step control sequence is selected to permit initial location along the tether by distance traveled, followed by a vernier movement to reach the final desired constant net acceleration level. The components for the control system are identified and arranged in a block disgram configuration. The support subsystems are also identified. The sections were integrated to develop a procedure for the determination of crawler system performance requirements and the initial design of tether crawler systems.
In previous work the problem of injecting the Galileo and Ulysses spacecraft from low earth orbit into their respective interplanetary trajectories has been discussed for the single stage (Centaur) vehicle. The central issue, in the event of spherically distributed injection errors, is what happens to the vehicle? The difficulties addressed in this paper involve the multi-stage problem since both Galileo and Ulysses will be utilizing the two-stage IUS system. Ulysses will also include a third stage: the PAM-S. The solution is expressed in terms of probabilities for total percentage of escape, orbit decay and reentry trajectories. Analytic solutions are found for Hill's Equations of Relative Motion (more recently called Clohessy-Wiltshire Equations) for multi-stage injections. These solutions are interpreted geometrically on the injection sphere. The analytic-geometric models compare well with numerical solutions, provide insight into the behavior of trajectories mapped on the injection sphere and simplify the numerical two-dimensional search for trajectory families.
An existing anisotropic anodic plasma contactor model is extended to include neutral gas emissions. Within the framework of the contactor model, the inclusion of external ionization leads to an integrodifferential equation. By choosing the ratio of the contactor radius to the neutral gas-electron ionization mean free path at the contactor exit as a suitable expansion parameter, an analytical expression for the current-voltage characteristic can be found. The performance of this contactor model is examined via a dynamic simulation model developed for electrodynamic tethers. Detailed models of the geomagnetic field and the ionosphere are used to create a realistic environment. The highly nonlinear contactor voltagle-current characteristic is incorporated into a circuit equation, which includes radiation impedance, and is solved along the tether's orbit. Results show that it is more effective to fully ionize the contactor gas internally, than to partially ionize it externally, and that based on the specific power for the tether system, the optimum gas to use is argon. In addition, the effect of the radiation impedance on tether system performance is examined.
NASA is planning missions to the vicinity of the Sun-Earth L(sub 2) point, some involving a distributed system of telescope spacecraft, configured in a plane about a hub. Several sets of differential equations are written for the formation flight of such telescopes relative to the hub, with varying levels of fidelity. Effects are cast as additive perturbations to the circular restricted three-body problem, expanded in terms of the system distanced, to an accuracy of 10-20 m. These include Earth's orbital eccentricity, lunar motion, solar radiation pressure, and small thrusting forces. Simulations validating the expanded differential equations are presented.
A document presents a concept for relatively inexpensive delivery of propellant to a large fuel depot in low orbit around the Earth, for use in rockets destined for higher orbits, the Moon, and for remote planets. The propellant is expected to be at least 85 percent of the mass needed in low Earth orbit to support the NASA Exploration Vision. The concept calls for the use of many small ( 10 ton) spin-stabilized, multistage, solid-fuel rockets to each deliver 250 kg of propellant. Each rocket would be winched up to a balloon tethered above most of the atmospheric mass (optimal altitude 26 2 km). There, the rocket would be aimed slightly above the horizon, spun, dropped, and fired at a time chosen so that the rocket would arrive in orbit near the depot. Small thrusters on the payload (powered, for example, by boil-off gases from cryogenic propellants that make up the payload) would precess the spinning rocket, using data from a low-cost inertial sensor to correct for small aerodynamic and solid rocket nozzle misalignment torques on the spinning rocket; would manage the angle of attack and the final orbit insertion burn; and would be fired on command from the depot in response to observations of the trajectory of the payload so as to make small corrections to bring the payload into a rendezvous orbit and despin it for capture by the depot. The system is low-cost because the small rockets can be mass-produced using the same techniques as those to produce automobiles and low-cost munitions, and one or more can be launched from a U.S. territory on the equator (Baker or Jarvis Islands in the mid-Pacific) to the fuel depot on each orbit (every 90 minutes, e.g., any multiple of 6,000 per year).
The flight environment of a reentry vehicle is predicted from the numerical solution of fluid dynamics equations for the region between the blunt nose and the swept wings. The inviscid portion of the shock layer is modeled by the Euler equations, but the laminar viscous flow adjacent to the wall is modeled by the approximate parabolic Navier-Stokes equations. The approximations made to the axial gradients of pressure and diffusive fluxes enable the coupled inviscid and viscous equations to be solved efficiently along the body axis. The equilibrium air aftbody code contains significant improvements over its predecessor which only considers a perfect gas model and noncircular configurations. The inclusion of numerical damping either explicitly or implicitly has extended its capabilities for predicting flow field around a winged configuration at higher Machs and greater angles of attack. The results are obtained on the cylindrical coordinates and satisfactory for the Shuttle Orbiter at a free-stream Mach number of 22 and an angle of attack of 40 deg. Also discussed are the inviscid formulation and its application for wind-tunnel conditions.
The long-term evolution characteristics (and stability) of an orbit are best characterized using a mean element propagation of the perturbed two body variational equations of motion. The averaging process eliminates short period terms leaving only secular and long period effects. In this study, a non-traditional approach is taken that averages the variational equations using adaptive numerical techniques and then numerically integrating the resulting EOMs. Doing this avoids the Fourier series expansions and truncations required by the traditional analytic methods. The resultant numerical techniques can be easily adapted to propagations at most solar system bodies.
In this paper, we examine the effects of radiation pressure, Poynting-Robertson (PR) drag, and solar wind drag on dust grains trapped in mean motion resonances with the Sun and Jupiter in the restricted (negligible dust mass) three-body Problem. We especially examine the evolution of dust grains in the 1:1 resonance. As a first step, the Sun and Jupiter are idealized to both be in circular orbit about a common center of mass (circular restricted three-body problem). From the equation of motion of the dust particle in the rotating reference frame, the drag-induced time rate of change of its Jacobi "constant," C, is then derived and expressed in spherical coordinates. This new mathematical expression in spherical coordinates shows that C, in the 1:1 resonance, both oscillates and secularly increases with increasing time. The new expression gives rise to an easy understanding of how an orbit evolves when the radiation force and solar wind drag are included. All dust grain orbits are unstable in time when PR and solar wind drag are included in the Sun-Jupiter-dust system. Tadpole orbits evolve into horseshoe orbits; and these orbits continuously expand in size to lead to close encounters with Jupiter. Permanent trapping is impossible. Orbital evolutions of a dust grain trapped in the 1:1 resonance in the planar circular, an inclined case, an eccentric case, and the actual Sun-Jupiter case are numerically simulated and compared with each other and show grossly similar time behavior. Resonances other than 1:1 are also explored with the new expression. Stable exterior resonance trapping may be possible under certain conditions. One necessary condition for such a trap is derived. Trapping in interior resonances is shown to be always unstable.
A technique is described for improving the positioning accuracy of earth survey satellites. Tracking data from Landsat-1, taken in contiguous spans equal in length to the repeat period of 18 days, have been reduced in a least-squares orbit-determination scheme to estimate two pairs of lumped resonance coefficients. Seven separate data spans yielded resonance coefficient pairs consistent to better than 10% of their values, based on the GEM-1 gravity model. Constraint equations were developed. The resonance terms were embedded in the GEM-8 gravity field to yield further improvements in orbital fit to the tracking data, and tests suggest a 20 m position capability for Landsat orbits when the resonance model is adopted. The lumped terms may be used to improve the perturbation model for all Landsat orbits, and the approach is applicable to the orbit of any satellite whose repeat period is long enough for geopotential resonance to have a significant effect on its motion.
This study is an extension of previous work to evaluate the applicability of Design of Experiments (DOE)/Response Surface Methodology to on-orbit thermal analysis. The goal was to determine if the methodology could produce a Response Surface Equation (RSE) that predicted the thermal model temperature results within +/-10 F. An RSE is a polynomial expression that can then be used to predict temperatures for a defined range of factor combinations. Based on suggestions received from the previous work, this study used a model with simpler geometry, considered polynomials up to fifth order, and evaluated orbital temperature variations to establish a minimum and maximum temperature for each component. A simplified Outer Mold Line (OML) thermal model of the Orion spacecraft was used in this study. The factors chosen were the vehicle's Yaw, Pitch, and Roll (defining the on-orbit attitude), the Beta angle (restricted to positive beta angles from 0 to 75), and the environmental constants (varying from cold to hot). All factors were normalized from their native ranges to a non-dimensional range from -1.0 to 1.0. Twenty-three components from the OML were chosen and the minimum and maximum orbital temperatures were calculated for each to produce forty-six responses for the DOE model. A customized DOE case matrix of 145 analysis cases was developed which used analysis points at the factor corners, mid-points, and center. From this data set, RSE s were developed which consisted of cubic, quartic, and fifth order polynomials. The results presented are for the fifth order RSE. The RSE results were then evaluated for agreement with the analytical model predictions to produce a +/-3(sigma) error band. Forty of the 46 responses had a +/-3(sigma) value of 10 F or less. Encouraged by this initial success, two additional sets of verification cases were selected. One contained 20 cases, the other 50 cases. These cases were evaluated both with the fifth order RSE and with the analytical model. For the maximum temperature predictions, 12 of the 23 components had all predictions within +/-10 F and 17 were within +/-20 F. For the minimum temperature predictions, only 4 of the 23 components (the four radiator temperatures), were within the 10 F goal. The maximum temperature RSEs were then run through 59,049 screening cases. The RSE predictions were then filtered to find 55 cases that produced the hottest temperatures. These 55 cases were then analyzed using the thermal model and the results compared against the RSE predictions. As noted earlier, 12 of the 23 responses were within +/-10 F at 17 within +/-20 F. These results demonstrate that if properly formulated, an RSE can provide a reliable, fast temperature prediction. Despite this progress, additional work is needed to determine why the minimum temperatures responses and 6 of the hot temperature responses did not produce reliable RSEs. Recommend focus areas are the model itself (arithmetic vs. diffusion nodes) and seeking consultations with statistical application experts.
A method is presented in satellite altimetry that attempts to simultaneously determine the geoid and sea surface toography with minimum wavelengths of about 500 km and to reduce the radial orbit errors caused by geopotential uncertainties. The modeling of these errors is made using the linearized Lagrangian perturbation theory. Observation equations are developed using sea surface heights and crossover discrepancies as observables. A minimum variance solution with prior information can then provide estimates of parametrs representing the sea surface topography and corrections to the orbit. The potential of the method is demonstrated in a solution where simulated geopotential errors and the Levitus sea surface topography are used to generate the observables for a Seasat 3 day arc. The simulation results suggest that the method can be used to efficiently process real altimeter data.
The precession and nutation of the Earth's equator arise from solar, lunar, and planetary torques on the oblate Earth. The mean lunar orbit plane is nearly coincident with the ecliptic plane. A small tilt out of the ecliptic is caused by planetary perturbations and the Earth's gravitational harmonic J(sub2). These planetary perturbations on the lunar orbit result in torques on the oblate Earth which contribute to precession, obliquity rate, and nutation while the J(sub 2) perturbations contribute to precession and nutation. Small additional contributions to the secular rates arise from tidal effects and planetary torques on the Earth's bulge. The total correction to the obliquity rate is -0.024sec/century, it is an observable motion in space (the much larger conventional obliquity rate is wholly from the motion of the ecliptic, not the equator), and it is not present in the IAU-adopted expressions for the orientation of the Earth's equator. The J(sub2) effects have generally been allowed for in past nutation theories and some procession theories. For the planetary effect, the contributions to the 18.6 yr nutation are -0.03 mas (milliarcseconds) for the in-phase Delta phi plus out-of-phase contributions of 0.14 mas in Delta phi and -0.03 mas in Delta epsilon. The latter terms demonstrate that out-of-phase contributions can arise by means other than dissipation. The sum of the contributions to the precession rate is considered and the inferred value of the moment of inertia combination (C-A)/C, which is used to scale the coefficients in the nutation series, is evaluated. Using an updated value for the precession rate, the rigid body (C-A)/C = 0.003 273 763 4 which, in combination with a satellite-derived J(sub2), gives a normalized polar moment of inertia C/MR(exp2) = 0.330 700 7. The planetary contributions to the precession and obliquity rates are not constant for long times causing accelerations in both quantities. Acceleration in precession also arises from tides and changing J(sub2). Contributions from the improved theory, masses, ecliptic motion, and measured values of the precession rate and obliquity are combined to give expressions (polynomials in time) for precession, obliquity, and Greenwich Mean Sidereal Time.
The precession and nutation of the Earth's equator arise from solar, lunar, and planetary torques on the oblate Earth. The mean lunar orbit plane is nearly coincident with the ecliptic plane. A small tilt out of the ecliptic is caused by planetary perturbations and the Earth's gravitational harmonic J(sub 2). These planetary perturbations on the lunar orbit result in torques on the oblate Earth which contribute to precession, obliquity rate, and nutation while the J(sub 2) perturbations contribute to precession and nutation. Small additional contributions to the secular rates arise from tidal effects and planetary torques on the Earth's bulge. The total correction to the obliquity rate is -0.024 sec/century, it is an observable motion in space (the much larger conventional obliquity rate is wholly from the motion of the ecliptic, not the equator), and it is not present in the IAU-adopted expressions for the orientation of the Earth's equator. The effects have generally been allowed for in past nutation theories and some precession theories. For the planetary effect, the contributions to the 18.6 yr nutation are -0.03 mas (milliarcseconds) for the in-phase Delta(psi) plus out-of-phase contributions of 0.14 mas in Delta(psi) and -0.03 mas in Delta(sub epsilon). The latter terms demonstrate that out-of-phase contributions can arise by means other than dissipation. The sum of the contributions to the precession rate is considered and the inferred value of the moment of inertia combination (C-A)/C, which is used to scale the coefficients in the nutation series, is evaluated. Using an updated value for the precession rate, the rigid body (C-A)/C =0.003 273 763 4 which, in combination with a satellite-derived J(sub 2), gives a normalized polar moment of inertia C/MR(exp 2) = 0.330 700 7. The planetary contributions to the precession and obliquity rates are not constant for long times causing accelerations in both quantities. Acceleration in precession also arises from tides and changing J(sub 2) Contributions from the improved theory, masses, ecliptic motion, and measured values of the precession rate and obliquity are combined to give expressions (polynomials in time) for precession, obliquity, and Greenwich Mean Sidereal Time.
We present an improved lunar digital elevation model (DEM) covering latitudes within +/-60 deg, at a horizontal resolution of 512 pixels per degree ( approx.60 m at the equator) and a typical vertical accuracy approx.3 to 4 m. This DEM is constructed from approx.4.5 ×10(exp 9) geodetically-accurate topographic heights from the Lunar Orbiter Laser Altimeter (LOLA) onboard the Lunar Reconnaissance Orbiter, to which we co-registered 43,200 stereo-derived DEMs (each 1 deg×1 deg) from the SELENE Terrain Camera (TC) ( approx.10(exp 10) pixels total). After co-registration, approximately 90% of the TC DEMs show root-mean-square vertical residuals with the LOLA data of < 5 m compared to approx.50% prior to co-registration. We use the co-registered TC data to estimate and correct orbital and pointing geolocation errors from the LOLA altimetric profiles (typically amounting to < 10 m horizontally and < 1 m vertically). By combining both co-registered datasets, we obtain a near-global DEM with high geodetic accuracy, and without the need for surface interpolation. We evaluate the resulting LOLA + TC merged DEM (designated as "SLDEM2015") with particular attention to quantifying seams and crossover errors.
“Bumper” is a computer program for analyzing spacecraft micrometeoroid and orbital debris (MMOD) risk. Bumper was developed in the late-1980s and has been continuously maintained and used since. The user base has grown from a few government entities now include numerous commercial entities as well. The National Aeronautics and Space Administration (NASA) Johnson Space Center (JSC) Hypervelocity Impact Technology (HVIT) group is responsible for all aspects of the Bumper software. Bumper has been used to characterize MMOD risk on many spacecraft. All of the International Space Station (ISS) modules, visiting vehicles and numerous external components and systems have been analyzed. Bumper was used to analyze the Space Shuttle, Orion, and many space probes, telescopes and satellites. Bumper is also being used to analyze future spacecraft such as the Deep Space Gateway (DSG) and Mars Sample Return (MSR) missions. The Bumper Configuration Control Board (CCB) ensures that all changes to the code are approved, reviewed, and documented. The current Bumper version – “Bumper 3” – is a Fortran executable that utilizes a 64-bit architecture. Bumper has numerous features that make it a powerful tool for analyzing spacecraft MMOD risk. Bumper uses the latest orbital debris and meteoroidenvironment models. Bumper also has a large library of ballistic limit “damage” equations available that can be used for a wide variety of MMOD shielding configurations. Bumper can also handle large spacecraft finite elementmodels (FEMs) and conducts checks of the model. This paper introduces Bumper and the MMOD risk analysis process using a simplified cube-shaped spacecraft model
Future space structures will have a low mass density and high flexibility, with ground test dynamic behavior differing significantly from that in zero-G orbit. Attention is presently given to the vibration behavior of a beam deformed by its own weight; the results obtained by the differential equations for both the static and dynamic responses of a large, simply supported beam, which are derived and solved analytically, allow ground test experiment measurements to be used for orbital dynamic characteristics verification efforts.
satellite encounters. Over this period, there have been several papers describing the orbit determination process and performance up through 2016 [1-5]. In April of 2017, Cassini will enter its Grand Finale mission phase when it will traverse the gap between the D-ring and the Saturn atmosphere twenty-two times before plunging deep into the atmosphere to end the mission. The lack of targeted satellite encounters during this period necessitates updates to the nominal Cassini Orbit Determination (OD) process. This paper describes these planned adaptations for the operation of the Grand Finale. During the Equinox and Solstice Mission Phase (2008-2016), navigation analysis has been divided into segments focused on two particular targeted satellite encounters, called an “arc”. Maneuvers in an arc were usually targeted to encounter B-plane position and time, so the OD state and covariance were mapped forward to the B-plane of the encounter within the arc. Trajectory dispersions during the Grand Finale need instead to be mapped to equator crossings and targeted Cartesian positions. In addition, trajectory arcs have typically covered a few orbital revolutions (~2-8 weeks), in order to span the time between two encounters. However, the Grand Finale will encompass five months of time without an encounter which necessitates an adjusted arc strategy. A modified arc strategy was developed based on OD behavior during long multi-rev periods between encounters in the year leading up to the Grand Finale. The OD covariance study conducted for the Grand Finale mission phase will also be examined.