Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “computer Codes Computer Programs”

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.

At least 1,009 records · Page 56

Modeling of Cryogenic Heated-Tube Flow Boiling Experiments of Nitrogen and Methane with the Generalized Fluid System Simulation Program

Accurate modeling of cryogenic boiling heat transfer is vital for the development of extended-duration space missions. Such missions may require the transfer of cryogenic propellants from in-space storage depots or the cooling of nuclear reactors. Purdue University in collaboration with NASA has assembled a database of cryogenic flow boiling data points from heated-tube experiments dating back to 1959, which has been used to develop new flow boiling correlations specifically for cryogens. Computational models of several of these experiments have been constructed in the Generalized Fluid System Simulation Program (GFSSP), a network flow code developed at NASA’s Marshall Space Flight Center. The new Purdue-developed correlations cover the full boiling curve: onset of nucleate boiling, nucleate boiling, critical heat flux, and film boiling. These correlations have been coded into GFSSP user subroutines. The fluids modeled are methane and nitrogen. Predictions of wall temperature and heat transfer coefficient are presented and compared to the test data.

Multi-Phase Flow↗

Modeling of Cryogenic Heated-Tube Flow Boiling Experiments of Hydrogen and Helium With the Generalized Fluid System Simulation Program

Accurate modeling of cryogenic boiling heat transfer is vital for the development of extended-duration space missions. Such missions may require the transfer of cryogenic propellants from in-space storage depots or the cooling of nuclear reactors. Purdue University in collaboration with NASA has assembled a database of cryogenic flow boiling data points from heated-tube experiments dating back to 1959, which has been used to develop new flow boiling correlations specifically for cryogens. Computational models of several of these experiments have been constructed in the Generalized Fluid System Simulation Program (GFSSP), a network flow code developed at NASA’s Marshall Space Flight Center. The new Purdue-developed universal correlations cover the full boiling curve: onset of nucleate boiling, nucleate boiling, critical heat flux, and film boiling. These correlations have been coded into GFSSP user subroutines. The fluids modeled in this study are liquid hydrogen and liquid helium. Predictions of local wall temperature and pressure drop are presented and compared to the test data.

multi-phase flow↗

Users manual for program ADMIT: Admittance and pressure transfer function developed for use on a PC computer

The piping in a liquid rocket can assume complex configurations due to multiple tanks, multiple engines, and structures that must be piped around. The capability to handle some of these complex configurations have been incorporated into the ADMIT code. The capability to modify the input on line has been implemented. The configurations allowed include multiple tanks, multiple engines, the splitting of a pipe into unequal segments going to different (or the same) engines. This program will handle the following type elements: straight pipes, bends, inline accumulators, tuned stub accumulators, Helmholtz resonators, parallel resonators, pumps, split pipes, multiple tanks, and multiple engines.

Armstrong, Wilbur C.↗

As-Flown Shuttle Micro-Meteoroid Orbital Debris (MMOD) Assessment History

The final flight of the Space Shuttle program was completed on July 21st of 2011. The reusable nature of the Orbiter vehicles has provided NASA with a unique opportunity to inspect and sample spacecraft surfaces exposed to the MMOD environment in low Earth orbit. This paper will provide details of MMOD damages found on Orbiter surfaces after recent flights, as well as putting these damages in context by providing historical comparisons of recent damage to previous flights. The Bumper threat assessment computer code is used by NASA to determine spacecraft MMOD penetration risk. The tool was used before each mission to estimate pre-flight MMOD risk for the Shuttle Program Office using planned vehicles attitudes and exposure times. This paper summarizes the efforts of the authors to assess the expected number MMOD impacts of selected particles sizes of 52 shuttle missions ranging from STS-50 (1992) through STS-133 (2011) using as-flown attitude data. The missions encompass the Shuttle/Mir precursors to the International Space Station (ISS) as well as the ISS assembly flights. The paper includes a comparison of observed MMOD impact damage to the predicted number of impacts from Bumper code assessments. This data can be used to validate engineering models of the orbital debris and micrometeoroid environment, such NASA s ORDEM debris model and MEM meteoroid model.

Hyde, James L.↗

Prediction of sound radiated from different practical jet engine inlets

