Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Computational 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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
NASA NTRS1949

Investigation of Conditions of Titanium Carbonization - IV

In a previous paper, results are presented of accurate investigations of the processes of titanium carbonization and the succeeding titanium carbide decarbonization as related to the phenomenon of the graphitization of soot by heating at a constant temperature in atmospheres of pure hydrogen and carbon monoxide. These tests showed that the processes of titanium carbonization-decarbonization in an atmosphere of pure gases without nitrogen proceed in the same direction as the analogous processes under the conditions of the production furnace. In this case, however, the presence of admixtures of nitrogen changes the quantitative results of the decarbonization process. Thermodynamic computations confirming the results of previous tests conducted at atmospheric pressure and additional tests of titanium carbonization at lowered pressures are presented herein.

Meerson, G. A.↗
NASA NTRS1951

General method and thermodynamic tables for computation of equilibrium composition and temperature of chemical reactions

A rapidly convergent successive approximation process is described that simultaneously determines both composition and temperature resulting from a chemical reaction. This method is suitable for use with any set of reactants over the complete range of mixture ratios as long as the products of reaction are ideal gases. An approximate treatment of limited amounts of liquids and solids is also included. This method is particularly suited to problems having a large number of products of reaction and to problems that require determination of such properties as specific heat or velocity of sound of a dissociating mixture. The method presented is applicable to a wide variety of problems that include (1) combustion at constant pressure or volume; and (2) isentropic expansion to an assigned pressure, temperature, or Mach number. Tables of thermodynamic functions needed with this method are included for 42 substances for convenience in numerical computations.

Huff, Vearl N↗
DOE OSTI2024

Chemomechanics in alloy phase stability

We describe a first-principles statistical mechanics method to calculate the free energies of crystalline alloys that depend on temperature, composition, and strain. The approach relies on an extension of the alloy cluster expansion to include an explicit dependence on homogeneous strain in addition to site occupation variables that track the degree of chemical ordering. The method is applied to the Si-Ge binary alloy and is used to calculate free energies that describe phase stability under arbitrary epitaxial constraints. We find that while the incoherent phase diagram (in which coexisting phases are not affected by coherency constraints) hosts a miscibility gap, coherent phase equilibrium predicts ordering and negative enthalpies of mixing. Instead of chemical instability, the chemomechanical free energy exhibits instabilities along directions that couple the composition of the alloy with a volumetric strain order parameter. Furthermore, this has fundamental implications for phase field models of spinodal decomposition as it indicates the importance of gradient energy coefficients that couple gradients in composition with gradients in strain.

Materials Science↗
DOE OSTI2026

An Ab Initio Molecular Dynamics Study of Key Thermodynamic Input Parameters for Computer Simulation of U-6Nb Solidification

The key to metallic fuel development is the fabrication of uranium metal and alloys into fuel forms. U-Nb alloys are one of the best candidates for a metallic fuel alloy with high-temperature strength sufficient to support the core, acceptable nuclear properties, good fabricability, and compatibility with usable coolant media. Melt processing has been a key component of the metallic fuel cycle, and process models require thermophysical parameters at elevated temperatures, particularly above the melting temperatures, regarding which experimental data are scarce, for accurate simulations and process development. By means of ab initio density-functional theory (DFT) quantum molecular dynamics (QMD), we have calculated the main thermophysical parameters—the density, thermal expansion coefficient, specific heat, thermal conductivity, melting temperature, latent heat of fusion, and viscosity—used in the modeling of the U-6 wt.% Nb alloy casting. The melting temperature of the U-6 wt.% Nb alloy at ambient pressure is obtained by means of QMD simulations using the Z-method. The ambient volume change and latent heat of melting of U-6 wt.% Nb are also derived from QMD simulations in conjunction with analytical fitting for the energy and pressure. The thermal conductivity for the solid U-Nb alloy is calculated from the semi-classical Boltzmann transport equation combined with an estimate of the electron relaxation time obtained from DFT simulations.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Structure and Synthesizability of Iron–Sulfur Metal–Organic Frameworks

