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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
NASA NTRS1972

Tables of critical-flow functions and thermodynamic properties for methane and computational procedures for both methane and natural gas

Procedures for calculating the mass flow rate of methane and natural gas through nozzles are given, along with the FORTRAN 4 subroutines used to make these calculations. Three sets of independent variables are permitted in these routines. In addition to the plenum pressure and temperature, the third independent variable is either nozzle exit pressure, Mach number, or temperature. A critical-flow factor that becomes a convenient means for determining the mass flow rate of methane through critical-flow nozzles is tabulated. Other tables are included for nozzle throat velocity and critical pressure, density, and temperature ratios, along with some thermodynamic properties of methane, including compressibility factor, enthalpy, entropy, specific heat, specific-heat ratio, and speed of sound. These tabulations cover a temperature range from 120 to 600 K and pressures to 3 million N/sq m.

Johnson, R. C.↗
DOE OSTI2022

Environmental and Energy Research at the Texas Center for Advanced Scientific Computing and Modeling (CASCaM)

The funding enabled the creation of the (University of North) Texas Center for Advanced Scientific Computing and Modeling (CASCaM). Due to its success, CASCaM has continued to grow beyond the initial team focused on the biological and environmental research discussed herein, with investments by the University of North Texas in new faculty and researchers who contribute to CASCaM research. CASCaM is now a sustainable effort, with researchers across multiple fields including chemistry, physics, biology, mechanical engineering, materials science, and other fields, with multiple researchers engaged in a broad variety of research, including other DOE projects. Over 20 PI’s and ~100 researchers are engaged in CASCaM research. The overall goal of the research effort was to provide greater understanding about carbon capture, storage, and utilization. As part of this effort, the development of computational methods and identification of novel strategies to support the CO2 and CO studies were necessary. Three primary themes were investigated: carbon dioxide and CO catalysis; protein- and amino acid-based strategies towards carbon dioxide capture, storage, and utilization; and, finally, the methodology improvements for modeling carbon dioxide chemistry

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1981

Calculation of the thermodynamic properties of a mixture of gases as a function of temperature and pressure

The evaluation of the thermodynamic properties of a gas mixture can be performed using a generalized correlation which makes use of the second virial coefficient. This coefficient is based on statistical mechanics and is a function of temperature and composition, but not of pressure. The method provides results accurate to within 3 percent for gases which are nonpolar or only slightly polar. When applied to highly polar gases, errors of 5 to 10 percent may result. For gases which associate, even larger errors are possible. The sequences of calculations can be routinely programmed for a digital computer. The thermodynamic properties of a mixture of neon, argon and ethane were calculated by such a program. The result will be used for the design of the gas replenishment system for the Energetic Gamma Ray Experiment Telescope.

Colon, G.↗
DOE OSTI2022

Calculating adsorption isotherms using the two-phase thermodynamic method and molecular dynamics simulations

We describe the calculation of adsorption isotherms from molecular dynamics simulations based on the two-phase thermodynamic (2PT) model. The 2PT model developed for bulk fluid phases treats the gas-like components as hard spheres (HSs), which correctly recovers the limiting behaviors of unconfined fluids. We showed that this treatment, however, does not always lead to the correct zero-loading behavior in strongly confining systems. For methane adsorption into zeolite MFI, the HS reference state underestimates entropy by up to 20% at low loadings and leads to an order-of-magnitude increase in the adsorption onset pressure. To fix these issues, we propose the use of ideal adsorbed gas (IAG) as the gas reference model, the properties of which can be computed using the Widom insertion method on an empty adsorbent. We further describe three routes to compute adsorption isotherms from the Helmholtz free energy at different loadings. Comparing against established Monte Carlo (MC) methods, we found that the adsorption isotherms obtained using the IAG reference state agrees to within 40%, which corresponds to deviations of <5% in adsorption free energy. The isotherms calculated using the HS reference state underestimate the adsorption uptake at low to medium loadings in strongly confining systems, but its accuracy improves at higher loadings and as the pore size increases relative to the sorbate diameter. The methods described here provide an alternative approach for computing adsorption isotherms when MC simulations in an open ensemble are undesirable and enable a direct comparison of computed adsorption thermodynamics with experiments.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

H2OI95: A Stand-Alone Fortran Code for Evaluating the IAPWS-95 Equation-of-State Model for Water (Rev. 1)

