SEARCH · Engineering Papers

Results for “flow 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 253 records · Page 14

All sources NASA NTRS DOE 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 NTRS2004

Routines for Computing Pressure Drops in Venturis

A set of computer-program routines has been developed for calculating pressure drops and recoveries of flows through standard venturis, nozzle venturis, and orifices. Relative to prior methods used for such calculations, the method implemented by these routines offers greater accuracy because it involves fewer simplifying assumptions and is more generally applicable to wide ranges of flow conditions. These routines are based on conservation of momentum and energy equations for real nonideal fluids, the properties of which are calculated by curve-fitting subroutines based on empirical properties data. These routines are capable of representing cavitating, choked, non-cavitating, and unchoked flow conditions for liquids, gases, and supercritical fluids. For a computation of flow through a given venturi, nozzle venturi, or orifice, the routines determine which flow condition occurs: First, they calculate a throat pressure under the assumption that the flow is unchoked or non-cavitating, then they calculate the throat pressure under the assumption that the flow is choked or cavitating. The assumption that yields the higher throat pressure is selected as the correct one.

de Quay, Laurence↗

NASA NTRS1979

Application of the principle of similarity fluid mechanics

Possible applications of the principle of similarity to fluid mechanics is described and illustrated. In correlating thermophysical properties of fluids, the similarity principle transcends the traditional corresponding states principle. In fluid mechanics the similarity principle is useful in correlating flow processes that can be modeled adequately with one independent variable (i.e., one-dimensional flows). In this paper we explore the concept of transforming the conservation equations by combining similarity principles for thermophysical properties with those for fluid flow. We illustrate the usefulness of the procedure by applying such a transformation to calculate two phase critical mass flow through a nozzle.

Hendricks, R. C.↗

NASA NTRS1991

Plasma physics of accreting neutron stars

Plasma concepts and phenomena that are needed to understand X- and gamma-ray sources are discussed. The capture of material from the wind or from the atmosphere or envelope of a binary companion star is described and the resulting types of accretion flows discussed. The reasons for the formation of a magnetosphere around the neutron star are explained. The qualitative features of the magnetospheres of accreting neutron stars are then described and compared with the qualitative features of the geomagnetosphere. The conditions for stable flow and for angular and linear momentum conservation are explained in the context of accretion by magnetic neutron stars and applied to obtain rough estimates of the scale of the magnetosphere. Accretion from Keplerian disks is then considered in some detail. The radial structure of geometrically thin disk flows, the interaction of disk flows with the neutron star magnetosphere, and models of steady accretion from Keplerian disks are described. Accretion torques and the resulting changes in the spin frequencies of rotating neutron stars are considered. The predicted behavior is then compared with observations of accretion-powered pulsars. Magnetospheric processes that may accelerate particles to very high energies, producing GeV and, perhaps, TeV gamma-rays are discussed. Finally, the mechanisms that decelerate and eventually stop accreting plasma at the surfaces of strongly magnetic neutron stars are described.

Ghosh, Pranab↗

NASA NTRS1985

Flow structure capturing on overset patched meshes

The present paper describes the application of a simple, robust stable implicit approach to solutions of the conservative equations of gasdynamics on either composite or overset meshes to capture flow structures such as shocks with high resolution in complex geometric domains. Without requiring special flux conservative operators, but rather, interpolating conservative variable data at mesh boundaries, the well posed upwind method provides accurate and oscillation free solutions, even when shocks cross the interior patch boundaries. In three problems with flow complexities that require grid refinement, the paper demonstrates the capability to conveniently carry out for gasdynamics the adaptive refined meshing strategy in overset patches proposed by Berger and Oliger (1984), and it extends this technique to rapidly convergent implicit methods for the Euler and Navier-Stokes equations. The numerical experiments show concretely, in a realistic aerodynamic problem, the savings in mesh points (about an order of magnitude here in two dimensions) for similar accuracy that flow structure aligned adaptive patched meshing affords compared to uniform grid refinement.

Venkatapathy, E.↗

NASA NTRS1974

Comparison of approximate and numerical analyses of nonlinear combustion instability

At the present time, there are three general analytical techniques available to study problems of unsteady motions in rocket motors: linear stability analysis; approximate nonlinear analysis, founded on examining the behavior of coupled normal modes; and numerical calculations based on the conservation equations for one-dimensional flows. The last two yield the linear results as a limit. It is the main purpose of this paper to check the accuracy of the approximate analysis against the numerical analysis for some special cases. The results provide some justification for using the approximate analysis to study three-dimensional problems.

Culick, F. E. C.↗

NASA NTRS1976

Conduction of thermal energy in the neighborhood of the earth's bow shock

