Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Conservation equations”

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
NASA NTRS1986

Computation of nonstationary strong shock diffraction by curved surfaces

A two-dimensional, high resolution shock-capturing algorithm was used on a supercomputer to solve Eulerian gasdynamic equations in order to simulate nonstationary strong shock diffraction by a circular arc model in a shock tube. The hypersonic Mach shock wave was assumed to arrive at a high angle of incidence, and attention was given to the effect of varying values of the ratio of specific heats on the shock diffraction process. Details of the conservation equations of the numerical algorithm, written in curvilinear coordinates, are provided, and model output is illustrated with the results generated for a Mach shock encountering a 15 deg circular arc. The sample graphics include isopycnics, a shock surface density profile, and pressure and Mach number contours.

Yang, J. Y.↗
NASA NTRS1972

Numerical solutions of reactive fluid flows during postignition transients in hybrid rocket systems.

A computational method has been developed for the study of the post-ignition transients in hybrid rocket systems. The particular system chosen consisted of a gaseous oxidizer flowing within a tube of solid fuel, resulting in heterogeneous combustion. With the appropriate assumptions, two-dimensional, time-dependent conservation equations were derived for the reacting gas phase, and for the solid phase, in a cylindrical coordinate system. These were then programmed for numerical computation, using two implicit finite-difference schemes, the Lax-Wendroff scheme for the gas phase, and the Crank-Nicolson scheme for the solid phase. Appropriate initial and boundary conditions were represented, including heat and mass conservation at the interface between gas and solid. Initially, no attempt was made to relate the recession rate at the surface to the surface temperature, or to include heat transfer by radiation. A simple case was selected for preliminary calculations, with aluminum and oxygen as fuel and oxidizer, and aluminum oxide as the product.

Hung, W. S. Y.↗
NASA NTRS1983

Nonlinear truncation error analysis of finite difference schemes for the Euler equations

It is pointed out that, in general, dissipative finite difference integration schemes have been found to be quite robust when applied to the Euler equations of gas dynamics. The present investigation considers a modified equation analysis of both implicit and explicit finite difference techniques as applied to the Euler equations. The analysis is used to identify those error terms which contribute most to the observed solution errors. A technique for analytically removing the dominant error terms is demonstrated, resulting in a greatly improved solution for the explicit Lax-Wendroff schemes. It is shown that the nonlinear truncation errors are quite large and distributed quite differently for each of the three conservation equations as applied to a one-dimensional shock tube problem.

Klopfer, G. H.↗
NASA NTRS1973

Supersonic turbulent boundary-layer flows with tangential slot injection.

An implicit finite-difference method has been developed for the solution of the compressible boundary-layer equations. This method is applied to tangential slot injection into supersonic turbulent boundary-layer flows. In addition, the effects induced by the interaction between the boundary-layer displacement thickness and the external pressure field are considered. Three different eddy viscosity models have been used to specify the turbulent momentum exchange. One model depends on the species concentration profile, and the species conservation equation has been included in the system of governing partial differential equations. For air injected into air, the freestream and injected gases are treated as separate species which have common fluid properties. Calculations were made and results were compared with experimental data at freestream Mach numbers of 2.4 and 6.0 and with results of another finite-difference method. Good agreement was obtained for the reduction of wall skin friction with slot injection.

Miner, E. W.↗
NASA NTRS1997

Impact of Periodic Unsteadiness on Performance and Heat Load in Axial Flow Turbomachines

Results of an analytical and experimental investigation, directed at the understanding of the impact of periodic unsteadiness on the time-averaged flows in axial flow turbomachines, are presented. Analysis of available experimental data, from a large-scale rotating rig (LSRR) (low speed rig), shows that in the time-averaged axisymmetric equations the magnitude of the terms representing the effect of periodic unsteadiness (deterministic stresses) are as large or larger than those due to random unsteadiness (turbulence). Numerical experiments, conducted to highlight physical mechanisms associated with the migration of combustor generated hot-streaks in turbine rotors, indicated that the effect can be simulated by accounting for deterministic stress like terms in the time-averaged mass and energy conservation equations. The experimental portion of this program shows that the aerodynamic loss for the second stator in a 1-1/2 stage turbine are influenced by the axial spacing between the second stator leading edge and the rotor trailing edge. However, the axial spacing has little impact on the heat transfer coefficient. These performance changes are believed to be associated with the change in deterministic stress at the inlet to the second stator. Data were also acquired to quantify the impact of indexing the first stator relative to the second stator. For the range of parameters examined, this effect was found to be of the same order as the effect of axial spacing.

