Search NASA⌕ Search

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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

An experimental parametric study of geometric, Reynolds number, and ratio of specific heats effects in three-dimensional sidewall compression scramjet inlets at Mach 6

The effects of the decrease in the ratio of specific heats have been investigated in generic 3D sidewall compression scramjet inlets with leading-edge sweep angles of 30 and 70 degrees. The effects of a decrease in ratio of specific heats were seen by comparing data from two facilities in two test gases: in the Langley Mach 6 CF4 Tunnel in tetrafluoromethane and in the Langley 15-Inch Mach 6 Air Tunnel in perfect gas air. In addition to the simulated real gas effects, the parametric effects of cowl position, contraction ratio, leading-edge sweep, and Reynolds number were investigated in the 15-Inch Mach 6 Air Tunnel. The models were instrumented with a total of 45 static pressure orifices distributed on the sidewalls and baseplate. Surface streamline patterns were examined via oil flow, and schlieren videos were made of the external flow field. The results of these tests have significant implications to ground based testing of inlets in facilities which do not operate at flight enthalpies.

Holland, Scott D.↗
NASA NTRS1993

Higher-order viscous shock-layer solutions for high altitude flows

A higher-order viscous shock-layer (HVSL) method has been developed and is used to obtain physically consistent results under varying degrees of low-density conditions for perfect-gas and nonequilibrium flows past long slender bodies. This method is implemented with higher-order body and shock slip conditions. Detailed comparisons with the direct simulation Monte Carlo (DSMC) method and Navier-Stokes (NS) calculations establish HVSL as a computationally efficient code for the low-density calculations. Further, it is shown that the deficiencies in the standard viscous shock-layer (VSL) method when predicting low-density flows lie in the governing equations and not with the slip boundary conditions as considered by earlier researchers.

Gupta, Roop N.↗
NASA NTRS1993

Viscous hypersonic shock-shock interaction on a blunt body at high altitude

The shock interaction produced when an incident oblique shock impinges on a blunt body, such as an engine inlet cowl lip of a hypersonic vehicle, is investigated for high altitude flight conditions. A perfect gas, Navier-Stokes numerical simulation of this problem at various altitudes representing continuum through transitional conditions is performed using the modified flux vector splitting method of Steger and Warming (1979). Two series of solutions are produced. First, a number-of-shock-positions are studied at a particular altitude and Mach number. Second, given a fixed shock position and Mach number, the interaction is investigated at several altitudes ranging from continuum to transitional flow conditions. It is shown that the interaction becomes fundamentally different as the density is lowered, and its effect on the overheating problem is progressively diminished. The maximum stagnation point heating at the highest altitude is reached only when the incident shock misses the cowl lip completely, and any interaction with the cowl bow shock that does occur takes place downstream and thus has little effect on the conditions at the stagnation point.

Comeaux, Keith A.↗
NASA NTRS1993

A viscous shock-layer analysis of 2-D and axisymmetric flows

Results are obtained for cylindrical leading edges of proposed transatmospheric vehicles by employing a 2D viscous shock-layer code for nonequilibrium and perfect gas flows. The accuracy and efficiency of the planar code is verified through detailed comparisons with other predictions. It is found to be as accurate and robust as its axisymmetric counterpart. This study includes results for nose radii ranging from 0.01 to 2.0 ft and half-angles of 5 and 6 deg for both cylindrically-blunted wedges and spherically blunted cones (included for comparison). Some results are presented as a ratio of the noncatalytic to the corresponding fully catalytic heating value to illustrate the maximum potential for a heating reduction in dissociated nonequilibrium flows. This ratio and the individual heating rates are smaller for cylindrically-blunted wedges with small nose radii, relative to the spherically-blunted cones (at the same radius). Therefore, a larger potential exists for heating reduction in cylindrically-blunted as compared with the spherically-blunted surfaces for finite-rate chemistry.

Gupta, Roop N.↗
NASA NTRS1993

Two-layer convective heating prediction procedures and sensitivities for blunt body reentry vehicles

This paper provides a description of procedures typically used to predict convective heating rates to hypersonic reentry vehicles using the two-layer method. These procedures were used to compute the pitch-plane heating distributions to the Apollo geometry for a wind tunnel test case and for three flight cases. Both simple engineering methods and coupled inviscid/boundary layer solutions were used to predict the heating rates. The sensitivity of the heating results in the choice of metrics, pressure distributions, boundary layer edge conditions, and wall catalycity used in the heating analysis were evaluated. Streamline metrics, pressure distributions, and boundary layer edge properties were defined from perfect gas (wind tunnel case) and chemical equilibrium and nonequilibrium (flight cases) inviscid flow-field solutions. The results of this study indicated that the use of CFD-derived metrics and pressures provided better predictions of heating when compared to wind tunnel test data. The study also showed that modeling entropy layer swallowing and ionization had little effect on the heating predictions.

