Search NASA⌕ Search

SEARCH · Search NASA

Results for “Complete 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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

On the thermodynamic framework of generalized coupled thermoelastic-viscoplastic-damage modeling

A complete potential based framework using internal state variables is put forth for the derivation of reversible and irreversible constitutive equations. In this framework, the existence of the total (integrated) form of either the (Helmholtz) free energy or the (Gibbs) complementary free energy are assumed a priori. Two options for describing the flow and evolutionary equations are described, wherein option one (the fully coupled form) is shown to be over restrictive while the second option (the decoupled form) provides significant flexibility. As a consequence of the decoupled form, a new operator, i.e., the Compliance operator, is defined which provides a link between the assumed Gibb's and complementary dissipation potential and ensures a number of desirable numerical features, for example the symmetry of the resulting consistent tangent stiffness matrix. An important conclusion reached, is that although many theories in the literature do not conform to the general potential framework outlined, it is still possible in some cases, by slight modifications of the used forms, to restore the complete potential structure.

Arnold, S. M.↗
NASA NTRS1975

Steady-state solutions of a diffusive energy-balance climate model and their stability

A diffusive energy-balance climate model, governed by a nonlinear parabolic partial differential equation, was studied. Three positive steady-state solutions of this equation are found; they correspond to three possible climates of our planet: an interglacial (nearly identical to the present climate), a glacial, and a completely ice-covered earth. Models similar to the main one are considered, and the number of their steady states was determined. All the models have albedo continuously varying with latitude and temperature, and entirely diffusive horizontal heat transfer. The stability under small perturbations of the main model's climates was investigated. A stability criterion is derived, and its application shows that the present climate and the deep freeze are stable, whereas the model's glacial is unstable. The dependence was examined of the number of steady states and of their stability on the average solar radiation.

Ghil, M.↗
NASA NTRS1976

Climate stability for a Sellers-type model

We study a diffusive energy-balance climate model governed by a nonlinear parabolic partial differential equation. Three positive steady-state solutions of this equation are found; they correspond to three possible climates of our planet: an interglacial (nearly identical to the present climate), a glacial, and a completely ice-covered earth. We consider also models similar to the main one studied, and determine the number of their steady states. All the models have albedo continuously varying with latitude and temperature, and entirely diffusive horizontal heat transfer. The diffusion is taken to be nonlinear as well as linear. We investigate the stability under small perturbations of the main model's climates. A stability criterion is derived, and its application shows that the 'present climate' and the 'deep freeze' are stable, whereas the model's glacial is unstable. A variational principle is introduced to confirm the results of this stability analysis. For a sufficient decrease in solar radiation (about 2%) the glacial and interglacial solutions disappear, leaving the ice-covered earth as the only possible climate.

Ghil, M.↗
NASA NTRS1981

Baroclinic instability of a rotating Hadley cell

The stability of a thin fluid layer between two rotating plates which are subjected to a horizontal temperature gradient is investigated. The solution for the stationary basic state is obtained in a closed form. It is pointed out that this solution identifies Ekman and thermal layers adjacent to the plates and interior temperature and velocity fields which are almost linear functions of height. The stability of that basic state with respect to infinitesimal zonal waves is then analyzed via the solution of the complete viscous linear equations for the perturbations. The character of the growth rates is found to be similar to those of the classical baroclinic waves. It is also found that the region of stability depends on the Prandtl number, the vertical stratification parameter, and both the meridional and zonal wavenumbers. The flow is generally unstable for small enough Ekman numbers and for Rossby numbers less than 10.

Antar, B. N.↗
NASA NTRS1995

Computing supernova collapse to neutron stars and black holes

We present a new numerical code for spherical hydrodynamics in general relativity. The code can handle gravitational collapse to a neutron star or to a black hole without the appearance of singularities. Moreover, the variables and equations in the code are very similar to those appearing in traditional Lagrangian supernova codes. Any such existing code can thus be easily adapted to treat collapse where the final fate is uncertain and may be either a neutron star or a black hole. The code is based on the formulation of Hernandez & Misner, in which retarded time is used as coordinate. This prevents the computational grid from penetrating inside any black hole that may form. We present the equations and a complete finite difference scheme for the adiabatic evolution of a fluid that obeys a gamma-law equation of state. We summarize the results of several testbed calculations performed to check our code. We also give the transformation of the analytic Oppenheimer-Snyder solution for homogeneous dust collapse to our coordinate system.

