Engineering PapersSearch

Engineering topics

Graves, R. A., Jr.

Publications and source records attributed to Graves, R. A., Jr..

At least 19 records

Hypersonic aerothermodynamics status report

This survey paper gives an overview of NASA's Aerothermodynamics program. The purpose is to present the elements of, example results from, and projected technology needs for hypersonic fluid and thermal physics. The NASA program is focused on improving the fundamental understanding of aerodynamic and aerothermodynamic flow phenomena over hypersonic vehicles operating in the continuum, transitional, and free molecule flow regime. Vehicle design capabilities, computational fluid dynamics, computational chemistry, turbulence modeling, aerothermal loading, Shuttle Orbiter flight data analysis, and aerothermal facilities are discussed.

Graves, R. A., Jr.

Computational fluid dynamics - The coming revolution

The development of aerodynamic theory is traced from the days of Aristotle to the present, with the next stage in computational fluid dynamics dependent on superspeed computers for flow calculations. Additional attention is given to the history of numerical methods inherent in writing computer codes applicable to viscous and inviscid analyses for complex configurations. The advent of the superconducting Josephson junction is noted to place configurational demands on computer design to avoid limitations imposed by the speed of light, and a Japanese projection of a computer capable of several hundred billion operations/sec is mentioned. The NASA Numerical Aerodynamic Simulator is described, showing capabilities of a billion operations/sec with a memory of 240 million words using existing technology. Near-term advances in fluid dynamics are discussed.

Graves, R. A., Jr.

Complete viscous flowfield solutions about a blunt parabolic body in a supersonic stream

A hybrid computational technique which splits the flowfield into inviscid and viscous regions is used to investigate the complete flowfield about axisymmetric parabolic blunt bodies in a supersonic stream. The solutions are carried out on the CDC CYBER-203 computer which, with its extensive memory, allows for the use of a large number of finite-difference mesh points, allowing resolution of important flowfield features. A range of freestream Mach number of 2-5 and a range of Re number based on nose radius of 500-125,000 was run for a sonic corner body. Contour plots of density, pressure, and Mach number, velocity vector plots, and surface distributions of pressure, heat transfer, and shear stress are presented. Also, correlations of the downstream extent of the base recirculation region with Re number based on nose radius are given.

Weilmuenster, K. J.

Viscous compressible flow about blunt bodies using a numerically generated orthogonal coordinate system

A numerical solution to the Navier-Stokes equations was obtained for blunt axisymmetric entry bodies of arbitrary shape in supersonic flow. These equations are solved on a finite difference mesh obtained from a simple numerical technique which generates orthogonal coordinates between arbitrary boundaries. The governing equations are solved in time dependent form using Stetter's improved stability three step predictor corrector method. For the present application, the metric coefficients were obtained numerically using fourth order accurate, finite difference relations and proved to be totally reliable for the highly stretched mesh used to resolve the thin viscous boundary layer. Solutions are obtained for a range of blunt body nose shapes including concavities.

Graves, R. A., Jr.

Flowfield calculations past an entry probe with ablated nose shape

The effects of ablated nose shapes on the flowfield solutions are studied, using a time-dependent finite-difference method developed by Kumar, et al. (1979). Solutions are obtained for the laminar flow of a radiating mixture of H-He in chemical equilibrium past a blunt axisymmetric body at zero angle of attack. The freestream conditions correspond to a point on a typical Jovian entry trajectory, and the initial probe shape is a 45-deg half-angle spherically blunted cone. It is found that as nose bluntness increases, the following occur: in the nose region, shock standoff distances and radiative heating rates increase substantially; surface pressure level increases, but convective heating rates decrease.

Kumar, A.

User's guide for vectorized code EQUIL for calculating equilibrium chemistry on Control Data STAR-100 computer

A vectorized code, EQUIL, was developed for calculating the equilibrium chemistry of a reacting gas mixture on the Control Data STAR-100 computer. The code provides species mole fractions, mass fractions, and thermodynamic and transport properties of the mixture for given temperature, pressure, and elemental mass fractions. The code is set up for the electrons H, He, C, O, N system of elements. In all, 24 chemical species are included.

Kumar, A.

Application of a numerically generated orthogonal coordinate system to the solution of inviscid axisymmetric supersonic flow over blunt bodies

A numerically generated orthogonal coordinate system (with the body surface and shock wave as opposite boundaries) was applied with a time asymptotic method to obtain steady flow solutions for axisymmetric inviscid flow over several blunt bodies including spheres, paraboloids, ellipsoids, hyperboloids, hemisphere cylinders, spherically blunted cones, and a body with a concavity in the stagnation region. Comparisons with experimental data and with the results of other computational methods are discussed. The numerically generated orthogonal coordinate system is described and applications of the method to complex body shapes, particularly those with concave regions, are discussed.

Hamilton, H. H., II

Laminar and turbulent flow solutions with radiation and ablation injection for Jovian entry

Laminar and turbulent flowfield solutions with coupled carbon-phenolic mass injection are presented for the forebody of a probe entering a nominal Jupiter atmosphere. Solutions are obtained for a 35-deg hyperboloid and for a 45-deg spherically blunted cone using a time-dependent finite-difference method. The radiative heating rates for the coupled laminar flow are significantly reduced as compared to the corresponding no-blowing case; however, for the coupled turbulent flow, it is found that the surface radiative heating rates are substantially increased and often exceed the corresponding no-blowing values. Turbulence is found to have no effect on the surface radiative heating rates for the no-blowing solutions. The present results are compared with the other available solutions and some additional solutions are presented.

