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,045 records · Page 58

Exascale Monte Carlo R&D [Slides]

Overview of this presentation is (1) Exascale computing - different technologies, getting there; (2) high-performance proof-of-concept MCMini - features and results; and (3) OpenCL toolkit - Oatmeal (OpenCL Automatic Memory Allocation Library) - purpose and features. Despite driver issues, OpenCL seems like a good, hardware agnostic tool. MCMini demonstrates the possibility for GPGPU-based Monte Carlo methods - it shows great scaling for HPC application and algorithmic equivalence. Oatmeal provides a flexible framework to aid in the development of scientific OpenCL codes.

97 MATHEMATICS AND COMPUTING↗

Automated metadata--final project report

This report summarizes the work of the Automated Metadata, Provenance Cataloging, and Navigable Interfaces: Ensuring the Usefulness of Extreme-Scale Data Project (MPO Project) funded by the United States Department of Energy (DOE), Offices of Advanced Scientific Computing Research and Fusion Energy Sciences. Initially funded for three years starting in 2012, it was extended for 6 months with additional funding. The project was a collaboration between scientists at General Atomics, Lawrence Berkley National Laboratory (LBNL), and Massachusetts Institute of Technology (MIT). The group leveraged existing computer science technology where possible, and extended or created new capabilities where required. The MPO project was able to successfully create a suite of software tools that can be used by a scientific community to automatically document their scientific workflows. These tools were integrated into workflows for fusion energy and climate research illustrating the general applicability of the project’s toolkit. Feedback was very positive on the project’s toolkit and the value of such automatic workflow documentation to the scientific endeavor.

96 KNOWLEDGE MANAGEMENT AND PRESERVATION↗

Propagation And Damping Of High Harmonic Fast Waves And Electron Cyclotron Waves In The Nstx-U-Device

For the spherical torus and for conventional tokamaks to be viable means of producing fusion energy there are two important requirements. The first is to heat the plasma in the experimental devices to high temperatures; preferably, to fusion relevant temperatures. The second is to drive a current in the plasma so that the fusion device operates in a steady-state mode. A pulsed fusion reactor for energy production is not desirable. A favored means to achieve both of these requirements is through the use of microwaves, commonly referred to as radio frequency waves. This proposal was to understand observations from the National Spherical Torus Experiment Upgrade (NSTX-U) device on radio frequency wave heating of the confined plasma during operations with high harmonic fast waves. The proposed research was also directed toward assisting in the optimization of the experimental modes of operation so as to maximize the heating efficiency. A detailed analysis of the experiments has shown that about one-third of the input radio frequency power does not make it into the core region of the plasma where heating is required. The losses are believed to occur primarily in the edge region where the radio frequency power is initially coupled into the plasma. Using sophisticated computational codes and theoretical models, that have been developed for the NSTX-U experimental program, this research quantified the physical processes which occur in the edge plasma and lead to a decrease in the heating efficiency. The effect of edge turbulence on the scattering of radio frequency waves was studied for waves in the ion cyclotron and electron cyclotron range of frequencies. For high harmonic fast waves in the ion cyclotron range of frequencies, simulations were carried out using a full-wave solver that was extended to the vacuum vessel wall. These simulations successfully reproduced trends that were observed experimentally. This simulation capability was also used to study the interaction of high harmonic fast waves with the scrape-off layer plasma. More recently, a new technique has been developed for the calculation of radio frequency waves in toroidal geometry that enables the simultaneous incorporation of antenna geometry, plasma facing components, the scrape off-layer, and core propagation in the high harmonic fast wave regime. The effect of turbulent plasmas on the mode conversion of extraordinary electromagnetic waves in the electron cyclotron range of frequencies in the NSTX-U plasma has been studied theoretically. In the over-dense plasmas of NSTX-U, the traditional electromagnetic waves are not effective for delivering energy and momentum to the core plasma. However, electron Bernstein waves generated by mode conversion are efficient carriers of wave energy and momentum to the core. The mode conversion occurs in the turbulent edge region. A theory for mode conversion in such plasmas has been developed. The parameters that can optimize the mode conversion can be readily determined. Also, the model is useful for experimentalists to analyze observation on a shot-to-shot basis.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Propagation And Damping Of High Harmonic Fast Waves And Electron Cyclotron Waves In The Nstx-U-Device

