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 73 records · Page 4

Nonequilibrium ionization due to thermal diffusion and mass flows

Recent calculations of diffusion coefficients are used in the continuity equation to compute ion populations of carbon in the solar transition region. Thermal diffusion causes strong departures from ionization equilibrium in the region where the temperature gradient is steepest. Mass-conserving flows are also included in our calculations. These dominate over thermal diffusion depending on the magnitude of the flows and also lead to departures from ionization equilibrium. These results have important implications for the interpretation of EUV line emission.

Roussel-Dupre, R.↗

Coronal closed structures. IV - Hydrodynamical stability and response to heating perturbations

The response of magnetically confined atmospheres to perturbations in the temperature and density distribution, and the local heating rate by means of a one-dimensional time-dependent hydrodynamical code, which incorporates the full energy, momentum and mass conservation equations is studied. These studies extend the linear instability analysis of Habbal and Rosner (1979) into the finite-amplitude regime, and generalize the confined atmosphere models of Serio et al., to the time-dependent domain. The results show that closed coronal atmospheres are stable against finite-amplitude perturbations if the chromospheric response is taken into account; and observed correlated increases in coronal density and temperature can only be achieved under quiescent conditions by increasing the heat deposition rate relatively more in the chromosphere than in the corona.

Peres, G.↗

Closed coronal structures. V - Gasdynamic models of flaring loops and comparison with SMM observations

A time-dependent one-dimensional code incorporating energy, momentum and mass conservation equations, and taking the entire solar atmospheric structure into account, is used to investigate the hydrodynamic response of confined magnetic structures to strong heating perturbations. Model calculation results are compared with flare observations which include the light curves of spectral lines formed over a wide range of coronal flare temperatures, as well as determinations of Doppler shifts for the high temperature plasma. It is shown that the numerical simulation predictions are in good overall agreement with the observed flare coronal plasma evolution, correctly reproducing the temporal profile of X-ray spectral lines and their relative intensities. The predicted upflow velocities support the interpretation of the blueshifts as due to evaporation of chromospheric material.

Peres, G.↗

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

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

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

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

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

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

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

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

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

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

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

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↗

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

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

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↗