Kumar, A.

A simple numerical orthogonal coordinate generator for fluid dynamic applications

An application of a simple numerical technique which allows for the rapid construction of orthogonal coordinate systems about two dimensional and axisymmetric bodies is presented. The technique is based on a predictor corrector numerical method. It can be used to generate orthogonal meshes which have unequally spaced points in two directions. These orthogonal meshes in their transformed computational plane are, however, equally spaced so that the differencing for the metric coefficients and the fluid dynamic equation terms can be easily determined using equally spaced central finite differences. Solutions to the Navier-Stokes equations for flow over blunt bodies with reverse curvature are presented. The coupling of the time dependent fluid dynamic equations and the coordinate generator worked well with no undesirable effects noted.

Graves, R. A., Jr.

A vectorized code for calculating laminar and turbulent hypersonic flows about blunt axisymmetric bodies at zero and small angles of attack

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.

Kumar, A.

Laminar and turbulent flows over a spherically blunted cone with massive surface blowing

Numerical solutions are presented for the flow over a spherically blunted cone with massive surface blowing. Time-dependent viscous shock layer equations are used to describe the flowfield. The boundary conditions on the body surface include a prescribed blowing rate distribution. The governing equations are solved by a time-asymptotic finite-difference method. Results presented here are only for a perfect gas-type flow at zero angle of attack. Both laminar and turbulent flow solutions are obtained. It is found that the surface blowing smooths out the effect of the curvature discontinuity at the sphere-cone juncture point on the laminar flowfield and results in a negative pressure gradient over the body. The shock slope increases on the downstream portion of the body as the surface blowing rate is increased. The turbulent flow with surface blowing is found to redevelop a boundary-layer-like region near the surface. The effects of this boundary-layer-like region on the flowfield and heating rates are discussed.

Kumar, A.

Application of a numerical orthogonal coordinate generator to axisymmetric blunt bodies

An application of a simple numerical technique which allows for the rapid construction of orthogonal coordinate systems about axisymmetric blunt bodies is presented. This technique can generate orthogonal meshes which have unequally spaced points in two directions. Relations are given for the numerical generation of the metric coefficients. Body shapes ranging from simple analytical bodies to complex reverse curvature bodies are presented together with their orthogonal coordinate systems. The relatively good accuracy of the technique is shown in tabular data describing coordinate line slopes and metric coefficients. The predictor-corrector numerical method used to generate these results is both simple in concept and easy to program, so that the application of the technique should be broader than the results presented.

Graves, R. A., Jr.

Comparative study of the convergence rates of two numerical techniques

The paper examines the applicability of the three-step Stetter (1968) method to the problem of hypersonic viscous flow over a blunt axisymmetric body used for planetary entry probes at zero angle of attack. The flow-field results using the two-step finite-difference MacCormack (1969) method are reported by Kumar and Graves (1977). Only the computational efficiency of Stetter's method is compared with that of MacCormack's in terms of the iterative time steps and computing time required for the steady-state solution. Advantages of Stetter's method over MacCormack's are established.

Kumar, A.

Navier-Stokes solutions using Stetter's method

Compressible Navier-Stokes equations for quasi-one-dimensional flow in a converging-diverging nozzle have been solved using Stetter's three-step predictor-corrector technique. Particular emphasis is given to the minimum iterative step feature for steady-state solutions. It is found that for the nonviscous-dominated case, Stetter's method attains the steady-state solution in the fewest steps when compared to four other currently used techniques.

Graves, R. A., Jr.

Effects of small angle of attack on the radiating viscous shock layer solutions for Jovian entry

Results are presented for the laminar and turbulent flows of a reacting mixture of hydrogen-helium in chemical equilibrium over blunted cones of large half angles at small angles of attack. The freestream conditions correspond to a typical Jovian entry trajectory point. Radiative heat flux is calculated by using a realistic step-function model for the absorption coefficient. An investigation is also made to determine whether a blunted cone at small angles of attack can be replaced by equivalent bodies at zero angle of attack for the purpose of calculating the shock shapes, surface pressure, skin friction and heating rates.

Kumar, A.

Turbulent viscous shock layer solutions for Jovian entry at small angles of attack

The equations governing the laminar and turbulent flows of reacting gas mixtures in chemical equilibrium over axially symmetric blunt bodies at small angles of attack are developed and presented in the unsteady conservative form. Solutions are obtained in the planes of symmetry of the flow field for the typical conditions encountered by a probe entering the Jovian atmosphere at small angles of attack. The eddy-viscosity is approximated by a two-layer model. The shock wave is treated as a discontinuity across which the shock relations are used to compute the flow conditions behind the shock. A time-asymptotic finite-difference method is used to solve the equations. The zero angle of attack results for a hyperboloid are compared with the existing results but neither experimental nor computational results are available for comparison with the present results at angle of attack.

Jumar, A.

Numerical solution of the viscous hypersonic flow past blunted cones at angle of attack

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.

Kumar, A.