Search NASASearch

SEARCH · Search NASA

Results for “PERFECT GAS”

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.

At least 73 records · Page 4

Transition Effects on Heating in the Wake of a Blunt Body

A series of aerodynamic heating tests was conducted on a 70-deg sphere-cone planetary entry vehicle model in a Mach 10 perfect-gas wind tunnel at freestream Reynolds numbers based on diameter of 8.23x104 to 3.15x105. Surface heating distributions were determined from temperature time-histories measured on the model and on its support sting using thin-film resistance gages. The experimental heating data were compared to computations made using an axisymmetric/2D, laminar, perfect-gas Navier-Stokes solver. Agreement between computational and experimental heating distributions to within, or slightly greater than, the experimental uncertainty was obtained on the forebody and afterbody of the entry vehicle as well as on the sting upstream of the free-shear-layer reattachment point. However, the distributions began to diverge near the reattachment point, with the experimental heating becoming increasingly greater than the computed heating with distance downstream from the reattachment point. It was concluded that this divergence was due to transition of the wake free shear layer just upstream of the reattachment point on the sting.

Hollis, Brian R.

Combined LAURA-UPS hypersonic solution procedure

A combined solution procedure for hypersonic flowfields around blunted slender bodies was implemented using a thin-layer Navier-Stokes code (LAURA) in the nose region and a parabolized Navier-Stokes code (UPS) on the after body region. Perfect gas, equilibrium air, and non-equilibrium air solutions to sharp cones and a sharp wedge were obtained using UPS alone as a preliminary step. Surface heating rates are presented for two slender bodies with blunted noses, having used LAURA to provide a starting solution to UPS downstream of the sonic line. These are an 8 deg sphere-cone in Mach 5, perfect gas, laminar flow at 0 and 4 deg angles of attack and the Reentry F body at Mach 20, 80,000 ft equilibrium gas conditions for 0 and 0.14 deg angles of attack. The results indicate that this procedure is a timely and accurate method for obtaining aerothermodynamic predictions on slender hypersonic vehicles.

Wood, William A.

An approximate viscous shock layer technique for calculating chemically reacting hypersonic flows about blunt-nosed bodies

An approximate axisymmetric method was developed which can reliably calculate fully viscous hypersonic flows over blunt nosed bodies. By substituting Maslen's second order pressure expression for the normal momentum equation, a simplified form of the viscous shock layer (VSL) equations is obtained. This approach can solve both the subsonic and supersonic regions of the shock layer without a starting solution for the shock shape. The approach is applicable to perfect gas, equilibrium, and nonequilibrium flowfields. Since the method is fully viscous, the problems associated with a boundary layer solution with an inviscid layer solution are avoided. This procedure is significantly faster than the parabolized Navier-Stokes (PNS) or VSL solvers and would be useful in a preliminary design environment. Problems associated with a previously developed approximate VSL technique are addressed before extending the method to nonequilibrium calculations. Perfect gas (laminar and turbulent), equilibrium, and nonequilibrium solutions were generated for airflows over several analytic body shapes. Surface heat transfer, skin friction, and pressure predictions are comparable to VSL results. In addition, computed heating rates are in good agreement with experimental data. The present technique generates its own shock shape as part of its solution, and therefore could be used to provide more accurate initial shock shapes for higher order procedures which require starting solutions.

Cheatwood, F. Mcneil

High performance Solid Rocket Motor (SRM) submerged nozzle/combustion cavity flowfield assessment

Two and three dimensional internal flowfield solutions for critical points in the Space Shuttle solid rocket booster burn time were developed using the Lockheed Huntsville GIM/PAID Navier-Stokes solvers. These perfect gas, viscous solutions for the high performance motor characterize the flow in the aft segment and nozzle of the booster. Two dimensional axisymmetric solutions were developed at t = 20 and t = 85 sec motor burn times. The t = 85 sec solution indicates that the aft segment forward inhibitor stub produces vortices with are shed and convected downwards. A three dimensional 3.5 deg gimbaled nozzle flowfield solution was developed for the aft segment and nozzle at t = 9 sec motor burn time. This perfect gas, viscous analysis, provided a steady state solution for the core region and the flow through the nozzle, but indicated that unsteady flow exists in the region under the nozzle nose and near the flexible boot and nozzle/case joint. The flow in the nozzle/case joint region is characterized by low magnitude pressure waves which travel in the circumferential direction. From the two and three dimensional flowfield calculations presented it can be concluded that there is no evidence from these results that steady state gas dynamics is the primary mechanism resulting in the nozzle pocketing erosion experienced on SRM nozzles 8A or 17B. The steady state flowfield results indicate pocketing erosion is not directly initiated by a steady state gas dynamics phenomenon.

