Search NASA⌕ Search

SEARCH · Search NASA

Results for “Computer Codes”

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,603 records · Page 89

OMV/VTE variable thrust engine analysis

The objective of the present work is to develop a predictive CFD based analytical tool for the variable thrust engine (VTE) in the orbital maneuvering vehicle (OMV). This objective is being accomplished within the framework of the Los Alamos KIVA computer code for chemically reactive flows with sprays. For the OMV application, the main structure of KIVA is to be retained while reformulating many of the phenomenological submodels, enhancing some of the numerics, and adding more features. The analytical model consists of the general conservation equations for two phase reactive flows and of submodels for turbulence, chemical reactions, and bipropellant sprays. Tailoring this model to the OMV engine brings about the added complexities of combustion and flow processes that occur in a liquid hypergolic propellant rocket chamber. This report exposes the foundation upon which the analytical tool is being constructed and developed. Results from a cursory computational exercise involving the simulation of the flow and combustion processes in a hypothetical N2H4/N2O4 rocket engine thrust chamber is presented and discussed.

Larosillere, Louis↗

Newton-Krylov Method in Applications of Solving Two-phase Problems using Drift Flux Model

In this paper, Jacobian-free Newton-Krylov (JFNK) method is investigated in an application to implicitly solve two-phase flow problems using four-equation drift flux model. The closure models include the EPRI drift flux correlations and additional constitutive models to determine flow regimes, wall boiling and interfacial heat/mass transfer, two-phase flow wall friction, etc. Different from many traditional computer codes, fully implicit methods are used for the time integration. The resulted nonlinear discretized equation system is solved using the JFNK method. Expensive and error prone derivation and implementation of the analytical Jacobian matrix are avoided. Numerical results are successfully validated using existing experimental data on flow boiling under forced convection conditions in both pipe and rod bundle geometries.

42 ENGINEERING↗

OMV/VTE variable thrust engine analysis

The objective of the present work is to develop a predictive CFD based analytical tool for the Variable Thrust Engine (VTE) in the Orbital Maneuvering Vehicle (OMV). This objective is being accomplished within the framework of the Los Alamos KIVA computer code for chemically reactive flows with sprays. For the OMV application, the main structure of KIVA is to be retained while reformulating many of the phenomenological submodels, enhancing some of the numerics, and adding more features. The analytical model consists of the general conservation equations for two-phase reactive flows and of submodels for turbulence, chemical reactions, and bipropellant sprays. Tailoring this model to the OMV engine brings about the added complexities of combustion and flow processes that occur in a liquid hypergolic propellant rocket chamber. This report exposes the foundation upon which the analytical tool is being constructed and developed. Results from a cursory computational exercise involving the simulation of the flow and combustion processes in a hypothetical N2H4/N204 rocket engine thrust chamber is presented and discussed.

Larosillere, L.↗

SCEPTRE 2.6 User's Guide

Sandia’s Computational Engine for Particle Transport for Radiation Effects (SCEPTRE) is a computer code that solves the linear Boltzmann transport equation, particularly targeting coupled photon-electron problems. It uses unstructured finite element meshes in space, multigroup in energy, and discrete ordinates (S n ) or other methods in angle. SCEPTRE uses an xml-based input file to specify the problem. This report documents the options and syntax of that input file.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

MatLib-1.1: Nuclear Material Properties Library: Developed under NQA-1-2017

The U.S. Nuclear Regulatory Commission (NRC) uses the computer code Fuel Analysis under Steady-state and Transients (FAST) to model steady-state and transient fuel behavior to support regulatory decisions. FAST relies on a material properties library (MatLib) that contains the thermal and mechanical properties of the nuclear materials and coolants of interest to support the US commercial nuclear industry. MatLib contains properties for a variety of nuclear fuels, cladding and other structural materials, gases, and coolants. In this document, material property correlations for the materials contained within MatLib are presented and discussed. When available, comparisons are made between the material property correlations and available data. Additionally, uncertainties are quantified on the material properties, which is then used by the NRC to support uncertainty quantification for best-estimate plus uncertainty safety evaluation reviews. This document describes MatLib-1.1, which was updated from MatLib-1.0 to include the addition of Lua functionality and code changes for integration with FAST. The results of MatLib did not change between MatLib-1.0 and MatLib-1.1; however, the source code changes warranted incrementing the minor numbering of the software version number. This document is one of a series of documents on FAST; the other documents detail the models used by FAST as well as its integral assessment to experiments and commercial data.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