Sulfur-based metal–organic frameworks (MOFs) and coordination polymers (CPs) are an emerging class of hybrid materials that have received growing attention due to their magnetic, conductive, and catalytic properties with potential applications in electrocatalysis and energy storage. In this work, we report a high-throughput virtual screening protocol to predict the synthesizability of candidate metal–sulfur MOFs/CPs by computing the thermodynamically stable structures resulting from a particular combination of metal cluster, linker, cation, and synthetic conditions. Free energies are computed by using all-atom classical mechanical thermodynamic integration. Low-free-energy structures are refined using ab initio density functional theory, and pair distribution functions and powder X-ray diffraction patterns are calculated to complement and guide experimental structure determination. We validate the computational approach by retrospective predictions of the stable structure produced by experimental syntheses, and a subsequent screen predicts Fe 4 S 4 -BDT–TPP as a new thermodynamically stable one-dimensional (1D) CP comprising a redox-active Fe 4 S 4 cluster, a 1,4-benzenedithiolate (BDT) linker, and a tetraphenylphosphonium (TPP) countercation. Furthermore, this material is experimentally synthesized, and the 1D chain structure of the crystal is confirmed using microcrystal electron diffraction. The computational screening pipeline is generically transferable to neutral and ionic MOFs/CPs comprising arbitrary metal clusters, linkers, cations, and synthetic conditions, and we make it freely available as an open source tool to guide and accelerate the discovery and engineering of novel porous materials.

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

Homologous Ladder Cyclohexasilanes

Here, we report the synthesis of five new examples of ladder cyclohexasilanes, possessing up to three consecutive fused rings and differing in relative ring fusion configuration and side chain structure. By coupling a 1,4-dipotassiooligosilyl dianion to a cyclohexasilane, we obtained bi- and tricyclic ladder cyclohexasilanes. Combined experimental and theoretical studies suggested that annulation could favor the cis configuration under kinetic control, while the trans configuration predominates under thermodynamic control. Computational studies show that with each additional ring in the trans-diastereomeric series, the predicted onset of light absorption shifts to longer wavelengths.

36 MATERIALS SCIENCE↗
DOE OSTI2026

A modified scheil approach for nucleation-dependent solidification pathways

As-solidified microstructures of near-eutectic alloys often contain multiple primary phases that are not expected from equilibrium phase diagrams. Such microstructures are caused by cooling-rate-dependent solidification pathways, a factor not captured by the Scheil–Gulliver model or variations thereof. Here, we present a model and algorithm that incorporate the critical nucleation undercooling for each solid phase into the Scheil–Gulliver model. We hypothesize that the non-equilibrium microstructure formation is primarily governed by a nucleation-competition mechanism. This mechanism accounts for both stable/metastable phase selection and primary-phase formation within eutectic regions driven by asymmetric nucleation barriers. The model is validated against a hypereutectic Al-Fe alloy, where it successfully reproduces the observed microstructural constituents, revealing the key dependencies of solidification microstructure on nucleation kinetics. Applicability to multicomponent systems is demonstrated through a hypereutectic Al–Fe–Si ternary alloy, where the model successfully predicts divorced eutectic microstructures and the associated oscillatory solidification pathways along univariant lines. As a result, the proposed framework establishes a nucleation-dependent computational approach for interpreting and predicting solidification microstructures.

Alloy design↗
DOE OSTI2025

Corrosion of 316 L stainless steel under the natural circulation of molten NaCl-MgCl 2 salt

