SEARCH · Engineering Papers

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 793 records · Page 44

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Theoretical study of CO2 capture mechanisms of SrO and Sr(OH)2.nH2O (n=0,1,8)

In this presentation, the adsorption reactions of CO2 on SrO surfaces and its hydrated counterparts are examined. The predicted energies of CO2 adsorption on follows the order of increasing H2O content. The carbonation on the monohydrate surface is barrierless, whereas the CO2 reaction on the Sr(OH)2 and the octahydrate surfaces follows pathways with finite activation barriers. The thermodynamics of bulk reaction is also addressed.

computational science and engineering↗

DOE OSTI2021

Defying Thermodynamics: Stabilization of Alane Within Covalent Triazine Frameworks for Reversible Hydrogen Storage

The highly unfavorable thermodynamics of direct aluminum hydrogenation can be overcome by stabilizing alane within a nanoporous bipyridine-functionalized covalent triazine framework (AlH 3 @CTF-bipyridine). This material and the counterpart AlH 3 @CTF-biphenyl rapidly desorb H 2 between 95 and 154°C, with desorption complete at 250°C. Sieverts measurements, 27 Al MAS NMR and 27 Al{ 1 H} REDOR experiments, and computational spectroscopy reveal that AlH 3 @CTF-bipyridine dehydrogenation is reversible at 60°C under 700 bar hydrogen, >10 times lower pressure than that required to hydrogenate bulk aluminum. DFT calculations and EPR measurements support an unconventional mechanism whereby strong AlH 3 binding to bipyridine results in single-electron transfer to form AlH 2 (AlH 3 ) n clusters. Here the resulting size-dependent charge redistribution alters the dehydrogenation/rehydrogenation thermochemistry, suggesting a novel strategy to enable reversibility in high-capacity metal hydrides.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Defying Thermodynamics: Stabilization of Alane Within Covalent Triazine Frameworks for Reversible Hydrogen Storage

Metastable metal hydrides such as AlH 3 have many attractive features as hydrogen storage media, but generally require complex reaction schemes for regeneration following H 2 release. Here in this paper, we demonstrate that the highly unfavorable thermodynamics of direct aluminum hydrogenation can be overcome by stabilizing alane within a nanoporous bipyridine-functionalized Covalent Triazine Framework (AlH 3 @CTF-bipyridine). This material and the counterpart AlH 3 @CTF-biphenyl rapidly desorb H 2 between 95 and 154 °C, with desorption complete at 250 °C. Sieverts measurements, 27 Al MAS NMR and 27 Al{ 1 H} REDOR experiments, and computational spectroscopy reveal that AlH 3 @CTF-bipyridine dehydrogenation is reversible at 60 °C under 700 bar hydrogen, >10 times lower pressure than that required to hydrogenate bulk aluminum. DFT calculations and EPR measurements support an unconventional mechanism whereby strong AlH 3 binding to bipyridine results in single-electron transfer to form AlH 2 (AlH 3 ) n clusters. The resulting size-dependent charge redistribution alters the dehydrogenation/rehydrogenation thermochemistry, suggesting a novel strategy to enable reversibility in high-capacity metal hydrides.

Coordination Chemistry↗

DOE OSTI2026

Gradient-informed Hamiltonian Monte Carlo for multicomponent CALPHAD model optimization and uncertainty quantification

CALPHAD model parameter optimization is inherently challenging due to non-smooth objective functions, high-dimensional parameter spaces, and the need for uncertainty quantification (UQ). Traditional weighted nonlinear least squares approaches are computationally efficient but local, whereas black-box global optimizers and ensemble Markov Chain Monte Carlo (MCMC) methods provide broader exploration at substantial computational cost. The objective of this work is to combine the global exploration capability of gradient-informed Hamiltonian Monte Carlo – specifically the No-U-Turn Sampler (NUTS) – with local deterministic refinement using BFGS to efficiently optimize multicomponent CALPHAD models with minimal manual intervention. Analytic gradients are computed via the Jansson derivative framework. The methodology is demonstrated on the Cr—Fe binary system and extended to the Cr—Fe—Ni ternary system with 32 degrees of freedom. For Cr—Fe, NUTS achieves comparable or superior optimality relative to ensemble MCMC while requiring over an order-of-magnitude fewer likelihood evaluations. Parameter uncertainties are quantified through NUTS sampling and propagated to thermodynamic observables using local expansion, demonstrating a novel modular approach that combines binary and ternary parameter subsets without requiring global relaxation. These results establish gradient-informed exploration as a scalable strategy for multicomponent CALPHAD optimization and provide a practical route towards efficient higher-order database development with quantified uncertainty.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Photochemistry of Hypervalent Iodoazide Derivatives