235U Resolved Resonance Evaluation for Benchmark Calculations in the Intermediate Energy Region

The Working Party on International Nuclear Data Evaluation Co-operation (WPEC) subgroup 29 (SG 29) was established to investigate an issue with the 235U capture cross-section in the energy range 0.1 to 2.25 keV. The WPEC criticality calculation results indicated an overestimation of the 235U capture cross-section of 10% or more. To understand and solve the problem, a recommendation was made to perform new capture cross-section measurements followed by a resonance evaluation. Hence, time of-flight capture cross-section measurements were done at the Rensselaer Polytechnic Institute (RPI) and at the Los Alamos National Laboratory (LANL). These new measurements were used together with the computer code SAMMY to reevaluate the 235U resonance parameters in the energy range from thermal to 2.25 keV. The impact of the new evaluation in benchmark calculations was done for the critical benchmark sensitive to the 0.1 to 2.25 keV energy range. The purpose of this work is to describe the 235U SAMMY evaluation and present the critical benchmark results.

97 MATHEMATICS AND COMPUTING↗

DYCAST: A finite element program for the crash analysis of structures

DYCAST is a nonlinear structural dynamic finite element computer code developed for crash simulation. The element library contains stringers, beams, membrane skin triangles, plate bending triangles and spring elements. Changing stiffnesses in the structure are accounted for by plasticity and very large deflections. Material nonlinearities are accommodated by one of three options: elastic-perfectly plastic, elastic-linear hardening plastic, or elastic-nonlinear hardening plastic of the Ramberg-Osgood type. Geometric nonlinearities are handled in an updated Lagrangian formulation by reforming the structure into its deformed shape after small time increments while accumulating deformations, strains, and forces. The nonlinearities due to combined loadings are maintained, and stiffness variation due to structural failures are computed. Numerical time integrators available are fixed-step central difference, modified Adams, Newmark-beta, and Wilson-theta. The last three have a variable time step capability, which is controlled internally by a solution convergence error measure. Other features include: multiple time-load history tables to subject the structure to time dependent loading; gravity loading; initial pitch, roll, yaw, and translation of the structural model with respect to the global system; a bandwidth optimizer as a pre-processor; and deformed plots and graphics as post-processors.

Pifko, A. B.↗

An implementation of a chemical and thermal nonequilibrium flow solver on unstructured meshes and application to blunt bodies

This paper presents a nonequilibrium flow solver, implementation of the algorithm on unstructured meshes, and application to hypersonic flow past blunt bodies. Air is modeled as a mixture of five chemical species, namely O2, N2, O, NO, and N, having two temperatures namely translational and vibrational. The solution algorithm is a cell centered, point implicit upwind scheme that employs Roe's flux difference splitting technique. Implementation of this algorithm on unstructured meshes is described. The computer code is applied to solve Mach 15 flow with and without a Type IV shock interference on a cylindrical body of 2.5mm radius representing a cowl lip. Adaptively generated meshes are employed, and the meshes are refined several times until the solution exhibits detailed flow features and surface pressure and heat flux distributions. Effects of a catalytic wall on surface heat flux distribution are studied. For the Mach 15 Type IV shock interference flow, present results showed a peak heat flux of 544 MW/m2 for a fully catalytic wall and 431 MW/m(exp 2) for a noncatalytic wall. Some of the results are compared with available computational data.

Prabhu, Ramadas K.↗

DHC-6 Twin Otter Tailplane Airfoil Section Testing in the Ohio State University 7x10 Wind Tunnel

