Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Thermodynamics”

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

A thermodynamically consistent discretization of 1D thermal-fluid models using their metriplectic 4-bracket structure

Thermodynamically consistent models in continuum physics, i.e. models which satisfy the first and second laws of thermodynamics, may be expressed using the metriplectic formalism. In this work, we leverage the structures underlying this modeling formalism to preserve thermodynamic consistency in discretizations of a fluid model. The procedure relies (1) on ensuring that the spatial semi-discretization retains certain symmetries and degeneracies of the Poisson and metriplectic 4-brackets, and (2) on the use of an appropriate energy conserving time-stepping method. Here, the minimally simple yet nontrivial example of a one-dimensional thermal-fluid model is treated. It is found that preservation of the requisite symmetries and degeneracies of the 4-bracket is relatively simple to ensure in Galerkin spatial discretizations, suggesting a path forward for thermodynamically consistent discretizations of more complex fluid models using more specialized Galerkin methods.

Hamiltonian structure↗
DOE OSTI2020

Critical evaluation and thermodynamic modeling of the Pd–Sn system

The single crystalline material PdSn 4 , a homologue of the Dirac nodal arc semimetal PtSn 4 , is a promising candidate to search for new topological states with fascinating quantum physical properties. In the present work, the thermochemistry and phase diagram information of the Pd–Sn system published in the literature was collected and critically reviewed, two thermodynamic reassessments of the Pd–Sn system were carried out in the frame of the CALPHAD approach to optimize the Gibbs free energy of each phase presented in the system. The liquid phase was described using respectively the Bragg-Williams random mixing model and the associate solution model. All the reliable thermodynamic property and phase equilibria data can be generally described using the presently obtained thermodynamic descriptions of the Pd–Sn system. Issues related to the further improvement of the thermodynamic descriptions of the Pd–Sn system are discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Complementing a continuum thermodynamic approach to constitutive modeling with symbolic regression

Here, an interpretable machine learning method, physics-informed genetic programming-based symbolic regression (P-GPSR), is integrated into a continuum thermodynamic approach to developing constitutive models. The proposed strategy for combining a thermodynamic analysis with P-GPSR is demonstrated by generating a yield function for an idealized material with voids, i.e., the Gurson yield function. First, a thermodynamic-based analysis is used to derive model requirements that are exploited in a custom P-GPSR implementation as fitness criteria or are strongly enforced in the solution. The P-GPSR implementation improved accuracy, generalizability, and training time compared to the same GPSR code without physics-informed fitness criteria. The yield function generated through the P-GPSR framework is in the form of a composite function that describes a class of materials and is characteristically more interpretable than GPSR-derived equations. The physical significance of the input functions learned by P-GPSR within the composite function is acquired from the thermodynamic analysis. Fundamental explanations of why the implemented P-GPSR capabilities improve results over a conventional GPSR algorithm are provided.

42 ENGINEERING↗
DOE OSTI2024

Thermodynamically consistent Cahn–Hilliard–Navier–Stokes equations using the metriplectic dynamics formalism

Cahn–Hilliard–Navier–Stokes (CHNS) systems describe flows with two-phases, e.g., a liquid with bubbles. Obtaining constitutive relations for general dissipative processes for such systems, which are thermodynamically consistent, can be a challenge. We show how the metriplectic 4-bracket formalism (Morrison and Updike, 2024) achieves this in a straightforward, in fact algorithmic, manner. First, from the noncanonical Hamiltonian formulation for the ideal part of a CHNS system we obtain an appropriate Casimir to serve as the entropy in the metriplectic formalism that describes the dissipation (e.g. viscosity, heat conductivity and diffusion effects). General thermodynamics with the concentration variable and its thermodynamics conjugate, the chemical potential, are included. Having expressions for the Hamiltonian (energy), entropy, and Poisson bracket, we describe a procedure for obtaining a metriplectic 4-bracket that describes thermodynamically consistent dissipative effects. The 4-bracket formalism leads naturally to a general CHNS system that allows for anisotropic surface energy effects. Furthermore, this general CHNS system reduces to cases in the literature, to which we can compare.

Cahn–Hilliard↗
DOE OSTI2022

Overcoming significant challenges in extracting off-stoichiometric thermodynamics using the compound energy formalism through complementary use of experimental and first principles data: A case study of Ba 1-x Sr x FeO 3-δ