For the spherical torus and for conventional tokamaks to be viable means of producing fusion energy there are two important requirements. The first is to heat the plasma in the experimental devices to high temperatures; preferably, to fusion relevant temperatures. The second is to drive a current in the plasma so that the fusion device operates in a steady-state mode. A pulsed fusion reactor for energy production is not desirable. A favored means to achieve both of these requirements is through the use of microwaves, commonly referred to as radio frequency waves. This proposal was to understand observations from the National Spherical Torus Experiment Upgrade (NSTX-U) device on radio frequency wave heating of the confined plasma during operations with high harmonic fast waves. The proposed research was also directed toward assisting in the optimization of the experimental modes of operation so as to maximize the heating efficiency. A detailed analysis of the experiments has shown that about one-third of the input radio frequency power does not make it into the core region of the plasma where heating is required. The losses are believed to occur primarily in the edge region where the radio frequency power is initially coupled into the plasma. Using sophisticated computational codes and theoretical models, that have been developed for the NSTX-U experimental program, this research quantified the physical processes which occur in the edge plasma and lead to a decrease in the heating efficiency. The effect of edge turbulence on the scattering of radio frequency waves was studied for waves in the ion cyclotron and electron cyclotron range of frequencies. For high harmonic fast waves in the ion cyclotron range of frequencies, simulations were carried out using a full-wave solver that was extended to the vacuum vessel wall. These simulations successfully reproduced trends that were observed experimentally. This simulation capability was also used to study the interaction of high harmonic fast waves with the scrape-off layer plasma. More recently, a new technique has been developed for the calculation of radio frequency waves in toroidal geometry that enables the simultaneous incorporation of antenna geometry, plasma facing components, the scrape off-layer, and core propagation in the high harmonic fast wave regime. The effect of turbulent plasmas on the mode conversion of extraordinary electromagnetic waves in the electron cyclotron range of frequencies in the NSTX-U plasma has been studied theoretically. In the over-dense plasmas of NSTX-U, the traditional electromagnetic waves are not effective for delivering energy and momentum to the core plasma. However, electron Bernstein waves generated by mode conversion are efficient carriers of wave energy and momentum to the core. The mode conversion occurs in the turbulent edge region. A theory for mode conversion in such plasmas has been developed. The parameters that can optimize the mode conversion can be readily determined. Also, the model is useful for experimentalists to analyze observation on a shot-to-shot basis.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Radiological Safety Analysis Computer (RSAC) Program Version 7.4 Users’ Manual

The Radiological Safety Analysis Computer (RSAC) Program Version 7.4 (RSAC-7) is the newest version of the RSAC legacy code. RSAC-7 calculates the consequences of a release of radionuclides to the atmosphere. Users generates a fission product inventory from either reactor operating history or a nuclear criticality event. RSAC-7 models the effects of high-efficiency particulate air filters or other cleanup systems and calculates the decay and ingrowth during transport through processes, facilities, and the environment. Doses are calculated for inhalation, air immersion, ground surface, ingestion, and cloud gamma pathways. RSAC-7 is used as a tool to evaluate accident conditions in emergency response scenarios, radiological sabotage events, and safety basis accident consequences. This users’ manual contains the mathematical models and operating instructions for RSAC-7. Instructions, screens, and examples are provided to guide the user through the functions provided by RSAC-7. This program is designed for users who are familiar with radiological dose assessment methods.

73 - NUCLEAR PHYSICS AND RADIATION PHYSICS↗

ACCESS 1: Approximation Concepts Code for Efficient Structural Synthesis program documentation and user's guide

The program documentation and user's guide for the ACCESS-1 computer program is presented. ACCESS-1 is a research oriented program which implements a collection of approximation concepts to achieve excellent efficiency in structural synthesis. The finite element method is used for structural analysis and general mathematical programming algorithms are applied in the design optimization procedure. Implementation of the computer program, preparation of input data and basic program structure are described, and three illustrative examples are given.

Miura, H.↗

Aeroelastic analysis for helicopter rotors with blade appended pendulum vibration absorbers. Mathematical derivations and program user's manual

