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 757 records · Page 42
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Computational Insights into Phase Equilibria Between Wide-Gap Semiconductors and Contact Materials

Novel wide-band-gap semiconductors are needed for next-generation power electronics, but there is a gap between a promising material and a functional device. Finding stable (metal) contacts is one of the major challenges that is currently dealt with mainly via trial and error. Herein, we computationally investigate the thermochemistry and phase coexistence at the junction between three wide-gap semiconductors, ..beta..-Ga2O3, GeO2, and GaN, and possible contact materials. The pool of possible contacts includes 47 elemental metals and a set of 4 common, n-type transparent conducting oxides (ZnO, TiO2, SnO2, and In2O3). We use first-principles thermodynamics to model the Gibbs free energies of chemical reactions as a function of gas pressure (pO2/pN2) and equilibrium temperature. We deduce whether a semiconductor/contact interface will be stable at relevant conditions or a chemical reaction between them is to be expected, possibly influencing the long-term reliability and performance of devices. We generally find that most elemental metals tend to oxidize or nitridize and form various interface oxide/nitride layers. Exceptions include select late- and post-transition metals and, in the case of GaN, also the alkali metals, which are predicted to exhibit stable coexistence, although in many cases at relatively low gas partial pressures. Similar is true for the transparent conducting oxides, for which, in most cases, we predict a preference toward forming ternary oxides when in contact with ..beta..-Ga2O3 and GeO2. The only exception is SnO2, which we find to form stable contacts with both oxides. Finally, we show how the same approach can be used to predict gas partial pressure vs temperature phase diagrams to help direct synthesis of ternary compounds. We believe these results provide a valuable guidance in selecting contact materials to wide-gap semiconductors and suitable growth conditions.

contact materials↗
DOE OSTI2025

Computational discovery of a novel double transition metal nitride MXene and its applications as an anchoring and catalytic material in Li–Se batteries

Double transition metal (DTM) nitride MXenes offer enhanced electrical conductivity and tunable catalytic properties compared to conventional carbide-based MXenes. In this study, we employed first-principles density functional theory (DFT) calculations to discover and computationally validate a novel DTM nitride MXene, Nb 2 TiN 2 , derived from its MAX phase precursor and investigated its potential as an anchoring material (AM) for Li–Se batteries. This newly proposed MXene expands the compositional landscape of DTM nitrides and opens new avenues for functional material design. We performed a comprehensive analysis of the thermodynamic and electronic properties of Nb 2 TiAlN 2 , and the MAX phase precursor to Nb 2 TiN 2 to assess its structural stability and exfoliation potential. Exfoliation energy calculations confirmed the feasibility of synthesizing Nb 2 TiN 2 from Nb 2 TiAlN 2 . We then explored the functionalized form, Nb 2 TiN 2 S 2 , evaluating its capability to serve as an effective anchoring material (AM) in Li–Se batteries by analyzing the reaction mechanisms and kinetics of the selenium reduction reaction (SeRR). Our results indicate that Nb 2 TiN 2 S 2 exhibits a strong binding affinity for lithium polyselenides (Li 2 Se n ), effectively suppressing the shuttle effect. Gibbs free energy calculations for the rate-limiting step of the SeRR reveal favorable kinetics and reduced reaction barriers. Overall, this study provides a detailed evaluation of the structural and electronic properties of a newly proposed DTM nitride MXene and its S-functionalized derivative and the catalyzing effect of Nb 2 TiN 2 S 2 in accelerating the reaction kinetics in Li–Se batteries. These findings underscore the potential importance of the further exploration of MXenes to address current challenges in high-performance Li–Se batteries.

25 ENERGY STORAGE↗
DOE OSTI2022

Crack-Free 30% Chromium-Nickel Alloy Welding Products for Nuclear Service

