SEARCH · Engineering Papers

Results for “Thermodynamic functions”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 145 records · Page 8

All sources NASA NTRS DOE OSTI

NASA NTRS1987

Simplified curve fits for the thermodynamic properties of equilibrium air

New, improved curve fits for the thermodynamic properties of equilibrium air have been developed. The curve fits are for pressure, speed of sound, temperature, entropy, enthalpy, density, and internal energy. These curve fits can be readily incorporated into new or existing computational fluid dynamics codes if real gas effects are desired. The curve fits are constructed from Grabau-type transition functions to model the thermodynamic surfaces in a piecewise manner. The accuracies and continuity of these curve fits are substantially improved over those of previous curve fits. These improvements are due to the incorporation of a small number of additional terms in the approximating polynomials and careful choices of the transition functions. The ranges of validity of the new curve fits are temperatures up to 25 000 K and densities from 10 to the -7 to 10 to the 3d power amagats.

Srinivasan, S.↗

DOE OSTI2026

Computational Design of Interlayers for Thermally Stable Compositionally Graded Coatings on Nickel Alloys

To extend the service life of Ni-based superalloys, refractory metal coatings are often used. However, direct bonding between metals with dissimilar crystal structure promotes brittle intermetallic phase formation. This work presents a computational thermodynamic framework for high throughput design of functionally graded interlayers to suppress deleterious phases that may form at the interlayer. The Thermo-Calc software package was used to screen candidate metallic interlayer elements based on stability of solid-solution phases. Vanadium was identified as a promising interlayer due to its consistent suppression of intermetallic phases. Temperature-dependent phase diagram mapping between 600 and 1000 °C guided selection of a compositional pathway that significantly reduced intermetallic formation compared to directly joining the Ni-based and Nb refractory alloys. Time–temperature–transformation analysis was performed to assess whether equilibrium-predicted phases are kinetically accessible along regions of the graded path where non-solid-solution phases are not fully suppressed. The methodology was further applied to additional Ni-based alloy and coating systems, illustrating its transferability as an approach for rapid computational design of graded interlayers in dissimilar high-temperature materials.

36 MATERIALS SCIENCE↗

DOE OSTI2025

First-Principles Studies on Sc 2 RuZ (Z = Si, Ge, Sn) Inverse Heusler Alloys: Structural, Electronic, and Transport Properties

The continuous demand for efficient, nontoxic, and thermally stable materials for room-temperature energy conversion motivates the exploration of novel thermoelectric systems beyond the traditional magnetic Heusler alloys. While full and half-Heusler compounds, especially Co-, Ni-, and Mn-based systems, have demonstrated promising thermoelectric properties, their typically high operating temperatures and magnetic complexities limit their applicability in ambient thermal management. In this context, we investigate whether Sc-based inverse Heusler alloys can offer a viable nonmagnetic alternative with competitive thermoelectric performance. In this work, we perform a systematic first-principles study of the inverse Heusler compounds Sc 2 RuZ (Z = Si, Ge, Sn), focusing on their structural, electronic, mechanical, and thermodynamic-thermoelectric properties. Density Functional Theory (DFT) was employed to compute optimized lattice structures and band dispersion, while dynamical stability was assessed via phonon calculations. Thermoelectric transport coefficients, including Seebeck coefficient, electrical conductivity, and thermal conductivity, were estimated using the semiclassical Boltzmann transport theory within the constant relaxation time approximation. Our results show that all Sc 2 RuZ compounds are thermodynamically stable semiconductors with indirect band gaps of 0.12–0.16 eV and exhibit high elastic moduli, especially Sc 2 RuSn, which demonstrates superior stiffness and incompressibility. Importantly, all compounds display promising room-temperature thermoelectric characteristics, including high Seebeck coefficients and power factors. These findings reveal that Sc 2 RuZ alloys represent a rare class of stable, nonmagnetic inverse Heusler semiconductors with intrinsic thermoelectric potential at room temperature, unlike many existing Heusler systems optimized for spintronics or high-temperature operation. This work expands the known design space for Heusler-based thermoelectrics and offers a theoretical basis for experimental realization of efficient, low-temperature, nonmagnetic thermoelectric materials.

alloys↗

NASA NTRS1986

Simplified curve fits for the thermodynamic properties of equilibrium air

