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.
Boundary conditions prove critical. Conference paper describes algorithm for calculation of shocks in hypersonic flows of gases at chemical equilibrium. Although algorithm represents intermediate stage in development of reliable, accurate computer code for two-dimensional flow, research leading up to it contributes to understanding of what is needed to complete task.
Analytical and experimental techniques for modeling the aerothermodynamics of hypersonic flight are assessed, together with the problems which will be encountered in developing reusable hypersonic vehicles. Emphasis is placed on a numerical coupling between a nonequilibrium chemistry model and hypersonic flow kinematics. Finite difference and finite element descriptions of flow fields in which molecules encounter a shock wave and undergo various motion (and thereby energy) transformations are discussed. The effects of artificial smearing of the shock wave are considered in terms of the resulting effects on the distribution of the energies and chemical composition of the transition region. Results are provided from schlieren photographs of shock-tube experiments, Navier-Stokes calculations of axisymmetric flow over a conical body, and calculations using a spectral multidomain approach for chemically reacting flows.
The feasibility of using corrugated panels as a thermal protection system for an advanced space transportation vehicle was investigated. The study consisted of two major tasks: development of improved correlations for wind tunnel heat transfer and pressure data to yield design techniques, and application of the design techniques to determine if corrugated panels have application future aerospace vehicles. A single-stage-to-orbit vehicle was used to assess advantages and aerothermodynamic penalties associated with use of such panels. In the correlation task, experimental turbulent heat transfer and pressure data obtained on corrugation roughened surfaces during wind tunnel testing were analyzed and compared with flat plate data. The correlations and data comparisons included the effects of a large range of geometric, inviscid flow, internal boundary layer, and bulk boundary layer parameters in supersonic and hypersonic flow.
The general problem of calculating the flow fields associated with hypersonic airbreathing aircraft is presented. Unique aspects of hypersonic aircraft aerodynamics are introduced and their demands on computational fluid dynamics are outlined. Example calculations associated with inlet/forebody integration and hypersonic nozzle design are presented to illustrate the nature of the problems considered.
The general problem of calculating the flow fields associated with hypersonic airbreathing aircrafts is presented. Unique aspects of hypersonic airplane aerodynamics are introduced and their demands on CFD are outlined. Example calculations associated with inlet/forebody integration and hypersonic nozzle design are presented to illustrate the nature of the problems considered.
Significant improvement has been obtained in the numerical performance of a Monte Carlo scheme for the analysis of nonequilibrium gas dynamics through an implementation of the algorithm which takes advantage of vector hardware, as presently demonstrated through application to three different problems. These are (1) a 1D standing-shock wave; (2) the flow of an expanding gas through an axisymmetric nozzle; and (3) the hypersonic flow of Ar gas over a 3D wedge. Problem (3) is illustrative of the greatly increased number of molecules which the simulation may involve, thanks to improved algorithm performance.
This paper presents a nonequilibrium flow solver, implementation of the algorithm on unstructured meshes, and application to hypersonic flow past blunt bodies. Air is modeled as a mixture of five chemical species, namely O2, N2, O, NO, and N, having two temperatures namely translational and vibrational. The solution algorithm is a cell centered, point implicit upwind scheme that employs Roe's flux difference splitting technique. Implementation of this algorithm on unstructured meshes is described. The computer code is applied to solve Mach 15 flow with and without a Type IV shock interference on a cylindrical body of 2.5mm radius representing a cowl lip. Adaptively generated meshes are employed, and the meshes are refined several times until the solution exhibits detailed flow features and surface pressure and heat flux distributions. Effects of a catalytic wall on surface heat flux distribution are studied. For the Mach 15 Type IV shock interference flow, present results showed a peak heat flux of 544 MW/m2 for a fully catalytic wall and 431 MW/m(exp 2) for a noncatalytic wall. Some of the results are compared with available computational data.
Composite photographs of the surface flow and shadowgraphs of the shock wave pattern are presented that depict the hypersonic flow field of a typical delta wing space shuttle orbiter. Results from a wind tunnel test in air are given in side, oblique, and projected plan views for angles of attack from 0 to 60 deg. The tests were conducted at a Mach number of 7.4 and for Reynolds numbers based on body length of 6,000,000 and 9,000,000. The interrelationship is shown for the intersecting bow and wing leading edge waves with the surface flow for angles of attack for which: (1) the leading edge wave is attached, and (2) the leading edge wave is detached.
ven with great advances in computational techniques and computing power during recent decades, the modeling of unsteady separated flows, such as those encountered in the wake of a re-entry vehicle, continues to be one of the most challenging problems in CFD. Of most interest to the aerothermodynamics community is accurately predicting transient heating loads on the base of a blunt body, which would result in reduced uncertainties and safety margins when designing a re-entry vehicle. However, the prediction of heat transfer can vary widely depending on the turbulence model employed. Therefore, selecting a turbulence model which realistically captures as much of the flow physics as possible will result in improved results. Reynolds Averaged Navier Stokes (RANS) models have become increasingly popular due to their good performance with attached flows, and the relatively quick turnaround time to obtain results. However, RANS methods cannot accurately simulate unsteady separated wake flows, and running direct numerical simulation (DNS) on such complex flows is currently too computationally expensive. Large Eddy Simulation (LES) techniques allow for the computation of the large eddies, which contain most of the Reynolds stress, while modeling the smaller (subgrid) eddies. This results in models which are more computationally expensive than RANS methods, but not as prohibitive as DNS. By complimenting an LES approach with a RANS model, a hybrid LES/RANS method resolves the larger turbulent scales away from surfaces with LES, and switches to a RANS model inside boundary layers. As pointed out by Bertin et al., this type of hybrid approach has shown a lot of promise for predicting turbulent flows, but work is needed to verify that these models work well in hypersonic flows. The very limited amounts of flight and experimental data available presents an additional challenge for researchers. Recently, a joint study by NASA and CUBRC has focused on collecting heat transfer data on the backshell of a scaled model of the Orion Multi-Purpose Crew Vehicle (MPCV). Heat augmentation effects due to the presence of cavities and RCS jet firings were also investigated. The high quality data produced by this effort presents a new set of data which can be used to assess the performance of CFD methods. In this work, a hybrid LES/RANS model developed at North Carolina State University (NCSU) is used to simulate several runs from these experiments, and evaluate the performance of high fidelity methods as compared to more typical RANS models. .
Theoretical results are presented for the structure of the hypersonic flow field of a blunt slab delta wing at moderately high angle of attack. Special attention is devoted to the interaction between the boundary layer and the inviscid entropy layer. The results are compared with experimental data. The three-dimensional inviscid flow is computed numerically by a marching finite difference method. Attention is concentrated on the windward side of the delta wing, where detailed comparisons are made with the data for shock shape and surface pressure distributions. Surface streamlines are generated, and used in the boundary layer analysis. The three-dimensional laminar boundary layer is computed numerically using a specially-developed technique based on small cross-flow in streamline coordinates. In the rear sections of the wing the boundary layer decreases drastically in the spanwise direction, so that it is still submerged in the entropy layer at the centerline, but surpasses it near the leading edge. Predicted heat transfer distributions are compared with experimental data.
This paper describes the current status of the operating characteristics of the electric arc-driven shock-tube facility at Ames Research Center, focusing on its potential usefulness in the current and anticipated future applications. The paper specifically addresses the questions as to: (1) how well the behavior of the arc driver is understood and controlled, (2) how well the facility is equipped to test low-density, very-high-velocity nonequilibrium flow regimes, and (3) how closely the facility is expected to produce an equilibrium hypersonic flow when operated in shock-tunnel modes. For these issues, it is shown that: (1) a plasma kinetics model of the exploding wire closely describes the arc behavior in the driver, (2) the facility can produce a spectroscopically-clean flow in a low density regime with a shock velocity of 13 km/sec in air when used with an aluminum driven tube, and (3) when operated as a shock-tunnel, the high enthalpy flow in the test section is expected to deviate only slightly from the perfect equilibrium flow conditions at enthalpies corresponding to flight speeds of 5 km/sec or less.
A real gas, upwind, parabolized Navier-Stokes (PNS) code has been developed to compute the three-dimensional hypersonic flow of equilibrium air around various body shapes. 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 (1981) flux-difference splitting scheme which has been modified to account for real gas effects using the nearly exact approach of Vinokur and Liu (1988). Simplified curve fits are employed to obtain the thermodynamic and transport properties of equilibrium air. The new code has been validated by computing the M-infinity = 25 laminar flow of air over cones at various angles of attack. The results of these computations are compared with the results from a conventional centrally-differenced, real gas PNS code and the previous axisymmetric, upwind, real gas code. The agreement is excellent in all cases.
This paper describes several applications of the global velocity measurement technique known as particle image velocimetry to three different facilities at Langley Research Center. These three applications represent global velocity measurements in distinctly different flows, namely a reacting hydrogen - air flame, a low Reynolds number boundary layer, and hypersonic flow over a wedge. In addition, these applications illustrate the power of the instrument to obtain data such as vorticity and Reynolds stress. Finally, the ability to obtain particle dynamics data in supersonic/hypersonic flows similar to that which can be obtained with point measurement techniques is demonstrated.
Numerical predictions of the Mars Science Laboratory reaction control system jets interacting with a Mach 10 hypersonic flow are compared to experimental nitric oxide planar laser-induced fluorescence data. The steady Reynolds Averaged Navier Stokes equations using the Baldwin-Barth one-equation turbulence model were solved using the OVERFLOW code. The experimental fluorescence data used for comparison consists of qualitative two-dimensional visualization images, qualitative reconstructed three-dimensional flow structures, and quantitative two-dimensional distributions of streamwise velocity. Through modeling of the fluorescence signal equation, computational flow images were produced and directly compared to the qualitative fluorescence data.
Computational results are presented for hypersonic viscous flow past spinning sharp and blunt cones of angle of attack, obtained with a parabolic Navier-Stokes marching code. The code takes into account the asymmetries in the flowfield resulting from spinning motion and computes the asymmetric shock shape, cross-flow and streamwise shear, heat transfer, cross-flow separation, and vortex structure. The Magnus force and moments are also computed. Comparisons are made with other theoretical analyses based on boundary-layer and boundary-region equations, and an anomaly is discovered in the displacement thickness contribution to the Magnus force when compared with boundary-layer results.
Ablative materials are used as the Thermal Protection System (TPS) for space vehicles that enter a planetary atmosphere. During the entry phase, the vehicle experiences high rates of aerodynamic heating. These ablative materials, acting as a heat shield, protect the vehicle from the harsh environment and ensure a safe entry. Ablative materials reduce the incoming heat flux through various phenomena such as near-surface oxidation, pyrolysis chemical reaction, vaporization, and other erosive processes, collectively known as ablation. Spallation is one such ablative process, in which solid particles are ejected from the material into the flow field. The many conditions leading to spallation, however, are not well understood and cannot be described accurately by simple phenomenological models. However, the spallation phenomenon might impact the aerodynamic heating rates of re-entry vehicles using ablative heat shields. For example, the solid particles ejected from the surface can modify the near surface chemical processes and flow field. To investigate spallation effects, a code was developed to compute the dynamics of spalled particles. The code uses a finite-rate chemistry model to study the chemical interactions of the particles with the flow field. The spallation code is dynamically coupled to a CFD solver that models the time accurate hypersonic flow field around an ablative sample. An experimental investigation was also performed at the NASA Langley HYMETS facility. HYMETS is a 400 kW arcjet that can simulate hypersonic flight and Earth or Martian entry conditions. Using high speed cameras, an IR camera, NIR and UV-VIS spectrometers, and other diagnostic instruments, test data on various graphite, LI-900, FiberForm, and PICA test samples were obtained, using target heat flux conditions of 100, 200, and 400 W/cm 2 . From the test campaign, measurements of particle velocity, acceleration and trajectory are obtained.
Application of a general method for calculation of inviscid hypersonic flow fields is discussed. General considerations are analyzed along with the sonic corner and the stagnation region. It is concluded that the complications caused by the requirement for sonic flow at the rear corner and particularly by the uncertain position of the stagnation streamline lead to sufficient difficulties with convergence of iterations that a practical procedure is not likely to be found.
Numerical predictions of the Mars Science Laboratory (MSL) reaction control system (RCS) jets interacting with a Mach 10 hypersonic flow are compared to experimental nitric oxide (NO) planar laser-induced fluorescence (PLIF) data. The steady Reynolds Averaged Navier Stokes (RANS) equations using the Baldwin-Barth one-equation turbulence model were solved using the OVERFLOW code. The experimental PLIF data used for comparison consists of qualitative two-dimensional visualization images, qualitative reconstructed three-dimensional flow structures, and quantitative two-dimensional distributions of streamwise velocity. Through modeling of the PLIF signal equation, computational flow images (CFI) were produced and directly compared to the qualitative PLIF data. Post processing of the experimental and simulation data enabled the jet trajectory to be extracted for a more quantitative comparison. The two-dimensional velocity fields were reconstructed through interpolation of a series of single-component velocity profiles. Each distribution of single-component velocity was obtained using molecular tagging velocimetry (MTV). After validating the numerical model, the numerical solution was further examined to gain insight into hypersonic jet-in-crossflow interaction. Future NO-PLIF experiments are proposed based on this analysis.