Mathematical development is presented for the expanded capabilities of the United Technologies Research Center (UTRC) G400 Rotor Aeroelastic Analysis. This expanded analysis, G400PA, simulates the dynamics of teetered rotors, blade pendulum vibration absorbers and the higher harmonic excitations resulting from prescribed vibratory hub motions and higher harmonic blade pitch control. Formulations are also presented for calculating the rotor impedance matrix appropriate to these higher harmonic blade excitations. This impedance matrix and the associated vibratory hub loads are intended as the rotor blade characteristics elements for use in the Simplified Coupled Rotor/Fuselage Vibration Analysis (SIMVIB). Sections are included presenting updates to the development of the original G400 theory, and material appropriate to the user of the G400PA computer program. This material includes: (1) a general descriptionof the tructuring of the G400PA FORTRAN coding, (2) a detaild description of the required input data and other useful information for successfully running the program, and (3) a detailed description of the output results.

Bielawa, R. L.↗

A Numerical Study of the Effects of Curvature and Convergence on Dilution Jet Mixing

An analytical program was conducted to assemble and assess a three-dimensional turbulent viscous flow computer code capable of analyzing the flow field in the transition liners of small gas turbine engines. This code is of the TEACH type with hybrid numerics, and uses the power law and SIMPLER algorithms, an orthogonal curvilinear coordinate system, and an algebraic Reynolds stress turbulence model. The assessments performed in this study, consistent with results in the literature, showed that in its present form this code is capable of predicting trends and qualitative results. The assembled code was used to perform a numerical experiment to investigate the effects of curvature and convergence in the transition liner on the mixing of single and opposed rows of cool dilution jets injected into a hot mainstream flow.

Holdeman, J. D.↗

A numerical study of the effects of curvature and convergence on dilution jet mixing

An analytical program was conducted to assemble and assess a three-dimensional turbulent viscous flow computer code capable of analyzing the flow field in the transition liners of small gas turbine engines. This code is of the TEACH type with hybrid numerics, and uses the power law and SIMPLER algorithms, an orthogonal curvilinear coordinate system, and an algebraic Reynolds stress turbulence model. The assessments performed in this study, consistent with results in the literature, showed that in its present form this code is capable of predicting trends and qualitative results. The assembled code was used to perform a numerical experiment to investigate the effects of curvature and convergence in the transition liner on the mixing of single and opposed rows of cool dilution jets injected into a hot mainstream flow.

Holdeman, J. D.↗

Aeroelastic tailoring for oblique wing lateral trim

Composite material aeroelastic tailoring is presently explored as a means for the correction of the roll trim imbalance of oblique-wing aircraft configurations. The concept is demonstrated through the analysis of a realistic oblique wing by a static aeroelastic computational procedure encompassing the full potential transonic aerodynamic code FLO22 and a Ritz structural plate program that models the stiffness due to symmetrical-but-unbalanced composite wing skins. Results indicate that asymetric composite tailoring reduces the aileron deflection needed for roll equilibrium, and reduces control surface hinge moment and drag. Wing skin stresses are, however, very high.

Bohlmann, Jonathan D.↗

NASTRAN migration to UNIX

COSMIC/NASTRAN, as it is supported and maintained by COSMIC, runs on four main-frame computers - CDC, VAX, IBM and UNIVAC. COSMIC/NASTRAN on other computers, such as CRAY, AMDAHL, PRIME, CONVEX, etc., is available commercially from a number of third party organizations. All these computers, with their own one-of-a-kind operating systems, make NASTRAN machine dependent. The job control language (JCL), the file management, and the program execution procedure of these computers are vastly different, although 95 percent of NASTRAN source code was written in standard ANSI FORTRAN 77. The advantage of the UNIX operating system is that it has no machine boundary. UNIX is becoming widely used in many workstations, mini's, super-PC's, and even some main-frame computers. NASTRAN for the UNIX operating system is definitely the way to go in the future, and makes NASTRAN available to a host of computers, big and small. Since 1985, many NASTRAN improvements and enhancements were made to conform to the ANSI FORTRAN 77 standards. A major UNIX migration effort was incorporated into COSMIC NASTRAN 1990 release. As a pioneer work for the UNIX environment, a version of COSMIC 89 NASTRAN was officially released in October 1989 for DEC ULTRIX VAXstation 3100 (with VMS extensions). A COSMIC 90 NASTRAN version for DEC ULTRIX DECstation 3100 (with RISC) is planned for April 1990 release. Both workstations are UNIX based computers. The COSMIC 90 NASTRAN will be made available on a TK50 tape for the DEC ULTRIX workstations. Previously in 1988, an 88 NASTRAN version was tested successfully on a SiliconGraphics workstation.

