Search NASA⌕ Search

SEARCH · Search NASA

Results for “Computer code”

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,261 records · Page 70

The 1991 version of the plume impingement computer program. Volume 2: User's input guide

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.

Bender, Robert L.↗

A Test of Thick-Target Nonuniform Ionization as an Explanation for Breaks in Solar Flare Hard X-Ray Spectra

Solar nonthermal hard X-ray (HXR) flare spectra often cannot be fitted by a single power law, but rather require a downward break in the photon spectrum. A possible explanation for this spectral break is nonuniform ionization in the emission region. We have developed a computer code to calculate the photon spectrum from electrons with a power-law distribution injected into a thick-target in which the ionization decreases linearly from 100% to zero. We use the bremsstrahlung cross-section from Haug (1997), which closely approximates the full relativistic Bethe-Heitler cross-section, and compare photon spectra computed from this model with those obtained by Kontar, Brown and McArthur (2002), who used a step-function ionization model and the Kramers approximation to the cross-section. We find that for HXR spectra from a target with nonuniform ionization, the difference (Delta-gamma) between the power-law indexes above and below the break has an upper limit between approx.0.2 and 0.7 that depends on the power-law index delta of the injected electron distribution. A broken power-law spectrum with a. higher value of Delta-gamma cannot result from nonuniform ionization alone. The model is applied to spectra obtained around the peak times of 20 flares observed by the Ramaty High Energy Solar Spectroscopic Imager (RHESSI from 2002 to 2004 to determine whether thick-target nonuniform ionization can explain the measured spectral breaks. A Monte Carlo method is used to determine the uncertainties of the best-fit parameters, especially on Delta-gamma. We find that 15 of the 20 flare spectra require a downward spectral break and that at least 6 of these could not be explained by nonuniform ionization alone because they had values of Delta-gamma with less than a 2.5% probability of being consistent with the computed upper limits from the model. The remaining 9 flare spectra, based on this criterion, are consistent with the nonuniform ionization model.

Holman, gordon↗

Rotating Desk for Collaboration by Two Computer Programmers

A special-purpose desk has been designed to facilitate collaboration by two computer programmers sharing one desktop computer or computer terminal. The impetus for the design is a trend toward what is known in the software industry as extreme programming an approach intended to ensure high quality without sacrificing the quantity of computer code produced. Programmers working in pairs is a major feature of extreme programming. The present desk design minimizes the stress of the collaborative work environment. It supports both quality and work flow by making it unnecessary for programmers to get in each other s way. The desk (see figure) includes a rotating platform that supports a computer video monitor, keyboard, and mouse. The desk enables one programmer to work on the keyboard for any amount of time and then the other programmer to take over without breaking the train of thought. The rotating platform is supported by a turntable bearing that, in turn, is supported by a weighted base. The platform contains weights to improve its balance. The base includes a stand for a computer, and is shaped and dimensioned to provide adequate foot clearance for both users. The platform includes an adjustable stand for the monitor, a surface for the keyboard and mouse, and spaces for work papers, drinks, and snacks. The heights of the monitor, keyboard, and mouse are set to minimize stress. The platform can be rotated through an angle of 40 to give either user a straight-on view of the monitor and full access to the keyboard and mouse. Magnetic latches keep the platform preferentially at either of the two extremes of rotation. To switch between users, one simply grabs the edge of the platform and pulls it around. The magnetic latch is easily released, allowing the platform to rotate freely to the position of the other user

Riley, John Thomas↗

Numerical simulation of Jet-A combustion approximated by improved propane chemical kinetics

Through the effort devoted to the chemical kinetics for propane air combustion, three mechanisms are developed. The full mechanism consists of 131 reactions. This mechanism is used as a guide for the evaluation of other mechanisms, but because of the long expected cpu time, it is not to be incorporated into the computer code KIVA-II for actual simulation. Through the sensitivity analysis, a reduced mechanism of 45 reactions is produced. But the calculated results from the 45 reaction mechanism are always low in temperature. Some efforts are devoted to correct this situation and details are included in this report. A simplified mechanism of reactions is successfully improved and computed results are compared with experimental data. Contour plots of physical parameters and species concentrations and results for emission indices of CO and NOx are presented.

Ying, Shuh-Jing↗

Overview of aerothermodynamic loads definition study

The objective of the Aerothermodynamic Loads Definition Study is to develop methods of accurately predicting the operating environment in advanced Earth-to-Orbit (ETO) propulsion systems, such as the Space Shuttle Main Engine (SSME) powerhead. Development of time averaged and time dependent three dimensional viscous computer codes as well as experimental verification and engine diagnostic testing are considered to be essential in achieving that objective. Time-averaged, nonsteady, and transient operating loads must all be well defined in order to accurately predict powerhead life. Described here is work in unsteady heat flow analysis, improved modeling of preburner flow, turbulence modeling for turbomachinery, computation of three dimensional flow with heat transfer, and unsteady viscous multi-blade row turbine analysis.

Gaugler, Raymond E.↗

On the generation of secondary motion in circular to rectangular transition ducts

A computer code has been developed for the calculation of three-dimensional turbulent flow in circular to rectangular transition ducts. The governing equations are written in a general curvilinear coordinate system so that all singly-connected geometrical configurations can be accommodated. Turbulence closure is through the standard k-epsilon model or a Reynolds stress closure model. Computed results show qualitative agreement with experimental data. Pressure-induced secondary flow is the dominant mode through the transition, but turbulence-induced secondary flow becomes prominent some distance downstream of the transition.

Demuren, A. O.↗

Three-dimensional modeling of diesel engine intake flow, combustion and emissions-2

A three-dimensional computer code, KIVA, is being modified to include state-of-the-art submodels for diesel engine flow and combustion. Improved and/or new submodels which have already been implemented and previously reported are: wall heat transfer with unsteadiness and compressibility, laminar-turbulent characteristic time combustion with unburned HC and Zeldo'vich NO(x), and spray/wall impingement with rebounding and sliding drops. Progress on the implementation of improved spray drop drag and drop breakup models, the formulation and testing of a multistep kinetics ignition model, and preliminary soot modeling results are described. In addition, the use of a block structured version of KIVA to model the intake flow process is described. A grid generation scheme was developed for modeling realistic (complex) engine geometries, and computations were made of intake flow in the ports and combustion chamber of a two-intake-value engine. The research also involves the use of the code to assess the effects of subprocesses on diesel engine performance. The accuracy of the predictions is being tested by comparisons with engine experiments. To date, comparisons were made with measured engine cylinder pressure, temperature and heat flux data, and the model results are in good agreement with the experiments. Work is in progress that will allow validation of in-cylinder flow and soot formation predictions. An engine test facility is described that is being used to provide the needed validation data. Test results were obtained showing the effect of injection rate and split injections on engine performance and emissions.

Reitz, R. D.↗

HITCAN for actively cooled hot-composite thermostructural analysis

A computer code, high temperature composite analyzer (HITCAN), was developed to analyze/design hot metal matrix composite structures. HITCAN is a general purpose code for predicting the global structural and local stress-strain response of multilayered (arbitrarily oriented) metal matrix structures both at the constituent (fiber, matrix, and interphase) and the structural level, including the fabrication process effects. The thermomechanical properties of the constituents are considered to be nonlinearly dependent on several parameters, including temperature, stress, and stress rate. The computational procedure employs an incremental iterative nonlinear approach utilizing a multifactor-interaction material behavior model, i.e., the material properties are expressed in terms of a product of several factors that affect the properties. HITCAN structural analysis capabilities (static, load stepping - a multistep static analysis with material properties updated at each step, modal, and buckling) for cooled hot structures are demonstrated through a specific example problem.

Chamis, C. C.↗

A numerical method for prediction of compressible turbulent flows with closure models

A new computer code to solve the time averaged Navier-Stokes equations is developed. Many of the state-of-the-art numerical techniques and algorithms have been tested and implemented in the program in order to achieve a better numerical accuracy and code efficiency. Various turbulence models are tested for a wide range of flows. The initial focus has been on two-equation eddy-viscosity models, which are the most advanced available in current compressible flow codes. The long term goal will be to test Reynolds-Stress models and to explore their performance in the high Mach number range. Although testing and improvement of turbulence models for supersonic and hypersonic flows is the primary objective of this research, part of the effort has been devoted to analyzing the vortex breakdown phenomena using new computer programs. Some preliminary results on the breakdown of a vortex flow in a tube are reported. Present calculations are restricted to two dimensional flow geometry.

Huang, P. G.↗

HITCAN for actively cooled hot-composite thermostructural analysis

A computer code, high temperature composite analyzer (HITCAN), was developed to analyze/design hot metal matrix composite structures. HITCAN is a general purpose code for predicting the global structural and local stress-strain response of multilayered (arbitrarily oriented) metal matrix structures both at the constituent (fiber, matrix, and interphase) and the structural level, including the fabrication process effects. The thermomechanical properties of the constituents are considered to be nonlinearly dependent on several parameters, including temperature, stress, and stress rate. The computational procedure employs an incremental iterative nonlinear approach utilizing a multifactor-interaction material behavior model, i.e., the material properties are expressed in terms of a product of several factors that affect the properties. HITCAN structural analysis capabilities (static, load stepping - a multistep static analysis with material properties updated at each step, modal, and buckling) for cooled hot structures are demonstrated through a specific example problem.

Chamis, C. C.↗

Eulerian Applications Project - xRage: High Explosives Package [Slides]

At the most elementary level, the high explosive model in xRage (or many other hydro codes) computes the amount of material that changes from a high explosive reactant to a high explosive product during a time step and accounts for the energy added to the flow by this reaction. The actual chemical reactions that occur during this conversion are not modeled. These processes are very complicated and involve levels of spatial and temporal resolution that are generally not practical for an application code. They would also require a state variable representation that includes all of the chemical species that are produced by the reaction and account for multiple stage chemical reactions. Instead the explosive models we generally use are phenomenological and try to represent the most important physics that drive the explosive wave in a very simplified way. These models generally have a number of tunable parameters that must be calibrated to experiment.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

Validation of RELAP5-3D He-FUS3 Data

The RELAP5-3D computer code was validated using data from the HE FUS3 experimental facility. HE FUS3 is a helium-cooled, electrically-heated experimental facility that was used to generate experimental data for a benchmark exercise aimed at the validation of system transient analysis codes for very-high-temperature gas reactor applications. RELAP5-3D was validated using results from seven steady-state tests and one transient that simulated a loss-of-flow accident.

97 MATHEMATICS AND COMPUTING↗

Predicting multi-wall structural response to hypervelocity impact using the hull code

Previously, multi-wall structures have been analyzed extensively, primarily through experiment, as a means of increasing the meteoroid/space debris impact protection of spacecraft. As structural configurations become more varied, the number of tests required to characterize their response increases dramatically. As an alternative to experimental testing, numerical modeling of high-speed impact phenomena is often being used to predict the response of a variety of structural systems under different impact loading conditions. The results of comparing experimental tests to Hull Hydrodynamic Computer Code predictions are reported. Also, the results of a numerical parametric study of multi-wall structural response to hypervelocity cylindrical projectile impact are presented.

Schonberg, William P.↗

SAM Code Validation using the Compact Integral Effects Test (CIET) Experimental Data

The primary objective of this study is to validate the system analysis code, SAM, using experimental data from the Compact Integral Effects Test (CIET) facility. The System Analysis Module (SAM) is an advanced and modern system analysis tool being developed at Argonne National Laboratory under the U.S. DOE Office of Nuclear Energy’s Nuclear Energy Advanced Modeling and Simulation (NEAMS) program. Originally aiming at supporting sodium fast reactor (SFR) design, it is now extended to cover other non-light water reactor (non-LWR) concepts, such as high-temperature gas-cooled reactor (HTGR) designs, molten salt-fueled, and molten salt-cooled reactor designs. To support SAM code development for the wide range of non-LWR applications, it is of paramount importance to validate the code against experiments relevant to these reactor concepts. The Compact Integral Effects Test (CIET) facility, which was designed to reproduce the thermal-hydraulics response of fluoride salt-cooled high-temperature reactors (FHR) under both forced- and natural-circulation conditions, has been identified and selected as one of benchmark test facilities for SAM code validation. Among available CIET experiments, three sets of them with distinctive characteristics were selected for SAM code validation purpose. The three sets experiments are: 1) Power step change transient tests; 2) DHX-DRACS natural circulation tests; and 3) heater frequency response tests. For all tests, steady-state or transient, SAM-predicted results show very good agreement with experimental data. The successful validation of SAM against these selected CIET experiments demonstrates that the computer code is well suited for thermal-hydraulics analysis of FHR designs.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Probabilistic Structural Analysis Methods (PSAM) for select space propulsion system components

