SEARCH · Search NASA
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.
The purpose of this report is to document the theoretical models utilized by the computer code MassTran. This report will focus on the theoretical models used to analyze high Mach number, fully compressible, transonic flows in pipes and networks.
Coupling generator current to electron-beam loads and dynamic-inductance loads (z pinches) is greatly improved with low-inductance vacuum transmission lines. Typically, magnetically insulated transmission lines (MITLs) utilize nearly constant impedance geometries where the minimum impedance (anode-cathode gap) is limited by electron loss during the setup of magnetic insulation. Constant-impedance MITL designs are not optimized for low inductance (Zτ) when given a constraint on vacuum electron flow. We describe the design and the modeling of MITLs that have a nonconstant impedance. The baseline variable-impedance MITL designs are based on equilibrium vacuum-electron flow models [1] that are incorporated into the Screamer circuit code [2]. Such designs satisfy the MITL limits in the areal density of the loss current and have magnetically insulated flow that is stable everywhere along the MITL. Different MITL designs are then modeled analytically and with a 2-D particle-in-cell computer code in order to quantify subtle electron losses driven by changes to the vacuum impedance as a function of distance along the transmission line. We present a variable-impedance MITL design that delivers more current with a lower inductance than found in constant-impedance MITL designs.
Finite-element, transient heat transfer analyses were performed for the first-stage blades of the space shuttle main engine (SSME) high-pressure fuel turbopump. The analyses were based on test engine data provided by Rocketdyne. Heat transfer coefficients were predicted by performing a boundary-layer analysis at steady-state conditions with the STAN5 boundary-layer code. Two different peak-temperature overshoots were evaluated for the startup transient. Cutoff transient conditions were also analyzed. A reduced gas temperature profile based on actual thermocouple data was also considered. Transient heat transfer analyses were conducted with the MARC finite-element computer code.
A Navier-Stokes computer code was validated using a number of two- and three-dimensional configurations for both laminar and turbulent flows. The validation data covers a range of freestream Mach numbers from 3 to 14, includes wall pressures, velocity profiles, and skin friction. Nozzle flow fields computed for a generic scramjet nozzle from Mach 3 to 20, wall pressures, wall skin friction values, heat transfer values, and overall performance are presented. In addition, three-dimensional solutions obtained for two asymmetric, single expansion ramp nozzles at a pressure ratio of 10 consists of the internal expansion region in the converging/diverging sections and the external supersonic exhaust in a quiescent ambient environment. The fundamental characteristics that were captured successfully include expansion fans; Mach wave reflections; mixing layers; and nonsymmetrical, multiple inviscid cell, supersonic exhausts. Comparison with experimental data for wall pressure distributions at the center planes shows good agreement.
A Navier-Stokes computer code was validated using a number of two- and three-dimensional configurations for both laminar and turbulent flows. The validation data covers a range of freestream Mach numbers from 3 to 14, including wall pressures, velocity pressure, and skin friction. Nozzle flow fields computed for a generic scramjet nozzle from Mach 3 to 20, wall pressures, wall skin friction values, heat transfer values, and overall performance are presented. In addition, three-dimensional solutions obtained for two asymmetric, single expansion ramp nozzles at a pressure ratio of 10 consists of the internal expansion region in the converging/diverging sections and the external superonic exhaust in a quiescent ambient environment. The fundamental characteristics that were captured successfully include expansion fans; Mach wave reflections; mixing layers; and nonsymmetrical, multiple inviscid cell, supersonic exhausts. Comparison with experimental data for wall pressure distributions at the center planes shows good agreement.
A computer code that generates shock-free transonic compressor cascade shapes while taking into account viscosity effects is developed. The mathematical model for the inviscid flow field is the full potential equation. The Kutta-Joukowski condition is satisfied by varying the free stream angle at downstream infinity. A boundary fitted computational grid of C-type is generated using a sequence of conformal mapping and nonorthogonal coordinate stretching and shearing transformations. The inviscid calculation is performed sequentially on up to four consecutively refined grids thereby accelerating the convergence of the solution process. The full potential equation is solved using a finite area technique and rotated, type-dependent finite differencing. Artificial viscosity of the first order is added in a fully conservative form. Shock-free cascade airfoil shapes are obtained using the fictitious gas concept of Sobieczky and the method of characteristics in the rheograph plane. Viscosity effects are incorporated via a boundary layer displacement thickness. The integral boundary layer code is based on Rotta's turbulence model and assumes transition region of zero length.
The primary objective of this study is to validate the system analysis code, SAM, using experimental data from the Compact Integral Effects Test (CIET) experimental loop. SAM is a modern system analysis code being developed at Argonne National Laboratory for safety analysis of designs for advanced non-light water reactors (non-LWRs), such as sodium-cooled fast reactors, high-temperature gas-cooled reactors, and fluoride salt-cooled high-temperature reactors (FHRs). To support SAM code development for the wide range of non-LWR applications, it is of paramount importance to validate the code against experiments highly relevant to these reactor concepts. Additionally, the CIET facility, which was designed to reproduce the thermal-hydraulics response of FHRs under both forced- and natural-circulation conditions, has been identified and selected as one of the benchmark test facilities for SAM code validation. In this study, two sets of available CIET tests were selected for SAM code validation purposes, namely, power step change transient tests and steady-state natural-circulation tests. For all selected tests, 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.
ATS platforms are some of the largest, most complex, and most expensive computer systems installed in the United States at just a few major national laboratories. This milestone describes our recent efforts to procure, install, and test a machine called Vortex at Sandia National Laboratories that is compatible with the larger ATS platform Sierra at LLNL. In this milestone, we have 1) configured and procured a machine with similar hardware characteristics as Sierra ATS, 2) installed the machine, verified its physical hardware, and measured its baseline performance, and 3) demonstrated the machine’s compatibility with Sierra ATS, and capacity for useful development and testing of Sandia computer codes (such as SPARC), including uses such as nightly regression testing workloads.
This report documents the library of thermodynamic data used with the NASA Glenn computer program CEA (Chemical Equilibrium with Applications). This library, containing data for over 2000 solid, liquid, and gaseous chemical species for temperatures ranging from 200 to 20,000 K, is available for use with other computer codes as well. The data are expressed as least-squares coefficients to a seven-term functional form for C((sup o)(sub p)) (T) / R with integration constants for H (sup o) (T) / RT and S(sup o) (T) / R. The NASA Glenn computer program PAC (Properties and Coefficients) was used to calculate thermodynamic functions and to generate the least-squares coefficients. PAC input was taken from a variety of sources. A complete listing of the database is given along with a summary of thermodynamic properties at 0 and 298.15 K.
A methodology for designing active control laws in a computational aeroelasticity environment is given. The methodology involves employing a systems identification technique to develop an explicit state-space model for control law design from the output of a computational aeroelasticity code. The particular computational aeroelasticity code employed in this paper solves the transonic small disturbance aerodynamic equation using a time-accurate, finite-difference scheme. Linear structural dynamics equations are integrated simultaneously with the computational fluid dynamics equations to determine the time responses of the structure. These structural responses are employed as the input to a modern systems identification technique that determines the Markov parameters of an "equivalent linear system". The Eigensystem Realization Algorithm is then employed to develop an explicit state-space model of the equivalent linear system. The Linear Quadratic Guassian control law design technique is employed to design a control law. The computational aeroelasticity code is modified to accept control laws and perform closed-loop simulations. Flutter control of a rectangular wing model is chosen to demonstrate the methodology. Various cases are used to illustrate the usefulness of the methodology as the nonlinearity of the aeroelastic system is increased through increased angle-of-attack changes.
In current practice, it is often difficult to draw firm conclusions about turbulence model accuracy when performing multi-code CFD studies ostensibly using the same model because of inconsistencies in model formulation or implementation in different codes. This paper describes an effort to improve the consistency, verification, and validation of turbulence models within the aerospace community through a website database of verification and validation cases. Some of the variants of two widely-used turbulence models are described, and two independent computer codes (one structured and one unstructured) are used in conjunction with two specific versions of these models to demonstrate consistency with grid refinement for several representative problems. Naming conventions, implementation consistency, and thorough grid resolution studies are key factors necessary for success.
The main goal of this report is to advance the development of methodologies for scattering by airborne composite vehicles. Although the primary focus continues to be the development of a general purpose computer code for analyzing the entire structure as a single unit, a number of other tasks are also being pursued in parallel with this effort. One of these tasks discussed within is on new finite element formulations and mesh termination schemes. The goal here is to decrease computation time while retaining accuracy and geometric adaptability.The second task focuses on the application of wavelets to electromagnetics. Wavelet transformations are shown to be able to reduce a full matrix to a band matrix, thereby reducing the solutions memory requirements. Included within this document are two separate papers on finite element formulations and wavelets.
Massive black hole (MBH) binaries are found at the centers of most galaxies. MBH mergers trace galaxy mergers and are strong sources of gravitational waves. 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 presentation shows how a spacetime is constructed on a computer to build a simulation laboratory for binary black hole mergers. Focus is on the recent advances that that reveal these waveforms, and the potential for discoveries that arises when these sources are observed by LIGO and LISA.
A computer code based on a finite-element solution algorithm is developed to perform an analytical investigation on the turbulent mixing and reaction of hydrogen jets injected from multiple orifices transverse and parallel to a supersonic airstream. A laser optical cavity flow field was also analyzed to demonstrate the generality of the proposed model. Computational results provide a three-dimensional description of velocity, temperature, and species-concentration fields downstream of injection. Major conclusions are that the analysis has immediate utility in evaluating the mixing effectiveness of transverse H2 injection data since it has been tested in its ability to model this type of data and that turbulent mixing length theory, constant effective Prandtl number, and a Lewis number of unity provide reasonable agreement with transverse H2 injection data downstream of the near-injection region. Efforts are presently being directed toward using the code in modeling laser and scramjet data for a wide range of flow conditions.
The need for a validation technique for computational fluid dynamics (CFD) codes in STOVL applications has led to research efforts to apply infrared thermal imaging techniques to visualize gaseous flow fields. Specifically, a heated, free-jet test facility was constructed. The gaseous flow field of the jet exhaust was characterized using an infrared imaging technique in the 2 to 5.6 micron wavelength band as well as conventional pitot tube and thermocouple methods. These infrared images are compared to computer-generated images using the equations of radiative exchange based on the temperature distribution in the jet exhaust measured with the thermocouple traverses. Temperature and velocity measurement techniques, infrared imaging, and the computer model of the infrared imaging technique are presented and discussed. From the study, it is concluded that infrared imaging techniques coupled with the radiative exchange equations applied to CFD models are a valid method to qualitatively verify CFD codes used in STOVL applications.
Big Bang Nucleosynthesis (BBN), the process of creation of lightest elements in the early universe, is a highly robust, precise, and ultimately successful theory that forms one of the three pillars of the standard hot-Big-Bang cosmological model. Existing theoretical treatments of BBN and the associated computer codes are accurate and flexible, but are typically highly technical and opaque, and not suitable for pedagogical understanding of the BBN. Here we present BBN-simple - a from-scratch numerical calculation of the lightest element abundances pitched at an advanced undergraduate or beginning graduate level. We review the physics of the early universe relevant for BBN, provide information about the reaction rates, and discuss computational-mathematics background that is essential in setting up a BBN calculation. Here, we calculate the abundances of the principal nuclear species in a standard cosmological model, and find a reasonably good agreement with public precision-level BBN codes.
Electronic compensation of reflector surface distortion using array feed with individual amplitude and phase control of the array elements is becoming increasingly attractive because of the recent advances in monolithic microwave integrated circuit (MMIC) technology. An algorithm has been developed previously using the concept of focal plane conjugate field matching in the receive mode and a computer code has been generated that predicts the proper excitation coefficients for the elements of the reflector feed array to compensate the effects of reflector surface distortion. This paper presents the results of an experimental study to verify the above compensation algorithm and in general to demonstrate the effectiveness of the array feed compensation technique.
In an attempt to provide a set of accurate standard test problems for computational code developers, a series of inviscid airfoil test cases were chosen by the AGARD Working Group 07, a subpanel of the AGARD Fluid Dynamics Panel. The cases include three different airfoils at transonic to supersonic conditions. A large number of international experts responded in this effort with computations that have been contrasted for accuracy and consistency. This paper is a summary of the authors' contribution to this study. In particular, the important aspects of the solution process that made it possible to obtain the high level of accuracy needed in this study are stressed.