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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Subtle metastability of the layered magnetic topological insulator MnBi 2 Te 4 from weak interactions

Abstract Layered quantum materials can host interesting properties, including magnetic and topological, for which enormous computational predictions have been done. Their thermodynamic stability is much less visited computationally, which however determines the existence of materials and can be used to guide experimental synthesis. MnBi 2 Te 4 is one of such layered quantum materials that was predicted to be an intrinsic antiferromagnetic topological insulator, and later experimentally realized but in a thermodynamically metastable state. Here, using a combined first-principles-based approach that considers lattice, charge, and spin degrees of freedom, we investigate the metastability of MnBi 2 Te 4 by calculating the Helmholtz free energy for the reaction Bi 2 Te 3 + MnTe → MnBi 2 Te 4 . We identify a temperature range (~500–873 K) in which the compound is stable with respect to the competing binary phases, consistent with experimental observation. We validate the predictions by comparing the calculated specific heats contributed from different degrees of freedom with experimental results. Our findings indicate that the degrees of freedom responsible for the van der Waals interaction, lattice vibration, magnetic coupling, and nontrivial band topology in MnBi 2 Te 4 not only enable emergent phenomena but also play a crucial role in determining its thermodynamic stability. This conclusion lays the foundation for the future computational material synthesis of novel layered systems.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Geometric and topological entropies of sphere packing

Here, we present a statistical mechanical description of randomly packed spherical particles, where the average coordination number is treated as a macroscopic thermodynamic variable. The overall packing entropy is shown to have two contributions: geometric, reflecting statistical weights of individual configurations, and topological, which corresponds to the number of topologically distinct states. Both of them are computed in the thermodynamic limit for isostatic and weakly underconstrained packings in 2D and 3D. The theory generalizes concepts of granular and glassy configurational entropies for the case of nonjammed systems. It is directly applicable to sticky colloids and predicts an asymptotic phase behavior of sticky spheres in the limit of strong binding.

36 MATERIALS SCIENCE↗
NASA NTRS1973

Thermodynamic and transport properties of gaseous tetrafluoromethane in chemical equilibrium

Equations and in computer code are presented for the thermodynamic and transport properties of gaseous, undissociated tetrafluoromethane (CF4) in chemical equilibrium. The computer code calculates the thermodynamic and transport properties of CF4 when given any two of five thermodynamic variables (entropy, temperature, volume, pressure, and enthalpy). Equilibrium thermodynamic and transport property data are tabulated and pressure-enthalpy diagrams are presented.

Hunt, J. L.↗
NASA NTRS1976

Correlations for determining thermodynamic properties of hydrogen-helium gas mixtures at temperatures from 7,000 to 35,000 K

Simple relations for determining the enthalpy and temperature of hydrogen-helium gas mixtures were developed for hydrogen volumetric compositions from 1.0 to 0.7. These relations are expressed as a function of pressure and density and are valid for a range of temperatures from 7,000 to 35,000 K and pressures from 0.10 to 3.14 MPa. The proportionality constant and exponents in the correlation equations were determined for each gas composition by applying a linear least squares curve fit to a large number of thermodynamic calculations obtained from a detailed computer code. Although these simple relations yielded thermodynamic properties suitable for many engineering applications, their accuracy was improved significantly by evaluating the proportionality constants at postshock conditions and correlating these values as a function of the gas composition and the product of freestream velocity and shock angle. The resulting equations for the proportionality constants in terms of velocity and gas composition and the corresponding simple realtions for enthalpy and temperature were incorporated into a flow field computational scheme. Comparison was good between the thermodynamic properties determined from these relations and those obtained by using a detailed computer code to determine the properties. Thus, an appreciable savings in computer time was realized with no significant loss in accuracy.

Zoby, E. V.↗
NASA NTRS1986

Thermodynamic Calculations for Complex Chemical Mixtures

General computer program, CECTRP, developed for calculation of thermodynamic properties of complex mixtures with option to calculate transport properties of these mixtures. Free-energy minimization technique used in equilibrium calculation. Rigorous equations used in transport calculations. Program calculates equilibrium compositions and corresponding thermodynamic and transport properties of mixtures. CECTRP accommodates up to 24 reactants, 20 elements, and 600 products, 400 of which are condensed. Written in FORTRAN IV for any large computer system.

Mcbride, B. J.↗
DOE OSTI2022

Nonequilibrium statistical thermodynamics of multicomponent interfaces

