Search NASA⌕ Search

SEARCH · Search NASA

Results for “Entropy stability”

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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
NASA NTRS1979

Comparison between thermochemical and phase stability data for the quartz-coesite-stishovite transformations

Phase stability and elasticity data have been used to calculate the Gibbs free energy, enthalpy, and entropy changes at 298 K and 1 bar associated with the quartz-coesite and coesite-stishovite transformations in the system SiO2. For the quartz-coesite transformation, these changes disagree by a factor of two or three with those obtained by calorimetric techniques. The phase boundary for this transformation appears to be well determined by experiment; the discrepancy, therefore, suggests that the calorimetric data for coesite are in error. Although the calorimetric and phase stability data for the coesite-stishovite transformation yield the same transition pressure at 298 K, the phase-boundary slopes disagree by a factor of two. At present, it is not possible to determine which of the data are in error. Thus serious inconsistencies exist in the thermodynamic data for the polymorphic transformations of silica.

Weaver, J. S.↗
DOE OSTI2023

Role of Anions in Stabilizing the [Zn–Al] Layered Double Hydroxides: A Thermodynamic Study

Here, room-temperature acid solution calorimetry, high-temperature oxide melt solution calorimetry, and low-temperature heat capacity measurements were employed to calculate the thermodynamic stabilities of the [Zn–Al–X] layered double hydroxides (LDH) containing different anions (X = Cl – , CO 3 2– , and SO 4 2– ). Cryogenic heat capacity measurements demonstrated a Schottky-type anomaly in the heat capacity of all three LDHs below 11 K. This anomaly is attributed to the tunneling of protons between adjacent oxygen atoms in the LDH interlayer as this creates an energy system similar to a two-level system modeled with a Schottky term. These heat capacity measurements were also used to determine vibrational entropies which, when combined with configurational entropies, provide standard entropies of these LDHs. Enthalpies of formation of LDHs from binary components were determined and combined with the entropies of formation to calculate Gibbs free energies. Based on these values, the order of stability is [Zn–Al–SO 4 ] > [Zn–Al–CO 3 ] > [Zn–Al–Cl]. This trend results from a combination of the interlayer spacing, amount of water in the interlayer, interactions among the interlayer species, and interactions between the metal hydroxide layer and the interlayer.

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

Materials Design Directions for Solar Thermochemical Water Splitting

The sustainable, economical production of molecular hydrogen is a crucial component of a net zero-greenhouse-gas-emissions future. Solar thermochemical water splitting (STWS) offers a renewable route to hydrogen with the potential to help decarbonize several industries, including transportation, manufacturing, mining, metals processing, and electricity generation, as well as provide sustainable hydrogen as a chemical feedstock. STWS uses high temperatures generated from concentrated sunlight or other sustainable means for high-temperature heat to produce hydrogen and oxygen from steam. For example, in its simplest form of a two-step thermochemical cycle, a redox-active metal oxide is heated to ≈1700-2000 K, driving off molecular oxygen while producing oxygen vacancies in the material. The reduced metal oxide then cools (ideally with the extracted heat recuperated for re-use) and, in a separate step, comes into contact with steam, which reacts with oxygen vacancies to produce molecular hydrogen while recovering the original state of the metal oxide. Despite its promising use of the entire solar spectrum to split water thermochemically, the current estimated cost of hydrogen produced via STWS is ≈4-6× the U.S. Department of Energy (DOE) Hydrogen Shot target value of $1/kg. One contributing approach to bridging this cost gap is the design of new materials with improved thermodynamic properties to enable higher efficiencies. The state-of-the-art (SOA) redox-active metal oxide for STWS is ceria (CeO 2 ), due to its close to optimal, although too high, oxygen vacancy formation enthalpy and large configurational and electronic entropy of reduction. However, ceria requires high operating temperatures and its efficiency is insufficient. Therefore, efforts to increase the efficiency of STWS cycles have focused on further optimizing oxygen vacancy formation enthalpies and augmenting the reduction entropy via substitution or doping and materials discovery schemes. Examples of the latter include the perovskites BaCe 0.25 Mn 0.75 O 3 and (Ca,Ce)(Ti,Mn)O 3 . These efforts and others have revealed intuitive chemical principles for the efficient and systematic design of more effective materials, such as the strong correlation between the enthalpies of crystal bond dissociation and solid-state cation reduction with the enthalpy of oxygen vacancy formation, as well as configurational entropy augmentation via the coexistence of two or more redox-active cation sublattices. The purpose of this chapter is to prepare the reader with an up-to-date account of STWS redox-active materials, both the SOA and promising newcomers, as well as to provide chemically intuitive strategies for improving their cycle efficiencies through materials design – in conjunction with ongoing efforts in reactor engineering and gas separations – to reach the cost points for commercial viability. First, we will introduce the thermodynamics of STWS using a two-step, metal-oxide, thermochemical cycle with economics in mind. We also will compare the pros and cons of processes that do or do not involve phase changes. Second, we will describe the qualities that make ceria the SOA STWS redox-active material, as well as its limitations. Third, we will survey some of the most promising candidates to date in the search for materials to supplant ceria, emphasizing the post-ternary, metal-oxide-perovskite alloys. Lastly, we will enumerate and discuss the following materials design directions for STWS redox-active materials: crystal reduction potentials as a proxy for oxygen vacancy formation enthalpies, engineering the electronic and configurational entropy of reduction via f-shells and simultaneous redox, and vetting materials stability via temperature-dependent phase diagrams and melting-point prediction.

