Engineering Papers⌕ Search

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.

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

At least 559 records · Page 31
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Multilevel Computational Studies Reveal the Importance of Axial Ligand for Oxygen Reduction Reaction on Fe–N–C Materials

The systematic improvement of Fe–N–C materials for fuel cell applications has proven challenging, due in part to an incomplete atomistic understanding of the oxygen reduction reaction (ORR) under electrochemical conditions. Herein, a multilevel computational approach, which combines ab initio molecular dynamics simulations and constant potential density functional theory calculations, is used to assess proton-coupled electron transfer (PCET) processes and adsorption thermodynamics of key ORR intermediates. These calculations indicate that the potential-limiting step for ORR on Fe–N–C materials is the formation of the Fe III –OOH intermediate. They also show that an active site model with a water molecule axially ligated to the iron center throughout the catalytic cycle produces results that are consistent with the experimental measurements. In particular, reliable prediction of the ORR onset potential and the Fe(III/II) redox potential associated with the conversion of Fe III –OH to Fe II and desorbed H 2 O requires an axial H 2 O co-adsorbed to the iron center. The observation of a five-coordinate rather than four-coordinate active site has significant implications for the thermodynamics and mechanism of ORR. Importantly, these findings highlight the importance of solvent–substrate interactions and surface charge effects for understanding the PCET reaction mechanisms and transition-metal redox couples under realistic electrochemical conditions.

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

Evaluation of two-particle properties within finite-temperature self-consistent one-particle Green’s function methods: Theory and application to GW and GF2

One-particle Green’s function methods can model molecular and solid spectra at zero or non-zero temperatures. One-particle Green’s functions directly provide electronic energies and one-particle properties, such as dipole moment. However, the evaluation of two-particle properties, such as $\langle$S 2 $\rangle$ and $\langle$N 2 $\rangle$, can be challenging because they require a solution of the computationally expensive Bethe–Salpeter equation to find two-particle Green’s functions. We demonstrate that the solution of the Bethe–Salpeter equation can be completely avoided. Applying the thermodynamic Hellmann–Feynman theorem to self-consistent one-particle Green’s function methods, we derive expressions for two-particle density matrices in a general case and provide explicit expressions for GF2 and GW methods. Such density matrices can be decomposed into an antisymmetrized product of correlated one-electron density matrices and the two-particle electronic cumulant of the density matrix. Cumulant expressions reveal a deviation from ensemble representability for GW, explaining its known deficiencies. We analyze the temperature dependence of $\langle$S 2 $\rangle$ and $\langle$N 2 $\rangle$ for a set of small closed-shell systems. Interestingly, both GF2 and GW show a non-zero spin contamination and a non-zero fluctuation of the number of particles for closed-shell systems at the zero-temperature limit.

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

Biological Electron Transfer and Catalysis EFRC

The focus of the Biological Electron Transfer and Catalysis (BETCy) EFRC research is elucidating mechanisms of conversion of electrochemical potential into chemical bond energy and is organized into two integrated Thrusts: Electron Bifurcation Catalysis and Electron Transfer Kinetics. We are developing a collective knowledge of metalloenzymes as models for redox reactions by applying physical science and computational tools to characterize biochemical reactions catalyzed by multi- subunit enzymes harboring arrays of iron-sulfur clusters and flavin cofactors. Understanding these mechanisms is central to overcoming the thermodynamic barriers that currently limit production of reduced products and fuels. Biological systems have elegant strategies for converting electrochemical potential energy into chemical bond energy (e.g., C-H, H-H, and N-H) stored in reduced compounds that can serve as advanced biofuels. One significant limitation for the production of highly reduced compounds is that their production in natural and industrial processes relies on low oxidation-reduction potential chemicals as feedstocks.

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

Weak scaling of the contact distance between two fluctuating interfaces with system size