Baumgarte, Thomas W.↗
DOE OSTI2020

Spin Transitions and Compressibility of ε-Fe 7 N 3 and γ'-Fe 4 N: Implications for Iron Alloys in Terrestrial Planet Cores

Iron nitrides are possible constituents of the cores of Earth and other terrestrial planets. Pressure-induced magnetic changes in iron nitrides and effects on compressibility remain poorly understood. In this work, we report synchrotron X-ray emission spectroscopy (XES) and X-ray diffraction (XRD) results for ε-Fe 7 N 3 and γ'-Fe 4 N up to 60 GPa at 300 K. The XES spectra reveal completion of high- to low-spin transition in ε-Fe 7 N 3 and γ'-Fe 4 N at 43 and 34 GPa, respectively. The completion of the spin transition induces stiffening in bulk modulus of ε-Fe 7 N 3 by 22% at ~40 GPa, but has no resolvable effect on the compression behavior of γ'-Fe 4 N. Fitting pressure-volume data to the Birch-Murnaghan equation of state yields V 0 = 83.29 ± 0.03 (Å 3 ), K 0 = 232 ± 9 GPa, K 0 ' = 4.1 ± 0.5 for nonmagnetic ε-Fe 7 N 3 above the spin transition completion pressure, and V 0 = 54.82 ± 0.02 (Å 3 ), K 0 = 152 ± 2 GPa, K 0 ' = 4.0 ± 0.1 for γ'-Fe 4 N over the studied pressure range. By reexamining evidence for spin transition and effects on compressibility of other candidate components of terrestrial planet cores, Fe 3 S, Fe 3 P, Fe 7 C 3 , and Fe 3 C based on previous XES and XRD measurements, we located the completion of high- to low-spin transition at ~67, 38, 50, and 30 GPa at 300 K, respectively. The completion of spin transitions of Fe 3 S, Fe 3 P, and Fe 3 C induces elastic stiffening, whereas that of Fe 7 C 3 induces elastic softening. Changes in compressibility at completion of spin transitions in iron-light element alloys may influence the properties of Earth's and planetary cores.

58 GEOSCIENCES↗
DOE OSTI2022

Benchmark calculations of infinite neutron matter with realistic two- and three-nucleon potentials

In this study, we present the equation of state of infinite neutron matter as obtained from highly realistic Hamiltonians that include nucleon-nucleon and three-nucleon coordinate-space potentials. We benchmark three independent many-body methods: Brueckner-Bethe-Goldstone (BBG), Fermi hypernetted chain/single-operator chain (FHNC/SOC), and auxiliary-field diffusion Monte Carlo (AFDMC). We find them to provide similar equations of state when the Argonne v 18 and the Argonne v' 6 nucleon-nucleon potentials are used in combination with the Urbana IX three-body force. Only at densities larger than about 1.5 the nuclear saturation density ( ρ 0 =0.16 fm ⁻3 ) the FHNC/SOC energies are appreciably lower than the other two approaches. The AFDMC calculations carried out with all of the Norfolk potentials fitted to reproduce the experimental trinucleon ground-state energies and nd doublet scattering length yield unphysically bound neutron matter, associated with the formation of neutron droplets. Including tritium β decay in the fitting procedure, as in the second family of Norfolk potentials, mitigates but does not completely resolve this problem. An excellent agreement between the BBG and AFDMC results is found for the subset of Norfolk interactions that do not make neutron-matter collapse, while the FHNC/SOC equations of state are moderately softer.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1984

Kalman filter modeling

