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.

SourceResearch connectionFrom yearThrough year

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

At least 217 records · Page 12

All sources NASA NTRS DOE OSTI

NASA NTRS1991

Numerical investigation of an internal layer in turbulent flow over a curved hill

A numerical investigation of incompressible and compressible turbulent flows over strongly curved surfaces is presented. The turbulent flow equations are solved by a pressure based Navier-Stokes equations solver. In the method, the conservation of mass equation is replaced by a pressure correction equation applicable for both compressible and incompressible flows. The turbulence is described by a multiple time scale turbulence model supplemented with a near-wall turbulence model. The numerical results show that the internal layer is a strong turbulence field which is developed beneath the external boundary layer and is located very close to the wall. The development of the internal layer is attributed to the enormous mean flow strain rate caused by the streamline curvature. The external boundary layer flow responds rather slowly to the streamline curvature. Thus, the turbulence field of the forward corner of the curved hill is characterized by two turbulence fields interacting with each other. The turbulence intensity of the internal layer is much stronger than that of the external boundary layer, so that the development of a new boundary layer in the downstream region of the curved hill depends mostly on the internal layer. These numerical results are in good agreement with the measured data, and show that the turbulence model can resolve the turbulence field subjected to the strong streamline curvature.

Kim, S.-W.↗

NASA NTRS1991

Comparison of SMAC, PISO, and iterative time-advancing schemes for unsteady flows

Calculations of unsteady flows using a simplified marker and cell (SMAC), a pressure implicit splitting of operators (PSIO), and an iterative time advancing scheme (ITA) are presented. A partial differential equation for incremental pressure is used in each time advancing scheme. Example flows considered are a polar cavity flow starting from rest and self-sustained oscillating flows over a circular and a square cylinder. For a large time step size, the SMAC and ITA are more strongly convergent and yield more accurate results than PSIO. The SMAC is the most efficient computationally. For a small time step size, the three time advancing schemes yield equally accurate Strouhal numbers. The capability of each time advancing scheme to accurately resolve unsteady flows is attributed to the use of new pressure correction algorithm that can strongly enforce the conservation of mass. The numerical results show that the low frequency of the vortex shedding is caused by the growth time of each vortex shed into the wake region.

Kim, Sang-Wook↗

NASA NTRS1991

A mathematical model of the maximum power density attainable in an alkaline hydrogen/oxygen fuel cell

A mathematical model of a hydrogen/oxygen alkaline fuel cell is presented that can be used to predict the polarization behavior under various power loads. The major limitations to achieving high power densities are indicated and methods to increase the maximum attainable power density are suggested. The alkaline fuel cell model describes the phenomena occurring in the solid, liquid, and gaseous phases of the anode, separator, and cathode regions based on porous electrode theory applied to three phases. Fundamental equations of chemical engineering that describe conservation of mass and charge, species transport, and kinetic phenomena are used to develop the model by treating all phases as a homogeneous continuum.

Kimble, Michael C.↗

NASA NTRS1990

Heating of minor ions by the coronal slow shock

The coronal slow shock has been predicted to exist embedded in large coronal holes at 4-10 solar radii. In this paper, a three-fluid model is used to study the jumps in minor ion properties across a slow shock such as the coronal slow shock. The jump conditions are formulated in the de Hoffmann-Teller frame of reference. The Rankine-Hugoniot solution determines the MHD flow and the magnetic field across the shocks. For each minor ion species, the fluid equations for the conservation of mass, momentum, and energy can be solved to determine the velocity and the temperature of the ions across the shock. Also obtained is a similarity solution for heavy ions. The results show that, on the downstream side of the slow shock, the ion temperatures are nearly proportional to the ion masses for He, O, Si, and Fe, in agreement with observed ion temperatures in the inner solar wind. This indicates that the possibly existing coronal slow shock can be responsible for the observed heating of minor ions in the solar wind.

Whang, Y. C.↗

NASA NTRS1991

A hybrid numerical method for three-dimensional spatially-developing free-shear flows