Freeman, J. A.

PARC Navier-Stokes code upgrade and validation for high speed aeroheating predictions

Applications of the PARC full Navier-Stokes code for hypersonic flowfield and aeroheating predictions around blunt bodies such as the Aeroassist Flight Experiment (AFE) and Aeroassisted Orbital Transfer Vehicle (AOTV) are evaluated. Two-dimensional/axisymmetric and three-dimensional perfect gas versions of the code were upgraded and tested against benchmark wind tunnel cases of hemisphere-cylinder, three-dimensional AFE forebody, and axisymmetric AFE and AOTV aerobrake/wake flowfields. PARC calculations are in good agreement with experimental data and results of similar computer codes. Difficulties encountered in flowfield and heat transfer predictions due to effects of grid density, boundary conditions such as singular stagnation line axis and artificial dissipation terms are presented together with subsequent improvements made to the code. The experience gained with the perfect gas code is being currently utilized in applications of an equilibrium air real gas PARC version developed at REMTECH.

Liver, Peter A.

Blunt body near wake flow field at Mach 6

Tests were conducted in a Mach 6 flow to examine the reattachment process of an axisymmetric free shear layer associated with the near wake of a 70 deg. half angle, spherically blunted cone with a cylindrical after body. Model angle of incidence was fixed at 0 deg. and free-stream Reynolds numbers based on body diameter ranged from 0.5 x 10(exp 6) to 4 x 10(exp 6). The sensitivity of wake shear layer transition on reattachment heating was investigated. The present perfect gas study was designed to compliment results obtained previously in facilities capable of producing real gas effects. The instrumented blunted cone model was designed primarily for testing in high enthalpy hypervelocity shock tunnels in both this country and abroad but was amenable for testing in conventional hypersonic blowdown wind tunnels as well. Surface heating rates were inferred from temperature - time histories from coaxial surface thermocouples on the model forebody and thin film resistance gages along the model base and cylindrical after body. General flow feature (bow shock, wake shear layer, and recompression shock) locations were visually identified by schlieren photography. Mean shear layer position and growth were determined from intrusive pitot pressure surveys. In addition, wake surveys with a constant temperature hot-wire anemometer were utilized to qualitatively characterize the state of the shear layer prior to reattachment. Experimental results were compared to laminar perfect gas predictions provided by a 3-D Navier Stokes code (NSHYP). Shear layer impingement on the instrumented cylindrical after body resulted in a localized heating maximum that was 21 to 29 percent of the forebody stagnation point heating. Peak heating resulting from the reattaching shear layer was found to be a factor of 2 higher than laminar predictions, which suggested a transitional shear layer. Schlieren flow visualization and fluctuating voltage time histories and spectra from the hot wire surveys across the shear layer substantiate this observation. The sensitivity of surface heating to forebody roughness was characterized for a reattaching shear layer. For example, at R(sub infinity), d = 4 x 10(exp 6), when the shear layer was transitional, the magnitude of peak heating from shear layer impingement was reduced by approximately 24 percent when transition grit was applied to the forebody. The spatial location of the local peak, however, remained unchanged.

Horvath, Thomas J.

Calculation of real-gas effects on blunt-body trim angles

The effect of vibrational excitation and dissociation at high temperatures on the trim angle of attack of a blunt lifting body is calculated for a nonequilibrium flow regime in air using a CFD technique. The vibrational-electronic temperature and the species densities are calculated assuming the flow to be in a nonequilibrium state. The forebody flow of a two-dimensional blunt body of the shape of the Apollo Command Module at a finite angle of attack is calculated. The results show that the pitching moment around a reference point is larger and the trim angle of attack is smaller for a reacting gas than for a perfect gas. The calculated shift in the trim angle due to the real-gas effect is of the same order as that seen during the Apollo flights.

