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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

Evaluation of Rotordynamic Coefficients of Look-through Labyrinths by Means of a Three Volume Bulk Flow Model

To describe the compressible, turbulent flow in a labyrinth seal, a three volume bulk flow model is presented. The conservation equations for mass, momentum, and energy are established in every control volume. A perturbation analysis is performed, yielding zeroth order equations for centric rotor position and first order equations describing the flow field for small rotor motions around the seal center. The equations are integrated numerically. From perturbation pressure, the forces on the shaft and the dynamic coefficients are calculated.

Nordmann, R.↗
NASA NTRS1992

Development of a 3-D upwind PNS code for chemically reacting hypersonic flowfields

Two new parabolized Navier-Stokes (PNS) codes were developed to compute the three-dimensional, viscous, chemically reacting flow of air around hypersonic vehicles such as the National Aero-Space Plane (NASP). The first code (TONIC) solves the gas dynamic and species conservation equations in a fully coupled manner using an implicit, approximately-factored, central-difference algorithm. This code was upgraded to include shock fitting and the capability of computing the flow around complex body shapes. The revised TONIC code was validated by computing the chemically-reacting (M(sub infinity) = 25.3) flow around a 10 deg half-angle cone at various angles of attack and the Ames All-Body model at 0 deg angle of attack. The results of these calculations were in good agreement with the results from the UPS code. One of the major drawbacks of the TONIC code is that the central-differencing of fluxes across interior flowfield discontinuities tends to introduce errors into the solution in the form of local flow property oscillations. The second code (UPS), originally developed for a perfect gas, has been extended to permit either perfect gas, equilibrium air, or nonequilibrium air computations. The code solves the PNS equations using a finite-volume, upwind TVD method based on Roe's approximate Riemann solver that was modified to account for real gas effects. The dissipation term associated with this algorithm is sufficiently adaptive to flow conditions that, even when attempting to capture very strong shock waves, no additional smoothing is required. For nonequilibrium calculations, the code solves the fluid dynamic and species continuity equations in a loosely-coupled manner. This code was used to calculate the hypersonic, laminar flow of chemically reacting air over cones at various angles of attack. In addition, the flow around the McDonnel Douglas generic option blended-wing-body was computed and comparisons were made between the perfect gas, equilibrium air, and the nonequilibrium air results.

Tannehill, J. C.↗
NASA NTRS1990

Computer program to solve two-dimensional shock-wave interference problems with an equilibrium chemically reacting air model

The computer program EASI, an acronym for Equilibrium Air Shock Interference, was developed to calculate the inviscid flowfield, the maximum surface pressure, and the maximum heat flux produced by six shock wave interference patterns on a 2-D, cylindrical configuration. Thermodynamic properties of the inviscid flowfield are determined using either an 11-specie, 7-reaction equilibrium chemically reacting air model or a calorically perfect air model. The inviscid flowfield is solved using the integral form of the conservation equations. Surface heating calculations at the impingement point for the equilibrium chemically reacting air model use variable transport properties and specific heat. However, for the calorically perfect air model, heating rate calculations use a constant Prandtl number. Sample calculations of the six shock wave interference patterns, a listing of the computer program, and flowcharts of the programming logic are included.

Glass, Christopher E.↗
NASA NTRS1992

Computation of the flow field in a centrifugal impeller with splitter blades

To support the design effort of the Space Transportation Main Engine (STME) Fuel Pump Stage, viscous flow calculations were performed in a centrifugal impeller with splitter blades. These calculations were carried out with a Navier-Stokes solver (MINT), which employs a linearized block-implicit Alternating Direction Implicit (ADI) procedure to iteratively solve a finite difference form of the system of conservation equations of mass, momentum, and energy in body-fitted coordinates. A computational grid was generated algebraically for the 'channel' between two main blades of the impeller and extended both upstream of the impeller inlet and downstream of the impeller exit so that the appropriate boundary conditions could be applied. The results of the calculations show that although the overall level of flow distortion near the impeller exit is not very large, there is a noticeable difference between the flow patterns in the two 'passages' (one passage between the pressure side of the splitter blade and the suction side of the next full blade).

Dejong, Frederik J.↗
NASA NTRS1992

Combustion instability analysis for liquid propellant rocket engines

