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 307 records · Page 17

All sources NASA NTRS DOE OSTI

NASA NTRS1985

SSME Turbopump Turbine Computations

A two-dimensional viscous code was developed to be used in the prediction of the flow in the SSME high-pressure turbopump blade passages. The rotor viscous code (RVC) employs a four-step Runge-Kutta scheme to solve the two-dimensional, thin-layer Navier-Stokes equations. The Baldwin-Lomax eddy-viscosity model is used for these turbulent flow calculations. A viable method was developed to use the relative exit conditions from an upstream blade row as the inlet conditions to the next blade row. The blade loading diagrams are compared with the meridional values obtained from an in-house quasithree-dimensional inviscid code. Periodic boundary conditions are imposed on a body-fitted C-grid computed by using the GRAPE GRids about Airfoils using Poisson's Equation (GRAPE) code. Total pressure, total temperature, and flow angle are specified at the inlet. The upstream-running Riemann invariant is extrapolated from the interior. Static pressure is specified at the exit such that mass flow is conserved from blade row to blade row, and the conservative variables are extrapolated from the interior. For viscous flows the noslip condition is imposed at the wall. The normal momentum equation gives the pressure at the wall. The density at the wall is obtained from the wall total temperature.

Jorgenson, P. G. E.↗

NASA NTRS1985

Time-dependent nonlinear cosmic ray shocks confirming abstract

Numerical studies of time dependent cosmic ray shock structures in planar geometry are interesting because analytical time-independent solutions are available which include the non-linear reactions on the plasma flow. A feature of these time asymptotic solutions is that for higher Mach numbers (M approximately 5) and for a low cosmic ray upstream pressure the solution is not uniquely determined by the usual conservation laws of mass, momentum and energy. These numerical solutions clearly indicate that much work needs to be done before we understand shock acceleration as a time dependent process. The slowness of the process is possibly due to the fact that there is a diffusive flux into the downstream region in addition to the usual advective losses. Analytic investigations of thie phenomenon are required.

Dorfi, E. A.↗

NASA NTRS1984

Numerical simulations of loops heated to solar flare temperatures. III - Asymmetrical heating

A numerical model is defined for asymmetric full solar flare loop heating and comparisons are made with observational data. The Dynamic Flux Tube Model is used to describe the heating process in terms of one-dimensional, two fluid conservation equations of mass, energy and momentum. An adaptive grid allows for the downward movement of the transition region caused by an advancing conduction front. A loop 20,000 km long is considered, along with a flare heating system and the hydrodynamic evolution of the loop. The model was applied to generating line profiles and spatial X-ray and UV line distributions, which were compared with SMM, P78-1 and Hintori data for Fe, Ca and Mg spectra. Little agreement was obtained, and it is suggested that flares be treated as multi-loop phenomena. Finally, it is concluded that chromospheric evaporation is not an effective mechanism for generating the soft X-ray bursts associated with flares.

Cheng, C.-C.↗

NASA NTRS1984

Dynamical relations for left ventricular ejection - Flow rate, momentum, force and impulse

An investigation was carried out to quantitatively evaluate left ventricular volume flow rate, momentum, force and impulse derived from application of conservation principles for mass and momentum of blood within the ventricle during the ejection phase. An automated digital image processing system was developed and applied to left ventricular angiograms which are computer processed and analyzed frame by frame to determine the dynamical relations by numerical methods. The initial experience with force and impulse has indicated that neither quantity seemed to be a sensitive indicator of coronary artery disease as evaluated by qualitative angiography for the particular patient group studied. Utilization of the dynamical relations in evaluating human left ventricular performance requires improved means of measurement and interpretation of clinical studies.

Back, L. H.↗

NASA NTRS1986

Transfer function analysis of thermospheric perturbations

Applying perturbation theory, a spectral model in terms of vectors spherical harmonics (Legendre polynomials) is used to describe the short term thermospheric perturbations originating in the auroral regions. The source may be Joule heating, particle precipitation or ExB ion drift-momentum coupling. A multiconstituent atmosphere is considered, allowing for the collisional momentum exchange between species including Ar, O2, N2, O, He and H. The coupled equations of energy, mass and momentum conservation are solved simultaneously for the major species N2 and O. Applying homogeneous boundary conditions, the integration is carred out from the Earth's surface up to 700 km. In the analysis, the spherical harmonics are treated as eigenfunctions, assuming that the Earth's rotation (and prevailing circulation) do not significantly affect perturbations with periods which are typically much less than one day. Under these simplifying assumptions, and given a particular source distribution in the vertical, a two dimensional transfer function is constructed to describe the three dimensional response of the atmosphere. In the order of increasing horizontal wave numbers (order of polynomials), this transfer function reveals five components. To compile the transfer function, the numerical computations are very time consuming (about 100 hours on a VAX for one particular vertical source distribution). However, given the transfer function, the atmospheric response in space and time (using Fourier integral representation) can be constructed with a few seconds of a central processing unit. This model is applied in a case study of wind and temperature measurements on the Dynamics Explorer B, which show features characteristic of a ringlike excitation source in the auroral oval. The data can be interpreted as gravity waves which are focused (and amplified) in the polar region and then are reflected to propagate toward lower latitudes.