Existing computer codes for calculating the far field radiation patterns surrounding various practical jet engine inlet configurations under different excitation conditions were upgraded. The computer codes were refined and expanded so that they are now more efficient computationally by a factor of about three and they are now capable of producing accurate results up to nondimensional wave numbers of twenty. Computer programs were also developed to help generate accurate geometrical representations of the inlets to be investigated. This data is required as input for the computer programs which calculate the sound fields. This new geometry generating computer program considerably reduces the time required to generate the input data which was one of the most time consuming steps in the process. The results of sample runs using the NASA-Lewis QCSEE inlet are presented and comparison of run times and accuracy are made between the old and upgraded computer codes. The overall accuracy of the computations is determined by comparison of the results of the computations with simple source solutions.

Zinn, B. T.↗

Engine structures modeling software system: Computer code. User's manual

ESMOSS is a specialized software system for the construction of geometric descriptive and discrete analytical models of engine parts, components and substructures which can be transferred to finite element analysis programs such as NASTRAN. The software architecture of ESMOSS is designed in modular form with a central executive module through which the user controls and directs the development of the analytical model. Modules consist of a geometric shape generator, a library of discretization procedures, interfacing modules to join both geometric and discrete models, a deck generator to produce input for NASTRAN and a 'recipe' processor which generates geometric models from parametric definitions. ESMOSS can be executed both in interactive and batch modes. Interactive mode is considered to be the default mode and that mode will be assumed in the discussion in this document unless stated otherwise.

Source record↗

Rapid Calculation of Spacecraft Trajectories Using Efficient Taylor Series Integration

A variable-order, variable-step Taylor series integration algorithm was implemented in NASA Glenn's SNAP (Spacecraft N-body Analysis Program) code. SNAP is a high-fidelity trajectory propagation program that can propagate the trajectory of a spacecraft about virtually any body in the solar system. The Taylor series algorithm's very high order accuracy and excellent stability properties lead to large reductions in computer time relative to the code's existing 8th order Runge-Kutta scheme. Head-to-head comparison on near-Earth, lunar, Mars, and Europa missions showed that Taylor series integration is 15.8 times faster than Runge- Kutta on average, and is more accurate. These speedups were obtained for calculations involving central body, other body, thrust, and drag forces. Similar speedups have been obtained for calculations that include J2 spherical harmonic for central body gravitation. The algorithm includes a step size selection method that directly calculates the step size and never requires a repeat step. High-order Taylor series integration algorithms have been shown to provide major reductions in computer time over conventional integration methods in numerous scientific applications. The objective here was to directly implement Taylor series integration in an existing trajectory analysis code and demonstrate that large reductions in computer time (order of magnitude) could be achieved while simultaneously maintaining high accuracy. This software greatly accelerates the calculation of spacecraft trajectories. At each time level, the spacecraft position, velocity, and mass are expanded in a high-order Taylor series whose coefficients are obtained through efficient differentiation arithmetic. This makes it possible to take very large time steps at minimal cost, resulting in large savings in computer time. The Taylor series algorithm is implemented primarily through three subroutines: (1) a driver routine that automatically introduces auxiliary variables and sets up initial conditions and integrates; (2) a routine that calculates system reduced derivatives using recurrence relations for quotients and products; and (3) a routine that determines the step size and sums the series. The order of accuracy used in a trajectory calculation is arbitrary and can be set by the user. The algorithm directly calculates the motion of other planetary bodies and does not require ephemeris files (except to start the calculation). The code also runs with Taylor series and Runge-Kutta used interchangeably for different phases of a mission.

Scott, James R.↗

Wind-tunnel preflight test program for aeroassist flight experiment