The formulation of appropriate state-space models for Kalman filtering applications is studied. The so-called model is completely specified by four matrix parameters and the initial conditions of the recursive equations. Once these are determined, the die is cast, and the way in which the measurements are weighted is determined foreverafter. Thus, finding a model that fits the physical situation at hand is all important. Also, it is often the most difficult aspect of designing a Kalman filter. Formulation of discrete state models from the spectral density and ARMA random process descriptions is discussed. Finally, it is pointed out that many common processes encountered in applied work (such as band-limited white noise) simply do not lend themselves very well to Kalman filter modeling.

Brown, R. G.↗
DOE OSTI2024

Testing α -attractor quintessential inflation against CMB and low-redshift data

Due to universality and attractor properties, α-attractor quintessential inflation establishes direct relations between inflationary observables such as the scalar tilt n s and the tensor-to-scalar ratio r, and late-time dark energy equation of state parameters w 0 and w a . In this work, we examine three different physically motivated regimes, considering complete freedom in the parameter α, models inspired by supergravity where α takes on values up to α = 7/3, and Starobinsky inflation (α = 1). We investigate the consistency and constraints imposed by Cosmic Microwave Background measurements from the Planck satellite, B-mode polarization data from the BICEP/Keck collaboration, and low-redshift observations. Additionally, we consider small-scale CMB measurements released by the Atacama Cosmology Telescope, which give results approaching the Harrison– Zel’dovich spectrum (n s ≈ 1). Here α-attractors lead to an improved fit over $\Lambda$CDM. For the large-scale CMB measurements, α ≳ 2 models can provide equally good fits as $\Lambda$CDM

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1976

Stability of laser heated flows

A local stability analysis is utilized to determine the stability of disturbances generated at each point along a nozzle of variable area ratio for a one-dimensional flow heated by laser radiation entering from the upstream direction. The governing equations for the quasi-one-dimensional flow without viscous dissipation, diffusion, and thermal conduction but including radiative heat transfer are given. The governing equations are combined to yield a relationship which governs the Mach number variation through the nozzle. The complete steady-state solution can be calculated from knowledge of the Mach number profile, the inlet conditions, and the laser power. The local stability analysis permits obtaining contour (or contours) of neutral stability. Solutions have been obtained for various nozzle configurations, but only one set of example calculations is presented. The results obtained indicate that the analysis serves as an important indicator as to where potential absorption wave phenomena may be initiated.

Wu, P. K. S.↗
NASA NTRS1981

Theory of dense hydrogen

Equations of state for molecular and metallic hydrogen are calculated in an accurate and comparable manner by using the density-functional method applied to static lattices. By isolating the electronic energies in this way, attention can be drawn to the relative importance of the protonic degrees of freedom. From the equation-of-state results it can be concluded that a remnant molecular pairing is preferred in a band-overlap metallic state, complete dissociation occurring only at very high densities. An accurate determination of the pressure needed to achieve complete dissociation is shown to require a self-consistent treatment of electron and proton degrees of freedom.

Chakravarty, S.↗
NASA NTRS1993

Trajectory optimization based on differential inclusion

A method for generating finite-dimensional approximations to the solutions of optimal control problems is introduced. By employing a description of the dynamical system in terms of its attainable sets in favor of using differential equations, the controls are completely eliminated from the system model. Besides reducing the dimensionality of the discretized problem compared to state-of-the-art collocation methods, this approach also alleviates the search for initial guesses from where standard gradient search methods are able to converge. The mechanics of the new method are illustrated on a simple double integrator problem. The performance of the new algorithm is demonstrated on a 1-D rocket ascent problem ('Goddard Problem') in presence of a dynamic pressure constraint.

Seywald, Hans↗
NASA NTRS1993

Trajectory optimization based on differential inclusion

A method for generating finite-dimensional approximations to the solutions of optimal control problems is introduced. By employing a description of the dynamical system in terms of its attainable sets in favor of using differential equations, the controls are completely eliminated from the system model. Besides reducing the dimensionality of the discretized problem compared to state-of-the-art collocation methods, this approach also alleviates the search for initial guesses from where standard gradient search methods are able to converge. The mechanics of the new method are illustrated on a simple double integrator problem. The performance of the new algorithm is demonstrated on a I-D rocket ascent problem ('Goddard Problem') in presence of a dynamic pressure constraint.

Seywald, Hans↗
NASA NTRS1993

