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.
This paper describes the development and testing of a blunt entry body wth a high degree of transonic dynamic stability and one hypersonic trim point. Because of the inadequacy of CFD techniques to predict dynamic characteristics of blunt vehicles in the transonic flight regime, a series of tests was performed.
A relatively simple method was developed previously (authors, 1973) for calculating laminar, transitional, and turbulent heating rates on three-dimensional bodies in hypersonic flows. This method was shown to yield reasonably accurate results for laminar heating on blunted circular and elliptical cones and an earlier version of the space shuttle vehicle. As the boundary layer along the surface grows, more and more of the inviscid-flow mass is entrained into the boundary layer, and the streamlines which passed through the nearly normal portion of the bow shock wave are 'swallowed' by the boundary layer. This phenomenon is often referred to as entropy-layer or streamline swallowing, and it can have a significant effect on the calculated heating rates. An approximate, yet simple, method for including the effects of entropy-layer swallowing in the heating-rate calculations is given.
A description is given of the computer code for predicting viscous shock-layer flows over nonanalytic blunt bodies (Program VISLNABB) for hypersonic, low Reynolds number flows. Four specific and one general body geometries are considered. In addition to sphere-cones, cylinder wedges and geometries defined in tabular form, options for hyperboloids and paraboloids are included. Details of the theory and results are included in a separate engineering report. The program, subroutines, variables in common, and input and output data are described. Listings of the program code, output data for a sample case, and the input data for this sample case are included.
Supersonic flow about cone calculated accurately. Report presents theoretical basis of computer code solving parabolized Navier-Stokes equations of supersonic and hypersonic flow. For increased accuracy in resolution of details of strong aerodynamic shocks, code incorporates implicit, finite-volume, upwind numerical-integration scheme. Performs well in numerical simulations of flows around simple bodies.
An integrated trajectory/control analysis algorithm has been used to generate trajectories and desired control strategies for two different hypersonic air-breathing vehicle models and orbit targets. Both models used cubic spline curve fit tabulated winged-cone accelerator vehicle representations. Near-fuel-optimal, horizontal takeoff trajectories, imposing a dynamic pressure limit of 1000 psf, were developed. The first model analysis case involved a polar orbit and included the dynamic effects of using elevons to maintain longitudinal trim. Analysis results indicated problems with the adequacy of the propulsion model and highlighted dynamic pressure/altitude instabilities when using vehicle angle of attack as a control variable. Also, the magnitude of computed elevon deflections to maintain trim suggested a need for alternative pitch moment management strategies. The second analysis case was reformulated to use vehicle pitch attitude relative to the local vertical as the control variable. A new, more realistic, air-breathing propulsion model was incorporated. Pitch trim calculations were dropped and an equatorial orbit was specified. Changes in flight characteristics due to the new propulsion model have been identified. Flight regimes demanding rapid attitude changes have been noted. Also, some issues that would affect design of closed-loop controllers were ascertained.
An experimental investigation of distributed sand-grain surface roughness effects on hypersonic boundary-layer transition and turbulent heating has been performed. Two representative entry vehicle geometries, a sphere-cone aeroshell and a spherical-cap aeroshell, were considered. Cast ceramic models of each geometry were fabricated with distributed sand-grain roughness patterns of different heights that simulated an ablated thermal protection system. Wind tunnel testing was performed at Mach 6 over a range of Reynolds numbers sufficient to produce laminar, transitional, and turbulent flow. Aeroheating and boundary-layer transition onset data were obtained using global phosphor thermography. The experimental heating data are presented herein, as are comparisons to laminar and turbulent smooth-wall heat transfer distributions from computational flow field simulations. The boundary-layer transition data were found to correlate with a functional representation developed in prior roughness studies, although the data scatter was greater owing to the height variability of distributed sand-grain roughness.
A class of vehicles for a mission to Mars are analyzed for aerodynamic characteristics using advanced Computational Fluid Dynamics (CFD). The general configuration is a modified cone-conical-frustum geometry where the nose radius has a large influence on the flowfield. Inviscid-compressible flow using the Euler equations and viscous-compressible flow using the thin-layer Navier-Stokes equations is applied to the configuration. The surface modeling, grid generation and application of state-of-the-art CFD software are described. The effects of nose radius, angle of attack, and hypersonic velocity on the flight characteristics of the vehicle are discussed. The numerical simulations demonstrate the merits of the inviscid and viscous software. Results are compared with wind tunnel experiments.
In its present state, the Amplification Factor Transport (AFT) model provides an efficient alternative for high-fidelity stability methods for complex geometries subject to low-speed conditions considered in the industrial-production CFD environment. Its performance has not yet been documented for the typical transition scenarios in hypersonic conditions. The goal of the present study is to apply the AFT model to a yawed- and flared-cone configuration, for which the laminar flows respectively support the crossflow and second-mode instability mechanisms, and to compare the amplification results against high-fidelity results computed with Linear Parabolized Stability Equations (LPSE). High-resolution viscous-flow solutions are provided with OVERFLOW and COFFE. For the yawed-cone case, the AFT model yields the largest amplification in a region of the flow where the expected crossflow instability is not amplified (in the leeward symmetry plane) according to LPSE. For the flared-cone case with a wall-temperature slightly lower than the adiabatic value, the AFT model yields extremely large amplification factors (in excess of 100), exceeding the computed LPSE amplification of second-mode instability results by a factor 14. Upon considering a cold-wall condition, the model shows much smaller amplification factors, which is the opposite with respect to the expected behavior of the second mode. The AFT model does not adequately capture the crossflow and second-mode instability mechanisms in its present state.
Measurements of waves initiated by wind-tunnel sounds are obtained for both supersonic and hypersonic speed. Flat-plate experiments carried out at Mach numbers of 4.5, 3.7, and 2.4 under the condition of laminar tunnel-wall boundary layers whereby a low background fluctuation level prevailed are outlined, as well as flat-plate experiments with sound-induced waves at supersonic speeds. Experimental studies of 'natural' wave growth on sharp cones at zero angle of attack and Mach numbers near 8 are discussed. Attention is focused on experimental research needed on boundary-layer receptivity to eddy Mach-wave radiation, hypersonic speed, and the stability-transition relationship in zero-pressure-gradient layers.
Two new parabolized Navier-Stokes (PNS) codes were developed to compute the three-dimensional, viscous, chemically reacting flow of air around hypersonic vehicles such as the National Aero-Space Plane (NASP). The first code (TONIC) solves the gas dynamic and species conservation equations in a fully coupled manner using an implicit, approximately-factored, central-difference algorithm. This code was upgraded to include shock fitting and the capability of computing the flow around complex body shapes. The revised TONIC code was validated by computing the chemically-reacting (M(sub infinity) = 25.3) flow around a 10 deg half-angle cone at various angles of attack and the Ames All-Body model at 0 deg angle of attack. The results of these calculations were in good agreement with the results from the UPS code. One of the major drawbacks of the TONIC code is that the central-differencing of fluxes across interior flowfield discontinuities tends to introduce errors into the solution in the form of local flow property oscillations. The second code (UPS), originally developed for a perfect gas, has been extended to permit either perfect gas, equilibrium air, or nonequilibrium air computations. The code solves the PNS equations using a finite-volume, upwind TVD method based on Roe's approximate Riemann solver that was modified to account for real gas effects. The dissipation term associated with this algorithm is sufficiently adaptive to flow conditions that, even when attempting to capture very strong shock waves, no additional smoothing is required. For nonequilibrium calculations, the code solves the fluid dynamic and species continuity equations in a loosely-coupled manner. This code was used to calculate the hypersonic, laminar flow of chemically reacting air over cones at various angles of attack. In addition, the flow around the McDonnel Douglas generic option blended-wing-body was computed and comparisons were made between the perfect gas, equilibrium air, and the nonequilibrium air results.
Viscous and bluntness induced pressures for flat plates, wedges, cylinders and cones at Mach 41
Boundary-layer receptivity and stability of Mach 6 flow over smooth and rough 7 half-angle sharp-tipped cones are numerically investigated. The receptivity of the boundary layer to slow acoustic disturbances, fast acoustic disturbances, and vortical disturbances are considered. The effects of two-dimensional isolated and distributed roughness on the receptivity and stability are also simulated. The results show that the instability waves are generated in the leading edge region and that the boundary layer is much more receptive to slow acoustic waves than to the fast waves. Vortical disturbances also generate unstable second modes, however the receptivity coefficients are smaller than that of the slow acoustic wave. An isolated two-dimensional roughness element of height h/delta =1/4 did not produce any difference in the receptivity or in the stability of the boundary layer. Distributed roughness elements produced a small decrease in the receptivity coefficient and also stabilized the boundary layer by small amounts.
Phenolic Impregnated Carbon Ablator (PICA) gained heritage during the stardust mission and is the baseline material for missions to Mars. The phenolic phase of PICA is friable, and therefore the surface is treated with a siloxane resin (NuSil) to mitigate the shedding of particulate matter. Historically, material response models have not accounted for the presence of this coating. Therefore, a test campaign was executed at the Hypersonic Materials Environmental Test System (HyMETS) to investigate the performance of NuSil under multiple heating rates and gas composition. Sphere-cone models were instrumented with thermocouples to measure the in-depth temperature response, and surface temperatures were monitored with a pyrometer. Furthermore, emission spectrometers were used to collect spectra from species emitting at the surface of the test article and the post-shock stagnation region. Post-test analysis suggests that the coating decomposes according to multiple stages. Inspection of silicon emission data reveals a rapid rise and decline in intensity within the first couple of seconds of testing, which is attributed to pyrolysis of NuSil. Pyrolysis is followed by a stagnation in silicon emission and surface temperature which are attributed to the formation of a thin layer of oxidation-resistant silicon oxycarbide. Moreover, the presence of the silicon oxycarbide layer significantly reduces the in-depth thermal response when compared to the baseline material. Eventually, the silicon oxycarbide layer decomposes via high-temperature reactions, which expose the underlying char layer to a highly reactive flow.
Phenolic Impregnated Carbon Ablator (PICA) is NASA’s baseline thermal protection system for missions to Mars. PICA is friable; therefore, NuSil ® CV-1144-0 is applied to the surface of flight hardware to mitigate the contamination of spacecraft sensors. NuSil is composed of a polysiloxane resin system that can transform into an oxidation-resistant coating upon heating. Therefore, a pathfinder campaign was conducted at the Hypersonic Materials Environmental Test System (HyMETS) at the NASA Langley Research Center to investigate the material response of instrumented sphere-cone models made from PICA coated with NuSil ® . Analysis of the post-test data is presented and suggests that NuSil ® impacts the material response of test articles subjected to oxidizing atmospheres, as evidenced by a reduction of the measured surface temperature and in-depth temperature response. Finally, a four-step ablation mechanism is proposed based on observations from high-speed video, optical emission spectroscopy, thermal analysis, and post-test surface characterization.
The fluid-dynamic characteristics of flat plates, 5 deg and 10 deg wedges, and 5 deg and 10 deg cones have been investigated at Mach numbers from 16.3 to 23.9 in helium flow. The flat-plate results are for a leading-edge Reynolds number range of 584 to 19,500 and show that the induced pressure distribution is essentially linear with the hypersonic viscous interaction parameter bar X within the scope of this investigation. It is also shown that the rate at which the induced pressure varies with bar X is a linear function of the leading-edge Reynolds number. The wedge and cone results show that as the flow-deflection angle increases, the induced-pressure effects decrease and the measured pressures approach those predicted by inviscid shock theory.
Nonslender bodies of revolution with minimum pressure drag, noting derivation of minimizing curves for Newtonian flow
The Aeroassist Flight Experiment (AFE) vehicle is a 4.26 m-diameter raked and blunted elliptic cone that is intended to obtain aerodynamic and aerothermodynamic data on blunt bodies at velocities near 10 km/s and altitudes above 75 km. Previous wind tunnel tests of this configuration have been at Mach numbers of 10 or lower. The present test program consisted of five test conditions at Mach numbers between 15 and 20 and simulated altitudes from 75 km to 97 km. Static longitudinal aerodynamic coefficients, the lift-to-drag ratio, and the trim angle-of-attack were obtained for angles-of-attack between 7 deg and 27 deg. The effects of angle-of-attack and shock Reynolds number, Re2, have been studied. These results indicate the AFE configuration is longitudinally stable. The variation of normal force slope, CNalpha and lift-to-drag slope, (L/D)alpha as well as the trim angle-of-attack with Re2 are documented. The trim angle-of-attack increases with decreasing values of shock Reynolds number (increasing altitude) while the lift-to-drag ratio at the trimmed conditions is relatively constant.
Boundary-layer receptivity and stability of Mach 6 flows over smooth and rough seven-degree half-angle sharp-tipped cones are numerically investigated. The receptivity of the boundary layer to slow acoustic disturbances, fast acoustic disturbances, and vortical disturbances is considered. The effects of three-dimensional isolated roughness on the receptivity and stability are also simulated. The results for the smooth cone show that the instability waves are generated in the leading edge region and that the boundary layer is much more receptive to slow acoustic waves than to the fast acoustic waves. Vortical disturbances also generate unstable second modes, however the receptivity coefficients are smaller than that of the slow acoustic wave. Distributed roughness elements located near the nose region decreased the receptivity of the second mode generated by the slow acoustic wave by a small amount. Roughness elements distributed across the continuous spectrum increased the receptivity of the second mode generated by the slow and fast acoustic waves and the vorticity wave. The largest increase occurred for the vorticity wave. Roughness elements distributed across the synchronization point did not change the receptivity of the second modes generated by the acoustic waves. The receptivity of the second mode generated by the vorticity wave increased in this case, but the increase is lower than that occurred with the roughness elements located across the continuous spectrum. The simulations with an isolated roughness element showed that the second mode waves generated by the acoustic disturbances are not influenced by the small roughness element. Due to the interaction, a three-dimensional wave is generated. However, the amplitude is orders of magnitude smaller than the two-dimensional wave.