Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Cubic 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 19 records

Efficient High-Pressure State Equations

A method is presented for a relatively accurate, noniterative, computationally efficient calculation of high-pressure fluid-mixture equations of state, especially targeted to gas turbines and rocket engines. Pressures above I bar and temperatures above 100 K are addressed The method is based on curve fitting an effective reference state relative to departure functions formed using the Peng-Robinson cubic state equation Fit parameters for H2, O2, N2, propane, methane, n-heptane, and methanol are given.

Harstad, Kenneth G.↗

Efficient High Pressure MixtureState Equations

A method is presented for an accurate noniterative, computationally efficient calculation of high pressure fluid mixture equations of state, especially targeted to gas turbines and rocket engines. Pressures above 1 bar and temperatures above 100 K are addressed. The method is based on curve fitting an effective reference state relative to departure funcitons formed using the Peng-Robinson cubic state equation. Fit parameters for H(sub 2), O(sub 2), N(sub 2), propane, n-heptane and methanol are given.

thermodynamics state equations noniterative mixtur↗

Computational Analyses of Pressurization in Cryogenic Tanks

A comprehensive numerical framework utilizing multi-element unstructured CFD and rigorous real fluid property routines has been developed to carry out analyses of propellant tank and delivery systems at NASA SSC. Traditionally CFD modeling of pressurization and mixing in cryogenic tanks has been difficult primarily because the fluids in the tank co-exist in different sub-critical and supercritical states with largely varying properties that have to be accurately accounted for in order to predict the correct mixing and phase change between the ullage and the propellant. For example, during tank pressurization under some circumstances, rapid mixing of relatively warm pressurant gas with cryogenic propellant can lead to rapid densification of the gas and loss of pressure in the tank. This phenomenon can cause serious problems during testing because of the resulting decrease in propellant flow rate. With proper physical models implemented, CFD can model the coupling between the propellant and pressurant including heat transfer and phase change effects and accurately capture the complex physics in the evolving flowfields. This holds the promise of allowing the specification of operational conditions and procedures that could minimize the undesirable mixing and heat transfer inherent in propellant tank operation. In our modeling framework, we incorporated two different approaches to real fluids modeling: (a) the first approach is based on the HBMS model developed by Hirschfelder, Beuler, McGee and Sutton and (b) the second approach is based on a cubic equation of state developed by Soave, Redlich and Kwong (SRK). Both approaches cover fluid properties and property variation spanning sub-critical gas and liquid states as well as the supercritical states. Both models were rigorously tested and properties for common fluids such as oxygen, nitrogen, hydrogen etc were compared against NIST data in both the sub-critical as well as supercritical regimes.

Ahuja, Vineet↗

Electrostriction of Near-Critical SF6 in Microgravity

We used interferometry to measure the electric-field-induced (i.e., "electrostrictive") increase of the density of sulfur hexafluoride near its critical point. The results at the temperatures T(sub c) + 5 mK, T(sub c) + 10 mK, and T(sub c) + 30 mK with T(sub c) = 319 K agree with a calculation based on the Clausius-Mossotti relation and the restricted cubic model equation of state. These measurements were performed in microgravity so that the electrostrictive density changes (less than or equal to 3.5% in this work) would not be complicated by stratification of the fluid's density induced by the Earth's gravity.

Zimmerli, G. A.↗

Electrostriction of a Near-Critical Fluid in Microgravity

We used interferometry to measure the electrostrictive increase of the density of sulfur hexaflouride (SF6) near its critical point. The results at three temperatures (T(sub c) + 5.0 mK, T(sub c) + 10.0 mK, T(sub c) + 30.0 mK with T(sub c) = 319 K) agree with a calculation based on the Clausius-Mossotti relation and the restricted cubic model equation of state. To measure electrostriction, an inhomogeneous electric field (less than or equal to 26 kV/cm) was applied to the SF6 sample by charging a fine wire that passed through it. These measurements were performed in microgravity so that the small electrostrictive density changes (less than or equal to 3.5 % in this work) would not be masked by the larger stratification of the fluid's density induced by the Earth's gravity. The predicted shifts of the critical temperature and density resulting from the electric field were too small to detect.

Zimmerli, Gregory A.↗