We report prior research in the development of 30% chromium-nickel alloy nuclear welding wires has resulted in the resolution of primary water stress corrosion cracking (PWSCC) and ductility dip cracking (DDC) as well as improvement in solidification cracking (SC) resistance. The resolution of DDC exhibits some Laves phase, which has a negative effect on SC resistance. In this study, the use of an alternate carbide former, tantalum (Ta), in combination with niobium (Nb) was researched. Three heats of recently designed Filler Metal 52MSS-Ta (i.e., HV1648, HV1673A, and VX131WXW) were melted, fabricated, and systematically studied. DDC and SC were evaluated with thermodynamic modeling using the Scheil solidification simulation model, two types of varestraint tests, and strain-to-fracture (STF) testing. The varestraint and STF test results showed an improved SC resistance with reduced Laves phase and concurrent excellent DDC resistance. Optimized compositions with low Laves phase also exhibited high threshold strain values (TSVs) in the STF test. VX131WXW — which contains 2.81 wt-% Ta, 0.6 wt-% Nb, and 6 wt-% iron (Fe) - exhibited a TSV of 24%. Thermo-Calc computed the Laves phase to be 0.24% for VX131WXW compared to 0.06% in HV1673A. This difference in Laves phase resulted in the lower SC resistance of VX131WXW compared to HV1673A when measured with longitudinal varestraint testing. The maximum crack distance for HV1673A was about 0.6 mm while that of heat VX131WXW was about 1.0 mm. The typical diluted weld deposit made with VX131WXW was also resistant to PWSCC due to the chromium content exceeding 24%. These simultaneous results mark progress toward crack-free welds and provide direction for further optimization of Ta-containing filler metals.

30% chromium-nickel alloy↗
DOE OSTI2024

Fuel‐Driven π‐Conjugated Superstructures to Form Transient Conductive Hydrogels

Despite advances in creating dissipative materials with transient properties, such as hydrogels and active droplets, their application remains confined to temporal changes in structural properties. Developing out-of-equilibrium materials whose electronic functions are parameterized by a chemical reaction cycle is challenging. Yet, this class of materials is required to construct biomimetic materials. In contrast to traditional chemical reaction cycles that exploit molecularly dissolved building blocks at thermodynamic equilibrium, we show that fiber structures derived from reactive naphthalene diimide (NDI) building blocks can be used as resting states to form far-from-equilibrium conductive hydrogels after the addition of chemical fuels. Upon fueling the NDI-derived fibers, a dual-component activation and deactivation pathway is deduced by kinetic analysis and is absent when using a molecularly dissolved resting state. Investigating the solid-state morphologies of the structures formed throughout the fuel-driven reaction cycle using cryo-EM reveals that the resting thermodynamic fibers evolve to transient thicker fibrils and layered superstructures. We show that the transient redox-active hydrogels exhibit a nearly threefold increase in electrical conductivity upon fuel consumption before reverting to their original value over hours. In conclusion, these far-from-equilibrium materials are potential candidates in applications such as programmable biorobotics and chemical computing.

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

Solvent-mediated contaminant removal from plastic waste using thermodynamic modeling

Plastics recycling is hindered by the compositional complexity of plastic waste, which can include numerous polymer components as well as low concentrations of additives and non-intentionally added substances. These latter small-molecule species, which we collectively refer to as contaminants, can harm human health and will build up in recycled plastic causing environmental and downstream processing challenges if not removed. In this work, we present molecular modeling approaches using the COnductor-like Screening MOdel for Real Solvents (COSMO-RS) to guide the selection of solvents that are capable of removing targeted contaminants from plastic waste. By considering the thermodynamic partitioning of contaminant species between a solvent phase and polymer phase, we identify guidelines for solvent selection to promote either the low-temperature extraction of contaminants from plastic waste or the removal of contaminants as part of a dissolution-based plastics recycling process. We present four case studies to illustrate the application of the computational approach to the removal of brominated flame retardants, phthalates, and selected perfluoroalkyl substances, and compare to both literature and newly collected experimental data to illustrate model prediction accuracy. Furthermore, the case studies highlight the capability of the modeling approach to help design recycling processes that explicitly account for contaminant removal, thereby increasing product purity during dissolution-based recycling or facilitating chemical recycling of contaminant-free plastics.

