Search NASA⌕ Search

SEARCH · Search NASA

Results for “Non-Equilibrium Flow”

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 109 records · Page 6

Upwind differencing and LU factorization for chemical non-equilibrium Navier-Stokes equations

By means of either the Roe or the Van Leer flux-splittings for inviscid terms, in conjunction with central differencing for viscous terms in the explicit operator and the Steger-Warming splitting and lower-upper approximate factorization for the implicit operator, the present, robust upwind method for solving the chemical nonequilibrium Navier-Stokes equations yields formulas for finite-volume discretization in general coordinates. Numerical tests in the illustrative cases of a hypersonic blunt body, a ramped duct, divergent nozzle flows, and shock wave/boundary layer interactions, establish the method's efficiency.

Shuen, Jian-Shun↗

Analysis of Flow From Arc-Jet Spectra

Materials testing is carried out at the JSC Atmospheric Reentry Materials and Structures Facility. A flow diagnostics program is under development to characterize the energy distribution in arc-jet flows used to simulate atmospheric entry. Heat transfer to materials depends on flow properties which includes the composition of and energy distribution among the atoms, ions, molecules, and molecular ions which make up the flow. This project involves analysis of shock layer characteristics from the radiation emitted and experimentally collected from near the front of the shock to near the surface of the material. The analysis has yielded information on relative populations of neutral molecules and molecular ions within the layer. In determining non-equilibrium temperatures within the layer, some insight into the spectral constants used to compute radiative emission has been gained.

Blackwell, H. E.↗

Hypervelocity Heat-Transfer Measurements in an Expansion Tube

A series of experiments has been conducted in the NASA HYPULSE Expansion Tube, in both CO2 and air test gases, in order to obtain data for comparison with computational results and to assess the capability for performing hypervelocity heat-transfer studies in this facility. Heat-transfer measurements were made in both test gases on 70 deg sphere-cone models and on hemisphere models of various radii. HYPULSE freestream flow conditions in these test gases were found to be repeatable to within 3-10%, and aerothermodynamic test times of 150 microsec in CO2 and 125 microsec in air were identified. Heat-transfer measurement uncertainty was estimated to be 10-15%. Comparisons were made with computational results from the non-equilibrium Navier-Stokes solver NEQ2D. Measured and computed heat-transfer rates agreed to within 10% on the hemispheres and on the sphere-cone forebodies, and to within 10% in CO2 and 25% in air on the afterbodies and stings of the sphere-cone models.

Hollis, Brian R.↗

Numerical simulations of non-equilibrium shock layers with efficient implicit schemes

Current and future calculations of nonequilibrium shock layers require the use of a very large number of equations, due to a multiplicity of chemical species, excited states, and internal energy modes. The computational cost associated with the use of standard implicit methods becomes prohibitive; it is therefore desirable to examine the potential of several methods and determine if any can be projected to be more efficient and accurate for large systems of equations. Here, the performance of several implicit schemes on several simple practical examples of reacting flows is examined. The Euler equations are solved by three different implicit methods, and two methods of coupling between the fluid dynamics and the chemistry are studied. Several cases of stiffness are considered and both 1D and 2D examples are computed.

Cambier, Jean-Luc↗

Fuel/oxidizer-rich high-pressure preburners

The analyses, designs, fabrication, and cold-flow acceptance testing of LOX/RP-1 preburner components required for a high-pressure staged-combustion rocket engine are discussed. Separate designs of injectors, combustion chambers, turbine simulators, and hot-gas mixing devices are provided for fuel-rich and oxidizer-rich operation. The fuel-rich design addresses the problem of non-equilibrium LOX/RP-1 combustion. The development and use of a pseudo-kinetic combustion model for predicting operating efficiency, physical properties of the combustion products, and the potential for generating solid carbon is presented. The oxygen-rich design addresses the design criteria for the prevention of metal ignition. This is accomplished by the selection of materials and the generation of well-mixed gases. The combining of unique propellant injector element designs with secondary mixing devices is predicted to be the best approach.

Schoenman, L.↗

Heating analysis for a Lunar Transfer Vehicle at near-equilibrium flow conditions

A heating analysis for a 15.2 m diameter Lunar Transfer Vehicle (LTV) at 0 and 10.6 deg angle of attack for a nominal trajectory through the earth's atmosphere is described. The analysis utilizes the Langley Aerothermodynamic Upwind Relaxation Algorithm (LAURA) with thin-layer, Navier-Stokes, thermochemical nonequilibrium options. Radiative heating levels are calculated using the Langley Optimized RAdiative Nonequilibrium (LORAN) and the Non-EQuilibrium AIr Radiation (NEQAIR) codes. At peak heating, the shock layer is substantially in equilibrium. Comprehensive spatial and spectral grid convergence studies have been implemented to quantify grid effects on the convective and radiative heating levels. Axisymmetric tests including the coupled effects of radiative energy transfer show negligible change to the convective heating and a 20 percent reduction in the radiative heating.