A pair of flat parallel surfaces, each freely diffusing along the direction of their separation, will eventually come into contact. If the shapes of these surfaces also fluctuate, then contact will occur when their centers-of-mass remain separated by a nonzero distance ℓ. An example of such a situation is the motion of interfaces between two phases at conditions of thermodynamic coexistence, and in particular the annihilation of domain wall pairs under periodic boundary conditions. Here we present a general approach to calculate the probability distribution of the contact distance ℓ and determine how its most likely value ℓ* depends on the surfaces' lateral size L. Using the Edward-Wilkinson equation as a model for interfaces, we demonstrate that ℓ* scales weakly with system size, i.e., the dependence of ℓ* on L for both (1+1)- and (2+1)-dimensional interfaces is such that lim L→∞ (ℓ*/L) = 0. In particular, for (2+1)-dimensional interfaces ℓ* is an algebraic function of logL, a result that is confirmed by computer simulations of slab-shaped domains formed under periodic boundary conditions. Overall, this weak scaling implies that such domains remain topologically intact until ℓ becomes very small compared to the lateral size of the interface, contradicting expectations from equilibrium thermodynamics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

An–imidophosphorane (An = U–Pu) bond covalency and proton-coupled electron transfer thermodynamics driven by orbital energy matching

A series of mid-actinide (An = U–Pu) tetrahomoleptic complexes supported by highly electron-donating imidophosphorane ligands, NPC ([NP t Bu(pyrr) 2 ] − , where t Bu = C(CH 3 ) 3 ; pyrr = pyrrolidinyl = N(C 4 H 8 )), are systematically investigated computationally and experimentally to elucidate the nature of actinide–ligand (An–L) covalency across the An 3+/4+/5+ oxidation states. Trends in An–L bonding and redox properties for these complexes, together with their protonated counterparts, are examined using orbital-, electron density-, and energy-decomposition-based methods. This integrated approach reveals progressively improved energy matching between α-spin An 5f and N im 2p orbitals with increasing atomic number and oxidation state, becoming particularly pronounced in the ligand-dominant π-bonding orbitals of An 4+ and An 5+ . In contrast to the An 3+ species, the enhanced An 5f π contributions in the higher-valent counterparts drive the increase in An–N im covalency for later An, thereby inverting the covalency trend to U < Np < Pu. Redistribution of electron density towards the An and N im atomic basins due to the growing energy-matching assisted covalency correlates with higher pKa values and increased N im –H bond dissociation free energies in protonated An 4+ complexes. Electron density at Nim in An 4+ shows a linear correlation with the p K a values calculated via the Bordwell equation. Calculations predict a cathodic shift of 0.84–1.00 V in the redox couples upon protonation, a trend validated when experimentally accessible. These findings demonstrate an increasing role of covalency driven by orbital energy matching from U to Pu in tuning the thermodynamic driving force for proton-coupled electron transfer in the An 5+ species.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2023

An autonomous laboratory for the accelerated synthesis of novel materials

To close the gap between the rates of computational screening and experimental realization of novel materials, we introduce the A-Lab, an autonomous laboratory for the solid-state synthesis of inorganic powders. This platform uses computations, historical data from the literature, machine learning (ML) and active learning to plan and interpret the outcomes of experiments performed using robotics. Over 17 days of continuous operation, the A-Lab realized 41 novel compounds from a set of 58 targets including a variety of oxides and phosphates that were identified using large-scale ab initio phase-stability data from the Materials Project and Google DeepMind. Synthesis recipes were proposed by natural-language models trained on the literature and optimized using an active-learning approach grounded in thermodynamics. Analysis of the failed syntheses provides direct and actionable suggestions to improve current techniques for materials screening and synthesis design. The high success rate demonstrates the effectiveness of artificial-intelligence-driven platforms for autonomous materials discovery and motivates further integration of computations, historical knowledge and robotics.

36 MATERIALS SCIENCE↗
DOE OSTI2024

The Pele Simulation Suite for Reacting Flows at Exascale

