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 programming manual for the two-dimensional Boeing Plastic Analysis Capability for Engines (BOPACE) program written in FORTRAN IV is presented. The BOPACE flow logic is summarized along with subroutines and input/output files. Emphasis is placed on linear equation solver routines written as an independent package.
Two plate elements, the linear strain triangular membrane element CTRIM6 and the higher order plate bending element CTRPLT1, were added to NASTRAN Level 16.0. The theoretical formulation, programming details, and bulk data information pertaining to the addition of these elements are discussed. Sample problems illustrating the use of these elements are presented.
A general version of Lagrange's equations, including auxiliary nonholonomic, rheonomic conditions of constraint, is used in the dynamic simulation and stability analysis of interconnected flexible bodies. Modeling of the nonlinear flexible/rigid dynamic coupling effects, the interaction forces/torques, and the elastic deformation effects is discussed. A digital computer program is developed to obtain time-domain solution for the nonlinear response of systems represented as a collection of individual bodies, numerical linearization of system-governing equations, time-domain solution for the perturbation response about a nominal state, and a frequency-domain stability analysis corresponding to the linearization. The digital simulation code is employed to study the dynamic behavior of a typical satellite and a spacecraft with deployable experiment booms.
Subroutine PLTMG is a FORTRAN program for solving self-adjoint elliptic boundary value problems in general regions of R-squared. It is based on a piecewise linear triangle finite element method, an adaptive grid refinement procedure, and a multi-level iterative method to solve the resulting sets of linear equations. This paper describes the method and presents some numerical results and comparisons.
A closed form solution is presented for the coplanar transfer between nearby circular orbits for spacecraft using their own blow down propulsion system. The decaying thrust is applied along the local horizontal and the linearized equation of motion in the orbital elements formulation are used to describe the transfer which consists of a thrust-coast-relight program. Sensitivity partials are also presented analytically in order to study the effect of maneuver execution errors and various other parameters affecting the blow down propulsion system characteristics on the transfer. This strategy is applied to study the transfer of the TOPEX spacecraft from the Shuttle parking orbit to its final operational orbit.
The demand for cartographic products at scales of 1:25,000 to 1:250,000 continues to increase throughout the world in order to meet requirements associated with: the survey and management of natural resources; environmental planning; and the establishment of geo-referenced data bases. A satellite system involving the use of multispectral linear array (MLA) sensors designed to meet cartographic requirements in terms of the completeness of detail and geometric accuracy standards associated with mapping programs offers great promise for rapidly providing the data with which to produce four types of map products/and data. These include topographic maps; digital terrain information (x,y,z coordinates; thematic mapping; and image maps. Each type of cartographic product/data is briefly considered in relation to current needs.
Some exactly soluble families of models for the chemical evolution of the Galaxy are presented. The parameters considered include gas mass, the age-metallicity relation, the star mass vs. metallicity, the age distribution, and the mean age of dwarfs. A short BASIC program for calculating these parameters is given. The calculation of metallicity gradients, nuclear cosmochronology, and extinct radioactivities is addressed. An especially simple, mathematically linear model is recommended as a standard model of galaxies with truncated infall due to its internal consistency and compact display of the physical effects of the parameters.
The design, implementation, and performance evaluation of a message passing facility (MPF) for shared memory multiprocessors are presented. The MPF is based on a message passing model conceptually similar to conversations. Participants (parallel processors) can enter or leave a conversation at any time. The message passing primitives for this model are implemented as a portable library of C function calls. The MPF is currently operational on a Sequent Balance 21000, and several parallel applications were developed and tested. Several simple benchmark programs are presented to establish interprocess communication performance for common patterns of interprocess communication. Finally, performance figures are presented for two parallel applications, linear systems solution, and iterative solution of partial differential equations.
An efficient finite element program for 3-D inelastic analysis of gas turbine hot section components was constructed and validated. A novel mixed iterative solution strategy is derived from the augmented Hu-Washizu variational principle in order to nodally interpolate coordinates, displacements, deformation, strains, stresses and material properties. A series of increasingly sophisticated material models incorporated in MHOST include elasticity, secant plasticity, infinitesimal and finite deformation plasticity, creep and unified viscoplastic constitutive model proposed by Walker. A library of high performance elements is built into this computer program utilizing the concepts of selective reduced integrations and independent strain interpolations. A family of efficient solution algorithms is implemented in MHOST for linear and nonlinear equation solution including the classical Newton-Raphson, modified, quasi and secant Newton methods with optional line search and the conjugate gradient method.
The Primary basis for heat transfer analysis of turbine airfoils is experimental data obtained in linear cascades. These data were very valuable in identifying the major heat transfer and fluid flow features of a turbine airfoil. The first program objective is to obtain a detailed set of heat transfer coefficients along the midspan of a stator and a rotor in a rotating turbine stage. The data are to be compared to some standard analysis of blade boundary layer heat transfer which is in use today. A second program objective is to obtain a detailed set of heat transfer coefficients along the midspan of a stator located in the wake of an upstream turbine stage.
This paper summarizes the development of the methods and a computer program to compute the probability of instability of a dynamic system than can be represented by a system of second-order ordinary linear differential equations. Two instability criteria based upon the roots of the characteristics equation or Routh-Hurwitz test functions are investigated. Computational methods based on system reliability analysis methods and importance sampling concepts are proposed to perform efficient probabilistic analysis. Numerical examples are provided to demonstrate the methods.
Development of design tools to furnish optimal acoustic environments for lightweight aircraft demands the ability to simulate the acoustic system on a workstation. In order to form an effective mathematical model of the phenomena at hand, we have begun by studying the propagation of acoustic waves inside closed spherical shells. Using a fully-coupled fluid-structure interaction model based upon variational principles, we have written a finite element analysis program and are in the process of examining several test cases. Future investigations are planned to increase model accuracy by incorporating non-linear and viscous effects.
In July of 1999 two linear aerospike rocket engines will power the first flight of NASA's X-33 advanced technology demonstrator. A successful X-33 flight test program will validate the aerospike nozzle concept, a key technical feature of Lockheed Martin's VentureStar(trademark) reusable launch vehicle. The aerospike received serious consideration for NASA's current space shuttle, but was eventually rejected in 1969 in favor of high chamber pressure bell engines, in part because of perceived technical risk. The aerospike engine (discussed below) has several performance advantages over conventional bell engines. However, these performance advantages are difficult to validate by ground test. The space shuttle, a multibillion dollar program intended to provide all of NASA's future space lift could not afford the gamble of choosing a potentially superior though unproven aerospike engine over a conventional bell engine. The X-33 demonstrator provides an opportunity to prove the aerospike's performance advantage in flight before commiting to an operational vehicle.
The frozen-flux hypothesis for the Earth's liquid core assumes that convective terms dominate diffusive terms in the induction equation governing the behaviour of the magnetic field at the surface of the core. While highly plausible on the basis of estimates of physical parameters, the hypothesis has been questioned. To study this hypothesis, we improve the method which tests the consistency of magnetic observations with the hypothesis by constructing simple, flux-conserving core-field models fitting the data at pairs of epochs. We introduce a new approach that fixes the patch configurations at each of the two epochs before inversion, so that each configuration is consistent with its respective data set but possesses the same patch topology. We expand upon the inversion algorithm, using quadratic programming to maintain the proper flux sign within patches; the modelling calculations are also extended to include data types that depend non-linearly on the model. Every test of a hypothesis depends on the characterization of the observational uncertainties; we undertake a thorough review of this question. For main-field models, the primary source of uncertainty comes from the crustal field. We base our analysis on statistical models of the crustal magnetization, adjusted to bring it into better conformity with our data set. The noise model permits us to take into account the correlations between the measurements and requires that a different weighting be given to horizontal and vertical components. It also indicates that the observations should be fit more closely than has been the practice heretofore. We apply the revised method to Magsat data from 1980 and survey and observatory data from 1915.5, two data sets believed to be particularly difficult to reconcile with the frozen-flux hypothesis. We compute a pair of simple, flux-conserving models that fit the averaged data from each epoch. We therefore conclude that present knowledge of the geomagnetic fields of 1980 and 1915.5 is consistent with the frozen-flux hypothesis.
An experimental program at Lewis Research Center investigates acoustical-mode combustion instability in liquid propellant rockets. One phase of this program is concerned with nitrogen tetroxide and a 50-50 fuel blend of hydrazine-UDMH. Effects on combustor stability and performance by variations in injection velocities, impingement angle and distance, amd thrust per element of triplet injectors were studied in a 10.77-inch- diameter cylindrical combustor at chamber pressures of 100 and 300 psia. Nominal thrust levels were 6700 and 20 000 pounds. Stability rating was accomplished by the use of various size RDX explosive charges to generate tangential pressure disturbances. The injection velocity effect correlated by using a ratio V(sub o)/V(sub f). Stability correlation with V(sub o)/V(sub f) resulted in a sharply humped curve with maximum stability at a ratio of 1.2. Performance decreased with increase in V(sub o)/V(sub f) after a maximum at 0.7, which agreed fairly well with the "uniform mixture ratio distribution" criteria of the Jet Propulsion Laboratory. Variation of impingement angle from 38 deg to 120 deg indicated that performance correlated with the absolute velocity ratio V(sub o)/V(sub f), but stability correlation appeared to be best with a velocity ratio using the radial component of the fuel injection velocity V(sub o)/V(sub fr). Impingement distances from 0.5 to 1.0 inch appeared to have only slight effect on stability and performance. Investigation of injector thrust per element for V(sub o)/V(sub fr) ranging from 0.8 to 2.8 yielded convex maxima curves. Acoustic liners were also tested to improve a liner design computer program, to obtain design criteria, and to develop flight-type liners. A marginally stable and spontaneously unstable combustor were used to evaluate linear configurations. All liners were effective in the marginally stable configuration, and only liners with large theoretical absorptivity were successful in the spontaneously unstable combustor.
The shared memory Multi-Level Parallelism (MLP) technique, developed last year at NASA Ames has been very successful in dramatically improving the performance of important NASA CFD codes. This new and very simple parallel programming technique was first inserted into the OVERFLOW production CFD code in FY 1998. The OVERFLOW-MLP code's parallel performance scaled linearly to 256 CPUs on the NASA Ames 256 CPU Origin 2000 system (steger). Overall performance exceeded 20.1 GFLOP/s, or about 4.5x the performance of a dedicated 16 CPU C90 system. All of this was achieved without any major modification to the original vector based code. The OVERFLOW-MLP code is now in production on the inhouse Origin systems as well as being used offsite at commercial aerospace companies. Partially as a result of this work, NASA Ames has purchased a new 512 CPU Origin 2000 system to further test the limits of parallel performance for NASA codes of interest. This paper presents the performance obtained from the latest optimization efforts on this machine for the LAURA-MLP and OVERFLOW-MLP codes. The Langley Aerothermodynamics Upwind Relaxation Algorithm (LAURA) code is a key simulation tool in the development of the next generation shuttle, interplanetary reentry vehicles, and nearly all "X" plane development. This code sustains about 4-5 GFLOP/s on a dedicated 16 CPU C90. At this rate, expected workloads would require over 100 C90 CPU years of computing over the next few calendar years. It is not feasible to expect that this would be affordable or available to the user community. Dramatic performance gains on cheaper systems are needed. This code is expected to be perhaps the largest consumer of NASA Ames compute cycles per run in the coming year.The OVERFLOW CFD code is extensively used in the government and commercial aerospace communities to evaluate new aircraft designs. It is one of the largest consumers of NASA supercomputing cycles and large simulations of highly resolved full aircraft are routinely undertaken. Typical large problems might require 100s of Cray C90 CPU hours to complete. The dramatic performance gains with the 256 CPU steger system are exciting. Obtaining results in hours instead of months is revolutionizing the way in which aircraft manufacturers are looking at future aircraft simulation work. Figure 2 below is a current state of the art plot of OVERFLOW-MLP performance on the 512 CPU Lomax system. As can be seen, the chart indicates that OVERFLOW-MLP continues to scale linearly with CPU count up to 512 CPUs on a large 35 million point full aircraft RANS simulation. At this point performance is such that a fully converged simulation of 2500 time steps is completed in less than 2 hours of elapsed time. Further work over the next few weeks will improve the performance of this code even further.The LAURA code has been converted to the MLP format as well. This code is currently being optimized for the 512 CPU system. Performance statistics indicate that the goal of 100 GFLOP/s will be achieved by year's end. This amounts to 20x the 16 CPU C90 result and strongly demonstrates the viability of the new parallel systems rapidly solving very large simulations in a production environment.
This report reviews six classes of models that are used for studying transportation network topologies. The report is motivated by two main questions. First, what can the "new science" of complex networks (scale-free, small-world networks) contribute to our understanding of transport network structure, compared to more traditional methods? Second, how can geographic information systems (GIS) contribute to studying transport networks? The report defines terms that can be used to classify different kinds of models by their function, composition, mechanism, spatial and temporal dimensions, certainty, linearity, and resolution. Six broad classes of models for analyzing transport network topologies are then explored: GIS; static graph theory; complex networks; mathematical programming; simulation; and agent-based modeling. Each class of models is defined and classified according to the attributes introduced earlier. The paper identifies some typical types of research questions about network structure that have been addressed by each class of model in the literature.
A method is developed to determine the flow field of a body of revolution in separated flow. The computer was used to integrate various solutions and solution properties of the sub-flow fields which made up the entire flow field without resorting to a finite difference solution to the complete Navier-Stokes equations. The technique entails the use of the unsteady cross flow analogy and a new solution to the two-dimensional unsteady separated flow problem based upon an unsteady, discrete-vorticity wake. Data for the forces and moments on aerodynamic bodies at low speeds and high angle of attack (outside the range of linear inviscid theories) such that the flow is substantially separated are produced which compare well with experimental data. In addition, three dimensional steady separated regions and wake vortex patterns are determined. The computer program developed to perform the numerical calculations is described.