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.
NASTRAN is a large, comprehensive, nonproprietary, general purpose finite element computer code for structural analysis. The Twelfth Users' Colloquim provides some comprehensive papers on the application of finite element methods in engineering, comparisons with other approaches, unique applications, pre and post processing or auxiliary programs, and new methods of analysis with NASTRAN.
An analytical comparison of the predicted combustion flowfields for the Phase 1 and Phase 2+ fuel preburner configurations using an advanced CFD combustion code is currently underway at Rocketdyne. The Phase 2+ injector has a modified element design. The redesigned element provides for improved atomization through a shift in the injector element/face bleed flow split while maintaining an acceptable pressure drop. Cold-flow tests of the new design have been encouraging. As a complement to those tests, the Advanced Rocket Injector Combustor Code (ARICC) was selected to conduct a detailed combustion flowfield simulation. ARICC models two-dimensional (or axisymmetric), transient, turbulent, two-phase mixing and combustion flowfields. It is derived from a member of the Los Alamos ICEd-ALE family of codes. Unique features of ARICC include explicit representation of the coaxial LOX jet and hydrogen gas flows, distributed atomization processes, droplet breakup, and supercritical vaporation processes. Predictions of temperature and OH concentration profiles for the Phase 1 and Phase 2+ injector flowfields indicate a trend toward more uniform temperature distributions and shorter flame lengths with the Phase 2+ design.
The hot section facility at the Lewis Research Center was used to demonstrate the capability of instruments to make required measurements of boundary conditions of the flow field and heat transfer processes in the hostile environment of the turbine. The results of thermal scaling tests show that low temperature and pressure rig tests give optimistic estimates of the thermal performance of a cooling design for high pressure and temperature application. The results of measuring heat transfer coefficients on turbine vane airfoils through dynamic data analysis show good comparison with measurements from steady state heat flux gauges. In addition, the data trends are predicted by the STAN5 boundary layer code. However, the magnitude of the experimental data was not predicted by the analysis, particularly in laminar and transitional regions near the leading edge. The infrared photography system was shown capable of providing detailed surface thermal gradients and secondary flow features on a turbine vane and endwell.
Experimental results were obtained for an Eppler 387 airfoil in the Langley Low Turbulence Pressure Tunnel. The tests were conducted over a Mach number range from 0.03 to 0.13 and a chord Reynolds number range for 60,000 to 460,000. Lift and pitching moment data were obtained from airfoil surface pressure measurements and drag data for wake surveys. Oil flow visualization was used to determine laminar separation and turbulent reattachment locations. Comparisons of these results with data on the Eppler 387 airfoil from two other facilities as well as the Eppler airfoil code are included.
An experimental and computational study has been performed on a close rotor-blade/vortex interaction. Surface pressure data was obtained from a rotor operating close to the tip-vortex from an upstream wing in a wind tunnel. Data was obtained for a wide range of blade-vortex proximities, orientations, and blade-tip Mach numbers (up to the transonic regime). A numerical model of these interactions was constructed using the unsteady, three-dimensional, full-potential rotor code called FPR. The model employed an undistorted full-field representation of the measured vortex. This simple model gave excellent comparisons with the data for a wide range of conditions, including parallel head-on interactions. Computational studies have also been performed on the manner of vortex representation and the influence of vortex-core size.
In this paper a method is described for predicting the compressible, free-wake, flow about a lifting rotor-body configuration. The method is an extension of a unique vorticity embedded full-potential method used to calculate free-wake rotor hover performance. An unusual feature of this method is that it obviates the requirement for multiple grids to treat the rotor-body problem. The approach used to treat the body is similar to that used to include the rotor wake in the full potential calculation. The body is modeled as a structured circulation sheet and the strength of this sheet is determined in an iterative manner. Initially the method is tested to compute the flow past simple isolated bodies like cylinders and spheres. After a comparison of these simple computations with exact solutions this procedure is included into the HELIX-I, free-wake rotor code, to compute the flow around a rotor mounted on a large whirl tower. The effects of the tower on rotor wake geometry and load distribution are presented.
The reflectance of snow-covered surfaces in Minnesota is analyzed using Landsat-5 Thematic Mapper satellite data. Calculations are performed for satellite-derived reflectances integrated over the spectral region (0.45-0.9 micron). Corrections are applied for atmospheric effects and integrated reflectances (R1) are compared over agricultural and forested areas and over a lake using TM scenes acquired in November 1984 and January 1985. Integrated reflectances are then mapped and inter- and intra-scene comparisons of surface reflectance are compared. Temporal analysis of reflectance changes can then be performed rapidly and efficiently using color-coded images. It is noted that the average R1 within the November 1984 subscene was 0.429 + or - 0.176, whereas R1 within the January 1985 subscene was 0.669 + or - 0.236.
A higher-order viscous shock-layer (HVSL) method has been developed and is used to obtain physically consistent results under varying degrees of low-density conditions for perfect-gas and nonequilibrium flows past long slender bodies. This method is implemented with higher-order body and shock slip conditions. Detailed comparisons with the direct simulation Monte Carlo (DSMC) method and Navier-Stokes (NS) calculations establish HVSL as a computationally efficient code for the low-density calculations. Further, it is shown that the deficiencies in the standard viscous shock-layer (VSL) method when predicting low-density flows lie in the governing equations and not with the slip boundary conditions as considered by earlier researchers.
Design of the thermal protection system for any hypersonic flight vehicle requires determination of both the peak temperatures over the surface and the heating-rate history along the flight profile. In this paper, the process used to generate the aerothermal environments required for the X-34 Testbed Technology Demonstrator thermal protection system design is described as it has evolved from a relatively simplistic approach based on engineering methods applied to critical areas to one of detailed analyses over the entire vehicle. A brief description of the trajectory development leading to the selection of the thermal protection system design trajectory is included. Comparisons of engineering heating predictions with wind-tunnel test data and with results obtained using a Navier-Stokes flowfield code and an inviscid/boundary layer method are shown. Good agreement is demonstrated among all these methods for both the ground-test condition and the peak heating flight condition. Finally, the detailed analysis using engineering methods to interpolate the surface-heating-rate results from the inviscid/boundary layer method to predict the required thermal environments is described and results presented.
An aerodynamic database has been generated for the Mars Smart Lander Shelf-All configuration using computational fluid dynamics (CFD) simulations. Three different CFD codes, USM3D and FELISA, based on unstructured grid technology and LAURA, an established and validated structured CFD code, were used. As part of this database development, the results for the Mars continuum were validated with experimental data and comparisons made where applicable. The validation of USM3D and LAURA with the Unitary experimental data, the use of intermediate LAURA check analyses, as well as the validation of FELISA with the Mach 6 CF(sub 4) experimental data provided a higher confidence in the ability for CFD to provide aerodynamic data in order to determine the static trim characteristics for longitudinal stability. The analyses of the noncontinuum regime showed the existence of multiple trim angles of attack that can be unstable or stable trim points. This information is needed to design guidance controller throughout the trajectory.
Detailed aeroheating information is critical to the successful design of a thermal protection system (TPS) for an aerospace vehicle. This report describes NASA Langley Research Center's (LaRC) two-color relative-intensity phosphor thermography method and the IHEAT software package which is used for the efficient data reduction and analysis of the phosphor image data. Development of theory is provided for a new weighted two-color relative-intensity fluorescence theory for quantitatively determining surface temperatures on hypersonic wind tunnel models; an improved application of the one-dimensional conduction theory for use in determining global heating mappings; and extrapolation of wind tunnel data to flight surface temperatures. The phosphor methodology at LaRC is presented including descriptions of phosphor model fabrication, test facilities and phosphor video acquisition systems. A discussion of the calibration procedures, data reduction and data analysis is given. Estimates of the total uncertainties (with a 95% confidence level) associated with the phosphor technique are shown to be approximately 8 to 10 percent in the Langley's 31-Inch Mach 10 Tunnel and 7 to 10 percent in the 20-Inch Mach 6 Tunnel. A comparison with thin-film measurements using two-inch radius hemispheres shows the phosphor data to be within 7 percent of thin-film measurements and to agree even better with predictions via a LATCH computational fluid dynamics solution (CFD). Good agreement between phosphor data and LAURA CFD computations on the forebody of a vertical takeoff/vertical lander configuration at four angles of attack is also shown. In addition, a comparison is given between Mach 6 phosphor data and laminar and turbulent solutions generated using the LAURA, GASP and LATCH CFD codes. Finally, the extrapolation method developed in this report is applied to the X-34 configuration with good agreement between the phosphor extrapolation and LAURA flight surface temperature predictions. The phosphor process outlined in the paper is believed to provide the aerothermodynamic community with a valuable capability for rapidly obtaining (4 to 5 weeks) detailed heating information needed in TPS design.
Characteristics of albedo, or backscatter current, providing a 'background' for calorimeter experiments in high energy cosmic rays are analyzed. The comparison of experimental data obtained in the flights of the ATIC spectrometer is made with simulations performed using the GEANT 3.21 code. The influence of the backscatter on charge resolution in the ATIC experiment is discussed.
During the period 3-10 February 1986 a series of major solar flares occurred on the Sun and several intense geomagnetic storms took place on the Earth. To examine the causality between the solar activity and the geomagnetic activity in this period, a magnetohydrodynamic (MHD) numerical simulation was performed using a 2 1/2 -D numerical code. In that period of February 1986, the Japanese spacecraft Sakigake was at 0.84 AU, 57 deg west of the Earth. Besides the in-situ measurements of the interplanetary plasma, Sakigake also provided Doppler scintillation observations. Comparisons between the results of the MHD simulation and the measurements made by the spacecraft Sakigake are presented.
Numerical optimization was used in conjunction with an inviscid, full potential equation, transonic flow analysis computer code to design an upper surface contour for a conventional airfoil to improve its supercritical performance. The modified airfoil was tested in a compressible flow wind tunnel. The modified airfoil's performance was evaluated by comparison with test data for the baseline airfoil and for an airfoil developed by optimization of leading edge of the baseline airfoil. While the leading edge modification performed as expected, the upper surface re-design did not produce all of the expected performance improvements. Theoretical solutions computed using a full potential, transonic airfoil code corrected for viscosity were compared to experimental data for the baseline airfoil and the upper surface modification. These correlations showed that the theory predicted the aerodynamics of the baseline airfoil fairly well, but failed to accurately compute drag characteristics for the upper surface modification.
In early 1991, the National Aeronautics and Space Administration's (NASA's) Office of Space Flight commissioned the Aeronautics and Space Engineering Board (ASEB) of the National Research Council (NRC) to investigate the adequacy of the current process by which NASA develops and verifies changes and updates to the Space Shuttle flight software. The Committee for Review of Oversight Mechanisms for Space Shuttle Flight Software Processes was convened in Jan. 1992 to accomplish the following tasks: (1) review the entire flight software development process from the initial requirements definition phase to final implementation, including object code build and final machine loading; (2) review and critique NASA's independent verification and validation process and mechanisms, including NASA's established software development and testing standards; (3) determine the acceptability and adequacy of the complete flight software development process, including the embedded validation and verification processes through comparison with (1) generally accepted industry practices, and (2) generally accepted Department of Defense and/or other government practices (comparing NASA's program with organizations and projects having similar volumes of software development, software maturity, complexity, criticality, lines of code, and national standards); (4) consider whether independent verification and validation should continue. An overview of the study, independent verification and validation of critical software, and the Space Shuttle flight software development process are addressed. Findings and recommendations are presented.
Predicted shock-layer emission signatures of the Japanese Hayabusa capsule during its reentry are presented for comparison with flight measurements made during an airborne observation mission using NASA s DC-8 Airborne Laboratory. For each altitude, lines of sight were extracted from flow field solutions computed using an inhouse high-fidelity CFD code, DPLR, at 11 points along the flight trajectory of the capsule. These lines of sight were used as inputs for the line-by-line radiation code NEQAIR, and emission spectra of the air plasma were computed in the wavelength range from 300 nm to 1600 nm, a range which covers all of the different experiments onboard the DC-8. In addition, the computed flow field solutions were post-processed with the material thermal response code FIAT, and the resulting surface temperatures of the heat shield were used to generate thermal emission spectra based on Planck radiation. Both spectra were summed and integrated over the flow field. The resulting emission at each trajectory point was propagated to the DC-8 position and transformed into incident irradiance. Comparisons with experimental data are shown.
An experimental method was developed to determine the droplet impingement characteristics on 2-D and 3-D bodies. The experimental results provide the essential droplet impingement data required to validate water droplet trajectory codes, which are used in the analysis of aircraft icing. A body, whose water droplet impingement characteristics are required, is covered at strategic locations by thin strips of moisture absorbing (blotter) paper, and is exposed to an air stream containing a water dye solution spray cloud. Water droplet impingement data are extracted from the dyed blotter strips by measuring the optical reflectance of the dye deposit on the strips, using an automated reflectometer. Models tested include a 4-inch diameter cylinder, a NACA 652015 airfoil section, a MS(1)-0317 supercritical airfoil section, three simulated ice shapes, an axisymmetric inlet and a Boeing 737-300 inlet model. Detailed descriptions of the dye tracer technique, instrumentation, data reduction method and the results obtained are presented. Analytical predictions of collection efficiency characteristics for most test configurations are included for comparison.
The feasibility of making quantitative nonintrusive NO concentration ([NO]) measurements in nonpremixed flames has been assessed by obtaining laser-induced fluorescence (LIF) measurements of [NO] in counterflow diffusion flames at atmospheric and higher pressures. Comparisons at atmospheric pressure between laser-saturated fluorescence (LSF) and linear LIF measurements in four diluted ethane-air counterflow diffusion flames with strain rates from 5 to 48/s yielded excellent agreement from fuel-lean to moderately fuel-rich conditions, thus indicating the utility of a model-based quenching correction technique, which was then extended to higher pressures. Quantitative LIF measurements of [NO] in three diluted methane-air counterflow diffusion flames with strain rates from 5 to 35/s were compared with OPPDIF model predictions using the GRI (version 2.11) chemical kinetic mechanism. The comparisons revealed that the GRI mechanism underpredicts prompt-NO by 30-50% at atmospheric pressure. Based on these measurements, a modified reaction rate coefficient for the prompt-NO initiation reaction was proposed which causes the predictions to match experimental data. Temperature measurements using thin filament pyrometry (TFP) in conjunction with a new calibration method utilizing a near-adiabatic H2-air Hencken burner gave very good comparisons with model predictions in these counterflow diffusion flames. Quantitative LIF measurements of [NO] were also obtained in four methane-air counterflow partially-premixed flames with fuel-side equivalence ratios (phi(sub B)) of 1.45, 1.6, 1.8 and 2.0. The measurements were in excellent agreement with model predictions when accounting for radiative heat loss. Spatial separation between regions dominated by the prompt and thermal NO mechanisms was observed in the phi(sub B) = 1.45 flame. The modified rate coefficient proposed earlier for the prompt-NO initiation reaction improved agreement between code predictions and measurements in the region where prompt-NO dominates. Finally, LIF measurements of NO were obtained in counterflow diffusion flames at 2 to 5 atm. Comparisons between [NO] measurements and predictions show that the GRI mechanism underpredicts prompt-NO by a factor of two to three at all pressures. In general, the results indicate a need for refinement of the CH chemistry, especially the pressure-dependent CH formation and destruction reactions.