The photochemical properties and reaction mechanisms of a series of hypervalent iodoazide compounds (R–IN3) were investigated, with substituents (−CH 3 , −H, −CF 3 ) tuning the electronic density of the phenyl ring. With the help of UV irradiation, ultrafast time-resolved spectroscopy, and density functional theory (DFT) calculations, we elucidated the mechanisms of azide radical (N 3 • ) release and its subsequent reactivity. UV–vis spectroscopy reveals that the photoconversion rates follow the trend CH3 > H > CF3, aligning with DFT-calculated ΔG values for ring-opening transitions. Homolytic cleavage of the I–N bond is identified as the dominant pathway for N3• generation upon UV irradiation, occurring within 400 fs. The released azide radicals react with the solvent molecules, while the iodo radical R–I C • fragments undergo a thermodynamically uphill lactone ring-opening (RO) reaction and subsequent hydrogen atom abstraction from the solvent to form carboxylic acids R–I–COOH, as validated by NMR and IR spectroscopy. The study also highlights the role of substituents in influencing reaction kinetics and intermediate stability, with electron-donating groups accelerating the release of the N 3 • species. This work bridges experimental observations with computational predictions, offering a foundation for future advancements in azide-based reactions and materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Neural network representations of multiphase Equations of State

Abstract Equations of State model relations between thermodynamic variables and are ubiquitous in scientific modelling, appearing in modern day applications ranging from Astrophysics to Climate Science. The three desired properties of a general Equation of State model are adherence to the Laws of Thermodynamics, incorporation of phase transitions, and multiscale accuracy. Analytic models that adhere to all three are hard to develop and cumbersome to work with, often resulting in sacrificing one of these elements for the sake of efficiency. In this work, two deep-learning methods are proposed that provably satisfy the first and second conditions on a large-enough region of thermodynamic variable space. The first is based on learning the generating function (thermodynamic potential) while the second is based on structure-preserving, symplectic neural networks, respectively allowing modifications near or on phase transition regions. They can be used either “from scratch” to learn a full Equation of State, or in conjunction with a pre-existing consistent model, functioning as a modification that better adheres to experimental data. We formulate the theory and provide several computational examples to justify both approaches, highlighting their advantages and shortcomings.

Science & Technology - Other Topics↗

DOE OSTI2025

Ba 1−x Sr x FeO 3−δ as an improved oxygen storage material for chemical looping air separation: a computational and experimental study

Chemical looping air separation (CLAS) is a promising technology to generate oxygen-rich gas streams to enable efficient carbon dioxide capture during fossil fuel combustion or gasification. CLAS relies on the capture and release of oxygen from the atmosphere using the redox properties of an oxygen-selective solid oxide carrier. This study investigates the redox characteristics of Ba 1−x Sr x FeO 3−δ (0.0 ≤ x ≤ 0.417, 0.0 ≤ δ ≤ 0.5) using a combination of density functional theory (DFT) calculations and experimental verification using X-ray diffraction, thermogravimetric analysis, and oxygen-temperature-programmed desorption. The DFT computed energies of the Ba 1−x Sr x FeO 3−δ perovskites reveal a composition-dependent transition from hexagonal to cubic phases as the Sr-concentration or oxygen vacancy concentration increases. Oxygen vacancy formation energies of the cubic perovskites are found to be lower than those of their hexagonal counterparts. A low oxygen diffusion barrier of ∼1 eV combined with the thermodynamic preference of Ba 1−x Sr x FeO 3−δ compositions that form in a cubic phase suggests them as promising candidates for oxygen storage applications. The experimental results corroborate this finding by identifying Ba 0.75 Sr 0.25 FeO 3−δ in the cubic phase as an optimal composition offering low-temperature oxygen storage capacities comparable to that of the state-of-the-art Sr 0.75 Ca 0.25 FeO 3−δ perovskite oxygen storage material at 325 °C and 350 °C.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Descriptors for Cu facets for CO2 reduction reaction activity

