Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “code comparisons”

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,171 records · Page 65
All sourcesNASA NTRSDOE OSTI
NASA NTRS1997

Analysis of Emission Spectra from Arc-jet Shock Layer Flows

This paper reports computational comparisons with experimental studies of a nonequilibrium blunt body shock layer in a high enthalpy arc-jet wind tunnel at NASA Ames Research Center. The primary objective of this work is to investigate the existence of a thermochemical equilibrium region in the shock layer. The existence of such an equilibrium region is of interest for following reasons: (1) to understand the equilibration process behind the shock in an arc-jet flow environment; (2) to interpret measured surface heat transfer data for purpose of determining surface catalytic efficiency, and (3) to determine the total enthalpy from the spectroscopic measurements. The paper will present an analysis of the experimental data obtained in the arc-jet wind tunnel. Experimental data includes measurements of emission spectra of radiation emanating from a shock layer formed in front of a 6-inch flat-faced cylinder. The measurements, obtained using a two dimensional CCD camera mounted on a spectrograph, provide spatially resolved spectra along the stagnation streamline of the model. Computational analysis includes simulation of nonequilibrium flow in the arc-jet facility (flow in the conical nozzle and shock layer in front of a flat-faced cylinder) using 2-D/axisymmetric Navier-Stokes codes and prediction of the radiation spectra from the axisymmetric flowfield using NEQAIR radiation code. Various line-of-sight averaged flow properties such as vibrational and rotational temperatures, species number densities within the shock layer are deduced from the experimental spectra. Comparison of the computed and experimental line-of-sight averaged flow properties provides assessment of thermochemical equilibration processes in an arc-jet shock layer.

Gokcen, T.↗
NASA NTRS2003

Investigation of the Flow Physics Driving Stall-Side Flutter in Advanced Forward Swept Fan Designs

Flutter-free operation of advanced transonic fan designs continues to be a challenging task for the designers of aircraft engines. In order to meet the demands of increased performance and lighter weight, these modern fan designs usually feature low-aspect ratio shroudless rotor blade designs that make the task of achieving adequate flutter margin even more challenging for the aeroelastician. This is especially true for advanced forward swept designs that encompass an entirely new design space compared to previous experience. Fortunately, advances in unsteady computational fluid dynamic (CFD) techniques over the past decade now provide an analysis capability that can be used to quantitatively assess the aeroelastic characteristics of these next generation fans during the design cycle. For aeroelastic applications, Mississippi State University and NASA Glenn Research Center have developed the CFD code TURBO-AE. This code is a time-accurate three-dimensional Euler/Navier-Stokes unsteady flow solver developed for axial-flow turbomachinery that can model multiple blade rows undergoing harmonic oscillations with arbitrary interblade phase angles, i.e., nodal diameter patterns. Details of the code can be found in Chen et al. (1993, 1994), Bakhle et al. (1997, 1998), and Srivastava et al. (1999). To assess aeroelastic stability, the work-per-cycle from TURBO-AE is converted to the critical damping ratio since this value is more physically meaningful, with both the unsteady normal pressure and viscous shear forces included in the work-per-cycle calculation. If the total damping (aerodynamic plus mechanical) is negative, then the blade is unstable since it extracts energy from the flow field over the vibration cycle. TURBO-AE is an integral part of an aeroelastic design system being developed at Honeywell Engines, Systems & Services for flutter and forced response predictions, with test cases from development rig and engine tests being used to validate its predictive capability. A recent experimental program (Sanders et al., 2002) was aimed at providing the necessary unsteady aerodynamic and vibratory response data needed to validate TURBO-AE for fan flutter predictions. A comparison of numerical TURBO-AE simulations with the benchmark flutter data is given in Sanders et al. (2003), with the data used to guide the validation of the code and define best practices for performing accurate unsteady simulations. The agreement between the analyses and the predictions was quite remarkable, demonstrating the ability of the analysis to accurately model the unsteady flow processes driving stall-side flutter.

Sanders, Albert J.↗
NASA NTRS1978

A computational model for the prediction of jet entrainment in the vicinity of nozzle boattails (the BOAT code)

