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 1,009 records · Page 56
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Reorganization Energies for Interfacial Proton-Coupled Electron Transfer to a Water Oxidation Catalyst

The reorganization energy (λ) for interfacial electron transfer (ET) and proton-coupled electron transfer (PCET) from a conductive metal oxide (In 2 O 3 :Sn, ITO) to a surface-bound water oxidation catalyst were extracted from kinetic data measured as a function of the thermodynamic driving force. Visible light excitation resulted in rapid excited-state injection (k inj > 10 8 s -1 ) to the ITO, which photo-initiated the two interfacial reactions of interest. The rate constants for both reactions increased with the driving force, −ΔG°, to a saturating limit, k max , with rate constants consistently larger for ET than for PCET. Marcus-Gerischer analysis of the kinetic data provided the reorganization energy for interfacial PCET (0.90 ± 0.02 eV) and ET (0.40 ± 0.02 eV), respectively. The magnitude of k max for PCET was found to decrease with pH, behavior that was absent for ET. Both the decrease in k max and the larger reorganization energy for an unwanted competing PCET reaction from the ITO to the oxidized catalyst showcases a significant kinetic advantage for driving solar water oxidation at high pH. Computational analysis revealed a larger inner-sphere reorganization energy contribution for PCET than for ET arising from a more significant change in the Ru–O bond length for PCET. Extending Marcus-Gerischer theory to PCET by including the excited electron-proton vibronic states and the proton donor-acceptor motion provided an apparent reorganization energy of 1.01 eV. Furthermore, this study demonstrates that the Marcus-Gerischer theory initially developed for ET can be reliably extended to PCET for quantifying and interpreting reorganization energies observed experimentally.

Catalysts Kinetic parameters↗
DOE OSTI2024

Divergent Bimetallic Mechanisms in Copper(II)-Mediated C–C, N–N, and O–O Oxidative Coupling Reactions

Copper-catalyzed aerobic oxidative coupling of diaryl imines provides a route for conversion of ammonia to hydrazine. Here the present study uses experimental and density functional theory computational methods to investigate the mechanism of N–N bond formation, and the data support a mechanism involving bimolecular coupling of Cu-coordinated iminyl radicals. Computational analysis is extended to Cu II -mediated C–C, N–N, and O–O coupling reactions involved in the formation of cyanogen (NC–CN) from HCN, 1,3-butadiyne from ethyne (i.e., Glaser coupling), hydrazine from ammonia, and hydrogen peroxide from water. The results reveal two different mechanistic pathways. Heteroatom ligands with an uncoordinated lone pair (iminyl, NH 2 , OH) undergo charge transfer to Cu II , generating ligand-centered radicals that undergo facile bimolecular radical–radical coupling. Ligands lacking a lone pair (CN and CCH) form bridged binuclear diamond-core structures that undergo C–C coupling. This mechanistic bifurcation is rationalized by analysis of spin densities in key intermediates and transition states, as well as multiconfigurational calculations. Radical–radical coupling is especially favorable for N–N coupling owing to energetically favorable charge transfer in the intermediate and thermodynamically favorable product formation.

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

CALPHAD Uncertainty Quantification and TDBX

CALPHAD uncertainty quantification (UQ) is the foundation of materials design with quantified confidence. We report a framework and software packages to enable CALPHAD UQ assessment and calculation using commercial CALPHAD software (Thermo-Calc). This Bayesian inference framework is coupled with a Markov chain Monte Carlo algorithm to establish uncertainty traces with a given thermodynamic database file (TDB) and corresponding experimental data points. This general framework is demonstrated with the Ni–Cr binary system. The algorithm is firstly validated on synthetic data with known ground truth. Then it is applied to real experimental data to generate posterior traces. We develop a file format named TDBX, which provides a single source of truth by combining the original TDB content and the traces for each assessed Gibbs energy parameter. CALPHAD UQ calculations are performed based on the TDBX file, from which uncertainties for phase boundaries, enthalpy curves, and solidification range are collected as examples of basic design parameters. This TDBX file with corresponding scripts are made open-source. Finally, the combination of CALPHAD UQ assessments and calculations connected by TDBX supports uncertainty-assisted modeling, enabling the integrated application of modern design with uncertainty methodologies to computational materials design.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Investigating two-dimensional adjoint QCD on the lattice

