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 127 records · Page 7

Extrapolation From Wind Tunnel to Flight: Shuttle Orbiter Aerodynamics

The paper reviews a combined numerical and experimental activity on the Shuttle Orbiter, first performed at NASA Langley within the Orbiter Experiment (OEX) and subsequently at ESA, as part of the AGARD FDP WG 18 activities. The study at Langley was undertaken to resolve the pitch up anomaly observed during the entry of the first flight of the Shuttle Orbiter. The present paper will focus on real gas effects on aerodynamics and not on heating. The facilities used at NASA Langley were the 15-in. Mach 6, the 20-in, Mach 6, the 31-in. Mach 10 and the 20-in. Mach 6 CF4 facility. The paper focuses on the high Mach, high altitude portion of the first entry of the Shuttle where the vehicle exhibited a nose-up pitching moment relative to pre-flight prediction of (AC,) = 0.03. In order to study the relative contribution of compressibility, viscous interaction and real gas effects on basic body pitching moment and flap efficiency, an experimental study was undertaken to examine the effects of Mach, Reynolds and ratio of specific heats at NASA. At high Mach, a decrease of gamma occurs in the shock layer due to high temperature effects. The primary effect of this lower specific heat ratio is a decrease of the pressure on the aft windward expansion surface of the Orbiter causing the nose-up pitching moment. Testing in the heavy gas, Mach 6 CF4 tunnel, gave a good simulation of high temperature effects. The facilities used at ESA were the 1m Mach 10 at ONERA Modane, the 0.7 m hot shot F4 at ONERA Le Fauga and the 0.88 m piston driven shock tube HEG at DLR Goettingen. Encouraging good force measurements were obtained in the F4 facility on the Orbiter configuration. Testing of the same model in the perfect gas Mach 10 S4 Modane facility was performed so as to have "reference" conditions. When one compares the P4 and S4 test results, the data suggests that the Orbiter "pitch up" is due to real gas effects. In addition, pressure measurements, performed on the aft portion of the windward side of the Halis configuration in HEG and P4, confirm that the pitch up is mainly attributed to a reduction of pressure due to a local decrease in gamma.

J. Muylaert↗

Numerical simulation of reentry flow around the Shuttle Orbiter including real gas effects

The flight environment of a reentry vehicle is predicted from the numerical solution of fluid dynamics equations for the region between the blunt nose and the swept wings. The inviscid portion of the shock layer is modeled by the Euler equations, but the laminar viscous flow adjacent to the wall is modeled by the approximate parabolic Navier-Stokes equations. The approximations made to the axial gradients of pressure and diffusive fluxes enable the coupled inviscid and viscous equations to be solved efficiently along the body axis. The equilibrium air aftbody code contains significant improvements over its predecessor which only considers a perfect gas model and noncircular configurations. The inclusion of numerical damping either explicitly or implicitly has extended its capabilities for predicting flow field around a winged configuration at higher Machs and greater angles of attack. The results are obtained on the cylindrical coordinates and satisfactory for the Shuttle Orbiter at a free-stream Mach number of 22 and an angle of attack of 40 deg. Also discussed are the inviscid formulation and its application for wind-tunnel conditions.

Li, C. P.↗

Compressible laminar boundary layer with real gas effects for flight conditions to M(e) = 8 and T(o) = 2500 K

One of the problems for the hypersonic flights is the lack of accurate methods for predicting the boundary layer parameters at high Mach numbers and temperatures. At high Mach numbers and temperatures, the real gas effects should be included in calculations of the compressible laminar and turbulent boundary layers. By introducing the Mangler-Levy-Lees transformation, the compressible laminar boundary layer equations were transformed into two ordinary differential equations. To solve these ordinary differential equations for the flow over a flat plate, the Runge-Kutta method of order 5 and 6 was used with a secant method of solving the simultaneous nonlinear equations for the boundary layer conditions. Instead of using the perfect gas assumptions, the properties of air in the temperature range of 100 to 2500 K at M(e) = 8 were used for the calculations. The calculations were carried out for free flight conditions. The real gas effects were significant for the conditions of high freestream Mach numbers, and the boundary layer parameters from the real gas solutions were quite different from the existing formulas and solutions.

Nagamatsu, H. T.↗

Simulation of Oxygen Disintegration and Mixing With Hydrogen or Helium at Supercritical Pressure

