ALBERT

Description: This paper describes the computational work performed on the simulation of a 16-in shock-tunnel facility. The numerical problems encountered during the computation of these flows are discussed along with the validity of some approximations used, notably concerning the reduction of the problem into problems of smaller dimensionality. Quasi-1D simulations can be used to help design experiments, or to better understanding the characteristics of the facility. An application to the design of a nonintrusive diagnostic is shown. The multidimensional flow transients computed include the shock reflection at the end of the driven tube, the shock propagation down the nozzle, and the breaking of the main diaphragm.

Description: A new upwind, parabolized Navier-Stokes (PNS) code has been developed to compute the hypersonic, viscous, chemically reacting flow around two-dimensional or axisymmetric bodies. The new code is an extension of the upwind (perfect gas) PNS code of Lawrence et al. (1986). The upwind algorithm is based on Roe's flux-difference splitting scheme which has been modified to account for real gas effects. The algorithm solves the gas dynamic and species continuity equations in a 'loosely' coupled manner. The new code has been validated by computing the laminar flow (at free stream Mach number 25) of chemically reacting air over a wedge and a cone. The results of these computations are compared with the results from a centrally-differenced, fully coupled, nonequilibrium PNS code. The agreement is excellent, except in the vicinity of the shock wave where the present code exhibits superior shock capturing capabilities.

Description: Hypersonic wake flows behind the Aeroassist Flight Experiment (AFE) geometry are analyzed using two Navier-Stokes flow solvers. Many of the AFE wake features observed in ballistic-range shadowgraphs are simulated using a simple, two-dimensional semicylinder geometry at moderate angles of attack. At free-stream conditions corresponding to a Hypersonic Free Flight Facility (HFFF) AFE experiment, the three-dimensional base flow for the AFE geometry is computed using an ideal-gas, Navier-Stokes solver. The computed results agree reasonably well with the shadowgraphs taken at the HFFF. An ideal-gas and a nonequilibrium Navier-Stokes solver have been coupled and applied to the complete flow around the AFE vehicle at the free-stream conditions corresponding to a nomial trajectory point. Limitations of the coupled ideal-gas and nonequilibrium solution are discussed. The nonequilibrium base flow solution is analyzed for the wake radiation and the radiation profiles along various lines of sight are compared. Finally, the wake unsteadiness is predicted using experimental correlations and the numerical solutions. An adaptive grid code, SAGE, has been used in all the simulations to enhance the solution accuracy. The grid adaptation is found to be necessary in obtaining base flow solutions with accurate flow features.

Description: Experiments are being conducted in the NASA Ames Hypervelocity Free Flight Aerodynamic Facility to quantify the effects on turbulent convective heat transfer of surface roughness representative of a new class of 3D woven thermal protection system mRough-wall turbulent heat transfer measurements were obtained on ballistic-range models in hypersonic flight in the NASA Ames Hypervelocity Free Flight Aerodynamic Facility. Each model had three different surface textures on segments of the conic frustum: smooth wall, sand roughness, and a pattern roughness, thus providing smooth-wall and sand-roughness reference data for each test. The pattern roughness was representative of a woven thermal protection system material developed by NASA's Heatshield for Extreme Entry Environment Technology project. The tests were conducted at launch speeds of 3.2 km/s in air at 0.15 atm. Roughness Reynolds numbers, k+, ranged for 12 to 70 for the sand roughness, and as high as 200 for the pattern roughness. Boundary-layer parameters required for calculating k+ were evaluated using computational fluid dynamics simulations. The effects of pattern roughness are generally characterized by an equivalent sand roughness determined with a correlation developed from experimental data obtained on specifically-designed roughness patterns that do not necessarily resemble real TPS materials. Two sand roughness correlations were examined: Dirling and van Rij, et al. Both gave good agreement with the measured heat-flux augmentation for the two larger pattern roughness heights tested, but not for the smallest height tested. It has yet to be determined whether this difference is due to limitations in the experimental approach, or due to limits in the correlations used. Future experiments are planned that will include roughness patterns more like those used in developing the equivalent sand roughness correlations.aterials being developed by NASA's Heatshield for Extreme Entry Environment Technology (HEEET) project. Data were simultaneously obtained on sand-grain roughened surfaces and smooth surfaces, which can be compared with previously obtained data. Results are presented in this extended abstract for one roughness pattern. The full paper will include results from three roughness patterns representing virgin HEEET, nominal turbulent ablated HEEET, and twice the roughness of nominal turbulent ablated HEEET. Results will be used to compare with commonly used equivalent sand grain roughness correlations.

