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 163 records · Page 9

All sources NASA NTRS DOE OSTI

NASA NTRS1983

A method for solving the transonic full-potential equation for general configurations

A method is developed for solving the full-potential equation for two-dimensional and axisymmetric flow which retains the grid and boundary condition simplicity of the transonic small-disturbance codes. The method is based on a finite-volume formulation of the mass conservation equation in a Cartesian coordinate system, and is an extension of the method of Purvis and Burkhalter (1979). This finite-volume approach, combined with the simple boundary treatment, is shown to result in a highly robust method applicable to a wide range of geometries and flow conditions. The accuracy of the method is demonstrated for general geometries in two-dimensional and axisymmetric flows. The use of this method results in significant gains in convergence rate over the vertical-line over-relaxation scheme by incorporating an AF2-type algorithm (Ballhaus et al., 1978). It is suggested that the simplicity of this method shold allow a relatively easy extension to complex geometries in three-dimensional flows, and complex two-dimensional configurations such as multielement airfoils should be amenable to this method.

Wedan, B.↗

NASA NTRS1984

A computational method for viscous incompressible flows

An implicit, finite-difference procedure for numerically solving viscous incompressible flows is presented. The pressure-field solution is based on the pseudocompressibility method in which a time-derivative pressure term is introduced into the mass-conservation equation to form a set of hyperbolic equations. The pressure-wave propagation and the spreading of the viscous effect is investigated using simple test problems. Computed results for external and internal flows are presented to verify the present method which has proved to be very robust in simulating incompressible flows.

Kwak, D.↗

NASA NTRS1984

A solution procedure for three-dimensional incompressible Navier-Stokes equation and its application

An implicit, finite-difference procedure is presented for numerically solving viscous incompressible flows. For convenience of applying the present method to three-dimensional problems, primitive variables, namely the pressure and velocities, are used. One of the major difficulties in solving incompressible flows that use primitive variables is caused by the pressure field solution method which is used as a mapping procedure to obtain a divergence-free velocity field. The present method is designed to accelerate the pressure-field solution procedure. This is achieved by the method of pseudocompressibility in which the time derivative pressure term is introduced into the mass conservation equation. The pressure wave propagation and the spreading of the viscous effect is investigated using simple test problems. The present study clarifies physical and numerical characteristics of the pseudo-compressible approach in simulating incompressible flows. Computed results for external and internal flows are presented to verify the present procedure. The present algorithm has been shown to be very robust and accurate if the selection of the pseudo-compressibility parameter has been made according to the guidelines given.

Kwak, D.↗

NASA NTRS1984

An incompressible Navier-Stokes flow solver in three-dimensional curvilinear coordinate systems using primitive variables

An implicit, finite-difference computer code has been developed to solve the incompressible Navier-Stokes equations in a three-dimensional, curvilinear coordinate system. The pressure-field solution is based on the pseudo compressibility approach in which the time derivative pressure term is introduced into the mass conservation equation to form a set of hyperbolic equations. The solution procedure employs an implicit, approximate factorization scheme. The Reynolds stresses, that are uncoupled from the implicit scheme, are lagged by one time-step to facilitate implementing various levels of the turbulence model. Test problems for external and internal flows are computed, and the results are compared with existing experimental data. The application of this technique for general three-dimensional problems is then demonstrated.

Kwak, D.↗

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 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

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.↗