Theoretical and Experimental Investigation of the Effect of Yaw on Heat Transfer to Circular Cylinders in Hypersonic Flow
Explore the source record for details and available documents.
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.
Explore the source record for details and available documents.
A procedure for calculating three-dimensional steady and nonsteady supersonic flows with the method of characteristics is developed and discussed. An approximate method is deduced from the characteristics method and shown to be of practical value at high supersonic speeds.
The quality of reentry simulation for Shuttle, HERMES, Sanger, and NASP systematically suffers from the strong non-equilibrium of rotational and vibrational temperature due to the rapid acceleration of the test gas in the nozzle. Therefore the determination of temperatures is necessary and, if possible, preferable by a non-intrusive technique. The specific interests of this review are optical techniques such as electron beam fluorescence, laser-induced fluorescence, and coherent anti-Stokes Raman scattering. The capabilities available for local measurements with temporal resolution and quantitative accuracy are discussed for velocity, temperature, density, species concentrations, and fluctuations due to turbulence. The applicability of these methods of measurement is presented and discussed for the coming topic in aerothermodynamics: experimental techniques of hot gases in high enthalpy flows.
A roadmap for computational fluid dynamics (CFD) code validation is developed. The elements of the roadmap are consistent with air-breathing vehicle design requirements and related to the important flow path components: forebody, inlet, combustor, and nozzle. Building block and benchmark validation experiments are identified along with their test conditions and measurements. Based on an evaluation criteria, recommendations for an initial CFD validation data base are given and gaps identified where future experiments would provide the needed validation data.
Critical issues concerning the modeling of low density hypervelocity flows where thermochemical nonequilibrium effects are pronounced are discussed. Emphasis is on the development of validated analysis tools, and the activity in the NASA Ames Research Center's Aerothermodynamics Branch is described. Inherent in the process is a strong synergism between ground test and real gas computational fluid dynamics (CFD). Approaches to develop and/or enhance phenomenological models and incorporate them into computational flowfield simulation codes are discussed. These models were partially validated with experimental data for flows where the gas temperature is raised (compressive flows). Expanding flows, where temperatures drop, however, exhibit somewhat different behavior. Experimental data for these expanding flow conditions is sparse and reliance must be made on intuition and guidance from computational chemistry to model transport processes under these conditions. Ground based experimental studies used to provide necessary data for model development and validation are described. Included are the performance characteristics of high enthalpy flow facilities, such as shock tubes and ballistic ranges.
The development of an algorithm for the solution of the compressible Euler equations at high Mach numbers on unstructured tetrahedral meshes is described. The basic algorithm is constructed in the form of a central difference scheme plus an explicit added artificial viscosity based upon fourth order differences of the solution. The stability of the solution in the vicinity of strong gradients is preserved by the incorporation of an additional artificial viscosity based upon a second order difference. Higher order accuracy is regained by using the ideas of flux corrected transport to limit the amount of added viscosity. The solution is advanced to steady state by means of an explicit multi-stage time-stepping method. The computational efficiency of the complete process is improved by incorporating an unstructured multigrid acceleration procedure. A number of flows of practical interest are analyzed to demonstrate the numerical performance of the proposed approach.
A detailed comparison of numerical results obtained by solving the Euler equations for the inviscid flow over delta wings in reentry configurations is presented. The investigation involves a side-by-side comparison between independently developed upwind Euler solvers at VKI (M3D) and NASA Johnson Space Center (E3D) using identical grids. In both solvers the governing equations are integrated by means of time marching finite volume shock capturing methods, based on a cell centered upwind evaluation of the cell face fluxes and nonlinear limiters. High resolution schemes are obtained via MUSCL characteristic variable extrapolation to ensure total variation diminishing (TVD) properties and therefore monotonic discontinuity capturing. Further comparisons are made with several published results by other authors based on both upwind and central discretizations. The present results contribute in asserting the high resolution upwind TVD schemes as the most reliable numerical technique to handle the strong discontinuities typical of high speed flows. The performances of the two upwind solvers are satisfactory and the numerical results in good agreement. However, the important issue of reaching a grid converged solution for the present complex three-dimensional inviscid flow problems was not achieved.
Three solution algorithms, explicit underrelaxation, point implicit, and lower upper symmetric Gauss-Seidel (LUSGS), are used to compute nonequilibrium flow around the Apollo 4 return capsule at 62 km altitude. By varying the Mach number, the efficiency and robustness of the solution algorithms were tested for different levels of chemical stiffness. The performance of the solution algorithms degraded as the Mach number and stiffness of the flow increased. At Mach 15, 23, and 30, the LUSGS method produces an eight order of magnitude drop in the L2 norm of the energy residual in 1/3 to 1/2 the Cray C-90 computer time as compared to the point implicit and explicit under-relaxation methods. The explicit under-relaxation algorithm experienced convergence difficulties at Mach 23 and above. At Mach 40 the performance of the LUSGS algorithm deteriorates to the point it is out-performed by the point implicit method. The effects of the viscous terms are investigated. Grid dependency questions are explored.
Three solution algorithms, explicit under-relaxation, point implicit, and lower-upper symmetric Gauss-Seidel, are used to compute nonequilibrium flow around the Apollo 4 return capsule at the 62-km altitude point in its descent trajectory. By varying the Mach number, the efficiency and robustness of the solution algorithms were tested for different levels of chemical stiffness.The performance of the solution algorithms degraded as the Mach number and stiffness of the flow increased. At Mach 15 and 30, the lower-upper symmetric Gauss-Seidel method produces an eight order of magnitude drop in the energy residual in one-third to one-half the Cray C-90 computer time as compared to the point implicit and explicit under-relaxation methods. The explicit under-relaxation algorithm experienced convergence difficulties at Mach 30 and above. At Mach 40 the performance of the lower-upper symmetric Gauss-Seidel algorithm deteriorates to the point that it is out performed by the point implicit method. The effects of the viscous terms are investigated. Grid dependency questions are explored.
An investigation of the effects of variation of leading-edge sweep and surface inclination on the flow over blunt flat plates was conducted at Mach numbers of 4 and 5.7 at free-stream Reynolds numbers per inch of 6,600 and 20,000, respectively. Surface pressures were measured on a flat plate blunted by a semicylindrical leading edge over a range of sweep angles from 0 deg to 60 deg and a range of surface inclinations from -10 deg to +10 deg. The surface pressures were predicted within an average error of +/- 8 percent by a combination of blast-wave and boundary-layer theory extended herein to include effects of sweep and surface inclination. This combination applied equally well to similar data of other investigations. The local Reynolds number per inch was found to be lower than the free-stream Reynolds number per inch. The reduction in local Reynolds number was mitigated by increasing the sweep of the leading edge. Boundary-layer thickness and shock-wave shape were changed little by the sweep of the leading edge.
An uncoupled Computational Fluid Dynamics-Direct Simulation Monte Carlo (CFD-DSMC) technique is developed and applied to provide solutions for continuum jets interacting with rarefied external flows. The technique is based on a correlation of the appropriate Bird breakdown parameter for a transitional-rarefied condition that defines a surface within which the continuum solution is unaffected by the external flow-jet interaction. The method is applied to two problems to assess and demonstrate its validity; one of a jet interaction in the transitional-rarefied flow regime and the other in the moderately rarefied regime. Results show that the appropriate Bird breakdown surface for uncoupling the continuum and non-continuum solutions is a function of a non-dimensional parameter relating the momentum flux and collisionality between the two interacting flows. The correlation is exploited for the simulation of a jet interaction modeled for an experimental condition in the transitional-rarefied flow regime and the validity of the correlation is demonstrated. The uncoupled technique is also applied to an aerobraking flight condition for the Mars Global Surveyor spacecraft with attitude control system jet interaction. Aerodynamic yawing moment coefficients for cases without and with jet interaction at various angles-of-attack were predicted, and results from the present method compare well with values published previously. The flow field and surface properties are analyzed in some detail to describe the mechanism by which the jet interaction affects the aerodynamics.
A new technique for the control of grid line spacing and intersection angles of a structured volume grid, using elliptic partial differential equations (PDEs) is presented. Existing structured grid generation algorithms make use of source term hybridization to provide control of grid lines, imposing orthogonality implicitly at the boundary and explicitly on the interior of the domain. A bridging function between the two types of grid line control is typically used to blend the different orthogonality formulations. It is shown that utilizing such a bridging function with source term hybridization can result in the excessive use of computational resources and diminishes robustness. A new approach, Anisotropic Lagrange Based Trans-Finite Interpolation (ALBTFI), is offered as a replacement to source term hybridization. The ALBTFI technique captures the essence of the desired grid controls while improving the convergence rate of the elliptic PDEs when compared with source term hybridization. Grid generation on a blunt cone and a Shuttle Orbiter is used to demonstrate and assess the ALBTFI technique, which is shown to be as much as 50% faster, more robust, and produces higher quality grids than source term hybridization.
A structured-grid quality measure is proposed, combining three traditional measurements: intersection angles, stretching, and curvature. Quality assesses whether the grid generated provides the best possible tradeoffs in grid stretching and skewness that enable accurate flow predictions, whereas the grid density is assumed to be a constraint imposed by the available computational resources and the desired resolution of the flow field. The usefulness of this quality measure is assessed by comparing heat transfer predictions from grid convergence studies for grids of varying quality in the range of [0.6-0.8] on an 8'half-angle sphere-cone, at laminar, perfect gas, Mach 10 wind tunnel conditions.
This paper presents computational results obtained with the direct simulation Monte Carlo (DSMC) method for several biconic test cases in which shock interactions and flow separation-reattachment are key features of the flow. Recent ground-based experiments have been performed for several biconic configurations, and surface heating rate and pressure measurements have been proposed for code validation studies. The present focus is to expand on the current validating activities for a relatively new DSMC code called DS2V that Bird (second author) has developed. Comparisons with experiments and other computations help clarify the agreement currently being achieved between computations and experiments and to identify the range of measurement variability of the proposed validation data when benchmarked with respect to the current computations. For the test cases with significant vibrational nonequilibrium, the effect of the vibrational energy surface accommodation on heating and other quantities is demonstrated.
A new procedure seeks to combine the thin-layer Navier-Stokes solver LAURA with the parabolized Navier-Stokes solver UPS for the aerothermodynamic solution of chemically-reacting air flow fields. The interface protocol is presented and the method is applied to two slender, blunted shapes. Both axisymmetric and three-dimensional solutions are included with surface pressure and heat transfer comparisons between the present method and previously published results. The case of Mach 25 flow over an axisymmetric six degree sphere-cone with a non-catalytic wall is considered to 100 nose radii. A stability bound on the marching step size was observed with this case and is attributed to chemistry effects resulting from the non-catalytic wall boundary condition. A second case with Mach 28 flow over a sphere-cone-cylinder-flare configuration is computed at both two and five degree angles of attack with a fully-catalytic wall. Surface pressures are seen to be within five percent with the present method compared to the baseline LAURA solution and heat transfers are within 10 percent. The effect of grid resolution is investigated in both the radial and streamwise directions. The procedure demonstrates significant, order of magnitude reductions in solution time and required memory for the three-dimensional case in comparison to an all thin-layer Navier-Stokes solution.
A capability to calculate surface heating rates has been incorporated in an approximate three-dimensional inviscid technique. Surface streamlines are calculated from the inviscid solution, and the axisymmetric analog is then used along with a set of approximate convective-heating equations to compute the surface heat transfer. The method is applied to blunted axisymmetric and three-dimensional ellipsoidal cones at angle of attack for the laminar flow of a perfect gas. The method is also applicable to turbulent and equilibrium-air conditions. The present technique predicts surface heating rates that compare favorably with experimental (ground-test and flight) data and numerical solutions of the Navier-Stokes (NS) and viscous shock-layer (VSL) equations. The new technique represents a significant improvement over current engineering aerothermal methods with only a modest increase in computational effort.
This paper considers the streamline-upwind Petrov/Galerkin (SUPG) method applied to the compressible Euler and Navier-Stokes equations in conservation-variable form. The spatial discretization, including a modified approach for interpolating the inviscid flux terms in the SUPG finite element formulation, is briefly reviewed. Of particular interest is the behavior of the shock capturing operator, which is required to regularize the scheme in the presence of strong, shock-induced gradients. A standard shock capturing operator which has been widely used in previous studies by several authors is presented and discussed. Specific modifications are then made to this standard operator which are designed to produce a more physically consistent discretization in the presence of strong shock waves. The actual implementation of the term in a finite dimensional approximation is also discussed. The behavior of the standard and modified scheme is then compared for several supersonic/hypersonic flows. The modified shock capturing operator is found to preserve enthalpy in the inviscid portion of the flowfield substantially better than the standard operator.
Presentation topics include background and motivation; physical modeling including governing equations and thermochemistry; finite element formulation; results of inviscid thermal nonequilibrium chemically reacting flow and viscous thermal equilibrium chemical reacting flow; and near-term effort.