Nonequilibrium interfacial thermodynamics has important implications for crucial biological, physical, and industrial-scale transport processes. Here, we discuss a theory of local equilibrium for multiphase multicomponent interfaces that builds upon the “sharp” interface concept first introduced by Gibbs, allowing for a description of nonequilibrium interfacial processes such as those arising in evaporation, condensation, adsorption, etc. By requiring that the thermodynamics be insensitive to the precise location of the dividing surface, one can identify conditions for local equilibrium and develop methods for measuring the values of intensive variables at the interface. We then use extensive, high-precision nonequilibrium molecular dynamics (NEMD) simulations to verify the theory and establish the validity of the local equilibrium hypothesis. In particular, we demonstrate that equilibrium equations of state are also valid out of equilibrium, and can be used to determine interfacial temperature and chemical potential(s) that are consistent with nonequilibrium generalizations of the Clapeyron and Gibbs adsorption equations. We also show, for example, that, far from equilibrium, temperature or chemical potential differences need not be uniform across an interface and may instead exhibit pronounced discontinuities. However, even in these circumstances, we demonstrate that the local equilibrium hypothesis and its implications remain valid. These results provide a thermodynamic foundation and computational tools for studying or revisiting a wide variety of interfacial transport phenomena.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Trajectory Ensemble Methods Provide Single-Molecule Statistics for Quantum Dynamical Systems

The emergence of experiments capable of probing quantum dynamics at the single-molecule level requires the development of new theoretical tools capable of simulating and analyzing these dynamics beyond an ensemble-averaged description. In this article, we present an efficient method for sampling and simulating the dynamics of the individual quantum systems that make up an ensemble and apply it to study the nonequilibrium dynamics of the ubiquitous spin-boson model. We generate an ensemble of single-system trajectories, and we analyze this trajectory ensemble using tools from classical statistical mechanics. Our results demonstrate that the dynamics of quantum coherence is highly heterogeneous at the single-system level due to variations in the initial bath configuration, which significantly affects the transient exchange of coherence between the system and its bath. Here, we observe that single systems tend to retain coherence over time scales longer than that of the ensemble. We also compute a novel thermodynamic entanglement entropy that quantifies a thermodynamic driving force favoring system–bath entanglement.

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

Thermodynamic and kinetic hydricity of transition metal hydrides

The prevalence of transition metal-mediated hydride transfer reactions in chemical synthesis, catalysis, and biology has inspired the development of methods for characterizing the reactivity of transition metal hydride complexes. Thermodynamic hydricity represents the free energy required for heterolytic cleavage of the metal–hydride bond to release a free hydride ion, H - , as determined through equilibrium measurements and thermochemical cycles. Kinetic hydricity represents the rate of hydride transfer from one species to another, as measured through kinetic analysis. This tutorial review describes the common methods for experimental and computational determination of thermodynamic and kinetic hydricity, including advice on best practices and precautions to help avoid pitfalls. The influence of solvation on hydricity is emphasized, including opportunities and challenges arising from comparisons across several different solvents. Connections between thermodynamic and kinetic hydricity are discussed, and opportunities for utilizing these connections to rationally improve catalytic processes involving hydride transfer are highlighted.

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

Computational thermomechanics for crystalline rock. Part II: Chemo-damage-plasticity and healing in strongly anisotropic polycrystals

We present a thermal–mechanical–chemical-phase field model that captures the multi-physical coupling effects of precipitation creeping, crystal plasticity, anisotropic fracture, and crack healing in polycrystalline rock at various temperature and strain-rate regimes. This model is solved via a fast Fourier transfer solver with an operator-split algorithm to update displacement, temperature and phase field, and chemical concentration incrementally. In nuclear waste disposal in salt formation, brine inside the crystal salt may migrate along the grain boundary and cracks due to the gradient of interfacial energy and pressure. This migration has a significant implication on the permeability evolution, creep deformation, and crack healing within rock salt but is difficult to incorporate implicitly via effective medium theories compared with computational homogenization. As such, we introduce a thermodynamic framework and a corresponding computational implementation that explicitly captures the brine diffusion along the grain boundary and crack at the grain scale. Meanwhile, the anisotropic fracture and healing are captured via a high-order phase field that represents the regularized crack region in which a newly derived non-monotonic driving force is used to capture the fracture and healing due to the solution–precipitation. Numerical examples are presented to demonstrate the capacity of the thermodynamic framework to capture the multiphysics material behaviors of rock salt.

42 ENGINEERING↗
DOE OSTI2023

Towards Quantum Computing Phase Diagrams of Gauge Theories with Thermal Pure Quantum States

