SEARCH · Search NASA

Results for “high enthalpy”

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 217 records · Page 12

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Wide-ranging predictions of new stable compounds powered by recommendation engines

The computational search for new stable inorganic compounds is faster than ever, thanks to high-throughput density functional theory (DFT). However, stable compound searches remain highly expensive because of the enormous search space and the cost of DFT calculations. To aid these searches, recommendation engines have been developed. We conduct a systematic comparison of the performance of previously developed recommendation engines, specifically ones based on elemental substitution, data mining, and neural network prediction of formation enthalpy. After identifying ways to improve the recommendation engines, we find the neural network to be superior at recommending stable Heusler compounds. Armed with improved recommendation engines, we identify tens of thousands of compounds that are stable at zero temperature and pressure, now available in the Open Quantum Materials Database. We summarize this diverse pool of compounds, including the elusive mixed anion compounds, and two of their many applications: thermoelectricity and solar thermochemical fuel production.

Science & Technology - Other Topics↗

DOE OSTI2021

Ab Initio Modeling of fcc Fe-Co-Cr-Ni High Entropy Alloys with Full Composition Range

The Fe-Co-Cr-Ni quaternary system has been studied extensively in the past decade, not only because of the superior properties achieved like high tensile ductility and fracture toughness, but also it is the foundation for the development of FeCrCoNi-based HEAs. However, most of the investigations are mainly focusing on the equiatomic and semi-equiatomic compositions. The physical properties of non-equiatomic compositions are barely explored. In the current work, the ab initio approach is adopted to predict the lattice parameter, structure stabilities, elastic properties, and enthalpy of formations of fcc Fe-Cr-Co-Ni HEA single crystals with the special quasirandom structure (SQS) method. To expand the design of Fe-Co-Cr-Ni HEAs into a full composition range, the current simulations cover all the binaries, ternaries, and the whole quaternary system. The 2-D composition-property contour diagrams are created to show the physical properties in the ternary and quaternary systems. Furthermore, our investigation shows that the fcc alloys' lattice parameter increases with the addition of Fe and Cr and decreases with Ni and Co. In addition to the enthalpy of formation at o K, the shape distortion rate is proposed to be another essential parameter to evaluate crystal stability. The results show the addition of Cr would destabilize the fcc lattice and cause symmetry breaking. Finally, the result of Pugh’s ratio shows that most of the fcc alloys show ductile behavior (P > 1.75), especially for the alloys around Fe 3 Ni, FeCo, and Fe 2 NiCo (P > 2.5). The brittle behavior (P < 1.75) is located around Fe 3 Cr and Fe 2 CoCr. Meanwhile, Cr 2 NiCo and Cr 2 FeNiCo are considered promising compositions to be tried and verified experimentally because of the higher bulk modulus and moderate shear modulus.

36 MATERIALS SCIENCE↗

DOE OSTI2025

An extension of the localized artificial diffusivity method for immiscible and high density ratio flows

The localized artificial diffusivity (LAD) method is widely regarded as the preferred multi-material regularization scheme for the compact finite difference method, because it is conservative, easy to implement, and generally robust for a wide range of multi-material problems. However, traditional LAD methods face significant challenges when applied to flows with large density ratios and when maintaining thermodynamic equilibrium across material interfaces. These limitations arise from the formulation of the artificial diffusivity flux and the reliance on enthalpy diffusion for interface regularization. Additionally, traditional LAD methods struggle to ensure stability under large density ratio conditions, fail to maintain a finite interface thickness, and are therefore unsuitable for modeling immiscible interfaces. Here, in this work, we discuss the origins of these issues in traditional LAD methods and propose modifications which enable the simulation of large density ratio and immiscible flows. The proposed method targets the artificial diffusion fluxes at gradients and ringing in the volume fraction, rather than the mass fraction in traditional methods, to consistently regularize large density ratio interfaces. Furthermore, the proposed method introduces an artificial bulk density diffusion term to enforce equilibrium conditions across interfaces. To address the challenge of modeling immiscible flows, a conservative diffuse interface term is incorporated into the formulation to ensure a finite interface thickness. Specific consideration is taken in the design of the method to ensure that these crucial properties are maintained for N -material flows. The effectiveness of the proposed method is demonstrated through a series of canonical test cases, and its accuracy is validated by comparison with experimental data on micro-bubble collapse in water. These results highlight the method’s robustness and its ability to overcome the limitations of traditional LAD approaches.

