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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Thermal evolution of neutron stars with internal frictional heating

It has been suggested that the frictional interaction of neutron superfluids with normal matter in the inner crust of neutron stars dissipates rotational energy of superfluids and generates heat. Incorporating a general formula of internal heating into the detailed numerical codes of thermal evolution, we examine the effects of the internal heating on thermal evolution of neutron stars. We find that when a very stiff equation of state is used, it takes as long as about 2 x 10 exp 4 yr for the interior of a neutron star to reach the isothermal state, even if a strong heat source is placed in a thin layer of the inner crust. This time scale reduces to a few hundred years or less when medium to soft equations of state are used. A neutron star cools by neutrino emissions during the earlier stages referred to as the neutrino cooling era, while during the later photon cooling era it cools primarily by emission of photons from its surface. We show that heating rates expected in the current superfluid-crust interaction model can greatly increase the surface temperature in the photon cooling era, significantly changing the thermal evolution of relatively old neutron stars if a stiff equation of state is adopted.

Umeda, Hideyuki↗
NASA NTRS1991

Electrostatic potential jump across fast-mode collisionless shocks

The electrostatic potential jump across fast-mode collisionless shocks is examined by comparing published observations, hybrid simulations, and a simple model, in order to better characterize its dependence on the various shock parameters. In all three, it is assumed that the electrons can be described by an isotropic power-law equation of state. The observations show that the cross-shock potential jump correlates well with the shock strength but shows very little correlation with other shock parameters. Assuming that the electrons obey an isotropic power law equation of state, the correlation of the potential jump with the shock strength follows naturally from the increased shock compression and an apparent dependence of the power law exponent on the Mach number which the observations indicate. It is found that including a Mach number dependence for the power law exponent in the electron equation of state in the simple model produces a potential jump which better fits the observations. On the basis of the simulation results and theoretical estimates of the cross-shock potential, it is discussed how the cross-shock potential might be expected to depend on the other shock parameters.

Mandt, M. E.↗
NASA NTRS1980

Structure and evolution of Uranus and Neptune

Three-layer interior models of Uranus and Neptune with central rocky cores, mantles of water, methane, and ammonia (the 'ices'), and outer envelopes primarily composed of hydrogen and helium are presented. The models incorporate a new H2O equation of state based on experimental data which is considerably 'softer' than previous H2O equations of state. Corrections for interior temperatures approximately 5000 K are included in the models, and the thermal evolution of both planets is investigated using recent heat flow measurements. It is found that the evolutionary considerations are consistent with gravitational field data in supporting models with approximately solar abundances of 'ice' and 'rock'. Evolutionary considerations indicate that initial temperatures and luminosities for Uranus and Neptune were not substantially higher than the present value. Both planets apparently have relatively small approximately 1-2 earth masses) hydrogen-helium envelopes, with Neptune's envelope smaller than Uranus'. A monotonic trend is evident among the Jovian planets: all have central rock-ice cores of approximately 15 earth masses, but with hydrogen-helium envelopes which decrease in mass from Jupiter to Saturn to Uranus to Neptune.

Hubbard, W. B.↗
NASA NTRS2013

The W-W02 Oxygen Fugacity Buffer at High Pressures and Temperatures: Implications for f02 Buffering and Metal-silicate Partitioning