Mayr, H. G.↗

NASA NTRS1985

Finite element algorithms for compressible flow computation on a supercomputer

Methods of applying computational fluid dynamics techniques to engineering problems are discussed. The linking factors between theoretical and applied research are the shape and weight functions which produce the Galerkin finite element schemes. The flow is unsteady, viscous, heat conducting and compressible, and steady-state flows are the asymptotic limit of unsteady flows. All flows are described by the time-averaged Navier-Stokes equations (NSE) with mass and energy conservation. Local curvilinear intrinsic coordinates are applied to discretize the NSE in arbitrary geometric domains. A modified weighted residuals approach defines the discrete analogs of the physical systems modeled and a general interpolants method is used to derive families of numerical models, both implicit and explicit, finite difference and finite elemnt, from a single point of departure. Time-dependent element approximations are achieved with a progressive assembly of generalized elements method. Hyperbolic steady-state Euler algorithms and quasi-parabolic are spatial marching algorithms for solving the discretized equations. Application of the techniques is illustrated through calculation of the compressible flow around a hypersonic flight vehicle at 100,000 ft at zero angle of attack at speeds of Mach 4-10.

Spradley, L. W.↗

NASA NTRS1986

Fluid flow and fuel-air mixing in a motored two-dimensional Wankel rotary engine

The implicit-factored method of Beam and Warming was employed to obtain numerical solutions to the conservation equations of mass, species, momentum, and energy to study the unsteady, multidimensional flow and mixing of fuel and air inside the combustion chambers of a two-dimensional Wankel rotary engine under motored conditions. The effects of the following engine design and operating parameters on fluid flow and fuel-air mixing during the intake and compression cycles were studied: engine speed, angle of gaseous fuel injection during compression cycle, and speed of the fuel leaving fuel injector.

Shih, T. I.-P.↗

NASA NTRS1986

Initiation of breakdown in slender compressible vortices

The onset of vortex breakdown in compressible flows is investigated analytically for the case in which the flow is axially symmetric, the vortex is isolated, its axis is parallel to the main flow, and the vortex radius is small compared to the breakdown length. The conservation equations for mass, momentum, and energy are formulated and solved numerically using a finite-difference scheme, as described by Krause (1985); numerical results are presented in graphs and briefly characterized.

Krause, E.↗

NASA NTRS1986

Initiation of breakdown in slender compressible vortices

The initiation of the breakdown process for axially symmetric compressible flows is investigated using a numerical solution of the conservation equations for mass, momentum, and energy. The vortex is isolated, with its axis parallel to the direction of the main stream, and the core radius is small compared to the breakdown length. Computations for several flowfields indicate that the breakdown of the solution is shifted further downstream with increasing Mach number until breakdown is no longer observed. In the subsonic case, the influence of the initial temperature distribution on the breakdown length of the solution is more pronounced than in the supersonic case, with heating of the core enhancing breakdown, and cooling delaying it. The breakdown of the solution is seen to always occur for nonvanishing axial velocity components.

Krause, E.↗

NASA NTRS1988

The computation of hypersonic ionized flows in chemical and thermal nonequlibrium

A numerical method to compute a two-dimensional hypersonic flowfield that is ionized and in thermochemical nonequilibrium has been developed. Such a flowfield is described by coupled time-dependent partial differential equations for the conservation of species mass, mass-average momentum, vibrational energy of each diatomic species, electron energy, and total mass-averaged energy. The steady-state solution to these fully coupled equations is obtained using an implicit Gauss-Seidel line relaxation technique. The computed electron densities in the flowfield compare well with experimental results.

Maccormack, Robert W.↗

NASA NTRS1988

Numerical modeling of chemical vapor deposition (CVD) in a horizontal reactor

In the present numerical prediction of the deposition rate of silicon from silane in a CVD process, the conservation equations for mass, momentum, energy, and chemical species are solved on a staggered grid using the SIMPLE algorithm, while the rate of chemical reactions in the gas phase and on the susceptor surface is obtained from an Arrhenius rate equation. Predicted deposition rates as a function of position along the susceptor with and without the gas phase chemical reaction are compared with the available experimental and numerical data; agreement is excellent except at the leading edge of the susceptor, where the deposition rate is overpredicted.

Sheikholeslami, M. Z.↗

NASA NTRS1989

Rotordynamic coefficients for labyrinth seals calculated by means of a finite difference technique

The compressible, turbulent, time dependent and three dimensional flow in a labyrinth seal can be described by the Navier-Stokes equations in conjunction with a turbulence model. Additionally, equations for mass and energy conservation and an equation of state are required. To solve these equations, a perturbation analysis is performed yielding zeroth order equations for centric shaft position and first order equations describing the flow field for small motions around the seal center. For numerical solution a finite difference method is applied to the zeroth and first order equations resulting in leakage and dynamic seal coefficients respectively.