A proposed Aeroassist Flight Experiment (AFE) will simulate return from geosynchronous orbit to low-earth orbit to provide fundamental information about the hypersonic flow conditions surrounding aeroassisted orbital transfer vehicles (AOTV's). Future AOTV design must rely heavily on computational fluid dynamic (CFD) computer codes now under development, since the entry conditions are beyond the simulation capabilities of ground-based test facilities. However, existing codes and best available ground-based test facilities must provide for design of the AFE. The AFE ground-based test program provides data for calibration of existing and developing pertinent CFD codes. Confidence in code prediction capability increases with ability to predict forces and moments, pressure and heat-transfer distributions, shock shapes, and surface streamline directions. The ground-based test program provides these data over a wide range of hypersonic test parameters. This paper describes the test program, models, facilities, and some representative results.

Wells, William L.↗

Demonstration of the range over which the Langley Research Center digital computer charring ablation program (CHAP) can be used with confidence: Comparisons of CHAP predictions and test data for three ablation materials

Comparisons of ablation calculations with the charring ablation computer code and ablation test data are presented over a wide range of environmental conditions in air for three materials: low-density nylon phenolic, Avcoat 5026-39HC/G, and a filled silicon elastomer. Heat fluxes considered range from over 500 Btu/sq ft-sec to less than 50 Btu/sq ft-sec. Pressures range from 0.5 atm to .004 atm. Enthalpies range from about 2000 Btu/lb to 18000 Btu/lb. Predictions of recession, pyrolysis penetration, and thermocouple responses are considered. Recession predictions for nylon phenolic are good as steady state is approached, but strongly transient cases are underpredicted. Pyrolysis penetrations and thermocouple responses are very well predicted. Recession amounts for Avcoat and silicone elastomer are less well predicted, although high heat flux cases near steady state are fairly satisfactory. Pyrolysis penetrations and thermocouple responses are very well predicted.

Moyer, C. B.↗

Computer simulation of dust grain evolution

The latest results are reported from a Monte Carlo code that is being developed at NASA Ames. The goal of this program, is to derive from the observed and presumed properties of the interstellar medium (ISM) the following information: (1) the size spectrum of interstellar dust; (2) the chemical structure of interstellar dust; (3) interstellar abundances; and (4) the lifetime of a dust grain in the ISM. Presently this study is restricted to refractory interstellar material, i.e., the formation and destruction of ices are not included in the program. The program is embedded in an analytic solution for the bulk chemical evolution of a two-phase interstellar medium in which stars are born in molecular clouds, but new nucleosynthesis products and stellar return are entered into a complementary intercloud medium. The well-mixed matter of each interstellar phase is repeatedly cycled stochastically through the complementary phase and back. Refractory dust is created by thermal condensation as stellar matter flows away from sites of nucleosynthesis such as novae and supernovae and/or from the matter returned from evolved intermediate stars. The history of each particle is traced by standard Monte Carlo techniques as it is sputtered and fragmented by supernova shock waves in the intercloud medium. It also accretes an amorphous mantle of gaseous refractory atoms when its local medium joins with the molecular cloud medium. Finally it encounters the possibility of astration (destruction by star formation) within the molecular clouds.

Liffman, K.↗

Rapid Aeroelastic Analysis of Blade Flutter in Turbomachines

The LINFLUX-AE computer code predicts flutter and forced responses of blades and vanes in turbomachines under subsonic, transonic, and supersonic flow conditions. The code solves the Euler equations of unsteady flow in a blade passage under the assumption that the blades vibrate harmonically at small amplitudes. The steady-state nonlinear Euler equations are solved by a separate program, then equations for unsteady flow components are obtained through linearization around the steady-state solution. A structural-dynamics analysis (see figure) is performed to determine the frequencies and mode shapes of blade vibrations, a preprocessor interpolates mode shapes from the structural-dynamics mesh onto the LINFLUX computational-fluid-dynamics mesh, and an interface code is used to convert the steady-state flow solution to a form required by LINFLUX. Then LINFLUX solves the linearized equations in the frequency domain to calculate the unsteady aerodynamic pressure distribution for a given vibration mode, frequency, and interblade phase angle. A post-processor uses the unsteady pressures to calculate generalized aerodynamic forces, response amplitudes, and eigenvalues (which determine the flutter frequency and damping). In comparison with the TURBO-AE aeroelastic-analysis code, which solves the equations in the time domain, LINFLUX-AE is 6 to 7 times faster.

Trudell, J. J.↗

Parallel Computation of Unsteady Flows on a Network of Workstations

Parallel computation of unsteady flows requires significant computational resources. The utilization of a network of workstations seems an efficient solution to the problem where large problems can be treated at a reasonable cost. This approach requires the solution of several problems: 1) the partitioning and distribution of the problem over a network of workstation, 2) efficient communication tools, 3) managing the system efficiently for a given problem. Of course, there is the question of the efficiency of any given numerical algorithm to such a computing system. NPARC code was chosen as a sample for the application. For the explicit version of the NPARC code both two- and three-dimensional problems were studied. Again both steady and unsteady problems were investigated. The issues studied as a part of the research program were: 1) how to distribute the data between the workstations, 2) how to compute and how to communicate at each node efficiently, 3) how to balance the load distribution. In the following, a summary of these activities is presented. Details of the work have been presented and published as referenced.

