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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
NASA NTRS1986

New integration techniques for chemical kinetic rate equations. I - Efficiency comparison

A comparison of the efficiency of several recently developed numerical techniques for solving chemical kinetic rate equations is presented. The solution procedures examined include two general-purpose codes, EPISODE and LSODE, developed as multipurpose differential equation solvers, and three specialzed codes, CHEMEQ, CREK1D, and GCKP84, developed specifically for chemical kinetics. The efficiency comparison is made by applying these codes to two practical combustion kinetics problems. Both problems describe adiabatic, constant-pressure, gas-phase chemical reactions and include all three combustion regimes: induction, heat release, and equilibration. The comparison shows that LSODE is the fastest routine currently available for solving chemical kinetic rate equations. An important finding is that an iterative solution of the algebraic enthalpy conservation equation for temperature can be significantly faster than evaluation of the temperature by integration of its time derivative. Significant increases in computational speed are realized by updating the reaction rate constants only when the temperature change exceeds an amount Delta-T that is problem dependent. An approximate expression for the automatic evaluation of Delta-T is presented and is shown to result in increased computational speed.

Radhakrishnan, K.↗
NASA NTRS2005

CFD Code Validation of Wall Heat Fluxes for a G02/GH2 Single Element Combustor

This paper puts forth the case for the need for improved injector design tools to meet NASA s Vision for Space Exploration goals. Requirements for this improved tool are outlined and discussed. The potential for Computational Fluid Dynamics (CFD) to meet these requirements is noted along with its current shortcomings, especially relative to demonstrated solution accuracy. The concept of verification and validation is introduced as the primary process for building and quantifying the confidence necessary for CFD to be useful as an injector design tool. The verification and validation process is considered in the context of the Marshall Space Flight Center (MSFC) Combustion Devices CFD Simulation Capability Roadmap via the Simulation Readiness Level (SRL) concept. The portion of the validation process which demonstrates the ability of a CFD code to simulate heat fluxes to a rocket engine combustor wall is the focus of the current effort. The FDNS and Loci-CHEM codes are used to simulate a shear coaxial single element G02/GH2 injector experiment. The experiment was conducted a t a chamber pressure of 750 psia using hot propellants from preburners. A measured wall temperature profile is used as a boundary condition to facilitate the calculations. Converged solutions, obtained from both codes by using wall functions with the K-E turbulence model and integrating to the wall using Mentor s baseline turbulence model, are compared to the experimental data. The initial solutions from both codes revealed significant issues with the wall function implementation associated with the recirculation zone between the shear coaxial jet and the chamber wall. The FDNS solution with a corrected implementation shows marked improvement in overall character and level of comparison to the data. With the FDNS code, integrating to the wall with Mentor s baseline turbulence model actually produce a degraded solution when compared to the wall function solution with the K--E model. The Loci-CHEM solution, produced by integrating to the wall with Mentor s baseline turbulence model, matches both the heat flux rise rate in the near injector region and the peak heat flux level very well. However, it moderately over predicts the heat fluxes downstream of the reattachment point. The Loci-CHEM solution achieved by integrating to the wall with Mentor s baseline turbulence model was clearly superior to the other solutions produced in this effort.

Lin, Jeff↗
NASA NTRS1980

Drag reduction effects in turbulent boundary layers over wavy walls

Two dimensional incompressible flow over wavy surfaces are analyzed numerically by spectral methods. Algorithms for periodic flows (Fourier modes in the periodic flow direction and Chebycheff modes in the normal direction), and inflow-outflow boundary conditions (Chebycheff modes used in both directions) are described. Results obtained using both codes are reported for laminar flows. Comparisons with known theoretical and experimental results are made.

Balasubramanian, R.↗
NASA NTRS1983

Design of the stift time and frequency transfer microwave ground terminal

The Satellite Time and Frequency Transfer System is intended to provide, simultaneously, global time comparisons at the subnanosecond level and frequency comparisons to better than 1 part in 10 to the 14th power. It utilizes an orbiting hydrogen maser clock and frequency standard that communicates, via microwave links, time and frequency information to Earth terminals operated by hydrogen masers controlling local clocks. A two-way microwave link, to and from the space vehicle provides Doppler information used to cancel the Doppler shifts in a one-way link from the spaceborne oscillator. Pseudo-random noise (PRN) modulation in the two-way link also provides range information to cancel the range delay in the PRN time transfer between space and Earth. The pseudo-random noise modulation system for time difference measurement and its incorporation in the Doppler cancellation system for frequency comparison is explained. The particular PRN code sequence selected and an analysis of the system is discussed.