Nordmann, R.↗

NASA NTRS1990

Two-dimensional Euler and Navier-Stokes Time accurate simulations of fan rotor flows

Two numerical methods are presented which describe the unsteady flow field in the blade-to-blade plane of an axial fan rotor. These methods solve the compressible, time-dependent, Euler and the compressible, turbulent, time-dependent, Navier-Stokes conservation equations for mass, momentum, and energy. The Navier-Stokes equations are written in Favre-averaged form and are closed with an approximate two-equation turbulence model with low Reynolds number and compressibility effects included. The unsteady aerodynamic component is obtained by superposing inflow or outflow unsteadiness to the steady conditions through time-dependent boundary conditions. The integration in space is performed by using a finite volume scheme, and the integration in time is performed by using k-stage Runge-Kutta schemes, k = 2,5. The numerical integration algorithm allows the reduction of the computational cost of an unsteady simulation involving high frequency disturbances in both CPU time and memory requirements. Less than 200 sec of CPU time are required to advance the Euler equations in a computational grid made up of about 2000 grid during 10,000 time steps on a CRAY Y-MP computer, with a required memory of less than 0.3 megawords.

Boretti, A. A.↗

NASA NTRS1990

Stability of stationary barotropic modons by Lyapunov's direct method

A new Liapunov stability condition is formulated for the shallow-water equations, using a gage-variable formalism. This sufficient condition is derived for the class of perturbations that conserve the total mass. It is weaker than existing stability criteria, i.e., it applies to a wider class of flows. Formal stability to infinitesimally small perturbations of arbitrary shape is obtained for two classes of large-scale geophysical flows: pseudo-eastward flow with constant shear, and localized coherent structures of modon type.

Sakuma, H.↗

NASA NTRS1990

New bound on neutrino dipole moments from globular-cluster stars

Neutrino dipole moments mu(nu) would increase the core mass of red giants at the helium flash by delta(Mc) = 0.015 solar mass x mu(nu)/10 to the -12th muB (where muB is the Bohr magneton) because of enhanced neutrino losses. Existing measurements of the bolometric magnitudes of the brightest red giants in 26 globular clusters, number counts of horizontal-branch stars and red giants in 15 globular clusters, and statistical parallax determinations of field RR Lyr luminosities yield delta(Mc) = 0.009 + or - 0.012 solar mass, so that conservatively mu(nu) is less than 3 x 10 to the -12th muB.

Raffelt, Georg G.↗

NASA NTRS1990

Energy conserving and potential-enstrophy dissipating schemes for the shallow water equations

To incorporate potential enstrophy dissipation into discrete shallow water equations with no or arbitrarily small energy dissipation, a family of finite-difference schemes have been derived with which potential enstrophy is guaranteed to decrease while energy is conserved (when the mass flux is nondivergent and time is continuous). Among this family of schemes, there is a member that minimizes the spurious impact of infinite potential vorticities associated with infinitesimal fluid depth. The scheme is, therefore, useful for problems in which the free surface may intersect with the lower boundary.

Arakawa, Akio↗

NASA NTRS1991

Simulation of mixing in the quick quench region of a rich burn-quick quench mix-lean burn combustor

A computer program was developed to study the mixing process in the quick quench region of a rich burn-quick quench mix-lean burn combustor. The computer program developed was based on the density-weighted, ensemble-averaged conservation equations of mass, momentum (full compressible Navier-Stokes), total energy, and species, closed by a k-epsilon turbulence model with wall functions. The combustion process was modeled by a two-step global reaction mechanism, and NO(x) formation was modeled by the Zeldovich mechanism. The formulation employed in the computer program and the essence of the numerical method of solution are described. Some results obtained for nonreacting and reacting flows with different main-flow to dilution-jet momentum flux ratios are also presented.

Shih, Tom I.-P.↗

NASA NTRS1991

Effect of ohmic, mass-transfer, and kinetic resistances on linear-sweep voltammetry in a cylindrical-pore electrode

A means is presented to account for the effect of ohmic, mass-transfer, and kinetic resistances on linear-sweep voltammograms by modeling a pore in a porous matrix as a cylindrical-pore electrode, and solving the mass and charge conservation equations in the context of this geometry for the simply redox reaction O + ne(-) yield R where both O and R are soluble species. Both analytical and numerical techniques are used to solve the governing equations. The calculated peak currents and potentials are correlated by empirical formulas to the measurable parameters: sweep rate, concentration of the redox species, diffusion coefficient, conductivity of the electrolyte, and pore dimensions. Using the correlations, a methodology is established for determining if the redox reaction kinetics are irreversible or reversible (Nernstian). If the reaction is irreversible, it is shown how the standard rate constant and the transfer coefficient may be extracted from linear-sweep voltammetry data, or, if the reaction is reversible, how the number of electrons transferred may be deduced.

Weidner, John W.↗