The simulation of high-pressure turbulent flows, where the pressure, p, is larger than the critical value, p(sub c), for the species under consideration, is relevant to a wide array of propulsion systems, e.g. gas turbine, diesel, and liquid rocket engines. Most turbulence models, however, have been developed for atmospheric-p turbulent flows. The difference between atmospheric-p and supercritical-p turbulence is that, in the former situation, the coupling between dynamics and thermodynamics is moderate to negligible, but for the latter it is very significant, and can dominate the flow characteristics. The reason for this stems from the mathematical form of the equation of state (EOS), which is the perfect-gas EOS in the former case, and the real-gas EOS in the latter case. For flows at supercritical pressure, p, the large eddy simulation (LES) equations consist of the differential conservation equations coupled with a real-gas EOS. The equations use transport properties that depend on the thermodynamic variables. Compared to previous LES models, the differential equations contain not only the subgrid scale (SGS) fluxes, but also new SGS terms, each denoted as a correction. These additional terms, typically assumed null for atmospheric pressure flows, stem from filtering the differential governing equations, and represent differences between a filtered term and the same term computed as a function of the filtered flow field. In particular, the energy equation contains a heat-flux correction (q-correction) that is the difference between the filtered divergence of the heat flux and the divergence of the heat flux computed as a function of the filtered flow field. In a previous study, there was only partial success in modeling the q-correction term, but in this innovation, success has been achieved by using a different modeling approach. This analysis, based on a temporal mixing layer Direct Numerical Simulation database, shows that the focus in modeling the q-correction should be on reconstructing the primitive variable gradients rather than their coefficients, and proposes the approximate deconvolution model (ADM) as an effective means of flow field reconstruction for LES heat flux calculation. Further, results for a study conducted for temporal mixing layers initially containing oxygen in the lower stream, and hydrogen or helium in the upper stream, show that, for any LES, including SGS-flux models (constant-coefficient Gradient or Scale-Similarity models, dynamic-coefficient Smagorinsky/Yoshizawa or mixed Smagorinsky/Yoshizawa/Gradient models), the inclusion of the q-correction in the LES leads to the theoretical maximum reduction of the SGS heat-flux difference. The remaining error in modeling this new subgrid term is thus irreducible.

Bellan, Josette↗

An Experimental Parametric Study of Geometric, Reynolds Number, and Ratio of Specific Heats Effects in Three-Dimensional Sidewall Compression Scramjet Inlets at Mach 6

Since mission profiles for airbreathing hypersonic vehicles such as the National Aero-Space Plane include single-stage-to-orbit requirements, real gas effects may become important with respect to engine performance. The effects of the decrease in the ratio of specific heats have been investigated in generic three-dimensional 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 (where gamma=1.22) and in the Langley 15-Inch Mach 6 Air Tunnel in perfect gas air (where gamma=1.4). 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.↗

Turbulent Aeroheating Testing of Mars Science Laboratory Entry Vehicle

An experimental investigation of turbulent aeroheating on the Mars Science Laboratory entry vehicle heat shield has been conducted in the Arnold Engineering Development Center Hypervelocity Wind Tunnel No. 9. Testing was performed on a 6-in. (0.1524 m) diameter MSL model in pure N2 gas in the tunnel's Mach 8 and Mach 10 nozzles at free stream Reynolds numbers of 4.1 x 10(exp 6)/ft to 49 x 10(exp 6)/ft (1.3 x 10(exp 7)/m to 19 x 10(exp 6/ft) and 1.2 x 10(exp 6)/ft to 19 x 10(exp 6)/ft (0.39 x 10(exp 7)/m to 62 x 10(exp 7)/m), respectively. These conditions were sufficient to span the regime of boundary-layer flow from completely laminar to fully-developed turbulent flow over the entire forebody. A supporting aeroheating test was also conducted in the Langley Research Center 20-Inch Mach 6 Air Tunnel at free stream Reynolds number of 1 x 10(exp 6)/ft to 7 x 10(exp 6)/ft (0.36 x 10(exp 7)/m to 2.2 x 10(exp 7)/m) in order to help corroborate the Tunnel 9 results. A complementary computational fluid dynamics study was conducted in parallel to the wind tunnel testing. Laminar and turbulent predictions were generated for the wind tunnel test conditions and comparisons were performed with the data for the purpose of helping to define uncertainty margins on predictions for aeroheating environments during entry into the Martian atmosphere. Data from both wind tunnel tests and comparisons with the predictions are presented herein. It was concluded from these comparisons that for perfect-gas conditions, the computational tools could predict fully-laminar or fully-turbulent heating conditions to within 12% or better of the experimental data.

