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.
Mathematical models and orbital mechanics for space shuttle guidance, control, and navigation during deorbit and atmospheric operations
Existing analysis of Orbiter-RMS-Payload kinematics were surveyed, including equations dealing with the two body kinematics in the presence of a massless RMS and compares analytical explicit solutions with numerical solutions. For the following operational phases of the RMS numerical demonstration, problems are provided: (1) payload capture; (2) payload stowage and removal from cargo bay; and (3) payload deployment. The equation of motion provided accounted for RMS control forces and torque moments and could be extended to RMS flexibility and control loop simulation without increasing the degrees of freedom of the two body system.
A large, free-flying space observatory is planned for (1) locating and characterizing gamma ray sources, in particular peculiar objects such as pulsars; (2) studying diffuse gamma ray emission, both from our galaxy and from extragalactic space, (3) investigating other galaxies, especially radio galaxies, Seyfert galaxies and QSOs; and (4) understanding the origin of gamma ray bursts. The spacecraft will weight about 10,000 kg, with about 6,500 kg available for scientific instruments, providing about 400 watts of power to the instruments, and about 12,000 bits per second of telemetry. The spacecraft which will be launched by the space shuttle into a 400 to 500 km circular orbit, inclined at 28.5 deg to the equator, will be three-axis stabilized, with orientation capability of about 20 arc minutes and aspect determination of 2 arc minutes. Absolute timing of 0.1 millisecond will be provided by the spacecraft. The instruments will be co-aligned along one axis of the spacecraft to optimize the spectral coverage of a particular region. A two-year mission is planned with launch in 1984. Possible retrieval by the shuttle at the end of the mission is considered.
A study is made of the collisionless tearing-mode stability properties of a field-reversed plasma layer whose temperature distribution is anisotropic. The plasma is confined by its self magnetic field with no external field. A kinetic description is used for both ions and electrons. The effects of the axis-crossing and nonaxis-crossing orbits are discussed. It is found that the conventional technique of matching the inner and outer asymptotic solutions at the electron inner-region is inadequate for the anisotropic case. An intermediate region in which the axis-crossing ion orbits are important is identified. The eigenvalue equation is solved using both analytic approximations and numerical methods to obtain the eigenmode structure and the linear dispersion relation. Previously announced in STAR as N84-14920
Ten years ago computer simulations of the flow field around the Space Shuttle Orbiter were limited to inviscid calculations for the windward side of the forebody and viscous calculations for selected 2-D problems. Advances in computer hardware and numerical methods during the past ten years made it possible to calculate viscous flow over the complete orbiter configuration at angle of attack. The equations solved are the Reynolds-averaged, Navier-Stokes equation, simplified by either the thin-layer or parabolized approximation. An algebraic eddy viscosity model is used for turbulent flow. The free stream is assumed to be a perfect gas for wind tunnel conditions and a real gas in thermodynamic equilibrium for flight conditions. Four examples of recent computer simulations are presented. Flow field results include oil flow patterns on the surface and Mach number contours, isobars, and cross-flow velocity vectors in the shock layer.
Ten years ago computer simulations of the flow field around the Space Shuttle Orbiter were limited to inviscid calculations for the windward side of the forebody and viscous calculations for selected 2-D problems. Advances in computer hardware and numerical methods during the past ten years made it possible to calculate viscous flow over the complete orbiter configuration at angle of attack. The equations solved are the Reynolds-averaged, Navier-Stokes equation, simplified by either the thin-layer or parabolized approximation. An algebraic eddy viscosity model is used for turbulent flow. The free stream is assumed to be a perfect gas for wind tunnel conditions and a real gas in thermodynamic equilibrium for flight conditions. Four examples of recent computer simulations are presented. Flow field results include oil flow patterns on the surface and Mach number contours, isobars, and cross-flow velocity vectors in the shock layer.
BUMPERII is modular computer program implementing numerical solution technique to calculate probability of no penetration (PNP) of spacecraft by orbiting artificial debris or impacts of meteoroids. Equation used to calculate PNP based on Poisson-distribution model applicable to similar analysis of smaller craft, but reflects more rigorous mathematical modeling of size, shape, orientation, and other characteristics of spacecraft that affect impacts, and necessary for PNP analysis of such larger structures as components of space station. Written in FORTRAN 77.
In International Terrestrial Reference Frame (ITRF) 2005 and ITRF2008, the approach for the construction of solutions by the IERS has been for individual analysis centers of each technique to process geodetic tracking data, and for each technique to develop a solution (or contribution) that is integrated into the final ITRF solution by careful combination of the technique solutions. The connections between the geodetic networks are realized by the application of local ties. In an alternate approach, we may assure processing homogeneity by creating normal equations for different techniques with the same orbit determination software, using identically derived algorithms. Another derivative of this approach is to realize the ties between the techniques using satellites tracked with multiple techniques; in effect tieing the networks together using satellite dynamics. In this solution, we develop a time series and a set of cumulative solutions from Satellite Laser Ranging (SLR) & Doppler Orbitography and Radio-positioning Integrated by Satellite (DORIS) based on homogeneous processing with the NASA GEODYN precise orbit determination suite of programs, where we jointly combine weekly the SLR data to Lageos1, Lageos2, Starlette, and Stella with the DORIS data from SPOT2-SPOT5, as well as satellites that utilize both techniques (TOPEX/Poseidon, Envisat, Jason-2). We discuss the modeling that is applied including upgrades implemented since the submission of the GSC ITRF2008 contributions for IDS. Firstly, we compare the SLR-only solutions comprising four geodetic satellites with the standard approach of utilizing only Lageos1 & Lageos2. Secondly, we evaluate the impact on the DORIS coordinates of the joint analysis with the SLR data.
Traditional azimuth-elevation ("az-el") stages are used to point a variety of devices ranging from large optical telescopes to tank guns. Such a stage typically has an elevation stage having a horizontal axis mounted on an azimuth stage with a vertical axis. Both stages are often motorized. Optical communications transceivers often require two-axis motorized control, as when the communications link is between a ground station and an aircraft or satellite. In such applications, the traditional azimuth-elevation stage has two important drawbacks: a gimbal lock exclusion zone at zenith and susceptibility to pointing errors caused by backlash. Az-el stages in which the azimuth stage cannot rotate a full 360deg have the additional drawback of an azimuth exclusion zone. The diagonal-axes stage described here mitigates or eliminates all of these problems. Instead of one vertical axis and one horizontal axis, a diagonal-axes stage has two horizontal axes, both oriented at 45 to the trajectory of the target. For example, a ground station located on the equator tracking a satellite with an equatorial orbit would have one axis parallel to northeast and southwest, and the other axis parallel to northwest and southeast. The diagonal-axes stage is considerably less vulnerable to backlash. If it is correctly oriented, its axes rotate in only one direction during an overhead pass by a satellite. As a result, the effects of backlash may be inherently eliminated. If the gravity-induced torque on either axis changes during the pass, then backlash may become important during the part of the pass where the gravity torque, instead of opposing the motion of the stage, pushes the stage in the direction of motion. This can result in the loss of gear-to-gear contact in one or more stages of the gear reduction mechanism. In this case, a preload spring used to eliminate backlash need only be sufficiently strong to overcome the gravity torque, i.e. it need not overcome friction in the gear train. The diagonal-axes stage is not backlash-free for arbitrary target trajectories such as an aircraft might execute. If properly oriented for any particular satellite, however, it is backlash-free for all passes of that satellite, which will trace out parallel paths on the sky, and for all passes of any other satellite that are perpendicular to the first. It will also be backlash-free for some fraction of other satellite trajectories.
Sulfur dioxide (SO2) is a major air pollutant that contributes to acid rain and aerosol formation (e.g., sulfates), adversely affects the environment and human health, and explosive volcanic SO2 emissions can impact climate. The majority of SO2 emissions are related to anthropogenic processes (e.g., fossil fuel burning, metal ore smelting operations), although natural processes such as volcanic eruptions and degassing also play an important role as anthropogenic SO2 emissions continue to decline. Generally, the most interest in volcanoes occurs during major eruptions. We will focus on comparing volcanic SO2 outgassing that occurs on an almost daily basis from lesser known volcanoes using satellite data. At NASA’s Global Sulfur Dioxide Monitoring Home page (https://so2.gsfc.nasa.gov/), we have been posting daily SO2 maps from 40 volcanic and industrial regions around the world using measurements from three satellite instruments; the Ozone Monitoring Instrument (OMI) onboard NASA’s Earth Observing System Aura satellite, the Ozone Monitoring and Profiler Suite (OMPS) onboard the NASA-NOAA Suomi National Polar-orbiting Partnership (NPP) satellite, and the TROPOspheric Monitoring Instrument (TROPOMI) onboard the ESA/Copernicus Sentinel-5 Precursor satellite. These instruments in low Earth sun-synchronous polar orbits with 1:30-2:00 pm equator crossing local time provide daily SO2 maps at different spatial resolutions: 13 x 24 km2, 50 x 50 km2 and 5.5 x 3.5 km2 for OMI, OMPS and TROPOMI respectively. Data from OMI are available since October 2004 (partial coverage since 2008), from OMPS since 2012 and from TROPOMI since 2018. We will present comparative SO2 mass time-series (see Hunga-Tonga plot) and statistical analyses of recent eruptions that have data from all the instruments.
NASA uses six-degrees-of-freedom (6-DOF) simulations tools to design, test, develop Guidance Navigation and Control (GN&C) software, and certify vehicle performance prior to flight. Therefore, it is critical that the 6-DOF tools used for vehicle design and certification are validated. The focus of this work is the vali-dation of the NASA Marshall Space Flight Center 6-DOF “GeneraLized Aero-space Simulation in Simulink” (GLASS) framework tool. The GLASS tool framework is currently used to support NASA GN&C insight for the Human Landing System (HLS) project, simulating vehicle dynamics during lunar descent and as-cent. The GLASS framework utilizes the off-the-shelf Mathworks ® Simscape Multibody® toolbox to model vehicle multi-body dynamics. NASA’s Engineering and Safety Center (NESC) provides a set of 6-DOF simulation verification “check cases” that are available to any user needing to verify 6-DOF tools. The check cases contain seventeen atmospheric and twenty-six orbital test scenarios are provided to validate equations of motion, environmental models (e.g., atmosphere, gravitation, and geodesy) and tool propagators. This paper compares GLASS 6-DOF simulation results against the NESC check cases’ results via simulation-to-simulation comparisons. The comparisons demonstrate that GLASS simulation results are “in family” with the outputs of the applicable NASA NESC check cases and verifies the GLASS core framework dynamics and the correct implementation of the check case scenario models.
NASA uses six-degrees-of-freedom (6-DOF) simulations tools to design, test, develop Guidance Navigation and Control (GN&C) software, and certify vehicle performance prior to flight. Therefore, it is critical that the 6-DOF tools used for vehicle design and certification are validated. The focus of this work is the validation of the NASA Marshall Space Flight Center 6-DOF “GeneraLized Aerospace Simulation in Simulink” (GLASS) framework tool. The GLASS tool framework is currently used to support NASA GN&C insight for the Human Landing System (HLS) project, simulating vehicle dynamics during lunar descent and ascent. The GLASS framework utilizes the off-the-shelf Mathworks (R) Simscape (TM) Multibody (TM) toolbox to model vehicle multi-body dynamics. NASA’s Engineering and Safety Center (NESC) provides a set of 6-DOF simulation verification “check cases” that are available to any user needing to verify 6-DOF tools. The check cases contain seventeen atmospheric and twenty-six orbital test scenarios are provided to validate equations of motion, environmental models (e.g., atmosphere, gravitation, and geodesy) and tool propagators. This paper compares GLASS 6-DOF simulation results against the NESC check cases’ results via simulation-to-simulation comparisons. The comparison demonstrates that GLASS simulation results are “in family” with the outputs of the applicable NASA NESC check-cases and verify the GLASS core framework dynamics and the implementation of the check case scenario models.
The motion of a satellite in orbit, subject to atmospheric force and the motion of a reentry vehicle are governed by gravitational and aerodynamic forces. This suggests the derivation of a uniform set of equations applicable to both cases. For the case of satellite motion, by a proper transformation and by the method of averaging, a technique appropriate for long duration flight, the classical nonlinear differential equation describing the contraction of the major axis is derived. A rigorous analytic solution is used to integrate this equation with a high degree of accuracy, using Poincare's method of small parameters and Lagrange's expansion to explicitly express the major axis as a function of the eccentricity. The solution is uniformly valid for moderate and small eccentricities. For highly eccentric orbits, the asymptotic equation is derived directly from the general equation. Numerical solutions were generated to display the accuracy of the analytic theory.
The back-scattering of outgas contamination near an orbiting spacecraft due to intermolecular collisions was analyzed. Analytical tools were developed for making reasonably accurate quantitative estimates of the outgas contamination return flux, given a knowledge of the pertinent spacecraft and orbit conditions. Two basic collision mechanisms were considered: (1) collisions involving only outgas molecules (self-scattering) and (2) collisions between outgas molecules and molecules in the ambient atmosphere (ambient-scattering). For simplicity, the geometry was idealized to a uniformly outgassing sphere and to a disk oriented normal to the freestream. The method of solution involved an integration of an approximation of the Boltzmann kinetic equation known as the BGK (or Krook) model equation. Results were obtained in the form of simple equations relating outgas return flux to spacecraft and orbit parameters. Results were compared with previous analyses based on more simplistic models of the collision processes.
In order to study and control the attitude of a spacecraft, it is necessary to understand the natural motion of a body in orbit. Assuming a spacecraft to be a rigid body, dynamics describes the complete motion of the vehicle by the translational and rotational motion of the body. The Simulink Attitude Analysis Model applies the equations of rigid body motion to the study of a spacecraft?s attitude in orbit. Using a TCP/IP connection, Matlab reads the values of the Remote Manipulator System (RMS) hand controllers and passes them to Simulink as specified torque and impulse profiles. Simulink then uses the governing kinematic and dynamic equations of a rigid body in low earth orbit (LE0) to plot the attitude response of a spacecraft for five seconds given known applied torques and impulses, and constant principal moments of inertia.
Distributed parameter models of the Solar Array Flight Experiment, the Mini-MAST truss, and Space Station Freedom assembly are discussed. The distributed parameter approach takes advantage of (1) the relatively small number of model parameters associated with partial differential equation models of structural dynamics, (2) maximum-likelihood estimation using both prelaunch and on-orbit test data, (3) the inclusion of control system dynamics in the same equations, and (4) the incremental growth of the structural configurations. Maximum-likelihood parameter estimates for distributed parameter models were based on static compliance test results and frequency response measurements. Because the Space Station Freedom does not yet exist, the NASA Mini-MAST truss was used to test the procedure of modeling and parameter estimation. The resulting distributed parameter model of the Mini-MAST truss successfully demonstrated the approach taken. The computer program PDEMOD enables any configuration that can be represented by a network of flexible beam elements and rigid bodies to be remodeled.
Within the context of the restricted problem of three bodies the effects caused by varying the mass ratio of the primaries and the eccentricity of their orbits, upon periodic orbits of the infinitesimal mass which are numerical continuations of circular orbits in the ordinary problem of two bodies are shown. A recursive power series technique is used to numerically integrate the equations of motion as well as the first variational equations in order to generate a two parameter family of perodic orbits and identify the linear stability characteristics. Seven such families are investigated with equally spaced mass ratios from 0.0 to 1.0 and eccentricities of the orbits of the primaries in a range 0.0 to 0.6. Stable orbits are associated with large distances of the infinitesimal mass from the perturbing primary, nearly circular motion of the primaries, and small mass ratios of the primaries. Unstable orbits for the infinitesimal mass are associated with small distances from the perturbing primary, highly elliptic orbits of the primaries and large mass ratios.
A mathematical model is developed to predict the dynamics of the proposed orbiting Spacecraft Control Laboratory Experiment (SCOLE) during the station keeping phase. The equations of motion are derived using a Newton-Euler formulation. The model includes the effects of gravity, flexibility, and orbital dynamics. The control is assumed to be provided to the system through the Shuttle's three torquers, and through six actuators located by pairs at two points on the mast and at the mass center of the reflector. The modal shape functions are derived using the fourth order beam equation. The generic mode equations are derived to account for the effects of the control forces on the modal shape and frequencies. The equations are linearized about a nominal equilibrium position. The linear regulator theory is used to derive control laws for both the linear model of the rigidized SCOLE as well as that of the actual SCOLE including the first four flexible modes. The control strategy previously derived for the linear model of the rigidized SCOLE is applied to the nonlinear model of the same configuration of the system and preliminary single axis slewing maneuvers conducted. The results obtained confirm the applicability of the intuitive and appealing two-stage control strategy which would slew the SCOLE system, as if rigid to its desired position and then concentrate on damping out the residual flexible motions.