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

The Principle of Energetic Consistency: Application to the Shallow-Water Equations

If the complete state of the earth's atmosphere (e.g., pressure, temperature, winds and humidity, everywhere throughout the atmosphere) were known at any particular initial time, then solving the equations that govern the dynamical behavior of the atmosphere would give the complete state at all subsequent times. Part of the difficulty of weather prediction is that the governing equations can only be solved approximately, which is what weather prediction models do. But weather forecasts would still be far from perfect even if the equations could be solved exactly, because the atmospheric state is not and cannot be known completely at any initial forecast time. Rather, the initial state for a weather forecast can only be estimated from incomplete observations taken near the initial time, through a process known as data assimilation. Weather prediction models carry out their computations on a grid of points covering the earth's atmosphere. The formulation of these models is guided by a mathematical convergence theory which guarantees that, given the exact initial state, the model solution approaches the exact solution of the governing equations as the computational grid is made more fine. For the data assimilation process, however, there does not yet exist a convergence theory. This book chapter represents an effort to begin establishing a convergence theory for data assimilation methods. The main result, which is called the principle of energetic consistency, provides a necessary condition that a convergent method must satisfy. Current methods violate this principle, as shown in earlier work of the author, and therefore are not convergent. The principle is illustrated by showing how to apply it as a simple test of convergence for proposed methods.

Cohn, Stephen E.↗

Airfoil Design Using a Coupled Euler and Integral Boundary Layer Method with Adjoint Based Sensitivities

The objective of this paper is to present a control theory approach for the design of airfoils in the presence of viscous compressible flows. A coupled system of the integral boundary layer and the Euler equations is solved to provide rapid flow simulations. An adjunct approach consistent with the complete coupled state equations is employed to obtain the sensitivities needed to drive a numerical optimization algorithm. Design to target pressure distribution is demonstrated on an RAE 2822 airfoil at transonic speed.

Edwards, S.↗

Impact jetting of water ice, with application to the accretion of icy planetesimals and Pluto

Jetting can occur during oblique impacts of water-ice bodies at relative velocities as low as 500 m/sec, because of the low Hugoniot elastic limit and high compressibility of ice compared to rock. In jetted ice, incipient melting, complete melting, and incipient vaporization occur, upon release to low pressure, at impact velocities of 1.3, 2.0, and 2.7 km/sec, respectively, much less than the 3.4, 4.4, and 5.3 km/sec, required in head-on collisions. Uncertainties in the shock equation-of-state may allow complete melting during jetting of relative velocities as low as 1.2 km/sec. Because jet speeds exceed impact speeds during the accretion of icy bodies greater than a few 100 km in radius, there may be a significant loss of icy material. Thus, jetting during a Charon-forming collision (and not vaporization) may account for Pluto-Charon's relatively large rock/ice ratio, should the C/O ratio of the solar nebula turn out to be too low to sufficiently raise the rock/ice ratio of outer solar nebula condensates by formation of noncondensable CO.

Mckinnon, William B.↗

Sound generation by a stenosis in a pipe

This paper presents a computational aeroacoustics study of sound generated by low Mach number flow through a closure (stenosis) in a circular pipe. The time-dependent incompressible flowfield in the pipe is first computed utilizing a vorticity/stream function formulation. The known velocity field is then utilized to determine the incompressible fluctuating pressure in the pipe. This hydrodynamic pressure field is then utilized to compute a hydrodynamic density perturbation to the constant incompressible density through the equation of state. Knowledge of this complete hydrodynamic field is then employed as the source of the resulting sound radiation. This tripartite technique allows separation of the compressible fluctuations from the purely hydrodynamic fluctuations. Results of the analysis are compared with experimental measurements of sound radiated by such a flow.

Hardin, J. C.↗

Momentum transfer in the one-dimensional impact of spaced plates

The two-dimensional hydrocode CSQ III was used to calculate the fraction of momentum transferred from a flyer plate to a target of two spaced plates. The effect of the vaporization phase transition, as calculated with the ANEOS analytical complete three-phase equation of state, was estimated. Application of these results to the protection of spacecraft from meteoroids and orbital debris is discussed.

Abbott, R. D.↗

Development of OPTIMOD - The new Space Shuttle Orbiter/Payload Thermal Integration Model

The development of the Orbiter/Payload Thermal Integration Model (OPTIMOD), a simplified Space Shuttle Orbiter thermal math model to be used for payload thermal integration analysis, is described. The Orbiter lower fuselage math model has been reduced to the minimum level of detail consistent with accuracy requirements. Instead, an 'environment radiation sensor', a simple closed-end cylinder oriented in a fixed position with respect to the Orbiter, is introduced. Equations are presented which relate the cylinder surface heating rates to external fuselage steady-state temperatures. These equations were derived from a regression correlation analysis of steady-state Orbiter thermal solutions for 2372 points on orbit, representing a complete range of vehicle-earth-sun orientations. An exponential equation is developed to relate the steady-state data to time-variant temperatures for the particular transient attitude timeline being simulated. The accuracy of OPTIMOD is shown to be excellent in the payload bay and to compare favorably with previous models in Orbiter elements outside the bay.