Gnoffo, Peter A.↗

The nonlinear evolution of inviscid Goertler vortices in three-dimensional boundary layers

The nonlinear development of inviscid Gortler vortices in a three-dimensional boundary layer is considered. We do not follow the classical approach of weakly nonlinear stability problems and consider a mode which has just become unstable. Instead we extend the method of Blackaby, Dando, and Hall (1992), which considered the closely related nonlinear development of disturbances in stratified shear flows. The Gortler modes we consider are initially fast growing and we assume, following others, that boundary-layer spreading results in them evolving in a linear fashion until they reach a stage where their amplitudes are large enough and their growth rates have diminished sufficiently so that amplitude equations can be derived using weakly nonlinear and non-equilibrium critical-layer theories. From the work of Blackaby, Dando and Hall (1993) is apparent, given the range of parameters for the Gortler problem, that there are three possible nonlinear integro-differential evolution equations for the disturbance amplitude. These are a cubic due to viscous effects, a cubic which corresponds to the novel mechanism investigated in this previous paper, and a quintic. In this paper we shall concentrate on the two cubic integro-differential equations and in particular, on the one due to the novel mechanism as this will be the first to affect a disturbance. It is found that the consideration of a spatial evolution problem as opposed to temporal (as was considered in Blackaby, Dando, and Hall, 1992) causes a number of significant changes to the evolution equations.

Blackaby, Nicholas↗

Innovative techniques for the production of energetic radicals for lunar processing including cold plasma processing of local planetary ores

Hydrogen reduction of ilmenite has been studied by a number of investigators as a potential means for recovery of oxygen from lunar soil. Interest in this process has always rested with the simplicity of the flow diagram and the utilization of established technology. Effective utilization of hydrogen in the reduction process at temperatures of 1200 C and below has always been disappointing and, as such, has led other investigators to focus attention on other systems. Effective utilization of hydrogen in the reduction of ilmenite can be significantly enhanced in the presence of a non-equilibrium hydrogen plasma. Ilmenite at solid specimen temperatures of 600 C to 970 C were reacted in a hydrogen plasma. Those experiments revealed that hydrogen utilization can be significantly enhanced. At a specimen temperature of 850 C the fraction of H2 reacted was 24 percent compared to the 7 percent theoretical limit calculated with thermodynamic theory for the same temperature. An added advantage for a hydrogen plasma involves further reduction of TiO2. Reduction of the iron oxide in ilmenite yields TiO2 and metallic iron as by products. Titanium forms a number of oxides including TiO, Ti2O3, Ti3O5 and the Magneli oxides (Ti4O7 to Ti50O99). In conventional processing of ilmenite with hydrogen it is possible to reduce TiO2 to Ti7O13 within approximately an hour, but with poor utilization of hydrogen on the order of one mole of H2 per thousand. In the cold or non-equilibrium plasma TiO2 can be rapidly reduced to Ti2O3 with hydrogen utilization exceeding 10 percent. Based on design considerations of the plasma reactor greater utilization of the hydrogen in the reduction of TiO2 is possible.

Bullard, D.↗

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.↗

A Technique for Determining the Nozzle-Flow Properties of Air in an Equilibrium, Nonequilibrium, or Frozen State

One of the problems associated with the development of high-enthalpy hypersonic test facilities hos been the determination of nozzle flow properties. This note presents a technique for determining the flow properties of air in an equilibrium, non-equilibrium, or frozen state using two test section measurements. A knowledge of two stagnation properties is required before the technique can be applied. The stagnation pressure can be measured, and reference 1 offers a method for determining stagnation enthalpy. Results of the calculations using the method of reference 1 are presented in reference 2 in a chart (chart 21) which can be readily used to determine the stagnation enthalpy from measurements of stagnation pressure, mass flow of air through the tunnel, and throat diameter.

Hayman, Lovick O.↗

RTE: A computer code for Rocket Thermal Evaluation

The numerical model for a rocket thermal analysis code (RTE) is discussed. RTE is a comprehensive thermal analysis code for thermal analysis of regeneratively cooled rocket engines. The input to the code consists of the composition of fuel/oxidant mixture and flow rates, chamber pressure, coolant temperature and pressure. dimensions of the engine, materials and the number of nodes in different parts of the engine. The code allows for temperature variation in axial, radial and circumferential directions. By implementing an iterative scheme, it provides nodal temperature distribution, rates of heat transfer, hot gas and coolant thermal and transport properties. The fuel/oxidant mixture ratio can be varied along the thrust chamber. This feature allows the user to incorporate a non-equilibrium model or an energy release model for the hot-gas-side. The user has the option of bypassing the hot-gas-side calculations and directly inputting the gas-side fluxes. This feature is used to link RTE to a boundary layer module for the hot-gas-side heat flux calculations.

