Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “mass conservation”

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 487 records · Page 27
All sourcesNASA NTRSDOE OSTI
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

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 NTRS1994

Direct numerical simulation of sheared turbulent flow

The summer assignment to study sheared turbulent flow was divided into three phases which were: (1) literature survey, (2) computational familiarization, and (3) pilot computational studies. The governing equations of fluid dynamics or Navier-Stokes equations describe the velocity, pressure, and density as functions of position and time. In principle, when combined with conservation equations for mass, energy, and thermodynamic state of the fluid a determinate system could be obtained. In practice the Navier-Stokes equations have not been solved due to the nonlinear nature and complexity of these equations. Consequently, the importance of experiments in gaining insight for understanding the physics of the problem has been an ongoing process. Reasonable computer simulations of the problem have occured as the computational speed and storage of computers has evolved. The importance of the microstructure of the turbulence dictates the need for high resolution grids in extracting solutions which contain the physical mechanisms which are essential to a successful simulation. The recognized breakthrough occurred as a result of the pioneering work of Orzag and Patterson in which the Navier-Stokes equations were solved numerically utilizing a time saving toggling technique between physical and wave space, known as a spectral method. An equally analytically unsolvable problem, containing the same quasi-chaotic nature as turbulence, is known as the three body problem which was studied computationally as a first step this summer. This study was followed by computations of a two dimensional (2D) free shear layer.

Harris, Vascar G.↗
NASA NTRS1996

Thermohydrodynamic Analysis of Cryogenic Liquid Turbulent Flow Fluid Film Bearings

Computational programs developed for the thermal analysis of tilting and flexure-pad hybrid bearings, and the unsteady flow and transient response of a point mass rotor supported on fluid film bearings are described. The motion of a cryogenic liquid on the thin film annular region of a fluid film bearing is described by a set of mass and momentum conservation, and energy transport equations for the turbulent bulk-flow velocities and pressure, and accompanied by thermophysical state equations for evaluation of the fluid material properties. Zeroth-order equations describe the fluid flow field for a journal static equilibrium position, while first-order (linear) equations govern the fluid flow for small amplitude-journal center translational motions. Solution to the zeroth-order flow field equations provides the bearing flow rate, load capacity, drag torque and temperature rise. Solution to the first-order equations determines the rotordynamic force coefficients due to journal radial motions.

Luis San Andres↗
NASA NTRS1995

Eight-moment approximation solar wind models

Heat conduction from the corona is important in the solar wind energy budget. Until now all hydrodynamic solar wind models have been using the collisionally dominated gas approximation for the heat conductive flux. Observations of the solar wind show particle distribution functions which deviate significantly from a Maxwellian, and it is clear that the solar wind plasma is far from collisionally dominated. We have developed a numerical model for the solar wind which solves the full equation for the heat conductive flux together with the conservation equations for mass, momentum, and energy. The equations are obtained by taking moments of the Boltzmann equation, using an 8-moment approximation for the distribution function. For low-density solar winds the 8-moment approximation models give results which differ significantly from the results obtained in models assuming the gas to be collisionally dominated. The two models give more or less the same results in high density solar winds.

Olsen, Espen Lyngdal↗
NASA NTRS1996

Impeller leakage flow modeling for mechanical vibration control