Russell, David J.↗

Dynamic simulation of coronal mass ejections

A model is developed for the formation and propagation through the lower corona of the loop-like coronal transients in which mass is ejected from near the solar surface to the outer corona. It is assumed that the initial state for the transient is a coronal streamer. The initial state for the streamer is a polytropic, hydrodynamic solution to the steady-state radial equation of motion coupled with a force-free dipole magnetic field. The numerical solution of the complete time-dependent equations then gradually approaches a stationary coronal streamer configuration. The streamer configuration becomes the initial state for the coronal transient. The streamer and transient simulations are performed completely independent of each other. The transient is created by a sudden increase in the pressure at the base of the closed-field region in the streamer configuration. Both coronal streamers and coronal transients are calculated for values of the plasma beta (the ratio of thermal to magnetic pressure) varying from 0.1 to 100.

Steinolfson, R. S.↗

Performance analysis of flexible aircraft with active control

Small perturbation equations of motion of a flexible aircraft with an active control technology (ACT) system were developed to evaluate the stability and performance of the controlled aircraft. The total aircraft system was formulated in state vector format and the system of equations was completed with fully unsteady and low frequency aerodynamics for arbitrary, complex configurations based on a potential aerodynamic method. The ACT system equations were incorporated in the digital computer program FCAP (Flight Control Analysis Program) which can be used for the analysis of complete aircraft configurations, including control system, with either low frequency or fully unsteady aerodynamics. The application of classical performance analyses including frequency response, poles and zeros, mean square response, and time response in FCAP in state vector format was discussed.

Noll, R. B.↗

Response analysis of flexible aircraft with active control

The small-perturbation equations of motion of a flexible aircraft with an active control technology (ACT) system are developed to evaluate the stability and performance of the controlled aircraft. The total aircraft system is formulated in state vector format and the system of equations is completed with fully unsteady and low-frequency aerodynamics for arbitrary, complex configurations based on a potential aerodynamic method. The ACT system equations are incorporated in the digital computer program FCAP (Flight Control Analysis Program) which can be used for the analysis of complete aircraft configurations, including control system, with either low-frequency or fully unsteady aerodynamics. The application of classical performance analyses including frequency response, poles and zeros, mean-square response, and time response in FCAP in state vector format is discussed.

Noll, R. B.↗

Equation-of-state spinning fluids in the Einstein-Cartan theory

The relativistic fluid equations may be completed in two physically distinct methods. One method assumes the mass, rho, (or particle number) is conserved, while the other method assumes an equation of state of the form P = P(rho). A variational principle for the mass conservation method both with and without an intrinsic spin for the fluid was constructed earlier (Ray and Smalley, 1982 and 1983). A variational principle for the fluid described by an equation of state both with and without spin is formulated. In all cases the variational principle is set in the Einstein-Cartan metric-torsion U4 geometry. The results for general relativity follow as a special case.

Ray, John R.↗

Semidirect calculation of steady two- and three-dimensional flows

This paper describes a semidirect method for rapidly solving steady-state viscous flows described by the complete Navier-Stokes equations. The current results are for two-dimensional incompressible flows in general channels at arbitrary Reynolds numbers, but work in progress on compressible and three-dimensional flows is also described. The basic concept of semidirect methods is to use the recently developed fast (direct, or noniterative) linear solvers to solve linearized equations, which are then iterated to solve the nonlinearity. The method used here is an extension of the Split NOS method (Roache, 1975)

Roache, P. J.↗

Hypersonic shock structure with Burnett terms in the viscous stress and heat flux

The continuum Navier-Stokes and Burnett equations are solved for one-dimensional shock structure in various monatomic gases. A new numerical method is employed which utilizes the complete time-dependent continuum equations and obtains the steady-state shock structure by allowing the system to relax from arbitrary initial conditions. Included is discussion of numerical difficulties encountered when solving the Burnett equations. Continuum solutions are compared to those obtained utilizing the Direct Simulation Monte Carlo method. Shock solutions are obtained for a hard sphere gas and for argon from Mach 1.3 to Mach 50. Solutions for a Maxwellian gas are obtained from Mach 1.3 to Mach 3.8. It is shown that the Burnett equations yield shock structure solutions in much closer agreement to both Monte Carlo and experimental results than do the Navier-Stokes equations. Shock density thickness, density asymmetry, and density-temperature separation are all more accurately predicted by the Burnett equations than by the Navier-Stokes equations.

Chapman, Dean R.↗

Stochastic time-optimal control problems

Two types of stochastic time-optimal controls in a one-dimensional setting are considered. Multidimensional problems, in the case of complete state information available and the system modeled by stochastic differential equations, are studied under the formulation of minimizing the expected transient-response time. The necessary condition of optimality is the satisfaction for the value function of a parabolic partial differential equation with boundary conditions. The sufficient condition of optimality is also provided, based on Dynkin's formula. Finally, three examples are given.

Zhang, W.↗

Comparison of Burnett, super-Burnett and Monte Carlo solutions for hypersonic shock structure