The Rankine-Hugoniot equations for MHD shocks are generalized by the addition of a term to the energy conservation equation representing a nonzero heat flow in the plasma in the neighborhood of the shock. This generalization is found to be compatible with the assumption of infinite electrical conductivity. The effects of plasma waves in this treatment are of the order of the reciprocal Alfvenic Mach number squared and hence are neglected. The effect of alpha particles in the solar wind is discussed. Seven crossings of the earth's bow shock by Explorer 35 in lunar orbit are analyzed. Sufficient data are available so that the determination of a dimensionless parameter, psi, characterizing the heat-flow difference across the bow shock is possible. The values of psi indicate energy-flux densities due to heat flow which are a nonnegligible fraction of the total energy flux. Two possible interpretations of psi are discussed.

Hohlfeld, R. G.↗

NASA NTRS1980

Status of linear boundary-layer stability and the e to the nth method, with emphasis on swept-wing applications

The-state-of-the-art for the application of linear stability theory and the e to the nth power method for transition prediction and laminar flow control design are summarized, with analyses of previously published low disturbance, swept wing data presented. For any set of transition data with similar stream distrubance levels and spectra, the e to the nth power method for estimating the beginning of transition works reasonably well; however, the value of n can vary significantly, depending upon variations in disturbance field or receptivity. Where disturbance levels are high, the values of n are appreciably below the usual average value of 9 to 10 obtained for relatively low disturbance levels. It is recommended that the design of laminar flow control systems be based on conservative estimates of n and that, in considering the values of n obtained from different analytical approaches or investigations, the designer explore the various assumptions which entered into the analyses.

Hefner, J. N.↗

NASA NTRS1982

An implicit solution of the three-dimensional Navier-Stokes equations for an airfoil spanning a wind tunnel

An implicit finite-difference algorithm is developed for the numerical solution of the incompressible three dimensional Navier-Stokes equations in the non-conservative primitive-variable formulation. The flow field about an airfoil spanning a wind-tunnel is computed. The coordinate system is generated by an extension of the two dimensional body-fitted coordinate generation techniques of Thompson, as well as that of Sorenson, into three dimensions. Two dimensional grids are stacked along a spanwise coordinate defined by a simple analytical function. A Poisson pressure equation for advancing the pressure in time is arrived at by performing a divergence operation on the momentum equations. The pressure at each time-step is calculated on the assumption that continuity be unconditionally satisfied. An eddy viscosity coefficient, computed according to the algebraic turbulence formulation of Baldwin and Lomax, simulates the effects of turbulence.

Moitra, A.↗

NASA NTRS1982

Numerical study of terrain-induced mesoscale motions in a mixed layer

Numerical integrations using a potential enstrophy-conserving scheme are presented for the flow within a mixed layer over hilly terrain using the hydrostatic shallow-water equations with a quadratic drag law. The mesoscale area treated is 150 km on a side; cyclic lateral boundary conditions are used. It is found that for the idealized conditions treated (no surface heating, no entrainment and no pressure adjustments aloft), the topography quickly induces a steady state flow pattern by means of surface friction. Unsteadiness does not occur unless a surface-friction Reynolds number is greater than approximately 100. Effects of varying the Rossby number, Froude number and terrain-height parameter are examined.

Han, Y.-J.↗

NASA NTRS1989

Numerical simulation of axisymmetric turbulent flow in combustors and diffusers

Numerical studies of turbulent flow in an axisymmetric 45-deg-expansion combustor and bifurcated diffuser are presented. The Navier-Stokes equations incorporating a k-epsilon model were solved in a nonorthogonal curvilinear coordinate system. A zonal-grid method, where the flow field was divided into several subsections, was developed. This approach permitted different computational schemes to be used in the various zones. In addition, grid generation was made a more simple task. Boundary overlap and interpolating techniques were used, and an adjustment of the flow variables was required to assure conservation of mass flux. Three finite-differencing methods (hybrid, quadratic upwind, and skew upwind) were used to represent the convection terms. Results were compared with existing experimental data. In general, good agreement between predicted and measured values was obtained.

Yung, Chain-Nan↗

NASA NTRS1990

Flux Jacobian Matrices For Equilibrium Real Gases

Improved formulation includes generalized Roe average and extension to three dimensions. Flux Jacobian matrices derived for use in numerical solutions of conservation-law differential equations of inviscid flows of ideal gases extended to real gases. Real-gas formulation of these matrices retains simplifying assumptions of thermodynamic and chemical equilibrium, but adds effects of vibrational excitation, dissociation, and ionization of gas molecules via general equation of state.

Vinokur, Marcel↗

NASA NTRS1990

Evolution and interaction of two- and three-dimensional instability waves

The present work is concerned with the nonlinear interactions of high-amplitude two- and three-dimensional instability waves in jets. A jet is excited by multifrequency axisymmetric and helical modes. Formulations are presented here for the nonlinear development of both the amplitudes and phases of the instability waves. The development of these waves are also coupled with the development of the mean flow and the background turbulence. In formulating this model it is assumed that each flow component can be characterized by conservations equations supplemented by closure models.

Mankbadi, R. R.↗

NASA NTRS1989

Numerical simulation of nonlinear development of instability waves

