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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 415 records · Page 23

All sources NASA NTRS DOE OSTI

NASA NTRS1993

On the structure of gaseous confined laminar diffusion flames: Numerical investigation

The structure and characteristics of gaseous confined laminar diffusion flames are investigated by numerically solving the time-dependent two-dimensional axisymmetric conservation equations. The numerical model accounts for the important chemical and physical processes involved, including axial diffusion, viscous effects, radial convection, and finite-rate chemistry. The numerical results clearly show that the flame has a finite thickness and leakage of fuel vapor into the flame zone is possible. The effect of heat release is found to induce some radial flow. Predicted flame shape and dimensions are compared to the classical Burke-Schumann flame. The numerically calculated flame is observed to be about 15 percent taller and 5 percent narrower than that of the Burke-Schumann solution under the same conditions.

Mawid, M. A.↗

NASA NTRS1992

A k-omega-multivariate beta PDF for supersonic combustion

In an attempt to study the interaction between combustion and turbulence in supersonic flows, an assumed PDF has been employed. This makes it possible to calculate the time average of the chemical source terms that appear in the species conservation equations. In order to determine the averages indicated in an equation, two transport equations, one for the temperature (enthalpy) variance and one for Q, are required. Model equations are formulated for such quantities. The turbulent time scale controls the evolution. An algebraic model similar to that used by Eklund et al was used in an attempt to predict the recent measurements of Cheng et al. Predictions were satisfactory before ignition but were less satisfactory after ignition. One of the reasons for this behavior is the inadequacy of the algebraic turbulence model employed. Because of this, the objective of this work is to develop a k-omega model to remedy the situation.

Alexopoulos, G. A.↗

NASA NTRS1993

Structure of confined laminar spray diffusion flames: Numerical investigation

The structure of confined laminar spray diffusion flames is investigated numerically by solving the gas-phase conservation equations for mass species, continuity, momentum, and energy and the liquid-phase equations for droplet position, velocity, size, and temperature. A one-step global reaction scheme along with six equilibrium reactions are employed to model the flame chemistry. Monodisperse as well as polydisperse sprays are considered. The numerical results demonstrate that liquid spray flames substantially differ from gaseous flames in their structure, i.e., temperature, concentration, and velocity fields, shape, and dimensions under the same conditions. Spray flames are predicted to be taller and narrower than their counterpart gaseous ones and their shapes are almost cylindrical. This is in agreement with experimental observations. The numerical computations also show that the use of the equilibrium reactions with the one-step reaction scheme decreases the flame temperature compared to the one-step reaction scheme without the equilibrium reactions and more importantly increases the surface area of the flame zone due to a phenomenon termed 'equilibrium broadening.' The spray flames also possess a finite thickness with minimal overlap of the fuel and oxygen species. A case for which a fuel-mixture consisting of 20 to 80 percent gas-liquid by mass is introduced into the combustor is also investigated and compared with predictions using only gaseous or liquid fuel.

Mawid, M. A.↗

NASA NTRS1992

A second-order bulk boundary-layer model

Bulk mass-flux models represent the large eddies that are primarily responsible for the turbulent fluxes in the planetary boundary layer as convective circulations, with an associated convective mass flux. In order for such models to be useful, it is necessary to determine the fractional area covered by rising motion in the convective circulations. This fraction can be used as an estimate of the cloud amount, under certain conditions. 'Matching' conditions have been developed that relate the convective mass flux to the ventilation and entrainment mass fluxes. These are based on conservation equations for the scalar means and variances in the entrainment and ventilation layers. Methods are presented to determine both the fractional area covered by rising motion and the convective mass flux. The requirement of variance balance is used to relax the 'well-mixed' assumption. The vertical structures of the mean state and the turbulent fluxes are determined analytically. Several aspects of this simple model's formulation are evaluated using results from large-eddy simulations.

Randall, David A.↗

NASA NTRS1992

Three-dimensional analytical model for the mixed layer depth

