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 145 records · Page 8

Comparative study of high-resolution shock-capturing schemes for a real gas

Recently developed second-order explicit shock-capturing methods, in conjunction with generalized flux-vector splittings, and a generalized approximate Riemann solver for a real gas are studied. The comparisons are made on different one-dimensional Riemann (shock-tube) problems for equilibrium air with various ranges of Mach numbers, densities and pressures. Six different Riemann problems are considered. These tests provide a check on the validity of the generalized formulas, since theoretical prediction of their properties appears to be difficult because of the non-analytical form of the state equation. The numerical results in the supersonic and low-hypersonic regimes indicate that these produce good shock-capturing capability and that the shock resolution is only slightly affected by the state equation of equilibrium air. The difference in shock resolution between the various methods varies slightly from one Riemann problem to the other, but the overall accuracy is very similar. For the one-dimensional case, the relative efficiency in terms of operation count for the different methods is within 30%. The main difference between the methods lies in their versatility in being extended to multidimensional problems with efficient implicit solution procedures.

Montagne, J.-L.↗

Comparative study of high-resolution shock-capturing schemes for a real gas

Recently developed second-order explicit shock-capturing methods, in conjunction with generalized flux-vector splittings, and a generalized approximate Riemann solver for a real gas are studied. The comparisons are made on different one-dimensional Riemann (shock-tube) problems for equilibrium air with various ranges of Mach numbers, densities and pressures. Six different Riemann problems are considered. These tests provide a check on the validity of the generalized formulas, since theoretical prediction of their properties appears to be difficult because of the non-analytical form of the state equation. The numerical results in the supersonic and low-hypersonic regimes indicate that these produce good shock-capturing capability and that the shock resolution is only slightly affected by the state equation of equilibrium air. The difference in shock resolution between the various methods varies slightly from one Riemann problem to the other, but the overall accuracy is very similar. For the one-dimensional case, the relative efficiency in terms of operation count for the different methods is within 30 percent. The main difference between the methods lies in their versatility in being extended to multidimensional problems with efficient implicit solution procedures.

Montagne, J.-L.↗

Automatic determination of fault effects on aircraft functionality

The problem of determining the behavior of physical systems subsequent to the occurrence of malfunctions is discussed. It is established that while it was reasonable to assume that the most important fault behavior modes of primitive components and simple subsystems could be known and predicted, interactions within composite systems reached levels of complexity that precluded the use of traditional rule-based expert system techniques. Reasoning from first principles, i.e., on the basis of causal models of the physical system, was required. The first question that arises is, of course, how the causal information required for such reasoning should be represented. The bond graphs presented here occupy a position intermediate between qualitative and quantitative models, allowing the automatic derivation of Kuipers-like qualitative constraint models as well as state equations. Their most salient feature, however, is that entities corresponding to components and interactions in the physical system are explicitly represented in the bond graph model, thus permitting systematic model updates to reflect malfunctions. Researchers show how this is done, as well as presenting a number of techniques for obtaining qualitative information from the state equations derivable from bond graph models. One insight is the fact that one of the most important advantages of the bond graph ontology is the highly systematic approach to model construction it imposes on the modeler, who is forced to classify the relevant physical entities into a small number of categories, and to look for two highly specific types of interactions among them. The systematic nature of bond graph model construction facilitates the process to the point where the guidelines are sufficiently specific to be followed by modelers who are not domain experts. As a result, models of a given system constructed by different modelers will have extensive similarities. Researchers conclude by pointing out that the ease of updating bond graph models to reflect malfunctions is a manifestation of the systematic nature of bond graph construction, and the regularity of the relationship between bond graph models and physical reality.

Feyock, Stefan↗

Comparison of some optimal control methods for the design of turbine blades

This paper attempts a comparative study of some numerical methods for the optimal control design of turbine blades whose vibration characteristics are approximated by Timoshenko beam idealizations with shear and incorporating simple boundary conditions. The blade was synthesized using the following methods: (1) conjugate gradient minimization of the system Hamiltonian in function space incorporating penalty function transformations, (2) projection operator methods in a function space which includes the frequencies of vibration and the control function, (3) epsilon-technique penalty function transformation resulting in a highly nonlinear programming problem, (4) finite difference discretization of the state equations again resulting in a nonlinear program, (5) second variation methods with complex state differential equations to include damping effects resulting in systems of inhomogeneous matrix Riccatti equations some of which are stiff, (6) quasi-linear methods based on iterative linearization of the state and adjoint equation. The paper includes a discussion of some substantial computational difficulties encountered in the implementation of these techniques together with a resume of work presently in progress using a differential dynamic programming approach.

Desilva, B. M. E.↗

Method of Characteristics Calculations and Computer Code for Materials with Arbitrary Equations of State and Using Orthogonal Polynomial Least Square Surface Fits

A numerical scheme using the method of characteristics to calculate the flow properties and pressures behind decaying shock waves for materials under hypervelocity impact is developed. Time-consuming double interpolation subroutines are replaced by a technique based on orthogonal polynomial least square surface fits. Typical calculated results are given and compared with the double interpolation results. The complete computer program is included.

