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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Direct numerical simulation of supercritical oxy-methane mixing layers with CO 2 substituted counterparts

Direct numerical simulations (DNS) of temporally developing, three dimensional, CH 4 /CO 2 , CH 4 /O 2 and CO 2 /O 2 mixing layers, are conducted at a supercritical pressure of 300 atm. To effectively model the supercritical regime, the employed formulation includes the compressible form of the governing equations, the cubic Peng-Robinson equation of state and a generalized formulation for heat and mass flux vectors derived from non-equilibrium thermodynamics and fluctuation theory. A linear inviscid stability analysis is also performed for each case, to determine its most unstable wavelength. Flow visualizations reveal the presence of high density gradient magnitude regions for all three mixing layers, with conditional averages indicating increased presence of heavier fluid species within these regions. No significant departures are observed from perfect gas behavior, with compressibility factors very close to unity for all three mixing cases. Applicability of presumed probability density function methods (PDF) is examined for the three supercritical mixing layers. An a priori analysis is also conducted to investigate various simplifying assumptions employed in modeling various subgrid scale (SGS) flux models. Two additional terms are identified in the large eddy simulations (LES) equations, the gradient of SGS contribution of pressure in the momentum equation and the gradient of SGS contribution of heat flux in energy equation, whose magnitudes are similar and comparable with their respective resolved terms. The performance of the scale similarity model to represent these additional terms is investigated. Lastly, the performance of Smagorinsky, gradient and scale similarity models are also investigated.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Dominant effect of anharmonicity on the equation of state and thermal conductivity of MgO under extreme conditions

Understanding heat transfer in the earth mantle region is of great scientific interest but highly challenging due to dissimilar conditions from our daily lives. Thermal conductivity of highly pressured and thermally activated materials can be completely different from that at ambient condition. In this study we calculate the equation of state and thermal conductivity of cubic crystalline MgO, a major component of the mantle, in such an environment. It is shown that the material properties are not accurately captured unless rigorous treatment of anharmonicity of phonon vibrational motion is explicitly applied. To predict renormalized phonon dispersion relations, interatomic force constants of high-order expansion are required, especially at elevated temperature. The anharmonicity plays an increasingly substantial role for regulating volume expansion as temperature rises, because the Gibbs free energy of MgO considerably deviates from what was predicted by the quasiharmonic approximation. We find that the lattice thermal conductivity of MgO and the heat transfer of the earth's mantle are reliably estimated only when both the temperature-dependent renormalization of phonon frequencies and four-phonon scattering effects are taken into consideration. We map the characteristics at Geotherm into contour diagrams in which equilibrium volumes and the thermal conductivities are described as functions of temperature and pressure. Furthermore, the thermal conductivity values of MgO in this study are smaller than those in previous theoretical studies due to the additional incorporation of four-phonon scatterings in the calculations. Specifically, deviation between the values in previous studies and in this study becomes larger as thermodynamic variables reach those near the core-mantle boundary. For instance, thermal conductivities of 52.15 and 23.67 Wm -1 K -1 are estimated at ambient and extreme conditions (T = 4000K and P = 135 GPa), respectively.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Density-Functional Tight-Binding Parameters for Bulk Zirconium: A Case Study for Repulsive Potentials

Density-functional tight-binding (DFTB) parameters are presented for the simulation of the bulk phases of zirconium. Electronic parameters were obtained using a band structure fitting strategy, while two-center repulsive potentials were created by particle swarm optimization. As objective functions for the repulsive potential fitting, we employed the Birch–Murnaghan equations of state for hexagonal close-packed (HCP), body-centered cubic (BCC) and ω phases of Zr from density-functional theory (DFT). When fractional atomic coordinates are not allowed to change in the generation of the equation-of-state curves, long-range repulsive DFTB potentials are able to almost perfectly reproduce equilibrium structures, relative DFT energies of the bulk phases, and bulk moduli. However, the same potentials lead to artifacts in the DFTB potential energy surfaces when atom positions in the unit cell are allowed to fully relax during the change of unit cell parameters. Finally, conventional short-range repulsive DFTB potentials, while inferior in their ability to reproduce DFT bulk energetics, are able to correctly reproduce the qualitative shape of the DFT potential energy surfaces, including the location of global minima, and can therefore be considered more transferable.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1985

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

Melting line and thermal equation of state of fcc-cobalt: A combined experimental and computational approach

The melting line of cobalt has been investigated both experimentally, using synchrotron X-ray diffraction coupled with laser-heated diamond anvil cells, and theoretically, using ab initio simulations. Over the investigated pressure and temperature range – between 30 and 100 GPa and from ambient temperature up to 4000 K – the hexagonal close-packed structure, stable at ambient conditions, is replaced at high temperature by the face-centered cubic structure, observed stable till melting. The melting temperatures obtained by the two methods are in remarkable agreement and the melting line can be well described by a Simon–Glatzel equation of the form T$_m$ = 1768(K)(P(GPa)/35.62+1) 0.64 . Finally, from the obtained results it was possible to determine a thermal equation of state for the cubic face-centered phase of Co.

