Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Coded computing”

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,315 records · Page 73
All sourcesNASA NTRSDOE OSTI
NASA NTRS1995

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.↗
DOE OSTI2023

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↗
DOE OSTI2022

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↗
NASA NTRS1987

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.↗
NASA NTRS1994

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.↗
NASA NTRS2000

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↗
NASA NTRS1989

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.↗
NASA NTRS1990

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.↗
NASA NTRS1990

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↗
NASA NTRS1995

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↗
DOE OSTI2024

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↗
NASA NTRS1989

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↗
NASA NTRS1986

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↗
NASA NTRS1986

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.↗
NASA NTRS1986

Numerical studies of motion and decay of vortex filaments

A computational code is developed for the integro-differential equations governing the motion of the centerlines of vortex filaments submerged in a background potential flow. These equations, which are derived from the method of matched asymptotic analysis, include the effect of decaying large-magnitude circumferential and axial velocity components in the vortical cores. Numerical examples are presented to assess the effect of large axial velocity and of nonsimilar initial profiles in vortical cores. The initial configurations of the filaments are chosen so as to fulfill the basic assumption of asymptotic analysis, which is the effective vortical core size is much smaller than all other length scales in the flowfield, e.g., the radius of curvature and interfilament distance. The computations are continued until the basic assumption is no longer valid, that is, when the merging or intersection of filaments have begun. Various types of local or global merging or intersection of filaments are classified and demonstrated by numerical examples.

Liu, C. H.↗
DOE OSTI2022

Stochastic Modeling in a Multimaterial Continuum Mixture Shock Physics Code

Stochastic modelling approaches are presented to capture random effects at multiple time and length scales. Random processes that occur at the microscale produce nondeterministic effects at the macroscale. Here we present three stochastic modeling approaches that describe random processes at microscopic length scales and map these processes to the macroscopic length scale. The first stochastic modeling approach is based upon a particle based numerical technique to solve a Stochastic Differential Equation (SDE) using an arbitrary diffusion process to capture random processes at the microstructural level. The second approach prescribes a Probability Density Function (PDF) for the drift and diffusion of the random variable derived using the forward and backward Kolmogorov equations. This method requires mean and drift evolution PDF transport equations. The third approach is the coupling of multiple random variables which are dependent on each other. The relationship of the PDFs and a coupling function, known as a copula, produces a Joint Probability Density Function (JPDF). These stochastic modeling approaches are implemented into a Multiple Component (MC) shock physics computational code and used to model statistical fracture and reactive flow applications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1984

Direct simulations of chemically reacting turbulent mixing layers

The report presents the results of direct numerical simulations of chemically reacting turbulent mixing layers. The work consists of two parts: (1) the development and testing of a spectral numerical computer code that treats the diffusion reaction equations; and (2) the simulation of a series of cases of chemical reactions occurring on mixing layers. The reaction considered is a binary, irreversible reaction with no heat release. The reacting species are nonpremixed. The results of the numerical tests indicate that the high accuracy of the spectral methods observed for rigid body rotation are also obtained when diffusion, reaction, and more complex flows are considered. In the simulations, the effects of vortex rollup and smaller scale turbulence on the overall reaction rates are investigated. The simulation results are found to be in approximate agreement with similarity theory. Comparisons of simulation results with certain modeling hypotheses indicate limitations in these hypotheses. The nondimensional product thickness computed from the simulations is compared with laboratory values and is found to be in reasonable agreement, especially since there are no adjustable constants in the method.

Riley, J. J.↗
NASA NTRS1982

A simple method to compute the change in earth-atmosphere radiative balance due to a stratospheric aerosol layer

A computer code was developed in terms of a three-layer model for the earth-atmosphere system, using a two-stream approximation for the troposphere and stratosphere. The analysis was limited to variable atmosphere loading by solar radiation over an unperturbed section of the atmosphere. The scattering atmosphere above a Lambertian ground layer was considered in order to derive the planar albedo and the spherical albedo. Attention was given to the influence of the aerosol optical thickness in the stratosphere, the single scattering albedo and asymmetry factor, and the sublayer albedo. Calculations were performed of the zonal albedo and the planetary radiation balance, taking into account a stratospheric aerosol layer containing H2SO4 droplets and volcanic ash. The resulting ground temperature disturbance was computed using a Budyko (1969) climate model. Local decreases in the albedo in the summer were observed in high latitudes, implying a heating effect of the aerosol. An accompanying energy loss of 23-27 W/sq m was projected, which translates to surface temperature decreases of either 1.1 and 0.45 C, respectively, for background and volcanic aerosols.

Lenoble, J.↗