Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Equation of State”

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 343 records · Page 19

Design and application of squeeze film dampers for turbomachinery stabilization

The steady-state transient response of the squeeze film damper bearing was investigated. Both the steady-state and transient equations for the hydrodynamic bearing forces are derived; the steady-state equations were used to determine the damper equivalent stiffness and damping coefficients. These coefficients are used to find the damper configuration which will provide the optimum support characteristics based on a stability analysis of the rotor-bearing system. The effects of end seals and cavitated fluid film are included. The transient analysis of rotor-bearing systems was conducted by coupling the damping and rotor equations and integrating forward in time. The effects of unbalance, cavitation, and retainer springs are included. Methods of determining the stability of a rotor-bearing system under the influence of aerodynamic forces and internal shaft friction are discussed.

Gunter, E. J.↗

Implications of tachyon-like matter for superdense stars.

Derivation of a new equation of state of superdense matter by treating superdense matter as a perfect, degenerate tachyon gas. Model calculations for superdense stars based on this equation of state are presented. By appropriately choosing a certain parameter, dynamical stability can be achieved for arbitrarily large central densities. Also, a somewhat larger than usual value for the maximum mass is obtained.

Bhatia, M. S.↗

Acceleration to a steady state for the Euler equations

A multi-stage Runge-Kutta method is analyzed for solving the Euler equations exterior to an airfoil. Highly subsonic, transonic and supersonic flows are evaluated. Various techniques for accelerating the convergence to a steady state are introduced and analyzed.

Turkel, E.↗

Acceleration to a steady state for the Euler equations

A multistage Runge-Kutta method is analyzed for solving the Euler equations exterior to an airfoil. Highly subsonic, transonic and supersonic flows are evaluated. Various techniques for accelerating the convergence to a steady state are introduced and analyzed.

Turkel, E.↗

An alternative approach to the numerical simulation of steady inviscid flow

A numerical procedure for the efficient simulation of steady inviscid flow is described and its utility is demonstrated. The method is uniformly valid for application in the subsonic, transonic and supersonic flow regimes. It does not rely on the introduction of additional assumptions beyond those necessary to obtain the Euler equations from the Navier-Stokes equations, nor does it make use of a time-asymptotic solution of the unsteady equations of motion. Application of the herein-defined surrogate equation technique allows the formulation of stable, fully-conservative, type-dependent finite difference equations for use in obtaining numerical solutions to systems of first-order partial differential equations, such as the steady-state Euler equations or their various approximations. Computational results are presented for the full Euler equations used to simulate rotational subsonic flow and for the transonic small disturbance equations. For the latter case, a computational efficiency greater than that obtained by means of the standard perturbation potential approach is indicated.

Johnson, G. M.↗

The origin of the moon and the single-impact hypothesis. I

One of the newer ideas regarding the origin of the moon is concerned with a single-impact hypothesis. It is pointed out that this theory has the advantage of overcoming most of the difficulties with the classical theories. The angular momentum of the earth-moon system can easily be obtained by varying the initial conditions of the impact. A series of three-dimensional numerical simulations of the collision between the earth and an object of about 1/10 its mass is presented. Different impact velocities, impact parameters, and initial internal energies are considered. Attention is given to assumptions, the equation of state, numerical techniques utilizing the momentum equation and the energy conservation equation, tests, and initial conditions and units.

Benz, W.↗

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

One-way nesting for a primitive equation ocean model