Oxygen fugacity (fO2) controls multivalent phase equilibria and partitioning of redox-sensitive elements, and it is important to understand this thermodynamic parameter in experimental and natural systems. The coexistence of a metal and its oxide at equilibrium constitutes an oxygen buffer which can be used to control or calculate fO2 in high pressure experiments. Application of 1-bar buffers to high pressure conditions can lead to inaccuracies in fO2 calculations because of unconstrained pressure dependencies. Extending fO2 buffers to pressures and temperatures corresponding to the Earth's deep interior requires precise determinations of the difference in volume (Delta) V) between the buffer phases. Synchrotron x-ray diffraction data were obtained using diamond anvil cells (DAC) and a multi anvil press (MAP) to measure unit cell volumes of W and WO2 at pressures and temperatures up to 70 GPa and 2300 K. These data were fitted to Birch-Murnaghan 3rd-order thermal equations of state using a thermal pressure approach; parameters for W are KT = 306 GPa, KT' = 4.06, and αKT = 0.00417 GPa K-1. Two structural phase transitions were observed for WO2 at 4 and 32 GPa with structures in P21/c, Pnma and C2/c space groups. Equations of state were fitted for these phases over their respective pressure ranges yielding the parameters KT = 190, 213, 300 GPa, KT' = 4.24, 5.17, 4 (fixed), and αKT = 0.00506, 0.00419, 0.00467 GPa K-1 for the P21/c, Pnma and C2/c phases, respectively. The W-WO2 buffer (WWO) was extended to high pressure by inverting the W and WO2 equations of state to obtain phase volumes at discrete pressures (1-bar to 100 GPa, 1 GPa increments) along isotherms (300 to 3000K, 100 K increments). The slope of the absolute fO2 of the WWO buffer is positive with increasing temperature up to approximately 70 GPa and is negative above this pressure. The slope is positive along isotherms from 1000 to 3000K with increasing pressure up to at least 100 GPa. The WWO buffer is at a higher fO2 than the IW buffer at pressures lower than 40 GPa, and the magnitude of this difference decreases at higher pressures. This qualitatively indicates an increasingly lithophile character for W at higher pressures. The WWO buffer was quantitatively applied to W metal-silicate partitioning by using the WWO-IW buffer difference in combination with literature data on W metal-silicate partitioning to model the exchange coefficient (KD) for the Fe-W exchange reaction. This approach captures the pressure dependence of W metal-silicate partitioning using the WWO-IW buffer difference and models the activities of the components in the silicate and metallic phases using an expression of the Gibbs excess energy of mixing. Calculation of KD along a peridotite liquidus predicts a decrease in W siderophility at higher pressures that supports the qualitative behavior predicted by the WWO-IW buffer difference, and agrees with findings of others. Comparing the competing effects of temperature and pressure on W metal-silicate partitioning, our results indicate that pressure exerts a greater effect.

Shofner, G. A.↗
NASA NTRS2020

Thermal Analysis, Compressibility, and Decomposition of Synthetic Bastnäsite-(La) to Lanthanum Oxyfluoride

Understanding basic material properties of rare earth element (REE) bearing minerals such as their phase stability and equations of state can assist in understanding how economically viable deposits might form. Bastnäsite is the most commonly mined REE bearing mineral. We synthesized the lanthanum-fluoride end member, bastnäsite-(La) (LaCO3F), and investigated its thermal behavior and decomposition products from 298 K to 1173 K under ambient pressure conditions through thermogravimetric analysis, differential scanning calorimetry, evolved gas analysis, and high temperature powder X-ray diffraction. We also investigated the compressibility of bastnäsite-(La) via single crystal X-ray diffraction in diamond anvil cells at an ambient temperature up to 11.3 GPa and from 4.9 GPa to 7.7 GPa up to 673 K. At ambient pressure, bastnäsite-(La) was stable up to 598 K in air, where it decomposed into CO2 and tetragonal γ-LaOF. Above 948 K, cubic α-LaOF is stable. High temperature X-ray diffraction data were used to fit the Fei thermal equation of state and the thermal expansion coefficient α(sub 298) for all three materials. Bastnäsite-(La) was fit from 298 K to 723 K with V(sub 0) = 439.82 Å(exp 3), α(sub 298) = 4.32 × 10(exp -5) K(exp -1), a(sub 0) = −1.68 × 10(exp -5) K(exp -1), a(sub 1) = 8.34 × 10(exp -8) K(exp -1), and a(sub 2) = 3.126 K(exp -1). Tetragonal γ-LaOF was fit from 723 K to 948 K with V(sub 0) = 96.51 Å(exp 3), α(sub 298) = 2.95×10(exp -4) K(exp -1), a(sub 0) = −2.41×10(exp -5) K(exp -1), a(sub 1) = 2.42×10(exp -7) K(exp -1), and a(sub 2) = 41.147 K(exp -1). Cubic α-LaOF was fit from 973 K to 1123 K with V(sub 0) = 190.71 Å(exp 3), α(sub 298) = −1.12×10(exp -5) K(exp -1), a(sub 0) = 2.36×10(exp -4) K(exp -1), a(sub 1) = −1.73 × 10(exp -7) K(exp -1), and a(sub 2) = −17.362 K(exp -1). An ambient temperature third order Birch–Murnaghan equation of state was fit with V(sub 0) = 439.82 Å(exp 3), K(sub 0) = 105 GPa, and K’ = 5.58.