Artificial diffusivity↗

DOE OSTI2023

Fundamentals of SrCoO 3-δ Based Oxygen Deficient Perovskites as Cathodes for Solid Oxide Fuel Cells

The local structure and oxygen stoichiometry in oxide systems have a profound impact on oxygen electrocatalysis encountered in metal-air batteries and solid oxide fuel cells (SOFCs). However, this knowledge is often obtained under conditions different from the real working conditions of the material, resulting in misinterpretation and misunderstanding. This PhD dissertation aims to obtain the structure and oxygen-stoichiometry information of a class of perovskite oxides under their real working conditions in solid oxide fuel cells. Several perovskite oxides were selected for the study: Sr 0.9 Y 0.1 CoO 3-δ (SYC10), Sr 0.9 Y 0.3 CoO 3-δ (SYC30), SrCo 0.9 Nb 0.1 O 3-δ (SCN10) and SrCo 0.9 Ta0.1O 3-δ (SCT10). The local crystal structure and oxygen stoichiometry of these materials were systematically characterized with in-situ neutron diffraction (ND). The oxygen stoichiometry was also measured by thermogravimetric analysis (TGA) and iodometric titration methods. To establish the correlation among structure, oxygen stoichiometry and oxygen reduction reaction (ORR) activity, high temperature electrochemical impedance spectroscopy, oxygen permeation and single cell testing were also performed. In chapter 2, SYC10 is described to have a more symmetrical structure and $V_o^{··}$ distribution, and higher $V_o^{··}$ concentration than SYC30. Molecular orbital energy analysis based on the local structure of the ORR-active Co1-polyhedra indicates SYC30 has a higher Fermi level relative to O-2p energy level in the active Co1-polyhedra, and thus a higher motional enthalpy for $V_o^{··}$ migration. However, the ORR activity of Y-doped SrCoO 3-δ (SCO) are found uncompetitive with other popular catalysts. In chapter 3 and 4, systematic structural studies are presented for Nb and Ta-doped SCO, namely SCN10 and SCT10, respectively. The methodology of unrevealing the structure-activity relationships in chapter 2 is also found applicable to SYC10, SCN10 and SCT10 as presented in chapter 5. The electrochemical tests indicate that SCT10 is the best catalyst due to its excellent thermal and electrochemical stability. In chapter 6, SCT10 is evaluated in a real solid oxide fuel cell in the form of nanoscaled SCT10 layer coating on the commercial (La 0.6 Sr 0.4 ) 0.95 Co 0.2 Fe 0.8 O 3-δ (LSCF). The bilayer structured cathode shows excellent properties of coarsening-resistant, Sr-segregation-free and high ORR-activity.

Yang, Tianrang↗

DOE OSTI2024

Influence of Carboxymethyl Cellulose as a Thickening Agent for Glauber’s Salt-Based Low Temperature PCM