The phase diagram of strong interactions in nature at finite temperature and chemical potential remains largely theoretically unexplored due to inadequacy of Monte-Carlo–based computational techniques in overcoming a sign problem. Quantum computing offers a sign-problem-free approach, but evaluating thermal expectation values is generally resource intensive on quantum computers. To facilitate thermodynamic studies of gauge theories, we propose a generalization of the thermal-pure-quantum-state formulation of statistical mechanics applied to constrained gauge-theory dynamics, and numerically demonstrate that the phase diagram of a simple low-dimensional gauge theory is robustly determined using this approach, including mapping a chiral phase transition in the model at finite temperature and chemical potential. Quantum algorithms, resource requirements, and algorithmic and hardware error analysis are further discussed to motivate future implementations. Thermal pure quantum states, therefore, may present a suitable candidate for efficient thermal simulations of gauge theories in the era of quantum computing.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Accelerating Thermochemical Equilibrium Calculations for Nuclear Reactor Applications

Thermochemical properties play a key role in modeling and simulation of several key phenomena in nuclear reactors. There has been an increasing interest in incorporating CALPHAD-based formulations in multiphysics simulations including for Molten Salt Reactors where knowledge of phase evolution of the salt and the chemical potentials of various elements are of utmost importance in source term analyses and redox control. However, the size of such simulations is often limited by the high computational cost of full thermodynamic equilibrium calculations. This work discusses the current efforts aimed at accelerating thermochemical equilibrium calculations for multiphysics simulations performed using the open-source finite element / finite volume code Multiphysics Object Oriented Simulation Environment (MOOSE) [1]. While several methods have been proposed for accelerating phase equilibrium calculations [2], most focus on relatively small systems and often rely on a- priori knowledge of the state-space of the system. Nuclear materials, however, are often multi-component systems owing to the evolution of composition under irradiation and an approach based on a-priori mapping of phase diagram is often not enough. This work is aimed at demonstrating an on-the-fly surrogate modeling framework that uses active learning to reduce the number of full equilibrium calculations that must be performed. By combining with efficient coupling approaches, the surrogate framework helps in reducing the computational cost of thermodynamic equilibrium informed multiphysics simulations of nuclear materials. The performance is benchmarked against full coupling with the thermochemistry library Thermochimica [3]. This work uses a machine learning based approach for constructing surrogate models to predict the stable phases in a multicomponent system. The surrogates were constructed using neural networks and Gaussian process classification. In this work, we compare the relative performance of the two methods. We also demonstrate the use of caching previous calculations by interpolating the values from nearest neighbors. References [1] Lindsay, A.D., et al. "2.0 – MOOSE: Enabling massively parallel multiphysics simulation", SoftwareX, 20 (2022): 101202. [2] Roos, W.A. and Zietsman J.H. "Accelerating complex chemical equilibrium calculations – A Review", Calphad, 77 (2022): 102380. [3] Piro, M.H.A., et al. "The thermochemistry library Thermochimica", Computational Materials Science, 67 (2013): 266-272.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Experimental and DFT Calculated IR Spectra of Guests in Zeolites: Acyclic Olefins and Host-Guest Interactions

We performed experimental and periodic density functional theory (DFT) IR spectroscopy to investigate the adsorption of acyclic olefins over both acidic and non-acidic zeolites. Two conjugated polyenes, 2,4-dimethyl-1,3-pentadiene (I) and 2,6-dimethyl-2,4,6-octatriene (II) were studied to probe organic intermediates that can form during methanol conversion, and that can lead to deactivating species known collectively as “coke.” We computed vibrational spectra using zeolite-adsorbed and gas-phase models for both neutral and protonated forms of I and II, and compared these DFT results to diffuse reflectance IR Fourier transform (DRIFT) spectra of zeolite-guest systems. Our experimental and computational results are precise enough to pinpoint that the gauche s-cis conformation of species I predominates during adsorption over de-aluminated zeolite beta. Computed zeolite-adsorbed spectra of the protonated species I and II best represent the DRIFT spectra obtained after the adsorption of the olefins on HMOR at 20 ?, with computed bands at 1543 and 1562 cm–1 for molecules I+ and II+, respectively, attributed to the allylic stretching mode, ?(C=C–C+). These computed band frequencies are within 6 cm–1 of experimental data and confirm that the interaction between neutral acyclic olefins and acidic zeolites leads to protonation of the olefin. Comparison of computed spectra of the protonated species in the gas phase to those in the zeolite indicates that the electrostatic interaction between alkenyl and alkadienyl cations and negative zeolite framework does not significantly impact the position of the allylic stretching bands. These results highlight that computed spectroscopy and thermodynamics coupled with experimental spectra can be used for elucidating complex mixtures in zeolites, and certain spectral features of adsorbed olefins can be accurately modeled by gas-phase calculations. ACKNOWLEDGMENTS BM and SMA were supported by the National Science Foundation under Award # CBET-1512442. Acknowledgment is made to the Donors of the American Chemical Society Petroleum Research Fund for support (or partial support) of this research (#589978-ND5, EDH and FCJ). MDB was supported by the US Department of Energy (DOE), Office of Science, Office of Basic Energy Sciences (BES), Division of Materials Science and Engineering. CJM was supported by the DOE, Office of Science, BES, Division of Chemical Sciences, Geosciences, and Biosciences. Pacific Northwest National Laboratory (PNNL) is a multiprogram national laboratory operated by Battelle for the U.S. DOE. Computing resources were generously provided through PNNL’s institutional computing (PIC).