Description: Grid generation issues relating to the simulation of the X33 aerothermal environment using the GASP code are explored. Required grid densities and normal grid stretching are discussed with regards to predicting the fluid dynamic and heating environments with the desired accuracy. The generation of volume grids is explored and includes discussions of structured grid generation packages such as GRIDGEN, GRIDPRO and HYPGEN. Volume grid manipulation techniques for obtaining desired outer boundary and grid clustering using the OUTBOUND code are examined. The generation of the surface grid with the required surface grid with the required surface grid topology is also discussed. Utilizing grids without singular axes is explored as a method of avoiding numerical difficulties at the singular line.

Description: Modern Computational Fluid Dynamics (CFD) techniques were used to compute aerodynamic forces and moments of the Space Shuttle Orbiter in specific portions of contingency abort trajectory space. The trajectory space covers a Mach number range of 3.5-15, an angle-of-attack range of 20-60 degrees, an altitude range of 100-190 kft, and several different settings of the control surfaces (elevons, body flap, and speed brake). While approximately 40 cases have been computed, only a sampling of the results is presented here. The computed results, in general, are in good agreement with the Orbiter Operational Aerodynamic Data Book (OADB) data (i.e., within the uncertainty bands) for almost all the cases. However, in a limited number of high angle-of-attack cases (at Mach 15), there are significant differences between the computed results, especially the vehicle pitching moment, and the OADB data. A preliminary analysis of the data from the CFD simulations at Mach 15 shows that these differences can be attributed to real-gas/Mach number effects.

Description: Results of flow simulations of turbulent shock wave boundary layer interaction experiments performed in the LENS-II tunnel at CUBRC.

Description: Plasma diagnostic measurement campaigns in the NASA Ames Interaction Heating Facility (IHF) have been conducted over the last several years with a view towards characterizing the flow in the arc jet facility by providing data necessary for modeling and simulation. Optical emission spectroscopy has been used in the plenum and in the free jet of the nozzle. Radiation incident over a probe surface has also been measured using radiometry. Plenum measurements have shown distinct radial profiles of temperature over a range of operating conditions. For cases where large amounts of cold air are added radially to the main arc-heated stream, the temperature profiles are higher by as much as 1500 K than the profiles assumed in flow simulations. Optical measurements perpendicular to the flow direction in the free jet showed significant contributions to the molecule emission through inverse pre-dissociation, thus allowing determination of atom number densities from molecular emission. This has been preliminarily demonstrated with the N2 1st Positive System. Despite the use of older rate coefficients, the resulting atom densities are reasonable and surprisingly close to flow predictions.

Description: Results of simulations of flow of an arc-heated stream around a 14-inch diameter 45 sphere-cone configuration are presented. Computations are first benchmarked against pressure and heat flux measurements made using copper slug calorimeters of different shapes and sizes. The influence of catalycity of copper on computed results is investigated. Good agreements between predictions and measurements are obtained by assuming the copper slug to be partially catalytic to atomic recombination. With total enthalpy estimates obtained from these preliminary computations, calculations are then performed for the test article, with the nozzle and test article considered as an integrated whole the same procedure adopted for calorimeter simulations. The resulting heat fluxes at select points on the test article (points at which fully instrumented plugs were placed) are used in material thermal response code calculations. Predicted time histories of temperature are compared against thermocouple data from the instrumented plugs, and recession determined. Good agreement is obtained for in-depth thermocouples.

Description: The Arc Jet Facilities at NASA Ames Research Center generate test streams with enthalpies ranging from 5 MJ/kg to 25 MJ/kg. The present work describes a rigorous method, based on equilibrium thermodynamics, for calculating the bulk enthalpy of the flow produced in two of these facilities. The motivation for this work is to determine a dimensionally-correct formula for calculating the bulk enthalpy that is at least as accurate as the conventional formulas that are currently used. Unlike previous methods, the new method accounts for the amount of argon that is present in the flow. Comparisons are made with bulk enthalpies computed from an energy balance method. An analysis of primary facility operating parameters and their associated uncertainties is presented in order to further validate the enthalpy calculations reported herein.