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.
This viewgraph presentation provides information on support sources available for the automatic parallelization of computer program. CAPTools, a support tool developed at the University of Greenwich, transforms, with user guidance, existing sequential Fortran code into parallel message passing code. Comparison routines are then run for debugging purposes, in essence, ensuring that the code transformation was accurate.
I am an intern in the Turbine Branch of the Turbomachinery and Propulsion Systems Division. The division is primarily concerned with experimental and computational methods of calculating heat transfer effects of turbine blades during operation in jet engines and land-based power systems. These include modeling flow in internal cooling passages and film cooling, as well as calculating heat flux and peak temperatures to ensure safe and efficient operation. The branch is research-oriented, emphasizing the development of tools that may be used by gas turbine designers in industry. The branch has been developing a computational fluid dynamics (CFD) and heat transfer code called GlennHT to achieve the computational end of this analysis. The code was originally written in FORTRAN 77 and run on Silicon Graphics machines. However the code has been rewritten and compiled in FORTRAN 90 to take advantage of more modem computer memory systems. In addition the branch has made a switch in system architectures from SGI's to Linux PC's. The newly modified code therefore needs to be tested and validated. This is the primary goal of my internship. To validate the GlennHT code, it must be run using benchmark fluid mechanics and heat transfer test cases, for which there are either analytical solutions or widely accepted experimental data. From the solutions generated by the code, comparisons can be made to the correct solutions to establish the accuracy of the code. To design and create these test cases, there are many steps and programs that must be used. Before a test case can be run, pre-processing steps must be accomplished. These include generating a grid to describe the geometry, using a software package called GridPro. Also various files required by the GlennHT code must be created including a boundary condition file, a file for multi-processor computing, and a file to describe problem and algorithm parameters. A good deal of this internship will be to become familiar with these programs and the structure of the GlennHT code. Additional information is included in the original extended abstract.
A vortex lattice code, CAMRAD II, and a Reynolds-Averaged Navier-Stoke code, OVERFLOW-D2, were used to predict the aerodynamic performance of a two-bladed horizontal axis wind turbine. All computations were compared with experimental data that was collected at the NASA Ames Research Center 80- by 120-Foot Wind Tunnel. Computations were performed for both axial as well as yawed operating conditions. Various stall delay models and dynamics stall models were used by the CAMRAD II code. Comparisons between the experimental data and computed aerodynamic loads show that the OVERFLOW-D2 code can accurately predict the power and spanwise loading of a wind turbine rotor.
A computer code is described which yields accurate solutions for a broad range of laminar, nonsimilar boundary layers, providing the inviscid flow field is known. The boundary layer may be subject to mass injection for perfect-gas, nonreacting flows. If no mass injection is present, the code can be used with either perfect-gas or real-gas thermodynamic models. Solutions, ranging from two-dimensional similarity solutions to solutions for the boundary layer on the Space Shuttle Orbiter during reentry conditions, have been obtained with the code. Comparisons of these solutions, and others, with solutions presented in the literature; and with solutions obtained from other codes, demonstrate the accuracy of the present code.
A code-to-code comparison has been performed for high-fidelity simulations ofNASA’s six-passenger quadrotor air taxi con- cept vehicle. Themultidisciplinary simulations combine comprehensive rotorcraft dynamics with high-fidelity fluid dynamics obtained from an unsteady Navier-Stokes computational fluid dynamics code. The results are then compared with those in the literature. A quantitative comparison of converged trim angles has been performed and normal force, chord force, and pitching-moment coefficients arepresented for qualitativecomparison.
Four experiments utilizing passive detectors (P0006, P0004, A0015, M0004) were flown on LDEF to study the radiation environment. These experiments have been summarized in a companion paper (Benton et al., 1996). One of the experimental goals was to measure LET spectra at different locations and shielding depths with plastic nuclear track detectors (PNTD). It was found that the LET spectra extended well above the LET cutoff imposed by the geomagnetic field on GCR particle penetration into LEO. The high LET particles detected were mostly short-range (range < 2000 m), indicating that they were secondaries produced locally within the PNTD. The presence of these high LET particle fluences is important for the determination of dose equivalent because of the high Quality Factors (Q) involved. A relatively small fraction of particle fluence can contribute a large fraction of dose equivalent. Short-range, inelastic secondary particles produced by trapped protons in the South Atlantic Anomaly (SAA) were found to be a major contributor to the LET spectra above 100 keV/micrometer. The LET spectra were found to extend beyond the approximately 137 keV/micrometer relativistic GCR Fe peak to over 1000 keV/micrometer. The high LET tail of the LET spectra was measured in CR-39 and polycarbonate PNTDs using different techniques. GCR made a relatively modest contribution to the LET spectra as compared to the contributions from short-range secondary particles and stopping protons. LET spectra intercomparisons were made between LDEF measurements and exposures to 154 MeV accelerated proton beams. The similarities support the role of nuclear interactions by trapped protons as the major source of secondary particles in the PNTDs. Also techniques were employed to reduce the range cutoff for detection of the short-range secondaries to approximately 1 micrometer, so that essentially all secondary particles were included in the LET spectra. This has allowed a more realistic assessment of secondary contribution to dose equivalent. Comparisons of measured and calculated LET spectra have been made that demonstrate the need for more accurate modeling of secondary particles in radiation transport codes. Comparisons include preliminary calculations in which attempts have been made to include secondary particles.
Examined is sensitivity of separation extent, wall pressure and heating to variation of primary input flow parameters, such as Mach and Reynolds numbers and shock strength, for 2D and Axisymmetric Hypersonic Shock Wave Turbulent Boundary Layer interactions obtained by Navier-Stokes methods using the SST turbulence model. Baseline parametric sensitivity response is provided in part by comparison with vetted experiments, and in part through updated correlations based on free interaction theory concepts. A recent database compilation of hypersonic 2D shock-wave/turbulent boundary layer experiments extensively used in a prior related uncertainty analysis provides the foundation for this updated correlation approach, as well as for more conventional validation. The primary CFD method for this work is DPLR, one of NASA's real-gas aerothermodynamic production RANS codes. Comparisons are also made with CFL3D, one of NASA's mature perfect-gas RANS codes. Deficiencies in predicted separation response of RANS/SST solutions to parametric variations of test conditions are summarized, along with recommendations as to future turbulence approach.
Predicted XV-15 exterior surface acoustic pressures are compared with previously published experimental data. Surface acoustic pressure transducers were concentrated near the tip-path-plane of the rotor in airplane mode. The comparison emphasized cruise conditions which are of interest for tiltrotor interior noise - level flight for speeds ranging from 72 m/s to 113 m/s. The predictions were produced by components of the NASA Langley Tiltrotor Aeroacoustic Code (TRAC) system of computer codes. Comparisons between measurements and predictions were made in both the time and frequency domains, as well as overall sound pressure levels. In general, the predictions replicated the measured data well. Discrepancies between measurements and predictions were noted. Some of the discrepancies were due to poor correlation of the measured data with the rotor tach signal. In other cases limitations of the predictive methodology have been indicated.
A comparison of the NPARC, PAB, and WIND (previously known as NASTD) Navier-Stokes solvers is made for two flow cases with turbulent mixing as the dominant flow characteristic, a two-dimensional ejector nozzle and a Mach 1.5 elliptic jet. The objective of the work is to determine if comparable predictions of nozzle flows can be obtained from different Navier-Stokes codes employed in a multiple site research program. A single computational grid was constructed for each of the two flows and used for all of the Navier-Stokes solvers. In addition, similar k-e based turbulence models were employed in each code, and boundary conditions were specified as similarly as possible across the codes. Comparisons of mass flow rates, velocity profiles, and turbulence model quantities are made between the computations and experimental data. The computational cost of obtaining converged solutions with each of the codes is also documented. Results indicate that all of the codes provided similar predictions for the two nozzle flows. Agreement of the Navier-Stokes calculations with experimental data was good for the ejector nozzle. However, for the Mach 1.5 elliptic jet, the calculations were unable to accurately capture the development of the three dimensional elliptic mixing layer.
Thermal analysis in both simple and complex models can require calculating the propagation of radiant energy to and from multiple surfaces. This can be accomplished through simple estimation techniques or complex computationally intense computer modeling simulations. Currently there are a variety of computer analysis techniques used to simulate the propagation of radiant energy, each having advantages and disadvantages. The major objective of this effort was to compare two ray tracing radiation propagation analysis programs (NEVADA and TSS) Net Energy Verification and Determination Analyzer and Thermal Synthesizer System with experimental data. Results from a non-flowing, electrically heated test rig was used to verify the calculated radiant energy propagation from a nozzle geometry that represents an aircraft propulsion nozzle system. In general the programs produced comparable overall results, and results slightly higher then the experimental data. Upon inspection of individual radiation interchange factors, differences were evident and would have been magnified if a more radical model temperature profile was analyzed. Bidirectional reflectivity data (BRDF) was not used do to modeling limitations in TSS. For code comparison purposes, this nozzle geometry represents only one case for one set of analysis conditions. Since each computer code has advantages and disadvantages bases on scope, requirements, and desired accuracy, the usefulness of this single case study may be limiting.
Experimental longitudinal and lateral-directional aerodynamic characteristics were obtained for the Pegasus and Pegasus XL configurations over a Mach number range from 1.6 to 6 and angles of attack from -4 to +24 degrees. Angle of sideslip was varied from -6 to +6 degrees, and control surfaces were deflected to obtain elevon, aileron, and rudder effectiveness. Experimental data for the Pegasus configuration are compared with engineering code predictions performed by Nielsen Engineering & Research, Inc. (NEAR) in the aerodynamic design of the Pegasus vehicle, and with results from the Aerodynamic Preliminary Analysis System (APAS) code. Comparisons of experimental results are also made with longitudinal flight data from Flight #2 of the Pegasus vehicle. Results show that the longitudinal aerodynamic characteristics of the Pegasus and Pegasus XL configurations are similar, having the same lift-curve slope and drag levels across the Mach number range. Both configurations are longitudinally stable, with stability decreasing towards neutral levels as Mach number increases. Directional stability is negative at moderate to high angles of attack due to separated flow over the vertical tail. Dihedral effect is positive for both configurations, but is reduced 30-50 percent for the Pegasus XL configuration because of the horizontal tail anhedral. Predicted longitudinal characteristics and both longitudinal and lateral-directional control effectiveness are generally in good agreement with experiment. Due to the complex leeside flowfield, lateral-directional characteristics are not as well predicted by the engineering codes. Experiment and flight data are in good agreement across the Mach number range.
Finite difference methods for unsteady transonic flow frequency use simplified equations in which certain of the time dependent terms are omitted from the governing equations. Kernel functions are derived for two dimensional subsonic flow, and provide accurate solutions of the linearized potential equation with the same time dependent terms omitted. These solutions make possible a direct evaluation of the finite difference codes for the linear problem. Calculations with two of these low frequency kernel functions verify the accuracy of the LTRAN2 and HYTRAN2 finite difference codes. Comparisons of the low frequency kernel function results with the Possio kernel function solution of the complete linear equations indicate the adequacy of the HYTRAN approximation for frequencies in the range of interest for flutter calculations.
The goal of the work is to use multiple codes and multiple configurations to provide an assessment of the capability of RANS solvers to predict circulation control dual radius airfoil performance and also to identify key issues associated with the computational predictions of these configurations that can result in discrepancies in the predicted solutions. Solutions were obtained for the Georgia Tech Research Institute (GTRI) dual radius circulation control airfoil and the General Aviation Circulation Control (GACC) dual radius airfoil. For the GTRI-DR airfoil, two-dimensional structured and unstructured grid computations predicted the experimental trend in sectional lift variation with blowing coefficient very well. Good code to code comparisons between the chordwise surface pressure coefficients and the solution streamtraces also indicated that the detailed flow characteristics were matched between the computations. For the GACC-DR airfoil, two-dimensional structured and unstructured grid computations predicted the sectional lift and chordwise pressure distributions accurately at the no blowing condition. However at a moderate blowing coefficient, although the code to code variation was small, the differences between the computations and experiment were significant. Computations were made to investigate the sensitivity of the sectional lift and pressure distributions to some of the experimental and computational parameters, but none of these could entirely account for the differences in the experimental and computational results. Thus, CFD may indeed be adequate as a prediction tool for dual radius CC flows, but limited and difficult to obtain two-dimensional experimental data prevents a confident assessment at this time.
Experimental data were obtained on an optimally contoured nozzle with an area ratio of 1025:1 and on a truncated version of this nozzle with an area ratio of 440:1. The nozzles were tested with gaseous hydrogen and liquid oxygen propellants at combustion chamber pressures of 1800 to 2400 psia and mixture ratios of 3.89 to 6.15. This report compares the experimental performance, heat transfer, and boundary layer total pressure measurements with theoretical predictions of the current Joint Army, Navy, NASA, Air Force (JANNAF) developed methodology. This methodology makes use of the Two-Dimensional Kinetics (TDK) nozzle performance code. Comparisons of the TDK-predicted performance to experimentally attained thrust performance indicated that both the vacuum thrust coefficient and the vacuum specific impulse values were approximately 2.0-percent higher than the turbulent prediction for the 1025:1 configurations, and approximately 0.25-percent higher than the turbulent prediction for the 440:1 configuration. Nozzle wall temperatures were measured on the outside of a thin-walled heat sink nozzle during the test fittings. Nozzle heat fluxes were calculated front the time histories of these temperatures and compared with predictions made with the TDK code. The heat flux values were overpredicted for all cases. The results range from nearly 100 percent at an area ratio of 50 to only approximately 3 percent at an area ratio of 975. Values of the integral of the heat flux as a function of nozzle surface area were also calculated. Comparisons of the experiment with analyses of the heat flux and the heat rate per axial length also show that the experimental values were lower than the predicted value. Three boundary layer rakes mounted on the nozzle exit were used for boundary layer measurements. This arrangement allowed total pressure measurements to be obtained at 14 different distances from the nozzle wall. A comparison of boundary layer total pressure profiles and analytical predictions show good agreement for the first 0.5 in. from the nozzle wall; but the further into the core flow that measurements were taken, the more that TDK overpredicted the boundary layer thickness.
Spectrally and spatially resolved radiance has been measured in the Electric Arc Shock Tube (EAST) facility, with the aim of improving fundamental understanding of high enthalpy flows in pure nitrogen. These tests provide data to inform models used for simulations of high speed flight in nitrogen rich atmospheres, such as Earth or Titan. The experiments presented in this paper cover conditions from approximately 6 km/s to 11 km/s at an initial pressure of 0.2 Torr. A wide range of physics, with different degrees of non-equilibrium and nitrogen dissociation, are covered. The EAST data are presented in different formats for analysis and comparisons. These formats include the spectral radiance at equilibrium (where appropriate), the spatial dependence of radiance over defined wavelength ranges and the mean non-equilibrium spectral radiance (the so-called "spectral non-equilibrium metric"). All the information needed to simulate each experimental trace, including free-stream conditions, shock time of arrival (i.e. x-t) relation, and the spectral and spatial resolution functions, are provided. Equilibrium radiation calculations are shown as a reference. It is the intention of this paper to motivate code comparisons benchmarked against this data set.
As NASA pursues the exploration mission, advanced propulsion for the next generation of spacecraft will be needed. These new propulsion systems will require higher performance and increased durability, despite current limitations on materials. A break-through technology is needed in the thrust chamber. In this paper the idea of using a porous metallic foam is examined for its potential cooling enhancement capabilities. The goal is to increase the chamber wall cooling without creating an additional pressure drop penalty. A feasibility study based on experiments at laboratory-scale conditions was performed and analysis at rocket conditions is underway. In the experiment, heat transfer and pressure drop data were collected using air as the coolant in a copper or nickel foam filled annular channel. The foam-channel performance was evaluated based on comparison with conventional microchannel cooling passages under equal pressure drop conditions. The heat transfer enhancement of the foam channel over the microchannel ranges from 130% to 172%. The enhancement is relatively independent of the pressure drop and increases with decreasing pore size. A direct numerical simulation model of the foam heat exchange has been built. The model is based on the actual metal foam microstructure of thin ligaments (0.2- 0.3 mm in diameter) that form a network of interconnected open-cells. The cell dimension is around 2 mm. The numerical model was built using the FLUENT CFD code. Comparison of the pressure drop results predicted by the current model with those experimental data of Leong and Jin [8] shows favorable comparisons. Pressure drop predictions have been made using hydrogen as a coolant at typical rocket conditions. Conjugate heat transfer analysis using the foam filled channel is planned for the future.
A second international AIAA Drag Prediction Workshop (DPW-II) was organized and held in Orlando Florida on June 21-22, 2003. The primary purpose was to inves- tigate the code-to-code uncertainty. address the sensitivity of the drag prediction to grid size and quantify the uncertainty in predicting nacelle/pylon drag increments at a transonic cruise condition. This paper presents an in-depth analysis of the DPW-II computational results from three state-of-the-art unstructured grid Navier-Stokes flow solvers exercised on similar families of tetrahedral grids. The flow solvers are USM3D - a tetrahedral cell-centered upwind solver. FUN3D - a tetrahedral node-centered upwind solver, and NSU3D - a general element node-centered central-differenced solver. For the wingbody, the total drag predicted for a constant-lift transonic cruise condition showed a decrease in code-to-code variation with grid refinement as expected. For the same flight condition, the wing/body/nacelle/pylon total drag and the nacelle/pylon drag increment predicted showed an increase in code-to-code variation with grid refinement. Although the range in total drag for the wingbody fine grids was only 5 counts, a code-to-code comparison of surface pressures and surface restricted streamlines indicated that the three solvers were not all converging to the same flow solutions- different shock locations and separation patterns were evident. Similarly, the wing/body/nacelle/pylon solutions did not appear to be converging to the same flow solutions. Overall, grid refinement did not consistently improve the correlation with experimental data for either the wingbody or the wing/body/nacelle pylon configuration. Although the absolute values of total drag predicted by two of the solvers for the medium and fine grids did not compare well with the experiment, the incremental drag predictions were within plus or minus 3 counts of the experimental data. The correlation with experimental incremental drag was not significantly changed by specifying transition. Although the sources of code-to-code variation in force and moment predictions for the three unstructured grid codes have not yet been identified, the current study reinforces the necessity of applying multiple codes to the same application to assess uncertainty.
Detailed spectrally and spatially resolved radiance has been measured in the Electric Arc Shock Tube at NASA Ames Research Center for conditions relevant to Titan entry, with varying atmospheric composition, free-stream density (equivalently, altitude) and shock velocity. The test campaign measured radiation at velocities from 4.7 km/s to 8 km/s and free-stream pressures of 0.1, 0.28 and 0.47 Torr with a variety of compositions. Radiances measured in this work are substantially larger compared to that reported both in past EAST test campaigns and in other shock tube facilities. Depending on the metric used for comparison, the discrepancy can be as high as an order of magnitude. Due to the difference with previously reported data, a substantial effort was undertaken to provide confidence in the new results. The present work provides a new benchmark set of data to replace those published in previous studies. The effect of gas impurities identified in previous shock tube studies was also examined by testing in pure N2 and deliberate addition of air to the CH4/N2 mixtures. Furthermore, a test campaign in pure N2 was also conducted with the aim of providing data for improving fundamental understanding of high enthalpy flows containing N2, such as high-speed entries into Earth or Titan. These experiments cover conditions from approximately 6 km/s to 11 km/s at an initial pressure of 0.2 Torr. It is the intention of this paper to motivate code comparisons benchmarked against this data set.