36 MATERIALS SCIENCE↗
NASA NTRS1984

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

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

Refined room-temperature equation of state of Bi up to 260 GPa

At room temperature, bismuth undergoes several structural transitions with increasing pressure before taking on a body-centered cubic (bcc) phase at approximately 8 GPa. The bcc structure is stable to the highest measured pressure and its simplicity, along with its high compressibility and atomic number, makes it an enticing choice as a pressure calibrant. Here, we present three data sets on the compression of bismuth in a diamond anvil cell in a neon pressure medium, up to a maximum pressure of about 260 GPa. The use of a soft pressure medium reduces deviatoric stress when compared to previous work. With an expanded pressure range, a higher point density, and a decreased uniaxial stress component, we are able to provide more reliable equation of state parameters. We also conduct density functional theory electronic-structure calculations that confirm that the bcc phase is energetically favored at high pressure.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS2020

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

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

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

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

Ab Initio Phase Diagram of Chromium to 2.5 TPa

Chromium possesses remarkable physical properties such as hardness and corrosion resistance. Chromium is also a very important geophysical material as it is assumed that lighter Cr isotopes were dissolved in the Earth’s molten core during the planet’s formation, which makes Cr one of the main constituents of the Earth’s core. Unfortunately, Cr has remained one of the least studied 3d transition metals. In a very recent combined experimental and theoretical study (Anzellini et al., Scientific Reports, 2022), the equation of state and melting curve of chromium were studied to 150 GPa, and it was determined that the ambient body-centered cubic (bcc) phase of crystalline Cr remains stable in the whole pressure range considered. However, the importance of the knowledge of the physical properties of Cr, specifically its phase diagram, necessitates further study of Cr to higher pressure. In this work, using a suite of ab initio quantum molecular dynamics (QMD) simulations based on the Z methodology which combines both direct Z method for the simulation of melting curves and inverse Z method for the calculation of solid–solid phase transition boundaries, we obtain the theoretical phase diagram of Cr to 2.5 TPa. We calculate the melting curves of the two solid phases that are present on its phase diagram, namely, the lower-pressure bcc and the higher-pressure hexagonal close-packed (hcp) ones, and obtain the equation for the bcc-hcp solid–solid phase transition boundary. We also obtain the thermal equations of state of both bcc-Cr and hcp-Cr, which are in excellent agreement with both experimental data and QMD simulations. We argue that 2180 K as the value of the ambient melting point of Cr which is offered by several public web resources (“Wikipedia,” “WebElements,” “It’s Elemental,” etc.) is most likely incorrect and should be replaced with 2135 K, found in most experimental studies as well as in the present theoretical work.

equation of state↗
DOE OSTI2021

Ab initio anharmonic thermodynamic properties of cubic Ca Si O 3 perovskite

We present an ab initio study of the thermodynamic properties of cubic CaSiO 3 perovskite (CaPv) over the pressure and temperature range of the Earth's lower mantle. We compute the anharmonic phonon dispersions throughout the Brillouin zone by utilizing the phonon quasiparticle approach, which characterizes the intrinsic temperature dependence of phonon frequencies and, in principle, captures full anharmonicity. Such temperature-dependent phonon dispersions are used to calculate ab initio free energy in the thermodynamic limit (N → ∞) within the framework of the phonon gas model. Accurate free energy calculations enable us to investigate cubic CaPv's thermodynamic properties, e.g., thermal expansivity, Grüneisen parameter, bulk modulus, heat capacity and thermal equation of state, where anharmonic effects are demonstrated. Furthermore, the present methodology provides an important theoretical approach for exploring phase boundaries, thermodynamic, and thermoelastic properties of strongly anharmonic materials at high pressures and temperatures.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Characterization of the high-pressure and high-temperature phase diagram and equation of state of chromium

The high-pressure and high-temperature phase diagram of chromium has been investigated both experimentally (in situ), using a laser-heated diamond-anvil cell technique coupled with synchrotron powder X-ray diffraction, and theoretically, using ab initio density-functional theory simulations. In the pressure–temperature range covered experimentally (up to 90 GPa and 4500 K, respectively) only the solid body-centred-cubic and liquid phases of chromium have been observed. Experiments and computer calculations give melting curves in agreement with each other that can both be described by the Simon–Glatzel equation $T_{m}(P) = 2136K (1 + P/25.9)^{0.41}$ In addition, a quasi-hydrostatic equation of state at ambient temperature has been experimentally characterized up to 131 GPa and compared with the present simulations. Both methods give very similar third-order Birch–Murnaghan equations of state with bulk moduli of 182–185 GPa and respective pressure derivatives of 4.74–5.15. According to the present calculations, the obtained melting curve and equation of state are valid up to at least 815 GPa, at which pressure the melting temperature is 9310 K. Finally, from the obtained results, it was possible to determine a thermal equation of state of chromium valid up to 65 GPa and 2100 K.

36 MATERIALS SCIENCE↗
NASA NTRS1982

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