Computation screening is crucial for designing efficient electrochemical catalysts for carbon dioxide (CO2R) reduction to valuable hydrocarbons and oxygenates. Herein, leveraging density functional theory calculations of the CO adsorption energy ΔE_CO on seventeen Cu terminations, we discover a strong linear correlation between ΔE_CO and the recently experimentally measured CO2R electrochemical currents (ACS Catal. 2022, 12, 11, 6578–6588). Examining the ab initio thermodynamics of the early critical intermediates CO*, COH*, and CHO*, we find that CO* → CHO* is the thermodynamically preferred step, and notably shows a volcano trend with the experimental currents where the maximum CO2R current corresponds to the moderate CHO* formation energy. Importantly, we show that increasing the step and kink density of the Cu termination not only enhances CO adsorption strength but also modulates the CO* → CHO* pathway, as respectively exemplified in the (941) and (741) facets. We also explain why (741) is exceptional with high CO2R activity as measured experimentally due to its relatively low activity toward the hydrogen evolution reaction compared with the other Cu surfaces. Beyond the general CO adsorption energy that only shows a linear trend with CO2R activity, we show that the reaction CO* → CHO* free energy is a descriptor that displays a volcano relationship with the overall CO2R activity on Cu facets.

machine learning↗

DOE OSTI2026

Ab Initio Many Body Quantum Embedding and Local Correlation in Crystalline Materials using Interpolative Separable Density Fitting

We present an efficient implementation of ab initio many-body quantum embedding and local correlation methods for infinite periodic systems through translational symmetry adapted interpolative separable density fitting, an approach which reduces the scaling of the calculations to only linear with the number of k-points. Employing this methodology, we compute correlated ground-state coupled cluster energies within density matrix embedding and local natural orbital correlation frameworks for both weakly and strongly correlated solids, using up to 1000 k-points. By extrapolating the local correlation domains and k-point sampling we further obtain estimates of the full coupled cluster with singles, doubles, and perturbative triples ground-state energies in the thermodynamic limit.

Chemical Physics (physics.chem-ph)↗

DOE OSTI2022

Modeling of Supercritical CO2 Shell-and-Tube Heat Exchangers Under Extreme Conditions: Part II: Heat Exchanger Model

Abstract Heat exchangers play a critical role in supercritical CO2 Brayton cycles by providing necessary waste heat recovery. Supercritical CO2 thermal cycles potentially achieve higher energy density and thermal efficiency operating at elevated temperatures and pressures. Accurate and computationally efficient estimation of heat exchanger performance metrics at these conditions is important for the design and optimization of sCO2 systems and thermal cycles. In this paper (Part II), a computationally efficient and accurate numerical model is developed to predict the performance of shell-and-tube heat exchangers (STHXs). Highly accurate correlations reported in Part I of this study are utilized to improve the accuracy of performance predictions, and the concept of volume averaging is used to abstract the geometry and reduce computation time. The numerical model is validated by comparison with computational fluid dynamics (CFD) simulations and provides high accuracy and significantly lower computation time compared to existing numerical models. A preliminary optimization study is conducted and the advantage of using supercritical CO2 as a working fluid for energy systems is demonstrated.

Engineering↗

DOE OSTI2023

$\mathrm{DRAGON}$: A multi-GPU orbital-free density functional theory molecular dynamics simulation package for modeling of warm dense matter

As progress in electronic structure theoretical methods is made, ab initio molecular dynamics (MD) based on orbital-free density functional theory (OF-DFT) is becoming increasingly more successful at substituting the traditional, very accurate but computationally costly Kohn–Sham (KS) approach for simulations of matter at the challenging warm dense matter (WDM) regime. However, despite the significant cost alleviation of eliminating the dependence on the KS orbitals, OF-DFT MD runs require ~10 2 to 10 3 CPU cores running for days, or even weeks, for simulations of systems comprised of 10 2 to 10 3 atoms, depending on thermodynamic conditions. Here we present DRAGON, a multi-GPU OF-DFT MD code for fast and efficient simulations of WDM. With a relatively small allocation of resources (4 to 8 GPU devices) it can provide an order of magnitude speedup for simulations containing $\mathscr{O}$(10 4 ) atoms and target systems composed of $\mathscr{O}$(10 5 ) atoms at conditions within the WDM regime, which is currently outside the capabilities of CPU codes.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2023

Dependence of predicted bulk properties of hexagonal hydroxyapatite on exchange-correlation functional