Penfield, H.↗
NASA NTRS1987

A comparison of five benchmarks

Five benchmark programs were obtained and run on the NASA Lewis CRAY X-MP/24. A comparison was made between the programs codes and between the methods for calculating performance figures. Several multitasking jobs were run to gain experience in how parallel performance is measured.

Huss, Janice E.↗
NASA NTRS1988

An initial investigation into methods of computing transonic aerodynamics sensitivity coefficients

Based upon the program outlined in the modified proposal, the initial effort has concentrated on developing the quasi-analytical approach for determining aerodynamic influence coefficients. In addition, in order to keep these initial efforts computationally efficient and straightforward, it was decided to first consider only two-dimensional flow and to model the problem using the transonic small perturbation equation with surface boundary conditions applied along the slit. The small perturbation equation was selected instead of the full potential equation because it permits the capturing of important nonlinear physical phenomena while being by comparison easy to handle from a coding standpoint. Consequently, the initial objectives were: (1) investigate the feasibility of using the quasi-analytical approach for determining aerodynamic sensitivity coefficients in the transonic regime; (2) determine the accuracy of the quasi-analytical approach by comparing results to those computed using finite difference techniques; and (3) investigate efficient methods for solving the large set of equations associated with the quasi-analytical approach.

Carlson, Leland A.↗
NASA NTRS1991

Computation of atmospheric cooling rates by exact and approximate methods

Infrared fluxes and cooling rates for several standard model atmospheres, with and without water vapor, carbon dioxide, and ozone, have been calculated using a line-by-line method at 0.01/cm resolution. The sensitivity of the results to the vertical integration scheme and to the model for water vapor continuum absorption is shown. Comparison with similar calculations performed at NOAA/GFDL shows agreement to within 0.5 W/sq m in fluxes at various levels and 0.05 K/d in cooling rates. Comparison with a fast, parameterized radiation code used in climate models reveals a worst case difference, when all gases are included, of 3.7 W/sq m in flux; cooling rate differences are 0.1 K/d or less when integrated over a substantial layer with point differences as large as 0.3 K/d.

Ridgway, William L.↗
NASA NTRS1992

Simulation of real-gas effects on pressure distributions for aeroassist flight experiment vehicle and comparison with prediction

Pressure distributions measured on a 60 degree half-angle elliptic cone, raked off at an angle of 73 degrees from the cone centerline and having an ellipsoid nose (ellipticity equal to 2.0 in the symmetry plane) are presented for angles of attack from -10 degrees to 10 degrees. The high normal shock density ratio aspect of a real gas was simulated by testing in Mach 6 air and CF sub 4 (density ratio equal to 5.25 and 12.0, respectively). The effects of Reynolds number, angle of attack, and normal shock density ratio on these measurements are examined, and comparisons with a three dimensional Euler code known as HALIS are made. A significant effect of density ratio on pressure distributions on the cone section of the configuration was observed; the magnitude of this effect decreased with increasing angle of attack. The effect of Reynolds number on pressure distributions was negligible for forebody pressure distributions, but a measurable effect was noted on base pressures. In general, the HALIS code accurately predicted the measured pressure distributions in air and CF sub 4.

Micol, John R.↗
NASA NTRS1992

Advanced ballistic range technology

Experimental interferograms, schlieren, and shadow-graphs are routinely used for validating ideal- and real-gas computational fluid dynamics (CFD) codes. Interferograms are used for quantitative comparisons of integrated density. Schlieren and shadowgraphs are used for comparing the locations of flow field features; the types of flowfield features recorded in these images include shocks, shear layers, and expansion fans. By constructing images from flow field simulations, one-to-one comparisons of computation and experiment are possible. During the period from 1 May 1991 to 31 Jan. 1992, the emphasis of this research was on developing techniques for constructing these images from flow field simulations.

Yates, Leslie A.↗
NASA NTRS1992

Application of computational fluid dynamics to the design of the Space Transportation Main Engine subscale nozzle

