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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
NASA NTRS1990

High speed turbulent reacting flows: DNS and LES

Work on understanding the mechanisms of mixing and reaction in high speed turbulent reacting flows was continued. Efforts, in particular, were concentrated on taking advantage of modern computational methods to simulate high speed turbulent flows. In doing so, two methodologies were used: large eddy simulations (LES) and direct numerical simulations (DNS). In the work related with LES the objective is to study the behavior of the probability density functions (pdfs) of scalar properties within the subgrid in reacting turbulent flows. The data base obtained by DNS for a detailed study of the pdf characteristics within the subgrid was used. Simulations are performed for flows under various initializations to include the effects of compressibility on mixing and chemical reactions. In the work related with DNS, a two-dimensional temporally developing high speed mixing layer under the influence of a second-order non-equilibrium chemical reaction of the type A + B yields products + heat was considered. Simulations were performed with different magnitudes of the convective Mach numbers and with different chemical kinetic parameters for the purpose of examining the isolated effects of the compressibility and the heat released by the chemical reactions on the structure of the layer. A full compressible code was developed and utilized, so that the coupling between mixing and chemical reactions is captured in a realistic manner.

Givi, Peyman↗
NASA NTRS1991

Mixing and non-equilibrium chemical reaction in a compressible mixing layer

The effects of compressibility, chemical reaction exothermicity, and non-equilibrium chemical modeling in a reacting plane mixing layer were investigated by means of two dimensional direct numerical simulations. The chemical reaction was irreversible and second order of the type A + B yields Products + Heat. The general governing fluid equations of a compressible reacting flow field were solved by means of high order finite difference methods. Physical effects were then determined by examining the response of the mixing layer to variation of the relevant non-dimensionalized parameters. The simulations show that increased compressibility generally results in a suppressed mixing, and consequently a reduced chemical reaction conversion rate. Reaction heat release was found to enhance mixing at the initial stages of the layer growth, but had a stabilizing effect at later times. The increased stability manifested itself in the suppression or delay of the formation of large coherent structures within the flow. Calculations were performed for a constant rate chemical kinetics model and an Arrhenius type kinetic prototype. The choice of the model was shown to have an effect on the development of the flow. The Arrhenius model caused a greater temperature increase due to reaction than the constant kinetic model. This had the same effect as increasing the exothermicity of the reaction. Localized flame quenching was also observed when the Zeldovich number was relatively large.

Steinberger, Craig J.↗
NASA NTRS1974

Termination of the solar wind in the hot, partially ionized interstellar medium

Theoretical foundations for understanding the problem of the termination of the solar wind are reexamined in the light of most recent findings concerning the states of the solar wind and the local interstellar medium. The investigation suggests that a simple extention of Parker's (1961) analytical model provides a useful approximate description of the combined solar wind, interstellar wind plasma flowfield under conditions presently thought to occur. A linear perturbation solution exhibiting both the effects of photoionization and charge exchange is obtained for the supersonic solar wind. A numerical algorithm is described for computing moments of the non-equilibrium hydrogen distribution function and associated source terms for the MHD equations. Computed using the algorithm in conjunction with the extended Parker solution to approximate the plasma flowfield, profiles of hydrogen number density are given in the solar wind along the upstream and downstream axes of flow with respect to the direction of the interstellar wind. Predictions of solar Lyman-alpha backscatter intensities to be observed at 1 a.u. have been computed, in turn, from a set of such hydrogen number density profiles varied over assumed conditions of the interstellar wind.

Lombard, C. K.↗
NASA NTRS1992

Model free simulations of a high speed reacting mixing layer

The effects of compressibility, chemical reaction exothermicity and non-equilibrium chemical modeling in a combusting plane mixing layer were investigated by means of two-dimensional model free numerical simulations. It was shown that increased compressibility generally had a stabilizing effect, resulting in reduced mixing and chemical reaction conversion rate. The appearance of 'eddy shocklets' in the flow was observed at high convective Mach numbers. Reaction exothermicity was found to enhance mixing at the initial stages of the layer's growth, but had a stabilizing effect at later times. Calculations were performed for a constant rate chemical rate kinetics model and an Arrhenius type kinetics prototype. The Arrhenius model was found to cause a greater temperature increase due to reaction than the constant kinetics model. This had the same stabilizing effect as increasing the exothermicity of the reaction. Localized flame quenching was also observed when the Zeldovich number was relatively large.

Steinberger, Craig J.↗
NASA NTRS1992

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↗
NASA NTRS1992

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↗
NASA NTRS1981

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

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

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

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

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

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

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

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

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↗
NASA NTRS1976

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

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

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↗