New improved curve fits for the thermodynamic properties of equilibrium air were developed. The curve fits are for p = p(e,rho), a = a(e,rho), T = T(e,rho), s = s(e,rho), T = T(p,rho), h = h(p,rho), rho = rho(p,s), e = e(p,s) and a = a(p,s). These curve fits can be readily incorporated into new or existing Computational Fluid Dynamics (CFD) codes if real-gas effects are desired. The curve fits were constructed using Grabau-type transition functions to model the thermodynamic surfaces in a piecewise manner. The accuracies and continuity of these curve fits are substantially improved over those of previous curve fits appearing in NASA CR-2470. These improvements were due to the incorporation of a small number of additional terms in the approximating polynomials and careful choices of the transition functions. The ranges of validity of the new curve fits are temperatures up to 25,000 K and densities from 10 to the minus 7th to 100 amagats (rho/rho sub 0).

Srinivasan, S.↗

DOE OSTI2024

Enhancing the accuracy and generality of the Debye–Grüneisen Model: Optimizing the volume dependence for accurate predictions across varied compositions

In this work, we have introduced an optimized Debye-Grüneisen model that revolutionizes the determination of the Debye temperature and Grüneisen parameters. Unlike conventional methods, our model requires only the 0 K energy volume data for a material as input, eliminating the need to determine the bulk modulus and its pressure derivative, which often pose challenges due to numerical uncertainties. This unique feature sets our model apart from existing approaches and streamlines the process, enabling accurate predictions of thermal expansion behavior across various materials. To demonstrate its effectiveness, we showcase its excellent agreement with measured coefficients of thermal expansion (CTE) for the nickel-cobalt-chromium-aluminum-yttrium (Ni-Co-Cr-Al-Y) bond-coating system. Additionally, we apply our approach by conducting a high-throughput search for potential bond-coating materials among 90,000 compositions within the aluminum-cobalt-chromium-iron-nickel (Al-Co-Cr-Fe-Ni) system. From this extensive search, four compositions are synthesized, and the measured CTE values agree very well with theoretical predictions, hence validating our approach. In conclusion, the current optimized Debye-Grüneisen model combined with Density Functional Theory (DFT)-based thermodynamic database enables reliable and efficient high-throughput calculations of CTE of of a material without expensive phonon calculations.

Bond coating materials↗

DOE OSTI2025

Correlation effects on coupled electronic and structural properties of doped rare-earth trihydrides

Rare-earth trihydride (𝑅⁢H 3 ) compounds exhibit intriguing coupled electronic and structural properties as a function of doping, hydrogen vacancies, and thermodynamic conditions. Theoretical studies of these materials typically rely on density functional theory (DFT), including the use of small supercells that may underestimate strong correlation effects and structural distortions which in turn may influence their metallicity. Here, we elucidate the roles of lattice distortions and correlation effects on the electronic properties of pristine and doped 𝑅⁢H 3 compounds by adopting DFT +U and quantum Monte Carlo (QMC) methods. Linear-response constrained DFT (LR-cDFT) methods find Hubbard U ≈ 2 eV for 𝑅 𝑑 orbitals and U ≈ 6 eV for H 𝑠 ⁡/N 𝑝 /O 𝑝 orbitals. The small U on Lu 𝑑 orbitals is consistent with QMC calculations on LuH 3 and LuH 2.875 ⁢N 0.125 . In pure face-centered-cubic (FCC) 𝑅⁢H 3 (𝑅 = Lu,Y) compounds, neither DFT nor DFT +U with the self-consistently determined U is enough to create a band gap, however a supercell with hydrogen distortions creates a small gap whose magnitude increases when performing DFT +U with self-consistently determined U values. Furthermore, correlation effects, in turn, have a moderate influence on the coupled structural and electronic properties of doped 𝑅⁢H 3 compounds and may be important when considering the competition between structural distortions and superconductivity.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

NASA NTRS1973

Thermodynamic Properties of Nitrogen Including Liquid and Vapor Phases from 63K to 2000K with Pressures to 10,000 Bar

Tables of thermodynamic properties of nitrogen are presented for the liquid and vapor phases for temperatures from the freezing line to 2000K and pressures to 10,000 bar. The tables include values of density, internal energy, enthalpy, entropy, isochoric heat capacity, isobaric heat capacity velocity of sound, the isotherm derivative, and the isochor derivative. The thermodynamic property tables are based on an equation of state, P=P (p,T), which accurately represents liquid and gaseous nitrogen for the range of pressures and temperatures covered by the tables. Comparisons of property values calculated from the equation of state with measured values for P-p-T, heat capacity, enthalpy, latent heat, and velocity of sound are included to illustrate the agreement between the experimental data and the tables of properties presented here. The coefficients of the equation of state were determined by a weighted least squares fit to selected P-p-T data and, simultaneously, to isochoric heat capacity data determined by corresponding states analysis from oxygen data, and to data which define the phase equilibrium criteria for the saturated liquid and the saturated vapor. The vapor pressure equation, melting curve equation, and an equation to represent the ideal gas heat capacity are also presented. Estimates of the accuracy of the equation of state, the vapor pressure equation, and the ideal gas heat capacity equation are given. The equation of state, derivatives of the equation, and the integral functions for calculating derived thermodynamic properties are included.

