Engineering topics
Publications and source records attributed to Davy, W. C..
Measurements of upper and lower shear turning angles and the wake-neck location that were obtained in the NASA Ames' ballistic range are presented. The upper shear turning angles measured in the ballistic range at velocities of 1.8 and 4.9 km/s are compared to the wind tunnel results.
This paper describes the aerothermal environment in the afterbody region of a blunt entry body. Recent ground-based experiments and computational predictions of the afterbody flow structure and radiation are presented. The similarity between the flowfield structures observed in the ground-based experiments and that obtained by calculation is encouraging. Approximate calculations of the radiative heating rate to the base are presented. Many of the phenomena associated with the expanding flow at the corner and the formation of the wake neck, however, are not well understood and require further study. A flight experiment is described that would use spectral and total measurements of the wake radiation as a nonintrusive diagnostic method to provide insight into the thermodynamic state of the wake gas.
This paper describes the rationale for conducting the proposed radiation measurement in the afterbody region of the Aeroassist Flight Experiment (AFE) vehicle, the results of the calculations of expected radiation intensities performed to date, and the instrumentation for the measurement. The experiment, named the Afterbody Radiometry Experiment (ARE), is one of the experiments that will be carried on board the AFE vehicle. The paper collates the existing experimental data to show the possibility that there may be substantial radiative heating of the afterbody region of an Aeroassisted Space Transfer Vehicle, which necessitates an experiment of the kind proposed. Calculations of the radiative heating rate to the base of the AFE are, qualitatively, in fair agreement with previous experimental work, which indicated that the radiative heating of the base is about 5 percent of that at the stagnation point.
The High Frequency Radiometer (HFR) is the only instrument on the Aeroassist Flight Experiment (AFE) with sufficient temporal resolution to discern the frequency of unsteady wake oscillations. Determining both the frequency and amplitude of wake unsteadiness during AFE atmospheric entry is essential for reliably predicting the geometry and motion of the wake of future Aeroassisted Space Transfer Vehicles (ASTV). These parameters directly affect the location and size of the pay- load and the weight of the required afterbody heat protection. The purpose of the AFE is to validate the technologies required for the design of ASTVs, which will be used to exploit Earth-lunar space. This validation will be conducted at a combination of vehicle size, altitude, and velocity not obtainable in ground-based facilities. The AFE will provide the experimental flight data needed to improve our understanding of hypersonic-wake physics and to validate computational predictions of the aerodynamic and heating loads, including afterbody radiative heating loads, on an ASTV. Reliable prediction of ASTV wake flows will ensure that payloads are located within the shear-layer envelope and will determine the amount of thermal protection the payloads require. Specifically, understanding the temporal nature of the wake unsteadiness is important for two reasons. Most importantly, analysis of ground-based experiments suggests that wake unsteadiness results in a variation of as much as +/- 5 deg in the shear-flow turning angle. This angle must be reliably predicted to avoid shear-layer impingement on the vehicle afterbody, which would result in heating rates of about 10 W/sq cm, of the same order as on the forebody stagnation point. Secondly, the energy associated with wake unsteadiness will reduce the static enthalpy of the wake fluid and cause an error of as much as 30% in the amount of predicted wake radiative heating. Therefore, the HFR flight data, which will quantify the frequency and amplitude of the wake unsteadiness, are required for the verification of computational models of ASTV flowfields. Measurement of wake unsteadiness has been endorsed by the Peer Science Steering Group and the AFE Computational Fluid Dynamics (CFD) Working Group. The Peer Science Steering Group has stressed that a reliable measurement of the wake unsteadiness is fundamental because incorporating unsteadiness represents a substantial challenge to the CFD community, and a reliable measurement will raise confidence in the computer simulation. The AFE CFD Working Group has stated that 44 any data which would prove the existence and effects of unsteady flow would be extremely valuable.
Radiation heat transfer (RHT) from the wake of a hypersonic vehicle to its afterbody is evaluated from Gnoffo and Greene's (1987) calculated wake flowfield and the radiative properties of ionized high-temperature air with the calculated nonequilibrium composition. The 4.2-m aeroassisted flight experiment at an altitude of 75 km and velocity of 8900 m/s causes a 0.1-m-thick layer initially at T = 10,000 K and P = 1 kN/sq m to separate from the shoulder of the forebody heat shield and spread aft to form a wake at approximately T = 5000 K and P = 20 N/sq m. Gas in the separated flow region at approximately T = 3000 K and P = 10 N/sq m, recirculates about the afterbody. It is shown that the radiating layer, recirculating gases, and wake are optically thin for purposes of making engineering RHT calculations. Directional, spectral, and spatial variations of the radiation incident upon the afterbody are presented.
A forthcoming NASA flight experiment is described that provides an opportunity to obtain a large base of radiometric data for high-altitude, high-velocity thermochemically nonequilibrated-flow conditions. As a preliminary to the design of a radiometer for this experiment, an approximate method for predicting both equilibrium and nonequilibrium radiative surface fluxes is described. Spectral results for one trajectory state, a velocity of 10 km/sec at an altitude of 85 km, are presented. These results are then used to develop some of the instrument parameters that will be needed for designing of the three genre of radiometers that are proposed for this experiment.
The Conservative Supra Characteristic Method (CSCM), an implicit upwind Navier-Stokes algorithm, is extended to the numerical simulation of flows in chemical equilibrium. The resulting computer code known as Chemistry and Gasdynamics Implicit - Version 2 (CAG12) is described. First-order accurate results are presented for inviscid and viscous Mach 20 flows of air past a hemisphere-cylinder. The solution procedure captures the bow shock in a chemically reacting gas, a technique that is needed for simulating high altitude, rarefied flows. In an initial effort to validate the code, the inviscid results are compared with published gasdynamic and chemistry solutions and satisfactorily agreement is obtained.
Explore the source record for details and available documents.
Laminar, real gas hypersonic flowfields over a three dimensional configuration are computed using an unsteady, factored implicit scheme. Local chemical and thermodynamic properties are evaluated by an equilibrium composition method. Transport properties are obtained from individual species properties and application of a mixture rule. Numerical solutions are presented for an ideal gas and equilibrium air for free-stream Mach numbers of 13 and 15 and at various angles of attack. The effect of real gas is to decrease the shock-layer thickness resulting from decreased shock-layer temperatures and corresponding increased density. The combined effects of viscosity and real gas are to increase the subsonic layer near the wall.
Real gas, inviscid supersonic flow fields over a three-dimensional configuration are determined using a factored implicit algorithm. Air in chemical equilibrium is considered and its local thermodynamic properties are computed by an equilibrium composition method. Numerical solutions are presented for both real and ideal gases at three different Mach numbers and at two different altitudes. Selected results are illustrated by contour plots and are also tabulated for future reference. Results obtained compare well with existing tabulated numerical solutions and hence validate the solution technique.
Implicit methods developed by Beam and Warming (1978 and Briley and McDonald (1977) make it possible to overcome the hyperbolic stiffness of the conservative compressible Navier-Stokes equations in the fine wall region computational mesh for high Reynold's number flow. Certain difficulties related to the use of these methods could be overcome by employing an approach reported by Roe (1981). In the present investigation Roe's conceptual framework has been adopted for constructing globally conservative finite difference methods. A globally conservative upwind finite difference method (CSCM) consisting of both implicit interior point and boundary point equations is constructed from a new characteristics based flux difference splitting. It is found that the employed upwind eigenvector split scheme which combines fully coupled implicit interior point and boundary point approximations has the desired properties of robust stability and accuracy.
Material response solutions for the forebody heat shield on the candidate 310-kg Galileo Probe are presented. A charring material ablation analysis predicts thermochemical surface recession, insulation thickness, and total required heat shield mass. Benchmark shock layer solutions provide the imposed entry heating environments on the ablating surface. Heat shield sizing results are given for a nominal entry into modeled nominal and cool-heavy Jovian atmospheres, and for two heat-shield property models. The nominally designed heat shield requires a mass of at least 126 kg and would require an additional 13 kg to survive entry into the less probable cool-heavy atmosphere. The material-property model with a 30% surface reflectance reduces these mass requirements by as much as 16%.
The factored-implicit, gas-dynamic algorithm has been adapted to the numerical simulation of equilibrium reactive flows. Changes required in the perfect gas version of the algorithm are developed, and the method of coupling gas-dynamic and chemistry variables is discussed. A flow-field solution that approximates a Jovian entry case was obtained by this method and compared with the same solution obtained by HYVIS, a computer program much used for the study of planetary entry. Comparison of surface pressure distribution and stagnation line shock-layer profiles indicates that the two solutions agree well.
Solutions are presented for the aerothermal heating environments and the material thermal response for the forebody heatshield on the candidate 242 kg Galileo probe entering the modeled nominal and cold-dense Jovian atmospheres. In the flowfield analysis, a finite difference procedure was employed to obtain benchmark predictions of pressure, radiation and convective heating rates (both laminar and turbulent) and the corresponding wall blowing obtained under the steady state approximation. The fluxes over the probe flank were found to be in a range where spallation is an important mass loss mechanism. The predicted heating rates were also used as boundary conditions for a charring materials ablation which was used to predict thermochemical based surface recession, mass loss and bondline temperatures. The contingency factor of 30% currently employed by NASA was found to be insufficient for entry into the cold-dense atmosphere.
A new code for the simulation of full (forebody and base region) flowfields about bluff bodies in the hypersonic regime of severe planetary entry is described. The present 'maximally conservative, maximally differenced' formulation of the unsteady compressible Navier-Stokes equations for 2-D axisymmetric 3-D flow is contrasted for stability with previous formulations of Viviand, Kutler, et al, and Thomas and Lombard. Discrete metric relations peculiar to the axisymmetric finite volume formulation are presented along with a general discussion of their relations to and consequences of failure to close computational cells. A computational mesh of curvilinear coordinate topology singular in the flow regime is presented that permits aligned capturing of the major physical features of the complex flowfield.
Utilizing a series of existing computer codes, ablation experiments in the Giant Planet Facility are numerically simulated. Of primary importance is the simulation of the low Mach number shock layer that envelops the test model. The RASLE shock-layer code, used in the Jupiter entry probe heat-shield design, is adapted to the experimental conditions. RASLE predictions for radiative and convective heat fluxes are in good agreement with calorimeter measurements. In simulating carbonaceous ablation experiments, the RASLE code is coupled directly with the CMA material response code. For the graphite models, predicted and measured recessions agree very well. Predicted recession for the carbon phenolic models is 50% higher than that measured. This is the first time codes used for the Jupiter probe design have been compared with experiments.
The ablation of carbonaceous materials in a hydrogen-helium stream has been simulated using a charring materials ablator computer code. These results are compared with the first ablation data to be obtained from the Ames-NASA Giant Planet Pilot Facility. Test stream diagnostics and ablation effects on convective and radiative heat transfer are discussed since these parameters constitute important input data to the numerical simulation. Graphite ablation was predicted to within 10 to 20%, and carbon-phenolic somewhat less accurately.
The nonequilibrium axisymmetric stagnation point boundary layer over an ablating graphite surface is considered. The external stream is a high temperature mixture of hydrogen and helium. Variable thermodynamic and transport properties are assumed. Lennard-Jones potential model is used to calculate the transport coefficients of each species. Although the mixture rules for viscosity of the gas mixture are used, the weighting functions are more sophisticated than those commonly employed. For the conductivity of the mixture, generalized Wassiljewa coefficients are used. Seven species with 28 dissociation/recombination reactions are considered. Hansen's model for the dissociation rate constants is employed. The recombination rate constants are obtained by invoking detailed balance principles assisted by the JANAF thermodynamic data and the Hansen-Pearson thermodynamic data for C3.