The corrosion behavior of 316 L stainless steel (SS) was studied via the natural circulation of molten eutectic NaCl-MgCl 2 salt through a microloop. The post-corrosion tested 316 L SS microloop sections were characterized with microscopy techniques to determine the microstructural and microchemical changes that occurred at the alloy/salt interface. It was found that 316 L SS showed heterogeneous dissolution at the hot-leg, whereas deposition of corrosion products occurred at the cold-leg. For the first time, experimentally obtained molten salt flow-induced corrosion of 316 L SS results were combined with computational thermodynamic-kinetic models to validate the dissolution and deposition in terms of elemental compositional changes at the alloy/salt interface. The thermodynamic-kinetic modeling predicted that the heterogeneous dissolution of Cr and Fe from the hot-leg section of 316 L SS persisted throughout the salt circulation. The model also estimated that, despite Cr deposition starting earlier than Fe, the total redeposition of Fe is expected to be significantly greater than that of Cr over the circulation of salt. Furthermore, the modeling accurately predicted the subsurface enrichment of Mo which is attributed to the reduced Cr activity and the relatively higher diffusion rate of Mo within the alloy matrix. Here, the agreement between modeling and experimental results confirms that Fe chlorides dissolve at the hot-leg and subsequently deposit at the cold-leg due to activity changes driven by the thermal gradient. By contrast, Cr was not detected in the cold-leg deposits, which is attributed to its weaker temperature dependence on activity, limiting its redeposition under these conditions.

36 - MATERIALS SCIENCE↗
DOE OSTI2025

In-situ formation of stable interface towards Li-in anode for halide solid-state electrolyte

Halide-based solid-state electrolytes (SSEs) are promising candidates for next-generation all-solid-state lithium batteries (ASSLBs) due to their high ionic conductivity and chemical stability. However, their poor interfacial compatibility with lithium metal anode and Li-In alloy significantly hinder practical application due to the requirement for a protective interlayer. In this study, a novel approach to overcome this limitation is presented by introducing iron (Fe) doping into Li 3 InCl 6 (LIC), which enables direct and stable contact with lithium-indium (Li-In) metal without a protective interlayer. Thermodynamic and computational analyses identified Fe 3+ as a suitable dopant based on its similar reduction potential to In 3+ and structural compatibility within the halide lattice. The synthesized 10 at. % Fe-doped LIC exhibits high phase purity, retained ionic conductivity, and notably improved interfacial stability. Full-cell tests using Fe-LIC achieve over 300 cycles with 80 % capacity retention. At the same time, symmetric Li-In/ Fe-LIC/ Li-In cells sustain over 500 h of operation, representing the first reported long-term cycling of LIC-based ASSLB without a protective interlayer. In conclusion, this work establishes Fe doping as an effective strategy to stabilize halide SSEs of In system against Li-In alloy, thereby simplifying cell architecture and advancing the development of safer, high-performance halide-based solid-state electrolytes.

Halide-based solid-state electrolytes↗
DOE OSTI2025

A Machine Learning Framework for Modeling Ensemble Properties of Atomically Disordered Materials

Atomic disorder can strongly influence material properties such as charge transport, optical response, and catalytic activity. However, efficiently modeling these disorder effects remains challenging for first-principles methods due to the cost of sampling large configurational spaces and computing complex physical quantities. Recent advances of machine learning techniques, particularly graph neural networks (GNNs), has enabled the efficient and accurate predictions of complex material properties, offering promising tools for studying disordered systems. In this work, we present a general machine-learning-assisted computational framework that integrates equivariant GNNs with Monte Carlo simulations to compute the thermodynamic and ensemble-averaged functional properties of disordered materials. Using the surface-termination-disordered MXene monolayer Ti 3 C 2 T 2–x as a representative system, we find that electrical conductivity exhibits an emergent peak near the order–disorder phase transition temperature due to the interplay between electron scattering and doping. In contrast, optical conductivity remains largely insensitive to local atomic disorder and reflects the global surface chemical composition. These results highlight the role of atomic disorder in affecting material properties and demonstrate the potential of our approach for statistically modeling disorder effects in a wide range of materials such as high-entropy alloys and spin liquids.

MXene↗
DOE OSTI2026

Concerted Proton–Electron Transfer Minimizes Substituent Effects on Adsorbed Phthalocyanine Electrocatalysis

