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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
NASA NTRS1985

Theoretical and experimental comparison of vapor cavitation in dynamically loaded journal bearings

Vapor cavitation for a submerged journal bearing under dynamically loaded conditions was investigated. The observation of vapor cavitation in the laboratory was done by high-speed photography. It was found that vapor cavitation occurs when the tensile stress applied to the oil exceeded the tensile strength of the oil or the binding of the oil to the surface. The theoretical solution to the Reynolds equation is determined numerically using a moving boundary algorithm. This algorithm conserves mass throughout the computational domain including the region of cavitation and its boundaries. An alternating direction implicit (MDI) method is used to effect the time march. A rotor undergoing circular whirl was studied. Predicted cavitation behavior was analyzed by three-dimensional computer graphic movies. The formation, growth, and collapse of the bubble in response to the dynamic conditions is shown. For the same conditions of dynamic loading, the cavitation bubble was studied in the laboratory using high-speed photography.

Brewe, D. E.↗
NASA NTRS1986

The construction and use of divergence free vector expansions for incompressible fluid flow calculations

For incompressible fluids the law of mass conservation reduces to a constraint on the velocity vector, namely that it be divergence free. This constraint has long been a source of great difficulty to the numericist seeking to discretize the Navier-Stokes and Euler equations. A spectral method is discussed which overcomes this difficulty. Its efficacy is demonstrated on some simple problems. The velocity is approximated by a finite sum of divergence free vectors, each of which satisfies the same boundary conditions as the velocity. Projecting the governing equation onto the space of inviscid vector fields eliminates the pressure term and produces a set of ordinary differential equations that must be solved for the coefficents in the velocity. The pressure can then be recovered if it is needed.

Mhuiris, N. M. G.↗
NASA NTRS1985

Theoretical modeling of the vapor cavitation in dynamically loaded journal bearings

A theoretical investigation is made of the evolution of a vapor-bubble for a submerged journal bearing under dynamically loaded conditions using the Elrod algorithm. This method conserves mass throughout the computational domain. A comparison study is performed to determine some of the consequences of applying a nonconservative theory (pseudo-Gumbel BC) to a dynamic problem. A complete dynamic cycle of a journal whirling in a circular path is chosen for the basis of comparison. Significant differences are observed in the load components near the end of the cycle. Further, good agreement with experiment is found for stationary and nonstationary cavitation.

Brewe, D. E.↗
NASA NTRS1986

Effect of vibration amplitude on vapor cavitation in journal bearings

Computational movies were used to analyze the formation and collapse of vapor cavitation bubbles in a submerged journal bearing. The effect of vibration amplitude on vapor cavitation was studied for a journal undergoing circular whirl. The boundary conditions were implemented using Elrod's algorithm, which conserves mass flow through the cavitation bubble as well as through the oil-film region of the bearing. The vibration amplitudes for the different cases studied resulted in maximum eccentricity ratios ranging from 0.4 to 0.9. The minimum eccentricity ratio reached in each case was 0.1. For the least vibration amplitude studied in which the eccentricity ratio varied between 0.1 and 0.4, no vapor cavitation occurred. The largest vibration amplitude (i.e., eccentricity ratios of 0.1 to 0.9) resulted in vapor cavitation present 76 percent of one complete orbit.

Brewe, D. E.↗
NASA NTRS1986

An assessment of transport coefficient approximations used in galactic heavy ion shielding calculations

An energy-dependent, perturbation expansion solution for heavy ion transport in one dimension is used to perform depth-dose calculations for 670/MeV nucleon Ne-20 beams incident upon a thick water target. Comparisons of predictions obtained by using typical energy-independent approximations and those obtained with fully energy-dependent input parameters are made. It is found that the calculated doses are underestimated when the energy-independent input approximations are used. The major source of error, however, is the lack of charge and mass conservation in the Silberberg-Tsao fragmentation parameters.

Townsend, Lawrence W.↗
NASA NTRS1987