Naraghi, Mohammad H. N.↗

Experimental and analytical study of cryogenic propellant boiloff to develop and verify alternate pressurization concepts for Space Shuttle external tank using a scaled down tank

Self pressurization by propellant boiloff is experimentally studied as an alternate pressurization concept for the Space Shuttle external tank (ET). The experimental setup used in the study is an open flow system which is composed of a variable area test tank and a recovery tank. The vacuum jacketed test tank is geometrically similar to the external LOx tank for the Space Shuttle. It is equipped with instrumentation to measure the temperature and pressure histories within the liquid and vapor, and viewports to accommodate visual observations and Laser-Doppler Anemometry measurements of fluid velocities. A set of experiments were conducted using liquid Nitrogen to determine the temperature stratification in the liquid and vapor, and pressure histories of the vapor during sudden and continuous depressurization for various different boundary and initial conditions. The study also includes the development and calibration of a computer model to simulate the experiments. This model is a one-dimensional, multi-node type which assumes the liquid and the vapor to be under non-equilibrium conditions during the depressurization. It has been tested for a limited number of cases. The preliminary results indicate that the accuracy of the simulations is determined by the accuracy of the heat transfer coefficients for the vapor and the liquid at the interface which are taken to be the calibration parameters in the present model.

Akyuzlu, K. M.↗

High-temperature LDV seed particle development

The feasibility of developing a method for making monodisperse, unagglomerated spherical particles greater than 50 nm in diameter was demonstrated. Carbonaceous particles were made by pyrolyzing ethylene with a pulsed CO2 laser, thereby creating a non-equilibrium mixture of carbon, hydrogen, hydrocarbon vapors, and unpyrolyzed ethylene. Via a complex series of reactions, the carbon and hydrocarbon vapors quickly condensed into the spherical particles. By cooling and dispersing them in a supersonic expansion immediately after their creation, the hot newly-formed spheres were prevented from colliding and coalescing, thus preventing the problem of agglomeration which as plagued other investigators studying laser-simulated particle formation. The cold particles could be left suspended in the residual gases indefinitely without agglomerating. Their uniform sizes and unagglomerated nature were visualized by collecting the particles on filters that were subsequently examined using electron microscopy. It was found the mean particle size can be coarsely controlled by varying the initial ethylene pressure, and can be finely controlled by varying the fluence (energy/unit area) with which the laser irradiates the gas. The motivating application for this research was to manufacture particles that could be used as laser Doppler velocimetry (LDV) seeds in high-temperature high-speed flows. Though the particles made in this program will not evaporate until heated to about 3000 K, and thus could serve as LDV seeds in some applications, they are not ideal when the hot atmosphere is also oxidizing. In that situation, ceramic materials would be preferable. Research performed elsewhere has demonstrated that selected ceramic materials can be manufactured by laser pyrolysis of appropriate supply gases. It is anticipated that, when the same gases are used in conjunction with the rapid cooling technique, unagglomerated spherical ceramic particles can be made with little difficulty. Such particles would also be valuable to manufacturers of ceramic or abrasive products, and this technique may find its greatest commercial potential in those areas.

Frish, Michael B.↗

The nonlinear evolution of modes on unstable stratified shear layers

The nonlinear development of disturbances in stratified shear flows (having a local Richardson number of value less than one quarter) is considered. Such modes are initially fast growing but, like related studies, we assume that the viscous, non-parallel spreading of the shear layer results in them evolving in a linear fashion until they reach a position where their amplitudes are large enough and their growth rates have diminished sufficiently so that amplitude equations can be derived using weakly nonlinear and non-equilibrium critical-layer theories. Four different basic integro-differential amplitude equations are possible, including one due to a novel mechanism; the relevant choice of amplitude equation, at a particular instance, being dependent on the relative sizes of the disturbance amplitude, the growth rate of the disturbance, its wavenumber, and the viscosity of the fluid. This richness of choice of possible nonlinearities arises mathematically from the indicial Frobenius roots of the governing linear inviscid equation (the Taylor-Goldstein equation) not, in general, differing by an integer. The initial nonlinear evolution of a mode will be governed by an integro-differential amplitude equations with a cubic nonlinearity but the resulting significant increase in the size of the disturbance's amplitude leads on to the next stage of the evolution process where the evolution of the mode is governed by an integro-differential amplitude equations with a quintic nonlinearity. Continued growth of the disturbance amplitude is expected during this stage, resulting in the effects of nonlinearity spreading to outside the critical level, by which time the flow has become fully nonlinear.