08 HYDROGEN↗
NASA NTRS1992

A quasi-static magnetospheric convection model in two-dimensions

A self-consistent 2D model of subsonic plasma sheet convection is presented. Specifically, time sequences of static equilibrium solutions for the 2D magnetospheric magnetic field are constructed consistent with adiabatic convection. This model self-consistently includes a dipole field and a reasonable accounting for the effects of inner magnetospheric shielding. Starting from a relaxed magnetospheric equilibrium, the earthward convection of plasma sheet flux tubes results in the stretching of inner plasma sheet field lines, the development of a local minimum in the equatorial magnetic field Be in the near-earth plasma sheet, and an increasing lobe magnetic field. This evolution in time occurs independent of the specific magnetopause or far-tail boundary conditions, provided the plasma sheet is not at marginal interchange stability. This behavior results solely from the convection of flux tubes of increasing entropy into the near-earth plasma sheet. These results are discussed in the general context of earthward convection in earth's plasma sheet and in the specific context of magnetospheric substorms.

Erickson, G. M.↗
DOE OSTI2025

Dehydration of Methyl Lactate on Alkali Cation-Exchanged Faujasite: Effects of Metal Cation Identity and Water Pressure

Turnover rates for the catalytic dehydration of methyl lactate (ML) over ion-exchanged faujasite (FAU) catalysts depend on the identity of alkali metal cations (Na + , K + , Cs + ) and local solvation effects. Analysis of rate measurements and in situ infrared spectroscopy gives evidence that the reaction involves kinetically relevant dissociation of adsorbed methyl lactate upon alkali metal cations. This process involves concerted methyl transfer to the surface and dissociation of the alkali metal from the framework, which occurs at cationic active sites that remain predominantly unoccupied under relevant conditions (0.5–10 kPa ML, 0.5–15 kPa H 2 O, 563–583 K). Despite the mechanistic similarities, apparent activation enthalpies (ΔH app ‡ ) decrease linearly (47 kJ mol –1 from Na + to Cs + ) with ionization energy and cationic radius, and apparent activation entropies (ΔS app ‡ ) decrease 74 J mol –1 K –1 . These trends reflect electrostatic interactions that stabilize the cations to the anionic sites on the zeolite: stronger association between these charges leads to increasingly endothermic processes to displace the alkali metal to form a cationic methoxy and an intrapore metal lactate intermediate. Water physisorption measurements suggest alkali metal ions bind superstoichiometric quantities of water within FAU pores, and in situ infrared spectra suggest the concerted adsorption of ML requires reorganization of this water. Consequently, these processes introduce entropic gains that partially offset entropy losses associated with ML adsorption. Hence, turnover rates differ only by a factor of 2 among Na-, K-, and Cs-FAU at 573 K (ΔΔG app ‡ = 5 kJ mol –1 ). These findings demonstrate the interplay of alkali metal ions with zeolite active sites and intrapore water clusters for ML dehydration, indicating that these interactions can be leveraged to deliver optimal performance under different reaction conditions.

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

Entropy Splitting for High Order Numerical Simulation of Vortex Sound at Low Mach Numbers