HPOTP and HPFTP vibration test results have exhibited transient and steady characteristics which may be due to impeller leakage path (ILP) related forces. For example, an axial shift in the rotor could suddenly change the ILP clearances and lengths yielding dynamic coefficient and subsequent vibration changes. ILP models are more complicated than conventional-single component-annular seal models due to their radial flow component (coriolis and centrifugal acceleration), complex geometry (axial/radial clearance coupling), internal boundary (transition) flow conditions between mechanical components along the ILP and longer length, requiring moment as well as force coefficients. Flow coupling between mechanical components results from mass and energy conservation applied at their interfaces. Typical components along the ILP include an inlet seal, curved shroud, and an exit seal, which may be a stepped labyrinth type. Von Pragenau (MSFC) has modeled labyrinth seals as a series of plain annular seals for leakage and dynamic coefficient prediction. These multi-tooth components increase the total number of 'flow coupled' components in the ILP. Childs developed an analysis for an ILP consisting of a single, constant clearance shroud with an exit seal represented by a lumped flow-loss coefficient. This same geometry was later extended to include compressible flow. The objective of the current work is to: supply ILP leakage-force impedance-dynamic coefficient modeling software to MSFC engineers, base on incompressible/compressible bulk flow theory; design the software to model a generic geometry ILP described by a series of components lying along an arbitrarily directed path; validate the software by comparison to available test data, CFD and bulk models; and develop a hybrid CFD-bulk flow model of an ILP to improve modeling accuracy within practical run time constraints.

Palazzolo, Alan B.↗
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 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 NTRS2000

An Unsteady Long Bearing Squeeze Film Damper Model: Circular Centered Orbits - Part 1

This paper, the first of a two-part series, presents results of an unsteady rotordynamic analysis of a long-bearing squeeze film damper executing circular centered orbits using a fluid circuit approach. A series of nodes and branches represent the geometry of the flow circuit. The mass and momentum conservation equations are solved to predict the pressure distribution in the squeeze film. The motion of the bearing is simulated by the variation of geometry within the flow path. The modeling methodology is benchmarked against published experimental long-bearing squeeze film damper test results. The model provides good agreement with the experimental damping coefficient.

Schallhorn, P. A.↗
NASA NTRS2000

An Unsteady Long Bearing Squeeze Film Damper Model: Statically Eccentric Operation - Part 2

This paper, the second of a two-part series, presents results of an unsteady rotordynamic analysis of a long-bearing squeeze film damper executing orbits about an off center position using a fluid circuit approach. A series of nodes and branches represent the geometry of the flow circuit. The mass and momentum conservation equations are solved to predict the pressure distribution in the squeeze film. The motion of the bearing is simulated by the variation of geometry within the flow path. This effort represents the first modeling approach which allows for an arbitrary orbit size about an arbitrary position.

Schallhorn, P. A.↗
NASA NTRS2000

Oceanic Excitation of Polar Motion During 1992-1994

Two global ocean general circulation models have been used to compute the angular momentum of nontidal oceanic current and sea level height variations: (1) the Princeton Modular Ocean Model (MOM) having 22 vertical layers and a rigid lid, and (2) the Miami Isopycnic-Coordinate Ocean Model (MICOM) having 11 vertical layers with a mixed layer and a free surface. Both models were run on the same 2 degree longitude by 1 degree latitude grid spanning 80 S to 80 N latitude. Following a 10-year spin-up with climatological air-sea fluxes, both models were forced during 1992-1994 with daily wind and heat flux from the NCEP operational analysis and sea surface salinity restoring to Levitus climatology. After correcting for the effects of mass non-conservation in these models, they are shown to predict similar effects on polar motion excitation of sea level height variations, but the effect of currents predicted by MOM is much larger than that predicted by MICOM.

Gross, R.↗
NASA NTRS1995

Evolution of Vapor Bubbles Nucleation Sites in Low Gravity

