Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “equations 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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
NASA NTRS1984

Discrete approximation methods for parameter identification in delay systems

Approximation schemes for parameter identification problems in which the governing state equation is a linear functional differential equation of retarded type are constructed. The basis of the schemes is the replacement of the parameter identification problem having an infinite dimensional state equation by a sequence of approximating parameter identification problems in which the states are given by finite dimensional discrete difference equations. The difference equations are constructed using linear semigroup theory and rational function approximations to the exponential. Sufficient conditions are given for the convergence of solutions to the approximating problems, which can be obtained using conventional methods, to solutions to the original parameter identification problem. Finite difference and spline based schemes using Paderational function approximations to the exponential are constructed, and shown to satisfy the sufficient conditions for convergence. A discussion and analysis of numerical results obtained through the application of the schemes to several examples is included.

Rosen, I. G.↗
NASA NTRS1973

Aircraft Range Optimization Using Singular Perturbations

An approximate analytic solution is developed for the problem of maximizing the range of an aircraft for a fixed end state. The problem is formulated as a singular perturbation and solved by matched inner and outer asymptotic expansions and the minimum principle of Pontryagin. Cruise in the stratosphere, and on transition to and from cruise at constant Mach number are discussed. The state vector includes altitude, flight path angle, and mass. Specific fuel consumption becomes a linear function of power approximating that of the cruise values. Cruise represents the outer solution; altitude and flight path angle are constants, and only mass changes. Transitions between cruise and the specified initial and final conditions correspond to the inner solutions. The mass is constant and altitude and velocity vary. A solution is developed which is valid for cruise but which is not for the initial and final conditions. Transforming of the independent variable near the initial and final conditions result in solutions which are valid for the two inner solutions but not for cruise. The inner solutions can not be obtained without simplifying the state equations. The singular perturbation approach overcomes this difficulty. A quadratic approximation of the state equations is made. The resulting problem is solved analytically, and the two inner solutions are matched to the outer solution.

Oconnor, Joseph Taffe↗
NASA NTRS1979

Multilevel modeling of the space construction base

A mathematical model representing a flexible space vehicle in terms of interconnected discrete masses is recast into hierarchical form to facilitate the application of multilevel control techniques. The model's state equations are temporarily decoupled by utilization of a set of coordination equations. A hierarchical form of the model is constructed with the set of coordination equations at the top of the hierarchy and the set of decoupled scalar state equations in the bottom level.

Chichester, F. D.↗
DOE OSTI2022

New equations of state for saturated tuff based on the CHEST model

The CHemical Equilibrium for Saturated Tuff (CHEST) model was originally presented in Ref. 1, with numerical results depicted in Figure 1 and tabulated in Appendix B of that work. The results consisted of specific internal energy as a function of pressure and specific volume, E(P, V), spanning the domain P = 10 –4 – 10 4 kbar and V = 0.17 – 752 cm 3 /g (ρ = 1.33 × 10 –3 – 5.9 g/cm 3 ).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Reply to “Comment on ‘Reexamining the relation between the binding energy of finite nuclei and the equation of state of infinite nuclear matter' ”

In their comment to our paper [1], Bertsch and Stroberg [2] provide three criticisms. Two of these concern the interpretation of our dispersive optical model (DOM) results and their relation to the liquid drop model (LDM) parameters. The third criticism focuses on the potential systematic uncertainties on our DOM results associated with missing three-body contributions. Before addressing these critiques, we want to state that the key message of our paper remains whether or not the DOM results agree with the LDM predictions in the nuclear interior. Here, the key point is that the standard determination of the saturation energy from the LDM is not ideal since the total binding energy has a minimal contribution from the core of the nucleus as pointed out in Figs. 1–3 of our paper.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

GEOSH: Ideal Gas Chemical Equation of State