Thermodynamic properties↗

DOE OSTI2024

Atomic structure of different surface terminations of polycrystalline ZnPd

The intermetallic compound ZnPd has been found to have desirable characteristics as a catalyst for the steam reforming of methanol. The understanding of the surface structure of ZnPd is important to optimize its catalytic behavior. However, due to the lack of bulk single-crystal samples and the complexity of characterizing surface properties in the available polycrystalline samples using common experimental techniques, all previous surface science studies of this compound have been performed on surface alloy samples formed through thin-film deposition. In this study, we present findings on the chemical and atomic structure of the surfaces of bulk polycrystalline ZnPd studied by a variety of complementary experimental techniques, including scanning tunneling microscopy (STM), x-ray photoelectron spectroscopy (XPS), low energy electron microscopy (LEEM), photoemission electron microscopy (PEEM), and microspot low-energy electron diffraction ( μ -LEED). These experimental techniques, combined with density functional theory (DFT)-based thermodynamic calculations of surface free energy and detachment kinetics at the step edges, confirm that surfaces terminated by atomic layers composed of both Zn and Pd atoms are more stable than those terminated by only Zn or Pd layers. DFT calculations also demonstrate that the primary contribution to the tunneling current arises from Pd atoms, in agreement with the STM results. The formation of intermetallics at surfaces may contribute to the superior catalyst properties of ZnPd over Zn or Pd elemental counterparts. Published by the American Physical Society 2024

36 MATERIALS SCIENCE↗

DOE OSTI2026

Thermal exchange-correlation functionals: Capturing quantum electron behavior in warm, dense plasmas

We summarize and give perspective upon recent progress in developing non-empirical constraint-based thermal (i.e., free energy) exchange-correlation (XC) density functionals essential for accurate description of the quantum behavior of electrons in warm, dense plasmas. After delineating the critical role of ground-state functionals for zero-temperature, time-dependent DFT, we outline the underpinnings of local density approximation, generalized gradient approximation (GGA), and meta-GGA XC free-energy functionals. Two basic thermalization principles for upgrading ground-state XC functionals to successful thermal ones are emphasized. Then, we turn to a long-standing challenge, assessment of the accuracy of well-founded functionals. Unlike the ground state, there are a few exact results for large T and P. An exception is path integral Monte Carlo (PIMC) data for dense H/D and He plasmas. For those, we did ab initio molecular dynamics simulations under selected thermodynamic conditions employing five thermal XC functionals: two approximate thermal GGAs, fully thermal GGA, an approximate meta-GGA, and fully thermal meta-GGA. Comparisons with the PIMC data show that functionals thermalized by augmenting a non-thermal functional with a lower-level thermal contribution are inferior to functionals with thermal XC and spatial inhomogeneity effects taken into account at the same level of refinement. We believe this and similar evidence should be convincing to the high-energy density physics community of the necessity of use of proper thermal XC functionals in simulation studies of finite-temperature quantum effects in warm, dense plasmas.

Ab-initio molecular dynamics↗

DOE OSTI2025

Understanding CO 2 adsorption on the surfaces of SrO and its hydroxylated variants Sr(OH) 2 · n H 2 O ( n = 0, 1, 8)

Strontium oxide (SrO) is a promising material for CO 2 capture through a reversible cycle of carbonation and calcination, where SrO reacts with CO 2 to form SrCO 3 and can be regenerated by calcination. In the presence of moisture, SrO forms strontium hydroxide and its hydrates (Sr(OH) 2 ·nH 2 O). Here, this study, which employs density functional theory, investigates the CO 2 adsorption mechanism of these processes on various crystal surfaces, including SrO, Sr(OH) 2 , Sr(OH) 2 ·1H 2 O, and Sr(OH) 2 ·8H 2 O. A significant finding is that the interaction of CO 2 with these surfaces leads to carbonate/bicarbonate formation via electron transfer, with notable differences in CO 2 orientation and bond characteristics between SrO surfaces and its hydroxylated surfaces. To explore the effects of moisture on CO 2 adsorption, H 2 O adsorption on these stable surfaces was investigated. The results showed that H 2 O reacts with the SrO (100) surface to form hydroxyl (OH) groups while it bonds with the surfaces of Sr(OH) 2 ·nH 2 O (n = 0, 1, 8) with hydrogen bonding. A small amount of H 2 O can enhance CO2 adsorption while a large amount of H 2 O could decrease the capability of CO 2 adsorption.