We present our investigations of SU(N) adjoint QCD in two dimensions with one Majorana fermion on the lattice. We determine the relevant parameter range for the simulations with Wilson fermions and present results for Polyakov loop, chiral condensate, and string tension. In the theory with massive fermions, all observables we checked show qualitative agreement between numerical lattice data and theory, while the massless limit is more subtle since chiral and non-invertible symmetry of the continuum theory are explicitly broken by lattice regularization. In thermal compactification, we observe N perturbative vacua for the holonomy potential at high-T with instanton events connecting them, and a unique vacuum at low-T. At finite-N, this is a cross-over and it turns to a phase transition at large-N thermodynamic limit. In circle compactification with periodic boundary conditions, we observe a unique center-symmetric minimum at any radius. In continuum, the instantons in the thermal case carry zero modes (for even N) and indeed, in the lattice simulations, we observe that chiral condensate is dominated by instanton centers, where zero modes are localized. We present lattice results on the issue of confinement vs. screening in the theory and comment on the roles of chiral symmetry and non-invertible symmetry.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Studies of Low and Intermediate Temperature Oxidation of Propane up to 100 Atm in a Supercritical-Pressure Jet-Stirred Reactor

Here, the low and intermediate temperature oxidation of propane has been investigated by using a novel supercritical pressure jet stirred reactor (SP-JSR) with and without 20% CO 2 additions at fuel lean and rich conditions at 10 and 100 atm and 500–1000 K. The mole fractions of C 3 H 8 , O 2 , CO, CO 2 , CH 2 O, C 2 H 4 , CH 3 CHO, and C 3 H 6 were quantified by using a micro-gas chromatograph (µ-GC). The experiment showed that different from that of 10 atm, at 100 atm only a weak negative temperature coefficient (NTC) behavior was observed because of the significant shift of the intermediate temperature HO 2 chemistry to lower temperature. In addition, at 100 atm, existing models in literatures could successfully capture the onset temperatures of the low and intermediate chemistry, while under-predict the fuel oxidation quantitatively and fail to capture the NTC behavior between 650 and 780 K at both fuel lean and rich conditions. Similar discrepancy was observed in studies of n-butane and dimethyl ether (DME) oxidations in literatures, implying that there existed large uncertainties in hierarchy model development of fuels with low temperature chemistries at extremely high pressures. Reaction pathways and sensitivity analyses showed that RO 2 competing reactions through (P1) RO 2 = QOOH, (P2) RO 2 = C 3 H 6 + HO 2 , (P3) RO 2 + CH 2 O/HO 2 = RO 2 H + HCO / O 2 dominated the low and intermediate temperature chemistries, followed by HO 2 / H 2 O 2 chemistry at 100 atm, which differed from the dominant pathway through QOOH consumption reactions at lower pressures. Especially, P3 is a new pathway of RO 2 consumption at high pressures, which was not observed in importance at low pressures. Special attention should be paid to the accurate computations of n-C 3 H 7 O 2 / i-C 3 H 7 O 2 + CH 2 O and n-C 3 H 7 O 2 / i-C 3 H 7 O 2 + in the P3 pathway and n-C 3 H 7 O 2 / i-C 3 H 7 O 2 decomposition reactions in the P2 pathway at high pressures.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2022

Thermodynamics of high-pressure ice phases explored with atomistic simulations

Most experimentally known high-pressure ice phases have a body-centred cubic (bcc) oxygen lattice. Our large-scale molecular-dynamics simulations with a machine-learning potential indicate that, amongst these bcc ice phases, ices VII, VII' and X are the same thermodynamic phase under different conditions, whereas superionic ice VII" has a first-order phase boundary with ice VII'. Moreover, at about 300 GPa, the transformation between ice X and the Pbcm phase has a sharp structural change but no apparent activation barrier, whilst at higher pressures the barrier gradually increases. Our study thus clarifies the phase behaviour of the high-pressure ices and reveals peculiar solid–solid transition mechanisms not known in other systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Comparison of Equilibrium and Nonequilibrium Approaches for Relative Binding Free Energy Predictions