Prognostic numerical models for atmospheric and oceanic circulations require initial fields, boundary conditions, and forcing functions in addition to a consistent set of partial differential equations, including a state relation and equations expressing conservation of mass, momentum, and energy. Depending on the horizontal domain to be modeled, the horizontal boundary conditions are either physically obvious or extremely difficult to specify consistently. If the entire atmosphere is modeled, periodic horizontal boundary conditions are appropriate. On the other hand, the physical horizontal boundaries on the entire ocean are solid walls. Obviously, the normal velocity at a solid wall is zero while the specification of the tangential velocity depends on the mathematical treatment of the horizontal viscous terms. Limitations imposed by computer capacity and cost, as well as research interests, have led to the use of limited area models to study flows in the atmosphere and ocean. The limited area models do not have physical horizontal boundaries, merely numerical ones. Correctly determining these open boundary conditions for limited-area numerical models has both intrigued and frustrated numerical modelers for decades. One common approach is to use the closed or solid wall boundary conditions for a limited-area model. The argument given for this approach is that the boundary conditions affect flow near the walls but that none of these effects are propagated into the interior. Therefore, one chooses a big enough domain that the central region of interest is not corrupted by the boundary flow. Research in progress to model the North Atlantic circulation vividly illustrates the pitfalls of this approach. Two model runs are compared: (1) the southern boundary at 20S between latitudes 0 and 40W is artificially closed; and (2) the same boundary is specified as open with an inward transport of 15 Sv (determined from a global model with the same physics) uniformly spread across the boundary. A comparison of both runs is presented.

Blake, D. W.↗

Sources of uncertainty in direct seismological measurements of the solar helium abundance

The methods by which Dappen et al. (1988, 1990, 1991) and Dziembowski et al. (1990, 1991, 1992) recently obtained discrepant estimates of the helium abundance in the solar convection zone are compared. The aim of the investigation reported in this paper is to identify the main source of the discrepancy. Using as proxy data eigenfrequencies of a set of modes of a theoretical solar model, computed with the same physics as were the frequencies of a reference model with which these data are compared, the two methods yield similar results. Thus we ascertain that the principal source of the discrepancy is not in the inversions themselves, which yield essentially a measure of the variation of the adiabatic exponent gamma of the material in the He II ionization zone. Instead it is in the approximations adopted in the treatment of heavy elements in the equation of state used to relate the variation of gamma to chemical composition. We obtain acceptably consistent results when inverting solar data by two methods using the same equation of state. We attempt to identify the likely residual sources of uncertainty.

Kosovichev, A. G.↗

Second coefficient of viscosity in air

Acoustic attenuation measurements in air were analyzed in order to estimate the second coefficient of viscosity. Data over a temperature range of 11 C to 50 C and at relative humidities between 6 percent and 91 percent were used. This analysis showed that the second coefficient of viscosity varied between 1900 and 20,000 times larger than the dynamic or first coefficient of viscosity over the temperature and humidity range of the data. In addition, the data showed that the molecular relaxation effects, which are responsible for the magnitude of the second coefficient of viscosity, place severe limits on the use of time-independent, thermodynamic equations of state. Compressible flows containing large streamwise velocity gradients, like shock waves, which cause significant changes in particle properties to occur during time intervals shorter than hundredths of seconds, must be modeled using dynamic equations of state. The dynamic model approach is described briefly.

Ash, Robert L.↗

Experimental Shock Transformation of Gypsum to Anhydrite: A New Low Pressure Regime Shock Indicator