Bouslog, Stanley A.↗
NASA NTRS1993

Application of an engineering inviscid-boundary layer method to slender three-dimensional vehicle forebodies

An engineering inviscid-boundary layer method has been modified for application to slender three-dimensional (3-D) forebodies which are characteristic of transatmospheric vehicles. An improved shock description in the nose region has been added to the inviscid technique which allows the calculation of a wider range of body geometries. The modified engineering method is applied to the perfect gas solution over a slender 3-D configuration at angle of attack. The method predicts surface pressures and laminar heating rates on the windward side of the vehicle that compare favorably with numerical solutions of the thin-layer Navier-Stokes equations. These improvements extend the 3-D capabilities of the engineering method and significantly increase its design applications.

Riley, Christopher J.↗
NASA NTRS1992

Requirements for facilities and measurement techniques to support CFD development for hypersonic aircraft

The design of a hypersonic aircraft poses unique challenges to the engineering community. Problems with duplicating flight conditions in ground based facilities have made performance predictions risky. Computational fluid dynamics (CFD) has been proposed as an additional means of providing design data. At the present time, CFD codes are being validated based on sparse experimental data and then used to predict performance at flight conditions with generally unknown levels of uncertainty. This paper will discuss the facility and measurement techniques that are required to support CFD development for the design of hypersonic aircraft. Illustrations are given of recent success in combining experimental and direct numerical simulation in CFD model development and validation for hypersonic perfect gas flows.

Sellers, William L., III↗
NASA NTRS1944

The numerical solution of compressible fluid flow problems

Numerical methods have been developed for obtaining the steady, adiabatic flow field of a frictionless, perfect gas about arbitrary two-dimensional bodies. The solutions include the subsonic velocity regions, the supersonic velocity regions, and the transition compression shocks, if required. Furthermore, the rotational motion and entropy changes following shocks are taken into account. Extensive use is made of the relaxation method. In this report the details of the methods of solution are emphasized so as to permit others to solve similar problems. Solutions already obtained are mentioned only by way of illustrating the possibilities of the methods described. The methods can be applied directly to wind tunnel and free air tests of arbitrary airfoil shapes at subsonic, sonic, and supersonic speeds.

Emmons, Howard W↗
NASA NTRS1993

Computational techniques for high-speed flows with viscous and chemical effects

Algorithms for solving the Euler and the Navier-Stokes equations in conjunction with chemical kinetic equations are presented. The convective flux is estimated from a quasi one dimensional interpolation procedure. Shock, contact, and expansion waves and thermochemical nonequilibrium phenomena are captured by the Lax-Friedrichs technique. Relaxation techniques were developed to enhance their effectiveness in dealing with spatial and temporal stiffness associated with the physical problems. Both explicit and implicit smoothers were implemented into the standard multigrid time stepping method. Unsteady and steady scalar problems are discussed. A perfect gas and equilibrium air shock tube problem is investigated. Numerical schemes and techniques are compared for the problems of shock and boundary layer interaction and three dimensional viscous, nonequilibrium flow encompassing an aerobrake. The results are comparable in accuracy against other high order non-oscillatory techniques. The multigrid methods are assessed using a Mach 8 flow over a complete planar body, a sphere, and a blunt delta wing at 20 deg incidence. Applying an implicit multigrid method on a nested grid of 128 by 64 nodes the reduction factor is 0.25. The central processing unit reduction factor is 2.2 after both the single and multigrid Runge-Kutta solutions converged to machine zero on a grid of 37 by 41 by 73 nodes.

Li, C. P.↗
NASA NTRS1994

Combined LAURA-UPS solution procedure for chemically-reacting flows

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 flowfields. 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 noncatalytic 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 noncatalytic 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 and the nonequilibrium results are compared with a perfect gas solution, showing that while the surface pressure is relatively unchanged by the inclusion of reacting chemistry the nonequilibrium heating is 25 percent higher. The procedure demonstrates significant, order of magnitude reductions in solution time and required memory for the three dimensional case over an all thin-layer Navier-Stokes solution.

Wood, William A.↗
NASA NTRS1971

Shock wave structure using nonlinear model Boltzmann equations