Zhou, Panzheng [University of, Wisconsin, Madison,↗
DOE OSTI2023

Interaction models and configurational entropies of binary MoTa and the MoNbTaW high entropy alloy

We introduce a simplified method to model the interatomic interactions of high entropy alloys based on a lookup table of cluster energies. Furthermore, these interactions are employed in replica exchange Monte Carlo simulations with histogram analysis to obtain thermodynamic properties across a broad temperature range. Kikuchi's cluster variation method entropy formalism is applied to directly calculate entropy from statistics on short- and long-range chemical order, and we discuss the convergence of the entropy as clusters of differing size are included. A high temperature series expansion aids in our understanding of the convergence. Computer codes implementing these methods, and supporting data, are freely available on the internet.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Finite-temperature many-body perturbation theory for anharmonic vibrations: Recursions, algebraic reduction, second-quantized reduction, diagrammatic rules, linked-diagram theorem, finite-temperature self-consistent field, and general-order algorithm

A unified theory is presented for finite-temperature many-body perturbation expansions of the anharmonic vibrational contributions to thermodynamic functions, i.e., the free energy, internal energy, and entropy. The theory is diagrammatically size-consistent at any order, as ensured by the linked-diagram theorem proved in this study, and, thus, applicable to molecular gases and solids on an equal footing. It is also a basis-set-free formalism, just like its underlying Bose–Einstein theory, capable of summing anharmonic effects over an infinite number of states analytically. It is formulated by the Rayleigh–Schrödinger-style recursions, generating sum-over-states formulas for the perturbation series, which unambiguously converges at the finite-temperature vibrational full-configuration-interaction limits. Two strategies are introduced to reduce these sum-over-states formulas into compact sum-over-modes analytical formulas. One is a purely algebraic method that factorizes each many-mode thermal average into a product of one-mode thermal averages, which are then evaluated by the thermal Born–Huang rules. Canonical forms of these rules are proposed, dramatically expediting the reduction process. The other is finite-temperature normal-ordered second quantization, which is fully developed in this study, including a proof of thermal Wick’s theorem and the derivation of a normal-ordered vibrational Hamiltonian at finite temperature. The latter naturally defines a finite-temperature extension of size-extensive vibrational self-consistent field theory. These reduced formulas can be represented graphically as Feynman diagrams with resolvent lines, which include anomalous and renormalization diagrams. Two order-by-order and one general-order algorithms of computing these perturbation corrections are implemented and applied up to the eighth order. The results show no signs of Kohn–Luttinger-type nonconvergence.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2026

Direct CO 2 Reduction to CO with an Fe 4 S 4 -Based Coordination Polymer

Fe 4 S 4 clusters play essential roles in nature, classically in electron transport but increasingly in newly discovered reactivity or catalysis. These roles have spurred interest in developing synthetic Fe 4 S 4 systems and while several molecular and material systems built from Fe 4 S 4 clusters have been developed, comparatively few examples of synthetic Fe 4 S 4 cluster-based catalysts exist. Herein, we present the use of an Fe 4 S 4 -based coordination polymer as a catalyst for the direct and selective electroreduction of CO 2 to CO. Computational studies suggest that the reaction proceeds through CO 2 binding to a reduced Fe 4 S 4 cluster, followed by a series of protonation, reduction, and H 2 O loss steps to yield a CO-bound cluster that can finally exchange with CO 2 to restart the catalytic cycle. CO bound clusters are predicted to be thermodynamically stable, suggesting that carbonyl species might be off-cycle intermediates. Mechanistic CV studies as well as in situ studies by IR spectroscopy provide evidence for carbonyl-ligated clusters, supporting these compounds as unusual examples of small molecule binding to Fe 4 S 4 clusters. Finally, this work establishes Fe 4 S 4 cluster-based coordination polymers as direct electrocatalysts for CO 2 reduction and provides mechanistic insights into how these species mediate catalytic conversions of small molecules.

cluster chemistry↗
DOE OSTI2020

Prediction of structure and cation ordering in an ordered normal-inverse double spinel

Abstract Spinels represent an important class of technologically relevant materials, used in diverse applications ranging from dielectrics, sensors and energy materials. While solid solutions combining two “single spinels” have been explored in a number of past studies, no ordered “double” spinels have been reported. Based on our first principles computations, here we predict the existence of such a double spinel compound MgAlGaO 4 , formed by an equimolar mixing of MgAl 2 O 4 normal and MgGa 2 O 4 inverse spinels. After studying the details of its atomic and electronic structure, we use a cluster expansion based effective Hamiltonian approach with Monte Carlo simulations to study the thermodynamic behavior and cation distribution as a function of temperature. Our simulations provide strong evidence for short-ranged cation order in the double spinel structure, even at significantly elevated temperatures. Finally, an attempt was made to synthesize the predicted double spinel compound. Energy Dispersive X-ray Spectrometry and X-ray diffraction Rietveld refinements were performed to characterize the single-phase chemical composition and local configurational environments, which showed a favorable agreement with the theoretical predictions. These findings suggest that a much larger number of compounds can potentially be realized within this chemical space, opening new avenues for the design of spinel-structured materials with tailored functionality.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Laser Measurements and Modeling of Shock Tunnel Freestream Velocity and Multispecies Thermal Nonequilibrium

Coherent anti-Stokes Raman scattering (CARS) and nitric oxide molecular tagging velocimetry (NO-MTV) were used to characterize the freestream in Sandia’s Hypersonic Shock Tunnel (HST) using a burst-mode laser operated at 100 kHz. Experiments were performed at nominal freestream velocities of 3 and 4 km/s using both air and N 2 test gas. The CARS diagnostic provides nonequilibrium characterization of the flow by measuring vibrational and rotational temperatures of N 2 and O 2 , which are compared to NO temperatures from separate laser absorption experiments. Simultaneous, collinear freestream velocities were measured using NO MTV along with pitot pressure measurements. Furthermore, this extensive freestream dataset is compared to nonequilibrium CFD that is capable of modeling species-specific, vibrational temperatures throughout the nozzle expansion. Significant nonequilibrium between vibrational and rotational temperatures was measured at each flow condition. N 2 exhibits the most nonequilibrium, followed by O 2 and NO. The CFD model captures this trend, although it consistently overpredicts N 2 and O 2 vibrational temperatures. At 3 km/s, the modeled NO nonequilibrium is underpredicted, whereas it is overpredicted at 4 km/s. Good agreement is seen between CFD and the velocity and rotational temperature measurements. Experiments with water added to the test gas yielded no discernable difference in vibrational relaxation.

Aerodynamics↗
DOE OSTI2022

Tailoring critical Al concentration to form external Al 2 O 3 scale on Ni–Al alloys by computational approach

Abstract Nickel (Ni)‐based superalloys for high‐temperature applications are often designed to form a continuous and slow‐growing oxide scale by adding Al and Cr and other beneficial elements. In the present work, the critical Al concentration in Ni–Al alloys needed to establish an α‐Al 2 O 3 scale in contrast to internal oxide formation is predicted as a function of temperature by means of the CALPHAD approach coupled with models in the literature, which account for the thermodynamics and kinetics of oxidation. The present thermodynamic remodeling of the Ni–O system results in a better agreement with experimental data of oxygen solubility in Ni at high temperatures. The oxygen solubility is combined with kinetic parameters to determine oxygen permeability in Ni, and the critical Al concentration needed to establish an α‐Al 2 O 3 scale at a given exposure temperature. Good agreement is found with available experimental data for both oxygen permeability and critical Al concentration, indicating the capacity of the CALPHAD approach to tailor oxidation resistance for materials of interest using thermodynamic and kinetic knowledge.

Ross, Austin↗
DOE OSTI2026

Data-based filtered dissipation rate modelling for multi-modal turbulent combustion: evaluating a priori model generalizability

Manifold-based models offer a computationally efficient alternative to directly transporting the thermochemical state in computational simulations of turbulent reacting flows, projecting the high-dimensional thermochemical state-space onto a low-dimensional manifold. Recent efforts have yielded a manifold-based model applicable to multi-modal combustion, enabling reconstruction of the thermochemical state from solutions to two-dimensional manifold equations in mixture fraction and generalized progress variable that are parameterised by three scalar dissipation rates. In coarse-grained simulations such as Large Eddy Simulation (LES), closure of the multi-modal manifold equations and subfilter variances/covariance requires closure of three filtered scalar dissipation rates. Here, the present work adopts a data-based approach, providing closure for the three filtered scalar dissipation rates via deep neural networks (DNNs). High-fidelity datasets corresponding to an autoigniting n-dodecane jet flame and a bluff body swirl-stabilized confined lifted spray flame of two aviation fuels (Jet-A and C1) with different ignition propensities are leveraged to generate training data that spans a diverse range of thermodynamic conditions and combustion modes, including low- and high-temperature ignition regimes in addition to premixed and nonpremixed behaviour. A final DNN model is trained to enforce inherent physical constraints by learning nonlinear functional transformations of the three filtered scalar dissipation rates. The generalizability of this constrained DNN model is demonstrated a priori via conditional statistics evaluated on the lifted spray flame with C1–a configuration that had not been included in the training data. Excellent DNN agreement with conditional DNS statistics is observed, and integrated gradients are computed to identify the most sensitive input variables. The similarity of the marginal PDFs of the most informative input variables and outputs across configurations are quantified via the Wasserstein metric, demonstrating that data-based models may successfully generalize to unseen parametric conditions so long as the most informative input variables share similar distributions across training and testing datasets.

Data-based modelling↗
DOE OSTI2020

Multidimensional simulations of Mckenna-driven flow tube configuration: Investigating non-ideality in NO x formation flow tube experiments

Multidimensional simulations have been conducted to simulate atmospheric pressure, flat-flame/McKenna-burner-driven-flow tube experiments targeted to obtain NO x speciation data for predicting/analyzing syngas combustion emissions. In a prior work, we demonstrated the impacts of multidimensional transport on post flame region prediction departures from those assuming unidimensional flow/transport conditions. In this work, we develop and utilize a multidimensional laminar reacting flow solver to simulate the fully coupled flame and post flame regions to further elucidate the impacts of the earlier unidimensional modeling assumptions on interpreting post flame NO x experimental data. The model is used to simulate a lean, premixed syngas/air flame and its associated post flame regions within a cylindrical flow-tube-like arrangement. The combustion process takes place under atmospheric condition with trace amount of NO x seeding fed into the inlet gas stream. The spatial evolution of NO x species (NO and NO 2 ) in the flame and in the post-combustion zone suggests two distinct regions: 1) a region encompassing the flame structure itself; and 2) a post flame region in which the temperature decays due to both axial and radial transport processes. The predictions show that for the conditions studied, a pulsatile flow field exists due to the formation of an expanding and contracting recirculation zone in the outer periphery of the flow tube. By resolving the nature of the flow, the resulting time-averaged temperature and species concentrations show improved agreement with existing experimental measurements. The flow-field interaction results in radial inhomogeneities in the NO 2 profiles with the maximum concentration offset from the flow centerline. The location of the peak in NO 2 is coupled with radial temperature gradients from wall cooling effects and their significant influence on NO/NO 2 interconversion kinetics, producing notable NO 2 accumulation in regions near the wall. Geometrical configurations capable of suppressing/minimizing the pulsatile nature are also investigated and the results are compared. Other experimental configurations could be considered in parametric simulations to determine the optimal configuration that would minimize non-idealities in the observations. The work shows the value in performing such computations in advance of settling on a particular design for flow tube/flow reactor experiments.