A three-dimensional (3D) semianalytical mixed-layer thermal inertia model for the mixed layer depth is developed, which includes the horizontal advection and diffusion processes modeled previously by Yan et al. (1991) in a one-dimensional (1D) model. The 3D model, which is based on the 3D thermal energy conservation equation, can be used as a supplement to the 1D model in the area where advection is important and for a finer resolution grid prediction. The model can be easily forced by remotely sensed data and is much simpler to compute than numerical 3D mixed-layer turbulence closure models. The new model was verified using data from the AVHRR/multichannel sea surface temperature data set and the Cooperative Ocean Atmosphere Data Set for the North Pacific Ocean, showing good agreement.

Yan, Xiao-Hai↗

NASA NTRS1993

Diffusion effects on the helium abundance of the solar transition region and corona

The diffusion of helium in the solar transition region is studied by solving the mass and momentum conservation equations for a hydrogen-helium plasma given a representative temperature profile. Steady state solutions show that two distinct atmospheres may result. In cases where the thermal force on alpha-particles is balanced by the partial pressure gradient force, helium is the dominant coronal species. On the other hand, if it is the frictional force between protons and alpha-particles which balances the thermal force on alpha-particles then hydrogen is the major coronal component. In order to explore which of these solutions are attainable within reasonable time scales, the time-dependent equations are solved, starting from an initial state with a uniform helium abundance of 10 percent. The atmosphere as a whole is close to hydrostatic equilibrium, but due the thermal forces the individual elements are not. This force inbalance leads to a differential flow between species. It is found that this differential flow leads to a significant enhancement of the coronal helium abundance. Even for the relatively shallow temperature gradient used the helium abundance in the lower corona increases to 30 percent over a 24 hr period.

Hansteen, Viggo H.↗

NASA NTRS1993

De-Dopplerization of aircraft acoustic signals

A de-Dopplerization scheme is devised and applied to both a tonal noise source, a XV-15 aircraft, and a broadband noise source, a F-18 aircraft. The procedure is developed from the conservation equations of fluid mechanics and is described in the paper. The corrected time history is constructed by using linear interpolation in the measured time history. This was made possible by knowing the position history of the aircraft provided by radar and/or laser tracking. The XV-15 data established that the scheme can accurately account for Doppler frequency shifts. The F-18 data confirms what has been noted in prediction models and static tests pertaining to broadband shock associated noise. That is, the peak frequency increases and the peaks broaden toward the jet axis. Another issue addressed in the study is the influence of correcting for spherical spreading and Doppler amplitude on the spectral shape and overall sound pressure levels of the source. Results from this investigation confirm that the dominant noise source in high-speed jets is due to turbulent mixing.

Kelly, Jeffrey J.↗

NASA NTRS1993

Heat Transfer from Radiatively Heated Material in a Low Reynolds Number Microgravity Environment

A mathematical model of the transient three-dimensional heat transfer between a slowly moving ambient gas stream and a thermally thick or thin flat surface heated by external radiation in a microgravity environment is presented. The problem is motivated in part by fire safety issues in spacecraft. The gas phase is represented by variable property convection-diffusion energy and mass conservation equations valid at low Reynolds numbers. The absence of gravity and low Reynolds number together permit the flow to be represented by a self-consistent velocity potential determined by the ambient velocity and the thermal expansion in the gas. The solid exchanges energy with the gas by conduction/convection and with the surroundings by surface absorption and re-emission of radiation. Heat conduction in the solid is assumed to be one dimensional at each point on the surface as a consequence of the limited times (of order of 10 seconds) of interest in these simulations. Despite the apparent simplicity of the model, the results show a complex thermally induced flow near the heated surface. The thermal exchange between the gas and solid produces an outward sourcelike flow upstream of the center of the irradiated area and a sinklike flow downstream. The responses of the temperature fields and the associated flows to changes in the intensity of the external radiation and the ambient velocity are discussed.

Yamashita, H.↗

NASA NTRS1993

Evaporation of LOX under supercritical and subcritical conditions

The evaporation of LOX under supercritical and subcritical conditions was studied experimentally and theoretically. In experiments, the evaporation rate and surface temperature were measured for LOX strand vaporizing in helium environments at pressures ranging from 5 to 68 atmospheres. Gas sampling and chromatography analysis were also employed to profile the gas composition above the LOX surface for the purpose of model validation. A comprehensive theoretical model was formulated and solved numerically to simulate the evaporation process of LOX at high pressures. The model was based on the conservation equations of mass, momentum, energy, and species concentrations for a multicomponent system, with consideration of gravitational body force, solubility of ambient gases in liquid, and variable thermophysical properties. Good agreement between predictions and measured oxygen mole fraction profiles was obtained. The effect of pressure on the distribution of the Lewis number, as well as the effect of variable diffusion coefficient, were further examined to elucidate the high-pressure transport behavior exhibited in the LOX vaporization process.