Ice contaminated tailplane stall (ICTS) has been found to be responsible for 16 accidents with 139 fatalities over the last three decades, and is suspected to have played a role in other accidents and incidents. The need for fundamental research in this area has been recognized at three international conferences sponsored by the FAA since 1991. In order to conduct such research, a joint NASA/FAA Tailplane Icing Program was formed in 1994: the Ohio State University has played an important role in this effort. The program employs icing tunnel testing, dry wind tunnel testing, flight testing, and analysis using a six-degrees-of-freedom computer code tailored to this problem. A central goal is to quantify the effect of tailplane icing on aircraft stability and control to aid in the analysis of flight test procedures to identify aircraft susceptibility to ICTS. This report contains the results ot testing of a full scale 2D model of a tailplane section of NASA's Icing Research Aircraft, with and without ice shapes, in an Ohio State University 7 x 10 Low Speed wind tunnel in 1994. The results have been integrated into a comprehensive database of aerodynamic coefficients and stability and control derivatives that will permit detailed analysis of flight test results with the analytical computer program. The testing encompassed a full range of angles of attack and elevator deflections, as well as two velocities to evaluate Reynolds number effects. Lift, drag, pitching moment, and hinge moment coefficients were obtained. In addition. instrumentation for use during flight testing was verified to be effective, all components showing acceptable fidelity. Comparison of clean and iced airfoil results show the ice shapes causing a significant decrease in the magnitude of CLmax (from -1.3 to -0.64) and associated stall angle (from -18.6 deg to -8.2 deg). Furthermore, the ice shapes caused an increase in hinge moment coefficient of approximately 0.02, the change being markedly abrupt for one of the ice shapes. A noticeable effect of elevator deflection is that magnitude of the stall angle is decreased for negative (upward) elevator deflections. All these result are consistent with observed tailplane phenomena. and constitute an effective set of data for comprehensive analysis of ICTS

Hiltner, Dale↗

A combined Eulerian-Lagrangian two-phase flow analysis of SSME HPOTP nozzle plug trajectories. I - Methodology

A computer code has been developed for the analysis of SSME (Space Shuttle Main Engine) HPOTP (High Pressure Oxidizer Turbo Pump) nozzle plug trajectories in the turnaround duct downstream of the turbine. The algorithm is based on a combined Eulerian-Lagrangian analysis originally developed for the study of two-phase flows. The Lagrangian part of this analysis has been enhanced to include three-dimensional particle motion and the effect of particle-wall collisions in complex geometries (with a large number of boundaries). The sensitivity of the nozzle plug trajectories to a variety of parameters has been determined, via the qualitative analysis of a select number of computed trajectories. The results of extensive parametric studies have been reported in a companion paper.

De Jong, Frederik J.↗

Computer access security code system

A security code system for controlling access to computer and computer-controlled entry situations comprises a plurality of subsets of alpha-numeric characters disposed in random order in matrices of at least two dimensions forming theoretical rectangles, cubes, etc., such that when access is desired, at least one pair of previously unused character subsets not found in the same row or column of the matrix is chosen at random and transmitted by the computer. The proper response to gain access is transmittal of subsets which complete the rectangle, and/or a parallelepiped whose opposite corners were defined by first groups of code. Once used, subsets are not used again to absolutely defeat unauthorized access by eavesdropping, and the like.

Collins, Earl R., Jr.↗

METCAN simulation of candidate metal matrix composites for high temperature applications

The METCAN (Metal Matrix Composite Analyzer) computer code is used to simulate the nonlinear behavior of select metal matrix composites in order to assess their potential for high temperature structural applications. Material properties for seven composites are generated at a fiber volume ratio of 0.33 for two bonding conditions (a perfect bond and a weak interphase case) at various temperatures. A comparison of the two bonding conditions studied shows a general reduction in value of all properties (except CTE) for the weak interphase case from the perfect bond case. However, in the weak interphase case, the residual stresses that develop are considerably less than those that form in the perfect bond case. Results of the computational simulation indicate that among the metal matrix composites examined, SiC/NiAl is the best candidate for high temperature applications at the given fiber volume ratio.

Lee, Ho-Jun↗

Code to Optimize Load Sharing of Split-Torque Transmissions Applied to the Comanche Helicopter

Most helicopters now in service have a transmission with a planetary design. Studies have shown that some helicopters would be lighter and more reliable if they had a transmission with a split-torque design instead. However, a split-torque design has never been used by a U.S. helicopter manufacturer because there has been no proven method to ensure equal sharing of the load among the multiple load paths. The Sikorsky/Boeing team has chosen to use a split-torque transmission for the U.S. Army's Comanche helicopter, and Sikorsky Aircraft is designing and manufacturing the transmission. To help reduce the technical risk of fielding this helicopter, NASA and the Army have done the research jointly in cooperation with Sikorsky Aircraft. A theory was developed that equal load sharing could be achieved by proper configuration of the geartrain, and a computer code was completed in-house at the NASA Lewis Research Center to calculate this optimal configuration.

