Search NASA⌕ Search

SEARCH · Search NASA

Results for “thermodynamics of computing”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

A mathematical framework for thermodynamic computing with applications to chemical reaction networks

The widespread adoption of energy-intensive computing applications has led to a growing need for energy-efficient computing approaches. Thermodynamic computing offers a promising approach for low-energy computation by leveraging the intrinsic computational capabilities of physical, chemical, or biological systems. However, the mathematical foundations of thermodynamic computing require further development to fully realize the potential energy efficiencies, as well as to assess factors like noise and operational speed. In this paper, we establish a mathematical framework for utilizing thermodynamic processes to perform fundamental operations, including addition, subtraction, multiplication, and division. We highlight the use of chemical reactions as potential computational units and explore synthetic chemical and biochemical systems as practical implementations. Additionally, we demonstrate how these principles can be applied to solving complex mathematical problems, such as ordinary differential equations (ODEs) and suggest the necessary components to implement the thermodynamic computing framework using chemical reactions based in a microfluidic device. This work enhances our understanding of thermodynamic processes for natural computing as a basis for scalable, energy-efficient computation in paradigm disruptive next-generation systems.

Cannon, William R. [Pacific Northwest National Lab↗
DOE OSTI2026

Nonlinear thermodynamic computing out of equilibrium

We present the design for a thermodynamic computer that can perform arbitrary nonlinear calculations in or out of equilibrium. Simple thermodynamic circuits, fluctuating degrees of freedom in contact with a thermal bath and confined by a quartic potential, display an activity that is a nonlinear function of their input. Such circuits can therefore be regarded as thermodynamic neurons, and can serve as the building blocks of networked structures that act as thermodynamic neural networks, universal function approximators whose operation is powered by thermal fluctuations. We simulate a digital model of a thermodynamic neural network, and show that its parameters can be adjusted by genetic algorithm to perform nonlinear calculations at specified observation times, regardless of whether the system has attained thermal equilibrium. This work expands the field of thermodynamic computing beyond the regime of thermal equilibrium, enabling fully nonlinear computations, analogous to those performed by classical neural networks, at specified observation times.

Whitelam, Stephen [Lawrence Berkeley National Labo↗
DOE OSTI2026

Generative Thermodynamic Computing

Here, we introduce a generative modeling framework for thermodynamic computing, in which structured data are synthesized from noise by the natural time evolution of a physical system governed by Langevin dynamics. While conventional diffusion models use neural networks to perform denoising, here the information needed to generate structure from noise is encoded by the dynamics of a thermodynamic system. Training proceeds by maximizing the probability with which the computer generates the reverse of a noising trajectory, which ensures that the computer generates data with minimal heat emission. We demonstrate this framework within a digital simulation of a thermodynamic computer. If realized in analog hardware, such a system would function as a generative model that produces structured samples without the need for artificially injected noise or active control of denoising.

Whitelam, Stephen [Lawrence Berkeley National Labo↗
DOE OSTI2021

Computational thermodynamic study of SiC chemical vapor deposition from MTS-H 2

This study focuses on the computational thermodynamic analysis of the chemical vapor deposition (CVD) of SiC from the methyltrichlorosilane-hydrogen (MTS-H 2 ) using up-to-date thermodynamic databases. High-resolution computation has been performed with the fine intervals of temperature and pressure at the various H 2 /MTS ratios of interest to systematically investigate the deposition condition range (800 to 1600°C, 0 to 26 664 Pa, and H 2 /MTS ratios of 0.1 to 100) to guide experimental exploration. The influence of deposition parameters on the compositions and phase stabilities of the deposit and gas phase pertinent to vapor processing is elucidated. Low pressure and medium temperatures (1000 to 1400°C) are beneficial to reaching a higher SiC deposition efficiency and provide an optimal window for preparing a high-purity (>99 wt.% SiC) deposit. This optimal processing window expands significantly with an increasing H 2 /MTS ratio (<20). These results are supported by a number of previous theoretical and experimental observations. The mass fraction of SiC in deposit is proposed as an additional perspective to understand the discrepancy between thermodynamic calculation and experimental observation of pure CVD SiC at low H 2 /MTS ratios.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Leveraging computational thermodynamics to guide SiC-ZrC chemical vapor deposition process development

Here using the CALPHAD approach to understand zirconium carbide deposition, a series of phase equilibria were calculated from a custom thermodynamic database based on a literature source, and the equilibria were used to explore the potential chemical vapor deposition (CVD) processing space in the ZrCl 4 -CH 3 SiCl 3 -CH 4 -H 2 system as a function of pressure, temperature, and gas composition. Several gas ratios were considered. At a given ZrCl 4 :CH 3 SiCl 3 ratio within the range studied, the most important factor was found to be the ratios of CH 4 :ZrCl 4 , wherein the nature of the composition – carbide vs. silicide – could be controlled. A pure binary composition of ZrC and SiC is expected to form by increasing the initial amount of methane and decreasing the amount of hydrogen from values predicted purely based on thermodynamic equilibrium. Rietveld analysis of the x-ray diffractograms from corresponding experimental depositions confirmed that increasing the CH 4 :ZrCl 4 ratio increased the fraction of carbon-containing species (SiC, ZrC) and decreased the fraction of non-carbides (ZrSi, ZrSi 2 , etc.), as predicted from the CALPHAD results.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Computing Thermodynamic Properties of Fluids Augmented by Nanoconfinement: Application to Pressurized Methane

Nanoconfined fluids exhibit remarkably different thermodynamic behavior compared to the bulk phase. These confinement effects render predictions of thermodynamic quantities of nanoconfined fluids challenging. In particular, confinement creates a spatially varying density profile near the wall that is primarily responsible for adsorption and capillary condensation behavior. Significant fluctuations in thermodynamic quantities, inherent in such nanoscale systems, coupled to strong fluid–wall interactions give rise to this near-wall density profile. Empirical models have been proposed to explain and model these effects, yet no first-principles based formulation has been developed. We present a statistical mechanics framework that embeds such a coupling to describe the effect of the fluid–wall interaction in amplifying the near-wall density behavior for compressible gases at elevated pressures such as pressurized methane in confinement. We show that the proposed theory predicts accurately the adsorbed layer thickness as obtained with small-angle neutron scattering measurements. Furthermore, the predictions of density under confinement from the proposed theory are shown to be in excellent agreement with available experimental and atomistic simulations data for a range of temperatures for nanoconfined methane. While the framework is presented for evaluating the near-wall density, owing to its rigorous foundation in statistical mechanics, the proposed theory can also be generalized for predicting phase-transition and nonequilibrium transport of nanoconfined fluids.

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

Design of crack-free laser additive manufactured Inconel 939 alloy driven by computational thermodynamics method

Abstract This paper examined the effect of Si addition on the cracking resistance of Inconel 939 alloy after laser additive manufacturing (AM) process. With the help of CALculation of PHAse Diagrams (CALPHAD) software Thermo-Calc, the amounts of specific elements (C, B, and Zr) in liquid phase during solidification, cracking susceptibility coefficients (CSC) and cracking criterion based on $$\left| {{\text{d}}T/{\text{d}}f_{{\text{s}}}^{1/2} } \right|$$ d T / d f s 1 / 2 values ( T : solidification temperature, f s : mass fraction of solid during solidification) were evaluated as the indicators for composition optimization. It was found that CSC together with $$\left| {{\text{d}}T/{\text{d}}f_{{\text{s}}}^{1/2} } \right|$$ d T / d f s 1 / 2 values provided a better prediction for cracking resistance. Graphical abstract

Zeng, Congyuan (ORCID:0000000214764277)↗
DOE OSTI2024

Computational thermodynamics-guided alloy design and phase stability in CoCrFeMnNi-based medium-and high-entropy alloys: An experimental-theoretical study

A computational thermodynamics approach has been employed to design CoCrFeMnNi-based medium- and high-entropy alloys (M/HEAs) with systematically varied compositions (Co(( 80-X )/ 2 )Cr(( 80-X )/ 2 )Fe X Mn 10 Ni 10 with x = 30, 40, and 50 at.%) and phase stability. Since the formation of sigma phase, usually brittle and undesirable, is a common concern, when this class of alloys is subjected to elevated temperatures (600–1000 °C), predicting its formation becomes essential. Thus, its formation and the phase equilibria were studied using the CALPHAD method, and two empirical methods, namely, valence electron concentration (VEC) and paired sigma-forming element (PSFE). Isothermal aging treatments at 900–1100 °C for 20 h were performed, since CALPHAD and VEC/PSFE predictions diverged. Both prediction methods were compared with experimental characterization by a combination of scanning electron microscopy and high-energy synchrotron X-ray diffraction. In conclusion, the predictions from the VEC/PSFE and CALPHAD calculations (depending on the database used) were shown to be quite accurate.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Chemistry Informed Machine Learning-Based Heat Capacity Prediction of Solid Mixed Oxides

Knowing heat capacity is crucial for modeling temperature changes with the absorption and release of heat and for calculating the thermal energy storage capacity of oxide mixtures with energy applications. The current prediction methods (ab initio simulations, computational thermodynamics, and the Neumann–Kopp rule) are computationally expensive, not fully generalizable, or inaccurate. Machine learning has the potential of being fast, accurate, and generalizable, but it has been scarcely used to predict mixture properties, particularly for mixed oxides. Here, we demonstrate a method for the generalizable prediction of heat capacity of solid oxide pseudobinary mixtures using heat capacity data obtained from computational thermodynamics and descriptors from ab initio databases. Further, models trained through this workflow achieved an error (mean absolute error of 0.43 J mol –1 K –1 ) lower than the uncertainty in differential scanning calorimetry measurements, and the workflow can be extended to predict other properties derived from the Gibbs free energy and for higher-order oxide mixtures.

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

The Fundamental Creep Behavior Model of Gr.91 Alloy by ICME Approach (Final Report)

The current project mainly focused on the fundamental creep cracking mechanism of the Grade 91 system during the operation conditions in the advanced technologies FE power plants and build the link among Composition-Processing Parameters-Phase Stability-Microstructure-Creep Resistance. A model based on computational thermodynamics and diffusion kinetics will be developed to provide guidance on how to improve the creep resistance of the alloy system. The PI mainly investigate the Gr.91 base alloy and weldment with the Integrated Computational Materials Engineering (ICME) approach. The long-term goal of the proposed program is to develop a model for creep and fatigue resistant alloys with different elemental systems, compositions, and processing parameters of weldment and heat-treatment. In addition, new creep and fatigue resistant alloys are to be designed based on the predictions from the model. The computational as well as experimental results include three sections: 1) The compositional optimization and secondary phases evaluation regarding the creep resistance in Grade 91 steel through the CALPHAD approach. In this section, the formation of the critical secondary phases, M23C6, MX and Z-phase were well predicted based on many conditions like the temperature as well as compositions. Meanwhile, the role of different alloying elements was also considered as they will directly affect the formation of the secondary phases mentioned above based on the compositional evolution of the alloying elements. 2) An investigation of Creep Resistance in Grade 91 Steel through Computational Thermodynamics. We systematically considered the evolution of the four critical temperatures, Ac3, Ac1, the threshold of M23C6 and Z-phase, on basis of the computational approach. Meanwhile, the equilibrium cooling as well as Scheil simulations were also considered to perform further predictions based on the different cooling rates. 3) Creep lifetime and microstructural evolution of the Grade 91 steel. In this part, mainly experimental approach was performed on the real creep test from the beginning of the commercial Gr.91 alloys. The alloys were normalized, tempered, welded (Gleeble) and PWHT and finally subjected to the creep lifetime analyses. Finally, all the samples after each process were further characterized based on the OM, XRD and SEM technique to draw final conclusions.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2024