H2OI95 is a stand-alone Fortran code for evaluating the IAPWS-95 equation-of-state model for the thermodynamic properties of water. It also evaluates the corresponding thermochemical properties of water consistent with the CODATA recommendations. The IAPWS-95 model is based on a model equation for the dimensionless Helmholtz energy for which the primary variables are the inverse reduced temperature (τ = T cr /T ) and reduced density (δ = ρ/ρ cr ). Here T is the absolute temperature (K), ρ is density (kg/m³), and the subscript “ cr ” refers to the critical point of water (647.096 K and 22.064 MPa pressure in this model, for which ρ cr is 322 kg/m³). The code solves four basic types of problems, distinguished by the specified inputs: 1. Temperature (K) and density (ρ kg/m³) or reduced density (δ) 2. Temperature (K) and pressure (MPa). 3. Temperature (K) on the saturation (liquid-vapor equilibrium) curve 4. Pressure (MPa) on the saturation curve. Each type of problem is run using a corresponding input (text) file. All but the first type of problem require iteration. For example, to solve for desired temperature and pressure, the reduced density must be adjusted to give the desired pressure. Iteration is accomplished using the Newton-Raphson method, though the secant method is also used in solving the fourth type of problem. For the last three types of problems, H2OI95 has been used to conduct numerical studies of convergence and the problem of multiple numerical solutions, some of which are not physically valid. Obtaining valid results depends mainly on appropriate choice of starting values for the density. The default values used in H2OI95 appear to consistently lead to generally desired results. With modification (not addressed here), H2OI95 can be used to support SUPCRT92 and similar codes that compute chemical thermodynamic properties of species and reactions over a wide range of temperature and pressure (273.16-1273K and 0-1000 MPa). Version 1.1 of this software calculates an extended set of water properties, including some derivable solely from IAPWS-95, such as the isobaric thermal expansivity and the isothermal compressibility, and still other properties based on supplementary models sanctioned by IAPWS, including the dynamic viscosity, thermal conductivity, surface tension, and static dielectric constant, and various other properties (notably Debye-Hückel parameters) derivable from the dielectric constant and its partial derivatives with respect to temperature and pressure. A coding error in Version 1.0 was identified while testing the implementation of the viscosity model. This error is described in Appendix C. It mainly affects calculation of “higher order” properties such as compressibility and heat capacity near the critical point.

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

Homologous Ladder Cyclohexasilanes

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

36 MATERIALS SCIENCE↗
DOE OSTI2026

A modified scheil approach for nucleation-dependent solidification pathways

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

Alloy design↗
DOE OSTI2021

Dynamics of disordered mechanical systems with large connectivity, free probability theory, and quasi-Hermitian random matrices

Disordered mechanical systems with high connectivity represent a limit opposite to the more familiar case of disordered crystals. Individual ions in a crystal are subjected essentially to nearest-neighbor interactions. In contrast, the systems studied in this paper have all their degrees of freedom coupled to each other. Thus, the problem of linearized small oscillations of such systems involves two full positive-definite and non-commuting matrices, as opposed to the sparse matrices associated with disordered crystals. Consequently, the familiar methods for determining the averaged vibrational spectra of disordered crystals, introduced many years ago by Dyson and Schmidt, are inapplicable for highly connected disordered systems. In this paper we apply random matrix theory (RMT) to calculate the averaged vibrational spectra of such systems, in the limit of infinitely large system size. At the heart of our analysis lies a calculation of the average spectrum of the product of two positive definite random matrices by means of free probability theory techniques. We also show that this problem is intimately related with quasi-hermitian random matrix theory (QHRMT), which means that the ‘hamiltonian’ matrix is hermitian with respect to a non-trivial metric. This extends ordinary hermitian matrices, for which the metric is simply the unit matrix. The analytical results we obtain for the spectrum agree well with our numerical results. The latter also exhibit oscillations at the high-frequency band edge, which fit well the Airy kernel pattern. We also compute inverse participation ratios of the corresponding amplitude eigenvectors and demonstrate that they are all extended, in contrast with conventional disordered crystals. Finally, we compute the thermodynamic properties of the system from its spectrum of vibrations. In addition to matrix model analysis, we also study the vibrational spectra of various multi-segmented disordered pendula, as concrete realizations of highly connected mechanical systems. A universal feature of the density of vibration modes, common to both pendula and the matrix model, is that it tends to a non-zero constant at vanishing frequency.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

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

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