The shock behavior of gypsum is important in understanding the Cretaceous/Paleogene event and other terrestrial impacts that contain evaporite sediments in their targets (e.g., Mars Exploration Rover Spirit detected sulfate at Gusev crater, [1]). Most interest focuses on issues of devolatilization to quantify the production of SO2 to better understand its role in generating a temporary atmosphere and its effects on climate and biota [2,3]. Kondo and Ahrens [4] measured induced radiation emitted from single crystal gypsum shocked to 30 and 40 GPa. They observed greybody emission spectra corresponding to temperatures in the range of 3,000 to 4,000 K that are a factor of 2 to 10 times greater than calculated pressure-density energy equation of state temperatures (Hugoniot) and are high enough to melt gypsum. Chen et al. [5] reported results of shock experiments on anhydrite, gypsum, and mixtures of these phases with silica. Their observations indicated little or no devolatilization of anhydrite shocked to 42 GPa and that the fraction of sulfur, by mass, that degassed is approx.10(exp -2) of theoretical prediction. In another report of shock experiments on calcite, anhydrite, and gypsum, Badjukov et al. [6] observed only intensive plastic deformation in anhydrite shock loaded at 63 GPa, and gypsum converted to anhydrite when shock loaded at 56 GPa but have not experimentally shocked gypsum in a step-wise manner to constrain possible incipient transformation effects. Schmitt and Hornemann [7] shock loaded anhydrite and quartz to a peak pressure of 60 GPa and report the platy anhydrite grains were completely pseudomorphed by small crystallized anhydrite grains. However, no evidence of interaction between the two phases could be observed and they suggested that recrystallization of anhydrite grains is the result of a solid-state transformation. They concluded that significant decomposition of anhydrite requires shock pressures higher than 60 GPa. Gupta et al. [8] reanalyzed the calcite and anhydrite shock wave experiments of Yang [9] using improved equations of state of porous materials and vaporized products. They determined the pressures for incipient and complete vaporization to be 32.5 and 122 GPa for anhydrite GPa which is a factor of 2 to 3 lower than reported earlier by Yang [9]. These studies are not in agreement regarding the onset of sulfate decomposition and documentation of shock effects in gypsum is incomplete.

Bell, Mary S.↗

Use of Generalized Fluid System Simulation Program (GFSSP) for Teaching and Performing Senior Design Projects at the Educational Institutions

This paper describes the experience of the authors in using the Generalized Fluid System Simulation Program (GFSSP) in teaching Design of Thermal Systems class at University of Alabama in Huntsville. GFSSP is a finite volume based thermo-fluid system network analysis code, developed at NASA/Marshall Space Flight Center, and is extensively used in NASA, Department of Defense, and aerospace industries for propulsion system design, analysis, and performance evaluation. The educational version of GFSSP is freely available to all US higher education institutions. The main purpose of the paper is to illustrate the utilization of this user-friendly code for the thermal systems design and fluid engineering courses and to encourage the instructors to utilize the code for the class assignments as well as senior design projects. The need for a generalized computer program for thermofluid analysis in a flow network has been felt for a long time in aerospace industries. Designers of thermofluid systems often need to know pressures, temperatures, flow rates, concentrations, and heat transfer rates at different parts of a flow circuit for steady state or transient conditions. Such applications occur in propulsion systems for tank pressurization, internal flow analysis of rocket engine turbopumps, chilldown of cryogenic tanks and transfer lines, and many other applications of gas-liquid systems involving fluid transients and conjugate heat and mass transfer. Computer resource requirements to perform time-dependent, three-dimensional Navier-Stokes computational fluid dynamic (CFD) analysis of such systems are prohibitive and therefore are not practical. Available commercial codes are generally suitable for steady state, single-phase incompressible flow. Because of the proprietary nature of such codes, it is not possible to extend their capability to satisfy the above-mentioned needs. Therefore, the Generalized Fluid System Simulation Program (GFSSP1) has been developed at NASA Marshall Space Flight Center (MSFC) as a general fluid flow system solver capable of handling phase changes, compressibility, mixture thermodynamics and transient operations. It also includes the capability to model external body forces such as gravity and centrifugal effects in a complex flow network. The objectives of GFSSP development are: a) to develop a robust and efficient numerical algorithm to solve a system of equations describing a flow network containing phase changes, mixing, and rotation; and b) to implement the algorithm in a structured, easy-to-use computer program. The analysis of thermofluid dynamics in a complex network requires resolution of the system into fluid nodes and branches, and solid nodes and conductors as shown in Figure 1. Figure 1 shows a schematic and GFSSP flow circuit of a counter-flow heat exchanger. Hot nitrogen gas is flowing through a pipe, colder nitrogen is flowing counter to the hot stream in the annulus pipe and heat transfer occurs through metal tubes. The problem considered is to calculate flowrates and temperature distributions in both streams. GFSSP has a unique data structure, as shown in Figure 2, that allows constructing all possible arrangements of a flow network with no limit on the number of elements. The elements of a flow network are boundary nodes where pressure and temperature are specified, internal nodes where pressure and temperature are calculated, and branches where flowrates are calculated. For conjugate heat transfer problems, there are three additional elements: solid node, ambient node, and conductor. The solid and fluid nodes are connected with solid-fluid conductors. GFSSP solves the conservation equations of mass and energy, and equation of state in internal nodes to calculate pressure, temperature and resident mass. The momentum conservation equation is solved in branches to calculate flowrate. It also solves for energy conservation equations to calculate temperatures of solid nodes. The equations are coupled and nonlinear; therefore, they are solved by an iterative numerical scheme. GFSSP employs a unique numerical scheme known as simultaneous adjustment with successive substitution (SASS), which is a combination of Newton-Raphson and successive substitution methods. The mass and momentum conservation equations and the equation of state are solved by the Newton-Raphson method while the conservation of energy and species are solved by the successive substitution method. GFSSP is linked with two thermodynamic property programs, GASP2 and WASP3 and GASPAK4, that provide thermodynamic and thermophysical properties of selected fluids. Both programs cover a range of pressure and temperature that allows fluid properties to be evaluated for liquid, liquid-vapor (saturation), and vapor region. GASP and WASP provide properties of 12 fluids. GASPAK includes a library of 36 fluids. GFSSP has three major parts. The first part is the graphical user interface (GUI), visual thermofluid analyzer of systems and components (VTASC). VTASC allows users to create a flow circuit by a 'point and click' paradigm. It creates the GFSSP input file after the completion of the model building process. GFSSP's GUI provides the users a platform to build and run their models. It also allows post-processing of results. The network flow circuit is first built using three basic elements: boundary node, internal node, and branch.

