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.
New version of TEACH computer program developed specifically for use in analysis of subsonic, swirling, reacting, turbulent flow in cylindrical bluff-body research combustor. Combustor design widely used as research tool for mathematical modeling of gas turbines and for development of diagnostic instrumentation.
The ceramics analysis and reliability evaluation of structures (CARES) computer program is described. The primary function of the code is to calculate the fast-fracture reliability or failure probability of macro-scopically isotropic ceramic components. These components may be subjected to complex thermomechanical loadings, such as those found in heat engine applications. CARES uses results from MSC/NASTRAN or ANSYS finite-element analysis programs to evaluate how inherent surface and/or volume type flaws component reliability. CARES utilizes the Batdorf model and the two-parameter Weibull cumulative distribution function to describe the effects of multiaxial stress states on material strength. The principle of independent action (PIA) and the Weibull normal stress averaging models are also included. Weibull material strength parameters, the Batdorf crack density coefficient, and other related statistical quantities are estimated from four-point bend bar or uniform uniaxial tensile specimen fracture strength data. Parameter estimation can be performed for a single or multiple failure modes by using a least-squares analysis or a maximum likelihood method. Kolmogorov-Smirnov and Anderson-Darling goodness-to-fit-tests, 90 percent confidence intervals on the Weibull parameters, and Kanofsky-Srinivasan 90 percent confidence band values are also provided. Examples are provided to illustrate the various features of CARES.
NASA-Lewis has developed a public-domain computer program, designated 'Ceramic Analysis and Reliability Evaluation of Structures' (CARES) for calculating the fast-fracture reliability of macroscopically isotropic ceramic components subjected to the complex thermomechanical loadings typical of heat engines. The design methodology employed by CARES encompasses linear elastic fracture mechanics theory, extreme value statistics, and material microstructures; component integrity is conceived as a function of the entire field solution of the stresses, rather than being based solely on the most highly stressed point.
Extension of the high charge and energy (HZE) transport computer program HZETRN for angular transport of neutrons is considered. For this paper, only light ion transport, He4 and lighter, will be analyzed using a pure solar proton source. The angular transport calculator is the ANISN/PC program which is being controlled by the HZETRN program. The neutron flux values are compared for straight-ahead transport and angular transport in one dimension. The shield material is aluminum and the target material is water. The thickness of these materials is varied; however, only the largest model calculated is reported which is 50 gm/sq cm of aluminum and 100 gm/sq cm of water. The flux from the ANISN/PC calculation is about two orders of magnitude lower than the flux from HZETRN for very low energy neutrons. It is only a magnitude lower for the neutrons in the 10 to 20 MeV range in the aluminum and two orders lower in the water. The major reason for this difference is in the transport modes: straight-ahead versus angular. The angular treatment allows a longer path length than the straight-ahead approximation. Another reason is the different cross section sets used by the ANISN/PC-BUGLE-80 mode and the HZETRN mode. The next step is to investigate further the differences between the two codes and isolate the differences to just the angular versus straight-ahead transport mode. Then, create a better coupling between the angular neutron transport and the charged particle transport.
A FORTRAN computer program for anisotropic shallow-shell finite elements with variable curvature is described. A listing of the program is presented together with printed output for a sample case. Computation times and central memory requirements are given for several different elements. The program is based on a stiffness (displacement) finite-element model in which the fundamental unknowns consist of both the displacement and the rotation components of the reference surface of the shell. Two triangular and four quadrilateral elements are implemented in the program. The triangular elements have 6 or 10 nodes, and the quadrilateral elements have 4 or 8 nodes. Two of the quadrilateral elements have internal degrees of freedom associated with displacement modes which vanish along the edges of the elements (bubble modes). The triangular elements and the remaining two quadrilateral elements do not have bubble modes. The output from the program consists of arrays corresponding to the stiffness, the geometric stiffness, the consistent mass, and the consistent load matrices for individual elements. The integrals required for the generation of these arrays are evaluated by using symbolic (or analytic) integration in conjunction with certain group-theoretic techniques. The analytic expressions for the integrals are exact and were developed using the symbolic and algebraic manipulation language.
The results are presented of a study to produce a computer program to calculate laminar and turbulent boundary layer flows. The program is capable of calculating the following types of flow: (1) incompressible or compressible, (2) two dimensional or axisymmetric, and (3) flows with significant transverse curvature. Also, the program can handle a large variety of boundary conditions, such as blowing or suction, arbitrary temperature distributions and arbitrary wall heat fluxes. The program has been specialized to the calculation of equilibrium air flows and all of the thermodynamic and transport properties used are for air. For the turbulent transport properties, the eddy viscosity approach has been used. Although the eddy viscosity models are semi-empirical, the model employed in the program has corrections for pressure gradients, suction and blowing and compressibility. The basic method of approach is to put the equations of motion into a finite difference form and then solve them by use of a digital computer. The program is written in FORTRAN 4 and requires small amounts of computer time on most scientific machines. For example, most laminar flows can be calculated in less than one minute of machine time, while turbulent flows usually require three or four minutes.
A computer program that uses a three-dimensional geometric technique for fitting a smooth surface to the component parts of an aircraft configuration is presented. The resulting surface equations are useful in performing various kinds of calculations in which a three-dimensional mathematical description is necessary. Programs options may be used to compute information for three-view and orthographic projections of the configuration as well as cross-section plots at any orientation through the configuration. The aircraft geometry input section of the program may be easily replaced with a surface point description in a different form so that the program could be of use for any three-dimensional surface equations.
The NASA/MSFC multilayer diffusion models are discribed which are used in applying meteorological information to the estimation of toxic fuel hazards resulting from the launch of rocket vehicle and from accidental cold spills and leaks of toxic fuels. Background information, definitions of terms, description of the multilayer concept are presented along with formulas for determining the buoyant rise of hot exhaust clouds or plumes from conflagrations, and descriptions of the multilayer diffusion models. A brief description of the computer program is given, and sample problems and their solutions are included. Derivations of the cloud rise formulas, users instructions, and computer program output lists are also included.
Spherical roller bearings have typically been used in applications with speeds limited to about 5000 rpm and loads limited for operation at less than about 0.25 million DN. However, spherical roller bearings are now being designed for high load and high speed applications including aerospace applications. A computer program, SASHBEAN, was developed to provide an analytical tool to design, analyze, and predict the performance of high speed, single row, angular contact (including zero contact angle), spherical roller bearings. The material presented is the mathematical formulation and analytical methods used to develop computer program SASHBEAN. For a given set of operating conditions, the program calculates the bearings ring deflections (axial and radial), roller deflections, contact areas stresses, depth and magnitude of maximum shear stresses, axial thrust, rolling element and cage rotational speeds, lubrication parameters, fatigue lives, and rates of heat generation. Centrifugal forces and gyroscopic moments are fully considered. The program is also capable of performing steady-state and time-transient thermal analyses of the bearing system.
A viscous-inviscid interaction method to calculate the subsonic and transonic flow over nozzle afterbodies with supersonic jet exhausts was developed. The method iteratively combines a relaxation solution of the full potential equation for the inviscid external flow, a shock capturing-shock fitting inviscid jet solution, an integral boundary layer solution, a control volume method for treating separated flows, and an overlaid mixing layer solution. A computer program called RAXJET which incorporates the method, illustrates the predictive capabilities of the method by comparison with experimental data is described, a user's guide to the computer program is provided. The method accurately predicts afterbody pressures, drag, and flow field properties for attached and separated flows for which no shock induced separation occurs.
The Plume Impingement Program (PLIMP) is a computer code used to predict impact pressures, forces, moments, heating rates, and contamination on surfaces due to direct impingement flowfields. Typically, it has been used to analyze the effects of rocket exhaust plumes on nearby structures from ground level to the vacuum of space. The program normally uses flowfields generated by the MOC, RAMP2, SPF/2, or SFPGEN computer programs. It is capable of analyzing gaseous and gas/particle flows. A number of simple subshapes are available to model the surfaces of any structure. The original PLIMP program has been modified many times of the last 20 years. The theoretical bases for the referenced major changes, and additional undocumented changes and enhancements since 1988 are summarized in volume 1 of this report. This volume is the User's Input Guide and should be substituted for all previous guides when running the latest version of the program. This version can operate on VAX and UNIX machines with NCAR graphics ability.
Numerical methods and computer programs for two-dimensional and axisymmetric supersonic flow over blunt-nosed and flared bodies
This report presents a description of a computer program mechanized to perform the Paull and Unger process of simplifying incompletely specified sequential machines. An understanding of the process, as given in Ref. 3, is a prerequisite to the use of the techniques presented in this report. This process has specific application in the design of asynchronous digital machines and was used in the design of operational support equipment for the Mariner 1966 central computer and sequencer. A typical sequential machine design problem is presented to show where the Paull and Unger process has application. A description of the Paull and Unger process together with a description of the computer algorithms used to develop the program mechanization are presented. Several examples are used to clarify the Paull and Unger process and the computer algorithms. Program flow diagrams, program listings, and a program user operating procedures are included as appendixes.
A detailed description of the Langley computer program PLOTWD which plots and fairs experimental wind-tunnel data is presented. The program was written for use primarily on the Langley CDC computer and CALCOMP plotters. The fundamental operating features of the program are that the input data are read and written to a random-access file for use during program execution, that the data for a selected run can be sorted and edited to delete duplicate points, and that the data can be plotted and faired using tension splines, least-squares polynomial, or least-squares cubic-spline curves. The most noteworthy feature of the program is the simplicity of the user-supplied input requirements. Several subroutines are also included that can be used to draw grid lines, zero lines, axis scale values and lables, and legends. A detailed description of the program operational features and each sub-program are presented. The general application of the program is also discussed together with the input and output for two typical plot types. A listing of the program code, user-guide, and output description are presented in appendices. The program has been in use at Langley for several years and has proven to be both easy to use and versatile.
Future Orbit Transfer Vehicles (OTV) presently under consideration require rocket engines delivering a high specific impulse. This high performance can be obtained with large area ratio thrust chambers using oxygen with either hydrogen or hydrocarbon fuels. In the projected nozzles, the combustion products are expanded to low pressure and temperature levels at high Mach number, a domain which has not been experienced with existing rocket engines. In order to assist in the design and evaluation of OTV engines, modifications have recently been incorporated in the Two-Dimensional-Kinetics (TDK) computer program. These modifications include the treatment of shock waves induced by the wall curvature, and shocks which are attached to the wall. The reflection of a shock wave at the flow axis is also treated in an appropriate manner and multiple shock reflections are allowed. A Boundary Layer Module (BLM) has also been incorporated into the TDK computer program in order to provide an automated treatment of the rigorous JANNAF rocket engine performance prediction procedure.
Computer program to solve heat conduction equation in lunar surface for temperature-dependent thermal properties
This report describes HYDES, a hybrid computer program capable of simulating one-spool turbojet, two-spool turbojet, or two-spool turbofan engine dynamics. HYDES is also capable of simulating two- or three-stream turbofans with or without mixing of the exhaust streams. The program is intended to reduce the time required for implementing dynamic engine simulations. HYDES was developed for running on the Lewis Research Center's Electronic Associates (EAI) 690 Hybrid Computing System and satisfies the 16384-word core-size and hybrid-interface limits of that machine. The program could be modified for running on other computing systems. The use of HYDES to simulate a single-spool turbojet and a two-spool, two-stream turbofan engine is demonstrated. The form of the required input data is shown and samples of output listings (teletype) and transient plots (x-y plotter) are provided. HYDES is shown to be capable of performing both steady-state design and off-design analyses and transient analyses.
The governing differential equations of motion for a high speed cylindrical roller bearing are programmed for numerical solution and plotted output. This computer program has the capability of performing a two dimensional or three dimensional simulation. Two numerical solutions of the governing differential equations were obtained to simulate the motion of a roller bearing, one for the two dimensional equations of motion and one for the three dimensional equations of motion. Computer generated plots were obtained and present such data as roller/cage interaction forces, roller/race traction forces, roller/race relative slip velocities and cage angular speed over a nondimensional time equivalent to 1.2 revolutions of the inner race. Roller axial displacement, roller skew angle, and skew moment are also plotted for the three dimensional solution. The trajectory of the cage center is plotted for both the two dimensional and three dimensional solutions.