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 667 records · Page 37
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Data-Driven Approach for Controlled Icosahedral Boron- Rich Compound Growth

This final technical report summarizes the research accomplishments and research highlights at the end of the funding period. This project aimed to leverage existing and new computational data produced from first-principles and molecular dynamics simulations to understand the thermodynamic, mechanical, and electronic properties of icosahedral boron compounds. The goal is to achieve targeted material properties by controlling the synthesis routes of these boron-rich compounds.

36 MATERIALS SCIENCE↗
NASA NTRS1991

Theoretical and computational dynamics of a compressible flow

An introduction to the theoretical and computational fluid dynamics of a compressible fluid is presented. The general topics addressed include: thermodynamics and physical properties of compressible fluids; 1D flow of an inviscid compressible fluid; shock waves; fundamental equations of the dynamics of a compressible inviscid non-heat-conducting and radiating fluid, method of small perturbations, linearized theory; 2D subsonic steady potential flow; hodograph and rheograph methods, exact solutions of 2D insentropic steady flow equations, 2D steady transonic and hypersonic flows, method of characteristics, linearized theory of 3D potential flow, nonlinear theory of 3D compressibe flow, anisentropic (rotational) flow of inviscid compressible fluid, electromagnetogasdynamics, multiphase flows, flows of a compressible fluid with transport phenomena.

Pai, Shih-I↗
NASA NTRS2008

Multispectral Emission of the Sun during the First Whole Sun Month: Magnetohydrodynamic Simulations

We demonstrate that a three-dimensional magnetohydrodynamic (MHD) simulation of the corona can model its global plasma density and temperature structure with sufficient accuracy to reproduce many of the multispectral properties of the corona observed in extreme ultraviolet (EW) and X-ray emission. The key ingredient to this new type of global MHD model is the inclusion of energy transport processes (coronal heating, anisotropic thermal conduction, and radiative losses) in the energy equation. The calculation of these processes has previously been confined to one-dimensional loop models, idealized two-dimensional computations, and three-dimensional active region models. We refer to this as the thermodynamic MHD model, and we apply it to the time period of Carrington rotation 1913 (1996 August 22 to September 18). The form of the coronal heating term strongly affects the plasma density and temperature of the solutions. We perform our calculation for three different empirical heating models: (1) a heating function exponentially decreasing in radius; (2) the model of Schrijver et al.; and (3) a model reproducing the heating properties of the quiet Sun and active regions. We produce synthetic emission images from the density and temperature calculated with these three heating functions and quantitatively compare them with observations from E W Imaging Telescope on the Solar and Heliospheric Observatory and the soft X-ray telescope on Yohkoh. Although none of the heating models provide a perfect match, heating models 2 and 3 provide a reasonable match to the observations.

Lionello, Roberto↗
DOE OSTI2024

Defect Thermodynamics and Transport Properties of Perovskite and Fluorite Materials for Solid-Oxide and Proton Conducting Oxide Cells Evaluated Based on Density Functional Theory Modeling

Density functional theory based defect thermodynamic modeling was performed to determine the effect of humidity and H2/O2 gas pressure on various defect chemistry and transport properties of perovskite and fluorite oxides for solid-oxide and proton-conducting-oxide cell applications, with inclusion of the electronic-conducting oxides (as electrodes) and insulating oxides (as electrolytes). Automatic defect generation workflow and first-principles charged defect analysis were implemented on NETL Joule supercomputer for modeling defect equilibria and transport properties of insulating oxides as electrolytes in SOCs and proton-conducting ceramic cells. A GNU Octave defect model subroutines were developed to facilitate defect modeling of electronic conducting oxides in a wide range of operating conditions guided by modeling and experiments. The model includes the hydride defect formation reaction under reducing conditions and allows to incorporate nonstoichiometry effects on the defect thermodynamic parameters. The developed model serves as a platform to facilitate fundamental understanding of the defect thermodynamics in SOC oxide materials and can be used as a novel tool in computational materials screening for SOC and other energy applications.

Lee, Yueh-Lin↗
DOE OSTI2024

Comparative study of machine learning techniques for post-combustion carbon capture systems