This work is focused on a novel, promising low temperature phase change material (PCM), based on the eutectic Glauber’s salt composition. To allow phase transition within the refrigeration range of temperatures of +5 °C to +12 °C, combined with a high repeatability of melting–freezing processes, and minimized subcooling, the application of three variants of sodium carboxymethyl cellulose (Na-CMC) with distinct molecular weights (700,000, 250,000, and 90,000) is considered. The primary objective is to optimize the stabilization of this eutectic PCM formulation, while maintaining the desired enthalpy level. Preparation methods are refined to ensure repeatability in mixing components, thereby optimizing performance and stability. Additionally, the influence of Na-CMC molecular weight on stabilization is examined through differential scanning calorimetry (DSC), T-history, and rheology tests. The PCM formulation of interest builds upon prior research in which borax, ammonium chloride, and potassium chloride were used as additives to sodium sulfate decahydrate (Glauber’s salt), prioritizing environmentally responsible materials. The results reveal that CMC with molecular weights of 250 kg/mol and 90 kg/mol effectively stabilize the PCM without phase separation issues, slowing crystallization kinetics. Conversely, CMC of 700 kg/mol proved ineffective due to the disruption of gel formation at its low gel point, hindering higher concentrations. Calculations of ionic concentration indicate higher Na ion content in PCM stabilized with 90 kg/mol CMC, suggesting increased ionic interactions and gel strength. A tradeoff is discovered between the faster crystallization in lower molecular weight CMC and the higher concentration required, which increases the amount of inert material that does not participate in the phase transition. After thermal cycling, the best formulation had a latent heat of 130 J/g with no supercooling, demonstrating excellent performance. This work advances PCM’s reliability as a thermal energy storage solution for diverse applications and highlights the complex relationship between Na-CMC molecular weight and PCM stabilization.

25 ENERGY STORAGE↗

DOE OSTI2023

Thermodynamic properties and superconductivity of natural carrollite (CuCo 2 S 4 )

The thermodynamic properties of natural carrollite with a formula Cu 0.92 Co 2.07 S 4 are studied using low temperature heat capacity and high temperature oxidative solution calorimetry. The standard thermodynamic functions at 298.15 K are: C p,m °, S m °, H m °, and Φ m ° are 158.48 J∙K −1 ∙mol −1 , 176.33 J∙K −1 ∙mol −1 , 28.40 kJ∙mol −1 , and 81.07 kJ∙mol −1 with an estimated error of 1%. The enthalpies (Δ f H), entropies (Δ f S), and Gibbs energies of formation (Δ f G) from elements at 298.15 K are: −344.46 ± 12.87 kJ∙mol −1 , −45.22 ± 1.81 J∙K −1 ∙mol −1 , and −331.14 ± 1.89 J∙K −1 ∙mol −1 . The stability of carrollite relative to the elements is demonstrated with a negative Gibbs energy (ΔG r °) of formation for several temperatures between 0 and 300 K. Furthermore, a superconductivity transition, which has been observed previously for synthetic carrollite, is confirmed by both the low temperature heat capacity and magnetization measurements.

Chemistry↗

DOE OSTI2024

Thermodynamic modeling of CsF with LiF-NaF-KF for molten fluoride-fueled reactors

Gibbs energy models were developed to describe the thermochemical behavior of CsF in molten FLiNaK (46.5LiF-11.5NaF-42KF mol%), a proposed molten salt reactor (MSR) fuel solvent and coolant, as cesium is of concern due to its high radiotoxicity and volatility. Initially, it was necessary to obtain a more accurate Gibbs energy function for CsF which required fitting parameters to reported vapor pressures over condensed phase CsF. The pseudo-binary systems CsF-LiF, CsF-NaF and CsF-KF were then evaluated utilizing phase equilibria and enthalpy of mixing (Δ mix H) values, together with original differential scanning calorimetry (DSC) measurements performed for the CsF-LiF and CsF-KF systems. The CsF-LiF-NaF, CsF-LiF-KF and CsF-NaF-KF pseudo-ternary system representations were obtained by interpolation of the constituent pseudo-binary systems, with DSC measurements performed for the CsF-LiF-NaF system to corroborate the calculated liquidus temperature. Ultimately, the pseudo-ternary systems were interpolated to obtain Gibbs energy models for the pseudo-quaternary CsF-LiF-NaF-KF system, supported by DSC measurements at low CsF compositions (1–10 mol%), yielding computed equilibria and cesium-containing vapor pressures that compare favorably with reported values. In conclusion, the Molten Salt Thermal Properties Database – Thermochemical (MSTDB-TC) was subsequently expanded to include these Gibbs energy models allowing description of the thermochemical behavior of the CsF-LiF-NaF-KF system.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2026