Bastnäsite↗
NASA NTRS1984

Rate equations analysis of phase-locked semiconductor laser arrays under steady state conditions

Rate equations analysis of phase-locked semiconductor laser arrays has been carried out. It was found that for given (laser) current densities, the photon density distribution in the array elements is that particular one which maximizes the total photon density. The results of this analysis were then combined with the waveguide properties of the laser array waveguide, yielding a basic model of phase-locked diode laser arrays. This model explains the effects of the variation of the current combination through the array elements on its mode structure that were observed recently.

Katz, J.↗
NASA NTRS1995

New model atmospheres for very cool white dwarfs with mixed H/He and pure He compositions

New model atmosphere calculations for very cool white dwarfs with mixed H/He and pure He compositions are presented. The hydrogen-rich models incorporate improved cross section calculations of the collision-induced absorption by molecular hydrogen due to collisions with H2, H, and He. The effects associated with variations in the effective temperature (T(sub eff) greater than or equal to 4000 and less than or equal to 10,000 K), the surface gravity (log g greater than or equal to 7.5 and less than or equal to 9.5), and the chemical composition (N(He)/N(H) greater than or equal to O and less than or equal to 100) are investigated. Results from earlier calculations are confirmed qualitatively, but a more detailed comparison reveals large quantitative deviations. Cool white dwarfs with mixed H/He chemical compositions are shown to be easily recognizable from their predicted strong infrared flux deficiency. Pure helium model calculations are described as well. These include a modified version of the recently developed equation of state of D. Saumon and G. Chabrier. Nonideal effects brought about by various equations of state are explored in detail. For the purpose of this analysis, a model of pressure ionization based on an accurate description of the interactions in a mostly atomic helium fluid is developed. The effects of pressure ionization are shown to be the most important issue in the model calculations. A critical discussion of previous generations of pure helium model calculations is presented. Finally, broadband color indices are provided for the complete model grid.

Bergeron, P.↗
NASA NTRS1980

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 demonstrated. Application of the 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 and for the transonic 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.↗
NASA NTRS1980

A mathematical model for Vertical Attitude Takeoff and Landing (VATOL) aircraft simulation. Volume 2: Model equations and base aircraft data

Equations incorporated in a VATOL six degree of freedom off-line digital simulation program and data for the Vought SF-121 VATOL aircraft concept which served as the baseline for the development of this program are presented. The equations and data are intended to facilitate the development of a piloted VATOL simulation. The equation presentation format is to state the equations which define a particular model segment. Listings of constants required to quantify the model segment, input variables required to exercise the model segment, and output variables required by other model segments are included. In several instances a series of input or output variables are followed by a section number in parentheses which identifies the model segment of origination or termination of those variables.

Fortenbaugh, R. L.↗
NASA NTRS1981

Surrogate-equation technique for simulation of steady inviscid flow