Effect of shaft frequency on cavitation in a journal bearing for noncentered circular whirl

The effect of shaft frequency on the performance of a submerged journal undergoing noncentered circular whirl is examined. The main emphasis of the paper is on the behavior of the vapor cavitation bubble and its effect on the bearing performance as a function of frequency. A cavitation algorithm due to Elrod was implemented in a computer program which solves a time-dependent Reynolds equation. This algorithm automatically handles the boundary conditions by using a switch function and a control volume approach which conserves mass throughout the entire flow. The shaft frequencies in this investigation ranged from 0 rad/s (squeeze-film damper) to -104 rad/s (a case in which oil-whip condition was produced momentarily). For the particular vibration amplitude chosen in this investigation it was observed that vapor cavitation had an effect on the load components for the full range of shaft frequencies investigated.

Brewe, David E.↗
NASA NTRS1986

Nongeostrophic theory of zonally averaged circulation. I - Formulation

A nongeostrophic theory of zonally averaged circulation is formulated using the nonlinear primitive equations (mass conservation, thermodynamics, and zonal momentum) on a sphere. The relationship between the mean meridional circulation and diabatic heating rate is studied. Differences between results of nongeostropic theory and the geostrophic formulation concerning the role of eddy forcing of the diabatic circulation and the nonlinear nearly inviscid limit versus the geostrophic limit are discussed. Consideration is given to the Eliassen-Palm flux divergence, the Eliassen-Palm pseudodivergence, the nonacceleration theorem, and the nonlinear nongeostrophic Taylor relationship.

Tung, Ka Kit↗
NASA NTRS1987

A study of carbon monoxide distribution determinations for a global transport model

A system of programs which model the chemistry and transport of carbon monoxide and methane in the Earth's atmosphere on a global scale was installed onto the NASA-Langley central scientific computer network. This system, the GLOBAL system, consists of a user-friendly set of procedural files which allow for simplified pre-processing, execution, and post-processing for all program elements. The package includes procedures for obtaining the NMC meteorological data, calculating the vertical winds to satisfy mass conservation, determining the boundary layer, and executing the transport/chemical model for carbon monoxide. In addition, plotting, saving to tape, and reading from tape routines were developed. Final modifications to the subprograms for processing the input data for the transport/chemistry model have improved these data to more accurately reflect true atmospheric conditions. The transfer of the transport/chemistry model from the NCAR CRAY system to the NASA-Langley VPS-32 system was completed. The problems encountered during this process and their resolutions are discussed.

Peters, Leonard K.↗
NASA NTRS1988

Consistent boundary conditions for Reduced Navier-Stokes (RNS) scheme applied to 3-dimensional internal viscous flows

A consistent and efficient set of boundary conditions is developed for the multi-sweep space-marching pressure-elliptic Reduced Navier-Stokes (RNS) scheme as applied for 3-D internal viscous flow problems. No-slip boundary conditions are directly imposed on the solid walls. There is no iteration procedure required in the cross plane to ensure mass conservation across each marching plane. The finite difference equations forming the coefficient matrix are ordered such that the surface normal velocity is specified on all the solid walls; unlike external flows, a pressure boundary condition in the cross plane is not required. Since continuity is directly satisfied at all points in the flow domain, the first order momentum equations can be solved directly for the pressure without the need for a Poisson pressure correction equation. The procedure developed herein can also be applied with periodic boundary conditions. The analysis is given for general compressible flows. Incompressible flow solutions are obtained, for straight and curved ducts of square cross section, to validate the procedure. These solutions are used to demonstrate the applicability of the RNS scheme, with the improved boundary conditions for internal flows with strong interaction, as would be encountered in ducts and turbomachinery geometries.

Reddy, D. R.↗
NASA NTRS1988

An improved flux-split algorithm applied to hypersonic flows in chemical equilibrium

