Search NASA⌕ Search

SEARCH · Search NASA

Results for “predictive thermodynamics”

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 235 records · Page 13

Low-Pressure Electrochemical Synthesis of Complex High-Pressure Superconducting Superhydrides

There is great current interest in multicomponent superhydrides due to their unique quantum properties under pressure. A remarkable example is the ternary superhydride Li 2 MgH 16 computationally identified to have an unprecedented high superconducting critical temperature T c of ~470 K at 250 GPa. However, the very high synthesis pressures required remains a significant hurdle for detailed study and potential applications. In this Letter, we evaluate the feasibility of synthesizing ternary Li-Mg superhydrides by the recently proposed pressure-potential (P 2 ) method that uniquely combines electrochemistry and applied pressure to control synthesis and stability. Furthermore, the results indicate that it is possible to synthesize Li-Mg superhydrides at modest pressures by applying suitable electrode potentials. Using pressure alone, no Li-Mg ternary hydrides are predicted to be thermodynamically stable, but in the presence of electrode potentials, both Li 2 MgH 16 and Li 4 MgH 24 can be stabilized at modest pressures. Three polymorphs are predicted as ground states of Li 2 MgH 16 below 300 GPa, with transitions at 33 and 160 GPa. The highest pressure phase is superconducting, while the two at lower pressures are not. Our findings point out the potentially important role of the P 2 method in controlling phase stability of complex multicomponent superhydrides.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

SiC and Si3N4 Recession Due to SiO2 Scale Volatility Under Combustor Conditions