The development of a computational model (BOAT) for calculating nearfield jet entrainment, and its incorporation in an existing methodology for the prediction of nozzle boattail pressures, is discussed. The model accounts for the detailed turbulence and thermochemical processes occurring in the mixing layer formed between a jet exhaust and surrounding external stream while interfacing with the inviscid exhaust and external flowfield regions in an overlaid, interactive manner. The ability of the BOAT model to analyze simple free shear flows is assessed by comparisons with fundamental laboratory data. The overlaid procedure for incorporating variable pressures into BOAT and the entrainment correction employed to yield an effective plume boundary for the inviscid external flow are demonstrated. This is accomplished via application of BOAT in conjunction with the codes comprising the NASA/LRC patched viscous/inviscid methodology for determining nozzle boattail drag for subsonic/transonic external flows.

Dash, S. M.↗
NASA NTRS2005

Stirling Analysis Comparison of Commercial Versus High-Order Methods

Recently, three-dimensional Stirling engine simulations have been accomplished utilizing commercial Computational Fluid Dynamics software. The validations reported can be somewhat inconclusive due to the lack of precise time accurate experimental results from engines, export control/proprietary concerns, and the lack of variation in the methods utilized. The last issue may be addressed by solving the same flow problem with alternate methods. In this work, a comprehensive examination of the methods utilized in the commercial codes is compared with more recently developed high-order methods. Specifically, Lele's compact scheme and Dyson's Ultra Hi-Fi method will be compared with the SIMPLE and PISO methods currently employed in CFD-ACE, FLUENT, CFX, and STAR-CD (all commercial codes which can in theory solve a three-dimensional Stirling model with sliding interfaces and their moving grids limit the effective time accuracy). We will initially look at one-dimensional flows since the current standard practice is to design and optimize Stirling engines with empirically corrected friction and heat transfer coefficients in an overall one-dimensional model. This comparison provides an idea of the range in which commercial CFD software for modeling Stirling engines may be expected to provide accurate results. In addition, this work provides a framework for improving current one-dimensional analysis codes.

Dyson, Rodger W.↗
NASA NTRS2007

Stirling Analysis Comparison of Commercial vs. High-Order Methods

Recently, three-dimensional Stirling engine simulations have been accomplished utilizing commercial Computational Fluid Dynamics software. The validations reported can be somewhat inconclusive due to the lack of precise time accurate experimental results from engines, export control/ proprietary concerns, and the lack of variation in the methods utilized. The last issue may be addressed by solving the same flow problem with alternate methods. In this work, a comprehensive examination of the methods utilized in the commercial codes is compared with more recently developed high-order methods. Specifically, Lele's Compact scheme and Dyson s Ultra Hi-Fi method will be compared with the SIMPLE and PISO methods currently employed in CFD-ACE, FLUENT, CFX, and STAR-CD (all commercial codes which can in theory solve a three-dimensional Stirling model although sliding interfaces and their moving grids limit the effective time accuracy). We will initially look at one-dimensional flows since the current standard practice is to design and optimize Stirling engines with empirically corrected friction and heat transfer coefficients in an overall one-dimensional model. This comparison provides an idea of the range in which commercial CFD software for modeling Stirling engines may be expected to provide accurate results. In addition, this work provides a framework for improving current one-dimensional analysis codes.

Dyson, Rodger W.↗
NASA NTRS1984

Finite Element Prediction of Acoustic Scattering and Radiation from Submerged Elastic Structures

A finite element formulation is derived for the scattering and radiation of acoustic waves from submerged elastic structures. The formulation uses as fundamental unknowns the displacement in the structure and a velocity potential in the field. Symmetric coefficient matrices result. The outer boundary of the fluid region is terminated with an approximate local wave-absorbing boundary condition which assumes that outgoing waves are locally planar. The finite element model is capable of predicting only the near-field acoustic pressures. Far-field sound pressure levels may be determined by integrating the surface pressures and velocities over the wet boundary of the structure using the Helmholtz integral. Comparison of finite element results with analytic results show excellent agreement. The coupled fluid-structure problem may be solved with general purpose finite element codes by using an analogy between the equations of elasticity and the wave equation of linear acoustics.

Everstine, G. C.↗
NASA NTRS1984

Numerical simulation of flow through inlets/diffusers with terminal shocks