The present novel algorithm for 3D incompressible Navier-Stokes equations treats a domain which is (1) infinite in the vertical direction, using a mapped spectral method; (2) finite in the streamwise direction, using a classical Fourier method; and (3) homogeneous in the spanwise direction, using high-order compact finite differencing. A projection method is employed which ensures the exact conservation of mass, as well as the satisfaction of the boundary conditions at infinity. The novel aspects of these methods are noted, and the code they constitute is validated in light of several test cases. Results are presented for two- and three-dimensional mixing layers.

Buell, Jeffrey C.↗

NASA NTRS1991

FLIP MHD - A particle-in-cell method for magnetohydrodynamics

The fluid-implicit-particle, or 'FLIP' method presently extended to 2D and 3D MHD flow incorporates a Lagrangian field representation and yields a grid magnetic Reynolds number of up to 16 while preserving contact continuities that retain the Galilean invariance of the MHD flow equations. Analytical arguments and numerical examples demonstrate the conservation of mass, momentum, magnetic flux, and energy; 2D calculation results for the illustrative cases of contact discontinuity convection, Rayleigh-Taylor unstable flow.

Brackbill, J. U.↗

NASA NTRS1991

One-way nesting for a primitive equation ocean model

Prognostic numerical models for atmospheric and oceanic circulations require initial fields, boundary conditions, and forcing functions in addition to a consistent set of partial differential equations, including a state relation and equations expressing conservation of mass, momentum, and energy. Depending on the horizontal domain to be modeled, the horizontal boundary conditions are either physically obvious or extremely difficult to specify consistently. If the entire atmosphere is modeled, periodic horizontal boundary conditions are appropriate. On the other hand, the physical horizontal boundaries on the entire ocean are solid walls. Obviously, the normal velocity at a solid wall is zero while the specification of the tangential velocity depends on the mathematical treatment of the horizontal viscous terms. Limitations imposed by computer capacity and cost, as well as research interests, have led to the use of limited area models to study flows in the atmosphere and ocean. The limited area models do not have physical horizontal boundaries, merely numerical ones. Correctly determining these open boundary conditions for limited-area numerical models has both intrigued and frustrated numerical modelers for decades. One common approach is to use the closed or solid wall boundary conditions for a limited-area model. The argument given for this approach is that the boundary conditions affect flow near the walls but that none of these effects are propagated into the interior. Therefore, one chooses a big enough domain that the central region of interest is not corrupted by the boundary flow. Research in progress to model the North Atlantic circulation vividly illustrates the pitfalls of this approach. Two model runs are compared: (1) the southern boundary at 20S between latitudes 0 and 40W is artificially closed; and (2) the same boundary is specified as open with an inward transport of 15 Sv (determined from a global model with the same physics) uniformly spread across the boundary. A comparison of both runs is presented.

Blake, D. W.↗

NASA NTRS1991

A guide to the use of theoretical models of the solar nebula for the interpretation of the meteoritic record

Attempts to derive a theoretical framework for the interpretation of the meteoritic record have been frustrated by our incomplete understanding of the fundamental processes that controlled the evolution of the primitive solar nebula. Nevertheless, it is possible to develop qualitative models of the nebula that illuminate its dynamic character, as well as the roles of some key parameters. These models draw on the growing body of observational data on the properties of disks around young, solar-type stars, and are constructed by applying the results of known solutions of protostellar collapse problems; making simple assumptions about the radial variations of nebular variables; and imposing the integral constraints demanded by conservation of mass, angular momentum, and energy. The models so constructed are heuristic, rather than predictive; they are intended to help us think about the nebula in realistic ways, but they cannot provide a definitive description of conditions in the nebula.

Cassen, Pat↗

NASA NTRS1991

Ignition transient analysis of solid rocket motor

Measurement data on the performance of Space Shuttle Solid Rocket Motor show wide variations in the head-end pressure changes and the total thrust build-up during the ignition transient periods. To analyze the flow and thermal behavior in the tested solid rocket motors, a 1-dimensional, ideal gas flow model via the SIMPLE algorithm was developed. Numerical results showed that burning patterns in the star-shaped head-end segment of the propellant and the erosive burning rate are two important factors controlling the ignition transients. The objective of this study is to extend the model to include the effects of aluminum particle commonly used in solid propellants. To treat the effects of aluminum-oxide particles in the combustion gas, conservation of mass, momentum, and energy equations for the particles are added in the numerical formulation and integrated by an inter-phase-slip algorithm.