CFD analyses of the Space Transportation Main Engine film/dump cooled subscale nozzle are presented, with an emphasis on the timely impact of CFD in the design of the subscale nozzle secondary coolant system. Calculations were performed with the Generalized Aerodynamic Simulation Program (GASP), using a Baldwin-Lomas Turbulence model, and finite rate hydrogen-oxygen chemistry. Design iterations for both the secondary coolant cavity passage and the secondary coolant lip are presented. In addition, validation of the GASP chemistry and turbulence models by comparison with data and other CFD codes are presented for a hypersonic laminar separation corner, a backward facing step, and a 2D scramjet nozzle with hydrogen-oxygen kinetics.

Garrett, J. L.↗
NASA NTRS1993

A residual flexibility approach to multibody dynamics

Many complex systems can be modeled as a collection of interacting bodies, where the relative motion of the bodies may be large. The dynamics of such systems are simulated using multibody dynamic formulations. Many of these treat each body as a rigid component, but recently the flexibility of the components has been incorporated. This paper presents a residual flexibility formulation of the multibody dynamics problem. The formulation is very simple and offers great computational efficiency since it treats each body as a free structure in space, interacting with other bodies only through interface forces. Each body's accelerations can be solved independently, as can each set of interface forces. We have applied the technique successfully to several special applications, and the initial implementation in a general mechanisms code has given excellent results in comparison to a direct finite element representation of flexibility.

Blelloch, Paul A.↗
NASA NTRS1994

Proof test of the computer program BUCKY for plasticity problems

A theoretical equation describing the elastic-plastic deformation of a cantilever beam subject to a constant pressure is developed. The theoretical result is compared numerically to the computer program BUCKY for the case of an elastic-perfectly plastic specimen. It is shown that the theoretical and numerical results compare favorably in the plastic range. Comparisons are made to another research code to further validate the BUCKY results. This paper serves as a quality test for the computer program BUCKY developed at NASA Johnson Space Center.

Smith, James P.↗
NASA NTRS1995

High-resolution measurements of the K-alpha spectra of low-ionizationm species of iron: A new spectral signature of nonequilibrium ionization conditions in young supernova remnants

We present the first systematic laboratory measurements of high-resolution K-alpha spectra of intermediate ions of iron, Fe X-XVII. These lines are not produced in collisional equilibrium plasmas because of the relevant charge states cannot exist at the high electron temperatures required for appreciable excitation of the K-alpha transitions. However, they can provide excellent spectral diagnostics for nonequilibrium ionization conditions, such the ionizing plasmas of young supernova remnants. To facilitate the line identifications, we compare our spectra with theoretical atomic calculations performed using multiconfiguration parametric potential and Dirac-Fock atomic codes. Our measurements also allow direct comparison with time-dependent ionization balance calculations for ionizing plasmas, and good agreement is found.

Decaux, V.↗
NASA NTRS2002

Nonlinear Pressurization and Modal Analysis Procedure for Dynamic Modeling of Inflatable Structures

An introduction and set of guidelines for finite element dynamic modeling of nonrigidized inflatable structures is provided. A two-step approach is presented, involving 1) nonlinear static pressurization of the structure and updating of the stiffness matrix and 2) hear normal modes analysis using the updated stiffness. Advantages of this approach are that it provides physical realism in modeling of pressure stiffening, and it maintains the analytical convenience of a standard bear eigensolution once the stiffness has been modified. Demonstration of the approach is accomplished through the creation and test verification of an inflated cylinder model using a large commercial finite element code. Good frequency and mode shape comparisons are obtained with test data and previous modeling efforts, verifying the accuracy of the technique. Problems encountered in the application of the approach, as well as their solutions, are discussed in detail.

Smalley, Kurt B.↗
NASA NTRS2004

Fundamental Investigations of Airframe Noise

An extensive numerical and experimental study of airframe noise mechanisms associated with a subsonic high-lift system has been performed at NASA Langley Research Center (LaRC). Investigations involving both steady and unsteady computations and experiments on a small-scale, part-span flap model are presented. Both surface (steady and unsteady pressure measurements, hot films, oil flows, pressure sensitive paint) and off surface (5 hole-probe, particle-imaged velocimetry, laser velocimetry, laser light sheet measurements) were taken in the LaRC Quiet Flow Facility (QFF) and several hard-wall tunnels up to flight Reynolds number. Successful microphone array measurements were also taken providing both acoustic source maps on the model, and quantitative spectra. Critical directivity measurements were obtained in the QFF. NASA Langley unstructured and structured Reynolds- Averaged Navier-Stokes codes modeled the flap geometries excellent comparisons with surface and offsurface experimental data were obtained. Subsequently, these meanflow calculations were utilized in both linear stability and direct numerical simulations of the flap-edge flow field to calculate unsteady surface pressures and farfield acoustic spectra. Accurate calculations were critical in obtaining not only noise source characteristics, but shear layer correction data as well. Techniques utilized in these investigations as well as brief overviews of results will be given.