In this work, we present the Pele suite of software tools for compressible and incompressible reacting flows. The Pele suite leverages several different libraries, notably AMReX and SUNDIALS, to achieve performance portability on heterogeneous computing architectures across the supercomputing landscape. The Pele suite is comprised of PeleC, a compressible reacting flow block-structured adaptive mesh refinement solver, PeleLMeX, a low-Mach number reacting flow block-structured adaptive mesh refinement solver, Pele-Physics, a library for transport, thermodynamics, finite rate chemistry, soot, spray and radiation physics. The objective of this paper is (i) to present the code development efforts necessary to achieve highly effective and scalable applications for exascale machines and (ii) to detail the performance results of the Combustion-Pele project applications on Oak Ridge National Laboratory's Frontier. We show good weak and strong scaling results for both PeleC and PeleLMeX up to more than 50 billion cells on more than 4096 Frontier graphics processing unit nodes. We also present a capability demonstration simulation of a dual-fuel pulse compression ignition engine (six adaptive mesh refinement levels, and 60 billion cells or 2.1 trillion degrees of freedom) on Frontier, to date one of the largest simulations performed on the first exascale-class supercomputer.

adaptive mesh refinement↗
DOE OSTI2025

Dynamic Atomistic Polar Structure Underpins Ultrahigh Linear Electro-Optic Coefficient in Transparent Ferroelectric Ceramics

Transparent ferroelectrics with high linear electro-optic (EO) coefficients are critical for advanced electro-optical devices. However, achieving optical transparency in ferroelectric ceramics remains challenging due to visible light scattering caused by defects such as domain walls, grain boundaries, and pores. Here, we report the successful fabrication of transparent ferroelectric ceramics through innovative chemical composition design and an advanced two-step sintering process in the La-doped Pb(Mg 1/3 Nb 2/3 )O 3 –PbTiO 3 system. The optical transparency, which is near the theoretical upper limit, can be attributed to the wide band gap and the minimization of light scattering of defects. By minimizing porosity and engineering grain/domain sizes to differ significantly from the wavelengths of visible light, we suppress scattering, achieving optical transparency near the theoretical upper limit. Strikingly, these ceramics exhibit an ultrahigh linear EO coefficient of ∼1417 pm/V, over 65 times greater than that of LiNbO 3 single crystals, the current industry standard. We attribute this exceptional performance to dynamic atomistic polar structures within switchable, thermally stable domains, which enhance electronic polarization sensitivity. This mechanism is corroborated by dielectric spectroscopy, high-resolution transmission electron microscopy and simulation. Our findings offer insights into the design of cost-effective transparent materials with exceptional EO properties, paving the way for next-generation electro-optical devices.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Nanoengineering of non-aqueous liquid electrolyte solutions for future lithium metal batteries

Research and development of non-aqueous electrolyte solutions are essential for practical advancement towards the production of high-energy lithium metal batteries (LMBs). An ideal LMB electrolyte solution should enable highly efficient, uniform and prolonged lithium metal plating and stripping, preserve the electrodes’ electro(chemo)mechanical properties and ensure compatibility with all cell components. However, despite extensive research efforts, scientists have yet to achieve an electrolyte design that meets these requirements simultaneously. Here, by examining the nanoengineering aspects of various non-aqueous electrolyte solution designs, we elucidate the understanding of the nanoscale physicochemical and electrochemical processes taking place in LMBs, which are mainly governed by the thermodynamic and kinetic properties of the electrolyte system. We also explore emerging research directions and propose an accelerated, iterative framework that integrates nanoengineering principles with machine learning, high-throughput computation and experimentation to facilitate the development of next-generation non-aqueous electrolyte solutions for practical LMBs.

Weintz, Dominik↗
DOE OSTI2023

Chemical evolution in nitrogen shocked beyond the molecular stability limit

Evolution of nitrogen under shock compression up to 100 GPa is revisited via molecular dynamics simulations using a machine-learned interatomic potential. The model is shown to be capable of recovering the structure, dynamics, speciation, and kinetics in hot compressed liquid nitrogen predicted by first-principles molecular dynamics, as well as the measured principal shock Hugoniot and double shock experimental data, albeit without shock cooling. Our results indicate that a purely molecular dissociation description of nitrogen chemistry under shock compression provides an incomplete picture and that short oligomers form in non-negligible quantities. Finally, this suggests that classical models representing the shock dissociation of nitrogen as a transition to an atomic fluid need to be revised to include reversible polymerization effects.

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