Computational Discovery of Ultralow Thermal Conductivity in the Energy-Degenerate Polymorphic Crystal Family A 2 M 2 M’Q 4

Crystalline materials, characterized by their well-defined lattices, typically exhibit a unique global thermodynamic minimum for a specific composition. However, in this study, we discover a quaternary chalcogenide family, A 2 M 2 M’Q 4 (A: alkali metals; M: coinage metal; M’: transition or group-IVA metals; Q: chalcogens), that exhibits pervasive energy (near-)degeneracy. For a given composition, multiple structurally distinct polymorphs exist within a formation enthalpy window of only a few milli-electron volts per atom. We quantify this inherent structural flexibility using a dedicated descriptor, σ f : the standard deviation of formation enthalpies among degenerate (meta)stable polymorphs. The consistently low σf observed across the A 2 M 2 M’Q 4 family signifies a characteristically shallow and frustrated potential energy landscape, which drives pronounced lattice anharmonicity, marking these materials as prime candidates for ultralow lattice thermal conductivity (κ L ). Employing an advanced high-throughput computational framework that integrates thermodynamics, lattice dynamics, and thermal conductivity calculations, we screen 1215 A 2 M 2 M’Q 4 compounds, identifying 30 stable candidates with κ L < 0.5 W m –1 K –1 at 300 K. Among them, Rb 2 Ag 2 SnTe 4 and Rb 2 Au 2 HfTe 4 , two representatives from the IVA and TM subgroups, are predicted to show ultralow room-temperature κ L of 0.174 W m –1 K –1 and 0.295 W m –1 K –1 , respectively. A systematic analysis suggests that the nonbonding and antibonding states induced by “dual rattlers” are the origin of low thermal conductivity in these compounds. Our results position the A 2 M 2 M’Q 4 family as a rich source of intrinsic thermal insulators and suggest that polymorphic energy degeneracy may serve as a valuable signpost for identifying crystalline families with potential anharmonicity.

cations↗

DOE OSTI2025

Atomic-Scale Insights into Carbon Dissolution in α-, γ-, and θ-Al 2 O 3 : Phase-dependent Transport Dynamics from First-Principles Calculations

α-Al 2 O 3 exhibits superior carburizing corrosion resistance compared to metastable γ-Al 2 O 3 and θ-Al 2 O 3 phases in high-temperature CO 2 environments, yet its atomic-scale origins remain unclear. Using first-principles density functional theory, we systematically investigate carbon dissolution and diffusion in α-Al 2 O 3 , γ-Al 2 O 3 , and θ-Al 2 O 3 , including the effects of oxygen (O) and aluminum (Al) vacancies. Our results show that α-Al 2 O 3 consistently exhibits higher carbon solution enthalpies than γ-Al 2 O 3 and θ-Al 2 O 3 in both pristine and defective structures, indicating lower intrinsic carbon solubility in α-Al 2 O 3 . Vacancies significantly enhance carbon incorporation: O vacancies reduce solution enthalpy, while Al vacancies further amplify this effect, with a strong preference for carbon at Al vacancy sites. Carbon diffusion barriers are also highest in α-Al 2 O 3 , reflecting slower carbon mobility. Al vacancies increase diffusion barriers across all phases, while O vacancies raise barriers in α- and γ-Al 2 O 3 but slightly lower them in θ-Al 2 O 3 . These results reveal a dual mechanism behind the carburizing resistance of α-Al 2 O 3 : reduced carbon solubility and elevated diffusion barriers. Furthermore, this work provides atomic-scale insights to guide the design of alumina-based materials with improved carburizing resistance through phase selection and defect engineering.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Technoeconomic analysis of hydrogen storage using 1,4-butanediol (BDO)/γ-butyrolactone (GBL) as a stationary backup power system