The compound energy formalism (CEF) is a powerful framework to describe the thermodynamics of metal oxides as a function of off-stoichiometry, temperature, and composition. The thermodynamic properties are crucial materials design attributes in metal oxide-based oxygen-exchange chemical processes. Despite the richness of information an accurate CEF model provides, a method to determine a unique and accurate fit for oxygen exchange materials remains elusive. This contribution details a method for fitting the CEF model that overcomes the current fitting challenges through three innovations: 1) the combination of density functional theory calculations with experimental data decorrelates excess terms and delineates the enthalpic/entropic contributions to the Gibbs free energy; 2) a systematic determination of the important CEF model terms, removing thermodynamic predetermining human intervention; 3) a self-consistent solution of the starting oxygen offstoichiometry (δ 0 ) of thermogravimetric measurements. Thus, our method enables the reliable extraction of off-stoichiometric metal oxide thermodynamic properties and facilitates rapid materials compositional screening, and reliable process design of systems dependent on off-stoichiometric redox-active metal oxides. We apply this method to a Ba x Sr 1-x FeO 3-δ test case. We find by systematically examining the performance of the CEF model fit with and without each innovation that all three innovations are necessary for an accurate fit. We determined that reduction enthalpy is higher and more sensitive to off-stoichiometry when the Sr fraction is large (139.5 and 185.3 kJ/mol O 2 for SrFeO 3 at δ = 0 and δ = 0.5, respectively vs. nearly constant 83 kJ/mol O 2 for BaFeO 3 ). However, the reduction entropy is mostly insensitive to Sr fraction, but highly dependent on δ suggesting larger contributions in the non-configurational entropy parameters.

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

Experimental Investigation of Thermodynamic Stabilization in Boron Imidazolate Frameworks (BIFs) Synthesized by Mechanochemistry

This study experimentally explores the energetics for the formation of boron–imidazolate frameworks (BIFs), which are synthesized by mechanochemistry. The topologically similar frameworks employ the same tetratopic linker based on tetrakis(imidazolyl)boric acid but differ in the monovalent cation metal nodes. This permits assessment of the stabilizing effect of metal nodes in frameworks with sodalite (SOD) and diamondoid (dia) topologies. The enthalpy of formation from endmembers (metal oxide and linker), which define thermodynamic stability of the structures, has been determined by use of acid solution calorimetry. The results show that heavier metal atoms in the node promote greater energetic stabilization of denser structures. Overall, in BIFs the relation between cation descriptors (ionic radius and electronegativity) and thermodynamic stability depends on framework topology. Thermodynamic stability increases with the metallic character of the cation employed as the metal node, independent of the framework topology. In conclusion, the results suggest unifying aspects for thermodynamic stabilization across MOF systems.

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

Proton, Electron, and Hydrogen-Atom Transfer Thermodynamics of the Metal–Organic Framework, Ti-MIL-125, Are Intrinsically Correlated to the Structural Disorder

Interfacial charge transfer reactions involving protons and/or electrons are fundamental to heterogeneous catalysis and many other reactions relevant to energy, chemical, and biological sectors. Metal–organic frameworks (MOFs) with redox-active metal-oxo nodes have emerged as candidate materials to examine these reactions with near-atomic-level precision, given their crystalline nature. Here, we employed a colloidally stable, Ti-based MOF, Ti-MIL-125, with different crystal sizes to examine catalytically relevant charge transfer thermodynamics. The Ti 8 (μ 2 -O) 8 (μ 2 -OH) 4 nodes structurally mimic TiO 2 , which has shown some PCET reactivity toward reactions of H 2 , O 2 , and others. In this report, we have demonstrated that a change in crystal size induces different amounts of structural disorder to the Ti-oxo node, further changing the thermodynamics of proton/electron/hydrogen-atom transfer reactions. Using electrochemical open-circuit potential (E OCP ) measurements, we have determined that all crystallites undergo a 1H + /1e – redox reaction, which, given the stoichiometry, can be considered as a net H atom transfer (HAT) reaction. The thermodynamics of this HAT reaction, the Ti 3+ O–H bond dissociation free energy (BDFE), was dependent on the crystal size of the MOF, as the decrease in crystal size induced more structural disorder. Our computational calculations have indicated that this difference in BDFE is due to a local change in the geometry of Ti cations, rather than the commonly invoked defects, such as the “missing-linker” defect sites. Individual proton/electron transfer (PT/ET) thermodynamics were also highly dependent on the crystal sizes. These were probed using pK a or band gaps (E g ), respectively. These findings suggest that, particularly when MOFs are nanosized with a large amount of structural disorder, they should no longer be considered “true” single-site catalysts; this is an implicit, but widespread assumption within the MOF-based catalysis field. Implications of these findings will be contrasted with structurally similar metal oxides like TiO 2 and other redox-active MOFs.

Bond dissociation free energy↗
DOE OSTI2024