The continuum Navier-Stokes, Burnett, and super-Burnett equations are solved for one-dimensional shock structures in various monoatomic gases. Solutions for a hard sphere gas, argon, and a Maxwellian gas from Mach 1.3 to Mach 50 are obtained. A numerical method utilizing the complete time-dependent continuum equations and obtaining the steady-state shock structure by allowing the system to relax from arbitrary initial conditions is employed. Shock density, velocity, temperature, and entropy profiles are also obtained using the direct simulation Monte Carlo method, and these results are used as bases for comparison for continuum solution profiles. It is shown that the Burnett equations yield shock structure solutions in much closer agreement to both Monte Carlo and experimental results than do the Navier-Stokes equations. Solutions to the super-Burnett equations with coefficients as presently derived, however, are considered to be inferior to those of the Burnett equations.

Fiscko, Kurt A.↗

Transient and steady-state velocity of domain walls for a complete range of drive fields

Approximate analytic solutions for transient and steady-state 180 deg domain wall motion in bulk magnetic material are obtained from the dynamic torque equations with a Gilbert damping term. The results for the Walker region in which the transient solution approaches the familiar Walker steady-state solution are presented in a slightly new form for completeness. An analytic solution corresponding to larger drive fields predicts an oscillatory motion with an average value which decreases with drive field for reasonable values of the damping parameter. These results agree with those obtained by a computer solution of the torque equation and those obtained with the assumption of a very large anisotropy field.

Bourne, H. C., Jr.↗

Transient and steady-state velocity of domain walls for a complete range of drive fields

Approximate analytic solutions for transient and steady-state 180 deg domain-wall motion in bulk magnetic material are obtained from the dynamic torque equations with a Gilbert damping term. The results for the Walker region in which the transient solution approaches the familiar Walker steady-state solution are presented in a slightly new form for completeness. An analytic solution corresponding to larger drive fields predicts an oscillatory motion with an average value of the velocity which decreases with drive field for reasonable values of the damping parameter. These results agree with those obtained by others from a computer solution of the torque equation and those obtained by others with the assumption of a very large anisotropy field.

Bourne, H. C., Jr.↗

Fully Associative, Nonisothermal, Potential-Based Unified Viscoplastic Model for Titanium-Based Matrices

A number of titanium matrix composite (TMC) systems are currently being investigated for high-temperature air frame and propulsion system applications. As a result, numerous computational methodologies for predicting both deformation and life for this class of materials are under development. An integral part of these methodologies is an accurate and computationally efficient constitutive model for the metallic matrix constituent. Furthermore, because these systems are designed to operate at elevated temperatures, the required constitutive models must account for both time-dependent and time-independent deformations. To accomplish this, the NASA Lewis Research Center is employing a recently developed, complete, potential-based framework. This framework, which utilizes internal state variables, was put forth for the derivation of reversible and irreversible constitutive equations. The framework, and consequently the resulting constitutive model, is termed complete because the existence of the total (integrated) form of the Gibbs complementary free energy and complementary dissipation potentials are assumed a priori. The specific forms selected here for both the Gibbs and complementary dissipation potentials result in a fully associative, multiaxial, nonisothermal, unified viscoplastic model with nonlinear kinematic hardening. This model constitutes one of many models in the Generalized Viscoplasticity with Potential Structure (GVIPS) class of inelastic constitutive equations.

Source record↗

Relative Shock Effects in Mixed Powders of Calcite, Gypsum, and Quartz: A Calibration Scheme from Shock Experiments

The shock behavior of calcite and gypsum is important in understanding the Cretaceous/Tertiary event and other terrestrial impacts that contain evaporite sediments in their targets. Most interest focuses on issues of devolatilization to quantify the production of CO2 or SO2 to better understand their role in generating a temporary atmosphere and its effects on climate and biota [e.g., papers in 1,2,3,4]. Devolatilization of carbonate is also important because the dispersion and fragmentation of ejecta is strongly controlled by the expansion of large volumes of gas during the impact process as well [5,6]. Shock recovery experiments for calcite yield seemingly conflicting results: early experimental devolatilization studies [7,8,9] suggested that calcite was substantially outgassed at 30 GPa (> 50%). However, the recent petrographic work of [10,11,12] presented evidence that essentially intact calcite is recovered from 60 GPa experiments. [13] reported results of shock experiments on anhydrite, gypsum, and mixtures of those phases with silica. Their observations indicate 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 (preliminary) report of shock experiments on calcite, anhydrite, and gypsum, [14] observe calcite recrystallization when shock loaded at 61 GPa, only intensive plastic deformation in anhydrite shock loaded at 63 GPa, and gypsum converted to anhydrite when shock loaded at 56 GPa. [15] shock loaded anhydrite and quartz to a peak pressure of 60 GPa. All of the quartz grains were trans-formed to glass and 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 suggest that recrystallization of anhydrite grains is the result of a solid state transformation. [16] reanalyzed the calcite and anhydrite shock wave experiments of [17] 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 and 17.8 and 54.1 GPa for calcite, respectively, a factor of 2 to 3 lower than reported earlier by [17].

Bell, Mary S.↗