Macaraeg, M. G.↗
NASA NTRS1992

A High Angle of Attack Inviscid Shuttle Orbiter Computation

An accurate description of the aerothermal environment is required to minimize the weight of the Thermal Protection System required on the leeside of winged reentry vehicles. The inability of ground-test facilities to reproduce the high enthalpy, separated flow present during reentry flight conditions, coupled with the prohibitive expense of flight tests, leads to the use of an analytical method - namely Computational Fluid Dynamics (CFD) - to describe the flow. While the ultimate goal of this work is to accurately predict the leeside flow and its associated thermal environment, an essential and reasonable first step towards that goal is a comparison of pressure predictions by a code with wind-tunnel data. Until such CFD pressure predictions agree with wind-tunnel test cases, there is little hope of accurately predicting the thermal environment at flight conditions. Thus, the objective of this study is to compare the pressure distributions predicted by inviscid, perfect gas CFD to Shuttle Orbiter wind-tunnel data and to address any significant issues encountered during the computation. While flight data is available for the Shuttle Orbiter, a wind-tunnel case is chosen for this study to allow a tractable problem for preliminary investigation. A wind-tunnel case allows the perfect gas assumption for the flow chemistry. This provides a significant computational savings over a several species finite-rate chemistry model which would be necessary if high-temperature effects present at flight conditions were to be included. In addition, by concentrating on the surface pressures, the analysis need only consider inviscid flow for general evaluation of the code capability. This further reduces the computational expense due to the absence of viscous terms and the associated decrease in the number of points required for the computational grid. Previous computational efforts (such as STEIN and HALIS) have been directed toward the windward surface quantities, primarily due to restrictions in treating either the winged geometry or its associated subsonic regions at high angle of attack. The code used for this study, the LAURA (Langley Aerothermodynamic Upwind Relaxation Algorithm) code of Gnoffo, represents a state-of-the-art code for computing the flow over complex configurations at hypersonic speeds. In the study, the LAURA code is applied to a wind-tunnel condition to initiate the assessment of the code's ability to predict the flow over a relatively complex hypersonic vehicle at high angles of attack. This presentation is used to highlight the pertinent results of a more detailed investigation of the inviscid calculation over the Shuttle Orbiter with the LAURA code.

William L. Kleb↗
NASA NTRS2008

Vibration Response Models of a Stiffened Aluminum Plate Excited by a Shaker

Numerical models of structural-acoustic interactions are of interest to aircraft designers and the space program. This paper describes a comparison between two energy finite element codes, a statistical energy analysis code, a structural finite element code, and the experimentally measured response of a stiffened aluminum plate excited by a shaker. Different methods for modeling the stiffeners and the power input from the shaker are discussed. The results show that the energy codes (energy finite element and statistical energy analysis) accurately predicted the measured mean square velocity of the plate. In addition, predictions from an energy finite element code had the best spatial correlation with measured velocities. However, predictions from a considerably simpler, single subsystem, statistical energy analysis model also correlated well with the spatial velocity distribution. The results highlight a need for further work to understand the relationship between modeling assumptions and the prediction results.

Cabell, Randolph H.↗
NASA NTRS2008

Time-Dependent Erosion of Ion Optics

The accurate prediction of thruster life requires time-dependent erosion estimates for the ion optics assembly. Such information is critical to end-of-life mechanisms such as electron backstreaming. CEX2D was recently modified to handle time-dependent erosion, double ions, and multiple throttle conditions in a single run. The modified code is called "CEX2D-t". Comparisons of CEX2D-t results with LDT and ELT post-tests results show good agreement for both screen and accel grid erosion including important erosion features such as chamfering of the downstream end of the accel grid and reduced rate of accel grid aperture enlargement with time.

ion thrusters↗