Two-dimensional Euler and Navier-Stokes solutions of the flow through three inlet/diffuser configurations with terminal shock systems are reported. Calculations without bleed indicate that the terminal shock location is very sensitive to the outflow back pressure. For cases where there are no available experimental results, it becomes difficult to estimate the back pressure that will result in a terminal shock. Estimates based on quasi-one-dimensional analysis are not found adequate for complex two-dimensional flows. It is found that since the flow downstream of the terminal shock is subsonic, and what happens at the outflow boundary affects the flow inside the inlet, enough of the subsonic diffuser must be modeled to accurately predict the terminal shock region. The diffuser portion should be fairly long with the outflow boundary occurring in a region of more or less uniform flow to be able to prescribe a uniform back pressure. The third configuration studied was investigated with and without incorporating bleed in the code. It is found that the use of bleed stabilizes the shock location and allows solutions which without bleed result in unstarting of the inlet. Comparisons are made with available experimental data.

Talcott, N. A., Jr.↗
NASA NTRS1987

Computational fluid dynamics research at the United Technologies Research Center requiring supercomputers

An overview of research activities at the United Technologies Research Center (UTRC) in the area of Computational Fluid Dynamics (CFD) is presented. The requirement and use of various levels of computers, including supercomputers, for the CFD activities is described. Examples of CFD directed toward applications to helicopters, turbomachinery, heat exchangers, and the National Aerospace Plane are included. Helicopter rotor codes for the prediction of rotor and fuselage flow fields and airloads were developed with emphasis on rotor wake modeling. Airflow and airload predictions and comparisons with experimental data are presented. Examples are presented of recent parabolized Navier-Stokes and full Navier-Stokes solutions for hypersonic shock-wave/boundary layer interaction, and hydrogen/air supersonic combustion. In addition, other examples of CFD efforts in turbomachinery Navier-Stokes methodology and separated flow modeling are presented. A brief discussion of the 3-tier scientific computing environment is also presented, in which the researcher has access to workstations, mid-size computers, and supercomputers.

Landgrebe, Anton J.↗
NASA NTRS1988

Benchmarking and performance analysis of the CM-2

A suite of benchmarking routines testing communication, basic arithmetic operations, and selected kernel algorithms written in LISP and PARIS was developed for the CM-2. Experiment runs are automated via a software framework that sequences individual tests, allowing for unattended overnight operation. Multiple measurements are made and treated statistically to generate well-characterized results from the noisy values given by cm:time. The results obtained provide a comparison with similar, but less extensive, testing done on a CM-1. Tests were chosen to aid the algorithmist in constructing fast, efficient, and correct code on the CM-2, as well as gain insight into what performance criteria are needed when evaluating parallel processing machines.

Myers, David W.↗
NASA NTRS1993

Prediction of rocket plume radiative heating using backward Monte-Carlo method

A backward Monte-Carlo plume radiation code has been developed to predict rocket plume radiative heating to the rocket base region. This paper provides a description of this code and provides sample results. The code was used to predict radiative heating to various locations during test firings of 48-inch solid rocket motors at NASA Marshall Space Flight Center. Comparisons with test measurements are provided. Predictions of full scale sea level Redesigned Solid Rocket Motor (RSRM) and Advanced Solid Rocket Motor (ASRM) plume radiative heating to the Space Shuttle external tank (ET) dome center were also made. A comparison with the Development Flight Instrumentation (DFI) measurements is also provided.

Wang, K. C.↗
NASA NTRS1993

Aerodynamics of Shuttle Orbiter at high altitudes

The high-altitude/high-Knudsen number aerodynamics of the Shuttle Orbiter are computed from Low-Earth Orbit down to 100 km using three-dimensional direct simulation Monte Carlo and free molecule codes. Results are compared with Blanchard's latest Shuttle aerodynamic model, which is based on in-flight accelerometer measurements, and bridging formula models. Good comparison is observed, except for the normal force and pitching moment coefficients. The present results were obtained for a generic Shuttle geometry configuration corresponding to a zero deflection for all control surfaces.

Rault, Didier F. G.↗
NASA NTRS1993

Laser diagnostics for small rockets