A numerical procedure for the iterative solution of inviscid flow problems is described, and its utility for the calculation of steady subsonic and transonic flow fields is demonstrated. Application of the surrogate equation technique defined herein 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. Steady, two dimensional solutions to the Euler equations for both subsonic, rotational flow and supersonic flow and to the small disturbance equations for transonic flow are presented.

Johnson, G. M.↗
NASA NTRS2006

MEASURING THE MASS OF 4UO900-40 DYNAMICALLY

Accurate measurements of neutron star masses are needed to constrain the equation of state of neutron star matter - of importance to both particle physics and the astrophysics of neutron stars - and to identify the evolutionary track of the progenitor stars that form neutron stars. The best measured values of the mass of 4UO900-40 (= Vela XR-l), 1.86 +/- 0.16 Msun (Barziv et al. 2001) and 1.93 +/- 0.20 Msun (Abubekerov et al. 2004), make it a leading candidate for the most massive neutron star known. The direct relationship between the maximum mass of neutron stars and the equation of state of ultra-dense matter makes 4UO900-40 an important neutron star mass to determine accurately. The confidence interval on previous mass estimates, obtained from observations that include parameters determined by non-dynamical methods, are not small enough to significantly restrict possible equations of state. We describe here a purely dynamical method for determining the mass of 4UO900-40, an X-ray pulsar, using the reprocessed UV pulses emitted by its BO.5Ib companion. One can derive the instantaneous radial velocity of each component by simultaneous X-ray and UV observations at the two quadratures of the system. The Doppler shift caused by the primary's rotational velocity and the illumination pattern of the X-rays on the primary, two of the three principal contributors to the uncertainty on the derived mass of the neutron star, almost exactly cancel by symmetry in this method. A heuristic measurement of the mass of 4UO900-40 using observations obtained previously with the High Speed Photometer on HST is given in Appendix A.

Dolan, J. F.↗
NASA NTRS1994

Storage and recycling of water and carbon dioxide in the earth

The stabilities and properties of water- and carbon-bearing phases in the earth have been determined from phase equilibrium measurements, combined with new data on the equations of state of water, carbon dioxide, carbonates and hydrates. The data have then been used to predict the fate of calcite and hydrous phases in subducting oceanic lithosphere. From the compositions of MORB's one can estimate concentrations of water and carbon of around 200 ppm and 80 ppm respectively in the upper mantle. Lower mantle estimates are very uncertain, but 1900 ppm water and 2000 ppm C are plausible concentrations. Measurements of the density of supercritical water to 3 GPa demonstrate that this phase is less compressible than anticipated from the equations of state of Haar et al. or Saul and Wagner and is closer to predictions based on molecular dynamics simulations. Conversely, fugacity measurements on carbon dioxide to 7 GPa show that this fluid is more compressible than predicted from the MRK equation of state. The results imply that hydrates are relatively more stable and carbonates less stable at pressures greater than 5 GPa than would be predicted from simple extrapolation of the low pressure data. Nevertheless, carbonates remain extremely refractory phases within both the upper and lower mantle.

Wood, Bernard J.↗
NASA NTRS1993

Thermal expansivity, bulk modulus, and melting curve of H2O-ice VII to 20 GPa

Equation of state properties of ice VII and fluid H2O at high pressures and temperatures have been studied experimentally from 6 to 20 GPa and 300-700 K. The techniques involve direct measurements of the unit-cell volume of the solid using synchrotron X-ray diffraction with an externally heated diamond-anvil cell. The pressure dependencies of the volume and bulk modulus of ice VII at room temperature are in good agreement with previous synchrotron X-ray studies. The thermal expansivity was determined as a function of pressure and the results fit to a newly proposed phenomenological relation and to a Mie-Gruneisen equation of state formalism. The onset of melting of ice VII was determined directly by X-ray diffraction at a series of pressures and found to be in accord with previous volumetric determinations. Thermodynamic calculations based on the new data are performed to evaluate the range of validity of previously proposed equations of state for fluid water derived from static and shock-wave compression experiments and from simulations.