Liquid organic hydrogen carriers (LOHCs) are compounds that store and release hydrogen in stable forms at high density. While one-way carriers such as methanol and ammonia have gained attention, their economic advantages are often realized from their use as an export product and direct use as a fuel. Liquid carrier materials that can instead be cycled for energy storage have promise for stationary power applications. In particular, the reversible LOHC system 1,4-butanediol (BDO, H 2 -rich) and gamma-butyrolactone (GBL, H 2 -lean) has a lower enthalpy of dehydrogenation compared to conventional cyclic hydrocarbons and can utilize non-precious metal copper-based catalysts. Here, in this study, BDO/GBL system capital and operating expenses are characterized for vapor phase versus liquid phase hydrogenation and dehydrogenation in a 10 MW backup power application corresponding to sizing of Tier 2 datacenters as well as other critical infrastructure such as hospitals. Costs are benchmarked against two incumbent technologies: a well-established methylcyclohexane/toluene carrier system and compressed gas storage. BDO/GBL storage costs are found to differ substantially between operating modes, with liquid phase hydrogenation coupled with liquid phase dehydrogenation leading to the lowest LCOS of $\$$4.58/kg H 2 in the absence of byproduct formation. In this bounding case, LCOS for the BDO/GBL system is lower than for MCH/TOL ($\$$6.97/kg H 2 ) and compressed gas ($\$$8.48/kg H 2 at 170 bara and $\$$12.05/kg H 2 at 350 bara). However, escalating costs of carrier replacement due to byproduct formation (ranging from an added $\$$6–11/kg H 2 ) illustrate the need for highly selective catalysts to ensure BDO/GBL carrier viability.

BDO/GBL↗

DOE OSTI2022

Overcoming significant challenges in extracting off-stoichiometric thermodynamics using the compound energy formalism through complementary use of experimental and first principles data: A case study of Ba 1-x Sr x FeO 3-δ

The compound energy formalism (CEF) is a powerful framework to describe the thermodynamics of metal oxides as a function of off-stoichiometry, temperature, and composition. The thermodynamic properties are crucial materials design attributes in metal oxide-based oxygen-exchange chemical processes. Despite the richness of information an accurate CEF model provides, a method to determine a unique and accurate fit for oxygen exchange materials remains elusive. This contribution details a method for fitting the CEF model that overcomes the current fitting challenges through three innovations: 1) the combination of density functional theory calculations with experimental data decorrelates excess terms and delineates the enthalpic/entropic contributions to the Gibbs free energy; 2) a systematic determination of the important CEF model terms, removing thermodynamic predetermining human intervention; 3) a self-consistent solution of the starting oxygen offstoichiometry (δ 0 ) of thermogravimetric measurements. Thus, our method enables the reliable extraction of off-stoichiometric metal oxide thermodynamic properties and facilitates rapid materials compositional screening, and reliable process design of systems dependent on off-stoichiometric redox-active metal oxides. We apply this method to a Ba x Sr 1-x FeO 3-δ test case. We find by systematically examining the performance of the CEF model fit with and without each innovation that all three innovations are necessary for an accurate fit. We determined that reduction enthalpy is higher and more sensitive to off-stoichiometry when the Sr fraction is large (139.5 and 185.3 kJ/mol O 2 for SrFeO 3 at δ = 0 and δ = 0.5, respectively vs. nearly constant 83 kJ/mol O 2 for BaFeO 3 ). However, the reduction entropy is mostly insensitive to Sr fraction, but highly dependent on δ suggesting larger contributions in the non-configurational entropy parameters.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Process Engineering to Increase the Layered Phase Concentration in the Immediate Products of Flame Spray Pyrolysis