Pyrophosphate-free glycolysis in Clostridium thermocellum increases both thermodynamic driving force and ethanol titers

Background: Clostridium thermocellum is a promising candidate for production of cellulosic biofuels, however, its final product titer is too low for commercial application, and this may be due to thermodynamic limitations in glycolysis. Previous studies in this organism have revealed a metabolic bottleneck at the phosphofructokinase (PFK) reaction in glycolysis. In the wild-type organism, this reaction uses pyrophosphate (PP i ) as an energy cofactor, which is thermodynamically less favorable compared to reactions that use ATP as a cofactor. Previously we showed that replacing the PP i -linked PFK reaction with an ATP-linked reaction increased the thermodynamic driving force of glycolysis, but only had a local effect on intracellular metabolite concentrations, and did not affect final ethanol titer. Results: In this study, we substituted PP i -pfk with ATP-pfk, deleted the other PPi-requiring glycolytic gene pyruvate:phosphate dikinase (ppdk), and expressed a soluble pyrophosphatase (PPase) and pyruvate kinase (pyk) genes to engineer PP i -free glycolysis in C. thermocellum. We demonstrated a decrease in the reversibility of the PFK reaction, higher levels of lower glycolysis metabolites, and an increase in ethanol titer by an average of 38% (from 15.1 to 21.0 g/L) by using PP i -free glycolysis. Conclusions: By engineering PP i -free glycolysis in C. thermocellum, we achieved an increase in ethanol production. These results demonstrate that optimizing the thermodynamic landscape through metabolic engineering can enhance product titers. While further increases in ethanol titers are necessary for commercial application, this work represents a significant step toward engineering glycolysis in C. thermocellum to increase ethanol titers.

09 BIOMASS FUELS↗
DOE OSTI2025

Thermodynamic Cloud Phase Classifications Using Machine Learning at NSA and ANX

Vertically resolved thermodynamic cloud phase classifications are essential for studies of atmospheric cloud and precipitation processes. The Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) THERMOCLDPHASE Value-Added Product (VAP) uses a multi-sensor approach to classify thermodynamic cloud phase by combining lidar backscatter and depolarization, radar reflectivity, Doppler velocity, spectral width, microwave radiometer-derived liquid water path, and radiosonde temperature measurements. The measured voxels are classified as ice, snow, mixed-phase, liquid (cloud water), drizzle, rain, and liq_driz (liquid+drizzle). We use this product as the ground truth to train three machine learning (ML) models to predict the thermodynamic cloud phase from multi-sensor remote sensing measurements taken at the ARM North Slope of Alaska (NSA) observatory: a random forest (RF), a multilayer perceptron (MLP), and a convolutional neural network (CNN) with a U-Net architecture. Evaluations against the outputs of the THERMOCLDPHASE VAP with one year of data show that the CNN outperforms the other two models, achieving the highest test accuracy, F1-score, and mean Intersection over Union (IOU). Analysis of ML confidence scores shows ice, rain, and snow have higher confidence scores, followed by liquid, while mixed, drizzle, and liq_driz have lower scores. Feature importance analysis reveals that the mean Doppler velocity and vertically resolved temperature are the most influential datastreams for ML thermodynamic cloud phase predictions. The ML models’ generalization capacity is further evaluated by applying them at another Arctic ARM site in Norway using data taken during the ARM Cold-Air Outbreaks in the Marine Boundary Layer Experiment (COMBLE) field campaign. Finally, we evaluate the ML models’ response to simulated instrument outages and signal degradation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Synthesis, characterization, and thermodynamic study of selected K-based zeolites

Potassium-rich zeolites often occur in cementitious systems, as K{sup +} is widespread in various cementitious materials, such as Portland, blended and alkali-activated cements. The knowledge of their stability and of thermodynamic models for solid solutions with Na{sup +} and Ca{sup 2+} are critical to understand long-term development and durability in such cements. Completing previous studies on Na- and Ca-based zeolites, the current work aims to determine the thermodynamic data of 14 types of K-based zeolites, which could possibly form in cementitious systems. The zeolites were synthesized hydrothermally, exchanged with K{sup +}, and characterized thoroughly with respect to framework structures, elemental compositions, water contents, and bond variations. Their thermodynamic properties were derived from the experimental solubility data, which allowed establishing predominance diagrams in the K{sub 2}O-SiO{sub 2}-Al{sub 2}O{sub 3}-H{sub 2}O system. The K-based zeolites typically showed the lowest solubility between 0 and 100 °C, with the notable exception of Ca-gismondine and two Na-based zeolites: natrolite and Na-mordenite.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Determining all thermodynamic transport coefficients for an interacting large N quantum field theory