The nonlinear interactions of high amplitude instability waves in turbulent jets are described. In plane shear layers Riley and Metcalf (1980) and Monkewitz (1987) have shown that these interactions are dependent, among other parameters, on the phase-difference between the two instability waves. Therefore, here researchers consider the nonlinear development of both the amplitudes and the phase of the instability waves. The development of these waves are also coupled with the development of the mean flow and the background turbulence. In formulating this model it is assumed that each of the flow components can be characterized by conservation equations supplemented by closure models. Results for the interactions between the two instability waves under high-amplitude forcing at fundamental and subharmonic frequencies are presented here. Qualitative agreements are found between the present predictions and available experimental data.

Mankbadi, Reda R.↗

NASA NTRS1992

A subgrid mixing model for LES of nonpremixed turbulent combustion

A new sub-grid mixing model for use in large eddy simulations of turbulent combustion is presented and applied to a hydrogen-air diffusion flame. The sub-grid model is based on Kerstein's Linear Eddy Model (1988). The model is first used to predict the mixing of a conserved scalar in a turbulent shear flow. The model correctly predicts the behavior of the pdf of the scalar field and displays a nonmarching peak at the preferred mixture fraction as the shear layer is traversed. It is then illustrated how a reduced chemical mechanism can be implemented within the linear eddy subgrid model formulation. The model is used to predict NO formation in hydrogen-air diffusion flame using a reduced chemical mechanism involving nine reactive scalars.

Mcmurtry, P. A.↗

NASA NTRS1993

A review of chemically reactive turbulent flow mixing mechanisms and a new design for a low NO(x) combustor

A review of chemically reactive flow analysis is presented. Key studies are reviewed. The approaches to modeling chemically reactive flows that are discussed are PDF, conserved PDF, conserved scalar, stochastic, and mixing length. Turbulent mixing mechanisms are also discussed and a model is formulated that could be used for optimizing non-premixed gas reacting systems. Consequently, a re-examination of the NO(x) reduction concepts for nonpremixed reactants are suggested together with two new designs for low NO(x) subsonic combustor research.

Mcbeath, Giorgio↗

NASA NTRS2002

Towards Resolving the Crab Sigma-Problem: A Linear Accelerator?

Using the exact solution of the axisymmetric pulsar magnetosphere derived in a previous publication and the conservation laws of the associated MHD flow, we show that the Lorentz factor of the outflowing plasma increases linearly with distance from the light cylinder. Therefore, the ratio of the Poynting to particle energy flux, generically referred to as sigma, decreases inversely proportional to distance, from a large value (typically approx. greater than 10(exp 4)) near the light cylinder to sigma approx. = 1 at a transition distance R(sub trans). Beyond this distance the inertial effects of the outflowing plasma become important and the magnetic field geometry must deviate from the almost monopolar form it attains between R(sub lc), and R(sub trans). We anticipate that this is achieved by collimation of the poloidal field lines toward the rotation axis, ensuring that the magnetic field pressure in the equatorial region will fall-off faster than 1/R(sup 2) (R being the cylindrical radius). This leads both to a value sigma = a(sub s) much less than 1 at the nebular reverse shock at distance R(sub s) (R(sub s) much greater than R(sub trans)) and to a component of the flow perpendicular to the equatorial component, as required by observation. The presence of the strong shock at R = R(sub s) allows for the efficient conversion of kinetic energy into radiation. We speculate that the Crab pulsar is unique in requiring sigma(sub s) approx. = 3 x 10(exp -3) because of its small translational velocity, which allowed for the shock distance R(sub s) to grow to values much greater than R(sub trans).

Contopoulos, Ioannis↗

NASA NTRS2004

Variable Step Integration Coupled with the Method of Characteristics Solution for Water-Hammer Analysis, A Case Study

One-dimensional water-hammer modeling involves the solution of two coupled non-linear hyperbolic partial differential equations (PDEs). These equations result from applying the principles of conservation of mass and momentum to flow through a pipe, and usually the assumption that the speed at which pressure waves propagate through the pipe is constant. In order to solve these equations for the interested quantities (i.e. pressures and flow rates), they must first be converted to a system of ordinary differential equations (ODEs) by either approximating the spatial derivative terms with numerical techniques or using the Method of Characteristics (MOC). The MOC approach is ideal in that no numerical approximation errors are introduced in converting the original system of PDEs into an equivalent system of ODEs. Unfortunately this resulting system of ODEs is bound by a time step constraint so that when integrating the equations the solution can only be obtained at fixed time intervals. If the fluid system to be modeled also contains dynamic components (i.e. components that are best modeled by a system of ODEs), it may be necessary to take extremely small time steps during certain points of the model simulation in order to achieve stability and/or accuracy in the solution. Coupled together, the fixed time step constraint invoked by the MOC, and the occasional need for extremely small time steps in order to obtain stability and/or accuracy, can greatly increase simulation run times. As one solution to this problem, a method for combining variable step integration (VSI) algorithms with the MOC was developed for modeling water-hammer in systems with highly dynamic components. A case study is presented in which reverse flow through a dual-flapper check valve introduces a water-hammer event. The predicted pressure responses upstream of the check-valve are compared with test data.

Turpin, Jason B.↗