Park, Chul

Hypersonic blunt body computations including real gas effects

Various second-order explicit and implicit TVD shock-capturing methods, a generalization of Roe's approximate Riemann solver, and a generalized flux-vector splitting scheme are used to study two-dimensional hypersonic real-gas flows. Special attention is given to the identification of some of the elements and parameters which can affect the convergence rate for high Mach numbers or real gases, but have negligible effect for low Mach numbers, for cases involving steady-state inviscid blunt flows. Blunt body calculations at Mach numbers of greater than 15 are performed to treat real-gas effects, and impinging shock results are obtained to test the treatment of slip surfaces and complex structures. Even with the addition of improvements, the convergence rate of algorithms in the hypersonic flow regime is found to be generally slower for a real gas than for a perfect gas.

Montagne, J.-L.

Base Flow and Heat Transfer Characteristics of a Four-Nozzle Clustered Rocket Engine: Effect of Nozzle Pressure Ratio

The base pressure and heating characteristics of a four-nozzle clustered rocket configuration is studied numerically with the aid of OVERFLOW Navier-Stokes code. A pressure ratio (chamber pressure to freestream static pressure) range of 990 to 5,920 and a freestream Mach number range of 2.5 to 3.5 are studied. The qualitative trends of decreasing base pressure with increasing pressure ratio and increasing base heat flux with increasing pressure ratio are correctly predicted. However, the predictions for base pressure and base heat flux show deviations from the wind tunnel data. The differences in absolute values between the computation and the data are attributed to factors such as perfect gas (thermally and calorically perfect) assumption, turbulence model inaccuracies in the simulation, and lack of grid adaptation.

Nallasamy, R.

Outflow Boundary Conditions for Turbine-Integrated Rotating Detonation Combustors

This study examines outflow boundary conditions (BCs) in computational fluid dynamics (CFD) simulations of a transition duct with and without guide vanes that converts supersonic flow exiting a rotating detonation combustor (RDC) to subsonic flow to drive a turbine. Since the flow exiting the transition duct has swirling shock waves with significant spatial and temporal variations in pressure, temperature, and Mach number, imposing proper BCs poses a challenge. To ensure all swirling shock waves exit the transition duct without creating non-physical reflected waves at its outlet, this study examined three outflow BCs: (1) the average pressure imposed at the duct’s outlet, (2) a nonreflecting BC (NRBC) with a specified average pressure imposed at the duct’s outlet, (3) the average pressure imposed at the outlet of an extension duct made up of a buffer layer and a sponge layer. This study is based on the three-dimensional, unsteady density-weighted-ensemble-averaged continuity, Navier–Stokes, and energy equations for a thermally perfect gas closed by the realizable k–ε model and “enhanced” wall functions. The results obtained show that imposing an average pressure at the transition duct’s outlet produces spurious waves that degrade the physical meaningfulness of the solution. When the NRBC was applied, swirling shock waves exited the duct’s outlet without creating spurious waves. However, its usage requires the gas to be thermally, as well as calorically, perfect, which this study shows could be a concern. By imposing the average pressure at the outlet of an extension duct, the gas does not need to be calorically perfect. The results obtained show the effects of the sponge layer’s length and coarsening ratio on damping nonuniformities in non-physical reflected waves to ensure the flow exiting the transition duct’s outlet can do so as if there are no boundaries present and has the desired average pressure—even though the BC is applied at the extension duct’s outlet.

gas turbines

Calculation of stiffness and damping coefficients for elastically supported gas foil bearings

The stiffness and damping coefficients of an elastically supported gas foil bearing are calculated. A perfect gas is used as the lubricant, and its behavior is described by the Reynolds equation. The structural model consists only of an elastic foundation. The fluid equations and the structural equations are coupled. A perturbation method is used to obtain the linearized dynamic coefficient equations. A finite difference formulation has been developed to solve for the four stiffness and the four damping coefficients. The effect of the bearing compliance on the dynamic coefficients is discussed in this paper.

Peng, J.-P.

Part 1 of a Computational Study of a Drop-Laden Mixing Layer