Trajectory optimization based on differential inclusion

A method for generating finite-dimensional approximations to the solutions of optimal control problems is introduced. By employing a description of the dynamical system in terms of its attainable sets in favor of using differential equations, the controls are completely eliminated from the system model. Besides reducing the dimensionality of the discretized problem compared to state-of-the-art collocation methods, this approach also alleviates the search for initial guesses from where standard gradient search methods are able to converge. The mechanics of the new method are illustrated on a simple double integrator problem. The performance of the new algorithm is demonstrated on a 1-D rocket ascent problem (`Goddard Problem') in presence of a dynamic pressure constraint.

Seywald, Hans↗
DOE OSTI2020

Derivation of K-matrix reaction theory in a discrete basis formalism

The usual derivations of the S and K matrices for two-particle reactions proceed through the Lippmann–Schwinger equation with formal definitions of the incoming and outgoing scattering states. Here we present a simpler alternative derivation that is carried out completely in the Hamiltonian representation, using a discrete basis of configurations for the scattering channels as well as the quasi-bound configurations of the combined fragments. We use matrix algebra to derive an explicit expression for the K matrix in terms of the Hamiltonian of the internal states of the compound system and the coupling between the channels and the internal states. The formula for the K matrix includes explicitly a real dispersive shift matrix to the internal Hamiltonian that is easily computed in the formalism. That expression is applied to derive the usual form of the S matrix as a sum over poles in the complex energy plane. Some extensions and limitations of the discrete-basis Hamiltonian formalism are discussed in the concluding remarks and in the Appendix.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1983

A theory of nondilute spray evaporation based upon multiple drop interactions

Drop interactions, changes in drop temperature, and the evolution of the ambient conditions due to heating and evaporation were taken into account in a theory of drop evaporation. A sphere of influence is defined for each drop, with all droplets having equal separation distances. Global conservation and state equations for a droplet spray are formulated, together with equations for single drop evaporation in finite surroundings. A criterion is defined for complete evaporation and the dynamics of spray evaporation are discussed. The theory is noted to be applicable for both dilute and nondilute spray evaporation effects.

Bellan, J.↗
NASA NTRS1990

Closed-form solution of decomposable stochastic models

Markov and semi-Markov processes are increasingly being used in the modeling of complex reconfigurable systems (fault tolerant computers). The estimation of the reliability (or some measure of performance) of the system reduces to solving the process for its state probabilities. Such a model may exhibit numerous states and complicated transition distributions, contributing to an expensive and numerically delicate solution procedure. Thus, when a system exhibits a decomposition property, either structurally (autonomous subsystems), or behaviorally (component failure versus reconfiguration), it is desirable to exploit this decomposition in the reliability calculation. In interesting cases there can be failure states which arise from non-failure states of the subsystems. Equations are presented which allow the computation of failure probabilities of the total (combined) model without requiring a complete solution of the combined model. This material is presented within the context of closed-form functional representation of probabilities as utilized in the Symbolic Hierarchical Automated Reliability and Performance Evaluator (SHARPE) tool. The techniques adopted enable one to compute such probability functions for a much wider class of systems at a reduced computational cost. Several examples show how the method is used, especially in enhancing the versatility of the SHARPE tool.

Sjogren, Jon A.↗
NASA NTRS1984

Use of the Burton-prim-slichter Equation at High Growth Rates

The following must be assumed in order to derive the classic Burton-Prim-Slichter equation for segregation during crystal growth: isobaric, isothermal, isopotential or uncharged species, binary, planar interface, steady state, constant diffusion coefficient in fluid, no diffusion in crystal, no lateral convection within fluid film at interface with complete mixing beyond (stagnant film model), and either density, total concentration or partial molar volumes constant in the fluid phase. In addition, the effective distribution coefficient and the interfacial distribution coefficient must be defined properly. The velocity in the equation is the growth rate times a factor correcting for the difference in volumetric properties between crystal and fluid. The stagnant film thickness is found to be a function of freezing rate, with the precise functionality depending on the type of stirring occurring in the growth fluid.

Wilcox, W. R.↗