Alchemical relative binding free energy calculations have recently found important applications in drug optimization. A series of congeneric compounds are generated from a preidentified lead compound, and their relative binding affinities to a protein are assessed in order to optimize candidate drugs. While methods based on equilibrium thermodynamics have been extensively studied, an approach based on nonequilibrium methods has recently been reported together with claims of its superiority. However, these claims pay insufficient attention to the basis and reliability of both methods. Here we report a comparative study of the two approaches across a large data set, comprising more than 500 ligand transformations spanning in excess of 300 ligands binding to a set of 14 diverse protein targets. Ensemble methods are essential to quantify the uncertainty in these calculations, not only for the reasons already established in the equilibrium approach but also to ensure that the nonequilibrium calculations reside within their domain of validity. If and only if ensemble methods are applied, we find that the nonequilibrium method can achieve accuracy and precision comparable to those of the equilibrium approach. Compared to the equilibrium method, the nonequilibrium approach can reduce computational costs but introduces higher computational complexity and longer wall clock times. There are, however, cases where the standard length of a nonequilibrium transition is not sufficient, necessitating a complete rerun of the entire set of transitions. This significantly increases the computational cost and proves to be highly inconvenient during large-scale applications. Our findings provide a key set of recommendations that should be adopted for the reliable implementation of nonequilibrium approaches to relative binding free energy calculations in ligand-protein systems.

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

Nickel B-site substitution in bulk Sr 1-x Ca x FeO 3 perovskite oxygen carriers: Benefits and limitations

We report oxygen (O 2 ) storage materials often rely on the presence of cobalt (Co) to reduce the thermodynamic penalty and increase the kinetics necessary for efficient O 2 storage and release. In this work, we investigate nickel (Ni) as an alternative B-site dopant in Sr 1-x Ca x FeO 3 to identify Co-free carriers that still show improved kinetics at low temperatures. In fact, we show a substantial increase in the reversible O2 release rate through mild Ni B-site substitution (y = 0.06) in select Sr 1-x Ca x Fe 1-y Ni y O 3 systems at 400 to 500 °C, reaching 2.00 wt.% O 2 release up to approximately 75% faster than Ni-free systems. To explain the role of Ni in these systems, we use density functional theory to calculate the O 2 vacancy (V O ) formation energy from separate metal-oxygen (M-O) bonding and relaxation components. We computationally show elongated Ni-O bonds are directly responsible for the decrease in V O upon Ni substitution.

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

Theoretical Investigation of Delafossite-Cu 2 ZnSnO 4 as a Promising Photovoltaic Absorber

In the quest for efficient and cost-effective photovoltaic absorber materials beyond silicon, considerable attention has been directed toward exploring alternatives. One such material, zincblende-derived Cu2ZnSnS4 (CZTS), has shown promise due to its ideal band gap size and high absorption coefficient. However, challenges such as structural defects and secondary phase formation have hindered its development. In this study, we examine the potential of another compound, Cu2ZnSnO4 (CZTO), with a similar composition to CZTS as a promising alternative. Employing ab initio density function theory (DFT) calculations in combination with an evolutionary structure prediction algorithm, we identify that the crystalline phase of delafossite structure is the most stable among the 900 (meta)stable CZTO. Its thermodynamic stability at room temperature is also confirmed by the molecular dynamics study. Excitingly, this new phase of CZTO displays a direct band gap where the dipole-allowed transition occurs, making it a strong candidate for efficient light absorptions. Furthermore, the estimation of spectroscopic limited maximum efficiency (SLME) directly demonstrates the high potential of delafossite-CZTO as a photovoltaic absorber. Our numerical results suggest that delafossite-CZTO holds promise for future photovoltaic applications.

14 SOLAR ENERGY↗
DOE OSTI2022

express: Extensible, high-level workflows for swifter ab initio materials modeling