The fourth year of technical developments on the Numerical Evaluation of Stochastic Structures Under Stress (NESSUS) system for Probabilistic Structural Analysis Methods is summarized. The effort focused on the continued expansion of the Probabilistic Finite Element Method (PFEM) code, the implementation of the Probabilistic Boundary Element Method (PBEM), and the implementation of the Probabilistic Approximate Methods (PAppM) code. The principal focus for the PFEM code is the addition of a multilevel structural dynamics capability. The strategy includes probabilistic loads, treatment of material, geometry uncertainty, and full probabilistic variables. Enhancements are included for the Fast Probability Integration (FPI) algorithms and the addition of Monte Carlo simulation as an alternate. Work on the expert system and boundary element developments continues. The enhanced capability in the computer codes is validated by applications to a turbine blade and to an oxidizer duct.

Source record↗

Binary Black Holes, Gravitational Waves, and Numerical Relativity

The final merger of two black holes releases a tremendous amount of energy and is one of the brightest sources in the gravitational wave sky. Observing these sources with gravitational wave detectors requires that we know the radiation waveforms they emit. Since these mergers take place in regions of very strong gravitational fields, we need to solve Einstein's equations of general relativity on a computer in order to calculate these waveforms. For more than 30 years, scientists have tried to compute these waveforms using the methods of numerical relativity. The resulting computer codes have been plagued by instabilities, causing them to crash well before the black holes in the binary could complete even a single orbit. Recently this situation has changed dramatically, with a series of amazing breakthroughs. This talk will take you on this quest for the holy grail of numerical relativity, showing how a spacetime is constructed on a computer to build a simutation laboratory for binary black hole mergers. We will focus on the recent advances that are revealing these waveforms, and the dramatic new potential for discoveries that arises when these sources will be observed by LIGO and LISA.