36 - MATERIALS SCIENCE↗
DOE OSTI2025

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

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

Halide-based solid-state electrolytes↗
DOE OSTI2021

Distinct Mechanisms and Hydricities of Cp*Ir-Based CO 2 Hydrogenation Catalysts in Basic Water

Transition metal-catalyzed reversible hydrogenation of CO 2 to formate in aqueous solutions under ambient conditions is an attractive environmentally friendly strategy for the storage and transportation of H 2 in a liquid chemical form. Here, mechanistic details for the CO 2 hydrogenation by a series of pentamethylcyclopentadienyl iridium(III) (Cp*Ir) complexes with picolinamidate ligands are investigated through density functional theory and kinetic isotope effects (KIEs) calculations. Additionally, adopting a speciation approach at a pH of 8.2, CO 2 hydrogenation mechanisms of 37 distinct protonation states and conformers of Cp*Ir type complexes are investigated. This imposes a challenging test for the computational modeling in terms of providing a consistent correlation with experimental kinetics and KIE studies of multiple catalysts instead of a single complex. Overall, H 2 heterolysis to generate an iridium hydride intermediate was demonstrated to be the rate determining step, which proceeds via distinct pathways depending on the stereoelectronic properties of the ligand. A peculiar iridium dihydride route was also uncovered, which could be optimized to accelerate the H 2 heterolysis. We have further computed the thermodynamic and kinetic hydricities of 74 complexes with [Cp*Ir(L)(H)] q ( q = -1, 0, 1 depending on the ligand charge by deprotonation or protonation events) general formula. Hydricity values do not show any noticeable correlation with the thermodynamics (Δ G ) of [Cp*Ir(L)(HCO 2 – )] formation but exhibit linear correlation with the kinetics (Δ G ‡ ) of electrophilic CO 2 attack to iridium hydride species, especially when the charges and local structural effects of the metal hydrides are considered. Insights gained in this study on (i) the factors determining the preferred pathway for the rate-limiting H 2 heterolysis step, (ii) the correlation between thermodynamic hydricity and the kinetics of CO 2 insertion to iridium hydride species, and (iii) mechanistic analysis via combination of free energy profiles and KIE studies for a series of iridium catalysts will provide guidelines for the design of next-generation catalysts for the reversible H 2 storage through CO 2 utilization.

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

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

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

MXene↗
DOE OSTI2026

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

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

Catalysts↗
DOE OSTI2020

A new magnesium sheet alloy with high tensile properties and room-temperature formability

Lightweight sheet alloys with superior mechanical performance such as high strength, ductility and formability at room temperature (RT) are desirable for high volume automotive applications. However, ductility or formability of metallic alloys at RT are generally inversely related to strength, thereby making it difficult to optimize all three simultaneously. Here we design a new magnesium sheet alloy-ZAXME11100 (Mg-1.0Zn-1.0Al-0.5Ca-0.4Mn-0.2Ce, wt. pct.) via CALPHAD (CALculation of PHAse Diagram) modeling and experimental validation. This new sheet alloy offers an excellent RT formability with a high Index Erichsen (I.E.) value of 7.8 mm in a solution-treated condition (T4), due to its weak and split basal texture and fine grain structure. The new ZAXME 11100 alloy also shows a rapid age-hardening response during post-forming artificial aging treatment at 210 °C for 1hour (T6), resulting in a significant increase of yield strength from 159 MPa (T4) to 270 MPa (T6). The excellent combination of T4 ductility (31%), T4 formability (7.8 mm) and T6 yield strength (270 MPa) in this new magnesium alloy is comparable to that of common 6xxx series aluminum sheet alloys. Thus, this new magnesium sheet alloy is highly attractive for sheet applications in automotive and other industries.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Nonequilibrium thermodynamics of erasure with superconducting flux logic

We implement a thermal-fluctuation-driven logical bit reset on a superconducting flux logic cell. We show that the logical state of the system can be continuously monitored with only a small perturbation to the thermally activated dynamics at 500 mK. We use the trajectory information to derive a single-shot estimate of the work performed on the system per logical cycle. We acquire a sample of 10 5 erasure trajectories per protocol and show that the work histograms agree with both microscopic theory and global fluctuation theorems. The results demonstrate how to design and diagnose complex, high-speed, and thermodynamically efficient computing using superconducting technology.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