The multi-dimensional numerical model has been developed to analyze the nonlinear combustion instabilities in liquid-fueled engines. The present pressure-based approach can handle the implicit pressure-velocity coupling in a non-iterative way. The additional scalar conservation equations for the chemical species, the energy, and the turbulent transport quantities can be handled by the same predictor-corrector sequences. This method is time-accurate and it can be applicable to the all-speed, transient, multi-phase, and reacting flows. Special emphasis is given to the acoustic/vaporization interaction which may act as the crucial rate-controlling mechanism in the liquid-fueled rocket engines. The subcritical vaporization is modeled to account for the effects of variable thermophysical properties, non-unitary Lewis number in the gas-film, the Stefan flow effect, and the effect of transient liquid heating. The test cases include the one-dimenisonal fast transient non-reacting and reacting flows, and the multi-dimensional combustion instabilities encountered in the liquid-fueled rocket thrust chamber. The present numerical model successfully demonstrated the capability to simulate the fast transient spray-combusting flows in terms of the limiting-cycle amplitude phenomena, correspondence between combustion and acoustics, and the steep-fronted wave and flame propagation. The investigated parameters include the spray initial conditions, air-fuel mixture ratios, and the engine geometry. Stable and unstable operating conditions are found for the liquid-fueled combustors. Under certain conditions, the limiting cycle behavior of the combusting flowfields is obtained. The numerical results indicate that the spray vaporization processes play an important role in releasing thermal energy and driving the combustion instability.

Kim, Y. M.↗
NASA NTRS1991

The effect of surface radiation on flame spread in a quiescent, microgravity environment

In most theories of laminar flame spread over solid combustibles, radiation is neglected. However, when gas motion is completely absent, which may occur in the quiescent, microgravity environment of a spacecraft, the importance of radiation compared to convection is enhanced. Here, a theoretical model is developed of flame spread over a thin solid fuel into an opposing flow of oxidizer, the opposing flow being present in the quiescent environment in flame-fixed coordinates, including the effects of radiative heat transfer from the fuel surface. Numerical solutions to the conservation equations in the gas and the solid phase that describe the spreading flame are presented for a variety of ambient conditions and surface emittances. When surface radiation is significant, the solid surface acts as a heat sink, causing the flame to cool, shrink in size, and spread more slowly. Better agreement with experiment for net heat transfer to the surface is obtained when surface radiation is included. For low ambient oxygen levels, radiative effects lead to flame extinction. A radiation/conduction parameter is identified that adequately describes the importance of surface radiation.

Bhattacharjee, S.↗
NASA NTRS1991

Opposed Flow Flame Spread in Normal, Enhanced and Reduced Gravity

Experimental and theoretical aspects of opposed-flow flame spread over solid fuels are presented with emphasis on the microgravity environments of spacecraft. For high opposing flow velocities, spread rate decreases with increasing velocity eventually leading to flame blowoff due to kinetic effects. At low opposing flow velocities, where diffusional effects are slowed and radiation becomes important, flame spread rate increases with increasing flow velocity. Extinction at low velocities is due to radiative effects. Modeling efforts that include radiation, both solid surface radiation and gas-phase radiation, predict qualitatively the experimental trends observed. Computationally, gas-phase radiation is conveniently included in solution of the conservation equations by employing a Plank mean absorption coefficient, a fraction of radiation that is fed back to the surface, and a shape function that describes the radiative flux distribution along the surface.

Altenkirch, Robert A.↗
NASA NTRS1990

Global analyses of water vapor, cloud and precipitation derived from a diagnostic assimilation of SSM/I geophysical retrievals

An analytical approach is described for diagnostically assimilating moisture data from Special Sensor Microwave Imager (SSM/I) into a global analysis of water vapor, cloud content, and precipitation. In this method, 3D fields of wind and temperature values taken from ECMWF gridded analysis are used to drive moisture conservation equations with parameterized microphysical treatment of vapor, liquid, and ice; the evolving field of water vapor is periodically updated or constrained by SSM/I retrievals of precipitable water. Initial results indicate that this diagnostic model can produce realistic large-scale fields of cloud and precipitation. The resulting water vapor analyses agree well with SSM/I and have an additional advantage of being synoptic.

Robertson, Franklin R.↗
NASA NTRS1992

Improved nonequilibrium viscous shock-layer scheme for hypersonic blunt-body flowfields