Blackaby, Nicholas↗

Supersonic flow of chemically reacting gas-particle mixtures. Volume 1: A theoretical analysis and development of the numerical solution

A numerical solution for chemically reacting supersonic gas-particle flows in rocket nozzles and exhaust plumes was described. The gas-particle flow solution is fully coupled in that the effects of particle drag and heat transfer between the gas and particle phases are treated. Gas and particles exchange momentum via the drag exerted on the gas by the particles. Energy is exchanged between the phases via heat transfer (convection and/or radiation). Thermochemistry calculations (chemical equilibrium, frozen or chemical kinetics) were shown to be uncoupled from the flow solution and, as such, can be solved separately. The solution to the set of governing equations is obtained by utilizing the method of characteristics. The equations cast in characteristic form are shown to be formally the same for ideal, frozen, chemical equilibrium and chemical non-equilibrium reacting gas mixtures. The particle distribution is represented in the numerical solution by a finite distribution of particle sizes.

Penny, M. M.↗

A survey of upwind methods for flows with equilibrium and non-equilibrium chemistry and thermodynamics

Several versions of flux-vector split and flux-difference split algorithms were compared with regard to general applicability and complexity. Test computations were performed using curve-fit equilibrium air chemistry for an M = 5 high-temperature inviscid flow over a wedge, and an M = 24.5 inviscid flow over a blunt cylinder for test computations; for these cases, little difference in accuracy was found among the versions of the same flux-split algorithm. For flows with nonequilibrium chemistry, the effects of the thermodynamic model on the development of flux-vector split and flux-difference split algorithms were investigated using an equilibrium model, a general nonequilibrium model, and a simplified model based on vibrational relaxation. Several numerical examples are presented, including nonequilibrium air chemistry in a high-temperature shock tube and nonequilibrium hydrogen-air chemistry in a supersonic diffuser.

Grossman, B.↗

Upwind schemes and bifurcating solutions in real gas computations

The area of high speed flow is seeing a renewed interest due to advanced propulsion concepts such as the National Aerospace Plane (NASP), Space Shuttle, and future civil transport concepts. Upwind schemes to solve such flows have become increasingly popular in the last decade due to their excellent shock capturing properties. In the first part of this paper the authors present the extension of the Osher scheme to equilibrium and non-equilibrium gases. For simplicity, the source terms are treated explicitly. Computations based on the above scheme are presented to demonstrate the feasibility, accuracy and efficiency of the proposed scheme. One of the test problems is a Chapman-Jouguet detonation problem for which numerical solutions have been known to bifurcate into spurious weak detonation solutions on coarse grids. Results indicate that the numerical solution obtained depends both on the upwinding scheme used and the limiter employed to obtain second order accuracy. For example, the Osher scheme gives the correct CJ solution when the super-bee limiter is used, but gives the spurious solution when the Van Leer limiter is used. With the Roe scheme the spurious solution is obtained for all limiters.

Suresh, Ambady↗

Study and modeling of finite rate chemistry effects in turbulent non-premixed flames

The development of numerical models that reflect some of the most important features of turbulent reacting flows requires information about the behavior of key quantities in well defined combustion regimes. In turbulent flames, the coupling between turbulent and chemical processes is so strong that it is extremely difficult to isolate the role played by one individual physical phenomenon. Direct numerical simulation (hereafter DNS) allows us to study in detail the turbulence-chemistry interaction in some restricted but completely defined situations. Globally, non-premixed flames are controlled by two limiting regimes: the fast chemistry case, where the turbulent flame can be pictured as a random distribution of local chemical equilibrium problems; and the slow chemistry case, where the chemistry integrates in time the turbulent fluctuations. The Damkoehler number, ratio of a mechanical time scale to chemical time scale, is used to distinguish between these regimes. Today most of the industrial computer codes are able to perform predictions in the hypothesis of local equilibrium chemistry using a presumed shape for the probability density function (pdt) of the conserved scalar. However, the finite rate chemistry situation is of great interest because industrial burners usually generate regimes in which, at some points, the flame is undergoing local extinction or at least non-equilibrium situations. Moreover, this variety of situations strongly influences the production of pollutants. To quantify finite rate chemistry effect, the interaction between a non-premixed flame and a free decaying turbulence is studied using DNS. The attention is focused on the dynamic of extinction, and an attempt is made to quantify the effect of the reaction on the small scale mixing process. The unequal diffusivity effect is also addressed. Finally, a simple turbulent combustion model based on the DNS observations and tractable in real flow configurations is proposed.

Vervisch, Luc↗