A method of minimizing numerical errors, and improving nonlinear stability and accuracy associated with low Mach number computational aeroacoustics (CAA) is proposed. The method consists of two levels. From the governing equation level, we condition the Euler equations in two steps. The first step is to split the inviscid flux derivatives into a conservative and a non-conservative portion that satisfies a so called generalized energy estimate. This involves the symmetrization of the Euler equations via a transformation of variables that are functions of the physical entropy. Owing to the large disparity of acoustic and stagnation quantities in low Mach number aeroacoustics, the second step is to reformulate the split Euler equations in perturbation form with the new unknowns as the small changes of the conservative variables with respect to their large stagnation values. From the numerical scheme level, a stable sixth-order central interior scheme with a third-order boundary schemes that satisfies the discrete analogue of the integration-by-parts procedure used in the continuous energy estimate (summation-by-parts property) is employed.

Mueller, B.↗
DOE OSTI2021

The u 3 si 2 -H system

U 3 Si 2 is of interest to the nuclear industry as a candidate fuel material due to its high uranium density and high thermal conductivity. However, it has been observed to react with hydrogen, resulting in material decrepitation. As a result, it is important to understand the thermodynamics of the U 3 Si 2 -H system. In this study, the thermodynamics of the hydrogen absorption reaction of USi were determined experimentally using Sievert’s gas absorption and related to crystallographic evolution with hydrogen content using X-ray diffraction. Experimentally-determined thermodynamic parameters were compared with results from density functional theory modeling. Results from this study were also compared with those determined in previous work. Sievert’s gas absorption results were used to develop the pressure-composition-temperature (PCT) curves of the U 3 Si 2 -H system. It was found that the hydride phase exhibited a maximum stoichiometry between U 3 Si 2 H 1.8 and U 3 Si 2 H 2 . The two-phase region for hydride formation from U 3 Si 2 exhibited a miscibility gap with a critical temperature between 623 and 673 K, as calculated from the PCT curves. Analysis of the PCT curves also showed that both the enthalpy and entropy of the hydrogen absorption reaction increased with hydrogen content but were lower than the values for uranium trihydride formation from uranium metal. The enthalpy of reaction for hydrogen absorption was calculated to range between -86.9 and -94.8 kJ mol -1 , while the entropy of reaction was calculated to range between 101.9 and 138.8 J mol -1 K -1 . Furthermore, DFT modeling of the thermoydnamic stability of the USi hydride phases yielded a decomposition temperature of U 3 Si 2 H 2 of approximately 600 K, which was consistent with the experimental results. Similarly, the DFT-calculated enthalpy and entropy of reaction to form USiH were determined to be -106.5kJ mol and 121.8J mol -1 K -1 , respectively, which were both in close agreement with the experimentally-determined values.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

From phase decomposition to evaporation: A multi-modal evaluation of thermally degraded model lightweight high-entropy alloy

Lightweight high-entropy alloys (LHEAs) have the potential to replace conventional lightweight materials due to their superior mechanical properties and thermal stability. However, the thermal degradation pattern of LHEAs from phase decomposition to evaporation is not clear. We develop a new Al-based dual phase (FCC + HCP) LHEA—AlTi 0.45 CuZn, and further investigate its thermal degradation behavior for potential high-temperature structural applications. Using multimodal advanced characterization techniques such as differential scanning calorimetry/thermogravimetric analysis, scanning/transmission electron microscopy, and synchrotron X-ray diffraction/pair distribution function (XRD/PDF), a sequence of thermal degradation events beyond the thermal phase stability limit—between 250 and 360 °C—is observed. These include phase decomposition at ~360 °C, Zn evaporation at ~750 °C, and LHEA melting at 880 °C which results in ~25% cumulative weight loss. The formation of Al-Ti phase off the AlTi 0.45 CuZn matrix is due to the largest negative mixing enthalpy for Al-Ti than other binary pairs. Similarly, Zn evaporation from AlTi 0.45 CuZn LHEA is due to its faster evaporation rate than other constituent elements. The high-resolution synchrotron XRD and PDF results support the aforementioned observations; in addition, they reveal local atomic arrangements, local strain, and sluggish grain growth in the LHEA. Among other LHEAs of close density range (5.55 ≤ ρ ≤ 5.85 g/cc), the investigated LHEA exhibits outstanding nano-indentation hardness values due to the coupled grain size effect and HCP phase strengthening of the FCC matrix. As the search for LHEAs for lightweight applications grows, this study shows the potential use of AlTi 0.45 CuZn LHEA for structural applications even at elevated temperatures.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Computational Design of Cost-Effective High-Entropy Thermal/Environmental Barrier Coatings