When liquid is expelled by a vapor bubble growing at a nucleation site on a superheated surface, a thin microlayer underneath the bubble is left behind. It is evaporated from the free microlayer surface that provides for bubble growth. The average thickness of the microlayer determining the evaporation rate increases with time if the latter does not exceed a threshold value associated with the burn-out crisis. The bubble is described as a spherical segment with its flattened part adjoining the microlayer. This introduces two independent variables - the radius of the spherical part of the bubble surface and the polar angle that defines the relative area of the flattened part. They are to be found out from a set of two strongly nonlinear equations resulting from mass and momentum conservation laws. The first one depends on both microlayer thickness and nonmonotonously changing bubble base area. The second involves two major factors favoring bubble detachment - the buoyancy and a force due to the initial momentum of vapor input into the bubble. The former force depends on gravity whereas the latter one does not. It is why the limiting regimes of bubble evolution that correspond to normal or moderately reduced gravity and to microgravity feature drastically different properties. In the first case, the buoyancy dominates and the bubble evolves in such a manner as to become a full sphere at a moment that can be viewed as that of detachment. The detachment volume grows as gravity decreases. In the second case, the buoyancy is negligible and the bubble stays near the surface, while its volume continues to increase for a sufficiently long time. The findings are discussed in connection with experimental data obtained under different gravity conditions, some unpublished experiments being included. They help to understand why the pool boiling heat transfer coefficient frequently increases as gravity falls down and eventually vanishes.

Buyevich, Yu A.↗
NASA NTRS1998

Integration of Biological, Physical/Chemical and Energy Efficient Systems in the CELSS Antarctic Analog: Performance of Prototype Systems and Issues for Life Support

The Controlled Ecological Life Support System (CELSS) Antarctic Analog Project (CAAP) is a joint endeavor between the National Science Foundation, Office of Polar Programs (NSF-OPP), and the National Aeronautics and Space Administration (NASA). The fundamental objective is to develop, deploy, and operate a testbed of advanced life support technologies at the Amundsen-Scott South Pole Station that enable the objectives of both the NSF and NASA. The functions of food production, water purification, and waste treatment, recycle, and reduction provided by CAAP will improve the quality of life for the South Pole inhabitants, reduce logistics dependence, enhance safety, and minimize environmental impacts associated with human presence on the polar plateau. Because of the analogous technical, scientific, and mission features with Planetary missions, such as a mission to Mars, CAAP provides NASA with a method for validating technologies and overall approaches to supporting humans. Prototype systems for waste treatment, water recycle, resource recovery and crop production are being evaluated in a testbed at Ames Research Center. The combined performance of these biological and physical/chemical systems as an integrated function in support of the human habitat will be discussed. Overall system performance will be emphasized. The effectiveness and efficiency of component technologies will be discussed in the context of energy and mass flow within the system and contribution to achieving a mass and energy conservative system. Critical to the discussion are interfaces with habitat functions outside of the closed-loop life support: the ability of the system to satisfy the life support requirements of the habitat and the ability to define input requirements. The significance of analog functions in relation to future Mars habitats will be discussed.

Bubenheim, David L.↗
NASA NTRS2003

Regional-Scale Modeling at NASA Goddard Space Flight Center

Over the past decade, the Goddard Mesoscale Modeling and Dynamics Group has used a popular regional scale model, MM5, to study precipitation processes. Our group is making contributions to the MM5 by incorporating the following physical and numerical packages: improved Goddard cloud processes, a land processes model (Parameterization for Land-Atmosphere-Cloud Exchange - PLACE), efficient but sophisticated radiative processes, conservation of hydrometeor mass (water budget), four-dimensional data assimilation for rainfall, and better computational methods for trace gas transport. At NASA Goddard, the MM5 has been used to study: (1) the impact of initial conditions, assimilation of satellite-derived rainfall, and cumulus parameterizations on rapidly intensifying oceanic cyclones, hurricanes and typhoons, (2) the dynamic and thermodynamic processes associated with the development of narrow cold frontal rainbands, (3) regional climate and water cycles, (4) the impact of vertical transport by clouds and lightning on trace gas distributiodproduction associated with South and North American mesoscale convective systems, (5) the development of a westerly wind burst (WWB) that occurred during the TOGA COARE and the diurnal variation of precipitation in the tropics, (6) a Florida sea breeze convective event and a Mid-US flood event using a sophisticated land surface model, (7) the influence of soil heterogeneity on land surface energy balance in the southwest GCIP region, (8) explicit simulations (with 1.33 to 4 km horizontal resolution) of hurricanes Bob (1991) and Bonnie (1998), (9) a heavy precipitation event over Taiwan, and (10) to make real time forecasts for a major NASA field program. In this paper, the modifications and simulated cases will be described and discussed.