CO2 capture↗

DOE OSTI2026

MoNbTi-based Refractory Multi-principal Element Alloy System: A Review of the Thermodynamic Phase Predictions, Formed Microstructures, and Mechanical Properties as a Function of the Fabrication Methods (AM versus SPS versus VAM)

Abstract Multi-principal element alloys, particularly of refractory compositions, are an increasingly popular candidate for extreme-environment applications, including for next-generation nuclear reactors and in other industries, such as biomedical and aerospace, due to their high strength. The ability to achieve solid-solution microstructures provides these alloys with improved mechanical properties; however, these microstructures are heavily composition- and processing-dependent. In this review, the multi-principal element alloy MoNbTi-based system (with Zr, V, and Cr additions) was selected to compare the effects of elemental composition and processing route on the resulting microstructures and mechanical properties. The review provides insight into identifying the optimal alloy composition and processing route for a balance of desired microstructure and mechanical properties.

Krogh, Kara↗

DOE OSTI2025

Fundamental Insights into Cathode Stability: Linking Compositional Tuning and Local Coordination in Complex Metal Oxides under Aqueous Transformations

Compositional tuning of complex metal oxides in Li-ion battery materials influences their performance as well as their end-of-life behavior, in particular, the tendency to release toxic metal cations in aqueous solution. We modeled ternary variants of a parent LiCoO 2 delafossite structure by varying the metal identity and relative amounts. This yielded ten model formulations of Li(A 4/6 B 1/6 C 1/6 )O 2 , where the material is enriched with the A metal and doped with B and C, with Ni, Mn, Co, Fe, Al, V, and Ti as constituent metals. To assess their stability in aqueous conditions, metal release energetics were calculated using a combination of Density Functional Theory calculations and thermodynamics. Metal release in ternary oxides is dictated by subtle variations in the coordination environment of the leaving group. To identify governing chemical features across diverse compositions with varying local coordination environments, we leverage random forest regression and descriptor importance analysis. A key result is that metal–oxygen orbital hybridization, quantified using a projected density-of-states-derived descriptor, H d/p , provides a physically grounded measure of interaction strength that governs metal release energetics. This refined perspective goes beyond conventional oxidation state considerations and offers more robust insights for materials science. Finally, we model defect surface-bound O 2 dimer formation as a proxy for reactive oxygen species (ROS) generation. The results show that Ni-rich compositions more readily stabilize spin-polarized O 2 dimers, corroborating experimental reports of an increased ROS-driven biological response. In conclusion, our results establish a compositional and electronic basis for metal release and surface oxygen reactivity that form a rationale for complex metal oxide design principles.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Dissociation of gas-phase anisole induced by low-energy electron interactions: understanding patterns of aromatic bond cleavage

Abstract In view of elucidating the fragmentation patterns of aromatic systems induced by low-energy electron interactions, dissociative electron attachment (DEA) to gas-phase anisole was performed. Anionic fragments resulting from this DEA process were detected by a quadrupole mass spectrometer, and ion yields of those fragments as a function of incident electron energy were rendered. Our study showed the formation of CH 3 − , HCC − , and OCH 3 − fragments, suggesting that various dissociation channels proceed out of DEA to anisole. We employed density functional theory to compute thermodynamic threshold energies for each potential dissociation channel. Those theoretical calculations supported the prediction that the CH 3 − and OCH 3 − fragments form via mechanisms of single-bond cleavage; the HCC − fragments may form through two-, three-, or four-body dissociation channels that entail hydrogen transfers and the cleavage of multiple aromatic bonds. The experimental resonance energies that form the CH 3 − , HCC − , and OCH 3 − fragments were 6.0 eV, 5.8 and 9.7 eV, and 9.8 eV, respectively. Given the classification of anisole as a monosubstituted aromatic species, our results explain generalizable patterns of electron-mediated dissociation in aromatic systems.

Finley, Jacob (ORCID:0009000606963352)↗