One-step sputtering of MoSSe metastable phase as thin film and predicted thermodynamic stability by computational methods

Abstract We present the fabrication of a MoS 2−x Se x thin film from a co-sputtering process using MoS 2 and MoSe 2 commercial targets with 99.9% purity. The sputtering of the MoS 2 and MoSe 2 was carried out using a straight and low-cost magnetron radio frequency sputtering recipe to achieve a MoS 2−x Se x phase with x = 1 and sharp interface formation as confirmed by Raman spectroscopy, time-of-flight secondary ion mass spectroscopy, and cross-sectional scanning electron microscopy. The sulfur and selenium atoms prefer to distribute randomly at the octahedral geometry of molybdenum inside the MoS 2−x Se x thin film, indicated by a blue shift in the A 1g and E 1 g vibrational modes at 355 cm −1 and 255 cm −1 , respectively. This work is complemented by computing the thermodynamic stability of a MoS 2−x Se x phase whereby density functional theory up to a maximum selenium concentration of 33.33 at.% in both a Janus-like and random distribution. Although the Janus-like and the random structures are in the same metastable state, the Janus-like structure is hindered by an energy barrier below selenium concentrations of 8 at.%. This research highlights the potential of transition metal dichalcogenides in mixed phases and the need for further exploration employing low-energy, large-scale methods to improve the materials’ fabrication and target latent applications of such structures.