Manookian, Babgen↗
DOE OSTI2022

Quantification of Dead Lithium on Graphite Anode under Fast Charging Conditions

A series of computational and experimental studies were conducted to understand the onset of lithium plating and subsequent quantification of dead lithium on graphite electrodes in the design of fast charging batteries. The experiments include titration and relaxation studies for detecting initiation of lithium metal plating for various SOC and C-rates, which are compared against the thermodynamically consistent phase field computational results. The collaborative study on "model graphite electrode" with 2.18 mAh cm-2 nominal capacity at 25 degrees C demonstrates: (1) the macroscopic voltage response during relaxation studies indicate the reintercalation of plated lithium into the graphite anode; (2) for SOC below 60% and low C-Rates, there is no dead lithium; (3) for SOC between 60% to 80%, and C-Rates in the range of 4C-6C show dead lithium both in experiments and simulations; (4) at 100% SOC and 4C-6C rates, large amounts of dead lithium are observed. The study presented here allows us to evaluate the effects of the physical properties of the electrochemical system on plating and stripping kinetics and the amount of dead lithium on graphite electrodes, which determines the cell capacity loss under fast charge.

batteries↗
DOE OSTI2023

Combustion of 10–100 μm aluminum droplets in detonation products gases

Here we describe a two-phase model of combustion effects in aluminized high explosive (HE) charges. It is based on: (i) a Gas Dynamic Model of the expansion of the detonation product gases and their turbulent combustion with air; and (ii) a Heterogeneous Continuum Model of aluminum (Al) droplets and their combustion with the detonation product gases. Initial conditions are based on an analytical similarity solution for a cylindrical Chapman-Jouguet (CJ) detonation propagating at the CJ detonation velocity. The CJ jump conditions are computed at the thermodynamic equilibrium state by the Cheetah code, assuming the Al droplets are inert. We assume that the Al is 10 % of the charge mass and occurs as droplets at the CJ state. Different initial droplet diameters, ranging from 10 to 100 microns, are studied. A hydrodynamic combustion model based on large Damköhler numbers is employed in this study, -3σ(1 + 0.276 $\sqrt{Re}$)/($Kd^2_w$). It has a square-root dependence on the Reynolds number (Re) and inversed-squared dependence on the droplet diameter (d w ). The burnout time (t B ) of the Al droplets has a three-halves dependence on the droplet diameter, t B ~ ($d^0_w$) 3/2 . After burnout, the detonation products act as detonation products of the HE charge with active Al. They turbulently mix with air and form a combustion layer on the outer edge of the fireball. Details of the two-phase model, initial conditions and evolution of the flow field will be described.

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

Cantor-derived medium-entropy alloys: bridging the gap between traditional metallic and high-entropy alloys

The year 2004 marked the beginning of a new era in the design of metallic materials, as the concept of multiple principal component alloys, commonly known as High-Entropy Alloys (HEAs), was proposed by Cantor and Yeh. The unexpected single-phase microstructure, instead of the expected brittle intermetallic compounds, was attributed to the large entropy of mixing and immediately caught the attention of the scientific community. Today, HEAs are considered important advanced materials and a broad range of alloys using nominally the same design principle have been investigated. Despite that, the CrMnFeCoNi (Cantor) alloy stands out as the most successful HEA due to its outstanding mechanical properties and microstructure. In this scenario, variants of the Cantor alloy, named medium-entropy alloys (MEAs), are gaining significant interest as they display a better industrial potential than both HEAs and traditional alloys. These variants of the Cantor alloy with only three or four main elements result in 15 possible combinations. The microstructure of these alloys is discussed in terms of advanced characterization as well as thermodynamic parameters and computational simulation. Their phase stability is addressed over a wide range of temperatures and strain rates. The mechanical properties, especially the fracture toughness, of the CrFeCoNi and CrCoNi alloys have been reported to be even superior to those of the Cantor alloy and most modern engineering alloys. This is associated with the formation of a continuous sequence of strengthening mechanisms, including hierarchical twin networks, which serve to prolong the strain hardening. The present article reviews and critically assesses, for the first time, recent advances in these Cantor-derived MEAs.

36 MATERIALS SCIENCE↗