Flame-spray-pyrolysis (FSP) is a robust and scalable process to synthesize particles at the commodity-scale. FSP has been used to produce the precursor powders which were converted to the layered structure $(R\bar{3m} \text { phase})$ by a postannealing step in making nickel-rich cathode materials (NCMs). Theoretically, the high flame temperature (normally >1500 K) in FSP can provide adequate energy for the phase conversion from rock-salt to layered structures and potentially enables one-step synthesis. However, the high flame temperature is a critical issue to cause lithium loss and structural degradation, preventing the formation of the layered phase. Here, guided by the gaseous nucleation theory, we implemented several FSP processes with different solution recipes. The layered phase concentration in the as-burned products can be increased with the solution enthalpies. By adding a rapid quench step to suppress the lithium loss and phase degradation, the layered phase can be further increased. This work contributes new ideas to innovating process regarding the process efficiency and throughput of manufacturing cathode materials at a large scale.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Structural phase transition, s ± -wave pairing, and magnetic stripe order in bilayered superconductor La 3 Ni 2 O 7 under pressure

Motivated by the recently discovered high-T c superconductor La 3 Ni 2 O 7 , we comprehensively study this system using density functional theory and random phase approximation calculations. At low pressures, the Amam phase is stable, containing the Y 2- mode distortion from the Fmmm phase, while the Fmmm phase is unstable. Because of small differences in enthalpy and a considerable Y 2- mode amplitude, the two phases may coexist in the range between 10.6 and 14 GPa, beyond which the Fmmm phase dominates. In addition, the magnetic stripe-type spin order with wavevector (π, 0) was stable at the intermediate region. Pairing is induced in the s ± -wave channel due to partial nesting between the M = (π, π) centered pockets and portions of the Fermi surface centered at the X = (π, 0) and Y = (0, π) points. This resembles results for iron-based superconductors but has a fundamental difference with iron pnictides and selenides. Moreover, our present efforts also suggest La 3 Ni 2 O 7 is qualitatively different from infinite-layer nickelates and cuprate superconductors.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Effects of Charged Solute-Solvent Interaction on Reservoir Temperature during Subsurface CO2 Injection

A short-term side-effect of CO2 injection is a developing low-pH front that forms ahead of the bulk water injectant, due to differences in solute diffusivity. Observations of downhole well temperature show a reduction in aqueous-phase temperature with the arrival of a low-pH front, followed by a gradual rise in temperature upon the arrival of a high concentration of bicarbonate ion. In this work, we model aqueous-phase transient heat advection and diffusion, with the volumetric energy generation rate computed from solute-solvent interaction using the Helgeson–Kirkham–Flowers (HKF) model, which is based on the Born Solvation model, for computing specific molar heat capacity and the enthalpy of charged electrolytes. A computed injectant water temperature profile is shown to agree with the actual bottom hole sampled temperature acquired from sensors. The modeling of aqueous-phase temperature during subsurface injection simulation is important for the accurate modeling of mineral dissolution and precipitation because forward dissolution rates are governed by a temperature-dependent Arrhenius model.

58 GEOSCIENCES↗

DOE OSTI2025

Quantum mechanical dataset of 836k neutral closed-shell molecules with up to 5 heavy atoms from C, N, O, F, Si, P, S, Cl, Br

Abstract We introduce the Vector-QM24 (VQM24) dataset comprehensively covering all possible neutral closed-shell small organic and inorganic molecules with up to five heavy (p-block) atoms: C, N, O, F, Si, P, S, Cl, Br. All valid stoichiometries, Lewis-rule-consistent graphs, and stable conformers (identified via GFN2-xTB) were enumerated combinatorially, yielding 577k conformational isomers spanning 258k constitutional isomers and 5,599 unique stoichiometries. DFT (ωB97X-D3/cc-pVDZ) optimizations were performed for all, and diffusion quantum Monte Carlo (DMC@PBE0(ccECP/cc-pVQZ)) energies are provided for 10,793 lowest-energy conformers with up to 4 heavy atoms. VQM24 includes structures, vibrational modes, rotational constants, thermodynamic properties (Gibbs free energies, enthalpies, ZPVEs, entropies, heat capacities), and electronic properties such as atomization, electron interaction, exchange-correlation, dispersion energies, multipole moments (dipole to hexadecapole), alchemical potentials, Mulliken charges, and wavefunctions. Machine learning models of atomization energies on this dataset reveal significantly higher complexity than QM9, with none achieving chemical accuracy. VQM24 offers a rigorous, high-fidelity benchmark for evaluating quantum machine learning models.