Developing cost-effective thermal/environmental barrier coatings (TEBC) requires balance among various properties including low thermal conductivity, matching coefficient of thermal expansion (CTE), high thermal stability, high fracture toughness, and high recession resistance while being affordable. Low oxygen diffusivity is desirable as it can slow down oxygen transport to reach the underlying bond coating and hence delay oxidation of the bond coating. This project aims to design low-cost high-performance TEBC based on high entropy rare earth disilicates to protect SiC-based ceramic matrix composites from chemical and thermal attack for better performance of components in the hot section of gas turbine engines. To accelerate the TEBC design, we utilize first-principles density functional theory to predict key properties including phase stability, CTE, lattice thermal conductivity, temperature-dependent elastic constants, and oxygen diffusivity. Alloying elements including Yb, Y, Er, Eu, Gd, Lu, La, and Ce are considered, and modeling prediction are compared with available experimental results.

coefficient of thermal expansion↗
NASA NTRS2020

Toward Transition Modeling in a Hypersonic Boundary Layer at Flight Conditions

An accurate physics-based transition prediction method integrated with computational fluid dynamics (CFD) solvers is pursued for hypersonic boundary layer flows over slender hypersonic vehicles at flight conditions. The geometry and flow conditions are selected to match relevant trajectory locations from the ascent phase of the HIFiRE-1 flight experiment, namely, a 7-degree half-angle cone with 2.5 mm nose radius, freestream Mach numbers in the range of 3.8 – 5.5 and freestream unit Reynolds numbers in the range of 3.3 × 10(exp 6) – 21.4 × 10(exp 6) m(exp -1). Earlier research had shown that the onset of transition during the HIFiRE-1 flight experiment correlated with an amplification factor of N ≈ 13.5 for the planar Mack modes. However, to incorporate the N-factor correlations into a CFD code, we investigate surrogate models for disturbance amplification that avoid the direct computation of stability characteristics. A commonly used approach for low-speed flows is based on an a priori database of stability characteristics for locally similar profiles. However, the results presented in this paper demonstrate that the application of this approach to hypersonic boundary layers over blunt spherical nose-tip cones leads to large, unacceptable errors in the predictions of amplification factors, mainly due to its failure in accounting for the effects of the entropy layer on the boundary-layer profiles along the length of the model. We propose and demonstrate an alternate approach that employs the stability computations for a canonical set of blunt cone configurations to train a physics-informed convolutional neural network model that is shown to provide substantially improved transition predictions for hypersonic flow configurations with entropy-layer effects. Furthermore, the excellent performance of the neural network model is also confirmed for cone configurations with nose radius and half-angle values that do not correspond to those used to build the database. Finally, the convolutional neural network model is shown to outperform the linear stability calculations for underresolved basic states.

Pedro Paredes↗
DOE OSTI2023

High-entropy electrolytes for practical lithium metal batteries

Electrolyte engineering is crucial for improving battery performance, particularly for lithium metal batteries. Recent advances in electrolytes have greatly improved cyclability by enhancing electrochemical stability at the electrode interfaces, but concurrently achieving high ionic conductivity has remained challenging. Here we report an electrolyte design strategy for enhanced lithium metal batteries by increasing the molecular diversity in electrolytes, which essentially leads to high entropy electrolytes (HEEs). We find that in weakly solvating electrolytes, the entropy effect reduces ion clustering while preserving the characteristic anion-rich solvation structures, which is characterized by synchrotron-based X-ray scattering and molecular dynamics simulations. Electrolytes with smaller- sized clusters exhibit a 2-fold improvement in ionic conductivity compared to conventional weakly- solvating electrolytes, enabling stable cycling at high current densities up to 2C (6.2 mA cm -2 ) in anode- free LiNi 0.6 Mn 0.2 Co 0.2 (NMC622)||Cu pouch cells. In conclusion, the efficacy of the design strategy is verified by performance improvements in three disparate weakly solvating electrolyte systems.

25 ENERGY STORAGE↗
DOE OSTI2023

The stability of FeH x and hydrogen transport at Earth’s core mantle boundary