36 MATERIALS SCIENCE↗
DOE OSTI2026

ZENN: A thermodynamics-inspired computational framework for heterogeneous data–driven modeling

Traditional entropy-based methods—such as cross-entropy loss in classification problems—have long been essential tools for representing the information uncertainty and physical disorder in data and for developing artificial intelligence algorithms. However, the rapid growth of data across various domains has introduced new challenges, particularly the integration of heterogeneous datasets with intrinsic disparities. To address this, we introduce a zentropy-enhanced neural network (ZENN), extending zentropy theory into the data science domain via intrinsic entropy, enabling more effective learning from heterogeneous data sources. ZENN simultaneously learns both energy and intrinsic entropy components, capturing the underlying structure of multisource data. To support this, we redesign the neural network architecture to better reflect the intrinsic properties and variability inherent in diverse datasets. We demonstrate the effectiveness of ZENN on classification tasks and energy landscape reconstructions, showing its superior generalization capabilities and robustness-particularly in predicting high-order derivatives. In image and text classification tasks, ZENN demonstrates superior generalization by introducing a learnable temperature variable that models latent multisource heterogeneity, allowing it to surpass state-of-the-art models on CIFAR-10/100, BBC News, and AG News. As a practical application in materials science, we employ ZENN to reconstruct the Helmholtz energy landscape of Fe3Pt using data generated from density functional theory and capture key material behaviors, including negative thermal expansion and the critical point in the temperature–pressure space. Overall, this work presents a zentropy-grounded framework for data-driven machine learning, positioning ZENN as a versatile and robust approach for scientific problems involving complex, heterogeneous datasets.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Understanding the microstructural stability in a y’-strengthened Ni-Fe-Cr-Al-Ti alloy