Computational analysis of countercurrent flows in packed absorption columns, often used in solvent-based post-combustion carbon capture systems (CCSs), is challenging. Typically, computational fluid dynamics (CFD) approaches are used to simulate the interactions between a solvent, gas, and column's packing geometry while accounting for the thermodynamics, kinetics, heat, and mass transfer effects of the absorption process. These simulations can then be used explain a column's hydrodynamic characteristics and evaluate its CO 2 -capture efficiency. However, these approaches are computationally expensive, making it difficult to evaluate numerous designs and operating conditions to improve efficiency at industrial scales. In this work, we comprehensively explore the application of statistical ML methods, convolutional neural networks (CNNs), and graph neural networks (GNNs) to aid and accelerate the scale-up and design optimization of solvent-based post-combustion CCSs. We apply these methods to CFD datasets of countercurrent flows in absorption columns with structured packings characterized by several geometric parameters. We train models to use these parameters, inlet velocity conditions, and other model-specific representations of the column to estimate key determinants of CO 2 -capture efficiency without having to simulate additional CFD datasets. We also evaluate the impact of different input types on the accuracy and generalizability of each model. We discuss the strengths and limitations of each approach to further elucidate the role of CNNs, GNNs, and other machine learning approaches for CO 2 -capture property prediction and design optimization.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1987

Simplified curve fits for the thermodynamic properties of equilibrium air

New, improved curve fits for the thermodynamic properties of equilibrium air have been developed. The curve fits are for pressure, speed of sound, temperature, entropy, enthalpy, density, and internal energy. These curve fits can be readily incorporated into new or existing computational fluid dynamics codes if real gas effects are desired. The curve fits are constructed from Grabau-type transition functions to model the thermodynamic surfaces in a piecewise manner. The accuracies and continuity of these curve fits are substantially improved over those of previous curve fits. These improvements are due to the incorporation of a small number of additional terms in the approximating polynomials and careful choices of the transition functions. The ranges of validity of the new curve fits are temperatures up to 25 000 K and densities from 10 to the -7 to 10 to the 3d power amagats.

Srinivasan, S.↗
DOE OSTI2025

A phase-field fracture formulation for generalized standard materials: The interplay between thermomechanics and damage

Accurately modeling fracture of ductile materials poses open challenges in the field of computational mechanics due to the multiphysics nature of their failure processes. Integrating the interplay between thermodynamics and damage into ductile fracture models is vital for predicting critical failure modes. Here, in this paper, we develop a versatile phase-field (PF) framework for modeling ductile fracture, taking into account finite-strain elasto-plasticity. The framework stems from a variational formulation of constitutive relations for generalized standard materials (GSMs), whose response is described by a Helmholtz free energy and a dissipation pseudo-potential. Its variational structure is based on a minimum principle for a functional that expresses the sum of power densities for reversible and irreversible processes. By minimizing this functional with a constraint on a von Mises yield function, we derive the evolution equation for the equivalent plastic strain and an associative flow rule. This constrained optimization problem is analytically solved for a wide class of thermo-viscoplasticity models. The key innovations of the current work include (i) a cubic plastic degradation function that accounts for a non-vanishing damage-dependent yield stress, (ii) closed-form expressions of the Helmholtz free energy and dissipation pseudo-potential for three thermo-viscoplasticity models, (iii) an extended Johnson–Cook plasticity model with a nonlinear hardening law, and (iv) a plastic work heat source that depends on the plastic degradation function and a variable Taylor–Quinney (TQ) coefficient. The capabilities of the proposed framework are tested with the aid of four ductile fracture problems, including the Sandia Fracture Challenge. In each of these problems, we examine the evolution of relevant field variables such as the PF order parameter, the equivalent plastic strain, the temperature, and the internal power dissipation density, in addition to the overall structural response quantified by the force–displacement curve. These numerical studies demonstrate that the proposed framework effectively represents ductile fracture, yielding computational results that exhibit good agreement with experimental data.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Mathematical Models and Numerical Methods for High-Fidelity Simulation of Ignition of Reactive Mixtures by Nanosecond Plasma Discharges in Realistic Configurations

We present a newly developed framework for the numerical simulation of ignition of reactive mixtures using single or repeated nanosecond discharge pulses. The framework builds upon the AMReX library, using the existing compressible solver PeleC and low-Mach solver PeleLMeX and allowing for adaptive mesh refinement, complex geometries, and execution on next-generation high-performance computing (HPC) systems. High-fidelity elementary models are adopted for weakly-ionised plasma discharges with significant energy deposition, consistent with nanosecond discharge pulses, and then implemented in the solver. The treatment of non-thermal electrons and charged species, thermodynamics of non-equilbrium species, plasma kinetics, limiting time scales, and boundary conditions for charged species are discussed and addressed for computational efficiency. The framework is demonstrated for three relevant applications: single and multi-pulse discharges in air, single pulse ignition of an ethylene/air mixture, and a three-dimensional plasma discharge in air with temperature stratification. The successful application of the framework demonstrates the feasibility of high-fidelity simulation of ignition of air/hydrocarbon mixtures in three-dimensions with multiple discharge pulses.

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

Real-space observation of ergodicity transitions in artificial spin ice

