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.
Pocket-size plastic manual tape reader device aids in reading, interpreting, and correcting binary and octal coded punched tapes. The coded information is more easily read if the color of the back plate contrasts sharply with that of the tape.
Language problem areas are summarized and possible solutions are discussed. Areas covered are: (1) double complex, (2) arrays of work space in library subroutines, (3) calls from a library subprogram to user code, (4) decomposition and synthesis of floating-point numbers, (5) environment parameters, (6) precision function, and (7) control structures.
High chamber pressure expander cycles proposed for orbit transfer vehicles depend primarily on the heat energy transmitted from the combustion products through the thrust wall chamber wall. The heat transfer to the nozzle wall is affected by such variables as wall roughness, relamarization, and the presence of particles in the flow. Motor performance loss for these nozzles with thick boundary layers is inaccurate using the existing procedure coded BLIMPJ. Modifications and innovations to the code are examined. Updated routines are listed.
This paper examines current research efforts underway at Lewis Research Center which have focused on improving the structural durability of the SSME turbodrive system through more fully examining the aerothermodynamics of the flow system. This objective was achieved by developing three-dimensional viscous and inviscid computer codes to predict the pressure and temperature distributions in the engine flow passages, particularly within the SSME turbines and turnaround ducts.
Two research programs are described which directly relate to the application of natural laminar flow (NLF) technology to transonic transport-type wind planforms. Each involved using state-of-the-art computational methods to design three-dimensional wing contours which generate significant runs of favorable pressure gradients. The first program supported the Variable Sweep Transition Flight Experiment and involves design of a full-span glove which extends from the leading edge to the spoiler hinge line on the upper surface of an F-14 outer wing panel. Boundary-layer and static-pressure data will be measured on this design during the supporting wind-tunnel and flight tests. These data will then be analyzed and used to infer the relationship between crossflow and Tollmein-Schlichting disturbances on laminar boundary-layer transition. A wing was designed computationally for a corporate transport aircraft in the second program. The resulting wing design generated favorable pressure gradients from the leading edge aft to the mid-chord on both upper and lower surfaces at the cruise design point. Detailed descriptions of the computational design approach are presented along with the various constraints imposed on each of the designs. Wing surface pressure distributions, which support the design objective and were derived from transonic three-dimensional analyses codes, are also presented. Current status of each of the research programs is included in the summary.
An evolving system for CCD data acquisition is described. The electronics are subdivided into those attached to the dewar, a 'smart' controller near the dewar, and a computer connected by serial link to the smart controller. Software for the smart controller is in assembler code, while the data acquisition and on-line analysis software is written in 'C' and uses the UNIX operating system.
The General Interpolants Method (GIM) is a 3-D, time-dependent, hybrid procedure for generating numerical analogs of the conservation laws. This study is directed toward the development and application of the GIM computer code for fluid dynamic research applications as implemented for the Cyber 200 series of supercomputers. An elliptic and quasi-parabolic version of the GIM code are discussed. Turbulence models, algebraic and differential equations, were added to the basic viscous code. An equilibrium reacting chemistry model and an implicit finite difference scheme are also included.
A generalized one-dimensional computer code for analyzing the flow and heat transfer in the turbomachinery cooling passages was developed. This code is capable of handling rotating cooling passages with turbulators, 180 degree turns, pin fins, finned passages, by-pass flows, tip cap impingement flows, and flow branching. The code is an extension of a one-dimensional code developed by P. Meitner. In the subject code, correlations for both heat transfer coefficient and pressure loss computations were developed to model each of the above mentioned type of coolant passages. The code has the capability of independently computing the friction factor and heat transfer coefficient on each side of a rectangular passage. Either the mass flow at the inlet to the channel or the exit plane pressure can be specified. For a specified inlet total temperature, inlet total pressure, and exit static pressure, the code computers the flow rates through the main branch and the subbranches, flow through tip cap for impingement cooling, in addition to computing the coolant pressure, temperature, and heat transfer coefficient distribution in each coolant flow branch. Predictions from the subject code for both nonrotating and rotating passages agree well with experimental data. The code was used to analyze the cooling passage of a research cooled radial rotor.
A generalized one-dimensional computer code for analyzing the flow and heat transfer in the turbomachinery cooling passages was developed. This code is capable of handling rotating cooling passages with turbulators, 180 degree turns, pin fins, finned passages, by-pass flows, tip cap impingement flows, and flow branching. The code is an extension of a one-dimensional code developed by P. Meitner. In the subject code, correlations for both heat transfer coefficient and pressure loss computations were developed to model each of the above mentioned type of coolant passages. The code has the capability of independently computing the friction factor and heat transfer coefficient on each side of a rectangular passage. Either the mass flow at the inlet to the channel or the exit plane pressure can be specified. For a specified inlet total temperature, inlet total pressure, and exit static pressure, the code computers the flow rates through the main branch and the subbranches, flow through tip cap for impingement cooling, in addition to computing the coolant pressure, temperature, and heat transfer coefficient distribution in each coolant flow branch. Predictions from the subject code for both nonrotating and rotating passages agree well with experimental data. The code was used to analyze the cooling passage of a research cooled radial rotor.
Jet-A spray combustion has been evaluated in gas turbine combustion with the use of propane chemical kinetics as the first approximation for the chemical reactions. Here, the numerical solutions are obtained by using the KIVA-2 computer code. The KIVA-2 code is the most developed of the available multidimensional combustion computer programs for application of the in-cylinder combustion dynamics of internal combustion engines. The released version of KIVA-2 assumes that 12 chemical species are present; the code uses an Arrhenius kinetic-controlled combustion model governed by a four-step global chemical reaction and six equilibrium reactions. Researchers efforts involve the addition of Jet-A thermophysical properties and the implementation of detailed reaction mechanisms for propane oxidation. Three different detailed reaction mechanism models are considered. The first model consists of 131 reactions and 45 species. This is considered as the full mechanism which is developed through the study of chemical kinetics of propane combustion in an enclosed chamber. The full mechanism is evaluated by comparing calculated ignition delay times with available shock tube data. However, these detailed reactions occupy too much computer memory and CPU time for the computation. Therefore, it only serves as a benchmark case by which to evaluate other simplified models. Two possible simplified models were tested in the existing computer code KIVA-2 for the same conditions as used with the full mechanism. One model is obtained through a sensitivity analysis using LSENS, the general kinetics and sensitivity analysis program code of D. A. Bittker and K. Radhakrishnan. This model consists of 45 chemical reactions and 27 species. The other model is based on the work published by C. K. Westbrook and F. L. Dryer.
A simple technique is presented which forms a cubic grid model of a helicopter from an Aircraft Modeling Code (AMC) input file. The AMC input file defines the helicopter fuselage as a series of polygonal cross sections. The cubic grid model is used as an input to a Finite Difference Time Domain (FDTD) code to obtain predictions of antenna performance on a generic helicopter model. The predictions compare reasonably well with measured data.
Three-dimensional shielding effects on cosmic ray charged particle fluences were measured with plastic nuclear track detectors in the P0006 experiment on Long Duration Exposure Facility (LDEF). The azimuthal and polar angle distributions of the galactic cosmic ray particles (mostly relativistic iron) were measured in the main stack and in four side stacks of the P0006 experiment, located on the west end of the LDEF satellite. A shadowing effect of the shielding of the LDEF satellite is found. Total fluence of stopping protons was measured as a function of the position in the main and side stacks of the P0006 experiment. Location dependence of total track density is explained by the three-dimensional shielding model of the P0006 stack. These results can be used to validate 3D mass model and transport code calculations and also for predictions of the outer radiation environment for the Space Station Freedom.
A nonperturbative analytic solution of the high charge and energy (HZE) Green's function is used to implement a computer code for laboratory ion beam transport. The code is established to operate on the Langley Research Center nuclear fragmentation model used in engineering applications. Computational procedures are established to generate linear energy transfer (LET) distributions for a specified ion beam and target for comparison with experimental measurements. The code is highly efficient and compares well with the perturbation approximations.
A nonperturbative analytic solution of the high charge and energy (HZE) Green's function is used to implement a computer code for laboratory ion beam transport in multilayered materials. The code is established to operate on the Langley nuclear fragmentation model used in engineering applications. Computational procedures are established to generate linear energy transfer (LET) distributions for a specified ion beam and target for comparison with experimental measurements. The code was found to be highly efficient and compared well with the perturbation approximation.
A brief overview of the transient dynamics capabilities at Sandia National Laboratories, with an emphasis on recent new developments and current research is presented. In addition, the Sandia National Laboratories (SNL) Engineering Analysis Code Access System (SEACAS), which is a collection of structural and thermal codes and utilities used by analysts at SNL, is described. The SEACAS system includes pre- and post-processing codes, analysis codes, database translation codes, support libraries, Unix shell scripts for execution, and an installation system. SEACAS is used at SNL on a daily basis as a production, research, and development system for the engineering analysts and code developers. Over the past year, approximately 190 days of CPU time were used by SEACAS codes on jobs running from a few seconds up to two and one-half days of CPU time. SEACAS is running on several different systems at SNL including Cray Unicos, Hewlett Packard PH-UX, Digital Equipment Ultrix, and Sun SunOS. An overview of SEACAS, including a short description of the codes in the system, are presented. Abstracts and references for the codes are listed at the end of the report.
GMIS is a general-purpose simulator for testing ground system software. GMIS can be adapted to any mission to simulate changes in the data state maintained by the mission's computers. GMIS was developed in Code 522 NASA Goddard Space Flight Center. The acronym GMIS stands for GOTT Mission Independent Simulator, where GOTT is the Ground Operations Technology Testbed. Within GOTT, GMIS is used to provide simulated data to an installation of TPOCC - the Transportable Payload Operations Control Center. TPOCC was developed by Code 510 as a reusable control center. GOTT uses GMIS and TPOCC to test new technology and new operator procedures.
A computational simulation procedure is presented for nonlinear analyses which incorporates microstress redistribution due to progressive fracture in ceramic matrix composites. This procedure facilitates an accurate simulation of the stress-strain behavior of ceramic matrix composites up to failure. The nonlinearity in the material behavior is accounted for at the constituent (fiber/matrix/interphase) level. This computational procedure is a part of recent upgrades to CEMCAN (Ceramic Matrix Composite Analyzer) computer code. The fiber substructuring technique in CEMCAN is used to monitor the damage initiation and progression as the load increases. The room-temperature tensile stress-strain curves for SiC fiber reinforced reaction-bonded silicon nitride (RBSN) matrix unidirectional and angle-ply laminates are simulated and compared with experimentally observed stress-strain behavior. Comparison between the predicted stress/strain behavior and experimental stress/strain curves is good. Collectively the results demonstrate that CEMCAN computer code provides the user with an effective computational tool to simulate the behavior of ceramic matrix composites.
During the all sky survey (May 1991 - Nov. 1992) of the Compton Gamma Ray Observatory the Vela pulsar PSR0833-45 was in the field of view of the Energetic Gamma Ray Experiment Telescope (EGRET) in ten separate viewing periods. The pulsar was detected in each one. The average intensity from 100 MeV to 2 GeV was (7.8 +/- 1.0) x 10 (exp -6) photons cm(exp -2) s(exp -1), which indicates that the pulsar in the years 1991/92 was in a state comparable to the low fluxes observed in 1977-1980. No significant changes in intensity were detected during the EGRET observations. The total spectrum of PSR0833-45 measured by EGRET can be described by a power law with index -(1.70 +/- 0.02) over the range 30 MeV to 2 GeV. The extrapolation of this spectrum into the 3 to 30 MeV range agrees with the observations by COMPTEL. Above 2 to 4 GeV EGRET detects a strong spectral break. The lightcurves obtained show a familar structure in the phase histogram: two peaks separated by 0.424 +/- 0.002 in phase with considerable emission in the phase interval between the peaks. The first gamma ray peak maximum trails the single radio peak maximum by 10.54 +/- 0.09 ms (= 0.118 +/- 0.001 in phase). The widths of the emission peaks (FWHM) are 2.7 ms for the first peak (0.03 phase) and 4.1 ms for the second peak (0.05 phase). The widths are approximately constant below a GeV, but show a tendency to become narrower at higher energies. On Jul. 20 1991 a glitch of the Vela period was registered in monitor radio observations. No significant differences between the pre- and post-glitch gamma ray lightcurves were found. The statistics available for the Vela observations allow for a division of the lightcurve into eight phase intervals. The emission peak cores (central FWHM) with leading and trailing wings and two interval regions were defined and spectra were derived for all parts of the lightcurve. The energy spectra for the eight phase intervals show significant differences: the first peak (approximately E(exp -1.81 +/- 0.04)) is somewhat softer than the second peak (approximately E(exp -1.71 +/- 0.03)); the wings attached to each peak show softer spectra than the code of the peaks; the interval emission has the hardest spectrum (approximately E(exp -1.52 +/- 0.03)).