Chan, Gordon C.↗

Critical evaluation of Jet-A spray combustion using propane chemical kinetics in gas turbine combustion simulated by KIVA-II

Numerical solutions of the Jet-A spray combustion were obtained by means of the KIVA-II computer code after Jet-A properties were added to the 12 chemical species the program had initially contained. Three different reaction mechanism models are considered. The first model consists of 131 reactions and 45 species; it is evaluated by comparing calculated ignition delay times with available shock tube data, and it is used in the evaluation of the other two simplified models. The simplified mechanisms consider 45 reactions and 27 species and 5 reactions and 12 species, respectively. In the prediction of pollutants NOx and CO, the full mechanism of 131 reactions is considered to be more reliable. The numerical results indicate that the variation of the maximum flame temperature is within 20 percent as compared with that of the full mechanism of 131 reactions. The chemical compositions of major components such as C3H8, H2O, O2, CO2, and N2 are of the same order of magnitude. However, the concentrations of pollutants are quite different.

Nguyen, H. L.↗

Panel methods: An introduction

Panel methods are numerical schemes for solving (the Prandtl-Glauert equation) for linear, inviscid, irrotational flow about aircraft flying at subsonic or supersonic speeds. The tools at the panel-method user's disposal are (1) surface panels of source-doublet-vorticity distributions that can represent nearly arbitrary geometry, and (2) extremely versatile boundary condition capabilities that can frequently be used for creative modeling. Panel-method capabilities and limitations, basic concepts common to all panel-method codes, different choices that were made in the implementation of these concepts into working computer programs, and various modeling techniques involving boundary conditions, jump properties, and trailing wakes are discussed. An approach for extending the method to nonlinear transonic flow is also presented. Three appendices supplement the main test. In appendix 1, additional detail is provided on how the basic concepts are implemented into a specific computer program (PANAIR). In appendix 2, it is shown how to evaluate analytically the fundamental surface integral that arises in the expressions for influence-coefficients, and evaluate its jump property. In appendix 3, a simple example is used to illustrate the so-called finite part of the improper integrals.

Erickson, Larry L.↗

Human life support during interplanetary travel and domicile. II - Generic Modular Flow Schematic modeling

This paper describes the Generic Modular Flow Schematic (GMFS) architecture capable of encompassing all functional elements of a physical/chemical life support system (LSS). The GMFS can be implemented to synthesize, model, analyze, and quantitatively compare many configurations of LSSs, from a simple, completely open-loop to a very complex closed-loop. The GMFS model is coded in ASPEN, a state-of-the-art chemical process simulation program, to accurately compute the material, heat, and power flow quantities for every stream in each of the subsystem functional elements (SFEs) in the chosen configuration of a life support system. The GMFS approach integrates the various SFEs and subsystems in a hierarchical and modular fashion facilitating rapid substitutions and reconfiguration of a life support system. The comprehensive ASPEN material and energy balance output is transferred to a systems and technology assessment spreadsheet for rigorous system analysis and trade studies.

Farral, Joseph F.↗

Analysis of aeroelastic and resonance responses of a wind tunnel model support system

Vibrations of the National Transonic Facility model support system were analyzed using CAP-TSD (Computational Aeroelasticity Program-Transonic Small Disturbance), a three-dimensional transonic small disturbance potential code designed for aeroelastic analysis of complex configurations. The model support system was represented as a tunnel-spanning flexible wing whose structural properties were obtained from measured responses. Aeroelastic transients were calculated, assuming no structural damping, and analyzed to obtain modal stability characteristics. The results showed that there is a tendency toward a hump mode instability in one of the modes and that there are other weakly damped or unstable modes. The close proximity of the solid wind-tunnel diffuser side walls to the model support system was shown to have adverse effects on its aeroelastic stability characteristics. When the effects of the wind-tunnel walls were modeled, the analyses indicated the presence of wind-tunnel resonance modes which were not present when free-air boundary conditions were used.