42 ENGINEERING↗
DOE OSTI2022

Metal dissolution from first principles: Potential-dependent kinetics and charge transfer

Metal dissolution and its inverse process are integral to both corrosion and electrodeposition; however, many mechanistic details regarding the dissolution process are challenging to decipher. These include how ion dissolution kinetics and charge transfer are influenced by the competition between metal and solvent interactions under an electrode potential. Here, in this work, we introduce a computational framework based on density functional theory with grand-canonical treatment of electrons to directly predict the potential energy landscape for metal dissolution at a constant potential. Using aluminum as an example, we demonstrate that dissolution kinetics is governed by competing kinetics between two physical processes associated with metal–metal bond breaking and ion-migration within the electrical double layer, respectively. We identify a kinetic transition between regimes dominated by each of these processes and show that this transition depends on the operating electrode potential, among other key factors. It is further found that kinetics and thermodynamics of these processes can be described with a simple, one-parameter Marcus-theory-type model. Beyond offering new understanding of charge transfer during dissolution, our simulation protocol provides a recipe for directly predicting other important quantities in electrochemical reactions from first principles that are difficult to measure, such as the symmetry factor.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Thermodynamic assessment of lithium halide reciprocal salt systems for energy applications

Lithium halides are of interest because of their unique applications as fuel or coolants in molten salt reactors (MSRs) and electrolytes in thermally activated batteries. Because iodine is an important fission product for MSRs, the Molten Salt Thermal Properties Database-Thermochemical (MSTDB-TC) is being expanded to include relevant iodide reciprocal salt systems. In this work, we assess the thermochemical properties of lithium halide pseudobinary systems based on reported experimental phase diagrams, eutectic temperatures, and enthalpies of mixing. Within the framework of the modified quasi-chemical model in the quadruplet approximation (MQMQA), extrapolations to the related pseudoternary and pseudoquaternary systems have been performed. The pseudoquaternary representations of phase equilibria, enthalpy, entropy, and heat capacity are represented in the resultant models. Of interest for generating source terms for accident analysis is the ability to compute the vapor pressures of the dominant iodide species for the systems at any composition. Here, an example of such calculations is provided for the LiF-LiI system at x LiI = 0.01 and x LiI = 0.70, in this case indicating that the LiI vapor species are predominant.

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