Sharma, Om P.↗
NASA NTRS1995

Nearly sonic and transsonic convective motions in the solar atmosphere related to the solar wind origin

MHD equations are considered for the solar atmosphere. 15 different simplest MHD regimes are indicated for the momentum transport equation depending on the mutual binary interplay between 6 possible and locally dominant terms: non-stationarity and inhomogeneity of the flow, gas pressure, magnetic tensions, viscous and gravity forces. These regimes are delimited by five physically independent dimensionless parameters, for example, Strouhal, sonic Mach, alfvenic Mach, Reynolds and Froude numbers or their combinations. Another partially overlapping classification of the simplest regimes may be introduced based on the energy conservation equation. There are also 15 independent binary combinations between nonstationary and inhomogeneous convective terms, dissipative energy sinks and sources (viscous, heat-conductive, Joule and radiative ones) in the energy equation. More complicated regimes are considered with multiple dominated terms. All these MHD regimes play their important role somewhere in the solar atmosphere complicated by the tensor transport coefficients in the magnetically dominated regions of the upper atmosphere. Nearly sonic and transsonic nonstationary convective motions with ascending and descending flows are observed in the solar chromosphere. the transition region and the lower pans of the solar corona together with related horizontal velocity components. This convection represents a kind of the 'cocoonery' manufacturing nonstationary vortices generated here and partially connected to the photosphere and to the solar wind. The solar wind originates from this powerful transsonic muddle in the solar atmosphere as a tiny fraction of the streamlines which are temporarily getting detached from the 'cocoons; and going to the infinity. The topologically complicated instantaneous 'runaway surface' around the Sun, i.e., the surface which separates outgoing to the infinity streams from other finite flows in the solar atmosphere was not described in the literature and needs additional investigation. We conclude that a simple one-connected smooth and quasistationary 'critical surface' which is often supposed to be placed somewhere in the solar corona (say, at 4 solar radii) in reality does not exist. Observational and theoretical arguments favor instead of this a highly structured, disordered, patched and time-dependent sonic transition surface in the solar atmosphere permanently fluctuating at positional dependent heights starting sometimes from photospheric (and maybe subphotospheric) levels and more frequently from the chromosphenc network up to several solar radii in the solar corona.

Veselovsky, I. S.↗
NASA NTRS1984

A solution procedure for three-dimensional incompressible Navier-Stokes equation and its application

An implicit, finite-difference procedure is presented for numerically solving viscous incompressible flows. For convenience of applying the present method to three-dimensional problems, primitive variables, namely the pressure and velocities, are used. One of the major difficulties in solving incompressible flows that use primitive variables is caused by the pressure field solution method which is used as a mapping procedure to obtain a divergence-free velocity field. The present method is designed to accelerate the pressure-field solution procedure. This is achieved by the method of pseudocompressibility in which the time derivative pressure term is introduced into the mass conservation equation. The pressure wave propagation and the spreading of the viscous effect is investigated using simple test problems. The present study clarifies physical and numerical characteristics of the pseudo-compressible approach in simulating incompressible flows. Computed results for external and internal flows are presented to verify the present procedure. The present algorithm has been shown to be very robust and accurate if the selection of the pseudo-compressibility parameter has been made according to the guidelines given.

Kwak, D.↗
NASA NTRS1985

New integration techniques for chemical kinetic rate equations. 2: Accuracy comparison

A comparison of the accuracy of several techniques recently developed for solving stiff differential equations is presented. The techniques examined include two general purpose codes EEPISODE and LSODE developed for an arbitrary system of ordinary differential equations, and three specialized codes CHEMEQ, CREKID, and GCKP84 developed specifically to solve chemical kinetic rate equations. The accuracy comparisons are made by applying these solution procedures to two practical combustion kinetics problems. Both problems describe adiabatic, homogeneous, gas phase chemical reactions at constant pressure, and include all three combustion regimes: induction, heat release, and equilibration. The comparisons show that LSODE is the most efficient code - in the sense that it requires the least computational work to attain a specified accuracy level. An important finding is that an iterative solution of the algebraic enthalpy conservation equation for the temperature can be more accurate and efficient than computing the temperature by integrating its time derivative.

Radhakrishnan, K.↗
NASA NTRS1986

New integration techniques for chemical kinetic rate equations. II - Accuracy comparison

