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.
Potential sources of error in the use of heated surface gages for separated-skin-friction measurement are studied. Emphasis is placed on the interpretation of local skin-friction measurements in two- and three-dimensional separated turbulent-shear flows before they are used to test the validity of current and proposed computer codes.
Accurate approximations are presented for the self-similar structures of nonradiating blast waves with adiabatic ions, isothermal electrons, and equation ion and electron temperatures at the shock. The cases considered evolve in cavities with power law ambient densities (including the uniform density case) and have negligible external pressure. The results provide the early time asymptote for systems with shock heating of electrons and strong thermal conduction. In addition, they provide analytical results against which two fluid numerical hydrodynamic codes can be checked.
A method for numerically solving the three dimensional unsteady Euler equations using flux vector splitting is developed. The equations are cast in curvilinear coordinates and a finite volume discretization is used. An explicit upwind second-order predictor-corrector scheme is used to solve the discretized equations. The scheme is stable for a CFL number of two and local time stepping is used to accelerate convergence for steady-state problems. Characteristic variable boundary conditions are developed and used in the far field and at surfaces. No additional dissipation terms are included in the scheme. Numerical results are compared with results from an existing three dimensional Euler code and experimental data.
A method for numerically solving the three-dimensional unsteady Euler equations using flux vector splitting is developed. The equations are cast in curvilinear coordinates and a finite volume discretization is used. An explicit upwind second-order predictor-corrector scheme is used to solve the discretized equations. The scheme is stable for a CFL number of 2 and local time stepping is used to accelerate convergence for steady-state problems. Characteristic variable boundary conditions are developed and used in the far-field and at surfaces. No additional dissipation terms are included in the scheme. Numerical results are compared with results from an existing three-dimensional Euler code and experimental data.
The design of a very fast, automatic black-box code for homogeneous, gas-phase chemical kinetics problems requires an understanding of the physical and numerical sources of computational inefficiency. Some major sources reviewed in this report are stiffness of the governing ordinary differential equations (ODE's) and its detection, choice of appropriate method (i.e., integration algorithm plus step-size control strategy), nonphysical initial conditions, and too frequent evaluation of thermochemical and kinetic properties. Specific techniques are recommended (and some advised against) for improving or overcoming the identified problem areas. It is argued that, because reactive species increase exponentially with time during induction, and all species exhibit asymptotic, exponential decay with time during equilibration, exponential-fitted integration algorithms are inherently more accurate for kinetics modeling than classical, polynomial-interpolant methods for the same computational work. But current codes using the exponential-fitted method lack the sophisticated stepsize-control logic of existing black-box ODE solver codes, such as EPISODE and LSODE. The ultimate chemical kinetics code does not exist yet, but the general characteristics of such a code are becoming apparent.
Based on the RECERT Program findings, NASA-Langley funded a fatigue study of code-unacceptable welds. Usage curves were developed which were based on the structural integrity of the welds. The details of this study are presented in NASA CR-178114. The information presented is a condensation and reinterpretation of the information in NASA CR-178114. This condensation and reinterpretation generated usage curves for welds having: (1) indications 0.20 -inch deep by 0.40-inch long, and (2) indications 0.195-inch deep by 8.4-inches long. These curves were developed using the procedures used in formulating the design curves in Section VIII, Division 2 of the American Society of Mechanical Engineers Boiler and Pressure Vessel Code.
Although the number of possible applications of boundary-layer control is large, a discussion is given only of those that have received the most attention recently at NASA Langley Research Center to improve airfoil drag characteristics. This research concerns stabilizing the laminar boundary layer through geometric shaping (natural laminar flow, NLF) and active control involving the removal of a portion of the laminar boundary layer (laminar flow control, LFC) either through discrete slots or a perforated surface. At low Reynolds numbers, a combination of shaping and forced transition has been used to achieve the desired run of laminar flow and control of laminar separation. In the design of both natural laminar flow and laminar flow control airfoils and wings, boundary layer stability codes play an important role. A discussion of some recent stability calculations using both incompressible and compressible codes is given.
Theoretical models of sensory-information processing by the human brain are reviewed from a human-factors perspective, with a focus on their implications for aircraft and avionics design. The topics addressed include perception (signal detection and selection), linguistic factors in perception (context provision, logical reversals, absence of cues, and order reversals), mental models, and working and long-term memory. Particular attention is given to decision-making problems such as situation assessment, decision formulation, decision quality, selection of action, the speed-accuracy tradeoff, stimulus-response compatibility, stimulus sequencing, dual-task performance, task difficulty and structure, and factors affecting multiple task performance (processing modalities, codes, and stages).
The latest version of Global Reference Atmospheric Model (GRAM-99) is presented and discussed. GRAM-99 uses either (binary) Global Upper Air Climatic Atlas (GUACA) or (ASCII) Global Gridded Upper Air Statistics (GGUAS) CD-ROM data sets, for 0-27 km altitudes. As with earlier versions, GRAM-99 provides complete geographical and altitude coverage for each month of the year. GRAM-99 uses a specially-developed data set, based on Middle Atmosphere Program (MAP) data, for 20-120 km altitudes, and NASA's 1999 version Marshall Engineering Thermosphere (MET-99) model for heights above 90 km. Fairing techniques assure smooth transition in overlap height ranges (20-27 km and 90-120 km). GRAM-99 includes water vapor and 11 other atmospheric constituents (O3, N2O, CO, CH4, CO2, N2, O2, O, A, He and H). A variable-scale perturbation model provides both large-scale (wave) and small-scale (stochastic) deviations from mean values for thermodynamic variables and horizontal and vertical wind components. The small-scale perturbation model includes improvements in representing intermittency ("patchiness"). A major new feature is an option to substitute Range Reference Atmosphere (RRA) data for conventional GRAM climatology when a trajectory passes sufficiently near any RRA site. A complete user's guide for running the program, plus sample input and output, is provided. An example is provided for how to incorporate GRAM-99 as subroutines in other programs (e.g., trajectory codes).
This document presents Mars Global Reference Atmospheric Model 2001 Version (Mars-GRAM 2001) and its new features. As with the previous version (mars-2000), all parameterizations fro temperature, pressure, density, and winds versus height, latitude, longitude, time of day, and season (Ls) use input data tables from NASA Ames Mars General Circulation Model (MGCM) for the surface through 80-km altitude and the University of Arizona Mars Thermospheric General Circulation Model (MTGCM) for 80 to 70 km. Mars-GRAM 2001 is based on topography from the Mars Orbiter Laser Altimeter (MOLA) and includes new MGCM data at the topographic surface. A new auxiliary program allows Mars-GRAM output to be used to compute shortwave (solar) and longwave (thermal) radiation at the surface and top of atmosphere. This memorandum includes instructions on obtaining Mars-GRAN source code and data files and for running the program. It also provides sample input and output and an example for incorporating Mars-GRAM as an atmospheric subroutine in a trajectory code.
Electromagnetic software (and software in general) at one time was coded onto a series of cards. These card decks contained both code and input data.
This paper presents comparisons of seven propagation codes for predicting liner attenuation in ducts with flow. The selected codes span the spectrum of methods available (finite element, parabolic approximation, and pseudo-time domain) and are collectively representative of the state-of-art in the liner industry. These codes are included because they have two-dimensional and three-dimensional versions and can be exported to NASA's Columbia Supercomputer. The basic assumptions, governing differential equations, boundary conditions, and numerical methods underlying each code are briefly reviewed and an assessment is performed based on two predefined metrics. The two metrics used in the assessment are the accuracy of the predicted attenuation and the amount of wall clock time to predict the attenuation. The assessment is performed over a range of frequencies, mean flow rates, and grazing flow liner impedances commonly used in the liner industry. The primary conclusions of the study are (1) predicted attenuations are in good agreement for rigid wall ducts, (2) the majority of codes compare well to each other and to approximate results from mode theory for soft wall ducts, (3) most codes compare well to measured data on a statistical basis, (4) only the finite element codes with cubic Hermite polynomials capture extremely large attenuations, and (5) wall clock time increases by an order of magnitude or more are observed for a three-dimensional code relative to the corresponding two-dimensional version of the same code.
Figures of Merit (FoMs) and the FoM software provide a method for quantitatively evaluating the quality of a regolith simulant by comparing the simulant to a reference material. FoMs may be used for comparing a simulant to actual regolith material, specification by stating the value a simulant s FoMs must attain to be suitable for a given application and comparing simulants from different vendors or production runs. FoMs may even be used to compare different simulants to each other. A single FoM is conceptually an algorithm that computes a single number for quantifying the similarity or difference of a single characteristic of a simulant material and a reference material and provides a clear measure of how well a simulant and reference material match or compare. FoMs have been constructed to lie between zero and 1, with zero indicating a poor or no match and 1 indicating a perfect match. FoMs are defined for modal composition, particle size distribution, particle shape distribution, (aspect ratio and angularity), and density. This TM covers the mathematics, use, installation, and licensing for the existing FoM code in detail.
Multiple nozzle plume flow field is computed with a 3-D, Navier-Stokes solver. Numerical simulation is performed with a flux-split, two-factor, time asymptotic viscous flow solver of Ying and Steger. The two factor splitting provides a stable 3-D solution procedure under ideal-gas assumptions. An ad-hoc acceleration procedure that shows promise in improving the convergence rate by a factor of three for steady state problems is utilized. Computed solutions to generic problems at various altitude and flight conditions show flow field complexity and three-dimensional effects due to multiple nozzle jet interactions. Viscous, ideal gas solutions for the symmetric nozzle are compared with other numerical solutions.
Performance prediction of turbomachines is a significant part of aircraft propulsion design. In the conceptual design stage, there is an important need to quantify compressor and turbine aerodynamic performance and develop initial geometry parameters at the 2-D level prior to more extensive Computational Fluid Dynamics (CFD) analyses. The Object-oriented Turbomachinery Analysis Code (OTAC) is being developed to perform 2-D meridional flowthrough analysis of turbomachines using an implicit formulation of the governing equations to solve for the conditions at the exit of each blade row. OTAC is designed to perform meanline or streamline calculations; for streamline analyses simple radial equilibrium is used as a governing equation to solve for spanwise property variations. While the goal for OTAC is to allow simulation of physical effects and architectural features unavailable in other existing codes, it must first prove capable of performing calculations for conventional turbomachines.OTAC is being developed using the interpreted language features available in the Numerical Propulsion System Simulation (NPSS) code described by Claus et al (1991). Using the NPSS framework came with several distinct advantages, including access to the pre-existing NPSS thermodynamic property packages and the NPSS Newton-Raphson solver. The remaining objects necessary for OTAC were written in the NPSS framework interpreted language. These new objects form the core of OTAC and are the BladeRow, BladeSegment, TransitionSection, Expander, Reducer, and OTACstart Elements. The BladeRow and BladeSegment consumed the initial bulk of the development effort and required determining the equations applicable to flow through turbomachinery blade rows given specific assumptions about the nature of that flow. Once these objects were completed, OTAC was tested and found to agree with existing solutions from other codes; these tests included various meanline and streamline comparisons of axial compressors and turbines at design and off-design conditions.
Performance prediction of turbomachines is a significant part of aircraft propulsion design. In the conceptual design stage, there is an important need to quantify compressor and turbine aerodynamic performance and develop initial geometry parameters at the 2-D level prior to more extensive Computational Fluid Dynamics (CFD) analyses. The Object-oriented Turbomachinery Analysis Code (OTAC) is being developed to perform 2-D meridional flowthrough analysis of turbomachines using an implicit formulation of the governing equations to solve for the conditions at the exit of each blade row. OTAC is designed to perform meanline or streamline calculations; for streamline analyses simple radial equilibrium is used as a governing equation to solve for spanwise property variations. While the goal for OTAC is to allow simulation of physical effects and architectural features unavailable in other existing codes, it must first prove capable of performing calculations for conventional turbomachines. OTAC is being developed using the interpreted language features available in the Numerical Propulsion System Simulation (NPSS) code described by Claus et al (1991). Using the NPSS framework came with several distinct advantages, including access to the pre-existing NPSS thermodynamic property packages and the NPSS Newton-Raphson solver. The remaining objects necessary for OTAC were written in the NPSS framework interpreted language. These new objects form the core of OTAC and are the BladeRow, BladeSegment, TransitionSection, Expander, Reducer, and OTACstart Elements. The BladeRow and BladeSegment consumed the initial bulk of the development effort and required determining the equations applicable to flow through turbomachinery blade rows given specific assumptions about the nature of that flow. Once these objects were completed, OTAC was tested and found to agree with existing solutions from other codes; these tests included various meanline and streamline comparisons of axial compressors and turbines at design and off-design conditions.
An experimental research program was conducted in the NASA Lewis Research Center 10 x 10 ft supersonic wind tunnel. The 2-D inlet model was designed to study the Mach 3.0 to 5.0 speed range for an over-under turbojet plus ramjet propulsion system. The model was extensively instrumented to provide both analytical code validation data as well as inlet performance information. Support studies for the program include flow field predictions with both 3-D parabolized Navier-Stokes (PNS) and 3-D full Navier-Stokes (FNS) analytical codes. Analytical predictions and experimental results are compared.
of almost any popular CAD format, e.g. IGES, Parasolid, DXF, ACIS etc. The solid models are processed (simplified) and meshed in GiD(R), and then converted into WIPL-D Pro input file by simple Fortran or Matlab code. This algorithm allows the user to control the mesh of imported geometry, and to assign electric p~perties to metalic and dielectric surfaces. Implementation of the algorithm is demonstrated by examples obtained from the NASA Discovery mission, Phoenix Lander 2008. Results for radiation pattern of Phoenix Lander UHF relay antenna with effect of Martian surface, both simulated in WIPL-D Pro and measured, are shown for comparison.