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.
FORTRAN 4 program for use on IBM 7094 computer for calculating temperature field of shaft seals is described. Shaft seals are composed of basically axisymmetric bodies. Cylindrical coordinate system is used as basis for analysis. Program is quite general and can be applied to variety of axisymmetric body problems.
Review of some of the considerations pertinent to efficient orbit utilization in the design of domestic communications satellite systems. A strategy is developed to efficiently locate a heterogeneous system of satellites within the available arc and provide room for future growth. A practical design is illustrated, using a computer simulation model, for the placement of 25 satellites within 73% of the available arc employing frequency and polarization coordination techniques. A number of widely variable factors that influence satellite spacing are examined. These factors include such critical system elements as telephony and television interference noise limits, frequency plan coordination, polarization plan coordination, ground antenna diameter, signal protection ratio, and satellite station keeping.
Skewed alignment of two redundant conventional inertial measuring units permits nonambiguous detection and isolation of hard and soft failures in real time by an airborne computer. Accelerometer outputs and gimbal readouts are monitored periodically, and attitude rate and velocity error vectors are computed from these data. Magnitudes of these vectors provide failure detection, and projection of these error vectors onto the coordinate axes of the two clusters permits isolation. A detailed Monte Carlo simulation of one version of the mechanization as applied to Space Shuttle boost trajectories demonstrates effectiveness down to very low levels of inertial instrument performance failures. The results indicate that worst case overall navigation performance occurs when accelerometer failures are of the order of 20 sigma and gyro failures are about 100 sigma for conventional state-of-the-art IMU instruments.
An iterative numerical method has been developed for the calculation of steady, three-dimensional, viscous, compressible flow fields in centrifugal compressor impellers. The computer code, which embodies the method, solves the steady three dimensional, compressible Navier-Stokes equations in rotating, curvilinear coordinates. The solution takes place on blade-to-blade surfaces of revolution which move from the hub to the shroud during each iteration.
Equations have been derived that transform perspectively viewed planar surface coordinates, as seen in a photograph, into coordinates of the original plane surface. These transformation equations are developed in terms of nine geometric variables that define the photographic setup and are redefined in terms of eight parameters. The parameters are then treated as independent quantities that fully characterize the transformation and are expressed directly in terms of the four corner coordinates of a reference rectangle in the object plane and their coordinates as seen in a photograph. Vehicle position is determined by transforming the perspectively viewed coordinate position of a representative vehicle target into runway coordinates. Vehicle heading is determined from the runway coordinates of two vehicle target points. When the targets are elevated above the plane of the reference grid, the computation of the heading angle is unaffected; however, the computation of the target position may require adjustment of two parameters. Methods are given for adjusting the parameters for elevation and an example is included for both nonelevated and elevated target conditions.
Calculated dynamic blade loads were compared with measured loads over a range of yaw stiffnesses of the DOE/NASA Mod-O wind turbine to evaluate the performance of two versions of the MOSTAS computer code. The first version uses a time-averaged coefficient approximation in conjunction with a multi-blade coordinate transformation for two bladed rotors to solve the equations of motion by standard eigenanalysis. The second version accounts for periodic coefficients while solving the equations by a time history integration. A hypothetical three-degree of freedom dynamic model was investigated. The exact equations of motion of this model were solved using the Floquet-Lipunov method. The equations with time-averaged coefficients were solved by standard eigenanalysis.
A computer code is presented for predicting incident thermal radiation from defined plume gas properties in either axisymmetric or cylindrical coordinate systems. The radiation model is a statistical band model for exponential line strength distribution with Lorentz/Doppler line shapes for 5 gaseous species (H2O, CO2, CO, HCl and HF) and an appoximate (non-scattering) treatment of carbon particles. The Curtis-Godson approximation is used for inhomogeneous gases, but a subroutine is available for using Young's intuitive derivative method for H2O with Lorentz line shape and exponentially-tailed-inverse line strength distribution. The geometry model provides integration over a hemisphere with up to 6 individually oriented identical axisymmetric plumes, a single 3-D plume, Shading surfaces may be used in any of 7 shapes, and a conical limit may be defined for the plume to set individual line-of-signt limits. Intermediate coordinate systems may specified to simplify input of plumes and shading surfaces.
Results of an experimental and analytical study of the nonlinear response and failure characteristics of internally pressurized 4- to 16-ply-thick graphite-epoxy cylindrical panels are presented. Specimens with clamped boundaries simulating the skin between two frames and two stringers of a typical transport fuselage were tested to failure. Failure results of aluminum specimens are compared with the graphite-epoxy test results. The specimens failed at their edges where the local bending gradients and interlaminar stresses are maximum. STAGS nonlinear two-dimensional shell analysis computer code results are used to identify regions of the panels where the response is independent of the axial coordinate. A geometrically nonlinear one-dimensional cylindrical panel analysis was derived and used to determine panel response and interlaminar stresses. Inclusion of the geometric nonlinearity was essential for accurate prediction of panel response. The maximum stress failure criterion applied to the predicted tensile stress in the fiber direction agreed best with the experimentally determined first damage pressures.
The development of time-dependent numerical simulations of unsteady interactive flows of an aerodynamic nature is reviewed with emphasis on compressible flows at flight Reynolds numbers and noniterative schemes based on Navier-Stokes equations. The importance of writing the equations in strong conservation-law form for a generalized body-oriented coordinate system is pointed out. The discussion covers time and length scales and numerical methods currently in use. Some computed results are presented and compared with experimental data.
The Environment for Application Software Integration and Execution (EASIE) provides both a methodology and a set of software utility programs to ease the task of coordinating engineering design and analysis codes. This document provides necessary information for installing the EASIE software on a host computer system. The target host is a DEX VAX running VMS version 4; host dependencies are noted when appropriate. Relevant directories and individual files are identified, and compile/load/execute sequences are specified. In the case of the data management utilities, database management system (DBMS) specific features are described in an effort to assist the maintenance programmer in converting to a new DBMS. The document also describes a sample EASIE program directory structure to guide the program implementer in establishing his/her application dependent environment.
Progress is reported on the high-resolution three-dimensional numerical simulation of flows characteristic of the Antarctic wintertime stratosphere. The numerical model is a modified version of the Reading University sigma-coordinate used previously for tropospheric studies. Physical parameterizations are kept to a minimum in order to concentrate as much computing power as possible on simulating details of the dynamical processes. The major question addressed is whether the features observed in recent high-resolution two-dimensional simulations - namely: (1) the formation of a sharp edge to the vortex (seen in the potential vorticity field), (2) the survival of the polar vortex in a material entity, and (3) the formation of small-scale eddies rough the break-up of tongues of high potential vorticity drawn out from the polar vortex - are realized in three-dimensional simulations.
A fundamental problem which must be resolved in virtually all non-trivial robotic operations is the well-known inverse kinematic question. More specifically, most of the tasks which robots are called upon to perform are specified in Cartesian (x,y,z) space, such as simple tracking along one or more straight line paths or following a specified surfacer with compliant force sensors and/or visual feedback. In all cases, control is actually implemented through coordinated motion of the various links which comprise the manipulator; i.e., in link space. As a consequence, the control computer of every sophisticated anthropomorphic robot must contain provisions for solving the inverse kinematic problem which, in the case of simple, non-redundant position control, involves the determination of the first three link angles, theta sub 1, theta sub 2, and theta sub 3, which produce a desired wrist origin position P sub xw, P sub yw, and P sub zw at the end of link 3 relative to some fixed base frame. Researchers outline a new inverse kinematic solution and demonstrate its potential via some recent computer simulations. They also compare it to current inverse kinematic methods and outline some of the remaining problems which will be addressed in order to render it fully operational. Also discussed are a number of practical consequences of this technique beyond its obvious use in solving the inverse kinematic question.
General simulation equations are derived for the rigid body motion of slung load systems. These systems are viewed as consisting of several rigid bodies connected by straight-line cables or links. The suspension can be assumed to be elastic or inelastic, both cases being of interest in simulation and control studies. Equations for the general system are obtained via D'Alembert's principle and the introduction of generalized velocity coordinates. Three forms are obtained. Two of these generalize previous case-specific results for single helicopter systems with elastic or inelastic suspensions. The third is a new formulation for inelastic suspensions. It is derived from the elastic suspension equations by choosing the generalized coordinates so as to separate motion due to cable stretching from motion with invariant cable lengths. The result is computationally more efficient than the conventional formulation, and is readily integrated with the elastic suspension formulation and readily applied to the complex dual lift and multilift systems. Equations are derived for dual lift systems. Three proposed suspension arrangements can be integrated in a single equation set. The equations are given in terms of the natural vectors and matrices of three-dimensional rigid body mechanics and are tractable for both analysis and programming.
A new and general method is presented for determining pole orientations of unresolved solar system bodies that exhibit brightness variation with rotation. The accuracy of the method is limited only by the number of available apparitions, which imposes the degree where the Legendre expansion must be truncated. It is shown that power spectra of lightcurves provide enough information to extract pole coordinates from the data, independent of body shape and albedo variegation. The method has been tested using several computer-generated synthetic light-curves, and the accuracy is found to be about 5 deg in both the longitude and latitude of the pole, which is considerably better than that obtained from the standard triaxial ellipsoid model. Finally, the new spherical harmonics method has been applied to six well-observed asteroids: 16 Psyche, 22 Kalliope, 39 Laetitia, 44 Nysa, 45 Eugenia, and 349 Dembowska.
Detailed understanding of heat transfer and fluid flow is required for many aerospace thermal systems. These systems often include phase change and operate over a range of accelerations or effective gravitational fields. An approach to analyzing such systems is presented which requires the simultaneous solution of the conservation laws of energy, momentum, and mass, as well as an equation of state. The variable property form of the governing equations are developed in two-dimensional Cartesian coordinates for a Newtonian fluid. A numerical procedure for solving the governing equations is presented and implemented in a computer program. The Galerkin form of the finite element method is used to solve the spatial variation of the field variables, along with the implicit Crank-Nicolson time marching algorithm. Quadratic Langrangian elements are used for the internal energy and the two components of velocity. Linear Lagrangian elements are used for the pressure. The location of the solid/liquid interface as well as the temperatures are determined form the calculated internal energy and pressure. This approach is quite general in that it can describe heat transfer without phase change, phase change with a sharp interface, and phase change without an interface. Analytical results from this model are compared to those of other researchers studying transient conduction, convection, and phase change and are found to be in good agreement. The numerical procedure presented requires significant computer resources, but this is not unusual when compared to similar studies by other researchers. Several methods are suggested to reduce the computational times.
A time-accurate solution method for the incompressible Navier-Stokes equations in generalized moving coordinates is presented. A finite-volume discretization method that satisfies the geometric conservation laws for time-varying computational cells is used. The discrete equations are solved by a fractional-step solution procedure. The solution is second-order-accurate in space and first-order-accurate in time. The pressure and the volume fluxes are chosen as the unknowns to facilitate the formulation of a consistent Poisson equation and thus to obtain a robust Poisson solver with favorable convergence properties. The method is validated by comparing the solutions with other numerical and experimental results. Good agreement is obtained in all cases.
"Compressed symbology" emerging technology involving one- and two-dimensional arrays of surface depressions to form optically readable dots. Patterns more durable and denser than common bar codes. Convey identification data in binary form and read by optoelectric sensors. Computers and compressed-symbology engraving machines they control constitute subsystems of "paperless" hardware-tracking and -identification systems coordinating flows of both identifying information and identified parts themselves, along with ancillary information like work orders. Modifications of software expected to accelerate marking operations, eliminate need for trial or practice marking, and reduce incidence of errors.
The Center-TRACON Automation System (CTAS), an advanced new air traffic automation program, requires knowledge of spatial and temporal atmospheric conditions such as the wind speed and direction, the temperature and the pressure in order to accurately predict aircraft trajectories. Real-time atmospheric data is available in a grid format so that CTAS must interpolate between the grid points to estimate the atmospheric parameter values. The atmospheric data grid is generally not in the same coordinate system as that used by CTAS so that coordinate conversions are required. Both the interpolation and coordinate conversion processes can introduce errors into the atmospheric data and reduce interpolation accuracy. More accurate algorithms may be computationally expensive or may require a prohibitively large amount of data storage capacity so that trade-offs must be made between accuracy and the available computational and data storage resources. The atmospheric data acquisition and processing employed by CTAS will be outlined in this report. The effects of atmospheric data processing on CTAS trajectory prediction will also be analyzed, and several examples of the trajectory prediction process will be given.