Han, Samuel S.↗

NASA NTRS1991

Elliptic flow computation by low Reynolds number two-equation turbulence models

A detailed comparison of ten low-Reynolds-number k-epsilon models is carried out. The flow solver, based on an implicit approximate factorization method, is designed for incompressible, steady two-dimensional flows. The conservation of mass is enforced by the artificial compressibility approach and the computational domain is discretized using centered finite differences. The turbulence model predictions of the flow past a hill are compared with experiments at Re = 10 exp 6. The effects of the grid spacing together with the numerical efficiency of the various formulations are investigated. The results show that the models provide a satisfactory prediction of the flow field in the presence of a favorable pressure gradient, while the accuracy rapidly deteriorates when a strong adverse pressure gradient is encountered. A newly proposed model form that does not explicitly depend on the wall distance seems promising for application to complex geometries.

Michelassi, V.↗

NASA NTRS1991

A two-dimensional Euler solution for an unbladed jet engine configuration

A two dimensional nonaxisymmetric Euler solution in a geometric representation of a jet engine configuration without blades is presented. The domain, including internal and external flow, is covered with a multiblock grid. To construct the grid, a domain decomposition technique is used to subdivide the domain and smooth grids are dimensioned and placed in each block. The grid contains 44 blocks which cover the external field, the inlet, bypass duct, core duct and nozzle of the nonaxisymmetric engine configuration. The geometry is symmetric about the meanline of the hub, but the grid is not since there is no symmetry condition applied to the grid between the two halves. With a symmetric grid at zero angle of attack, the measures of the solution would cancel exactly. With an asymmetric grid, the solution will not necessarily be symmetric and the lift coefficient will not necessarily be zero. Thus, grid asymmetry can be exploited to verify the resolution of the solution. The solution may be verified on the basis of five theoretical quantities: conservation of mass and energy, deviation of the lift coefficient from zero, deviation of the drag coefficient from zero, deviations from constant entropy, and deviations in the pressure distributions over the symmetric surfaces of the components. This technique is suitable for obtaining numerical solutions in complex geometries and provides a foundation for complete engine throughflow calculations.

Stewart, Mark E. M.↗

NASA NTRS1992

Development of an improved model for runback water on aircraft surfaces

A computer simulation for 'running wet' and evaporative aircraft anti-icing systems is developed. The model is based on the analysis of the liquid water film which forms in the regions of direct impingement and, then, breaks up near the impingement limits into rivulets. The wetness factor distribution resulting from the film breakup and the rivulet configuration on the surface are predicted using a stability analysis theory and the laws of mass energy conservation. The solid structure is modeled as a multiple layer wall. The anti-icing system modeled is of the thermal type utilizing hot air and/or electrical heating elements embedded within the wall layers. Experimental observations revealing some of the basic physics of the water flow on the surface are presented. Detailed qualitative documentation of the tests are given. Several numerical examples are considered, and the effect of some of the involved parameters on the system performance are investigated.

Al-Khalil, Kamel M.↗

NASA NTRS1992

On-orbit cryogenic fluid transfer research at NASA Lewis Research Center

A summary of research into on-orbit cryogenic fluid transfer at the NASA Lewis Research Center (LRC) is presented. Variable test parameters and liquid injection configurations elucidated the conditions necessary for a successful transfer of liquid hydrogen by the no-vent fill method. The model is based on conservation of mass and a first-law energy balance for a control volume. The ullage, the bulk liquid, and the tank wall are each represented by a single node. The magnitude of the maximum receiver tank pressure was found to be dependent on the liquid inlet temperature, the inlet mass flow rate, and the initial temperature. A comparison of the test data and the analytical results for a no-vent fill test with a small receiver tank is presented.

Taylor, W. J.↗

NASA NTRS1991

Sensing and actuation in a liquid propellant rocket engine

