Sphere drag in near-free-molecule hypersonic flow
Spheres drag coefficient at hypersonic Mach numbers for near free molecular flow
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.
Spheres drag coefficient at hypersonic Mach numbers for near free molecular flow
The general three-dimensional direct simulation Monte Carlo method is used to study the hypersonic flow around two blunt wedges that intersect at a 90 deg angle. Results are obtained for the transitional flow regimes found at 85 and 100 km altitude with a reentry velocity of 7.5 km/s. The disturbance field in front of the double-wedge body is found to be larger than that produced by a single wedge. Surface pressures and flow densities are higher near the wedge intersection, whereas surface heating and shear streses are greater at locations removed from the corner. Results also show that three-dimensional flow structure occurs only near the wedge corner, and that the flow monotonically approaches the limiting two-dimensional wedge flow case in the spanwise direction.
An approximate axisymmetric method has been developed which can reliably calculate nonequilibrium fully viscous hypersonic flows over blunt-nosed bodies. By substituting Maslen's second-order pressure expression for the normal momentum equation, a simplified form of the viscous shock layer (VSL) equations is obtained. This approach can solve both the subsonic and supersonic regions of the shock layer without a starting solution for the shock shape. This procedure is significantly faster than the parabolized Navier-Stokes and VSL solvers and would be useful in a preliminary design environment. Solutions have been generated for air flows over several analytic body shapes. Surface heat transfer and pressure predictions are comparable to VSL results. Computed heating rates are in good agreement with experimental data. The present technique generates its own shock shape as part of its solution, and therefore could be used to provide more accurate initial shock shapes for higher-order procedures which require starting solutions.
The interaction between a swept shock wave and a laminar boundary layer was investigated experimentally in high-enthalpy hypersonic flow. The effect of high-temperature, real gas physics on the interaction was examined by conducting tests in air and helium. Heat transfer measurements were made on the surface of a flat plate and a shock-generating fin using thin-film resistance sensors for fin incidence angles of 0, 5, and 10 deg at Mach numbers of 6.9 in air and 7.2 in helium. The experiments were conducted in the NASA HYPULSE expansion tube, an impulse-type facility capable of generating high-enthalpy, high-velocity flow with freestream levels of dissociated species that are particularly low. The measurements indicate that the swept shock wave creates high local heat transfer levels in the interaction region, with the highest heating found in the strongest interaction. The maximum measured heating rates in the interaction are order of magnitude greater than laminar flat plate boundary layer heating levels at the same location.
We describe the extension of a 2-D simplified face-averaged nodal-gradient (F-ANG) method to 3-D and demonstrate that the 3-D simplified F-ANG method is accomplished by augmenting the nodecentered gradient least squares stencil. This augmented stencil F-ANG method is shown to result in advection and diffusion schemes that are stable for hexahedral, prismatic, pyramidal and tetrahedral cells without having to resort to cell-averaged nodal gradients. In addition, we describe the modifications to the augmented stencil required to support the use of wall function boundary conditions. Finally we describe a consistent, face-stencil based multi-dimensional limiter procedure (MLP), and show it to be fully consistent and compatible with the linearity-preserving unstructured- MUSCL (LP-U-MUSCL) scheme for all values of kappa. These methods and schema are implemented in the cell-centered finite-volume code VULCAN-CFD, which is then used to investigate whether the robustness improvements demonstrated in 2-D carry over to 3-D by computing hypersonic flows using mixed-element grids as well as highly adapted tetrahedral grids.
The general status of a grant to investigate the applications of quantum theory in atomic and molecular scattering problems in hypersonic flow is summarized. Abstracts of five articles and eleven full-length articles published or submitted for publication are included as attachments. The following topics are addressed in these articles: fragmentation of heavy ions (HZE particles); parameterization of absorption cross sections; light ion transport; emission of light fragments as an indicator of equilibrated populations; quantum mechanical, optical model methods for calculating cross sections for particle fragmentation by hydrogen; evaluation of NUCFRG2, the semi-empirical nuclear fragmentation database; investigation of the single- and double-ionization of He by proton and anti-proton collisions; Bose-Einstein condensation of nuclei; and a liquid drop model in HZE particle fragmentation by hydrogen.
The development of new hypersonic flight vehicles is limited by the physical understanding that may be obtained from ground test facilities. This has motivated the present development of a temporally and spatially resolved velocimetry measurement for Sandia National Laboratories (SNL) Hypersonic Wind Tunnel (HWT) using Femtosecond Laser Electronic Excitation Tagging (FLEET). First, a multi-line FLEET technique has been created for the first time and tested in a supersonic jet, allowing simultaneous measurements of velocities along multiple profiles in a flow. Secondly, two different approaches have been demonstrated for generating dotted FLEET lines. One employs a slit mask pattern focused into points to yield a dotted line, allowing for two- or three-component velocity measurements free of contamination between components. The other dotted-line approach is based upon an optical wedge array and yields a grid of points rather than a dotted line. Two successful FLEET measurement campaigns have been conducted in SNL’s HWT. The first effort established optimal diagnostic configurations in the hypersonic environment based on earlier benchtop reproductions, including validation of the use of a 267 nm beam to boost the measurement signal-to-noise ratio (SNR) with minimal risk of perturbing the flow and greater simplicity than a comparable resonant technique at 202 nm. The same FLEET system subsequently was reconstituted to demonstrate the ability to make velocimetry measurements of hypersonic turbulence in a realistic flow field. Mean velocity profiles and turbulence intensity profiles of the shear layer in the wake of a hypersonic cone model were measured at several different downstream stations, proving the viability of FLEET as a hypersonic diagnostic.
The finite element method was applied to compressible flows past complete blunt bodies in a first attempt to use the method to model hypersonic flow past re-entry vehicles. A Runge-Kutta time marching scheme was utilized and was found to adequately handle the bow shock and downstream wake. Adaptive refinement was successively applied to the flow field grid, and improvement in solution accuracy was observed. The results reported in this paper include the flow past a circular cylinder at Mach 3, and the flow in the plane of symmetry for the flow past a sphere at Mach 10.
The approximate axisymmetric method presented for accurately calculating the surface and flowfield properties of fully viscous hypersonic flow over blunt-nosed bodies incorporates the turbulence model of Cebeci-Smith (1970) and the equilibrium air tables of Hansen (1959). The method is faster than the parabolized Navier-Stokes or viscous shock layer solvers that it could replace for preliminary design determinations. Surface heat transfer and pressure predictions for the present method are comparable with the more accurate viscous shock layer method as well as flight test and wind tunnel data. A starting solution is not required.
This paper reviews techniques for transitional- and turbulent-flow measurements and describes current research in support of turbulence modeling. Special attention is given to the potential of applying hot wire and laser velocimeter to measuring turbulent fluctuations in hypersonic flow fields. The results of recent experiments conducted in two hypersonic wind tunnels are presented and compared with previous hot-wire turbulence measurements.
An analytical approximation is obtained to the solution of the equation describing the combined effects of bluntness and boundary-layer displacement on the hypersonic flow over a wedge for the case where the constant parameter in the governing equation, which is proportional to the angle of attack, is greater than zero. It is shown that the approximation has good physical accuracy over the entire range of conditions to which the theory is applicable for positive angle of attack.
An experimental study of surface pressure distributions on a family of blunt and sharp large angle cones was made in hypersonic flows of helium, air, and tetrafluoromethane. The effective isentropic exponents of these flows were 1.67, 1.40, and 1.12. Thus, the effect of large shock density ratios such as might be encountered during planetary entry because of real-gas effects could be studied by comparing results in tetrafluoromethane with those in air and helium. It was found that shock density ratio had a large effect on both shock shape and pressure distribution. The differences in pressure distribution indicate that for atmospheric flight at high speed where real-gas effects produce large shock density ratios, large-angle cone vehicles can be expected to experience different trim angles of attack, drag coefficient, and lift-drag ratios than those for ground tests in air wind tunnels.
As the capability of the space transportation vehicles (STV's) expand to meet the requirements for future space exploration and utilization, the effects of rarefied hypersonic flows will play a more significant role in defining the aerodynamic and aerothermodynamic performance of STV's. This is particularly true of the low lift/drag aeroassisted STV's where aerobraking occurs at relatively high altitudes and high velocity. Because of the limitations of the continuum description as expressed by the Navier-Stokes equations and the difficulties of solving the Boltzmann equation, the particle of molecular approach has been developed over the last three decades for modeling rarefied gas effects. The direct simulation Monte Carlo (DSMC) method of Bird is the most used method today for simulating rarefied flows. The DSMC method provides a direct physical simulation as opposed to a numerical solution of a set of model equations. This is accomplished by developing phenomenological models of the relevant physical events. The DSMC method accounts for translational, thermal, chemical, and radiative nonequilibrium effects. The general features of the DSMC method, the numerical requirements for obtaining meaningful results, the modeling used to simulate high temperature gas effects, and applications of the method to calculate the flow about an aeroassist flight experiment vehicle (AFE) are reviewed. The AFE simulates a geosynchronous return while entering the Earth's upper atmosphere at approximately 10 km/s. Results obtained using a general 3-D code are presented for the more rarefied portion of the atmospheric encounter (altitudes of 200 to 100 km) emphasizing surface, flowfield, and aerodynamic characteristics of the AFE. Finally, results obtained using axisymmetric and 1-D versions of the code are presented for lower altitude conditions.
Hypersonic viscous flow over spherically blunted cones of large half angle is computed at small angles of attack in the plane of symmetry of the flow field. Time-dependent viscous shock-layer equations in body-oriented coordinates are used to describe the flow field. The shock wave is treated as a discontinuity, across which the Rankine-Hugoniot relations are used to compute the flow conditions behind the shock. A time-marching second-order finite-difference method is used to solve the equations for a perfect gas. The local CFL (Courant-Friedrich-Lewy) time increment is used to advance the solution in time at each grid point. A fourth-order damping is used to damp the oscillations in the flow quantities. The numerical results of the present analysis for quantities such as shock standoff distance, surface-pressure distribution, and heating rates compare well with the existing theoretical and experimental results.
Hypersonic transitional flow predictions have been made using the parabolized Navier-Stokes equations with an algebraic transition/turbulence model by appropriately modulating the turbulent viscosity with the available intermittency functions for incompressible and compressible flows. A comparison between the predictions with and without a low Reynolds number correction has also been made. The predictions are compared with the available experimental data and with the theory over a range of Mach number. A simple Re(theta)/M(delta) criterion is shown to satisfactorily predict the meridional variation of the onset location of transition on a cone at a small angle of attack, whereas none of the correlations discussed can do that. Various available correlations are discussed vis-a-vis the predictions as to the locations of the onset and the end of transition.
Experimental studies have been conducted to examine slot film cooling effectiveness and the interaction between the cooling film and an incident planar shock wave in turbulent hypersonic flow. The experimental studies were conducted in the 48-inch shock tunnel at Calspan at a freestream Mach number of close to 6.4 and at a Reynolds number of 35 x 10(exp 6) based on the length of the model at the injection point. The Mach 2.3 planar wall jet was generated from 40 transverse nozzles (with heights of both 0.080 inch and 0.120 inch), producing a film that extended the full width of the model. The nozzles were operated at pressures and velocities close to matching the freestream, as well as at conditions where the nozzle flows were over- and under-expanded. A two-dimensional shock generator was used to generate oblique shocks that deflected the flow through total turnings of 11, 16, and 21 degrees; the flows impinged downstream of the nozzle exits. Detailed measurements of heat transfer and pressure were made both ahead and downstream of the injection station, with the greatest concentration of measurements in the regions of shock-wave/boundary layer interaction. The major objectives of these experimental studies were to explore the effectiveness of film cooling in the presence of regions of shock-wave/boundary layer interaction and, more specifically, to determine how boundary layer separation and the large recompression heating rates were modified by film cooling. Detailed distributions of heat transfer and pressure were obtained in the incident shock/wall-jet interaction region for a series of shock strengths and impingement positions for each of the two nozzle heights. Measurements were also made to examine the effects of nozzle lip thickness on cooling effectiveness. The major conclusion from these studies was that the effect of the cooling film could be readily dispersed by relatively weak incident shocks, so the peak heating in the recompression region was not significantly reduced by even the largest levels of film cooling. For the case studies in the absence of film cooling, the interaction regions were unseparated. However, adding film cooling resulted in regions of boundary layer separation induced in the film cooling layer -- the size of which regions first increased and then decreased with increased film cooling. Surprisingly, the size of the separated regions and the magnitude of the recompression heating were not strongly influenced by the thickness of the cooling film, nor by the point of shock impingement relative to the exit plane of the nozzles. The lip thickness was found to have little effect on cooling effectiveness. Measurements with and in the absence of shock interaction were compared with the results of earlier experimental studies and correlated in terms of the major parameters controlling these flows.
A user's guide is provided for a computer code which calculates the laminar and turbulent hypersonic flows about blunt axisymmetric bodies, such as spherically blunted cones, hyperboloids, etc., at zero and small angles of attack. The code is written in STAR FORTRAN language for the CDC-STAR-100 computer. Time-dependent, viscous-shock-layer-type equations are used to describe the flow field. These equations are solved by an explicit, two-step, time asymptotic, finite-difference method. For the turbulent flow, a two-layer, eddy-viscosity model is used. The code provides complete flow-field properties including shock location, surface pressure distribution, surface heating rates, and skin-friction coefficients. This report contains descriptions of the input and output, the listing of the program, and a sample flow-field solution.
Experimental studies have been conducted to examine slot film cooling effectiveness and the interaction between the cooling film and an incident planar shock wave in turbulent hypersonic flow. The experimental studies were conducted in the 48-inch shock tunnel at Calspan at a freestream Mach number of close to 6.4 and at a Reynolds number of 35 x 10(exp 6) based on the length of the model at the injection point. The Mach 2.3 planar wall jet was generated from 40 transverse nozzles (with heights of both 0.080 inch and 0.120 inch), producing a film that extended the full width of the model. The nozzles were operated at pressures and velocities close to matching the freestream, as well as at conditions where the nozzle flows were over- and under-expanded. A two-dimensional shock generator was used to generate oblique shocks that deflected the flow through total turnings of 11, 16, and 21 degrees; the flows impinged downstream of the nozzle exits. Detailed measurements of heat transfer and pressure were made both ahead and downstream of the injection station, with the greatest concentration of measurements in the regions of shock-wave/boundary layer interaction. The major objectives of these experimental studies were to explore the effectiveness of film cooling in the presence of regions of shock-wave/boundary layer interaction and, more specifically, to determine how boundary layer separation and the large recompression heating rates were modified by film cooling. Detailed distributions of heat transfer and pressure were obtained in the incident-shock/wall-jet interaction region for a series of shock strengths and impingement positions for each of the two nozzle heights. Measurements were also made to examine the effects of nozzle lip thickness on cooling effectiveness. The major conclusion from these studies was that the effect of the cooling film could be readily dispersed by relatively weak incident shocks, so the peak heating in the recompression region was not significantly reduced by even the largest levels of film cooling. For the case studies in the absence of film cooling, the interaction regions were unseparated. However, adding film cooling resulted in regions of boundary layer separation induced in the film cooling layer, the size of which regions first increased and then decreased with increased film cooling. Surprisingly, the size of the separated regions and the magnitude of the recompression heating were not strongly influenced by the thickness of the cooling film, nor by the point of shock impingement relative to the exit plane of the nozzles. The lip thickness was found to have little effect on cooling effectiveness. Measurements with and in the absence of shock interaction were compared with the results of earlier experimental studies and correlated in terms of the major parameters controlling these flows.