Centrella, Joan↗

Binary Black Holes, Gravitational Waves, and Numerical Relativity

The final merger of two black holes releases a tremendous amount of energy and is one of the brightest sources in the gravitational wave sky. Observing these sources with gravitational wave detectors requires that we know the radiation waveforms they emit. Since these mergers take place in regions of very strong gravitational fields, we need to solve Einstein's equations of general relativity on a computer in order to calculate these waveforms. For more than 30 years, scientists have tried to compute these waveforms using the methods of numerical relativity. The resulting computer codes have been plagued by instabilities. causing them to crash well before the black hole:, in the binary could complete even a single orbit. Recently this situation has changed dramatically, with a series of amazing breakthroughs. This talk will take you on this quest for the holy grail of numerical relativity, showing how a spacetime is constructed on a computer to build a simulation laboratory for binary black hole mergers. We will focus on the recent advances that are revealing these waveforms, and the dramatic new potential for discoveries that arises when these sources will be observed by LIGO and LISA.

Centrella, Joan↗

Binary Black Holes, Gravitational Waves, and Numerical Relativity

The final merger of two black holes releases a tremendous amount of energy and is one of the brightest sources in the gravitational wave sky. Observing these sources with gravitational wave detectors requires that we know the radiation waveforms they emit. Since these mergers take place in regions of very strong gravitational fields, we need to solve Einstein's equations of general relativity on a computer in order to calculate these waveforms. For more than 30 years, scientists have tried to compute these waveforms using the methods of numerical relativity. The resulting computer codes have been plagued by instabilities, causing them to crash well before the black holes in the binary could complete even a single orbit. Recently this situation has changed dramatically, with a series of amazing breakthroughs. This talk will take you on this quest for the holy grail of numerical relativity, showing how a spacetime is constructed on a computer to build a simulation laboratory for binary black hole mergers. We will focus on the recent advances that are revealing these waveforms, and the dramatic new potential for discoveries that arises when these sources will be observed by LIGO and LISA.

Centrella, Joan↗