We present the framework and methodology for the new Los Alamos National Laboratory (LANL) G as chemical E quation O f S tate at H igher temperatures code (GEOSH) which aims to accurately model the behavior of chemically complex gaseous mixtures in equilibrium. Assuming the ideal gas approximation, GEOSH leverages the recursive nature of the Saha ionization and molecular equations in order to eliminate the molecular and ionic degrees of freedom, thereby reducing the problem size to the number of atomic species plus one for the free electrons if ions are included. This approach allows the chemical species, both molecular and ionic, of the mixture to be expressed in terms of the abundances of the elemental species. As a result, the GEOSH framework achieves a reduction in computational expense, increased processing speed, and the capability to efficiently model large-scale chemical networks. This report provides the necessary physical background and theoretical foundations for the GEOSH code, accompanied by benchmarking studies.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

X-ray diffraction and equation of state of the C–S–H room-temperature superconductor

X-ray diffraction indicates that the structure of the recently discovered carbonaceous sulfur hydride (C–S–H) room-temperature superconductor is derived from previously established van der Waals compounds found in the H 2 S–H 2 and CH 4 –H 2 systems. Crystals of the superconducting phase were produced by a photochemical synthesis technique, leading to the superconducting critical temperature T c of 288 K at 267 GPa. X-ray diffraction patterns measured from 124 to 178 GPa, within the pressure range of the superconducting phase, are consistent with an orthorhombic structure derived from the Al 2 Cu-type determined for (H 2 S) 2 H 2 and (CH 4 ) 2 H 2 that differs from those predicted and observed for the S–H system at these pressures. The formation and stability of the C–S–H compound can be understood in terms of the close similarity in effective volumes of the H 2 S and CH 4 components, and denser carbon-bearing S–H phases may form at higher pressures. Furthermore, the results are crucial for understanding the very high superconducting T c found in the C–S–H system at megabar pressures.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1982

Integrating matrix solution of the hybrid state vector equations for beam vibration

A simple, versatile, and efficient computational technique has been developed for dynamic analysis of linear elastic beam and rod type of structures. Moreover, the method provides a rather general solution approach for two-point boundary value problems that are described by a single independent spatial variable. For structural problems, the method is implemented by a mixed state vector formulation of the differential equations, combined with an integrating matrix solution procedure. Highly accurate solutions are easily achieved with this approach. Example solutions are given for beam vibration problems including discontinuous stiffness and mass parameters, elastic restraint boundary conditions, concentrated inertia loading, and rigid body modes

Lehman, L. L.↗
NASA NTRS1974

Scaled parametric equation of state for oxygen in the critical region

Thermodynamic properties show an anomalous nonanalytic behavior in the critical region that cannot be represented by analytic equations used in engineering correlations. The formulation of scaling laws is discussed, taking into account an approach proposed by Schofield et al. (1969). A procedure for data analysis is described. The model considered has been found to give satisfactory results in tests conducted with a number of gases. The application of the described approach to density profiles is discussed. Attention is given to the model parameters, the critical exponents, and the coefficients for oxygen.

Levelt Sengers, J. M. H.↗
DOE OSTI2024

Diffusion-Limited Kinetics in Reactive Systems

A proper representation of chemical kinetics is vital to understanding, modeling, and optimizing many important chemical processes. In liquid and surface phases, where diffusion is slow, the rate at which the reactants diffuse together limits the overall rate of many elementary reactions. Commonly, the textbook Smoluchowski theory is utilized to estimate effective rate coefficients in the liquid phase. On surfaces, modelers commonly resort to much more complex and expensive Kinetic Monte Carlo (KMC) simulations. Here, in this study, we extend the Smoluchowski model to allow the diffusing species to undergo chemical reactions and derive analytical formulas for the diffusion-limited rate coefficients for 3D, 2D, and 2D/3D interface cases. With these equations, we are able to demonstrate that when species react faster than they diffuse they can react orders of magnitude faster than predicted by Smoluchowski theory, through what we term “the reactive transport effect”. We validate the derived steady-state equations against particle Monte Carlo (PMC) simulations, KMC simulations, and non-steady-state solutions. Furthermore, using PMC and KMC simulations, we propose corrections that agree with all limits and the computed data for the 2D and 2D/3D interface steady-state equations, accounting for unique limitations in the associated derived equations. Additionally, we derive equations to handle couplings between diffusion-limited rate coefficients in reaction networks. We believe these equations should make it possible to run much more accurate mean-field simulations of liquids, surfaces, and liquid–surface interfaces accounting for diffusion limitations and the reactive transport effect.