Hollis, Brian R.↗

A scalable compressible volume of fluid solver using a stratified flow model

Here, in this paper, a Volume of Fluid method which uses a Stratified Flow model for flux calculation is proposed to solve a compressible multiphase flow in a way that has high parallel efficiency. The solver is part of an open-source computing tool to study fuel entrainment and combustion in the application of paraffin-based hybrid rocket motors. This paper focuses on the compressible Volume of Fluid solver for the liquid and gaseous phases in this environment. This solver uses perfect gas and stiffened gas models for equation of state, with the ability to easily add other models as desired. Additionally a comparison of a Riemann solver versus an AUSM+up scheme is conducted, showing improvement with the later in the solution of multiphase flows. This solver successfully reproduces shock tube cases with the expected accuracy, as well as simulating two-dimensional shear and gravity-driven flows. Results also show that the code is able to utilize a Stratified Flow model to evaluate shear flow without the need for interface reconstruction or gradient calculations. Finally, it is shown that the solver has near-ideal scaling under strong scaling tests as well as good performance in static scaling, giving improved performance over current options and implying future advancement of high-performance computing efficiency for multiphase flow solvers.

97 MATHEMATICS AND COMPUTING↗

A finite scale model for shock structure

In this paper we explore the use of a finite scale model of fluid dynamics to predict the finite structure of a shock wave in a perfect gas. We begin by documenting the history and issues that have arisen when Navier–Stokes theory is applied to the shock structure problem, and continue by motivating the improvement that finite scale theory might provide from its representation of inviscid (anomalous) dissipation. Furthermore, our primary results include the formulation of a traveling wave equation from finite scale theory, an analysis of the solutions of that equation as regards shock width and monotonicity properties, and an estimation of the critical parameter of the theory from experimental data.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Finite scale theory: Predicting nature’s shocks

In this paper, we describe a continuum model that accurately reproduces the experimentally measured structure of physical shocks in a perfect gas. In this work, we begin by presenting a history of shock structure research, theoretical, experimental and numerical, to quantify the significant discrepancies between Navier–Stokes predictions and laboratory measurements. In our first main result, we discuss modifications that generalize the Chapman–Enskog approximation and lead to a new continuum model termed the Finite Scale Equations. In our second main result, we use this continuum model to calculate shock structure with results comparable in accuracy to numerical simulations based on the Boltzmann equation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Direct numerical simulation of supercritical oxy-methane mixing layers with CO 2 substituted counterparts

Direct numerical simulations (DNS) of temporally developing, three dimensional, CH 4 /CO 2 , CH 4 /O 2 and CO 2 /O 2 mixing layers, are conducted at a supercritical pressure of 300 atm. To effectively model the supercritical regime, the employed formulation includes the compressible form of the governing equations, the cubic Peng-Robinson equation of state and a generalized formulation for heat and mass flux vectors derived from non-equilibrium thermodynamics and fluctuation theory. A linear inviscid stability analysis is also performed for each case, to determine its most unstable wavelength. Flow visualizations reveal the presence of high density gradient magnitude regions for all three mixing layers, with conditional averages indicating increased presence of heavier fluid species within these regions. No significant departures are observed from perfect gas behavior, with compressibility factors very close to unity for all three mixing cases. Applicability of presumed probability density function methods (PDF) is examined for the three supercritical mixing layers. An a priori analysis is also conducted to investigate various simplifying assumptions employed in modeling various subgrid scale (SGS) flux models. Two additional terms are identified in the large eddy simulations (LES) equations, the gradient of SGS contribution of pressure in the momentum equation and the gradient of SGS contribution of heat flux in energy equation, whose magnitudes are similar and comparable with their respective resolved terms. The performance of the scale similarity model to represent these additional terms is investigated. Lastly, the performance of Smagorinsky, gradient and scale similarity models are also investigated.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Efficient GPU Implementation of Automatic Differentiation for Computational Fluid Dynamics