Source record↗

A New Inversion Routine to Produce Vertical Electron-Density Profiles from Ionospheric Topside-Sounder Data

Two software applications have been produced specifically for the analysis of some million digital topside ionograms produced by a recent analog-to-digital conversion effort of selected analog telemetry tapes from the Alouette-2, ISIS-1 and ISIS-2 satellites. One, TOPIST (TOPside Ionogram Scalar with True-height algorithm) from the University of Massachusetts Lowell, is designed for the automatic identification of the topside-ionogram ionospheric-reflection traces and their inversion into vertical electron-density profiles Ne(h). TOPIST also has the capability of manual intervention. The other application, from the Goddard Space Flight Center based on the FORTRAN code of John E. Jackson from the 1960s, is designed as an IDL-based interactive program for the scaling of selected digital topside-sounder ionograms. The Jackson code has also been modified, with some effort, so as to run on modern computers. This modification was motivated by the need to scale selected ionograms from the millions of Alouette/ISIS topside-sounder ionograms that only exist on 35-mm film. During this modification, it became evident that it would be more efficient to design a new code, based on the capabilities of present-day computers, than to continue to modify the old code. Such a new code has been produced and here we will describe its capabilities and compare Ne(h) profiles produced from it with those produced by the Jackson code. The concept of the new code is to assume an initial Ne(h) and derive a final Ne(h) through an iteration process that makes the resulting apparent-height profile fir the scaled values within a certain error range. The new code can be used on the X-, O-, and Z-mode traces. It does not assume any predefined profile shape between two contiguous points, like the exponential rule used in Jackson s program. Instead, Monotone Piecewise Cubic Interpolation is applied in the global profile to keep the monotone nature of the profile, which also ensures better smoothness in the final profile than in Jackson s program. The new code uses the complete refractive index expression for a cold collisionless plasma and can accommodate the IGRF, T96, and other geomagnetic field models.

Wang, Yongli↗

Neutron ages computed from experimental activation data

Computer program reduces time manually required to compute neutron age and to provide definite plan of procedural choices. Program allows convenient comparison of several fitting and error analysis procedures. Program code provides for error estimation of various integration options.

Semler, T. T.↗

Optimal low thrust geocentric transfer

A computer code which will rapidly calculate time-optimal low thrust transfers is being developed as a mission analysis tool. The final program will apply to NEP or SEP missions and will include a variety of environmental effects. The current program assumes constant acceleration. The oblateness effect and shadowing may be included. Detailed state and costate equations are given for the thrust effect, oblateness effect, and shadowing. A simple but adequate model yields analytical formulas for power degradation due to the Van Allen radiation belts for SEP missions. The program avoids the classical singularities by the use of equinoctial orbital elements. Kryloff-Bogoliuboff averaging is used to facilitate rapid calculation. Results for selected cases using the current program are given.

Edelbaum, T. N.↗

Numerical solution of the Navier-Stokes equations for arbitrary blunt bodies in supersonic flows

A time-dependent, two-dimensional Navier-Stokes code employing the body-fitted coordinate technique has been developed for supersonic flows past blunt bodies of arbitrary shape. The computer program is based on the finite-difference approximation of the compressible Navier-Stokes equations transformed to nonorthogonal curvilinear coordinates with the contravariant components of the velocity vector as dependent variables. The bow shock ahead of the body is obtained as part of the solution, by 'shock capturing'. Numerical solutions of the complete equations are presented in detail for free-stream Mach number 4.6, Reynolds number 10,000, and an isothermal wall temperature of 556 K for a circular cylinder with the free-stream outer boundaries forming a hyperbola in the front and a circular arc in the back.

Warsi, Z. U. A.↗

Solar Power Satellite antenna phase control system hardware simulation, phase 4. Volume 3: SOLARSIM users manual

Documentation of the interactive software package designed to predict the effect of certain electrical and mechanical imperfections on the performance of the SPS system is presented. The capabilities of the SOLARSIM program to quantify the spacetenna performance parameter values are described. The SOLARSIM package can compute the RMS pointing error, title effects, MPTX code tracking loop performance, MPTX carrier tracking loop performance, averaged power pattern, and the power transfer efficiency.

Lindsey, W. C.↗