Iron hydride in Earth’s interior can be formed by the reaction between hydrous minerals (water) and iron. Studying iron hydride improves our understanding of hydrogen transportation in Earth’s interior. Our high-pressure experiments found that face-centered cubic (fcc) FeH x (x≤1) is stable up to 165 GPa, and our ab initio molecular dynamics simulations predicted that fcc FeH x transforms to a superionic state under lower mantle conditions. In the superionic state, H-ions in fcc FeH become highly diffusive-like fluids with a high diffusion coefficient of ~3.7 × 10 -4 cm 2 s -1 , which is comparable to that in the liquid Fe-H phase. The densities and melting temperatures of fcc FeH x were systematically calculated. Similar to superionic ice, the extra entropy of diffusive H-ions increases the melting temperature of fcc FeH. Further, the wide stability field of fcc FeH enables hydrogen transport into the outer core to create a potential hydrogen reservoir in Earth’s interior, leaving oxygen-rich patches (ORP) above the core mantle boundary (CMB).

58 GEOSCIENCES↗
DOE OSTI2020

Efficient, Thermally Stable, and Mechanically Robust All-Polymer Solar Cells Consisting of the Same Benzodithiophene Unit-Based Polymer Acceptor and Donor with High Molecular Compatibility

All-polymer solar cells (all-PSCs) are a highly attractive class of photovoltaics for wearable and portable electronics due to their excellent morphological and mechanical stabilities. Recently, new types of polymer acceptors (P A s) consisting of non-fullerene small molecule acceptors (NFSMAs) with strong light absorption have been proposed to enhance the power conversion efficiency (PCE) of all-PSCs. However, polymerization of NFSMAs often reduces entropy of mixing in PSC blends and prevents the formation of intermixed blend domains required for efficient charge generation and morphological stability. One approach to increase compatibility in these systems is to design P A s that contain the same building blocks as their polymer donor (P D ) counterparts. Here, a series of NFSMA-based P A s [P(BDT2BOY5-X), (X = H, F, Cl)] are reported, by copolymerizing NFSMA (Y5-2BO) with benzodithiophene (BDT), a common donating unit in high-performance P D s such as PBDB-T. All-PSC blends composed of PBDB-T PD and P(BDT2BOY5-X) P A show enhanced molecular compatibility, resulting in excellent morphological and electronic properties. Specifically, PBDB-T:P(BDT2BOY5-Cl) all-PSC has a PCE of 11.12%, which is significantly higher than previous PBDB-T:Y5-2BO (7.02%) and PBDB-T:P(NDI2OD-T2) (6.00%) PSCs. Additionally, the increased compatibility of these all-PSCs greatly improves their thermal stability and mechanical robustness. For example, the crack onset strain (COS) and toughness of the PBDB-T:P(BDT2BOY5-Cl) blend are 15.9% and 3.24 MJ m -3 , respectively, in comparison to the PBDB-T:Y5-2BO blends at 2.21% and 0.32 MJ m -3 .

14 SOLAR ENERGY↗
DOE OSTI2022

Energetically consistent model reduction for metriplectic systems

The metriplectic formalism is useful for describing complete dynamical systems which conserve energy and produce entropy. This creates challenges for model reduction, as the elimination of high-frequency information will generally not preserve the metriplectic structure which governs long-term stability of the system. Based on proper orthogonal decomposition, a provably convergent metriplectic reduced-order model is formulated which is guaranteed to maintain the algebraic structure necessary for energy conservation and entropy formation. Further, numerical results on benchmark problems show that the proposed method is remarkably stable, leading to improved accuracy over long time scales at a moderate increase in cost over naive methods.

42 ENGINEERING↗
DOE OSTI2024

Halting Oxygen Evolution to Achieve Long Cycle Life in Sodium Layered Cathodes

Oxygen redox chemistries at high voltage have materialized as a revolutionary paradigm for cathodes with high-energy density; however, they are plagued by the challenges of labile oxygen loss and rapid degradations upon cycling, even after concerted endeavors from the research community. Here we propose a multi-concentration stratagem propelled by entropy reinforcement to enhance the electronic structure disorder (ESD) at high desodiation states for impeding undesired oxygen mobility and ensuring controlled oxygen activity, elucidated by density functional theory calculations. The increased disorder strengthens the reversible electrochemistry of lattice oxygen redox, leading to effectively suppressed P−O structural evolution and highly stable localized TMO 6 octahedral environments, as demonstrated by soft/hard X-ray absorption spectroscopy. Furthermore, through a comparative analysis of sodium-layered cathodes with different configuration entropy, we reveal that a high-entropy state induced by cationic disordering has the capacity to perturb cationic redox boundaries, significantly restraining the formation of detrimental O′3 phases. As a consequence, the high-voltage cycling stability has been greatly upgraded, up to 4.4 V versus Na + /Na, with an impressive 90.1 % capacity retention at 1 C over 100 cycles and 76.1 % capacity retention at 2 C over 300 cycles. In conclusion, the resilient oxygen redox, enabled through the control of ESD, broadens the horizons for entropy engineering and lays the foundation for advancements in high-energy, long-cycling, and safe batteries.