Yang, A. S.↗

NASA NTRS1990

Diabatic forcing and initialization with assimilation of cloud and rain water in a forecast model: Methodology

The focus of this part of the investigation is to find one or more general modeling techniques that will help reduce the time taken by numerical forecast models to initiate or spin-up precipitation processes and enhance storm intensity. If the conventional data base could explain the atmospheric mesoscale flow in detail, then much of our problem would be eliminated. But the data base is primarily synoptic scale, requiring that a solution must be sought either in nonconventional data, in methods to initialize mesoscale circulations, or in ways of retaining between forecasts the model generated mesoscale dynamics and precipitation fields. All three methods are investigated. The initialization and assimilation of explicit cloud and rainwater quantities computed from conservation equations in a mesoscale regional model are examined. The physical processes include condensation, evaporation, autoconversion, accretion, and the removal of rainwater by fallout. The question of how to initialize the explicit liquid water calculations in numerical models and how to retain information about precipitation processes during the 4-D assimilation cycle are important issues that are addressed. The explicit cloud calculations were purposely kept simple so that different initialization techniques can be easily and economically tested. Precipitation spin-up processes associated with three different types of weather phenomena are examined. Our findings show that diabatic initialization, or diabatic initialization in combination with a new diabatic forcing procedure, work effectively to enhance the spin-up of precipitation in a mesoscale numerical weather prediction forecast. Also, the retention of cloud and rain water during the analysis phase of the 4-D data assimilation procedure is shown to be valuable. Without detailed observations, the vertical placement of the diabatic heating remains a critical problem.

Raymond, William H.↗

NASA NTRS1993

Simplified, inverse, ejector design tool

A simple lumped parameter based inverse design tool has been developed which provides flow path geometry and entrainment estimates subject to operational, acoustic, and design constraints. These constraints are manifested through specification of primary mass flow rate or ejector thrust, fully-mixed exit velocity, and static pressure matching. Fundamentally, integral forms of the conservation equations coupled with the specified design constraints are combined to yield an easily invertible linear system in terms of the flow path cross-sectional areas. Entrainment is computed by back substitution. Initial comparison with experimental and analogous one-dimensional methods show good agreement. Thus, this simple inverse design code provides an analytically based, preliminary design tool with direct application to High Speed Civil Transport (HSCT) design studies.

Dechant, Lawrence J.↗

NASA NTRS1992

Modeling of hypersonic reacting flows

The effect of nonequilibrium thermochemical processes on the aerodynamic and heat transfer characteristics of a hypersonic vehicle is illustrated with examples. It is shown that the conventional method of predicting chemical reactions always predict that the flow is closer to equilibrium than it actually is, leading to incorrect predictions of aerodynamic characteristics of a vehicle. A method is presented for predicting chemical processes using a two-temperature kinetic model, which is developed on the basis of combining conservation equations for vibrational energy and for electron electronic energy, which are also derived. The validity of the two-temperature model is proven using two types of experimental data: radiation and shock shapes.

Park, Chul↗

NASA NTRS1994

Lox droplet vaporization in a supercritical forced convective environment

A systematic investigation has been conducted to study the effects of ambient flow conditions (i.e. pressure and velocity) on supercritical droplet gasification in a forced-convective environment. The model is based on the time-dependent conservation equations in axisymmetric coordinates, and accommodates thermodynamic nonidealities and transport anomalies. In addition, an efficient scheme for evaluating thermophysical properties over the entire range of fluid thermodynamic states is established. The analysis allows a thorough examination of droplet behavior during its entire lifetime, including transient gasification, dynamic deformation, and shattering. A parametric study of droplet vaporization rate in terms of ambient pressure and Reynolds number is also conducted.

Hsiao, Chia-Chun↗

NASA NTRS1994

A combined Eulerian-volume of fraction-Lagrangian method for atomization simulation