Equations describing the conservation of mass, momentum, and energy, and multicomponent transport and chemical reaction in a liquid propellant rocket (LPR) engine are given, along with appropriate boundary conditions. The formulation includes a detailed description of the two-phase (gas/liquid) flow in terms of a moving boundary problem, with droplets and vapor occupying time-dependent subdomains of the LPR chamber. Several sensor and actuator types are considered in detail, and their mathematical descriptions are incorporated into the differential equations and boundary conditions.

Bentsman, Joseph↗

NASA NTRS1992

Comparison of the SMAC, PISO, and iterative time-advancing schemes for unsteady flows

Calculations of unsteady flows using a simplified marker and cell (SMAC), a pressure implicit splitting of operators (PISO), and an iterative time advancing scheme (ITA) are presented. A partial differential equation for incremental pressure is used in each time advancing scheme. Example flows considered are a polar cavity flow starting from rest and self-sustained oscillating flows over a circular and a square cylinder. For a large time step size, the SMAC and ITA are more strongly convergent and yield more accurate results than PISO. The SMAC is the most efficient computationally. For a small time step size, the three time advancing schemes yield equally accurate Strouhal numbers. The capability of each time advancing scheme to accurately resolve unsteady flows is attributed to the use of new pressure correction algorithm that can strongly enforce the conservation of mass. The numerical results show that the low frequency of the vortex shedding is caused by the growth time of each vortex shed into the wake region.

Kim, S.-W.↗

NASA NTRS1993

A Fractional-Step Method Of Computing Incompressible Flow

Method of computing time-dependent flow of incompressible, viscous fluid involves numerical solution of Navier-Stokes equations on two- or three-dimensional computational grid based on generalized curvilinear coordinates. Equations of method derived in primitive-variable formulation. Dependent variables are pressure at center of each cell of computational grid and volume fluxes across faces of each cell. Volume fluxes replace Cartesian components of velocity; these fluxes correspond to contravariant components of velocity multiplied by volume of computational cell, in staggered grid. Choice of dependent variables enables simple extension of previously developed staggered-grid approach to generalized curvilinear coordinates and facilitates enforcement of conservation of mass.

Kwak, Dochan↗

NASA NTRS1993

Light ion flow in the nightside ionosphere of Venus

The flow characteristics of the light ions H(+) and He(+) have been studied in the midnight region of the ionosphere of Venus. Measurements of ion composition, electron and ion temperatures and magnetic fields by instruments onboard the Pioneer Venus Orbiter have been used in rite electron and ion equations of conservation of mass and momentum to derive the vertical flow velocities of H(+) and He(+). When average height profiles of the measured quantities were used, H(+) was found to flow upward, accelerating to speeds of almost 1 km/s at the ion-exobase. In a similar fashion, He(+) was found to flow downward into the neutral atmosphere where it is readily quenched by charge transfer reactions. The polarization electric field played an important role in forcing H(+) upward, but did not contribute enough to the He(+) force balance to produce upward flow. At the ion-exobase, the outward electric polarization force on H(+) was shown to be five times the gravitational force. Using an analogy with the terrestrial ion-exosphere, H(+) was inferred to flow upward into the ionotail of Venus and accelerate to escape speeds. A planet averaged escape flux of 1.4 x 10 exp 7/sq cm/s was calculated, which is comparable to hydrogen loss rates estimated by other investigators.

Hartle, R. E.↗

NASA NTRS1993

An Eulerian/Lagrangian coupling procedure for three-dimensional vortical flows

A coupled Eulerian/Lagrangian method is presented for the reduction of numerical diffusion observed in solutions of 3D vortical flows using standard Eulerian finite-volume time-marching procedures. A Lagrangian particle tracking method, added to the Eulerian time-marching procedure, provides a correction of the Eulerian solution. In turn, the Eulerian solution is used to integrate the Lagrangian state-vector along the particles trajectories. While the Eulerian solution ensures the conservation of mass and sets the pressure field, the particle markers describe accurately the convection properties and enhance the vorticity and entropy capturing capabilities of the Eulerian solver. The Eulerian/Lagrangian coupling strategies are discussed and the combined scheme is tested on a constant stagnation pressure flow in a 90 deg bend and on a swirling pipe flow. As the numerical diffusion is reduced when using the Lagrangian correction, a vorticity gradient augmentation is identified as a basic problem of this inviscid calculation.

Felici, Helene M.↗