Thermodynamic transport coefficients can be calculated directly from quantum field theory without requiring analytic continuation to real time. We determine all second-order thermodynamic transport coefficients for the uncharged N-component massless (critical) scalar field theory with quartic interaction in the large N limit, for any value of the coupling. We find that in the large N limit, all thermodynamic transport coefficients for the interacting theory can be expressed analytically in terms of the in-medium mass and sums over modified Bessel functions. We expect our technique to allow a similar determination of all thermodynamic transport coefficients for all theories that are solvable in the large N limit, including certain gauge theories.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Thermodynamic theory of crystal plasticity: Formulation and application to polycrystal fcc copper

We present a thermodynamic description of crystal plasticity. Our formulation is based on the Langer-Bouchbinder-Lookman thermodynamic dislocation theory (TDT), which asserts the fundamental importance of an effective temperature that describes the state of configurational disorder and therefore the dislocation density of the crystalline material. We extend the TDT description from isotropic plasticity to crystal plasticity with many slip systems. Finite-element simulations show favourable comparison with experiments on polycrystal fcc copper under uniaxial compression, tension, and simple shear. Here, the thermodynamic theory of crystal plasticity thus provides a thermodynamically consistent and physically rigorous description of dislocation motion in crystals. We also discuss new insights about the interaction of dislocations belonging to different slip systems.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Coupling computational thermodynamics with density-function-theory based calculations to design L12 precipitates in Fe Ni based alloys

Achieving a high-volume fraction of thermodynamically stable L12-type precipitates that are resistant to coarsening is of great importance for the development of low-cost FeNi based austenitic steels. With the aid of computational thermodynamics, this work designed two model alloys: Fe-37.4Ni-6.1Al-2.9Ti (FNAT) and Fe-45.2Ni-5.9Al-8.5Si (FNAS). Both alloys were designed to contain a similar amount of L1 2 precipitate in Fe-Ni matrix without forming other precipitates. Density-Function-Theory (DFT) calculation was coupled with computational thermodynamics to predict the critical radius at which the precipitates change shape from spherical to cuboidal. The calculation results suggest that critical radius for the FNAT alloy is about two orders of magnitude larger than that for the FeNiAlSi alloy. Phase stability and morphology of the L1 2 precipitates in these two alloys were experimentally investigated through X-ray diffraction, atom probe tomography, and scanning and transmission electron microscopy. The L1 2 precipitates in the FeNiAlSi system were found to be cuboidal and rod shaped, with much larger size than the spherical ones in the FeNiAlTi system, agreeing with the calculation results. This work suggested that coupling computational thermodynamics with DFT calculations can be reliably used to design L1 2 precipitates in FeNi based alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Self-Assembly of Metallo-Supramolecules under Kinetic or Thermodynamic Control: Characterization of Positional Isomers Using Scanning Tunneling Spectroscopy

Coordination-driven self-assembly has been extensively employed to construct a variety of discrete structures as a bottom-up strategy. However, mechanistic understanding regarding whether self-assembly is under kinetic or thermodynamic control is less explored. To date, such mechanistic investigation has been limited to distinct, assembled structures. It still remains a formidable challenge to study the kinetic and thermodynamic behavior of self-assembly systems with multiple assembled isomers due to the lack of characterization methods. Herein, we use a stepwise strategy which combined self-recognition and self-assembly processes to construct giant metallo-supramolecules with 8 positional isomers in solution. With the help of ultrahigh-vacuum, low-temperature scanning tunneling microscopy and scanning tunneling spectroscopy, we were able to unambiguously differentiate 14 isomers on the substrate which correspond to 8 isomers in solution. Through measurement of 162 structures, the experimental probability of each isomer was obtained and compared with the theoretical probability. Finally, such a comparison along with density functional theory (DFT) calculation suggested that although both kinetic and thermodynamic control existed in this self-assembly, the increased experimental probabilities of isomers compared to theoretical probabilities should be attributed to thermodynamic control.

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

Projecting Permafrost Thaw of Sub-Arctic Tundra With a Thermodynamic Model Calibrated to Site Measurements