Hydroxyapatite (HA) is the main component of human bones and teeth. HA has also been widely applied in various technological fields due to its unique properties. Reliable computational simulations for this critical material often require input parameters obtained from the first-principles density functional theory (DFT) calculations with appropriate exchange–correlation functionals. Previous DFT calculations are insufficient to verify the reliabilities of various functionals, particularly as they fail to assess predictions for multiple properties of the material. In this paper, we first select 18 different functionals to calculate geometric, elastic, electronic, and thermodynamic properties of hexagonal HA bulk crystal. We find that the results from optB86b-vdW and optB88-vdW functionals with dispersion corrections have the overall best agreement with available experimental data. Then, we choose optB88-vdW functional, as well as PBE functional without dispersion corrections as a comparison, to perform extensive first-principles DFT phonon calculations under the quasiharmonic approximation. Various thermodynamic properties (including phonon contributions to internal energy, entropy, and Helmholtz free energy) and thermal parameters (including thermal expansions of volume, thermal expansion coefficients, heat capacities, isothermal bulk moduli, etc.) versus temperature are consequently obtained. By comparing these quantities, we report that the results from optB88-vdW functional have significantly better agreement with available experimental data than those from PBE functional although the latter has been widely used in previous DFT calculations for HA-based material systems.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Stress-Dependent Chemo-Mechanical Performance of Amorphous Si Anodes for Li-Ion Batteries upon Lithiation

Alloying-type anodes are significantly governed by their chemo-mechanical performance during the electrochemical cycling. The reaction-induced huge volumetric change of these anodes may cause material degradation and failure under mechanical constraints. Here, we investigate the stress-dependent lithiation behavior of amorphous Si (a-Si) anodes using molecular dynamics simulations. It is indicated that a-Si anodes can sustain higher hydrostatic stress than biaxial/uniaxial ones without the occurrence of mechanical failure. Thermodynamic and electrochemical calculations demonstrate that although the lithiation procedure also affects the thermodynamic stability of a-Si anodes, it is mainly dominated by the external mean stresses. Compressive stress is confirmed to destabilize a-Si anodes and further trigger their capacity fading. Compared with our atomistic simulations, previous continuum models underestimate the open-cell potentials of a-Si anodes, due to their ignored large volumetric deformation at higher stresses and Li concentrations. Finally, this computational study provides the intensive atomic-level understanding of the stress-dependent lithiation behavior of a-Si anodes.

25 ENERGY STORAGE↗

DOE OSTI2021

Oxygenate Production from Plasma-Activated Reaction of CO 2 and Ethane

Upgrading ethane with CO 2 as a soft oxidant represents a desirable means of obtaining oxygenated hydrocarbons. This reaction is not thermodynamically feasible under mild conditions and has not been previously achieved as a one-step process. Non-thermal plasma was implemented as an alternative means of supplying energy to overcome activation barriers, leading to the production of alcohols, aldehydes, and acids as well as C 1 -C 5+ hydrocarbons under ambient pressure, with a maximum total oxygenate selectivity of 12%. A plasma chemical kinetic computational model was developed and found to be in good agreement with the experimental trends. Results from this study illustrate the potential to use plasma for the direct synthesis of value-added alcohols, acids, and aldehydes from ethane and CO 2 under mild conditions.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Proton-exchange induced reactivity in layered oxides for lithium-ion batteries

LiNi x Co y Mn 1-x-y O 2 (0 < x, y < 1, NCM) is the dominant positive material for the state-of-the-art lithium-ion batteries. However, the sensitivity of NCM materials to moisture makes their manufacturing, storage, transportation, electrode processing and recycling complicated. Although it is recognized that protons play a critical role in their structure stability and performance, proton exchange with Li + in NCM materials has not been well understood. Here, we employ advanced characterizations and computational studies to elucidate how protons intercalate into the layered structure of NCM, leading to the leaching of Li + and the formation of protonated NCM. It is found that protonation facilitates cation rearrangement and formation of impurity phases in NCM, significantly deteriorating structural stability. The adverse effects induced by protons become increasingly pronounced with a higher Ni content in NCM. Through a comprehensive investigation into the thermodynamics and kinetics of protonation, we discover that Li deficiencies in NCM materials can be resolved via solution process in the presence of Li + ions and controlled proton concentration. The underlying mechanism of relithiation is further explored through materials characterizations and kinetics modeling. This work provides crucial insights into controlling structural and compositional defects of Li-ion battery positive material in complicated processing environment.