Insight into pseudo-film-boiling in supercritical water under simplified interfacial conditions

The focus of this paper is the physics of supercritical water film boiling through the simulation of natural convection of supercritical water in a rectangular cavity. It has been widely believed that above a fluid's critical point, only one phase exists. However, recent research has indicated the existence of distinct gas-like and liquid-like regions separated by the Widom line, which is the locus of the maxima of the specific heat. Along the Widom line, density decreases sixfold, viscosity drops by a factor of two, while specific heat spikes by an order of magnitude. These variations, specifically in density and viscosity, produce a thick pseudo-interface and flow dynamics behavior akin to film-boiling, denoted here as pseudo-two-phase flow dynamics. When the cavity is heated from the bottom, a pseudo-film of gas-like fluid is formed at the heated wall, and the flow complexity evolves, characterized by the formation of multiple vortices. Results demonstrate the ability of modern conjugate heat transfer computational fluid dynamics to capture the pseudo two-phase regime, with the formation of a pseudo-film at the heated wall. Further, the boiling curve is studied, showing an increase in heat flux with an increase in the temperature differential between the wall and the pseudo-critical temperature. Furthermore, results also show a previously uncharacterized instability, where the convection cells in the cavity oscillate and can flip directions under certain conditions.

Computational fluid dynamics↗
DOE OSTI2020

Photonic Refrigeration from Time-Modulated Thermal Emission

We develop theoretical and computational formalisms to describe thermal radiation from temporally modulated systems. We demonstrate that such a modulation results in a photon-based active cooling mechanism. This mechanism has a high thermodynamic performance that can approach the Carnot limit. Our work points to exciting new avenues in active, time-modulated control of thermal emission for cooling and energy harvesting applications.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Feature Detection

Focal Area(s): This proposal aims to develop and evaluate statistical models and machine learning algorithms for detecting and tracking features in spatiotemporal remotely sensed data with uncertainty quantification. We focus a particular application on the detection of sea ice leads and ridges in the Arctic and use these key sea ice features for model calibration and to gain insight into the physics of sea ice thermodynamics and deformation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Laplace's method applied to functional integrals

In statistical thermodynamics, partition functions are required to evaluate the current state. It is typically impossible to analytically calculate the resulting configuration integrals, but it is sometimes possible to obtain approximations in the limit of steep potential energies. This approach directly leverages Laplace’s method for approximating integrals and has been quite successful in recent applications. Motivated by models with an infinite number of degrees of freedom, this approach should be extended to the statistical thermodynamics of fields. To that end, Laplace’s method is developed for functional integrals.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Cooperativity and Metal–Linker Dynamics in Spin Crossover Framework Fe(1,2,3-triazolate) 2

Cooperative interactions are responsible for the useful properties of spin crossover (SCO) materials–large hysteresis windows, critical temperatures near room temperature, and abrupt transitions–with hybrid framework materials exhibiting the greatest cooperativity and hysteresis of all SCO systems. However, little is known about the chemical origin of cooperativity in frameworks. Here, we present a combined experimental–computational approach for identifying the origin of cooperativity in the metal–organic framework (MOF) Fe(1,2,3-triazolate) 2 (Fe(TA) 2 ), which exhibits the largest known hysteresis window of all SCO materials and unusually high transition temperatures, as a roadmap for understanding the manipulation of SCO behavior in general. Variable-temperature vibrational spectroscopy provides evidence that “soft modes” associated with dynamic metal–linker bonding trigger the cooperative SCO transition. Thermodynamic analysis also confirms a cooperativity magnitude much larger than those of other SCO systems, while electron density calculations of Fe(TA) 2 support previous theoretical predictions that large cooperativity arises in materials where SCO produces considerable differences in metal–ligand bond polarities between different spin states. Taken together, this combined experimental–computational study provides a microscopic basis for understanding cooperative magnetism and highlights the important role of dynamic bonding in the functional behavior of framework materials.

36 MATERIALS SCIENCE↗
DOE OSTI2023