Chang, T. S.↗

Aluminum under high pressure: 1 equation of state

A curve of applied pressure P versus lattice constant a is calculated for single crystal aluminum. It results from an application of the method of structural expansions for deriving the energies of simple metals, a method known to give reasonable results for the elastic constants even at second order in the effective electron-ion interaction. The method is taken from Fermi surface analysis . This essentially experimental information is used to verify that the extant face-centered cubic structure remains the preferred crystalline phase up to the highest pressures considered. Arguments are given to suggest that the P versus a curve should have reasonable a priori accuracy, which may be improved if experimental data in the intermediate pressure region can be provided to refine the energy dependent pseudo-potential. At 3 megabars the lattice constant is reduced by only 22%; the ion cores at this pressure are still very well separated.

Friedli, C.↗

Lattice model theory of the equation of state covering the gas, liquid, and solid phases

The three stable states of matter and the corresponding phase transitions were obtained with a single model. Patterned after Lennard-Jones and Devonshires's theory, a simple cubic lattice model containing two fcc sublattices (alpha and beta) is adopted. The interatomic potential is taken to be the Lennard-Jones (6-12) potential. Employing the cluster variation method, the Weiss and the pair approximations on the lattice gas failed to give the correct phase diagrams. Hybrid approximations were devised to describe the lattice term in the free energy. A lattice vibration term corresponding to a free volume correction is included semi-phenomenologically. The combinations of the lattice part and the free volume part yield the three states and the proper phase diagrams. To determine the coexistence regions, the equalities of the pressure and Gibbs free energy per molecule of the coexisting phases were utilized. The ordered branch of the free energy gives rise to the solid phase while the disordered branch yields the gas and liquid phases. It is observed that the triple point and the critical point quantities, the phase diagrams and the coexistence regions plotted are in good agreement with the experimental values and graphs for argon.

Bonavito, N. L.↗

Aluminum under high pressure. I - Equation of state

A curve of applied pressure versus lattice constant is calculated for single-crystal aluminum. It results from an application of the method of structural expansions for deriving the energies of simple metals, a method known to give reasonable results for the elastic constants even at second order in the effective electron-ion interaction. The latter is taken from Fermi-surface analysis, and it is verified that the extant face-centered cubic structure remains the preferred crystalline phase up to the highest pressures considered. Arguments are given to suggest that the curve should have reasonable a priori accuracy and can admit possible improvement if experimental data in the intermediate-pressure region can be provided to refine the (in principle) energy-dependent pseudopotential. At three megabars, the lattice constant is reduced by only 22 per cent; the ion cores at this pressure are still very well separated.

Friedli, C.↗

On the contribution of intramolecular zero point energy to the equation of state of solid H2

Experimental evidence shows that the internal zero-point energy of the H2 molecule exhibits a relatively strong pressure dependence in the solid as well as changing considerably upon condensation. It is shown that these effects contribute about 6% to the total sublimation energy and to the pressure in the solid state. Methods to modify the ab initio isolated pair potential to account for these environmental effects are discussed.

Chandrasekharan, V.↗

Dynamic properties of mare basalts - Relation of equation of state to petrology

A comparison of shock compression and release adiabat data to 160 GPa for the 12063 low-titanium mare basalt with previous results for the high-titanium, 70215 mare basalt shows that, while the two rocks have very similar zero-pressure densities, the former is much less compressible above 100 GPa. This demonstrates that rocks with complex mineralogies can have significantly different dynamic properties despite a similarity of oxide components. The 10% higher ilmenite content of 70215 is suggested as an explanation for the contrasting dynamic properties of the two rocks. Because the 12063 high-pressure data presented lie midway between the 70215 data and those previously obtained for the anorthosite that is the dominant constituent of the lunar highlands, it is concluded that the behavior of lunar surface material under dynamic compression is terrane-dependent; where this dependence is (1) found to extend to different regions of the maria, and (2) may affect the determination of the relative ages of different rock units from areal crater densities.

Ahrens, T. J.↗

On the equation of state of solar wind ions derived from Helios measurements

The radial evolution of such adiabatic invariants as the ion magnetic moments are studied on the basis of solar wind ion velocity distribution observations made by the Helios spacecraft between 0.3 and 1 AU. Significant differences between proton and alpha-particle parameters in wind velocity dependence are noted, and adiabaticity is on the average found to be violated. This violation is interpreted as evidence that protons are heated perpendicular to the field in fast streams. The contribution of the differential streaming energy to the total internal energy of the ions is also considered, and average heliocentric radial profiles for the ion heat fluxes are presented together with an examination of a possible ion heat flux role in the supply of thermal energy during the winds radial expansion. The present findings suggest that wave-particle interactions, and/or Coulomb collisions, are required to explain solar wind ion thermal energy states and radial temperature profiles.

Marsch, E.↗