Whitlow, Woodrow, Jr.↗

Numerical thermal analyses of heat exchangers for the stirling engine application

The Regenerator, Cooler and Heater for the NASA Space Power Research Engine (SPRE) have been analyzed in detail for laminar, incompressible and oscillatory flow conditions. Each component has been analyzed independently and in detail with the regenerator being modeled as two-parallel-plates channel with a solid wall. The ends of the channel are exposed to two reservoir maintained at different temperature thus providing an axial temperature gradient along the channel. The cooler and heater components have been modeled as circular pipes with isothermal walls. Two different types of thermal boundary conditions have been investigated for the cooler and heater, namely, symmetric and asymmetric temperature inflow. In symmetric temperature inflow the flow enters the channel with the same temperature in throughout the velocity cycle whereas, in asymmetric temperature inflow the flow enters with a different temperature in each half cycle. The study was conducted over a wide range of Maximum Reynolds number (RE(max) varying from 75 to 60000, Valensi number (Va) from 2.5 to 800, and relative amplitude of fluid displacement (A(sub r) from 0.357 to 1.34. A two dimensional Finite volume method based on the SIMPLE algorithm was used to solve the governing partial differential equations. Post processing programs were developed to effectively describe the heat transfer mechanism under oscillatory flows. The computer code was validated by comparing with existing analytical solutions for oscillating flows. The thermal field have been studied with the help of temperature contour and three dimensional plots. The instantaneous friction factor, wall heat flux and heat transfer coefficient have been examined. It has been concluded that in general, the frictional factor and heat transfer coefficient are higher under oscillatory flow conditions when the Valensi number is high. Also, the thermal efficiency decreases for lower A(r) values. Further, the usual steady state definition for the heat transfer coefficient does not seem to be valid.

Kannapareddy, Mohan Raj↗

Computational Modeling of Multi-Phase/Multi-Species Flows with Applications to Liquid Rocket Engines

Accurate prediction of all physical phenomena in a combustion chamber is essential for better understanding of the system performance. Atomization, evaporation, combustion, chemical kinetics, and turbulence are those processes of great importance that need to be well understood. Processes involving the liquid phase in a combustion chamber will be further complicated under supercritical conditions. More advanced and accurate numerical techniques are required to extend our understanding of the above phenomena. A computer program for multi-species/multi-phase flow was developed for NASA/MSFC in 1992. This code, called Liquid Thrust Chamber Performance (LTCP) program takes an Eulerian- Eulerian approach and is based on the Total Variation Diminishing (TVD) technique with Lax-Friedrichs upwind method. Under the NASA/ASEE SFFP the LTCP code was used to predict the performance characteristics of several engines that were of particular interest to NASA. This code was also successful in a combustion detonation study. Converting the program to the PC platform was accomplished which extends usability and makes it available to a wider range of users. The Eulerian formulation of the liquid phase provides a suitable model that can be extended to include combustion modeling under supercritical conditions. The results have been compared against the ones of other codes and available measured data. The algorithm proved to be robust and efficient for problems with stiff source terms.

Navaz, Homayun K.↗

CFD Evaluation of Aviation Fuels for Commercial Supersonics Technology

An overview is provided of a CFD analysis of a Lean-Direct Injection (LDI) combustion concept for NASA’s Commercial Supersonic Transport (CST) program. The National Combustion Code (OpenNCC) was used to perform two-phase reacting flow computations with three different aviation fuels for UTRC’s Axially Controlled Stoichiometry (ACS) combustor at CST cruise conditions. All computations were performed with a consistent approach of mesh-generation, spray and chemical-kinetics modeling for the combustor. Emissions (EICO, EINOx) and performance characteristics were predicted with OpenNCC, to evaluate the impact of replacing Jet-A with RP2 and Iso-Paraffinic Kerosene (IPK) for supersonic combustor applications. The CFD predictions of flame structure for IPK fuel were significantly different in comparison to those of Jet-A and RP2 fuel. However, the CFD analysis predicted very similar EINOx emissions for Jet-A, RP-2 and IPK fuels. The predicted EINOx emissions for all three fuels also compared very well with measured experimental data.

gas turbine combustion↗