Majumdar, A. K.↗

Fluid dynamics of multiphase flow in volcanic environments on the Earth, Io, and other planets

Data was compiled on H2O, CO2, S, and SO2 and it was analyzed by various techniques to prepare phase diagrams relevant to conditions on the Earth, Io, Mars, and other planets. The use of Cray models for volcanic eruptions was examined. Attempts are being made to develop equations of state of the compounds listed above to be compatible with the Cray code so that they can be called as subroutines for the volcanic models. This is a major project, as the current code uses only the equation of state for a perfect gas, and requires that the atmosphere be the same composition as the volcanic gas. The analysis of the thermodynamic data is time consuming, and was complicated by changes from a VAX/VMS system to Macintosh's, Crays using CTSS, Crays using Unicos, and Sun workstations.

Kieffer, Susan Werner↗

Outer magnetospheric fluctuations and pulsar timing noise

The Cheng, Ho, and Ruderman (1986) outer-magnetosphere gap model was used to investigate the stability of Crab-type outer magnetosphere gaps for pulsars having the parameter (Omega-square B) similar to that of the Crab pulsar. The Lamb, Pines, and Shaham (1978) fluctuating magnetosphere noise model was applied to the Crab pulsar to examine the type of the equation of state that best describes the structure of the neutron star. The noise model was also applied to other pulsars, and the theoretical results were compared with observational data. The results of the comparison are consistent with the stiff equation of state, as suggested by the vortex creep model of the neutron star interior. The timing-noise observations also contribute to the evidence for the existence of superfluid in the core of the neutron star.

Cheng, K. S.↗

Rotating neutron star models and pulsars

Recently developed equations of state are used to construct fully relativistic models of rotating neutron stars. All of these stable neutron star models are bound. Consideration of binding energies suggests that very-low-mass neutron stars will not be formed in supernova explosions. Some criteria for equations of state are derived from a comparison with observations. Limits for the physical parameters of the Crab pulsar are obtained. In addition, we find that cosmic-ray protons cannot be produced by pulsars.

Boerner, G.↗