The structure of a strong plane shock wave in a monatomic rarefied perfect gas is one of the simplest problems able to be posed in kinetic theory, and one of the hardest to solve. Its simplicity lies in the absence of solid boundaries, geometrical complications, or internal molecular energy. Its difficulty arises from the great departure of the gas from equilibrium within the shock, which invalidates many of the techniques used successfully elsewhere in kinetic theory. In addition to this theoretical challenge, the modern development of ballistics and hypersonic flight has helped to stimulate extensive theoretical and experimental interest in the shock problem. The experimenters in turn have encountered great difficulties on account of the very small physical dimensions of shocks. In fact, until very recently indeed, any close comparisons of theoretical and experimental shock structure results have been rather unprofitable due to the inadequacies of both theory and experiment. During the last few years this situation has been appreciably improved by development of the Monte Carlo method. This allows idealized 'experiments' to be performed on large computers instead of in wind tunnels, using a known intermolecular force law. The most developed of these methods has been shown to be equivalent theoretically to the Boltzmann equation and to give results which agree extremely closely with measurements of high accuracy. Thus Monte Carlo results not only form the soundest basis for our present theoretical knowledge of shock wave structure, but, for purposes of developing other theories, can also be considered a very valuable experimental resource. However, such results remain very expensive to obtain. In this thesis we develop more economical kinetic theory methods for the approximate prediction of shock structure, and compare our results with those of the Monte Carlo method.

Segal, Ben Maurice↗
NASA NTRS1994

New Developments in the Method of Space-Time Conservation Element and Solution Element-Applications to Two-Dimensional Time-Marching Problems

A new numerical discretization method for solving conservation laws is being developed. This new approach differs substantially in both concept and methodology from the well-established methods, i.e., finite difference, finite volume, finite element, and spectral methods. It is motivated by several important physical/numerical considerations and designed to avoid several key limitations of the above traditional methods. As a result of the above considerations, a set of key principles for the design of numerical schemes was put forth in a previous report. These principles were used to construct several numerical schemes that model a 1-D time-dependent convection-diffusion equation. These schemes were then extended to solve the time-dependent Euler and Navier-Stokes equations of a perfect gas. It was shown that the above schemes compared favorably with the traditional schemes in simplicity, generality, and accuracy. In this report, the 2-D versions of the above schemes, except the Navier-Stokes solver, are constructed using the same set of design principles. Their constructions are simplified greatly by the use of a nontraditional space-time mesh. Its use results in the simplest stencil possible, i.e., a tetrahedron in a 3-D space-time with a vertex at the upper time level and other three at the lower time level. Because of the similarity in their design, each of the present 2-D solvers virtually shares with its 1-D counterpart the same fundamental characteristics. Moreover, it is shown that the present Euler solver is capable of generating highly accurate solutions for a famous 2-D shock reflection problem. Specifically, both the incident and the reflected shocks can be resolved by a single data point without the presence of numerical oscillations near the discontinuity.

Chang, Sin-Chung↗
NASA NTRS1993

Time-Accurate Unsteady Flow Simulations Supporting the SRM T+68-Second Pressure Spike Anomaly Investigation (STS-54B)

Time-accurate unsteady flow simulations are being performed supporting the SRM T+68sec pressure 'spike' anomaly investigation. The anomaly occurred in the RH SRM during the STS-54 flight (STS-54B) but not in the LH SRM (STS-54A) causing a momentary thrust mismatch approaching the allowable limit at that time into the flight. Full-motor internal flow simulations using the USA-2D axisymmetric code are in progress for the nominal propellant burn-back geometry and flow conditions at T+68-sec--Pc = 630 psi, gamma = 1.1381, T(sub c) = 6200 R, perfect gas without aluminum particulate. In a cooperative effort with other investigation team members, CFD-derived pressure loading on the NBR and castable inhibitors was used iteratively to obtain nominal deformed geometry of each inhibitor, and the deformed (bent back) inhibitor geometry was entered into this model. Deformed geometry was computed using structural finite-element models. A solution for the unsteady flow has been obtained for the nominal flow conditions (existing prior to the occurrence of the anomaly) showing sustained standing pressure oscillations at nominally 14.5 Hz in the motor IL acoustic mode that flight and static test data confirm to be normally present at this time. Average mass flow discharged from the nozzle was confirmed to be the nominal expected (9550 lbm/sec). The local inlet boundary condition is being perturbed at the location of the presumed reconstructed anomaly as identified by interior ballistics performance specialist team members. A time variation in local mass flow is used to simulate sudden increase in burning area due to localized propellant grain cracks. The solution will proceed to develop a pressure rise (proportional to total mass flow rate change squared). The volume-filling time constant (equivalent to 0.5 Hz) comes into play in shaping the rise rate of the developing pressure 'spike' as it propagates at the speed of sound in both directions to the motor head end and nozzle. The objectives of the present analysis are to: (1) capture the dynamic responses of the motor combustion gas flow to correlate with available low-frequency (less than 12.5 sample/sec) data and (2) observe the high-frequency (up to 50 Hz) characteristics of the response to determine any potentials for dynamic coupling.

Dougherty, N. S.↗
NASA NTRS1995

Recent improvements to and validation of the one dimensional NASA wave rotor model