The nonequilibrium viscous shock-layer (VSL) solution scheme is revisited to improve its solution accuracy in the stagnation region and also to minimize and control errors in the conservation of elemental mass. The stagnation-point solution is improved by using a second-order expansion for the normal velocity, and the elemental mass conservation is improved by directly imposing the element conservation equations as solution constraints. These modifications are such that the general structure and computational efficiency of the nonequilibrium VSL scheme is not affected. This revised nonequilibrium VSL scheme is used to study the Mach 20 flow over a 7-deg sphere-cone vehicle under 0- and 20-deg angle-of-attack conditions. Comparisons are made with the corresponding predictions of Navier-Stokes and parabolized Navier-Stokes solution schemes. The results of these tests show that the nonequilibrium blunt-body VSL scheme is indeed an accurate, fast, and extremely efficient means for generating the blunt-body flowfield over spherical nose tips at zero-to-large angles of attack.

Bhutta, Bilal A.↗
NASA NTRS1992

Self-pressurization of a spherical liquid hydrogen storage tank in a microgravity environment

Thermal stratification and self-pressurization of partially filled liquid hydrogen (LH2) storage tanks under microgravity condition is studied theoretically. A spherical tank is subjected to a uniform and constant wall heat flux. It is assumed that a vapor bubble is located in the tank center such that the liquid-vapor interface and tank wall form two concentric spheres. This vapor bubble represents an idealized configuration of a wetting fluid in microgravity conditions. Dimensionless mass and energy conservation equations for both vapor and liquid regions are numerically solved. Coordinate transformation is used to capture the interface location which changes due to liquid thermal expansion, vapor compression, and mass transfer at liquid-vapor interface. The effects of tank size, liquid fill level, and wall heat flux on the pressure rise and thermal stratification are studied. Liquid thermal expansion tends to cause vapor condensation and wall heat flux tends to cause liquid evaporation at the interface. The combined effects determine the direction of mass transfer at the interface. Liquid superheat increases with increasing wall heat flux and liquid fill level and approaches an asymptotic value.

Lin, C. S.↗
NASA NTRS1992

Low-to-high altitude predictions of three-dimensional ablative reentry flowfields

A new coupled-chemistry PNS scheme has been developed to study 3D hypersonic reentry flows over sphere cones at small angles of attack. A new quasi-steady ablation model for Teflon was used to predict ablative reentry flow fields under laminar, transitional, and turbulent conditions. This solution scheme uses a new two-step solution scheme for solving the species conservation equations under low-altitude conditions. As an example, the scheme is used to study hypersonic flow over a 6-deg sphere-cone configuration at 50 and 125 kft conditions under zero and nonzero angle of attack conditions.

Bhutta, Bilal A.↗
NASA NTRS1992

Variable gamma computations for a forward-firing Space Shuttle Main Engine

A procedure for solving an axisymmetric SSME nozzle exhausting with an equivalent mass flow of two nozzles is presented for the return-to-launch-site (RTLS) maneuver for an engine-out case. The algorithm incorporates a two-species variable gamma formulation along with diffusion terms associated with a species conservation equation. Coupling of this algorithm is accomplished explicitly to the Navier-Stokes PARC2DR and PARC3DR codes. Axisymmetric results indicate a very complicated flowfield with a bow shock extending to approximately 570 radii in front of the nozzle, and low temperatures on what would be an external tank (ET) aft of the nozzle.

Saladino, A. J.↗
NASA NTRS1992

Semi-coupled flow and thermal analysis of the field joint during rapid pressurization of the redesigned solid rocket motor

A transient, semi-coupled, multi-dimensional thermal and flow analysis methodology was developed to predict the thermal/gas dynamic conditions in the field joint region of the Redesigned Solid Rocket Motor (RSRM). Transient temperature response, pressure history, and combustion gas flow rate (within the field joint region), were of principle interest, in the course of this study. The thermal environment in the field joint was modeled using SINDA, a finite difference based thermal network analyzer. The combustion gas flow boundary condition was generated using the FLAP code; this code performs a transient, lumped-parameter, control volume analysis to solve the mass, momentum, and energy conservation equations. The FLAP computer code was modified to account for erosion of the NBR insulation material, following ignition. An independent grid sensitivity study was conducted to determine an appropriate grid distribution near the wall. The predicted results, obtained using an optimum grid distribution and computer generated flow boundary condition, were compared with subscale test data.