Science & Technology - Other Topics↗

DOE OSTI2020

Development of a CPU/GPU portable software library for Lagrangian–Eulerian simulations of liquid sprays

The Lagrangian–Eulerian method is widely used for simulations of fuel sprays in turbulent combustion because of the advantage of treating the spray droplets as discrete points. One challenge of the Lagrangian–Eulerian method is the intense computational requirement when tracking the large number of Lagrangian particles needed for high fidelity. We have developed a performance-portable library, Grit, to track the Lagrangian particles in parallel on central processing unit (CPU) and graphics processing unit (GPU) accelerated high performance computing (HPC) architectures. Grit is a C++ library which employs Message Passing Interface (MPI) for distributed memory parallelism and Kokkos programming model for on-node shared memory parallelism with performance portability across different architectures of GPUs and multi-core/manycore CPUs. The parallel algorithms, key parallel kernels, and their performances on the pre-exascale supercomputer, Summit, are presented. Grit is coupled with a direct numerical simulation (DNS) solver, S3D, for multiphase simulations. A conservative formulation has been developed and implemented in Grit for phase coupling with thermodynamic consistency. The formulation separates the conservation of mass, momentum and energy from the physical models to prevent accidental violation of conservation laws due to inconsistent models. The formulation also enforces consistent definitions of enthalpies of the fuel for both phases and the latent heat of evaporation. Finally, simulations of turbulent particle-laden flow and the evaporation of dilute turbulent spray jet are performed to verify the software implementation, and to demonstrate the scalability of Grit for large-scale multiphase simulations.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2022

Thermodynamics of the Nonrelativistic Free-Electron Fermi Gas in One, Two, and Three Dimensions from the Degenerate to the Nondegenerate Temperature Regime

The thermodynamic properties of a nonrelativistic free-electron Fermi gas is of fundamental interest in condensed matter physics. Properties previously studied in three-dimensions (3D) in the lowand high-temperature limits include the internal energy, heat capacity, zero-field magnetic spin susceptibility, and pressure. Here we report solutions for the temperature dependence spanning these two temperature regimes of the chemical potential, internal energy, magnetic susceptibility, and the heat capacity at constant volume in 1D, 2D, and 3D. Also calculated are the pressure, enthalpy, heat capacity at constant pressure, isothermal compressibility, and thermal expansion coefficient versus temperature in 2D and 3D. Of primary interest here are the detailed dimension-dependent crossovers of these properties between the degenerate and nondegenerate temperature regime, which are graphically illustrated for each of the above properties.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Development of Precipitation-Strengthened Al 0.8 NbTiVM (M = Co, Ni) Light-Weight Refractory High-Entropy Alloys

Single-phase solid-solution refractory high-entropy alloys (RHEAs) have been receiving significant attention due to their excellent mechanical properties and phase stability at elevated temperatures. Recently, many studies have been reported regarding the precipitation-enhanced alloy design strategy to further improve the mechanical properties of RHEAs at elevated temperatures. In this study, we attempted to develop precipitation-hardened light-weight RHEAs via addition of Ni or Co into Al 0.8 NbTiV HEA. The added elements were selected due to their smaller atomic radius and larger mixing enthalpy, which is known to stimulate the formation of precipitates. The addition of the Ni or Co leads to the formation of the sigma precipitates with homogeneous distribution. The formation and homogeneous distribution of sigma particles plays a critical role in improvement of yield strength. Furthermore, the Al 0.8 NbTiVM 0.2 (M = Co, Ni) HEAs show excellent specific yield strength compared to single-phase AlNbTiV and NbTiVZr RHEA alloys and conventional Ni-based superalloy (Inconel 718) at elevated temperatures.

36 MATERIALS SCIENCE↗