Ni-Fe-Cr-Al-Ti alloys, with Ni levels near 50 wt.%, have the potential to develop a microstructure consisting of a face-centered cubic ? matrix with the homogeneous precipitation of fine ordered ?’ precipitates similar to traditional Ni-based superalloys with significantly greater Ni content. Scanning electron microscopy (SEM), transmission electron microscopy (TEM), atom probe tomography (APT), and CALPHAD -based thermodynamic modeling were employed to understand the phase stabilities and microstructural evolution in an age-hardenable Ni-27Fe-18Cr-1Co-1.6Al-3.75Ti-1.2Mo-0.03C (wt%) alloy. The primary heat-treatment of solution annealing at 1121°C for 4h and age-hardening treatment at 760 °C for 16h resulted in a microstructure consisting of fine ?’ precipitates in an austenitic matrix along with grain boundary carbides, consistent with thermodynamic calculations. Long-term aging at 900 °C for 250h resulted in the coarsening of ?’ along with a change in the morphology of the precipitates from spherical to a more cuboidal shape. In addition, ? phase formation was observed concomitant with the partial dissolution of the ?’ phase. The ability of computational thermodynamic models to predict microstructural characteristics is discussed.

Gwalani, Bharat↗
DOE OSTI2024

Physics-coupled data-driven design of high-temperature alloys

We present a materials design loop, which streamlines physics-coupled machine learning (ML) surrogate models to discover new alloy chemistries with improved properties. The efficacy is demonstrated by discovering a high-temperature alumina-forming austenitic (AFA) stainless steel with enhanced creep, followed by experimental validation. The ML models have been trained using a well-curated, highly consistent experimental dataset augmented with synthetic microstructural features from a computational thermodynamic approach. We have populated a large number of hypothetical AFA alloys to explore the high-dimensional composition space and have predicted their creep properties by providing the same synthetic input features obtained from the trained ML models. Uncertainties from the ML training were taken as thresholds for truncating predicted results to identify alloys with improved or deteriorated creep. Individual elemental compositions have been determined via probability density distribution analysis from the group of alloys at the top and bottom of the predicted creep values for further virtual and experimental validations. In conclusion, we anticipate that this workflow can be applied to screen desired conditions, such as chemistry and processing parameters, in high-dimensional space through physics-guided data analytics.

Alloy design↗
DOE OSTI2021

A Computational Study of RNA Tetraloop Thermodynamics, Including Misfolded States

An important characteristic of RNA folding is the adoption of alternative configurations of similar stability, often referred to as misfolded configurations. These configurations are considered to compete with correctly folded configurations, although their rigorous thermodynamic and structural characterization remains elusive. Tetraloop motifs found in large ribozymes are ideal systems for an atomistically detailed computational quantification of folding free energy landscapes and the structural characterization of their constituent free energy basins, including nonnative states. In this work, we studied a group of closely related 10-mer tetraloops using a combined parallel tempering and metadynamics technique that allows a reliable sampling of the free energy landscapes, requiring only knowledge that the stem folds into a canonical A-RNA configuration. Here we isolated and analyzed unfolded, folded, and misfolded populations that correspond to different free energy basins. We identified a distinct misfolded state that has a stability very close to that of the correctly folded state. This misfolded state contains a predominant population that shares the same structural features across all tetraloops studied here and lacks the noncanonical A-G base pair in its loop portion. Further analysis performed with biased trajectories showed that although this competitive misfolded state is not an essential intermediate, it is visited in most of the transitions from unfolded to correctly folded states. Moreover, the tetraloops can transition from this misfolded state to the correctly folded state without requiring extensive unfolding.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