A numerical model has been developed at the NASA Lewis Research Center which can predict both the unsteady flow quantities within a wave rotor passage and the steady averaged flows in the ports. The model is based on the assumptions of one-dimensional, unsteady, perfect gas flow. The model assesses not only the dominant wave behavior, but the loss effects of finite passage opening time, leakage from the passage ends, viscosity, and heat transfer to and from the passages. The model operates in the rotor reference frame; however, until recently no account was made for the often significant effect of the rotor circumferential velocity component. The present model accounts for this by modifying the passage boundary conditions, allowing the internal computational scheme to remain the rotor reference frame, while quantities such as inlet duct stagnation properties may be specified in the fixed or absolute reference frame. Accurate modeling of this effect is critical to successful wave rotor analysis and design, particularly in off-design predictions where the flows in the inlet ducts are mismatched with the rotor passages and significant turning may take place (i.e., work is done on the gas). The relative simplicity of the model makes it useful for design and optimization, as well as analysis, of wave rotor cycles for many applications. This report, building on several earlier papers, describes the most recent modifications to the model. These include accounting for the relative/absolute transition at the passage boundaries and refinements to the viscous source term correlation which resulted from this accounting. Comparison of model predictions with measured data is then presented and discussed.

Paxson, Daniel E.↗
NASA NTRS1995

A numerical model for dynamic wave rotor analysis

A numerical model has been developed which can predict the dynamic (and steady state) performance of a wave rotor, given the geometry and time dependent boundary conditions. The one-dimensional, perfect gas, CFD based code tracks the gasdynamics in each of the wave rotor passages as they rotate past the various ducts. The model can operate both on and off-design, allowing dynamic behavior to be studied throughout the operating range of the wave rotor. The model accounts for several major loss mechanisms including finite passage opening time, fluid friction, heat transfer to and from the passage walls, and leakage to and from the passage ends. In addition, it can calculate the amount of work transferred to and from the fluid when the flow in the ducts is not aligned with the passages such as occurs in off-design operation. Since it is one-dimensional, the model runs reasonably fast on a typical workstation. This paper will describe the model and present the results of some transient calculations for a conceptual four port wave rotor designed as a topping cycle for a small gas turbine engine.

Paxson, D. E.↗
NASA NTRS1994

Penetration below a convective zone

Two-dimensional numerical simulations are used to investigate how fully compressible nonlinear convection penetrates into a stably stratified zone beneath a stellar convection zone. Estimates are obtained of the extent of penetration as the relative stability S of the stable to the unstable zone is varied over a broad range. The model deals with a perfect gas possessing a constant dynamic viscosity. The dynamics is dominated by downward-directed plumes which can extend far into the stable material and which can lead to the excitation of a broad spectrum of internal gravity waves in the lower stable zone. The convection is highly time dependent, with the close coupling between the lateral swaying of the plumes and the internal gravity waves they generate serving to modulate the strength of the convection. The depth of penetration delta, determined by the position where the time-averaged kinetic flux has its first zero in the stable layer, is controlled by a balance between the kinetic energy carried into the stable layer by the plumes and the buoyancy braking they experience there. A passive scalar is introduced into the unstable layer to evaluate the transport of chemical species downward. Such a tracer is effectively mixed within a few convective overturning times down to a depth of delta within the stable layer. Analytical estimates based on simple scaling laws are used to interpret the variation of delta with S, showing that it first involves an interval of adiabatic penetration if the local Peclet number of the convection exceeds unity, followed by a further thermal adjustment layer, the depths of each interval scaling in turn as S(exp -1) and S(exp -1/4). These estimates are in accord with the penetration results from the simulations.

Hurlburt, Neal E.↗
NASA NTRS1994

Data-parallel lower-upper relaxation method for reacting flows

The implicit lower-upper symmetric Gauss-Seidel (LU-SGS) method of Yoon and Jameson is modified for use on massively parallel computers. The method has been implemented on the Thinking Machines CM-5 and the MasPar MP-1 and MP-2, where large percentages of the theoretical peak floating point performance are obtained. It is shown that the new data-parallel LU relaxation method has better convergence properties than the original method for two different inviscid compressible flow simulations. The convergence is also improved for five-species reacting air computations. The performance of the method on various partitions of the CM-5 and on the MasPar computers is discussed. The new method shows promise for the efficient simulation of very large perfect gas and reacting flows.

Candler, Graham V.↗
NASA NTRS1995

Boundary layer blockage in expansion tube nozzles

The results of a first order perfect gas correction for the effects of the boundary layer formation within expansion tubes with nozzles are presented. The analytical model developed to describe the boundary layer formation within the expansion tube and an expansion nozzle located at the end of the acceleration tube is based on the Karman integral equations. The results of this analytical model are compared with experimental data from an expansion diffuser. The model provides a useful tool for the preliminary design of nozzles for such facilities.

Sudnitsin, Olga↗