In this work, we introduce an open-source Julia project, express, an extensible, lightweight, high-throughput, high-level workflow framework that aims to automate ab initio calculations for the materials science community. express is shipped with well-tested workflow templates, including structure optimization, equation of state (EOS) fitting, phonon spectrum (lattice dynamics) calculation, and thermodynamic property calculation in the framework of the quasi-harmonic approximation (QHA). It is designed to be highly modularized so that its components can be reused across various occasions, and customized workflows can be built on top of that. Users can also track the status of workflows in real-time, and rerun failed jobs thanks to the data lineage feature express provides. Finally, two working examples, i.e., all workflows applied to lime and akimotoite, are also presented in the code and this paper.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Velocity Measurements in Channel Gas Flows in the Slip Regime by means of Molecular Tagging Velocimetry

Direct measurements of the slip velocity in rarefied gas flows produced by local thermodynamic non-equilibrium at the wall represent crucial information for the validation of existing theoretical and numerical models. In this work, molecular tagging velocimetry (MTV) by direct phosphorescence is applied to argon and helium flows at low pressures in a 1-mm deep channel. MTV has provided accurate measurements of the molecular displacement of the gas at average pressures of the order of 1 kPa. To the best of our knowledge, this work reports the very first flow visualizations of a gas in a confined domain and in the slip flow regime, with Knudsen numbers up to 0.014. MTV is cross-validated with mass flowrate measurements by the constant volume technique. The two diagnostic methods are applied simultaneously, and the measurements in terms of average velocity at the test section are in good agreement. Moreover, preliminary results of the slip velocity at the wall are computed from the MTV data by means of a reconstruction method

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

Thermodynamics and transport in molten chloride salts and their mixtures

Relationship between thermophysical properties and phonon mean free path. Heat capacity, viscosity, and thermal conductivity in ionic liquids decrease as mean free path decreases and dynamics become less “solid-like” and collect motion diminishes.