Molecularly modified electrodes (MMEs) are potent electrocatalysts, but few principles exist for their rational design. Electrocatalysis by soluble molecules depends strongly on substituents that tune the catalyst redox potential (E 1/2 ), but it is unclear if this parameter similarly impacts MME catalysis. Herein, we employ the hydrogen evolution reaction (HER) as a test case for comparing carbon-adsorbed cobalt phthalocyanine (CoPc/C) and cobalt hexadecafluoro-phthalocyanine (CoFPc/C). By correlating HER activity and voltammetric data to total Co surface concentration across a wide range of catalyst loadings, we find that only 5–25% of adsorbed Co sites contribute to the Co(II/I) redox wave and that this subpopulation poorly correlates with catalytic activity. Instead, in the low-loading limit, catalytic activity correlates linearly with the majority Co(II/I)-silent Co population, revealing per-site turnover frequency (TOF) values for HER. Despite a 230 mV difference in Co(II/I) redox potentials, CoPc/C and CoFPc/C display TOF values differing by less than a factor of 3 when compared over a wide potential range. Mechanistic studies point to an inner-sphere concerted proton–electron transfer step as rate-determining, suggesting that the Co–H bond dissociation free energy (BDFE) rather than the Co(II/I) E 1/2 is thermodynamically relevant. Computational studies indicate that the fluoro-substituents lead to compensatory changes in Co(II/I) E 1/2 and Co(I) basicity, leaving the Co–H BDFE largely unchanged between CoPc and CoFPc and thereby manifesting in similar catalytic rates. Furthermore, these results highlight the limited effect of E 1/2 -tuning on MME catalytic activity and motivate the development of methods to directly alter active site–substrate BDFE.

Catalysts↗
DOE OSTI2026

Disorder-induced magnetoelastic behaviors of MnTexSbyBi1-x-y alloys

This dataset contains input and output files from density functional theory (DFT) simulations used to study the disorder-induced magnetoelastic behaviors of MnTexSbyBi1-x-y (0 ≤ x + y ≤ 1) alloys and their binary end members MnTe, MnSb, and MnBi. The alloys adopt the hexagonal NiAs-type (nickeline) structure and span ternary (MnTexSb1-x, MnTexBi1-x, MnBixSb1-x), and quaternary compositions across the full MnTe–MnSb–MnBi composition triangle. For each alloy composition, the dataset provides DFT calculations in three magnetic configurations: A-type antiferromagnetic (AFM), C-type AFM, and ferromagnetic (FM). Every magnetic configuration folder contains the fully relaxed crystal structure (CONTCAR), VASP input parameters (INCAR), and the main VASP output file (OUTCAR), from which total electronic energies, Mn magnetic moments, lattice parameters, and percent volume changes between magnetic states are extracted. These data are used to construct compositional phase diagrams, evaluate thermodynamic stability (formability), and map magnetoelastic responses across the alloy space. For A-type AFM and FM configurations, additional data are provided as follows: (i) FORCE_CONSTANTS and thermal_properties.yaml files at the top level of A-type_AFM/ and FM/ folders — present only for compositions marked with an asterisk (*) in Table I of the main text. These are derived from Phonopy finite-displacement calculations on full disordered 128-atom supercells and provide vibrational free energy, entropy (Svib)contribution from explicit disorder calculations. (Table I of the associated main manuscript) (ii) A VCA/ subfolder within A-type_AFM/ and FM/, containing FORCE_CONSTANTS and thermal_properties.yaml from Virtual Crystal Approximation phonon calculations (without spin-orbit coupling). VCA data are available for all compositions and are used to estimate vibrational contributions to the Gibbs free energy across the full composition space. (iii) A SOC/ subfolder containing CONTCAR, INCAR, and OUTCAR from spin-orbit coupling calculations, providing relativistic corrections to electronic energies and lattice parameters (Tables S2–S3 of the SM, and Table I of the main manuscript). (iv) A SOC/VCA/ subfolder containing FORCE_CONSTANTS and thermal_properties.yaml from VCA phonon calculations performed within the SOC framework, combining relativistic and vibrational thermodynamic corrections. The computed properties are used to map the AFM–FM magnetic crossover near MnTe0.75Sb0.25, demonstrate disorder- and spin-induced phonon broadening, identify a semiconductor-to-metal crossover, and quantify the pronounced magnetoelastic volume response near the magnetic phase boundary.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Dual dressed black holes as the end point of the charged superradiant instability in $\mathcal{N} = 4$ Yang Mills