Colossal Cryogenic Electro‐Optic Response Through Metastability in Strained BaTiO 3 Thin Films

The search for thin film electro-optic materials that can retain superior performance under cryogenic conditions has become critical for quantum computing. Barium titanate thin films show large linear electro-optic coefficients in the tetragonal phase at room temperature, which is severely degraded down to ≈200 pm V −1 in the rhombohedral phase at cryogenic temperatures. There is immense interest in manipulating these phase transformations and retaining superior electro-optic properties down to liquid helium temperature. Utilizing the thermodynamic theory of optical properties, a large low-temperature electro-optic response is designed by engineering the energetic competition between different ferroelectric phases, leading to a low-symmetry monoclinic phase with a massive electro-optic response. The existence of this phase is demonstrated in a strain-tuned BaTiO 3 thin film that exhibits a linear electro-optic coefficient of 2516 ± 100 pm V −1 at 5 K, which is an order of magnitude higher than the best reported performance thus far. Importantly, the electro-optic coefficient increases by 100 × during cooling, unlike the conventional films, where it degrades. Further, at the lowest temperature, significant higher order electro-optic responses also emerge. These results represent a new framework for designing materials with property enhancements by stabilizing highly tunable metastable phases with strain.

barium titanate↗
DOE OSTI2021

The statistical mechanics of near-BPS black holes

Due to the failure of thermodynamics for low temperature near-extremal black holes, it has long been conjectured that a 'thermodynamic mass gap' exists between an extremal black hole and the lightest near-extremal state. For non-supersymmetric near-extremal black holes in Einstein gravity with an AdS2 throat, no such gap was found. Rather, at that energy scale, the spectrum exhibits a continuum of states, up to non-perturbative corrections. In this paper, we compute the partition function of near-BPS black holes in supergravity where the emergent, broken, symmetry is PSU(1, 1|2). To reliably compute this partition function, we show that the gravitational path integral can be reduced to that of a $\mathcal{N}=4$ supersymmetric extension of the Schwarzian theory, which we define and exactly quantize. In contrast to the non-supersymmetric case, we find that black holes in supergravity have a mass gap and a large extremal black hole degeneracy consistent with the Bekenstein–Hawking area. Furthermore, our results verify a plethora of string theory conjectures, concerning the scale of the mass gap and the counting of extremal micro-states.

79 ASTRONOMY AND ASTROPHYSICS↗