Cockrell, C. [Nuclear Futures Institute, Bangor Un↗
DOE OSTI2020

Extended-magnetohydrodynamics in under-dense plasmas

Extended-magnetohydrodynamics (MHD) transports magnetic flux and electron energy in high-energy-density experiments, but individual transport effects remain unobserved experimentally. Two factors are responsible in defining the transport: electron temperature and electron current. Here, each electron energy transport term has a direct analog in magnetic flux transport. To measure the thermally driven transport of magnetic flux and electron energy, a simple experimental configuration is explored computationally using a laser-heated pre-magnetized under-dense plasma. Changes to the laser heating profile precipitate clear diagnostic signatures from the Nernst, cross-gradient-Nernst, anisotropic conduction, and Righi-Leduc heat-flow. With a wide operating parameter range, this configuration can be used in both small and large scale facilities to benchmark MHD and kinetic transport in collisional/semi-collisional, local/non-local, and magnetized/unmagnetized regimes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Integration Roadmap for Multi-Scale, Multi-Physics Mass Accountancy in Molten Salt Reactors

The objective of the integration roadmap is to describe the necessary physics required and a plan for modeling and simulation approach for predicting mass accountancy in molten salt reactors on an engineering scale. The thermophysical properties and the underlying thermodynamics are fundamental inputs. Therefore, the modeling will span length scales from first principles calculations to the engineering scale. The intention is to predict where material accumulates in a reactor core and loop and to understand perturbations on the systems level, for example the downstream effects from a turbine failure

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

Compositionally Tuning Electron Transfer from Photoexcited Core/Shell Quantum Dots via Cation Exchange

It is critical to find methods to control the thermodynamic driving force for photoexcited charge transfer from quantum dots (QDs) and explore how this affects charge transfer rates, since the efficiency of QD-based photovoltaic and photocatalysis technologies depends on both this rate and the associated energetic losses. In this work, we introduce a single-pot shell growth and Cu-catalyzed cation exchange method to synthesize Cd x Zn 1-x Se/Cd y Zn 1-y S QDs with tunable driving forces for electron transfer. Functionalizing them with two molecular electron acceptors—naphthalenediimide (NDI) and anthraquinone (AQ)—allowed us to probe nearly 1 eV of driving forces. For AQ, at lower driving forces, we find that higher Zn content results in a 130-fold increase of electron transfer rate constants. However, at higher driving forces electron transfer dynamics are unaltered. Here, the data are understood using an Auger-assisted electron transfer model and analyzed with computational work to determine approximate binding geometries of these electron acceptors. Our work provides a method to tune QD reducing power and produces useful metrics for optimizing QD charge transfer systems that maximize rates of electron transfer while minimizing energetic losses.

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

Efficacy of the radial pair potential approximation for molecular dynamics simulations of dense plasmas

Macroscopic simulations of dense plasmas rely on detailed microscopic information that can be computationally expensive and is difficult to verify experimentally. In this work, we delineate the accuracy boundary between microscale simulation methods by comparing Kohn–Sham density functional theory molecular dynamics (KS-MD) and radial pair potential molecular dynamics (RPP-MD) for a range of elements, temperature, and density. By extracting the optimal RPP from KS-MD data using force matching, we constrain its functional form and dismiss classes of potentials that assume a constant power law for small interparticle distances. Our results show excellent agreement between RPP-MD and KS-MD for multiple metrics of accuracy at temperatures of only a few electron volts. The use of RPPs offers orders of magnitude decrease in computational cost and indicates that three-body potentials are not required beyond temperatures of a few eV. Due to its efficiency, the validated RPP-MD provides an avenue for reducing errors due to finite-size effects that can be on the order of ∼20%.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Nanotomography for Quantitative 3D Particle Reconstruction

Particulates are ubiquitous across fuel cycle operations and carry critical information about particle formation, processing, and potential proliferation-related activities. Traditional analytical techniques, including micro-Raman spectroscopy and standard electron microscopy, are often limited in spatial resolution or dimensionality, particularly when used to examine metallic or submicron-scale features. Understanding particle morphology, phase distribution, and internal porosity is essential for constraining formation conditions, thermodynamic environments, and material transport behavior. In this report, we demonstrate the application of plasma focused ion beam nanotomography to reconstruct micron-scale particulates at nanoscale resolution. Using high-resolution backscattered electron imaging and Avizo software, we obtained 3D reconstructions that enabled quantitative analysis of particle morphology, phase composition, and internal voids. Representative examples include a Ta particle with a large central void and a composite particle with embedded tetrahedral crystalline structures. These reconstructions reveal structural and compositional details that are inaccessible through conventional 2D imaging. The results demonstrate that nanotomography provides both qualitative and quantitative insights into particle formation and behavior. Using nanotomography, porosity and phase distributions can be quantified to inform models of particle density, transport, and solidification conditions. Beyond technical insights, the workflow developed here establishes a transferable capability for analyzing heterogeneous particles and has potential applications in bulk materials studies via x-ray computed tomography or other volumetric imaging modalities. Ongoing efforts are focused on optimizing the workflow to process multiple particles simultaneously, increasing throughput and statistical robustness. Overall, this work illustrates the power of nanotomography as a tool for connecting particulate morphology to formation mechanisms, composition, and transport, thereby strengthening analytical capabilities for nuclear forensics, fuel cycle analysis, and related scientific investigations.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Modeling of Supercritical CO2 Shell-and-Tube Heat Exchangers Under Extreme Conditions. Part I: Correlation Development

Abstract High-temperature supercritical CO2 Brayton cycles are promising possibilities for future stationary power generation and hybrid electric propulsion applications. Heat exchangers are critical components in supercritical CO2 thermal cycles and require accurate correlations and comprehensive performance modeling under extreme temperatures and pressures. In this paper (Part I), new Colburn and friction factor correlations are developed to quantify shell-side heat transfer and friction characteristics of flow within heat exchangers in the shell-and-tube configuration. Using experimental and computational fluid dynamics (CFD) data sets from existing literature, multivariate regression analysis is conducted to achieve correlations that capture the effect of multiple critical geometric parameters. These correlations offer superior accuracy and versatility as compared to previous studies and predict the thermohydraulic performance of about 90% of the existing experimental and CFD data within ±15%. Supplementary thermohydraulic performance data are acquired from CFD simulations with supercritical CO2 as working fluid to validate the developed correlations and demonstrate its capability to be applied to supercrtical CO2 heat exchangers.

Engineering↗