NASA NTRS1985

Diurnal movements of cotton leaves expressed as thermodynamic work and entropy changes

It is pointed out that some important agricultural crops show heliotropic leaf movements. In these species, the proclivity of leaves to orient either perpendicularly or parallel or in some combination of these positions with respect to the sun is controlled by the leaf turgor and the availability of water. Such an orientational response is particularly noticeable for cotton. Schutt et al. (1985) have detailed leaf trajectories using three angles. The present investigation applies the three-angle representation to leaf trajectory mapping and to the calculation of the phase angle 'gamma' between the individual leaf normals and the solar direction. Using gamma, the thermodynamic work and entropy functions are evaluated and used to distinguish between the behavior of water-stressed and well watered cotton canopies.

Schutt, J. B.↗

NASA NTRS2003

Probabilistic Analysis of Solid Oxide Fuel Cell Based Hybrid Gas Turbine System

The emergence of fuel cell systems and hybrid fuel cell systems requires the evolution of analysis strategies for evaluating thermodynamic performance. A gas turbine thermodynamic cycle integrated with a fuel cell was computationally simulated and probabilistically evaluated in view of the several uncertainties in the thermodynamic performance parameters. Cumulative distribution functions and sensitivity factors were computed for the overall thermal efficiency and net specific power output due to the uncertainties in the thermodynamic random variables. These results can be used to quickly identify the most critical design variables in order to optimize the design and make it cost effective. The analysis leads to the selection of criteria for gas turbine performance.

Gorla, Rama S. R.↗

NASA NTRS1983

A Potential Function Derivation of a Constitutive Equation for Inelastic Material Response

Physical and thermodynamic concepts are used to develop a potential function for application to high temperature polycrystalline material response. Inherent in the formulation is a differential relationship between the potential function and constitutive equation in terms of the state variables. Integration of the differential relationship produces a state variable evolution equation that requires specification of the initial value of the state variable and its time derivative. It is shown that the initial loading rate, which is directly related to the initial hardening rate, can significantly influence subsequent material response. This effect is consistent with observed material behavior on the macroscopic and microscopic levels, and may explain the wide scatter in response often found in creep testing.

Stouffer, D. C.↗

NASA NTRS1992

Evolution of a terrestrial magma ocean: Thermodynamics, kinetics, rheology, convection, differentiation

The evolution of an initially totally molten magma ocean is constrained on the basis of analysis of various physical problems in the magma ocean. First of all an equilibrium thermodynamics of the magma ocean is developed in the melting temperature range. The equilibrium thermodynamical parameters are found as functions only of temperature and pressure and are used in the subsequent models of kinetics and convection. Kinematic processes determine the crystal size and also determine a non-equilibrium thermodynamics of the system. Rheology controls all dynamical regimes of the magma ocean. The thermal convection models for different rheological laws are developed for both the laminar convection and for turbulent convection in the case of equilibrium thermodynamics of the multiphase system. The evolution is estimated on the basis of all the above analysis.

Solomatov, V. S.↗

DOE OSTI2025

Insights into the Oxidation Mechanism of Vivianite to Metavivianite from First-Principles Calculations

Vivianite, a hydrous ferrous iron-phosphate mineral (Fe 3 (PO 4 ) 2 ·8H 2 O), readily oxidizes in contact with air yielding the less hydrous mixed-valent iron-phosphate mineral metavivianite. This topotactic transformation nominally occurs by oxidative dehydrogenation in which outgoing electrons from the iron sublattice are charge compensated by hydrolysis of structural water. However, the details of this internal charge balancing mechanism that allows the structure to remain electrostatically stable remain unknown. Here, in this study, we use density functional theory (DFT) calculations and ab initio thermodynamics (AIT) to evaluate the energetics of this process in terms of hydrogen release as a function of environmental variables such as partial pressure and temperature. The results show a thermodynamic driving force for vivianite phase transformation as oxidation progresses that is triggered by its rigid structure that has a limited accommodation for hydrogen vacancies. In contrast, metavivianite has a more flexible lattice and hydrogen bond network that stabilizes these hydrogen defects. Metavivianite is shown to be a stable intermediate for 66% residual Fe 2+ down to 33% Fe 2+ , below which other phases/structures such as santabarbarite should be more thermodynamically favorable. Our study provides a basis for experimental tests of our mechanistic findings and helps fill a basic knowledge gap about the solid-state process that defines how vivianite interacts with its surrounding environment.

Sassi, Michel [Pacific Northwest National Laborato↗