Ghaffarian, Benny↗
NASA NTRS1992

Magnetic fields and radiative shocks in protogalaxies and the origin of globular clusters

The paper examines the hypothesis that globular clusters formed from gravitational instability in dense sheets of gas produced behind radiative shocks inside protogalaxies, such as those produced by the collision of subgalactic mass fragments partaking of the virial motions within the protogalaxy, in order to determine the differences which result if a magnetic field is present in the preshock medium. The MHD conservation equations are solved along with rate equations for nonequilibrium ionization, recombination, molecular formation and dissociation, and the equations of radiative transfer for steady-state shocks of velocity 300 km/s in a gas of preshock densities of 0.1-1 cu cm, and magnetic field strengths of 0.1-6 micro-G. The magnetic field is found to limit the degree of postshock compression and, thereby, to reduce the level of external radiation flux required to suppress H2 formation and cooling.

Shapiro, Paul R.↗
NASA NTRS1992

Combustion of liquid-fuel droplets in supercritical conditions

A comprehensive analysis of liquid-fuel droplet combustion in both subcritical and supercritical environments has been conducted. The formulation is based on the complete conservation equations for both gas and liquid phases, and accommodates variable thermophysical properties, finite-rate chemical kinetics, and a full treatment of liquid-vapor phase equilibrium at the drop surface. The governing equations and associated interfacial boundary conditions are solved numerically using a fully coupled, implicit scheme with the dual time-stepping integration technique. The model is capable of treating the entire droplet history, including the transition from the subcritical to supercritical state. As a specific example, the combustion of n-pentane fuel droplets in air is studied for pressures in the range of 5-140 atm. Results indicate that the ambient gas pressure exerts significant control of droplet gasification and burning processes through its influence on fluid transport, gas-liquid interfacial thermodynamics, and chemical reactions. The droplet gasification rate increases progressively with pressure. However, the data for the overall burnout time exhibit a considerable change in the combustion mechanism at the critical pressure, mainly as a result of reduced mass diffusivity and latent heat of vaporization with increased pressure.

Shuen, J. S.↗
NASA NTRS1992

The structure of premixed particle-cloud flames

The structure of premixed flames propagating in combustible systems, containing uniformly distributed volatile fuel particles, in an oxidizing gas mixture, is analyzed. It is presumed that the fuel particles vaporize first to yield a gaseous fuel of known chemical structure, which is subsequently oxidized in the gas phase. The analysis is performed in the asymptotic limit, where the value of the characteristic Zeldovich number, based on the gas-phase oxidation of the gaseous fuel is large, and for values of phi(u) greater than or equal to 1.0, where phi(u) is the equivalence ratio based on the fuel available in the fuel particles. The structure of the flame is presumed to consist of a preheat vaporization zone where the rate of the gas-phase chemical reaction is small, a reaction zone where convection and the rate of vaporization of the fuel particles are small and a convection zone where diffusive terms in the conservation equations are small. For given values phi(u) the analysis yields results for the burning velocity and phi(g) where phi(g) is the effective equivalence ratio in the reaction zone. The analysis shows that even though phi(u) greater than or equal to 1.0, for certain cases the calculated value of phi(g) is less than unity. This prediction is in agreement with experimental observations.

Seshadri, K.↗
NASA NTRS1992

Properties of thermospheric gravity waves on earth, Venus and Mars

A spectral model with spherical harmonics and Fourier components that can simulate atmospheric perturbations in the global geometry of a multiconstituent atmosphere is presented. The boundaries are the planetary surface where the transport velocities vanish and the exobase where molecular heat conduction and viscosity dominate. The time consuming integration of the conservation equations is reduced to computing the transfer function (TF) which describes the dynamic properties of the medium divorced from the complexities in the temporal and horizontal variations of the excitation source. Given the TF, the atmospheric response to a chosen source distribution is then obtained in short order. Theoretical studies are presented to illuminate some properties of gravity waves on earth, Venus and Mars.

Mayr, H. G.↗
NASA NTRS1992

Calculating Shocks In Chemically Reacting Flows

Total-variation-diminishing dissipation operator results in crisp shocks. Report discusses method for numerical simulation of shocks in hypersonic, chemically reacting flows. Introduces mathematical model of gas of as many as three species for which species-conservation equations coupled with gas-dynamic equations.

Eberhardt, Scott↗