Two nonintrusive flowfield diagnostics based on spectrally-resolved elastic (Rayleigh) and inelastic (Raman) laser light scattering were developed for obtaining local flowfield measurements in low-thrust gaseous H2/O2 rocket engines. The objective is to provide an improved understanding of phenomena occurring in small chemical rockets in order to facilitate the development and validation of advanced computational fluid dynamics (CFD) models for analyzing engine performance. The laser Raman scattering diagnostic was developed to measure major polyatomic species number densities and rotational temperatures in the high-density flowfield region extending from the injector through the chamber throat. Initial application of the Raman scattering diagnostic provided O2 number density and rotational temperature measurements in the exit plane of a low area-ratio nozzle and in the combustion chamber of a two-dimensional, optically-accessible rocket engine. In the low-density nozzle exit plane region where the Raman signal is too weak, a Doppler-resolved laser Rayleigh scattering diagnostic was developed to obtain axial and radial mean gas velocities, and in certain cases, H2O translational temperature and number density. The results from these measurements were compared with theoretical predictions from the RPLUS CFD code for analyzing rocket engine performance. Initial conclusions indicate that a detailed and rigorous modeling of the injector is required in order to make direct comparisons between laser diagnostic measurements and CFD predictions at the local level.

Zupanc, Frank J.↗
NASA NTRS1994

Synthesized spectra of turbulent clouds

A comparison has been made of the velocity structure of non-star-forming dense molecular clouds and the output of a 512(exp 3) hydrodynamic code for compressible gas displaying features of turbulence. A remarkable similarity exists between the interstellar line profiles and the velocity profiles found by integrating, with an optically thin assumption, through subsections of the simulation. This similarity is in terms of both the structure of individual spectra, showing skewness, multiple peaks, and wings, and also the variations of the spectra from point to point. It is possible to see the origin of the various velocity substructures and how they correlate with regions of high density and velocity divergence, or regions of high vorticity, by deprojecting individual spectra of the simulation. It is found that the interstellar medium spectra are best represented by the epoch of the simulation in which the energy spectrum is essentially Kolmogorov-like and initially developed shocks have been largely converted into vorticity structures. Further, it is found that the vorticity peaks ascribed to intermittency contribute to the non-Gaussian features and wings of the line profile, but at an intensity level that is weak compared to the contribution of the bulk of the flow.

Falgarone, E.↗
NASA NTRS1998

Deterministic Multiaxial Creep and Creep Rupture Enhancements for CARES/Creep Integrated Design Code

High temperature and long duration applications of monolithic ceramics can place their failure mode in the creep rupture regime. A previous model advanced by the authors described a methodology by which the creep rupture life of a loaded component can be predicted. That model was based on the life fraction damage accumulation rule in association with the modified Monkman-Grant creep rupture criterion. However, that model did not take into account the deteriorating state of the material due to creep damage (e.g., cavitation) as time elapsed. In addition, the material creep parameters used in that life prediction methodology, were based on uniaxial creep curves displaying primary and secondary creep behavior, with no tertiary regime. The objective of this paper is to present a creep life prediction methodology based on a modified form of the Kachanov-Rabotnov continuum damage mechanics (CDM) theory. In this theory, the uniaxial creep rate is described in terms of sum, temperature, time, and the current state of material damage. This scalar damage state parameter is basically an abstract measure of the current state of material damage due to creep deformation. The damage rate is assumed to vary with stress, temperature, time, and the current state of damage itself. Multiaxial creep and creep rupture formulations of the CDM approach are presented in this paper. Parameter estimation methodologies based on nonlinear regression analysis are also described for both, isothermal constant stress states and anisothermal variable stress conditions This creep life prediction methodology was preliminarily added to the integrated design code CARES/Creep (Ceramics Analysis and Reliability Evaluation of Structures/Creep), which is a postprocessor program to commercially available finite element analysis (FEA) packages. Two examples, showing comparisons between experimental and predicted creep lives of ceramic specimens, are used to demonstrate the viability of Ns methodology and the CARES/Creep program.

Jadaan, Osama M.↗
NASA NTRS1998

Deterministic and Probabilistic Creep and Creep Rupture Enhancement to CARES/Creep: Multiaxial Creep Life Prediction of Ceramic Structures Using Continuum Damage Mechanics and the Finite Element Method