Tao, W.-K.↗
NASA NTRS2005

GRAVITATIONAL WAVE EXTRACTION FROM AN INSPIRALING CONFIGURATION OF MERGING BLACK HOLES

We present new techniques for evolving binary black hole systems which allow the accurate determination of gravitational waveforms directly from the wave zone region of the numerical simulations. Rather than excising the black hole interiors, our approach follows the "puncture" treatment of black holes, but utilizing a new gauge condition which allows the black holes to move successfully through the computational domain. We apply these techniques to an inspiraling binary, modeling the radiation generated during the final plunge and ringdown. We demonstrate convergence of the waveforms and and good conservation of mass-energy, with just over 3% of the system s mass converted to gravitational radiation.

Baker, John G.↗
NASA NTRS2002

Geostrophic balance with a full Coriolis Force: implications for low latitutde studies

In its standard form, geostrophic balance uses a partial representation of the Coriolis force. The resulting formation has a singularity at the equator, and violates mass and momentum conservation. When the horizontal projection of the planetary rotation vector is considered, the singularity at the equator disappears, continuity can be preserved, and quasigeostrophy can be formulated at planetary scale.

balance atmosphere geostrophy remote sensing wind↗
NASA NTRS2008

GFSSP Training Course Lectures

GFSSP has been extended to model conjugate heat transfer Fluid Solid Network Elements include: a) Fluid nodes and Flow Branches; b) Solid Nodes and Ambient Nodes; c) Conductors connecting Fluid-Solid, Solid-Solid and Solid-Ambient Nodes. Heat Conduction Equations are solved simultaneously with Fluid Conservation Equations for Mass, Momentum, Energy and Equation of State. The extended code was verified by comparing with analytical solution for simple conduction-convection problem The code was applied to model: a) Pressurization of Cryogenic Tank; b) Freezing and Thawing of Metal; c) Chilldown of Cryogenic Transfer Line; d) Boil-off from Cryogenic Tank.

Majumdar, Alok K.↗
NASA NTRS2010

Analytic Modeling of Pressurization and Cryogenic Propellant

An analytic model for pressurization and cryogenic propellant conditions during all mission phases of any liquid rocket based vehicle has been developed and validated. The model assumes the propellant tanks to be divided into five nodes and also implements an empirical correlation for liquid stratification if desired. The five nodes include a tank wall node exposed to ullage gas, an ullage gas node, a saturated propellant vapor node at the liquid-vapor interface, a liquid node, and a tank wall node exposed to liquid. The conservation equations of mass and energy are then applied across all the node boundaries and, with the use of perfect gas assumptions, explicit solutions for ullage and liquid conditions are derived. All fluid properties are updated real time using NIST Refprop.1 Further, mass transfer at the liquid-vapor interface is included in the form of evaporation, bulk boiling of liquid propellant, and condensation given the appropriate conditions for each. Model validation has proven highly successful against previous analytic models and various Saturn era test data and reasonably successful against more recent LH2 tank self pressurization ground test data. Finally, this model has been applied to numerous design iterations for the Altair Lunar Lander, Ares V Core Stage, and Ares V Earth Departure Stage in order to characterize Helium and autogenous pressurant requirements, propellant lost to evaporation and thermodynamic venting to maintain propellant conditions, and non-uniform tank draining in configurations utilizing multiple LH2 or LO2 propellant tanks. In conclusion, this model provides an accurate and efficient means of analyzing multiple design configurations for any cryogenic propellant tank in launch, low-acceleration coast, or in-space maneuvering and supplies the user with pressurization requirements, unusable propellants from evaporation and liquid stratification, and general ullage gas, liquid, and tank wall conditions as functions of time.

Corpening, Jeremy H.↗