This first of three reports on a computational study of a drop-laden temporal mixing layer presents the results of direct numerical simulations (DNS) of well-resolved flow fields and the derivation of the large-eddy simulation (LES) equations that would govern the larger scales of a turbulent flow field. The mixing layer consisted of two counterflowing gas streams, one of which was initially laden with evaporating liquid drops. The gas phase was composed of two perfect gas species, the carrier gas and the vapor emanating from the drops, and was computed in an Eulerian reference frame, whereas each drop was tracked individually in a Lagrangian manner. The flow perturbations that were initially imposed on the layer caused mixing and eventual transition to turbulence. The DNS database obtained included transitional states for layers with various liquid mass loadings. For the DNS, the gas-phase equations were the compressible Navier-Stokes equations for conservation of momentum and additional conservation equations for total energy and species mass. These equations included source terms representing the effect of the drops on the mass, momentum, and energy of the gas phase. From the DNS equations, the expression for the irreversible entropy production (dissipation) was derived and used to determine the dissipation due to the source terms. The LES equations were derived by spatially filtering the DNS set and the magnitudes of the terms were computed at transitional states, leading to a hierarchy of terms to guide simplification of the LES equations. It was concluded that effort should be devoted to the accurate modeling of both the subgridscale fluxes and the filtered source terms, which were the dominant unclosed terms appearing in the LES equations.

Okong'o, Nora A.

Parallelization of Lower-Upper Symmetric Gauss-Seidel Method for Chemically Reacting Flow

Development of technologies for exploration of the solar system has revived an interest in computational simulation of chemically reacting flows since planetary probe vehicles exhibit non-equilibrium phenomena during the atmospheric entry of a planet or a moon as well as the reentry to the Earth. Stability in combustion is essential for new propulsion systems. Numerical solution of real-gas flows often increases computational work by an order-of-magnitude compared to perfect gas flow partly because of the increased complexity of equations to solve. Recently, as part of Project Columbia, NASA has integrated a cluster of interconnected SGI Altix systems to provide a ten-fold increase in current supercomputing capacity that includes an SGI Origin system. Both the new and existing machines are based on cache coherent non-uniform memory access architecture. Lower-Upper Symmetric Gauss-Seidel (LU-SGS) relaxation method has been implemented into both perfect and real gas flow codes including Real-Gas Aerodynamic Simulator (RGAS). However, the vectorized RGAS code runs inefficiently on cache-based shared-memory machines such as SGI system. Parallelization of a Gauss-Seidel method is nontrivial due to its sequential nature. The LU-SGS method has been vectorized on an oblique plane in INS3D-LU code that has been one of the base codes for NAS Parallel benchmarks. The oblique plane has been called a hyperplane by computer scientists. It is straightforward to parallelize a Gauss-Seidel method by partitioning the hyperplanes once they are formed. Another way of parallelization is to schedule processors like a pipeline using software. Both hyperplane and pipeline methods have been implemented using openMP directives. The present paper reports the performance of the parallelized RGAS code on SGI Origin and Altix systems.

Yoon, Seokkwan

Turbulent viscous-shock-layer solutions with strong vorticity interaction

Numerical solutions of the viscous-shock-layer equations governing laminar and turbulent flows of a perfect gas and radiating and nonradiating mixtures of perfect gases in chemical equilibrium are presented for hypersonic flow over spherically blunted cones and hyperboloids. The results are compared with boundary-layer and inviscid flow-field solutions. The agreement with the inviscid flow field data is satisfactory. The agreement with boundary-layer solutions is good except in regions of strong vorticity interaction. In these flow regions, the viscous-shock-layer solutions appear to be more satisfactory than the boundary-layer solutions.

Anderson, E. C.

Vorticity interaction effects on blunt bodies

Numerical solutions of the viscous shock layer equations governing laminar and turbulent flows of a perfect gas and radiating and nonradiating mixtures of perfect gases in chemical equilibrium are presented for hypersonic flow over spherically blunted cones and hyperboloids. Turbulent properties are described in terms of the classical mixing length. Results are compared with boundary layer and inviscid flowfield solutions; agreement with inviscid flowfield data is satisfactory. Agreement with boundary layer solutions is good except in regions of strong vorticity interaction; in these flow regions, the viscous shock layer solutions appear to be more satisfactory than the boundary layer solutions. Boundary conditions suitable for hypersonic viscous shock layers are devised for an advanced turbulence theory.

Anderson, E. C.