High temperature and long duration applications of monolithic ceramics can place their failure mode in the creep rupture regime. A previous model advanced by the authors described a methodology by which the creep rupture life of a loaded component can be predicted. That model was based on the life fraction damage accumulation rule in association with the modified Monkman-Grant creep ripture criterion However, that model did not take into account the deteriorating state of the material due to creep damage (e.g., cavitation) as time elapsed. In addition, the material creep parameters used in that life prediction methodology, were based on uniaxial creep curves displaying primary and secondary creep behavior, with no tertiary regime. The objective of this paper is to present a creep life prediction methodology based on a modified form of the Kachanov-Rabotnov continuum damage mechanics (CDM) theory. In this theory, the uniaxial creep rate is described in terms of stress, temperature, time, and the current state of material damage. This scalar damage state parameter is basically an abstract measure of the current state of material damage due to creep deformation. The damage rate is assumed to vary with stress, temperature, time, and the current state of damage itself. Multiaxial creep and creep rupture formulations of the CDM approach are presented in this paper. Parameter estimation methodologies based on nonlinear regression analysis are also described for both, isothermal constant stress states and anisothermal variable stress conditions This creep life prediction methodology was preliminarily added to the integrated design code CARES/Creep (Ceramics Analysis and Reliability Evaluation of Structures/Creep), which is a postprocessor program to commercially available finite element analysis (FEA) packages. Two examples, showing comparisons between experimental and predicted creep lives of ceramic specimens, are used to demonstrate the viability of this methodology and the CARES/Creep program.

Jadaan, Osama M.↗
NASA NTRS1999

Influence of Ply Waviness on Fatigue Life of Tapered Composite Flexbeam Laminates

Nonlinear tapered flexbeam laminates, with significant ply waviness, were cut from a full-size composite rotor hub flexbeam. The specimens were tested under combined axial tension and cyclic bending loads. All of the specimens had wavy plies through the center and near the surfaces (termed marcelled areas), although for some of the specimens the surface marcels were very obvious, and for others they were much smaller. The specimens failed by first developing cracks through the marcels at the surfaces, and then delaminations grew from those cracks, in both directions. Delamination failure occurred in these specimens at significantly shorter fatigue lives than similar specimens without waviness, tested in ref. 2. A 2D finite element model was developed which closely approximated the flexbeam geometry, boundary conditions, and loading. In addition, the FE model duplicated the waviness observed in one of the test specimens. The model was analyzed using a geometrically nonlinear FE code. Modifications were made to the original model to reduce the amplitude of the marcels near the surfaces. The analysis was repeated for each modification. Comparisons of the interlaminar normal stresses, sigma(sub n), in the various models showed that under combined axial-tension and cyclic-bending loading, for marcels of the same aspect ratio, sigma(sub n) stresses increased as the distance along the taper, from thick to thin end, increased. For marcels of the same aspect ratio and at the same X-location along the taper, sigma(sub n) stresses decreased as the distance from the surface into the flexbeam interior increased. A technique was presented for determining the smallest acceptable marcel aspect ratio at various locations in the flexbeam.

Murri, Gretchen B.↗
NASA NTRS1999

Numerical Analysis of the Trailblazer Inlet Flowfield for Hypersonic Mach Numbers

A study of the Trailblazer vehicle inlet was conducted using the Global Air Sampling Program (GASP) code for flight Mach numbers ranging from 4-12. Both perfect gas and finite rate chemical analysis were performed with the intention of making detailed comparisons between the two results. Inlet performance was assessed using total pressure recovery and kinetic energy efficiency. These assessments were based upon a one-dimensional stream-thrust-average of the axisymmetric flowfield. Flow visualization utilized to examine the detailed shock structures internal to this mixed-compression inlet. Kinetic energy efficiency appeared to be the least sensitive to differences between the perfect gas and finite rate chemistry results. Total pressure recovery appeared to be the most sensitive discriminator between the perfect gas and finite rate chemistry results for flight Mach numbers above Mach 6. Adiabatic wall temperature was consistently overpredicted by the perfect gas model for flight Mach numbers above Mach 4. The predicted shock structures were noticeably different for Mach numbers from 6-12. At Mach 4, the perfect gas and finite rate chemistry models collapse to the same result.

Steffen, C. J., Jr.↗
NASA NTRS1997

A Numerical Simulation of a Normal Sonic Jet into a Hypersonic Cross-Flow

This study involves numerical modeling of a normal sonic jet injection into a hypersonic cross-flow. The numerical code used for simulation is GASP (General Aerodynamic Simulation Program.) First the numerical predictions are compared with well established solutions for compressible laminar flow. Then comparisons are made with non-injection test case measurements of surface pressure distributions. Good agreement with the measurements is observed. Currently comparisons are underway with the injection case. All the experimental data were generated at the Southampton University Light Piston Isentropic Compression Tube.

Jeffries, Damon K.↗