The tracking of free surfaces between liquid and gas phases and analysis of the interfacial phenomena between the two during the atomization and breakup process of a liquid fuel jet is modeled. Numerical modeling of liquid-jet atomization requires the resolution of different conservation equations. Detailed formulation and validation are presented for the confined dam broken problem, the water surface problem, the single droplet problem, a jet breakup problem, and the liquid column instability problem.

Seung, S. P.↗

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 NTRS1994

Effects of non-unity Lewis numbers in diffusion flames

The purpose of this work is to carry out direct numerical simulations of diffusion controlled combustion with non-unity Lewis numbers for the reactants and products, thus accounting for the differential diffusion effects of the temperature and concentration fields. We use a formulation based on combining the conservation equations in a way to eliminate the reaction terms similar to the method used by Burke and Schumann (1928) for unity Lewis numbers. We present calculations for an axisymmetric fuel jet and for a planar, time evolving mixing layer, leaving out the effects of thermal expansion and variations of the transport coefficients due to the heat release. Our results show that the front of the flame shifts toward the fuel or oxygen sides owing to the effect of the differential diffusion and that the location of maximum temperature may not coincide with the flame. The dependence of the distribution of the reaction products on their Lewis number has been investigated.

Linan, A.↗

NASA NTRS1993

A Generalized Eulerian-Lagrangian Analysis, with Application to Liquid Flows with Vapor Bubbles

Under a NASA MSFC SBIR Phase 2 effort an analysis has been developed for liquid flows with vapor bubbles such as those in liquid rocket engine components. The analysis is based on a combined Eulerian-Lagrangian technique, in which Eulerian conservation equations are solved for the liquid phase, while Lagrangian equations of motion are integrated in computational coordinates for the vapor phase. The novel aspect of the Lagrangian analysis developed under this effort is that it combines features of the so-called particle distribution approach with those of the so-called particle trajectory approach and can, in fact, be considered as a generalization of both of those traditional methods. The result of this generalization is a reduction in CPU time and memory requirements. Particle time step (stability) limitations have been eliminated by semi-implicit integration of the particle equations of motion (and, for certain applications, the particle temperature equation), although practical limitations remain in effect for reasons of accuracy. The analysis has been applied to the simulation of cavitating flow through a single-bladed section of a labyrinth seal. Models for the simulation of bubble formation and growth have been included, as well as models for bubble drag and heat transfer. The results indicate that bubble formation is more or less 'explosive'. for a given flow field, the number density of bubble nucleation sites is very sensitive to the vapor properties and the surface tension. The bubble motion, on the other hand, is much less sensitive to the properties, but is affected strongly by the local pressure gradients in the flow field. In situations where either the material properties or the flow field are not known with sufficient accuracy, parametric studies can be carried out rapidly to assess the effect of the important variables. Future work will include application of the analysis to cavitation in inducer flow fields.

Dejong, Frederik J.↗

NASA NTRS1994

Evolution of solitary density waves in stellar winds of early-type stars: A simple explanation of discrete absorption component behavior

We model the evolution of a density shell propagating through the stellar wind of an early-type star, in order to investigate the effects of such shells on UV P Cygni line profiles. Unlike previous treatments, we solve the mass, momentum, and energy conservation equations, using an explicit time-differencing scheme, and present a parametric study of the density, velocity, and temperature response. Under the assumed conditions, relatively large spatial scale, large-amplitude density shells propagate as stable waves through the supersonic portion of the wind. Their dynamical behavior appears to mimic propagating 'solitary waves,' and they are found to accelerate at the same rate as the underlying steady state stellar wind (i.e., the shell rides the wind). These hydrodynamically stable structures quantitatively reproduce the anomalous 'discrete absorption component' (DAC) behavior observed in the winds of luminous early-type stars, as illustrated by comparisons of model predictions to an extensive International Ultraviolet Explorer (IUE) time series of spectra of zeta Puppis (O4f). From these comparisons, we find no conclusive evidence indicative of DACs accelerating at a significantly slower rate than the underlying stellar wind, contrary to earlier reports. In addition, these density shells are found to be consistent within the constraints set by the IR observations. We conclude that the concept of propagating density shells should be seriously reconsidered as a possible explanation of the DAC phenomenon in early-type stars.

Waldron, Wayne L.↗