batteries↗

DOE OSTI2020

Predicting synthesizable multi-functional edge reconstructions in two-dimensional transition metal dichalcogenides

Abstract Two-dimensional (2D) transition metal dichalcogenides (TMDCs) have attracted tremendous interest as functional materials due to their exceptionally diverse and tunable properties, especially in their edges. In addition to the conventional armchair and zigzag edges common to hexagonal 2D materials, more complex edge reconstructions can be realized through careful control over the synthesis conditions. However, the whole family of synthesizable, reconstructed edges remains poorly studied. Here, we develop a computational approach integrating ensemble-generation, force-relaxation, and electronic-structure calculations to systematically and efficiently discover additional reconstructed edges and screen their functional properties. Using MoS 2 as a model system, we screened hundreds of edge-reconstruction to discover over 160 reconstructed edges to be more stable than the conventional ones. More excitingly, we discovered nine new synthesizable reconstructred edges with record thermodynamic stability, in addition to successfully reproducing three recently synthesized edges. We also find our predicted reconstructed edges to have multi-functional properties—they show near optimal hydrogen evolution activity over the conventional edges, ideal for catalyzing hydrogen-evolution reaction (HER) and also exhibit half-metallicity with a broad variation in magnetic moments, making them uniquely suitable for nanospintronic applications. Our work reveals the existence of a wide family of synthesizable, reconstructed edges in 2D TMDCs and opens a new materials-by-design paradigm of ‘intrinsic’ edge engineering multifunctionality in 2D materials.

Hu, Guoxiang (ORCID:0000000329428564)↗

DOE OSTI2020

A higher-order finite element reactive transport model for unstructured and fractured grids

Abstract This work presents a new reactive transport framework that combines a powerful geochemistry engine with advanced numerical methods for flow and transport in subsurface fractured porous media. Specifically, the PhreeqcRM interface (developed by the USGS) is used to take advantage of a large library of equilibrium and kinetic aqueous and fluid-rock reactions, which has been validated by numerous experiments and benchmark studies. Fluid flow is modeled by the Mixed Hybrid Finite Element (FE) method, which provides smooth velocity fields even in highly heterogenous formations with discrete fractures. A multilinear Discontinuous Galerkin FE method is used to solve the multicomponent transport problem. This method is locally mass conserving and its second order convergence significantly reduces numerical dispersion. In terms of thermodynamics, the aqueous phase is considered as a compressible fluid and its properties are derived from a Cubic Plus Association (CPA) equation of state. The new simulator is validated against several benchmark problems (involving, e.g., Fickian and Nernst-Planck diffusion, isotope fractionation, advection-dispersion transport, and rock-fluid reactions) before demonstrating the expanded capabilities offered by the underlying FE foundation, such as high computational efficiency, parallelizability, low numerical dispersion, unstructured 3D gridding, and discrete fraction modeling.

58 GEOSCIENCES↗

DOE OSTI2025

Thermodynamic modeling of countercurrent chemical looping reverse water gas shift process for redox material screening

The reverse water gas shift (RWGS) reaction is a key pathway for CO 2 utilization, particularly within Power-to-X process chains aimed at sustainable fuel and chemical production. Countercurrent chemical looping (CL-RWGS) using non-stoichiometric oxides can overcome equilibrium limitations of conventional RWGS reactors, enabling significantly higher CO 2 conversions. However, modeling the limiting performance of such systems is challenging due to their multiphase nature and coupled spatial and temporal variation in chemical composition. In this work, we present a discretized batch equilibrium model that simulates CL-RWGS reactors as a series of localized equilibrium exchanges between gas and solid elements. The model is numerically stable, computationally efficient, and free of kinetic source terms, making it well-suited for parametric studies and system-level integration. It is validated against established convection–diffusion models and shown to predict reasonable upper bounds on experimental results. Application of the model to a range of oxygen carrier materials identifies cerium–zirconium solid solutions, particularly Ce 0.80 Zr 0.20 O 2 , as a promising class offering superior oxygen storage characteristics compared to state-of-the-art La 0.6 Sr 0.4 FeO 3 . This framework provides a robust platform for materials screening, reactor sizing, and performance optimization in chemical looping systems. The model implementation is available as open-source software to support further research and development.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