An explicit, finite-difference, shock-capturing numerical algorithm is presented and applied to hypersonic flows assumed to be in thermochemical equilibrium. Real-gas chemistry is either loosely coupled to the gasdynamics by way of a Gibbs free energy minimization package or fully coupled using species mass conservation equations with finite-rate chemical reactions. A scheme is developed that maintains stability in the explicit, finite-rate formulation while allowing relatively high time steps. The codes use flux vector splitting to difference the inviscid fluxes and employ real-gas corrections to viscosity and thermal conductivity. Numerical results are compared against existing ballistic range and flight data. Flows about complex geometries are also computed.

Palmer, Grant↗
NASA NTRS1988

Igneous proceses and closed system evolution of the Tharsis region of Mars

A quantitative petrologic model for the evolution of the Tharsis region on Mars is presented, which is consistent with global gravity and topography data. It is demonstrated that it is possible to form and support the topographic relief of the Tharsis plateau by a closed-system mass-conservative nearly isostatic process involving generation of magmas from a mantle source region. Extrusion and/or intrusion (or underplating) of such magmas allows low-pressure solidification, with a consequent increase in volume relative to that which would be possible in the high-pressure source region, leading to elevated topology. The distribution of densities with depth obtained by the model is quantitatively consistent with the isostatic models of Sleep and Phillips (1979, 1985).

Finnerty, A. A.↗
NASA NTRS1989

The free-wake prediction of rotor hover performance using a vortex embedding method

A method is developed to predict the rotor hover performance. This method solves the compressible mass conservation equation much like current full potential codes and can therefore predict the transonic flows on a rotor. However, the newly developed approach also allows for the free convection of shed vorticity and permits the computation of the entire hover wake system. The method uses a vortex embedding scheme in potential flow and has been implemented in a computer code, HELIX -I. To predict power we implement a simple boundary layer and two different induced-drag integration schemes. The induced-drag is obtained from surface pressure integration and an energy flux integral. Comparisons between computations and experiment show good agreement for the prediction of power polars, surface pressure distribution, and tip vortex geometry.

Ramachandran, K.↗
NASA NTRS1989

The ultraviolet variability of HD 45166 (qWR + B8 V) - Evidence for stellar wind radiative instabilities

High-resolution IUE spectroscopy of the low-mass qWR star HD 45166, obtained at several epochs during 1980-1988, reveals a complex pattern of variability in its stellar wind emission lines, P Cygni profiles, and highly ionized photospheric absorption spectrum. The most intensive observations, obtained during a 36 hr continuous run in February 1988, show the presence of discrete absorption components (DACs) in the C IV 1550 resonance doublet. These DACs can blend to give the appearance of a well-developed P Cygni absorption profile seen at some other times, roughly doubling the column density of the ground state C(3+) ions in the wind. In the 1988 spectra, two principal features are seen, at about -950 km/s and -750 km/s, with column densities of about 14.0 and 14.7 dex/sq cm, respectively. Both features migrate in velocity with accelerations of about 140 cm/sq sec, implying a recurrence time scale of about 1.6 days and a lifetime of about 3 days for the DACs. A striking characteristic of the acceleration is its slowness. The pattern of DAC variability in HD 45166 is similar to that found for O-type stars and provides further evidence that the DACs are not mass-conserving features but rather that a given feature is formed by different material at different times. It is concluded that the C IV DAC variability observed in HD 45166 results from structural changes in the wind arising from radiative instabilities.

Willis, Allan J.↗
NASA NTRS1989

A complete spherical harmonic approach to luni-solar tides

In this work a spherical harmonic theory of ocean tides is presented. The theory is based on Laplace tide equations modified to include turbulence with constant eddy viscosity, linearized bottom friction, and oceanic loading and self-gravitation. Variable bathymetry is also treated in harmonic terms, and no-flow boundary conditions are applied at continental coastlines. The tide and boundary constraint equations are reduced to matrix form and solved by a weighted least-squares procedure. Five zonal luni-solar tides, ranging in period from 14 days to 18.6 yr, are investigated using the theory; such tides have typically been difficult to compute using traditional numerical approaches. The polar motion and changes in the length of day induced by these long-period tides are calculated. Tidal solutions are compared extensively with results from other tidal theories and from recent satellite and sea-level observations. The greatest limitation to accurate prediction of zonal tides - for any theory - appears to be the marginal failure of all tide theories to conserve mass globally; the use of additional mass constraints may be warranted.