Abstract Ever since its introduction by Ludwig Boltzmann, the ergodic hypothesis became a cornerstone analytical concept of equilibrium thermodynamics and complex dynamic processes. Examples of its relevance range from modeling decision-making processes in brain science to economic predictions. In condensed matter physics, ergodicity remains a concept largely investigated via theoretical and computational models. Here, we demonstrate the direct real-space observation of ergodicity transitions in a vertex-frustrated artificial spin ice. Using synchrotron-based photoemission electron microscopy we record thermally-driven moment fluctuations as a function of temperature, allowing us to directly observe transitions between ergodicity-breaking dynamics to system freezing, standing in contrast to simple trends observed for the temperature-dependent vertex populations, all while the entropy features arise as a function of temperature. These results highlight how a geometrically frustrated system, with thermodynamics strictly adhering to local ice-rule constraints, runs back-and-forth through periods of ergodicity-breaking dynamics. Ergodicity breaking and the emergence of memory is important for emergent computation, particularly in physical reservoir computing. Our work serves as further evidence of how fundamental laws of thermodynamics can be experimentally explored via real-space imaging.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Ab initio anharmonic thermodynamic properties of cubic Ca Si O 3 perovskite

We present an ab initio study of the thermodynamic properties of cubic CaSiO 3 perovskite (CaPv) over the pressure and temperature range of the Earth's lower mantle. We compute the anharmonic phonon dispersions throughout the Brillouin zone by utilizing the phonon quasiparticle approach, which characterizes the intrinsic temperature dependence of phonon frequencies and, in principle, captures full anharmonicity. Such temperature-dependent phonon dispersions are used to calculate ab initio free energy in the thermodynamic limit (N → ∞) within the framework of the phonon gas model. Accurate free energy calculations enable us to investigate cubic CaPv's thermodynamic properties, e.g., thermal expansivity, Grüneisen parameter, bulk modulus, heat capacity and thermal equation of state, where anharmonic effects are demonstrated. Furthermore, the present methodology provides an important theoretical approach for exploring phase boundaries, thermodynamic, and thermoelastic properties of strongly anharmonic materials at high pressures and temperatures.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Smaller Is Better: The Case for Lower-Order Iodoplumbate Species Dominating MAPbI 3 /Dimethylformamide Solutions

Here, using complementary experimental measurements and computational predictions of spectroscopic measurements (EXAFS, XANES, and UV–vis), we have determined the identity of the most stable iodoplumbate species in dilute lead halide perovskite precursor solutions. We have determined which species are most likely to be thermodynamically stable compared to others that are unstable or metastable. Condensed phase ab initio models were constructed, and the resulting ensembles were used to directly compare the computed signals to the experimental results of the EXAFS, XANES, and UV–vis spectra of PbI 2 :MAI in DMF. The results of this study suggest that only Pb 2+ , PbI + , and PbI 2 are dominant in the dilute lead perovskite precursor solutions as thermodynamically stable entities. Our interpretation of the relative stability of iodoplumbate species in solution, based on an analysis of EXAFS and XANES spectra, provides critically important new insight into the species most likely to be responsible for crystal nucleation and growth in these materials. This insight will have a significant consequence on the broad scientific community and will necessitate the reinterpretation of peaks in the UV–vis spectra of lead halide perovskite precursor solutions.

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

Cloud based tool for analysis of chemical kinetic mechanisms

