SEARCH · Engineering Papers
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 Basic Linear Algebra Subprograms (BLAS) library is a collection of 38 FORTRAN-callable routines for performing basic operations of numerical linear algebra. BLAS library is portable and efficient source of basic operations for designers of programs involving linear algebriac computations. BLAS library is supplied in portable FORTRAN and Assembler code versions for IBM 370, UNIVAC 1100 and CDC 6000 series computers.
This paper briefly describes three user-oriented computer codes developed at the Ohio State University ElectroScience Laboratory: namely, the On-Aircraft Antenna Code, the Basic Scattering Code, and the Reflector Antenna Code. Each code uses the Geometrical Theory of Diffraction (GTD). An example is used to describe the basic GTD approach for all three codes. A brief description is given for each code which includes sample results compared to measurements.
Large-scale integrated (LSI) devices fabricated with integrated injection logic (I2L) were studied by both heavy ion experiments and device analysis to determine the upset threshold, the charge collection volumes and ultimately the probability for upset in the galactic cosmic ray spectrum. The devices studied, the SEP9900A, the SBP9989 and the P-Code generator were fabricated by Texas Instruments, Inc.
The principal of multiple material layering is well known as an effective method of reducing heat transfer; however, thermal gradients can impose significant loads and lead to delamination and subsequent component failure. An analysis is developed and experimental data are discussed for the thermomechanical effects of multilayered materials on a heat sink substrate of cylindrical geometry subject to thermal cycling. The geometry is heated in cross-flow by a high-velocity flame and cooled in crossflow by ambient-temperature air from a critical flow orifice. Each layer of material possesses a threshold beyond which small changes in temperature or mechanical loading greatly influence the life in thermal cycling of the layered materials. Comparisons are made between the thermomechanical loads predicted by various numerical codes for the linear case and by a simplified analytic model.
Laminar heat-transfer rates were measured on spherically blunted, 13 degrees/F degrees on-axis and bent biconics (fore cone bent 7 degrees upward relative to aft cone) at hypersonic-hypervelocity flow conditions in the Langley Expansion Tube. Freestream velocities from 4.5 to 6.9 km/sec and Mach numbers from 6 to 9 were generated using helium, nitrogen, air, and carbon dioxide test gases, resulting in normal shock density ratios from 4 to 19. Angle of attack, referenced to the axis of the aft cone, was varied from zero to 20 degrees in 4 degree increments. The effect of nose bend, angle of attack, and real-gas phenomena on heating distributions are presented along with comparisons of measurement to prediction from a code which solves the three-dimensional 'parabolized Navier-Stokes' equations.
A computer code has been developed that calculates radiative properties of nonequilibrium air in the low-density regimes expected during the flight of aeroassisted, orbital transfer vehicles. From the given nonequilibrium thermodynamic state variables, the code calculates number densities of internal states and the accompanying emission and absorption characteristics. In addition, the code calculates the number density of the hypothetical gas in radiative equilibrium that produces the same radiation emission as the given nonequilibrium gas. Sample results are shown to demonstrate how the code predicts the nonlinear variation of radiation with density at low densities.
Transonic two-dimensional cascade flows have been analyzed using viscous/inviscid coupling concepts. A full potential cascade code is coupled with an inverse integral boundary layer/wake method that permits calculation of separated laminar or turbulent flow. The semi-inverse coupling method of Wigton converges slowly in the case of a strong shock in the region between the shock and the trailing edge. The location of a strong shock is not well predicted by the coupling method, which indicates the need for an entropy correction in the potential code or the inclusion of a shock-boundary layer interaction module.
A supplement to the data relating to the calculation of nonequilibrium radiation in flight regimes of aeroassisted orbital transfer vehicles contains the listings of the computer code NEQAIR (Nonequilibrium Air Radiation), its primary input data, and explanation of the user-supplied input variables. The user-supplied input variables are the thermodynamic variables of air at a given point, i.e., number densities of various chemical species, translational temperatures of heavy particles and electrons, and vibrational temperature. These thermodynamic variables do not necessarily have to be in thermodynamic equilibrium. The code calculates emission and absorption characteristics of air under these given conditions.
Detailed longitudinal and circumferential heating distributions were measured on a spherically blunted, 13 deg/7 deg biconic at angles of attack from 0 deg to 27 deg. The measurements were made in the Langley 20-inch Mach 6 tunnel and 31-inch Mach 10 tunnel to provide heating distributions at Mach 6 and 10 in air. The free-stream Reynolds number based on model length varied from 0.4 to 4.8 million at Mach 6 and 0.3 to 1.3 million at Mach 10. The basic trends in leeward heating are found to be described by a single parameter, the viscous interaction parameter. A decrease in this parameter results in a decrease in leeward heating for attached flow, an increase in heating for separated flow, and a forward movement of the separation region. The parabolized Navier-Stokes code for laminar flow qualitatively predicted the longitudinal and circumferential heating distributions.
Laminar heating distributions have been measured on a 1.9 percent scale model of a generic aeroassisted vehicle taking the shape of a spherically blunted, 13-deg/7-deg biconic whose forecone section is bent upward (by 7 deg) to furnish self-trim capability at a 20-deg angle-of-attack. The results thus obtained were compared with data gathered for a straight biconic. While no Reynolds number effect on heating was noted on the windward side of the forecone, the opposite was true of the leeward side, where a Reynolds number increase caused circumferential flow separation at lower angles of attack. Generally, windward heating was predicted to within 10 percent with a computer code solving the steady, three-dimensional parabolized Navier-Stokes equations.
A preliminary quantitative study of the compressible flowfield in a steady, nonreacting model SCRAMJET combustor using laser-induced iodine fluorescence (LIIF) is reported. Measurements of density, temperature, and velocity were conducted with the calibrated, nonintrusive, optical technique for two different combustor operating conditions. First, measurements were made in the supersonic flow over a rearward-facing step without transverse injection for comparison with calculated pressure profiles. The second configuration was staged injection behind the rearward-facing step at an injection dynamic pressure ratio of 1.06. These experimental results will be used to validate computational fluid dynamic (CFD) codes being developed to model supersonic combustor flowfields.
This paper presents the results of a recent flight experiment conducted to evaluate surface-mounted hot-film sensors designed to detect Tollmien-Schlichting instability in the boundary layer in flight. Flight data are presented to illustrate the design and operation of the sensors and to provide information on disturbance growth and transition mode. The results of calculations using a boundary-layer analysis code and boundary-layer linear stability analysis methods are included.
The multitasking of a finite-difference scheme using multiple overset meshes is described. In this chimera, or multiple overset mesh approach, a multiple body configuration is mapped using a major grid about the main component of the configuration, with minor overset meshes used to map each additional component. This type of code is well suited to multitasking. Both steady and unsteady two dimensional computations are run on parallel processors on a CRAY-X/MP 48, usually with one mesh per processor. Flow field results are compared with single processor results to demonstrate the feasibility of running multiple mesh codes on parallel processors and to show the increase in efficiency.
A calibrated, nonintrusive optical technique, laser-induced iodine fluorescence (LIIF) was used to quantify the steady, compressible flowfield of a nonreacting, supersonic combustor. The combustor was configured with single and staged, transverse-air injection into a supersonic-air freestream behind a rearward-facing step. Pressure, temperature, two-velocity components, and injectant mole fraction were measured with high spatial resolution in the three-dimensional flowfields. These experimental results provide a benchmark set of data for validation of computational fluid dynamic (CFD) codes being developed to model supersonic combustor flowfields.
The straightforward automatic-differentiation and the hand-differentiated incremental iterative methods are interwoven to produce a hybrid scheme that captures some of the strengths of each strategy. With this compromise, discrete aerodynamic sensitivity derivatives are calculated with the efficient incremental iterative solution algorithm of the original flow code. Moreover, the principal advantage of automatic differentiation is retained (i.e., all complicated source code for the derivative calculations is constructed quickly with accuracy). The basic equations for second-order sensitivity derivatives are presented; four methods are compared. Each scheme requires that large systems are solved first for the first-order derivatives and, in all but one method, for the first-order adjoint variables. Of these latter three schemes, two require no solutions of large systems thereafter. For the other two for which additional systems are solved, the equations and solution procedures are analogous to those for the first order derivatives. From a practical viewpoint, implementation of the second-order methods is feasible only with software tools such as automatic differentiation, because of the extreme complexity and large number of terms. First- and second-order sensitivities are calculated accurately for two airfoil problems, including a turbulent flow example; both geometric-shape and flow-condition design variables are considered. Several methods are tested; results are compared on the basis of accuracy, computational time, and computer memory. For first-order derivatives, the hybrid incremental iterative scheme obtained with automatic differentiation is competitive with the best hand-differentiated method; for six independent variables, it is at least two to four times faster than central finite differences and requires only 60 percent more memory than the original code; the performance is expected to improve further in the future.
We report new predictions for the EUV spectral emission of FeIX-FeXXIV, based on data now available from the Solar EUV Rocket Telescope and Spectrograph (SERTS) and the Extreme Ultraviolet Explorer (EUVE) spectrometers. The iron spectral emission model is the first result of a larger effort to revise the Raymond & Smith model and to update the atomic rates. We present here predicted emissivities for selected densities and temperatures applicable to various astrophysical plasmas. Comparisons of our predicted spectra with two recent observations provide important tests of the atomic data. They also test to some extent some basic assumptions of coronal emission codes: optically thin spectral lines and ionization equilibrium.
NASA is developing a 'tool box' that includes a number of advanced structural analysis computer codes which, taken together, represent the comprehensive fracture mechanics capability required to predict the onset of widespread fatigue damage. These structural analysis tools have complementary and specialized capabilities ranging from a finite-element-based stress-analysis code for two- and three-dimensional built-up structures with cracks to a fatigue and fracture analysis code that uses stress-intensity factors and material-property data found in 'look-up' tables or from equations. NASA is conducting critical experiments necessary to verify the predictive capabilities of the codes, and these tests represent a first step in the technology-validation and industry-acceptance processes. NASA has established cooperative programs with aircraft manufacturers to facilitate the comprehensive transfer of this technology by making these advanced structural analysis codes available to industry.
Observations regarding the use of advanced computational fluid dynamics (CFD) analysis, sensitivity analysis (SA), and design codes in gradient-based multidisciplinary design optimization (MDO) reflect our perception of the interactions required of CFD and our experience in recent aerodynamic design optimization studies using CFD. Sample results from these latter studies are summarized for conventional optimization (analysis - SA codes) and simultaneous analysis and design optimization (design code) using both Euler and Navier-Stokes flow approximations. The amount of computational resources required for aerodynamic design using CFD via analysis - SA codes is greater than that required for design codes. Thus, an MDO formulation that utilizes the more efficient design codes where possible is desired. However, in the aerovehicle MDO problem, the various disciplines that are involved have different design points in the flight envelope; therefore, CFD analysis - SA codes are required at the aerodynamic 'off design' points. The suggested MDO formulation is a hybrid multilevel optimization procedure that consists of both multipoint CFD analysis - SA codes and multipoint CFD design codes that perform suboptimizations.