Dickman, S. R.↗
NASA NTRS1990

Relativistic winds from pulsar and black hole magnetospheres

The present consideration of the extraction of the rotational energy stored in a Kerr black hole by a large-scale MHD wind as a possible source for AGN jet energy proceeds by viewing both the ingoing and outgoing winds as propagating toward asymptotic infinities and transporting the conserved mass, angular momentum, and energy fluxes established by the source region. It is judged that the field-line rotation rate or total voltage drop across the system must also be determined by the dissipative process that injects plasma onto the magnetospheric field lines.

Punsly, Brian↗
NASA NTRS1991

A dynamically adaptive multigrid algorithm for the incompressible Navier-Stokes equations: Validation and model problems

An algorithm is described for the solution of the laminar, incompressible Navier-Stokes equations. The basic algorithm is a multigrid based on a robust, box-based smoothing step. Its most important feature is the incorporation of automatic, dynamic mesh refinement. This algorithm supports generalized simple domains. The program is based on a standard staggered-grid formulation of the Navier-Stokes equations for robustness and efficiency. Special grid transfer operators were introduced at grid interfaces in the multigrid algorithm to ensure discrete mass conservation. Results are presented for three models: the driven-cavity, a backward-facing step, and a sudden expansion/contraction.

Thompson, C. P.↗
NASA NTRS1991

Ion dynamics in the Venus ionosphere

Measurement data on the ion velocity in the Venus ionosphere (mainly from the Pioneer Venus Orbiter Retarding Potential Analyzer) are summarized, and theoretical models developed to explain them are reviewed. Data and theoretical predictions are compared in extensive graphs and diagrams and discussed in detail. It is shown that the predominant flow is away from the subsolar point, at up to 3 km/sec in the terminator region. A model of axisymmetric flow based on momentum, energy, and mass conservation laws is found to reproduce the observed ion velocities at solar zenith angles less than about 140 deg, but not the high velocities and chaotic behavior seen near the antisolar point. Also discussed are significant differences between the flow above and below about 400 km and the effects of changes in the dynamic pressure of the solar wind.

Miller, K. L.↗
NASA NTRS1991

The turbomachine blading design using S2-S1 approach

The boundary conditions corresponding to the design problem when the blades being simulated by the bound vorticity distribution are presented. The 3D flow is analyzed by the two steps S2 - S1 approach. In the first step, the number of blades is supposed to be infinite, the vortex distribution is transformed into an axisymmetric one, so that the flow field can be analyzed in a meridional plane. The thickness distribution of the blade producing the flow channel striction is taken into account by the modification of metric tensor in the continuity equation. Using the meridional stream function to define the flow field, the mass conservation is satisfied automatically. The governing equation is deduced from the relation between the azimuthal component of the vorticity and the meridional velocity. The value of the azimuthal component of the vorticity is provided by the hub to shroud equilibrium condition. This step leads to the determination of the axisymmetric stream sheets as well as the approximate camber surface of the blade. In the second step, the finite number of blades is taken into account, the inverse problem corresponding to the blade to blade flow confined in each stream sheet is analyzed. The momentum equation implies that the free vortex of the absolute velocity must be tangential to the stream sheet. The governing equation for the blade to blade flow stream function is deduced from this condition. At the beginning, the upper and the lower surfaces of the blades are created from the camber surface obtained from the first step with the assigned thickness distribution. The bound vorticity distribution and the penetrating flux conservation applied on the presumed blade surface constitute the boundary conditions of the inverse problem. The detection of this flux leads to the rectification of the geometry of the blades.

Luu, T. S.↗