Charged Black holes in AdS_5 × S^5 A d S 5 × S 5 suffer from superradiant instabilities over a range of energies. Hairy black hole solutions (constructed within gauged supergravity) have previously been proposed as endpoints to this instability. We demonstrate that these hairy black holes are themselves unstable to the emission of large dual giant gravitons. We propose that the endpoint to this instability is given by Dual Dressed Black Holes (DDBH)s; configurations consisting of one, two, or three very large dual giant gravitons surrounding a core AdS A d S black hole with one, two, or three SO(6) S O ( 6 ) chemical potentials equal to unity. The dual giants each live at AdS A d S radial coordinates of order \sqrt{N} N and each carry charge of order N^2 N 2 . The large separation makes DDBHs a very weakly interacting mix of their components and allows for a simple computation of their thermodynamics. We conjecture that DDBHs dominate the phase diagram of \mathcal{N}=4 𝒩 = 4 Yang-Mills over a range of energies around the BPS plane, and provide an explicit construction of this phase diagram, briefly discussing the interplay with supersymmetry. We develop the quantum description of dual giants around black hole backgrounds and explicitly verify that DDBHs are stable to potential tunneling instabilities, precisely when the chemical potentials of the core black holes equal unity. We also construct the 10-dimensional DDBH bulk solutions.

Choi, Sunjin (ORCID:0000000283310199)↗
DOE OSTI2026

Harnessing Ocean Thermal Gradients Using Thermoelectric Based Submersibles for Ocean Power Applications

The urgent need for energy solutions in marine environments has accelerated the development of innovative technologies capable of leveraging natural resources for power generation. This study introduces a buoyancy-driven submersible system designed to harness ocean thermal gradients using thermoelectric generators (TEGs) and phase change materials (PCMs). The technology aims to provide autonomous power to offshore aquaculture farms, unmanned underwater vehicles (UUVs), offshore platform illumination, and ocean sensors, significantly reducing dependence on fossil fuels. Ocean thermal gradients, especially prevalent in mid-latitude regions, exhibit temperature differences between surface and deep waters ranging from 7 degrees Celsius to 30 degrees Celsius depending on seasonal variations. The proposed submersible technology utilizes TEGs to convert thermal energy from these gradients into electrical power, generating between 0.2 and 0.5 watts, while PCMs are employed to store and regulate this energy, ensuring a stable and continuous power supply. The buoyancy-driven mechanism of the submersible enhances its capability to navigate through varying depths, optimizing its exposure to different thermal gradients and maximizing energy harvesting. The performance of this submersible system is analyzed through detailed thermodynamic assessments and computational fluid dynamics (CFD) modeling focused on heat transfer. These analyzes consider real-world ocean temperature profiles and seek to refine the interaction between TEGs and PCMs to optimize energy extraction. The evaluation encompasses several key performance metrics, including power output and energy efficiency. Results confirm the potential of this innovative technology to provide a continuous and reliable power source for marine applications. By demonstrating the feasibility of using ocean thermal gradients for energy generation, this study contributes to the broader efforts of innovation in energy technologies for harsh, remote marine environments. The implementation of such promises is significant advancements in the autonomy of marine operations. The ongoing research will further investigate scalability ensuring that the system can be effectively adapted to various marine settings and operational demands.

16 TIDAL AND WAVE POWER↗