Many scientific and engineering applications require repeated calculations of derivatives of output functions with respect to input parameters. Automatic Differentiation (AD) is a method that automates derivative calculations and can significantly speed up code development. In Computational Fluid Dynamics (CFD), derivatives of flux functions with respect to state variables (Jacobian) are needed for efficient solutions of the nonlinear governing equations. AD of flux functions on graphics processing units (GPUs) is challenging as flux computations involve many intermediate variables that create high register pressure and require significant memory traffic because of the need to store the derivatives. This paper presents a forward-mode AD method based on multivariate dual numbers that addresses these challenges and simultaneously reduces the floating-point operation count. The dimension of the multivariate dual numbers is optimized for performance. The flux computations are restructured to minimize the number of temporary variables and reduce register pressure. For effective utilization of memory bandwidth, shared memory is used to store the local flux Jacobian. This AD implementation is compared with several other Jacobian implementations on an NVIDIA V100 GPU (V100). For three-dimensional perfect-gas compressible-flow equations implemented in a practical CFD code, the AD implementation of a flux Jacobian based on multivariate dual numbers of dimension 5 outperforms all other GPU AD implementations on V100. Its performance is comparable with the optimized hand-differentiated version. Finally, the implementation achieves 75% of the peak floating-point throughput and 61 % of the peak global device memory bandwidth usage.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

A Polytropic Approximation of Compressible Flow in Pipes With Friction

This paper demonstrates the usefulness of treating subsonic Fanno flow (adiabatic flow, with friction, of a perfect gas in a constant-area pipe) as a polytropic process. It is shown that the polytropic model allows an explicit equation for mass flow rate to be developed. The concept of the energy transfer ratio is used to develop a close approximation to the polytropic index. Explicit equations for mass flow rate and net expansion factor in terms of upstream properties and pressure ratio are developed for Fanno and isothermal flows. An approximation for choked flow is also presented. Here, the deviation of the results of this polytropic approximation from the values obtained from a traditional gas dynamics analysis of subsonic Fanno flow is quantified and discussed, and a typical design engineering problem is analyzed using the new method.

42 ENGINEERING↗

Modern infinitesimals and the entropy jump across an inviscid shock wave

Here, this article applies nonstandard analysis to study the generalized solutions of entropy and energy across one-dimensional shock waves in a compressible, inviscid, perfect gas. Nonstandard analysis is an area of modern mathematics that studies number systems that contain both infinitely large and infinitely small numbers. For an inviscid shock wave, it is assumed that the shock thickness occurs on an infinitesimal interval and that the jump functions for the field variables are smoothly defined on this interval. A weak converse to the existence of the entropy peak is derived and discussed. Generalized solutions of the Euler equations for entropy and energy are then derived for both theoretical and realistic normalized velocity profiles.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Symmetry Determining Equations of the Rankine-Hugoniot Equations for Variable Velocity Shock Waves

The “constant velocity piston” problem (Fig. 1), also known as the “piston problem,” is a standard model for a one dimensional, in our case linear, symmetric shock wave moving through an inviscid, perfect gas. The model can be divided into two regions - a perturbed section on the left and an unperturbed section on the right - by a moving shock wave moving left to right. Both the perturbed and unperturbed sections, i.e. the shocked and unshocked regions, respectively, obey the Eulerian conservation equations; however, at the exact location of the shock, there is a mathematical discontinuity not satisfied by the Euler equations. To ensure continuity and conservation of certain quantities when crossing between the unshocked and shocked regions, we evoke a series of equations derived from the Eulerian conservation equations, called the Rankine-Hugoniot equations, or “jump” equations as it is often referred to in the literature on the topic. The classical constant-velocity piston problem assumes the piston features a constant driving velocity (among many other willing suspensions of belief required in the pursuit of a first principles equation model); consequent to this assumption is a constant-velocity shock and a constant-velocity shocked flow state. However, using Lie Group Theory (LGT), also known as symmetry analysis, we can attempt to reinterpret the model with a shock wave of variable velocity in time and space. An extension of the model in this way opens up the possibility for obtaining new analytical solutions to the piston problem for certain shock velocity models. In this report, we use LGT to derive the symmetry determining equations (SDEs), whose solutions are Lie groups, which permit analytical solutions. In the future, we can then use the SDEs to define constraint equations on the shock velocity model and what the successive solutions to the Euler equations might be based off such constraints. This report is structured as follows: Section 2 provides a brief derivation of the Rankine-Hugoniot (“jump”) equations; Section 3 gives an overview of Lie group theory; Section 4 derives the SDEs of the jump equations; Section 5 derives the Euler conservation equations for fluids; and Section 6 presents concluding remarks and opportunities for future studies.

42 ENGINEERING↗

Galactic winds

Galactic wind as steady radial perfect gas flow from stars and central gravitating source

Axford, W. I.↗