Johnson, Matthew Sean↗
NASA NTRS1995

Determination of stability and control derivatives from the NASA F/A-18 HARV from flight data using the maximum likelihood method

This report is a compilation of PID (Proportional Integral Derivative) results for both longitudinal and lateral directional analysis that was completed during Fall 1994. It had earlier established that the maneuvers available for PID containing independent control surface inputs from OBES were not well suited for extracting the cross-coupling static (i.e., C(sub N beta)) or dynamic (i.e., C(sub Npf)) derivatives. This was due to the fact that these maneuvers were designed with the goal of minimizing any lateral directional motion during longitudinal maneuvers and vice-versa. This allows for greater simplification in the aerodynamic model as far as coupling between longitudinal and lateral directions is concerned. As a result, efforts were made to reanalyze this data and extract static and dynamic derivatives for the F/A-18 HARV (High Angle of Attack Research Vehicle) without the inclusion of the cross-coupling terms such that more accurate estimates of classical model terms could be acquired. Four longitudinal flights containing static PID maneuvers were examined. The classical state equations already available in pEst for alphadot, qdot and thetadot were used. Three lateral directional flights of PID static maneuvers were also examined. The classical state equations already available in pEst for betadot, p dot, rdot and phi dot were used. Enclosed with this document are the full set of longitudinal and lateral directional parameter estimate plots showing coefficient estimates along with Cramer-Rao bounds. In addition, a representative time history match for each type of meneuver tested at each angle of attack is also enclosed.

Napolitano, Marcello R.↗
NASA NTRS1973

Study of helium emissions from active solar regions

A theoretical study is made of the visible and UV line radiation of He I atoms and He II ions from a plane-parallel model flare layer. Codes were developed for the solution of the statistically steady state equation for a 30 level He I - II - III model, and the line and continuum transport equations. These codes are described and documented in the report along with sample solutions. Optical depths and some line intensities are presented for a 1000 km thick layer. Solutions of the steady state equations are presented for electron temperatures 10,000 to 50,000 K and electron densities 10 to the 10th power to 10 to the 14th power/cu cm.

Kulander, J. L.↗
NASA NTRS1987

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

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

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↗
NASA NTRS1977

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

Evidence of Free-Bound Transitions in Warm Dense Matter and Their Impact on Equation-of-State Measurements

Warm dense matter (WDM) is now routinely created and probed in laboratories around the world, providing unprecedented insights into conditions achieved in stellar atmospheres, planetary interiors, and inertial confinement fusion experiments. However, the interpretation of these experiments is often filtered through models with systematic errors that are difficult to quantify. Due to the simultaneous presence of quantum degeneracy and thermal excitation, transitions in which free electrons are de-excited into thermally unoccupied bound states transferring momentum and energy to a scattered x-ray photon become viable. Here we show that such free-bound transitions are a particular feature of WDM and vanish in the limits of cold and hot temperatures. The inclusion of these transitions into the analysis of recent X-ray Thomson Scattering experiments on WDM at the National Ignition Facility and the Linac Coherent Light Source is required to obtain a physically consistent temperature from the Chihara decomposition. This interpretation is corroborated by agreement with a recently developed model-free thermometry technique and presents an important step for precisely characterizing and understanding the complex WDM state of matter.

Boehme, M [HZDR Dresden]↗