Pressure and heat transfer on blunt curved plates with concave and convex surfaces in hypersonic flow
Pressure and heat transfer on blunt concave and convex plates in hypersonic 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.
Pressure and heat transfer on blunt concave and convex plates in hypersonic flow
Heat input influence on hypersonic flows past sphere and circular cylinder
For hypersonic flow with a shock wave, there is a similar solution consistent throughout the viscous and inviscid layers along a very slender three-fourths-power body of revolution The strong pressure interaction problem can then be treated by the method of similarity. Numerical calculations are performed in the viscous region with the edge pressure distribution known from the inviscid similar solutions. The compressible laminar boundary-layer equations are transformed into a system of ordinary differential equations. The resulting two-point boundary value problem is then solved by the Runge-Kutta method with a modified Newton's method for the corresponding boundary conditions. The effects of wall temperature, mass bleeding, and body transverse curvature are investigated. The induced pressure, displacement thickness, skin friction, and heat transfer due to the previously mentioned parameters are estimated and analyzed.
A direct solution approach for surface erosion in particle-laden hypersonic flows is extended for use in low-cost two-way coupled solutions of dilute gas-particle flows. The Trajectory Control Volume method, which uses a sparse set of probe particles to predict surface erosion distributions on general vehicles, is reformulated for the solution of source terms by mean trajectory subdivision and computing a flux differencing. The approach is verified successfully against a boundary layer solution and shown to agree well with experimental measurements. A representative Mars entry case, with conditions and geometry based on the ExoMars Schiaparelli capsule, is solved with the approach to study the impact of two-way coupling on surface heating and erosion. Results indicate that for realistic loading conditions, heating is largely unmodified compared to one-way coupled results at peak heating trajectory conditions, and no measureable difference is observed in the surface erosion rate. At exaggerated loading conditions high enough to observe coupling effects, the worst case collisional heating can increase heating by up to 60%.
A comprehensive study for the inviscid numerical calculation of the hypersonic flow past a class of elliptic-cone derived waveriders is presented. The theoretical background associated with hypersonic small-disturbance theory (HSDT) is reviewed. Several approximation formulas for the waverider compression surface are established. A CFD algorithm is used to calculate flow fields for the on-design case and a variety of off-design cases. The results are compared with HSDT, experiment, and other available CFD results. For the waverider shape used in previous investigations, the bow shock for the on-design condition stands off from the leading-edge tip of the waverider. It was found that this occurs because the tip was too thick according to the approximating shape formula that was used to describe the compression surface. When this was corrected, the bow shock became closer to attached as it should be. At Mach numbers greater than the design condition, a lambda-shock configuration develops near the tip of the compression surface. At negative angles of attack, other complicated shock patterns occur near the leading-edge tip. These heretofore unknown flow patterns show the power and utility of CFD for investigating novel hypersonic configurations such as waveriders.
Vortex-induced heating alleviation to lee side of slender wings in hypersonic flow by contouring leading edge planform
Shock tunnel drag measurements on sharp slender cones in near free molecule hypersonic flow in air and He
This report describes the research work undertaken at the Massachusetts Institute of Technology. The aim of this research is to identify effective algorithms and methodologies for the efficient and routine solution of hypersonic viscous flows about re-entry vehicles. For over ten years we have received support from NASA to develop unstructured mesh methods for Computational Fluid Dynamics. As a result of this effort a methodology based on the use, of unstructured adapted meshes of tetrahedra and finite volume flow solvers has been developed. A number of gridding algorithms flow solvers, and adaptive strategies have been proposed. The most successful algorithms developed from the basis of the unstructured mesh system FELISA. The FELISA system has been extensively for the analysis of transonic and hypersonic flows about complete vehicle configurations. The system is highly automatic and allows for the routine aerodynamic analysis of complex configurations starting from CAD data. The code has been parallelized and utilizes efficient solution algorithms. For hypersonic flows, a version of the, code which incorporates real gas effects, has been produced. One of the latest developments before the start of this grant was to extend the system to include viscous effects. This required the development of viscous generators, capable of generating the anisotropic grids required to represent boundary layers, and viscous flow solvers. In figures I and 2, we show some sample hypersonic viscous computations using the developed viscous generators and solvers. Although these initial results were encouraging, it became apparent that in order to develop a fully functional capability for viscous flows, several advances in gridding, solution accuracy, robustness and efficiency were required. As part of this research we have developed: 1) automatic meshing techniques and the corresponding computer codes have been delivered to NASA and implemented into the GridEx system, 2) a finite element algorithm for the solution of the viscous compressible flow equations which can solve flows all the way down to the incompressible limit and that can use higher order (quadratic) approximations leading to highly accurate answers, and 3) and iterative algebraic multigrid solution techniques.
The meaning of catalysis and its relation to aerodynamic heating in nonequilibrium hypersonic flows are discussed. The species equations are described and boundary conditions for them are derived for a multicomponent gas and for a binary gas. Slip effects are included for application of continuum methods to low-density flows. Measurement techniques for determining catalytic wall recombination rates are discussed. Among them are experiments carried out in arc jets as well as flow reactors. Diagnostic methods for determining the atom or molecule concentrations in the flow are included. Results are given for a number of materials of interest to the aerospace community, including glassy coatings such as the RCG coating of the Space Shuttle and for high temperature refractory metals such as coated niobium. Methods of calculating the heat flux to space vehicles in nonequilibrium flows are described. These methods are applied to the Space Shuttle, the planned Aeroassist Flight Experiment, and a hypersonic slender vehicle such as a transatmospheric vehicle.
Results of an experimental and numerical investigation of tangential swept slot injection (sweep angles of 22.5 and 45 deg) into a thick turbulent boundary layer at Mach 6 are presented. Film cooling effectiveness, skin friction, and flow structure downstream of the swept slot injection are investigated. The data are compared to that for unswept slots, and it is found that cooling effectiveness and skin-friction reductions are not significantly affected by sweeping the slot. Predictions of cooling effectiveness and skin friction obtained by a numerical finite-difference technique agree reasonably well with experimental surface variables. As in previous supersonic two-dimensional slot research, reduced mixing was found downstream of the slot lip in the present three-dimensional case.
Aerodynamic heating in hypersonic space vehicles is an important factor to be considered in their design. Therefore the designers of such vehicles need reliable heat transfer data in this respect for a successful design. Such data is usually produced by testing the models of hypersonic surfaces in wind tunnels. Most of the hypersonic test facilities at present are conventional blow-down tunnels whose run times are of the order of several seconds. The surface temperatures on such models are obtained using standard techniques such as thin-film resistance gages, thin-skin transient calorimeter gages and coaxial thermocouple or video acquisition systems such as phosphor thermography and infrared thermography. The data are usually reduced assuming that the model behaves like a semi-infinite solid (SIS) with constant properties and that heat transfer is by one-dimensional conduction only. This simplifying assumption may be valid in cases where models are thick, run-times short, and thermal diffusivities small. In many instances, however, when these conditions are not met, the assumption may lead to significant errors in the heat transfer results. The purpose of the present paper is to investigate this aspect. Specifically, the objectives are as follows: (1) to determine the limiting conditions under which a model can be considered a semi-infinite body; (2) to estimate the extent of errors involved in the reduction of the data if the models violate the assumption; and (3) to come up with correlation factors which when multiplied by the results obtained under the SIS assumption will provide the results under the actual conditions.
A differential equation governing the geometry of a two-dimensional ballute in hypersonic flow and its constraining boundary conditions are derived under idealized assumptions. By solving these equations, the shape of the ballute is determined over a range of conditions. Lift, drag, pitching moment, and the allowed limit of center-of-gravity location for stability (meta-center) are then calculated using Newtonian hypersonic approximation. It is shown that the meta-center occurs near the forward end because of compliance of the ballute membrane to the shock layer pressures, especially at low free-stream densities. In order for the vehicle employing the ballute to be stable at all densities, the center of gravity must be within approximately the forward 20 percent of overall length of the vehicle. However, typical flight trajectories of an aeroassisted orbital transfer vehicle employing the ballute for aerobraking show that the vehicle may be able to complete its atmospheric flight without tumbling provided that the center of gravity is located within the forward 43 percent of the vehicle length because of the relatively short duration of flight through the destabilizing low-density regime.
Steady hypersonic flow over blunt bodies with high rates of continuous mass and heat transfer from body surface
Hypersonic flows over straight and bent biconics are calculated for a range of freestream conditions in which the gas behind the shock is treated as either perfect or real. The Parabolized-Navier-Stokes (PNS) equations form the basis of the approximation scheme. Good comparisons with experimental data for pressures, forces and moments, heat transfer, and oil-flow patterns serve to validate the perfect-gas version of the code. Circumferential velocity vector plots further aid in the interpretation of leeside oil-flow patterns. A variable-effective-gamma (VEG) option is implemented for the real-gas calculations. Gamma, now defined as the ratio of enthalpy to internal energy, is determined from a locally valid linear relation in enthalpy and pressure at every mesh point which in turn is calculated from a benchmark equilibrium code. The VEG option is easily incorporated into a host PNS code because it uses the underlying perfect-gas structure of the code. Comparisons to experimental data for heat transfer and shock shape in high enthalpy air have been obtained using the VEG option.
Results of a numerical study of low-density hypersonic flow about cylindrically blunted wedges and spherically blunted cones with body half angles of 0, 5, and 10 deg are presented. Most of the transitional flow regime encountered during entry between the free molecule and continuum regimes is simulated for a reentry velocity of 7.5 km/s by including freestream conditions of 70 to 100 km. The bodies are at zero angle of incidence and have diffuse and finite catalytic surfaces. Translational, thermodynamic, and chemical nonequilibrium effects are considered in the numerical simulation by utilizing the direct simulation Monte Carlo (DSMC) method. The numerical simulations show that noncontinuum effects such as surface temperature jump, and velocity slip are evident for all cases considered. The onset of chemical dissociation occurs at a simulated altitude of 96 km for the two-dimensional configurations. Comparisons between the DSMC and continuum viscous shock-layer calculations highlight the significant difference in flowfield structure predicted by the two methods.
Implicit, finite-difference, shock-capturing algorithm calculates inviscid, hypersonic flows in chemical equilibrium. Implicit formulation chosen because overcomes limitation on mathematical stability encountered in explicit formulations. For dynamical portion of problem, Euler equations written in conservation-law form in Cartesian coordinate system for two-dimensional or axisymmetric flow. For chemical portion of problem, equilibrium state of gas at each point in computational grid determined by minimizing local Gibbs free energy, subject to local conservation of molecules, atoms, ions, and total enthalpy. Major advantage: resulting algorithm naturally stable and captures strong shocks without help of artificial-dissipation terms to damp out spurious numerical oscillations.
Report discusses mathematical models of turbulence used in numerical simulations of complicated viscous, hypersonic flows. Includes survey of essential features of models and their statuses in applications.
An approximate solution technique has been developed for three-dimensional, inviscid, hypersonic flows. The method uses Maslen's explicit pressure equation and the assumption of approximate stream surfaces in the shock layer. This approximation represents a simplification of Maslen's asymmetric method. The solution procedure involves iteratively changing the shock shape in the subsonic-transonic region until the correct body shape is obtained. Beyond this region, the shock surface is determined by using a marching procedure. Results are presented herein for a paraboloid and elliptic cone at angle of attack. Calculated surface pressure distributions, shock shapes, and property profiles are compared with experimental data and finite-difference solutions of the Euler equations. Comparisons of the results of the present method with experimental data and detailed predictions are very good. Since the present method provides a very rapid computational procedure, it can be used for parametric or preliminary design applications. One useful application would be to incorporate a heating procedure for aerothermal studies.