New High-Entropy Perovskite Oxides with Increased Reducibility and Stability for Thermochemical Hydrogen Generation

This project aims to design, synthesize, and test a transformative class of High-Entropy Perovskite Oxides (HEPOs) as redox oxides to enable thermochemical hydrogen generation with improved stability, kinetics, and efficiency. These developed HEPOs are expected to demonstrate improved kinetics with oxygen surface exchange coefficient ( k > 7.5×10 -4 cm/s) in Budget Period (BP) 1, retain its structural stability in a broad range of oxygen non-stoichiometry (Δδ > 0.15) at a low operating reduction temperature of T red < 1400°C in BP 2, and deliver a H 2 yield of over 400 µmol per gram of oxide and high stability with less than 20% degradation after at least 50 cycles in BP 3. This project is feasible due to the unique thermodynamic properties (simultaneously increased reducibility and phase stability) and kinetic characters (stability against particle coarsening and potentially enhanced oxygen transport and surface reaction kinetics) of such HEPOs, and it is enabled by a unique active learning computational design approach. Computational studies have been conducted to investigate the oxygen vacancy formation in complex perovskite systems. * Accurate prediction of V O .. concentration with disordered A-site cations in Fe-based high-entropy perovskite oxides * Combined MC/DFT computation elucidates the mechanism of Co preference on the redox due to the strain introduced by local distortion. In this project, we explored a large number (~150) of perovskite compositions, which are listed in Tables 2 – 4). * All perovskite specimens have been synthesized through a high-throughput high-energy ball milling process, followed by the conventional sintering process. * XRD, SEM/EDS and TGA were performed to confirm the crystal structure, phase homogeneity and oxygen non-stoichiometry for compositionally complex perovskite oxides (CCPOs). * 110 compositions show single-phase from XRD * Unusual aliovalent doping effects in medium-entropy perovskite compositions. * V-shape relation between Δδ vs. x (La 1-x Sr x )(Mn 1/3 Fe 1/3 Ti 1/3 )O 3 * The highest reported hydrogen production for the CCPOs made in this project ( T re = 1350 ºC 30 min, T Ox = 1100 ºC 30 min) * B-site mixing (La 0.8 Sr 0.2 )(Mn 0.2 Fe 0.2 Co 0.4 Al 0.2 )O 3 : 89.97 ± 2.73 mmol H2 /mol oxide (395 ± 10 μmol/g oxide ) (i) No phase transformation during reactions when Co molar ratio is less than 61% (ii) Balance between intrinsic kinetics (oxygen surface exchange) and thermodynamics (oxygen non-stoichiometry) (iii) Preference of Co identified by in-situ XPS * A-site mixing (La 1/6 Pr 1/6 Nd 1/6 Gd 1/6 Ba 1/6 Sr1/6)MnO 3 : 98.48 mmol H2 /mol oxid e (~415 μmol/g oxide )

08 HYDROGEN↗
DOE OSTI2025

Simultaneous inference of equation of state parameters and unknown data errors with uncertainty quantification via hierarchical Bayesian posterior maximization

Equations of state (EOSs) are a key component in running hydrodynamic simulations as they relate the thermodynamic states for the material. The Davis reactants EOS is commonly used for modeling high explosives (HEs), and the EOS model parameters are calibrated using material specific data. The calibrations are often performed with uncertainty quantification via Bayesian inference to account for uncertainty in the data and generate ensembles of likely parameters. However, there are relatively few HE data sets to use for calibration and many are historical and lack error information. In this work, we simultaneously calibrate the Davis reactants EOS model parameters and unknown data error terms for the high explosive PBX 9501. To quantify the uncertainty in the models and the data, we use a Bayesian framework for the calibration and compute the hierarchical Bayesian posterior distribution with both a posteriori maximization approach and Markov Chain Monte Carlo. In general, we find that, given our assumptions, the two approaches result in similar calibrated parameters, posterior covariance matrices, and insights about the parameters but that the posterior maximization requires far less computational resources.

97 MATHEMATICS AND COMPUTING↗