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 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
NASA NTRS1995

OMV/VTE variable thrust engine analysis

The objective of the present work is to develop a predictive CFD based analytical tool for the Variable Thrust Engine (VTE) in the Orbital Maneuvering Vehicle (OMV). This objective is being accomplished within the framework of the Los Alamos KIVA computer code for chemically reactive flows with sprays. For the OMV application, the main structure of KIVA is to be retained while reformulating many of the phenomenological submodels, enhancing some of the numerics, and adding more features. The analytical model consists of the general conservation equations for two-phase reactive flows and of submodels for turbulence, chemical reactions, and bipropellant sprays. Tailoring this model to the OMV engine brings about the added complexities of combustion and flow processes that occur in a liquid hypergolic propellant rocket chamber. This report exposes the foundation upon which the analytical tool is being constructed and developed. Results from a cursory computational exercise involving the simulation of the flow and combustion processes in a hypothetical N2H4/N204 rocket engine thrust chamber is presented and discussed.

Larosillere, L.↗
NASA NTRS1995

Anisotropic shock jump conditions: Theory and observations

The MHD Rankine-Hugoniot (RH) relations for shock waves in a collisionless plasma with bi-Maxwellian distribution functions are considered. While by introducing the pressure anisotropy parameter xi in the RH relations, the number of unknowns -- B, V, n, p and xi (a total of 9) -- becomes one more than the total number of the conservation equations, it is possible to use the observed quantities on both sides of the shock to study the anisotropy changes across the shock. A simple relation for the anisotropy change across the shock is derived as a function of the ratio of magnetic fields m(= B'/B), the shock normal angle theta(sub Bn) and the plasma beta and beta' (primes are downstream values). Since m and theta(sub Bn) can be determined accurately in observation, the reliability of the anisotropy change deduced is mostly dependent on the accuracy of the measurements beta and beta'. We have applied the results to six low-beta quasi-perpendicular (Q perpendicular) laminar bow shock crossings with temperature anisotropy measured in the magnetosheath. In the six test cases, it is found that the predicted pressure anisotropies agree well with those observed in the magnetosheath.

Chao, J. K.↗
NASA NTRS1989

OMV/VTE variable thrust engine analysis

The objective of the present work is to develop a predictive CFD based analytical tool for the variable thrust engine (VTE) in the orbital maneuvering vehicle (OMV). This objective is being accomplished within the framework of the Los Alamos KIVA computer code for chemically reactive flows with sprays. For the OMV application, the main structure of KIVA is to be retained while reformulating many of the phenomenological submodels, enhancing some of the numerics, and adding more features. The analytical model consists of the general conservation equations for two phase reactive flows and of submodels for turbulence, chemical reactions, and bipropellant sprays. Tailoring this model to the OMV engine brings about the added complexities of combustion and flow processes that occur in a liquid hypergolic propellant rocket chamber. This report exposes the foundation upon which the analytical tool is being constructed and developed. Results from a cursory computational exercise involving the simulation of the flow and combustion processes in a hypothetical N2H4/N2O4 rocket engine thrust chamber is presented and discussed.

Larosillere, Louis↗
NASA NTRS1996

Investigation of shock-induced combustion past blunt projectiles

A numerical study is conducted to simulate shock-induced combustion in premixed hydrogen-air mixtures at various free-stream conditions and parameters. Two-dimensional axisymmetric, reacting viscous flow over blunt projectiles is computed to study shock-induced combustion at Mach 5.11 and Mach 6.46 in hydrogen-air mixture. A seven-species, seven reactions finite rate hydrogen-air chemical reaction mechanism is used combined with a finite-difference, shock-fitting method to solve the complete set of Navier-Stokes and species conservation equations. The study has allowed an improved understanding of the physics of shock-induced combustion over blunt projectiles and the numerical results can now be explained more readily with one-dimensional wave-interaction model.

Ahuja, J. K.↗
NASA NTRS1995

Non-parabolic hydrodynamic formulations for the simulation of inhomogeneous semiconductor devices

Hydrodynamic models are becoming prevalent design tools for small scale devices and other devices in which high energy effects can dominate transport. Most current hydrodynamic models use a parabolic band approximation to obtain fairly simple conservation equations. Interest in accounting for band structure effects in hydrodynamic device simulation has begun to grow since parabolic models can not fully describe the transport in state of the art devices due to the distribution populating non-parabolic states within the band. This paper presents two different non-parabolic formulations of the hydrodynamic model suitable for the simulation of inhomogeneous semiconductor devices. The first formulation uses the Kane dispersion relationship (hk)(exp 2)/2m = W(1 + alpha(W)). The second formulation makes use of a power law ((hk)(exp 2)/2m = xW(sup y)) for the dispersion relation. Hydrodynamic models which use the first formulation rely on the binomial expansion to obtain moment equations with closed form coefficients. This limits the energy range over which the model is valid. The power law formulation readily produces closed form coefficients similar to those obtained using the parabolic band approximation. However, the fitting parameters (x,y) are only valid over a limited energy range. The physical significance of the band non-parabolicity is discussed as well as the advantages/disadvantages and approximations of the two non-parabolic models. A companion paper describes device simulations based on the three dispersion relationships: parabolic, Kane dispersion, and power low dispersion.

Smith, Arlynn W.↗
NASA NTRS1995

General Fluid System Simulation Program to Model Secondary Flows in Turbomachinery

The complexity and variety of turbomachinery flow circuits created a need for a general fluid system simulation program for test data anomaly resolution as well as design review. The objective of the paper is to present a computer program that has been developed to support Marshall Space Flight Center's turbomachinery internal flow analysis efforts. The computer program solves for the mass. energy and species conservation equation at each node and flow rate equation at each branch of the network by a novel numerical procedure which is a combination of both Newton-Ralphson and successive substitution method and uses a thermodynamic property program for computing real gas properties. A generalized, robust, modular, and 'user-friendly' computer program has been developed to model internal flow rates, pressures, temperatures, concentrations of gas mixtures and axial thrusts. The program can be used for any network for compressible and incompressible flows, choked flow, change of phase and gaseous mixturecs. The code has been validated by comparing the predictions with Space Shuttle Main Engine test data.

Majumdar, Alok K.↗
NASA NTRS1994

Numerical Computation of Flame Spread over a Thin Solid in Forced Concurrent Flow with Gas-phase Radiation

Excerpts from a paper describing the numerical examination of concurrent-flow flame spread over a thin solid in purely forced flow with gas-phase radiation are presented. The computational model solves the two-dimensional, elliptic, steady, and laminar conservation equations for mass, momentum, energy, and chemical species. Gas-phase combustion is modeled via a one-step, second order finite rate Arrhenius reaction. Gas-phase radiation considering gray non-scattering medium is solved by a S-N discrete ordinates method. A simplified solid phase treatment assumes a zeroth order pyrolysis relation and includes radiative interaction between the surface and the gas phase.

Jiang, Ching-Biau↗
NASA NTRS1995

Theory of the Effects of Small Gravitational Levels on Droplet Gasification

A mathematical model taking into account small (and constant) gravitational levels is developed for vaporization of an isolated liquid droplet suspended in a stagnant atmosphere. A goal of the present analysis is to see how small gravitational levels affect droplet gasification characteristics. Attention is focused upon determining the effects on gas-phase phenomena. The conservation equations arc normalized and nondimensionalized, and a small parameter that accounts for the effects of gravity is identified. This parameter is the square of the inverse of a Froude number based on the gravitational acceleration, the droplet radius, and a characteristic gas-phase velocity at the droplet surface. Asymptotic analyses are developed in terms of this parameter. In the analyses, different spatial regions are identified. Near a droplet, gravitational effects are negligible in the first approximation, and the flowfield is spherically symmetric to the leading order. Analysis shows, however, that outer zones exist where gravitational effects cannot be neglected; it is expected that a stagnation point will be present in an outer zone that is not present when gravity is totally absent. The leading order and higher-order differential equations for each zone are derived and solved. The solutions allow the effects of gravity on vaporization rates and temperature, velocity and species fields to be determined.

Beitelmal, A.↗
NASA NTRS1994

Numerical Simulation of Shock-Induced Combustion Past Blunt Bodies Using Shock-Fitting Technique

Two-dimensional axisymmetric, reacting viscous flow over blunt projectiles is computed to study shock induced combustion at Mach 5.11 and Mach 6.46 in hydrogen-air mixture. A finite-difference, shock-fitting method is used to solve the complete set of Navier Stokes and species conservation equations. In this approach, the bow shock represents a boundary of the computational domain and is treated as a discontinuity across which Rankine-Hugoniot conditions are applied. All interior details of the flow such as compression waves, reaction front, and the wall boundary layer are captured automatically in the solution. Since shock-fitting approach reduces the amount of artificial dissipation, all the intricate details of the flow are captured much more clearly than has been possible with the shock-capturing approach. This has allowed an improved understanding of the physics of shock-induced combustion over blunt projectiles and the numerical results can now be explained more readily with one dimensional wave-interaction model than before.

Ahuja, J. K.↗
NASA NTRS1998

Flow Simulation in Secondary Flow Passages of a Rocket Engine Turbopump

This paper presents application of a Generalized Fluid System Simulation Program (GFSSP) to model fluid flow in a very complex secondary flow circuit of a rocket engine turbopump, SIMPLEX. SIMPLEX was a low cost turbopump designed, manufactured and tested to demonstrate the ability to reduce the overall cost and design cycle time of turbomachinery. GFSSP is a general purpose computer program for analyzing flowrate, pressure, temperature and concentration distribution in a complex flow network. The program employs a finite volume formulation of mass, momentum and energy conservation equations in conjunction with thermodynamic equation of state of a real fluid. The secondary flow passages modeled include flow between the impeller shroud and housing, impeller back face and housing, flow through labyrinth seal and bearing, return line flow to inducer and flow adjacent to front and back face of the turbine. Axial load on the bearing are calculated from the predicted pressure on all axial faces of the rotor. The predicted pressure and temperature distributions compared well with the test data.

Schallhorn, Paul↗
NASA NTRS1979

Modeling of Radiative Heating in Base Region of Jovian Entry Probe

A theoretical model is derived to determine the average thermodynamic properties in the expanding region, recirculating region, recompression region, and neck region through application of one-dimensional conservation equations. Radiative transfer is calculated using spectrally detailed computer codes accounting for nonequilibrium. The results show that the most severe heating occurs immediately behind the frustum, and that the recompression and neck regions are the major sources of radiation that heats the base stagnation point, The radiation flux to the base point is slightly stronger with ablation than without. Its value is 0.11(43 P(sub b)/P(sub s))(sup 2) times that to the front stagnation point where the base pressure P(sub b) is defined as the average pressure in the recirculating region and P(sub s) is the front stagnation-point pressure. The time-integrated heat load to the base point is 18(43 P(sub b)/P(sub s))(sup 2) kJ/sq cm.

Park, Chul↗
NASA NTRS1998

Numerical Modeling of Pressurization of a Propellant Tank

An unsteady finite volume procedure has been developed to predict the history of pressure, temperature and mass flow rate of the pressurant and propellant during the expulsion of the propellant from a tank. The time dependent mass, momentum and energy conservation equations are solved at the ullage space. The model accounts for the change in the ullage volume due to expulsion of the propellant. It also accounts for the heat transfer from the tank wall and propellant to the ullage gas. The procedure was incorporated in the Generalized Fluid System Simulation Program (GFSSP). The results of several test cases were then compared with a published correlation of pressurant requirements for a given displacement of propellant. The agreement between the predictions and the correlation was found to be satisfactory.

Majumdar, Alok↗
NASA NTRS1998

Current and Future Critical Issues in Rocket Propulsion Systems

The objective of this research was to tackle several problems that are currently of great importance to NASA. In a liquid rocket engine several complex processes take place that are not thoroughly understood. Droplet evaporation, turbulence, finite rate chemistry, instability, and injection/atomization phenomena are some of the critical issues being encountered in a liquid rocket engine environment. Pulse Detonation Engines (PDE) performance, combustion chamber instability analysis, 60K motor flowfield pattern from hydrocarbon fuel combustion, and 3D flowfield analysis for the Combined Cycle engine were of special interest to NASA. During the summer of 1997, we made an attempt to generate computational results for all of the above problems and shed some light on understanding some of the complex physical phenomena. For this purpose, the Liquid Thrust Chamber Performance (LTCP) code, mainly designed for liquid rocket engine applications, was utilized. The following test cases were considered: (1) Characterization of a detonation wave in a Pulse Detonation Tube; (2) 60K Motor wall temperature studies; (3) Propagation of a pressure pulse in a combustion chamber (under single and two-phase flow conditions); (4) Transonic region flowfield analysis affected by viscous effects; (5) Exploring the viscous differences between a smooth and a corrugated wall; and (6) 3D thrust chamber flowfield analysis of the Combined Cycle engine. It was shown that the LTCP-2D and LTCP-3D codes are capable of solving complex and stiff conservation equations for gaseous and droplet phases in a very robust and efficient manner. These codes can be run on a workstation and personal computers (PC's).

Navaz, Homayun K.↗
NASA NTRS1999

Numerical Speed of Sound and its Application to Schemes for all Speeds

The concept of "numerical speed of sound" is proposed in the construction of numerical flux. It is shown that this variable is responsible for the accurate resolution of' discontinuities, such as contacts and shocks. Moreover, this concept can he readily extended to deal with low speed and multiphase flows. As a results, the numerical dissipation for low speed flows is scaled with the local fluid speed, rather than the sound speed. Hence, the accuracy is enhanced the correct solution recovered, and the convergence rate improved. We also emphasize the role of mass flux and analyze the behavior of this flux. Study of mass flux is important because the numerical diffusivity introduced in it can be identified. In addition, it is the term common to all conservation equations. We show calculated results for a wide variety of flows to validate the effectiveness of using the numerical speed of sound concept in constructing the numerical flux. We especially aim at achieving these two goals: (1) improving accuracy and (2) gaining convergence rates for all speed ranges. We find that while the performance at high speed range is maintained, the flux now has the capability of performing well even with the low: speed flows. Thanks to the new numerical speed of sound, the convergence is even enhanced for the flows outside of the low speed range. To realize the usefulness of the proposed method in engineering problems, we have also performed calculations for complex 3D turbulent flows and the results are in excellent agreement with data.

Liou, Meng-Sing↗
NASA NTRS1999

Numerical Prediction of Transient Axial Thrust and Internal Flows in a Rocket Engine Turbopump

This paper presents the application of the Generalized Fluid System Simulation Program (GFSSP) to model the time-dependent flow in a complex secondary flow circuit of the turbopump of the Fastrac engine currently under development at Marshall Space Flight Center. GFSSP is a general purpose computer program for analyzing steady-state and time-dependant flowrates, pressures, temperatures, and concentrations in a complex flow network. The program employs a finite volume formulation of mass, momentum and energy conservation equations in conjunction with the thermodynamic equation of state of real fluids. GFSSP was used to calculate the axial thrust and internal flow distribution of the Fastrac engine turbopump during the start and shut down transients. The models discussed in this paper use boundary conditions that were extracted from turbopump test data. The GFSSP predicted turbopump secondary flow passage pressures and temperatures were compared with actual measured values.

VanHooser, Katherine↗
NASA NTRS1999

Physically-Derived Dynamical Cores in Atmospheric General Circulation Models

The algorithm chosen to represent the advection in atmospheric models is often used as the primary attribute to classify the model. Meteorological models are generally classified as spectral or grid point, with the term grid point implying discretization using finite differences. These traditional approaches have a number of shortcomings that render them non-physical. That is, they provide approximate solutions to the conservation equations that do not obey the fundamental laws of physics. The most commonly discussed shortcomings are overshoots and undershoots which manifest themselves most overtly in the constituent continuity equation. For this reason many climate models have special algorithms to model water vapor advection. This talk focuses on the development of an atmospheric general circulation model which uses a consistent physically-based advection algorithm in all aspects of the model formulation. The shallow-water model is generalized to three dimensions and combined with the physics parameterizations of NCAR's Community Climate Model. The scientific motivation for the development is to increase the integrity of the underlying fluid dynamics so that the physics terms can be more effectively isolated, examined, and improved. The expected benefits of the new model are discussed and results from the initial integrations will be presented.

Rood, Richard B.↗
NASA NTRS1998

Development of Numerical Extended Hydrodynamics for Transition-Regime Non-Equilibrium Flows Encountered in Semiconductor Manufacturing Processes

Six months of funding was received for the proposed three year research program (funding for the period from March 1, 1997 to August 31, 1997). Although the official starting date for the project was March 1, 1997, no funding for the project was received until July 1997. In the funded research period, considerable progress was made on Phase I of the proposed research program. The initial research efforts concentrated on applying the 10-, 20-, and 35-moment Gaussian-based closures to a series of standard two-dimensional non-reacting single species test flow problems, such as the flat plate, couette, channel, and rearward facing step flows, and to some other two-dimensional flows having geometries similar to those encountered in chemical-vapor deposition (CVD) reactors. Eigensystem analyses for these systems for the case of two spatial dimensions was carried out and efficient formulations of approximate Riemann solvers have been formulated using these eigenstructures. Formulations to include rotational non-equilibrium effects into the moment closure models for the treatment of polyatomic gases were explored, as the original formulations of the closure models were developed strictly for gases composed of monatomic molecules. The development of a software library and computer code for solving relaxing hyperbolic systems in two spatial dimensions of the type arising from the closure models was also initiated. The software makes use of high-resolution upwind finite-volumes schemes, multi-stage point implicit time stepping, and automatic adaptive mesh refinement (AMR) to solve the governing conservation equations for the moment closures. The initial phase of the code development was completed and a numerical investigation of the solutions of the 10-moment closure model for the simple two-dimensional test cases mentioned above was initiated. Predictions of the 10-moment model were compared to available theoretical solutions and the results of direct-simulation Monte Carlo (DSMC) calculations. The first results of this study were presented at a meeting last year.

Groth, Clinton P. T.↗
NASA NTRS1999

Physically-Derived Dynamical Cores in Atmospheric General Circulation Models

The algorithm chosen to represent the advection in atmospheric models is often used as the primary attribute to classify the model. Meteorological models are generally classified as spectral or grid point, with the term grid point implying discretization using finite differences. These traditional approaches have a number of shortcomings that render them non-physical. That is, they provide approximate solutions to the conservation equations that do not obey the fundamental laws of physics. The most commonly discussed shortcomings are overshoots and undershoots which manifest themselves most overtly in the constituent continuity equation. For this reason many climate models have special algorithms to model water vapor advection. This talk focuses on the development of an atmospheric general circulation model which uses a consistent physically-based advection algorithm in all aspects of the model formulation. The shallow-water model of Lin and Rood (QJRMS, 1997) is generalized to three dimensions and combined with the physics parameterizations of NCAR's Community Climate Model. The scientific motivation for the development is to increase the integrity of the underlying fluid dynamics so that the physics terms can be more effectively isolated, examined, and improved. The expected benefits of the new model are discussed and results from the initial integrations will be presented.

Rood, Richard B.↗