Silicon carbide (SiC) and Si3N4 materials were tested in various turbine engine combustion environments chosen to represent either conventional fuel-lean or fuel-rich mixtures proposed for high-speed aircraft. Representative chemical vapor-deposited (CVD), sintered, and composite materials were evaluated by furnace and high-pressure burner rig exposures. Although protective SiO2 scales formed in all cases, the evidence presented supports a model based on paralinear growth kinetics (i.e., parabolic growth moderated simultaneously by linear volatilization). The volatility rate is dependent on temperature, moisture content, system pressure, and gas velocity. The burner tests were thus used to map SiO2 volatility (and SiC recession) over a range of temperatures, pressures, and velocities. The functional dependency of material recession (volatility) that emerged followed the form A[exp(-Q / RT)](P(sup x)v(sup y). These empirical relations were compared with rates predicted from the thermodynamics of volatile SiO and SiOxHy reaction products and a kinetic model of diffusion through a moving boundary layer. For typical combustion conditions, recession of 0.2 to 2 micrometers/hr is predicted at 1200 to 1400 C, far in excess of acceptable long-term limits.

Smialek, James L.↗

Nitromethane Decomposition via Automated Reaction Discovery and an Ab Initio Corrected Kinetic Model

In the explore the systematic construction of kinetic models from in silico reaction data for the decomposition of nitromethane. Our models are constructed in a computationally affordable manner by using reactions discovered through accelerated molecular dynamics simulations using the ReaxFF reactive force field. The reaction paths are then optimized to determine reaction rate parameters. We introduce a reaction barrier correction scheme that combines accurate thermochemical data from density functional theory with ReaxFF minimal energy paths. We validate our models across different thermodynamic regimes, showing predictions of gas phase CO and NO concentrations and high-pressure induction times that are similar to experimental data. The kinetic models are analyzed to find fundamental decomposition reactions in different thermodynamic regimes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Metal oxide candidates for thermochemical water splitting obtained with a generative diffusion model

Generative diffusion models (DMs) for inorganic crystalline materials are being actively investigated for their potential to expand the chemical and structural design spaces for known functional materials. Generative candidates are particularly useful for applications where few functional, let alone commercially viable, materials currently exist, such as metal oxides for thermochemical water-splitting, which have strict requirements for defect thermodynamics and host stability. Here, we critically examine generated metal oxides from the M ATTER G EN DM conditioned on select chemical systems for thermochemical water splitting applications. Perhaps most notably, we find that M ATTER G EN predicts a novel, thermodynamically stable, quinary metal oxide, Ba 2 SrInFeO 6 , although this compound represents an ordered and layered substitution within the same A 3 B 2 O 6 structural prototype as its two ternary end members. Detailed density functional theory calculations and spin configuration sampling for this material and its possible decomposition products—beyond what existed in M ATTER G EN training data—are required to quantitatively validate hull energy predictions and conclusions of stability. Furthermore, the material exhibits oxygen defect formation energies appropriate for thermochemical water splitting, warranting targeted investigation in an experimental validation campaign, along with other future M ATTER G EN candidates in this application space.

36 MATERIALS SCIENCE↗

The solubility of NdPO 4 and DyPO 4 and stability of Nd and Dy chloride and hydroxyl complexes as a function of pH and salinity to 450 °C

The mobility of rare earth elements (REE) in geological systems is often controlled by the stability of monazite and xenotime, which are important hosts for the light and heavy REE. These REE phosphates constitute important resources and provide information on ore formation conditions and timing due to their uses as geothermometers and geochronometers. While the thermodynamic properties for these minerals are well-established up to high temperature and pressure, the properties for the aqueous REE complexes are not well constrained with limited information up to 300 °C. In this study, the solubility of synthetic NdPO 4 and DyPO 4 endmembers were measured in hydrothermal sub- to supercritical NaCl-bearing aqueous solutions from 350 at saturated water vapor pressure to 450 °C at 700 bar. The speciation of Nd and Dy was determined in acidic to alkaline solutions (pH 25 °C values from 2 to 10), and indicates that the hydroxyl species REE(OH) 2 + and REE(OH) 3 0 predominate at low salinity (0.01 mol/kg NaCl) with some contribution of REE chloride species REECl 2+ and REECl 2 + in acidic fluids. The REE phosphate solubility measured in these experiments is up to two orders of magnitude higher than predicted using existing thermodynamic properties from literature for aqueous species extrapolated from lower temperature data to supercritical conditions. To address this discrepancy, the thermodynamic properties of the REE aqueous species were optimized using GEMSFITS to derive the formation constants for the REE chloride (β Cl ) and hydroxyl (β OH ) complexes in the studied temperature range. This study highlights the importance of the hydroxyl species REE(OH) 2 + and REE(OH) 3 0 over a wide range of pH and temperature. As a result, the revised thermodynamic properties for Dy and Nd chloride and hydroxyl species more accurately predict the solubility of NdPO 4 and DyPO 4 and provide insight into the speciation and mobility of the light and heavy REE in hydrothermal supercritical fluids in the crust.

58 GEOSCIENCES↗

Extreme Metastability of Diamond and its Transformation to the BC8 Post-Diamond Phase of Carbon

Diamond possesses exceptional physical properties due to its remarkably strong carbon–carbon bonding, leading to significant resilience to structural transformations at very high pressures and temperatures. Despite several experimental attempts, synthesis and recovery of the theoretically predicted post-diamond BC8 phase remains elusive. Through quantum-accurate multimillion atom molecular dynamics (MD) simulations, we have uncovered the extreme metastability of diamond at very high pressures, significantly exceeding its range of thermodynamic stability. We predict the post-diamond BC8 phase to be experimentally accessible only within a narrow high pressure–temperature region of the carbon phase diagram. The diamond to BC8 transformation proceeds through premelting followed by BC8 nucleation and growth in the metastable carbon liquid. Here we propose a double-shock compression pathway for BC8 synthesis, which is currently being explored in experiments at the National Ignition Facility.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

SiC and Si3N4 Recession Due to SiO2 Scale Volatility Under Combustor Conditions

SiC and Si3N4 materials were tested under various turbine engine combustion environments, chosen to represent either conventional fuel-lean or fuel-rich mixtures proposed for high speed aircraft. Representative CVD, sintered, and composite materials were evaluated in both furnace and high pressure burner rig exposure. While protective SiO2 scales form in all cases, evidence is presented to support paralinear growth kinetics, i.e. parabolic growth moderated simultaneously by linear volatilization. The volatility rate is dependent on temperature, moisture content, system pressure, and gas velocity. The burner tests were used to map SiO2 volatility (and SiC recession) over a range of temperature, pressure, and velocity. The functional dependency of material recession (volatility) that emerged followed the form: exp(-QIRT) * P(exp x) * v(exp y). These empirical relations were compared to rates predicted from the thermodynamics of volatile SiO and SiO(sub x)H(sub Y) reaction products and a kinetic model of diffusion through a moving, boundary layer. For typical combustion conditions, recession of 0.2 to 2 micron/h is predicted at 1200- 1400C, far in excess of acceptable long term limits.

Smialek, James L.↗

Mixed-Chalcogen 2D Silver Phenylchalcogenides (AgE 1– x E x Ph; E = S, Se, Te)

Alloying is a powerful strategy for tuning the electronic band structure and optical properties of semiconductors. Here, we investigate the thermodynamic stability and excitonic properties of mixed-chalcogen alloys of two-dimensional (2D) hybrid organic-inorganic silver phenylchalcogenides (AgEPh; E = S, Se, Te). Using a variety of structural and optical characterization techniques, we demonstrate that the AgSePh-AgTePh system forms homogeneous alloys (AgSe 1-x Te x Ph, 0 ≤ x ≤ 1) across all compositions, whereas the AgSPh-AgSePh and AgSPh-AgTePh systems exhibit distinct miscibility gaps. Density functional theory calculations reveal that chalcogen mixing is energetically unfavorable in all cases, but comparable in magnitude to the ideal entropy of mixing at room temperature. Because AgSePh and AgTePh have the same crystal structure (which is different from AgSPh), alloying is predicted to be thermodynamically preferred over phase separation in the case of AgSePh-AgTePh, whereas phase separation is predicted to be more favorable than alloying for both the AgSPh-AgSePh and AgSPh-AgTePh systems, in agreement with experimental observations. Homogeneous AgSe1-xTexPh alloys exhibit continuously tunable excitonic absorption resonances in the ultraviolet-visible range, while the emission spectrum reveals competition between exciton delocalization and self-trapping behavior. Altogether, these observations provide new insight into the thermodynamics of 2D silver phenylchalcogenides and the effect of lattice composition on electron-phonon interactions in 2D hybrid organic-inorganic semiconductors.

36 MATERIALS SCIENCE↗

A graph neural network-state predictive information bottleneck (GNN-SPIB) approach for learning molecular thermodynamics and kinetics

Molecular dynamics simulations offer detailed insights into atomic motions but face timescale limitations. Enhanced sampling methods have addressed these challenges but even with machine learning, they often rely on pre-selected expert-based features. Here, in this work, we present a Graph Neural Network-State Predictive Information Bottleneck (GNN-SPIB) framework, which combines graph neural networks and the state predictive information bottleneck to automatically learn low-dimensional representations directly from atomic coordinates. Tested on three benchmark systems, our approach predicts essential structural, thermodynamic and kinetic information for slow processes, demonstrating robustness across diverse systems. The method shows promise for complex systems, enabling effective enhanced sampling without requiring pre-defined reaction coordinates or input features.

Zou, Ziyue↗

TankSIM: A Cryogenic Tank Performance Prediction Program

Accurate prediction of the thermodynamic state of the cryogenic propellants in launch vehicle tanks is necessary for mission planning and successful execution. Cryogenic propellant storage and transfer in space environments requires that tank pressure be controlled. The pressure rise rate is determined by the complex interaction of external heat leak, fluid temperature stratification, and interfacial heat and mass transfer. If the required storage duration of a space mission is longer than the period in which the tank pressure reaches its allowable maximum, an appropriate pressure control method must be applied. Therefore, predictions of the pressurization rate and performance of pressure control techniques in cryogenic tanks are required for development of cryogenic fluid long-duration storage technology and planning of future space exploration missions. This paper describes an analytical tool, Tank System Integrated Model (TankSIM), which can be used for modeling pressure control and predicting the behavior of cryogenic propellant for long-term storage for future space missions. It is written in the FORTRAN 90 language and can be compiled with any Visual FORTRAN compiler. A thermodynamic vent system (TVS) is used to achieve tank pressure control. Utilizing TankSIM, the following processes can be modeled: tank self-pressurization, boiloff, ullage venting, and mixing. Details of the TankSIM program and comparisons of its predictions with test data for liquid hydrogen and liquid methane will be presented in the final paper.

Bolshinskiy, L. G.↗

Finding predictive models for singlet fission by machine learning

Singlet fission (SF), the conversion of one singlet exciton into two triplet excitons, could significantly enhance solar cell efficiency. Molecular crystals that undergo SF are scarce. Computational exploration may accelerate the discovery of SF materials. However, many-body perturbation theory (MBPT) calculations of the excitonic properties of molecular crystals are impractical for large-scale materials screening. We use the sure-independence-screening-and-sparsifying-operator (SISSO) machine-learning algorithm to generate computationally efficient models that can predict the MBPT thermodynamic driving force for SF for a dataset of 101 polycyclic aromatic hydrocarbons (PAH101). SISSO generates models by iteratively combining physical primary features. The best models are selected by linear regression with cross-validation. The SISSO models successfully predict the SF driving force with errors below 0.2 eV. Based on the cost, accuracy, and classification performance of SISSO models, we propose a hierarchical materials screening workflow. Three potential SF candidates are found in the PAH101 set.

36 MATERIALS SCIENCE↗

PNNL-Predictive-Phenomics/ThermoOpt

ThermoOpt: Thermodynamic optimization of metabolic models in dynamic environments through pathway confinement.

Johnson, Connah [Pacific Northwest National Labora↗

Predictions of nucleation theory applied to Ehrenfest thermodynamic transitions

A modified nucleation theory is used to determine a critical nucleus size and a critical activation-energy barrier for second-order Ehrenfest thermodynamic transitions as functions of the degree of undercooling, the interfacial energy, the heat-capacity difference, the specific volume of the transformed phase, and the equilibrium transition temperature. The customary approximations of nucleation theory are avoided by expanding the Gibbs free energy in a Maclaurin series and applying analytical thermodynamic expressions to evaluate the expansion coefficients. Nonlinear correction terms for first-order-transition calculations are derived, and numerical results are presented graphically for water and polystyrene as examples of first-order and quasi-second-order transitions, respectively.

Barker, R. E., Jr.↗

Real gas properties and Space Shuttle Main Engine fuel turbine performance prediction

The H2/H2O mixture thermodynamic and transport properties variations for the Space Shuttle Main Engine (SSME) fuel turbine over a range of temperatures and pressures are examined. The variation of molecular viscosity, specific heat at constant pressure, and Prandtl number for the hydrogen/steam mixture are fitted using polynominal relationships for future turbine performance use. The mixture property variations are calculated using GASP and WASP computer programs. The air equivalent performance of the SSME fuel turbine is computed.

Harloff, G. J.↗

Solid State Reduction Driven Synthesis of Mn Containing Multi-principal Component Alloys

In recent years, high entropy alloys (HEAs), also known as multi-principal component alloys (MPCAs) have emerged as a new and exciting class of materials. This paper reports on the solid state reduction synthesis of a series of CoFeNiMn-based MPCA compositions, starting from a mixture of the corresponding oxides. One of the aims of the study was to test whether the degree of reduction of MnO, a highly stable oxide, could be enhanced by tailoring the alloy composition. Specifically, the influence of Ni content was studied because Ni exhibits a significant negative enthalpy of mixing with Mn. High purity precursor powders of Co(OH) 2 , Fe 2 O 3 , MnO 2 , and NiO were milled and mixed using standard ceramic processing methods. Here, the nominal sample compositions (assuming complete oxide reduction) were (CoFeMn) x Ni (1-x) , for x = 0, 0.083, 0.166, and 0.25. The oxide samples were subjected to a series of isothermal reduction anneals in flowing 3 pct H 2 –Ar at 1100 °C. The resulting microstructures were characterized using scanning electron microscopy (SEM), X-ray energy dispersive spectroscopy (EDS) and X-ray diffraction (XRD). The composition of the resulting MPCAs was determined quantitatively using wavelength dispersive spectroscopy (WDS) in the electron microprobe. The study revealed that for each of the initial oxide compositions studied, it was possible to achieve an MPCA with ~ 25 at. pct Mn. These results were found to be consistent with the predictions of a thermodynamic model whereby a negative enthalpy of mixing (ΔH mix ), combined with a contribution from configurational entropy, can offset a positive free energy of reduction (ΔG red ). The incorporation of vibrational entropy into first principles calculations was found to have a significant effect on the predicted crystal structure of the MPCAs.

36 MATERIALS SCIENCE↗

Temperature effects on the universal equation of state of solids

Recently it has been argued based on theoretical calculations and experimental data that there is a universal form for the equation of state of solids. This observation was restricted to the range of temperatures and pressures such that there are no phase transitions. The use of this universal relation to estimate pressure-volume relations (i.e., isotherms) required three input parameters at each fixed temperature. It is shown that for many solids the input data needed to predict high temperature thermodynamical properties can be dramatically reduced. In particular, only four numbers are needed: (1) the zero pressure (P=0) isothermal bulk modulus; (2)it P=0 pressure derivative; (3) the P=0 volume; and (4) the P=0 thermal expansion; all evaluated at a single (reference) temperature. Explicit predictions are made for the high temperature isotherms, the thermal expansion as a function of temperature, and the temperature variation of the isothermal bulk modulus and its pressure derivative. These predictions are tested using experimental data for three representative solids: gold, sodium chloride, and xenon. Good agreement between theory and experiment is found.

Vinet, P.↗

Temperature effects on the universal equation of state of solids

Recently it has been argued based on theoretical calculations and experimental data that there is a universal form for the equation of state of solids. This observation was restricted to the range of temperatures and pressures such that there are no phase transitions. The use of this universal relation to estimate pressure-volume relations (i.e., isotherms) required three input parameters at each fixed temperature. It is shown that for many solids the input data needed to predict high temperature thermodynamical properties can be dramatically reduced. In particular, only four numbers are needed: (1) the zero pressure (P = 0) isothermal bulk modulus; (2) its P = 0 pressure derivative; (3) the P = 0 volume; and (4) the P = 0 thermal expansion; all evaluated at a single (reference) temperature. Explicit predictions are made for the high temperature isotherms, the thermal expansion as a function of temperature, and the temperature variation of the isothermal bulk modulus and its pressure derivative. These predictions are tested using experimental data for three representative solids: gold, sodium chloride, and xenon. Good agreement between theory and experiment is found.

Vinet, Pascal↗