Source record↗

SCEPTRE 2.7 User's Guide

Sandia’s Computational Engine for Particle Transport for Radiation Effects (SCEPTRE) is a computer code that solves the linear Boltzmann transport equation, particularly targeting coupled photon-electron problems. It uses unstructured finite element meshes in space, multigroup in energy, and discrete ordinates (Sn) or other methods in angle. SCEPTRE uses an xml-based input file to specify the problem. This report documents the options and syntax of that input file.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Euler solvers for transonic applications

The 1980s may well be called the Euler era of applied aerodynamics. Computer codes based on discrete approximations of the Euler equations are now routinely used to obtain solutions of transonic flow problems in which the effects of entropy and vorticity production are significant. Such codes can even predict separation from a sharp edge, owing to the inclusion of artificial dissipation, intended to lend numerical stability to the calculation but at the same time enforcing the Kutta condition. One effect not correctly predictable by Euler codes is the separation from a smooth surface, and neither is viscous drag; for these some form of the Navier-Stokes equation is needed. It, therefore, comes as no surprise to observe that the Navier-Stokes has already begun before Euler solutions were fully exploited. Moreover, most numerical developments for the Euler equations are now constrained by the requirement that the techniques introduced, notably artificial dissipation, must not interfere with the new physics added when going from an Euler to a full Navier-Stokes approximation. In order to appreciate the contributions of Euler solvers to the understanding of transonic aerodynamics, it is useful to review the components of these computational tools. Space discretization, time- or pseudo-time marching and boundary procedures, the essential constituents are discussed. The subject of grid generation and grid adaptation to the solution are touched upon only where relevant. A list of unanswered questions and an outlook for the future are covered.

Vanleer, Bram↗

Fast Reactor Methods Development Plan

Fast reactor design and safety analysis methods (computer codes) have been under development for approximately sixty years. During the first forty years, code development was driven by the immediate needs of operational facilities such as Experimental Breeder Reactors I and II and the Fast Flux Test Facility. Development was further driven by programmatic needs to support advanced reactor programs such as the Clinch River Breeder Reactor Project and the Advanced Liquid Metal Reactor Program. By the early 1990s, design and safety analysis codes had been well established to addressed the unique characteristics of fast reactors. During the subsequent twenty years, however, development was nearly stagnant as a result of U.S. Government policy changes that abruptly canceled nearly all of the fast reactor research in the United States. As a result, the state of fast reactor design and safety analysis methods has deteriorated.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Blade loss transient dynamics analysis, volume 2. Task 2: TETRA 2 user's manual

This is the user's manual for the TETRA 2 Computer Code, a program developed in the NASA-Lewis Blade Loss Program. TETRA 2 calculates a turbine engine's dynamic structural response from applied stimuli. The calculation options are: (1) transient response; and (2) steady state forced response. Based on the method of modal syntheses, the program allows the use of linear, as well as nonlinear connecting elements. Both transient and steady state options can include: flexible Bladed Disk Module, and Nonlinear Connecting Elements (including deadband, hardening/softening spring). The transient option has the additional capability to calculate response with a squeeze film bearing module. TETRA 2 output is summarized in a plotfile which permits post processing such as FFT or graphical animation with the proper software and computer equipment.

Black, Gerald↗

Application of a low order panel method to complex three-dimensional internal flow problems

An evaluation of the ability of a low order panel method to predict complex three-dimensional internal flow fields was made. The computer code VSAERO was used as a basis for the evaluation. Guidelines for modeling internal flow geometries were determined and the effects of varying the boundary conditions and the use of numerical approximations on the solutions accuracy were studied. Several test cases were run and the results were compared with theoretical or experimental results. Modeling an internal flow geometry as a closed box with normal velocities specified on an inlet and exit face provided accurate results and gave the user control over the boundary conditions. The values of the boundary conditions greatly influenced the amount of leakage an internal flow geometry suffered and could be adjusted to eliminate leakage. The use of the far-field approximation to reduce computation time influenced the accuracy of a solution and was coupled with the values of the boundary conditions needed to eliminate leakage. The error induced in the influence coefficients by using the far-field approximation was found to be dependent on the type of influence coefficient, the far-field radius, and the aspect ratio of the panels.

Ashby, D. L.↗