Fei, Yingwei↗
NASA NTRS1994

Merging neutron stars. 1. Initial results for coalescence of noncorotating systems

We present three-dimensional Newtonian simulations of the coalescence of two neutron stars, using a smoothed particle hydrodynamics (SPH) code. We begin the simulations with the two stars in a hard, circular binary, and have them spiral together as angular momentum is lost through gravitational radiation at the rate predicted by modeling the system as two point masses. We model the neutron stars as hard polytropes (gamma = 2.4) of equal mass, and investigate the effect of the initial spin of the two stars on the coalescence. The process of coalescence, from initial contact to the formation of an axially symmetric object, takes only a few orbital periods. Some of the material from the two neutron stars is shed, forming a thick disk around the central, coalesced object. The mass of this disk is dependent on the initial neutron star spins; higher spin rates result in greater mass loss and thus more massive disks. For spin rates that are most likely to be applicable to real systems, the central coalesced object has a mass of 2.4 solar mass, which is tantalizingly close to the maximum mass allowed by any neutron star equation of state for an object that is supported in part by rotation. Using a realistic nuclear equation of state, we estimate the temperature of the material after the coalescence. We find that the central object is at a temperature of approximately 10 MeV, while the disk is heated by shocks to a temperature of 2 to 4 MeV.

Davies, M. B.↗
NASA NTRS1988

An upwind-differencing scheme for the incompressible Navier-Stokes equations

The steady state incompressible Navier-Stokes equations in 2-D are solved numerically using the artificial compressibility formulation. The convective terms are upwind-differenced using a flux difference split approach that has uniformly high accuracy throughout the interior grid points. The viscous fluxes are differenced using second order accurate central differences. The numerical system of equations is solved using an implicit line relaxation scheme. Although the current study is limited to steady state problems, it is shown that this entire formulation can be used for solving unsteady problems. Characteristic boundary conditions are formulated and used in the solution procedure. The overall scheme is capable of being run at extremely large pseudotime steps, leading to fast convergence. Three test cases are presented to demonstrate the accuracy and robustness of the code. These are the flow in a square-driven cavity, flow over a backward facing step, and flow around a 2-D circular cylinder.

Rogers, Stuart E.↗
NASA NTRS1987

Constraints on the mass-radius relation for the neutron star in the X-ray burst source 4U/MXB 1820-30 located in the globular cluster NGC 6624

From the spectral evolution of seven X-ray bursts observed from 4U/MXB 1820-30 (Haberl et al., 1986) during the phases of photospheric expansion and cooling, constraints are derived on the mass-radius relation of the neutron star in this system. Taking into account the uncertainty in the source distance and the possible anisotropy of the burst emission. If during bursts with radius expansion the outer layers of the neutron star are of a hydrogen-poor composition, none of the current equations of state of high-density matter is excluded. For a hydrogen-rich composition, only very soft equations of state are allowed.

Van Paradijs, J.↗
NASA NTRS1972

The thermodynamic properties of oxygen and nitrogen

The development of a single equation of state for oxygen and nitrogen based on the thermodynamic properties of the gases is described. The coefficients of the equation of state were determined by simultaneous least squares fits to values of isochoric heat capacity and saturation density values used to define the criteria for phase equilibrium. Tables of data for the conditions of both gases are included.

Stewart, R. B.↗
NASA NTRS1983

Closed-form solutions for a class of optimal quadratic tracking problems

Closed-form solutions are derived for a class of tracking problems including a linear optimal regulator and a prefilter for a time-invariant plant. The solutions for the prefilter equation and state trajectory coupled by the Riccati equation are exponentially related to the stability matrix of the plant. A computational procedure is presented in recursive form when the desired output state dynamics is assumed linear and time-invariant. Several examples are given for illustration.

Turner, J. D.↗