A Galerkin method for the estimation of parameters in hybrid systems governing the vibration of flexible beams with tip bodies

An approximation scheme is developed for the identification of hybrid systems describing the transverse vibrations of flexible beams with attached tip bodies. In particular, problems involving the estimation of functional parameters are considered. The identification problem is formulated as a least squares fit to data subject to the coupled system of partial and ordinary differential equations describing the transverse displacement of the beam and the motion of the tip bodies respectively. A cubic spline-based Galerkin method applied to the state equations in weak form and the discretization of the admissible parameter space yield a sequence of approximating finite dimensional identification problems. It is shown that each of the approximating problems admits a solution and that from the resulting sequence of optimal solutions a convergent subsequence can be extracted, the limit of which is a solution to the original identification problem. The approximating identification problems can be solved using standard techniques and readily available software.

Banks, H. T.↗

Approximation methods for inverse problems involving the vibration of beams with tip bodies

Two cubic spline based approximation schemes for the estimation of structural parameters associated with the transverse vibration of flexible beams with tip appendages are outlined. The identification problem is formulated as a least squares fit to data subject to the system dynamics which are given by a hybrid system of coupled ordinary and partial differential equations. The first approximation scheme is based upon an abstract semigroup formulation of the state equation while a weak/variational form is the basis for the second. Cubic spline based subspaces together with a Rayleigh-Ritz-Galerkin approach were used to construct sequences of easily solved finite dimensional approximating identification problems. Convergence results are briefly discussed and a numerical example demonstrating the feasibility of the schemes and exhibiting their relative performance for purposes of comparison is provided.

Rosen, I. G.↗

The effect of an electric field on the morphological stability of the crystal-melt interface of a binary alloy. III - Weakly nonlinear theory

The effect of a constant electric current on the crystal-melt interface morphology during directional solidification at constant velocity of a binary alloy is considered. A linear temperature field is assumed, and thermoelectric effects and Joule heating are neglected; electromigration and differing electrical conductivities of crystal and melt are taken into account. A two-dimensional weakly nonlinear analysis is carried out to third order in the interface amplitude, resulting in a cubic amplitude equation that describes whether the bifurcation from the planar state is supercritical or subcritical. For wavelengths corresponding to the most dangerous mode of linear theory, the demarcation between supercritical and subcritical behavior is calculated as a function of processing conditions and material parameters. The bifurcation behavior is a sensitive function of the magnitude and direction of the electric current and of the electrical conductivity ratio.

Wheeler, A. A.↗

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↗

Spline-based distributed system identification with application to large space antennas

A parameter and state estimation technique for distributed models is demonstrated through the solution of a problem generic to large space antenna system identification. Assuming the position of the reflective surface of the maypole (hoop/column) antenna to be approximated by the static two-dimensional, stretched-membrane partial differential equation with variable-stiffness coefficient functions, a spline-based approximation procedure is described that estimates the shape and stiffness functions from data set observations. For given stiffness functions, the Galerkin projection with linear spline-based functions is applied to project the distributed problem onto a finite-dimensional subspace wherein algebraic equations exist for determining a static shape (state) prediction. The stiffness functions are then parameterized by cubic splines and the parameters estimated by an output error technique. Numerical results are presented for data descriptive of a 100-m-diameter maypole antenna.

Banks, H. T.↗

Weak cubic Langmuir turbulence

The cubically nonlinear Schroedinger equation model of Langmuir turbulence is solved in the weak turbulence limit. Steady-state power-law solutions for the energy spectra are found in arbitrary dimensionality. In one spatial dimension, the theory incorrectly predicts that no spectrum evolves in time. In three spatial dimensions, numerical solutions are obtained for the undriven, undamped, initial value problem and for the driven, damped, initial value problem.

Hansen, P. J.↗

Piezoresistive sensors for measurement of thermally-induced stresses in microelectronics

The general theory of conduction in a piezoresistive material subjected to temperature change is derived and the basic equations needed to design and apply test chips for thermal stress measurement are established. A general expression is obtained for the resistance changes which occur in an arbitrarily oriented 1D filamentary conductor subjected to general 3D stress states and arbitrary temperature changes. The derived formula is valid for crystals such as silicon which have cubic symmetry and diamond lattice structures. Basic equations are presented for the resistance changes experienced by stressed in-plane resistors fabricated on (100) silicon wafers. The established theory is applied to a simple two-element sensor rosette. A test chip containing this rosette configuration is designed and fabricated, and a suitable calibration system is constructed.