A comparison of the accuracy of several techniques recently developed for solving stiff differential equations is presented. The techniques examined include two general purpose codes EEPISODE and LSODE developed for an arbitrary system of ordinary differential equations, and three specialized codes CHEMEQ, CREKID, and GCKP84 developed specifically to solve chemical kinetic rate equations. The accuracy comparisons are made by applying these solution procedures to two practical combustion kinetics problems. Both problems describe adiabatic, homogeneous, gas phase chemical reactions at constant pressure, and include all three combustion regimes: induction heat release, and equilibration. The comparisons show that LSODE is the most efficient code - in the sense that it requires the least computational work to attain a specified accuracy level. An important finding is that an iterative solution of the algebraic enthalpy conservation equation for the temperature can be more accurate and efficient than computing the temperature by integrating its time derivative.

Radhakrishnan, K.↗
NASA NTRS2005

Incorporation of a Chemical Equilibrium Equation of State into LOCI-Chem

Renewed interest in development of advanced high-speed transport, reentry vehicles and propulsion systems has led to a resurgence of research into high speed aerodynamics. As this flow regime is typically dominated by hot reacting gaseous flow, efficient models for the characteristic chemical activity are necessary for accurate and cost effective analysis and design of aerodynamic vehicles that transit this regime. The LOCI-Chem code recently developed by Ed Luke at Mississippi State University for NASA/MSFC and used by NASA/MSFC and SSC represents an important step in providing an accurate, efficient computational tool for the simulation of reacting flows through the use of finite-rate kinetics [3]. Finite rate chemistry however, requires the solution of an additional N-1 species mass conservation equations with source terms involving reaction kinetics that are not fully understood. In the equilibrium limit, where the reaction rates approach infinity, these equations become very stiff. Through the use of the assumption of local chemical equilibrium the set of governing equations is reduced back to the usual gas dynamic equations, and thus requires less computation, while still allowing for the inclusion of reacting flow phenomenology. The incorporation of a chemical equilibrium equation of state module into the LOCI-Chem code was the primary objective of the current research. The major goals of the project were: (1) the development of a chemical equilibrium composition solver, and (2) the incorporation of chemical equilibrium solver into LOCI-Chem. Due to time and resource constraints, code optimization was not considered unless it was important to the proper functioning of the code.

Cox, Carey F.↗
NASA NTRS1976

On the energy deposition of photoelectrons in the atmosphere of Venus

Vertical components of photoelectron fluxes in the atmosphere of Venus are computed by solving an appropriate form of the Boltzmann equation in the cases where there is no flux of either photoelectrons or solar-wind particles across the ionopause and where photoelectrons are free to escape from the atmosphere. It is assumed that Venus has no magnetic field and that the atmosphere is composed of carbon dioxide, atomic oxygen, and helium. The results are plotted as a function of altitude for several energies in the range from 100 eV to a cutoff of the order of 1 eV. Heating rates for the two upper boundary conditions and the case of no spatial transport are determined which show that transport effects dominate at altitudes greater than about 200 km. Electron temperatures are calculated for the adopted model atmosphere and ionosphere by solving the pertinent conservation equation, and excitation rates are computed for the CO Cameron band as well as the CO2(+) A and B bands.

Mccormick, P. T.↗
NASA NTRS1993

A critical analysis of the accuracy of several numerical techniques for combustion kinetic rate equations

A detailed analysis of the accuracy of several techniques recently developed for integrating stiff ordinary differential equations is presented. The techniques include two general-purpose codes EPISODE and LSODE developed for an arbitrary system of ordinary differential equations, and three specialized codes CHEMEQ, CREK1D, and GCKP4 developed specifically to solve chemical kinetic rate equations. The accuracy study is made by application of these codes to two practical combustion kinetics problems. Both problems describe adiabatic, homogeneous, gas-phase chemical reactions at constant pressure, and include all three combustion regimes: induction, heat release, and equilibration. To illustrate the error variation in the different combustion regimes the species are divided into three types (reactants, intermediates, and products), and error versus time plots are presented for each species type and the temperature. These plots show that CHEMEQ is the most accurate code during induction and early heat release. During late heat release and equilibration, however, the other codes are more accurate. A single global quantity, a mean integrated root-mean-square error, that measures the average error incurred in solving the complete problem is used to compare the accuracy of the codes. Among the codes examined, LSODE is the most accurate for solving chemical kinetics problems. It is also the most efficient code, in the sense that it requires the least computational work to attain a specified accuracy level. An important finding is that use of the algebraic enthalpy conservation equation to compute the temperature can be more accurate and efficient than integrating the temperature differential equation.

Radhadrishnan, Krishnan↗
NASA NTRS1964

An inviscid model of the solar wind

Inviscid, thermal conductive, spherically symmetric solar wind model - model solution by integration of mass, momentum, and energy conservation equations

SOLAR WIND↗