Northern circumpolar permafrost thaw affects global carbon cycling, as large amounts of stored soil carbon becomes accessible to microbial breakdown under a warming climate. The magnitude of carbon release is linked to the extent of permafrost thaw, which is locally variable and controlled by soil thermodynamics. Soil thermodynamic properties, such as thermal diffusivity, govern the reactivity of the soil-atmosphere thermal gradient, and are controlled by soil composition and drainage. In order to project permafrost thaw for an Alaskan tundra experimental site, we used seven years of site data to calibrate a soil thermodynamic model using a data assimilation technique. The model reproduced seasonal and interannual temperature dynamics for shallow (5–40 cm) and deep soil layers (2–4 m), and simulations of seasonal thaw depth closely matched observed data. The model was then used to project permafrost thaw at the site to the year 2100 using climate forcing data for three future climate scenarios (RCP 4.5, 6.0, and 8.5). Minimal permafrost thawing occurred until mean annual air temperatures rose above the freezing point, after which we measured over a 1 m increase in thaw depth for every 1 °C rise in mean annual air temperature. Under no projected warming scenario was permafrost remaining in the upper 3 m of soil by 2100. We demonstrated an effective data assimilation method that optimizes parameterization of a soil thermodynamic model. The sensitivity of local permafrost to climate warming illustrates the vulnerability of sub-Arctic tundra ecosystems to significant and rapid soil thawing.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Impact of a New Sea Ice Thermodynamic Formulation in the CESM2 Sea Ice Component

Abstract The sea ice component of the Community Earth System Model version 2 (CESM2) contains new “mushy‐layer” physics that simulates prognostic salinity in the sea ice, with consequent modifications to sea ice thermodynamics and the treatment of melt ponds. The changes to the sea ice model and their influence on coupled model simulations are described here. Two simulations were performed to assess the changes in the vertical thermodynamics formulation with prognostic salinity compared to a constant salinity profile. Inclusion of the mushy layer thermodynamics of Turner et al. (2013, https://doi.org/10.1002/jgrc.20171 ) in a fully coupled Earth system model produces thicker and more extensive sea ice in the Arctic, with relatively unchanged sea ice in the Antarctic compared to simulations using a constant salinity profile. While this is consistent with the findings of uncoupled ice‐ocean model studies, the role of the frazil and congelation growth is more important in fully coupled simulations. Melt pond drainage is also an important contribution to simulated ice thickness differences as also found in the uncoupled simulations of Turner and Hunke (2015; https://doi.org/10.1002/2014JC010358 ). However, it is an interaction of the ponds and the snow fraction that impacts the surface albedo and hence the top melt. The changes in the thermodynamics and resulting ice state modify the ice‐ocean‐atmosphere fluxes with impacts on the atmosphere and ocean states, particularly temperature.

Bailey, David A.↗
DOE OSTI2020

Thermodynamics, dynamics, and structure of supercritical water at extreme conditions

Molecular dynamics (MD) simulations to understand the thermodynamic, dynamic, and structural changes in supercritical water across the Frenkel line and the melting line have been performed. The two-phase thermodynamic model [J. Phys. Chem. B(2010), 114(24), 8191-8198] and the velocity autocorrelation functions are used to locate the Frenkel line and to calculate the thermodynamic and dynamic properties. The Frenkel lines obtained from the two-phase thermodynamic model and the velocity autocorrelation criterion do not agree with each other. Structural characteristics and the translational diffusion dynamics of water suggest that this inconsistency could arise from the two oscillatory modes in water, which are associated with the bending of hydrogen bonds and intermolecular collisions inside the first coordination shell. Finally, the overall results lead us to conclude that the universality of the Frenkel line as a dynamic crossover line from rigid to nonrigid fluids is preserved in water.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

Thermodynamic guiding principles of high-capacity phase transformation materials for splitting H 2 O and CO 2 by thermochemical looping

Here, thermochemical looping splitting of water and carbon dioxide (CO 2 ) with greenhouse-gas-free (GHG-free) energy has the potential to help address the Gt-scale GHG emissions challenge. Reaction thermodynamics largely contributes to the main bottlenecks of cost reduction for thermochemical looping water/CO 2 splitting cycle. Here, we analyze thermodynamic driving forces in such cycles with two-phase ternary ferrites as model systems. We find that cation configurational entropy chiefly determines the change of partial molar entropy with oxygen stoichiometry. In addition, our phase diagram analysis accurately predicts the optimal Fe ratio for maximal water/CO 2 splitting capacity in thermal reduction and in chemical reduction based cycles, underlining the significance of phase boundary positions. With chemical reduction, >10% CO 2 conversion and high oxygen exchange capacity can both be achieved. Furthermore, our reduced Gibbs free energy model illustrates critical thermodynamic factors that influence the water/CO 2 splitting capacity. Our research reveals the thermodynamic driving forces underlying the unconventional high-capacity Fe-poor ferrites, further explained via phase diagrams of Fe–Co–O, Fe–Ni–O and Fe–Mg–O. Future materials improvements can be guided by our reduced Gibbs free energy model.

08 HYDROGEN↗