25 ENERGY STORAGE↗
DOE OSTI2026

Creep behavior of high-entropy alloys: A critical review

High-entropy alloys (HEAs) have been extensively investigated during the last two decades. While substantial progress has been made in understanding their phase stability, microstructure, and deformation mechanisms at room and cryogenic temperatures, the long-term creep behavior (>100 h) of HEAs at high temperatures (>0.6 T m , where T m is the melting temperature) remains relatively underexplored. This knowledge gap is critical, as many engineering applications, including those for power generation and propulsion, require materials with good creep resistance to maintain structural integrity over extended service lifetimes. This review provides a focused and critical assessment of the current understanding of high-temperature deformation and creep behavior of HEAs, with particular attention paid to face-centered cubic HEAs and body-centered cubic refractory HEAs. The underlying deformation mechanisms governing their creep response and the influence of phase stability at elevated temperatures are examined in detail. Recent studies reveal mechanistic differences between HEAs and conventional dilute alloys that do not always lead to improved creep resistance belying their initial promise. Based on these findings, we discuss the challenges in designing HEAs for high-temperature structural applications and outline future research directions that may lead to creep-resistant HEAs.

Creep↗
DOE OSTI2022

The Two-Energies Turbulence Scheme Coupled to the Assumed PDF Method

An update of the two-energy turbulence scheme is presented, the 2TE + APDF scheme. The original version of the two-energy scheme is able to successfully model shallow convection without the need of an additional parameterization for non-local fluxes. However, the performance of the two-energy scheme is worse in stratocumulus cases, where it tends to overestimate the erosion of the stable layers. We have identified the causes: the non-local stability parameter does not consider local stratification, the scheme lacks an internal parameter that could distinguish between a shallow convection regime and a stratocumulus regime, and it uses an inflexible turbulence length scale formulation. To alleviate this problem, we propose several modifications: an update of the stability parameter, a modified computation of the turbulence length scale, and the introduction of the entropy potential temperature to distinguish between a shallow convection and a stratocumulus regime. In addition, the two-energy scheme is coupled to a simplified assumed probability density function method in order to achieve a more universal representation of the cloudy regimes. The updated turbulence scheme is evaluated for several idealized cases and one selected real case in the ICOsahedral Nonhydrostatic (ICON) modeling framework. The results show that the updated scheme corrects the overmixing problem in the stratocumulus cases. The performance of the updated scheme is comparable to the operational setup, and can be thus used instead of the operational turbulence and shallow convection scheme in ICON. Additionally, the updated scheme improves the coupling with dynamics, which is beneficial for the modeling of coherent flow structures in the atmospheric boundary layer.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

A prediction of the thermodynamic, thermophysical, and mechanical properties of CrTaO 4 from first principles

It is reported that the self-forming CrTaO 4 oxide scale can protect refractory high-entropy alloys from oxidation, superior to Cr 2 O 3 . In this paper, the phase stability, mechanical, and thermal properties of three polymorphous phases of CrTaO 4 are systematically investigated from first-principles density functional theory calculations. The mechanical properties predicted using the strain–energy methods indicated that all three phases are mechanically stable. The temperature dependence of elastic constants and polycrystalline moduli of three phases demonstrated the thermal softening as temperature increase. The Helmholtz-free energies as a function of volume and temperature are derived from phonon dispersions within the quasi-harmonic approximation at six strained volumes. The calculated apparent bulk coefficients of thermal expansion of these three phases are evaluated, the highest value approximately 13.4× 10 –6 K –1 within a temperature range of 500–2000 K for the rutile I4 1 md phase. The lattice thermal conductivity calculated by the Debye–Callaway model suggested that the rutile type I4 1 md phase has the lowest value of approximately 2.1 W/m/K at 1800 K. The other two phases, C2/m and P2/c, exhibit higher values due to relatively lower Grüneisen parameters and larger phonon velocities. The melting point of CrTaO 4 is predicted to be between 1975 and 2449 K using ab initio molecular dynamics simulations. Finally, this work provides a comprehensive theoretical understanding of the thermodynamic, mechanical, and thermal properties for the new material CrTaO 4 and serves as an example of a viable computational design strategy for improved oxidation resistance of refractory alloys at high temperatures.

36 MATERIALS SCIENCE↗