Suhling, J. C.↗

Topics in strong Langmuir turbulence

Progress in two approaches to the study of strong Langmuir turbulence is reported. In two spatial dimensions, numerical solution of the Zakharov equations yields a steady state involving linear growth, linear damping, and a collection of coherent, long-lived entities which might loosely be called solitons. In one spatial dimension, a statistical theory is applied to the cubically nonlinear Schroedinger equation and is solved analytically in a special case.

Nicholson, D. R.↗

Topics in strong Langmuir turbulence

Progress in two approaches to the study of strong Langmuir turbulence is reported. In two spatial dimensions, numerical solution of the Zakharov equations yields a steady state involving linear growth, linear damping, and a collection of coherent, long-lived entities which might loosely be called solitons. In one spatial dimension, a statistical theory is applied to the cubically nonlinear Schroedinger equation and is solved analytically in a special case.

Nicholson, D. R.↗

Exact solution of the Lifshitz equations governing the growth of fluctuations in cosmology

The exact solution of the Lifshitz equations governing the cosmological evolution of an initial fluctuation is presented. Lifshitz results valid for squares of the sound velocity equal to zero and 1/3 are extended in closed form to any equation of state where the pressure equals the total energy density times the square of the sound velocity. The solutions embody all the results found previously for special cases of the square of the sound velocity. It is found that the growth of any initial fluctuation is only an exponential function of time with an exponent of not more than 4/3 and is insufficient to produce galaxies unless the initial fluctuation is very large. A possible way to produce very large initial fluctuations by modifying the equation of state by including gravitational interactions is also examined. It is found that a phase transition can occur at baryonic density of 1 nucleon per cubic Planck length or equivalently, at a time of about 10 to the -43rd power sec. At those early times, the masses allowed by causality requirements are too small to be of interest in galaxy formation.

Adams, P. J.↗

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

A Finite Element Theory for Predicting the Attenuation of Extended-Reacting Liners

A non-modal finite element theory for predicting the attenuation of an extended-reacting liner containing a porous facesheet and located in a no-flow duct is presented. The mathematical approach is to solve separate wave equations in the liner and duct airway and to couple these two solutions by invoking kinematic constraints at the facesheet that are consistent with a continuum theory of fluid motion. Given the liner intrinsic properties, a weak Galerkin finite element formulation with cubic polynomial basis functions is used as the basis for generating a discrete system of acoustic equations that are solved to obtain the coupled acoustic field. A state-of-the-art, asymmetric, parallel, sparse equation solver is implemented that allows tens of thousands of grid points to be analyzed. A grid refinement study is presented to show that the predicted attenuation converges. Excellent comparison of the numerically predicted attenuation to that of a mode theory (using a Haynes 25 metal foam liner) is used to validate the computational approach. Simulations are also presented for fifteen porous plate, extended-reacting liners. The construction of some of the porous plate liners suggest that they should behave as resonant liners while the construction of others suggest that they should behave as broadband attenuators. In each case the finite element theory is observed to predict the proper attenuation trend.

Watson, W. R.↗

Finite element computation of a viscous compressible free shear flow governed by the time dependent Navier-Stokes equations

A finite element algorithm for solution of fluid flow problems characterized by the two-dimensional compressible Navier-Stokes equations was developed. The program is intended for viscous compressible high speed flow; hence, primitive variables are utilized. The physical solution was approximated by trial functions which at a fixed time are piecewise cubic on triangular elements. The Galerkin technique was employed to determine the finite-element model equations. A leapfrog time integration is used for marching asymptotically from initial to steady state, with iterated integrals evaluated by numerical quadratures. The nonsymmetric linear systems of equations governing time transition from step-to-step are solved using a rather economical block iterative triangular decomposition scheme. The concept was applied to the numerical computation of a free shear flow. Numerical results of the finite-element method are in excellent agreement with those obtained from a finite difference solution of the same problem.

Cooke, C. H.↗