The size and complexity of chemical kinetic models used to simulate combustion of hydrocarbon fuels continue to increase with our improved understanding of the underlying physical mechanisms, along with the improved ability to utilize such models as the result of increased computational power and efficiency. As mechanisms grow beyond thousands of species and tens of thousands of reactions, it becomes increasingly difficult to manually check for errors and physical inconsistencies. In this work, we present several automated methods to check for such issues in the specification of chemical reaction models. First, we demonstrate how discontinuities in thermodynamic data can cause simulation difficulties manifested as long wall-clock times and convergence failures. To correct this type of problem, we describe an automated method for refitting thermodynamic parameters. We also outline several methods to check the timescales of reaction rate coefficients to ensure physical consistency. All the methods are made available through a web-based tool (https://combustiontools.llnl.gov) intended to aid mechanism developers and users to improve the accuracy and performance of fuel models used by the combustion community.

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

Chemically Separable Ruthenium Hydride Isomers with Distinct Hydride Transfer Properties

Transition metal hydrides are key intermediates in biological and industrial catalysis, where control over hydride donor ability (hydricity) is essential for optimizing reactivity. Herein, we report a series of isomeric ruthenium hydrides in which the ligand trans to the ruthenium hydride bond–either an N-heterocyclic carbene or pyridine–modulates the thermodynamic hydricity by up to 8 kcal mol –1 . These hydrides exhibited distinct reactivity toward CO 2 and various organic hydride acceptors, enabling detailed kinetic and mechanistic analysis. Stopped-flow experiments, isotopic labeling, and computational studies supported an outer-sphere hydride transfer mechanism. Linear free energy relationships (LFERs) and Marcus relationships were utilized to correlate the changes in thermodynamics and kinetics, revealing that ligands with strong trans effects and trans influence can disrupt conventional scaling relationships, which can be a powerful strategy for promoting new reactivity paradigms in hydride transfer reactions.

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

Unexpected Hydricity of a Bis-Carbene Iridium Complex Provides Insight into Electronic Structure Impacts on Hydride Donor Ability

The hydricity (ΔG° H– ), or hydride donor ability, of a transition metal complex is a thermodynamic parameter which can aid in the design and interpretation of various catalytic reactions that involve hydride transfer as a key step. In an attempt to generate a strong hydride donor, the bis-carbene ligand 3,3′-methylenebis(1-methyl-imidazol-2-ylidene) (“bis-mim”) was installed in an iridium hydride complex, [Cp*Ir(bis-mim)H] + . Experimental and computational studies show that [Cp*Ir(bis-mim)H] + is actually a relatively weak hydride donor, however. To understand why the complex is an unexpectedly weak hydride donor, experimental and computational studies probing the steric and electronic effects on hydricity were conducted. Steric factors had a minimal impact on thermodynamic hydricity but dampened kinetic hydricity. In conclusion, the poor thermodynamic hydride donor ability can be attributed to an electronic structure that results in relatively long Cp*−Ir bonds, an unusually high Ir–H BDFE and pK a values that have an outsized influence on thermochemical cycles for hydricity.

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

Time-reversal of an unknown quantum state

For decades, researchers have sought to understand how the irreversibility of the surrounding world emerges from the seemingly time-symmetric, fundamental laws of physics. Quantum mechanics conjectured a clue that final irreversibility is set by the measurement procedure and that the time-reversal requires complex conjugation of the wave function, which is overly complex to spontaneously appear in nature. Building on this Landau-Wigner conjecture, it became possible to demonstrate that time-reversal is exponentially improbable in a virgin nature and to design an algorithm artificially reversing a time arrow for a given quantum state on the IBM quantum computer. However, the implemented arrow-of-time reversal embraced only the known states initially disentangled from the thermodynamic reservoir. Here we develop a procedure for reversing the temporal evolution of an arbitrary unknown quantum state. This opens the route for general universal algorithms sending temporal evolution of an arbitrary system backward in time.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1991

Nonequilibrium radiation and chemistry models for aerocapture vehicle flowfields, volume 3

The computer programs developed to calculate the shock wave precursor and the method of using them are described. This method calculated the precursor flow field in a nitrogen gas including the effects of emission and absorption of radiation on the energy and composition of gas. The radiative transfer is calculated including the effects of absorption and emission through the line as well as the continuum process in the shock layer and through the continuum processes only in the precursor. The effects of local thermodynamic nonequilibrium in the shock layer and precursor regions are also included in the radiative transfer calculations. Three computer programs utilized by this computational scheme to calculate the precursor flow field solution for a given shock layer flow field are discussed.

Carlson, Leland A.↗
DOE OSTI2022

A Computational Framework to Accelerate the Discovery of Perovskites for Solar Thermochemical Hydrogen Production: Identification of Gd Perovskite Oxide Redox Mediators

A high-throughput computational framework to identify novel multinary perovskite redox mediators is presented, and this framework is applied to discover the Gd-containing perovskite oxide compositions Gd 2 BB'O 6 , GdA'B 2 O 6 , and GdA'BB'O 6 that split water. The computational scheme uses a sequence of empirical approaches to evaluate the stabilities, electronic properties, and oxygen vacancy thermodynamics of these materials, including contributions to the enthalpies and entropies of reduction, ΔH TR and ΔS TR . This scheme uses the machine-learned descriptor τ to identify compositions that are likely stable as perovskites, the bond valence method to estimate the magnitude and phase of BO 6 octahedral tilting and provide accurate initial estimates of perovskite geometries, and density functional theory including magnetic- and defect-sampling to predict STCH-relevant properties. Eighty-three promising STCH candidate perovskite oxides down-selected from 4392 Gd-containing compositions are reported, three of which are referred to experimental collaborators for characterization and exhibit STCH activity. Our results demonstrate that the high-throughput computational scheme described herein—which is used to evaluate Gd-containing compositions but can be applied to any multinary perovskite oxide compositional space(s) of interest—accelerates the discovery of novel STCH active redox mediators with reasonable computational expense.

36 MATERIALS SCIENCE↗