Search NASA⌕ Search

SEARCH · Search NASA

Results for “Thermodynamic Modeling”

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

The Physics behind Precipitation Onset Bias in CMIP6 Models: The Pseudo-Entrainment Diagnostic and Trade-Offs between Lapse Rate and Humidity

Conditional instability and the buoyancy of plumes drive moist convection but have a variety of representations in model convective schemes. Vertical thermodynamic structure information from Atmospheric Radiation Measurement (ARM) sites and reanalysis (ERA5), satellite-derived precipitation (TRMM3b42), and diagnostics relevant for plume buoyancy are used to assess climate models. Previous work has shown that CMIP6 models represent moist convective processes more accurately than their CMIP5 counterparts. However, certain biases in convective onset remain pervasive among generations of CMIP modeling efforts. Here we diagnose these biases in a cohort of nine CMIP6 models with sub-daily output, assessing conditional instability in profiles of equivalent potential temperature θ e and saturation equivalent potential temperature θ es in comparison to a plume model with different mixing assumptions. Most models capture qualitative aspects of the θ es vertical structure, including a substantial decrease with height in the lower free troposphere associated with the entrainment of subsaturated air. We define a “pseudo-entrainment” diagnostic that combines subsaturation and a θ es measure of conditional instability similar to what entrainment would produce under the small-buoyancy approximation. This captures the trade-off between larger θ es lapse rates (entrainment of dry air) and small subsaturation (permits positive buoyancy despite high entrainment). This pseudo-entrainment diagnostic is also a reasonable indicator of the critical value of integrated buoyancy for precipitation onset. Models with poor θ e /θ es structure (those using variants of the Tiedtke Scheme) or low entrainment runs of CAM5, and models with low subsaturation, such as NASA-GISS, lie outside the observational range in this diagnostic.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Phase field modeling of hydrogen release in nitrile-butadiene rubber composites after high-pressure hydrogen exposure

Experimental data showed that a simple diffusion model doesn’t suffice to explain the hydrogen (H 2 ) release kinetics in nitrile butadiene rubber (NBR) composites (one with no filler or plasticizer, one with plasticizer only, one with fillers and plasticizer, and one with fillers only) after exposed to hydrogen gas at a pressure of 27.6 MPa. In this work, we developed a phase field model that considers the effect of H 2 gas bubble evolution on H 2 release kinetics. With the model, we simulated the effect of thermodynamic and kinetics properties on the nucleation, growth, and shrinkage of gas bubbles and H 2 release kinetics during and after decompression. Here, the results demonstrated that 1) the model built upon certain thermodynamic and kinetics properties of a given material can well describe its H 2 release kinetics measured in experiments, and 2) the model can predict the gas bubble evolution, which is associated with the material property degradation and failures, with accurate thermodynamic and kinetics properties or H 2 release data from experiments.

08 HYDROGEN↗
DOE OSTI2024

Elucidating the complex interplay between thermodynamics, kinetics, and electrochemistry in battery electrodes through phase-field modeling

Abstract This article highlights applications of phase-field modeling to electrochemical systems, with a focus on battery electrodes. We first provide an overview on the physical processes involved in electrochemical systems and applications of the phase-field approach to understand the thermodynamic and kinetic mechanisms underlying these processes. We employ two examples to highlight how realistic thermodynamics and kinetics can naturally be incorporated into phase-field modeling of electrochemical processes. One is a composite battery cathode with an intercalation compound (Li x FePO 4 ) as the electrochemically active material, and the other is a displacement reaction compound (Li–Cu–TiS 2 ). With the input parameters mostly from atomistic calculations and experimental measurements, phase-field simulations allowed us to untangle the interactions among transport, reaction, electricity, chemistry, and thermodynamics that lead to highly complex evolution of the materials within battery electrodes. The implications of these observations for battery performance and degradation are discussed. Graphical abstract

Andrews, W. Beck (ORCID:0000000287824621)↗
DOE OSTI2024

Chemistry Informed Machine Learning-Based Heat Capacity Prediction of Solid Mixed Oxides

Knowing heat capacity is crucial for modeling temperature changes with the absorption and release of heat and for calculating the thermal energy storage capacity of oxide mixtures with energy applications. The current prediction methods (ab initio simulations, computational thermodynamics, and the Neumann–Kopp rule) are computationally expensive, not fully generalizable, or inaccurate. Machine learning has the potential of being fast, accurate, and generalizable, but it has been scarcely used to predict mixture properties, particularly for mixed oxides. Here, we demonstrate a method for the generalizable prediction of heat capacity of solid oxide pseudobinary mixtures using heat capacity data obtained from computational thermodynamics and descriptors from ab initio databases. Further, models trained through this workflow achieved an error (mean absolute error of 0.43 J mol –1 K –1 ) lower than the uncertainty in differential scanning calorimetry measurements, and the workflow can be extended to predict other properties derived from the Gibbs free energy and for higher-order oxide mixtures.

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

Discrete generative diffusion models without stochastic differential equations: A tensor network approach

Diffusion models (DMs) are a class of generative machine learning methods that sample a target distribution by transforming samples of a trivial (often Gaussian) distribution using a learned stochastic differential equation. In standard DMs, this is done by learning a “score function” that reverses the effect of adding diffusive noise to the distribution of interest. Here we consider the generalisation of DMs to lattice systems with discrete degrees of freedom, and where noise is added via Markov chain jump dynamics. We show how to use tensor networks (TNs) to efficiently define and sample such “discrete diffusion models” (DDMs) without explicitly having to solve a stochastic differential equation. We show the following: (i) by parametrising the data and evolution operators as TNs, the denoising dynamics can be represented exactly; (ii) the auto-regressive nature of TNs allows to generate samples efficiently and without bias; (iii) for sampling Boltzmann-like distributions, TNs allow to construct an efficient learning scheme that integrates well with Monte Carlo. We illustrate this approach to study the equilibrium of two models with non-trivial thermodynamics, the d = 1 constrained Fredkin chain and the d = 2 Ising model. Published by the American Physical Society 2025

Causer, Luke (ORCID:0000000194243473)↗
DOE OSTI2020

Three-dimensional modeling and simulations of single-crystal and bi-crystal titanium for high-strain-rate loading conditions

High purity single crystal titanium (Ti) under shock wave loading is modeled under both one-dimensional and three-dimensional cylindrical conditions. Cylinder sizes of 10 μm and 20 μm radius are both considered in order to assess influence of boundary conditions. A thermodynamically consistent single crystal model for application to shock conditions is presented. The model accounts for the coupled non-linear elastic, dislocation slip, deformation twinning, and structural phase transformation response of the titanium material. Plate impact experiment results using a copper flyer are used to compare against the simulations for crystals oriented in [0001] and [10$\overline1$1] crystallographic directions. The one-dimensional and threedimensional simulations of the two differently oriented single crystals indicate differences between the one-dimensional and three-dimensional representation, especially for the [10$\overline1$1] oriented single crystal. This orientation breaks the relative orientation symmetry between the crystal and cylinder which otherwise exists for the [0001] oriented single crystal. A significant amount of heterogeneity in the field response of the[10$\overline1$1] oriented simulation was demonstrated due to the highly coupled nature of the deformation. A bi-crystal model composed of both the [0001] and [10$\overline1$1] orientations with the boundary between the two along the axis of the cylinder is also considered for a cylinder model size of 10 μm. The results indicate a strong interaction between the two grains that affects the ω phase volume fraction achieved relative to the single crystal calculations.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Observing and Modeling the Sequential Pairwise Reactions that Drive Solid‐State Ceramic Synthesis

Abstract Solid‐state synthesis from powder precursors is the primary processing route to advanced multicomponent ceramic materials. Designing reaction conditions and precursors for ceramic synthesis can be a laborious, trial‐and‐error process, as heterogeneous mixtures of precursors often evolve through a complicated series of reaction intermediates. Here, ab initio thermodynamics is used to model which pair of precursors has the most reactive interface, enabling the understanding and anticipation of which non‐equilibrium intermediates form in the early stages of a solid‐state reaction. In situ X‐ray diffraction and in situ electron microscopy are then used to observe how these initial intermediates influence phase evolution in the synthesis of the classic high‐temperature superconductor YBa 2 Cu 3 O 6+ x (YBCO). The model developed herein rationalizes how the replacement of the traditional BaCO 3 precursor with BaO 2 redirects phase evolution through a low‐temperature eutectic melt, facilitating the formation of YBCO in 30 min instead of 12+ h. Precursor selection plays an important role in tuning the thermodynamics of interfacial reactions and emerges as an important design parameter in planning kinetically favorable synthesis pathways to complex ceramic materials.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Snow thermal conductivity and conductive flux in the Central Arctic: Estimates from observations and implications for models

During the Arctic winter, the conductive heat flux through the sea ice and snow balances the radiative and turbulent heat fluxes at the surface. Snow on sea ice is a thermal insulator that reduces the magnitude of the conductive flux. The thermal conductivity of snow, that is, how readily energy is conducted, is known to vary significantly in time and space from observations, but most forecast and climate models use a constant value. This work begins with a demonstration of the importance of snow thermal conductivity in a regional coupled forecast model. Varying snow thermal conductivity impacts the magnitudes of all surface fluxes, not just conduction, and their responses to atmospheric forcing. Given the importance of snow thermal conductivity in models, we use observations from sea ice mass balance buoys installed during the Multidisciplinary drifting Observatory for the Study of Arctic Climate expedition to derive the profiles of thermal conductivity, density, and conductive flux. From 13 sites, median snow thermal conductivity ranges from 0.33 W m −1 K −1 to 0.47 W m −1 K −1 with a median from all data of 0.39 W m −1 K −1 from October to February. In terms of surface energy budget closure, estimated conductive fluxes are generally smaller than the net atmospheric flux by as much as 20 W m −2 , but the average residual during winter is −6 W m −2 , which is within the uncertainties. The spatial variability of conductive heat flux is highest during clear and cold time periods. Higher surface temperature, which often occurs during cloudy conditions, and thicker snowpacks reduce temporal and spatial variability. These relationships are compared between observations and the coupled forecast model, emphasizing both the importance and challenge of describing thermodynamic parameters of snow cover for modeling the Arctic as a coupled system.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

A Novel Standard Gibbs Energy of Formation Model for High-Enthalpy Water Systems

This work provides an advanced standard molar Gibbs energy of formation model, informed by molecular statistical thermodynamics (MST), and validated by mineral solubility and ion association reactions. This new model aligns with experimental data within uncertainties and highlights the role of specific MST interactions around the critical point. This work will provide a reliable tool for the optimization of industrial processes operating at the edges of the supercritical domain.

Hall, Derek↗
DOE OSTI2025

Zavitsas’ hydration model for electrolytes only stable in the presence of another: NaAl(OH) 4 in aqueous NaOH solution

Many aqueous electrolytes are only stable in the presence of another electrolyte, such as electrolytes that are only soluble in strong acid or base. An example is sodium aluminate [NaAl(OH) 4 ], which is only stable in aqueous NaOH. This complicates developing thermodynamic parameters because it is difficult to separate the contributions of individual electrolytes in multicomponent solutions to measured bulk thermodynamic properties. The present study develops a method to determine the liquid phase parameters from solubility data for Zavitsas’ Hydration model, a model that incorporates hydration parameters into the activities of dissolved species. Furthermore, this study uses gibbsite [Al(OH) 3 ] solubility data in aqueous NaOH solution to develop Zavitsas' model parameters for NaAl(OH) 4 simultaneously with the equilibrium constants. A good fit of the solubility data was found, showing that Zavitsas’ model is effective for this system.

Bayer process↗
DOE OSTI2020

Quantitative assessment of eddy viscosity rans models for turbulent mixed convection in a differentially heated plane channel

Turbulent mixed convection between two vertical, infinite parallel plates at different temperatures is studied using various two-equation turbulence models. The numerical simulations are performed at a turbulent Reynolds number of Re τ = 150 and a Grashof number of Gr = 9.6 × 10 5 . Comparisons are made against the highly trusted DNS results. Consistent with the DNS approach, the current simulations are performed using constant properties and the Boussinesq approximation to predict the influence of buoyancy. Previous studies have provided assessments of two-equation turbulence models for various scenarios, but often rely on a qualitative “eye” test in order to determine the most appropriate model to predict a given flow. This study aims to provide a new form of quantitative assessment that accounts for both the physics captured by the turbulence model as well as the magnitude of the system response quantities (SRQ) using a modified symmetric mean absolute percent error (SMAPE) method. This method is designed to be approachable to researchers at any level and can be applied to system response quantities from multiple research fields. Uncertainty quantification is also performed to determine the discretization error for each turbulence model. Recommendations are made as to which turbulence models best capture the physics – hydrodynamically and thermally – using both local and global validation metrics. Lastly, a sensitivity analysis is performed on the damping functions used in the most accurate models. This underpins the potential of model developments and adjustments most worth pursuing for buoyant flows. Finally, this framework provides a more physics-based comparative analysis of the selected turbulence models.

42 ENGINEERING↗
DOE OSTI2021

Comparison of Conventional and Constrained Variational Methods for Computing Large-Scale Budgets and Forcing Fields

Analyses of atmospheric heat and moisture budgets serve as an effective tool to study convective characteristics over a region and to provide large-scale forcing fields for various modeling applications. This paper examines two popular methods for computing large-scale atmospheric budgets: the conventional budget method (CBM) using objectively gridded analyses based primarily on radiosonde data and the constrained variational analysis (CVA) approach which supplements vertical profiles of atmospheric fields with measurements at the top of the atmosphere and at the surface to conserve mass, water, energy, and momentum. Successful budget computations are dependent on accurate sampling and analyses of the thermodynamic state of the atmosphere and the divergence field associated with convection and the large-scale circulation that influences it. Utilizing analyses generated from data taken during Dynamics of the Madden-Julian Oscillation (DYNAMO) field campaign conducted over the central Indian Ocean from October to December 2011, we evaluate the merits of these budget approaches and examine their limitations. While many of the shortcomings of the CBM, in particular effects of sampling errors in sounding data, are effectively minimized with CVA, accurate large-scale diagnostics in CVA are dependent on reliable background fields and rainfall constraints. Finally, for the DYNAMO analyses examined, the operational model fields used as the CVA background state provided wind fields that accurately resolved the vertical structure of convection in the vicinity of Gan Island. However, biases in the model thermodynamic fields were somewhat amplified in CVA resulting in a convective environment much weaker than observed.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

High-Energy Collision of Quarks and Mesons in the Schwinger Model: From Tensor Networks to Circuit QED

With the aim of studying nonperturbative out-of-equilibrium dynamics of high-energy particle collisions on quantum simulators, we investigate the scattering dynamics of lattice quantum electrodynamics in 1 + 1 dimensions. Working in the bosonized formulation of the model and in the thermodynamic limit, we use uniform-matrix-product-state tensor networks to construct multiparticle wave-packet states, evolve them in time, and detect outgoing particles post collision. This facilitates the numerical simulation of scattering experiments in both confined and deconfined regimes of the model at different energies, giving rise to rich phenomenology, including inelastic production of quark and meson states, meson disintegration, and dynamical string formation and breaking. We obtain elastic and inelastic scattering cross sections, together with time-resolved momentum and position distributions of the outgoing particles. Furthermore, we propose an analog circuit-QED implementation of the scattering process that is native to the platform, requires minimal ingredients and approximations, and enables practical schemes for particle wave-packet preparation and evolution. Furthermore, this study highlights the role of classical and quantum simulation in enhancing our understanding of scattering processes in quantum field theories in real time.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Resolving the Heat of Trimethylaluminum and Water Atomic Layer Deposition Half-Reactions

We report that atomic layer deposition (ALD) is a surface synthesis technique that is characterized by self-limiting reactions between gas-phase precursors and a solid substrate. Although ALD processes have been demonstrated that span the periodic table, a greater understanding of the surface chemistry that affords ALD is necessary to enable greater precision, including area- and site-selective growth. We offer new insight into the thermodynamics and kinetics of the trimethylaluminum (TMA) and H 2 O ALD half-reactions with calibrated and time-resolved in situ pyroelectric calorimetry. The half-reactions produce 3.46 and 2.76 eV/Al heat, respectively, which is greater than the heat predicted by computational models based on crystalline Al 2 O 3 substrates and closely aligned with the heat predicted by standard heats of formation. The pyroelectric thin-film calorimeter offers submilisecond temporal resolution that uniquely and clearly resolves precursor delivery and reaction kinetics. Both half-reactions are observed to exhibit multiple kinetic rates, with average TMA half-reaction rates at least 2 orders of magnitude faster than the H 2 O half-reaction kinetics. Comparing the experimental heat with published computational literature and additional first-principles modeling highlights the need to refine our models and mechanistic understanding of even the most ubiquitous ALD reactions.

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

The Fundamental Creep Behavior Model of Gr.91 Alloy by ICME Approach (Final Report)

The current project mainly focused on the fundamental creep cracking mechanism of the Grade 91 system during the operation conditions in the advanced technologies FE power plants and build the link among Composition-Processing Parameters-Phase Stability-Microstructure-Creep Resistance. A model based on computational thermodynamics and diffusion kinetics will be developed to provide guidance on how to improve the creep resistance of the alloy system. The PI mainly investigate the Gr.91 base alloy and weldment with the Integrated Computational Materials Engineering (ICME) approach. The long-term goal of the proposed program is to develop a model for creep and fatigue resistant alloys with different elemental systems, compositions, and processing parameters of weldment and heat-treatment. In addition, new creep and fatigue resistant alloys are to be designed based on the predictions from the model. The computational as well as experimental results include three sections: 1) The compositional optimization and secondary phases evaluation regarding the creep resistance in Grade 91 steel through the CALPHAD approach. In this section, the formation of the critical secondary phases, M23C6, MX and Z-phase were well predicted based on many conditions like the temperature as well as compositions. Meanwhile, the role of different alloying elements was also considered as they will directly affect the formation of the secondary phases mentioned above based on the compositional evolution of the alloying elements. 2) An investigation of Creep Resistance in Grade 91 Steel through Computational Thermodynamics. We systematically considered the evolution of the four critical temperatures, Ac3, Ac1, the threshold of M23C6 and Z-phase, on basis of the computational approach. Meanwhile, the equilibrium cooling as well as Scheil simulations were also considered to perform further predictions based on the different cooling rates. 3) Creep lifetime and microstructural evolution of the Grade 91 steel. In this part, mainly experimental approach was performed on the real creep test from the beginning of the commercial Gr.91 alloys. The alloys were normalized, tempered, welded (Gleeble) and PWHT and finally subjected to the creep lifetime analyses. Finally, all the samples after each process were further characterized based on the OM, XRD and SEM technique to draw final conclusions.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2025

A generalizable machine learning-assisted fast Fourier transform algorithm to simulate the large strain phenomena in polycrystalline materials

Machine learning methods have shown initial promise in constitutive modeling for single crystals or homogenized polycrystals, delivering notable computational efficiency. However, existing machine learning-based constitutive models often lack generalizability, limiting their application across diverse boundary value problems. This study introduces a thermodynamics-informed artificial neural network model to accelerate rate-tangent crystal plasticity fast Fourier transform simulations for cross-scale deformation behaviors of polycrystals under complex loading. Our model integrates microstructural variability and local interactions effectively. To address local effects in each grain, we employ K-means clustering to group Gauss points within the microstructure into clusters assumed to be in similar mechanical states. This approach, based on self-clustering analysis, extends model scope from macroscopic stress response to the granular level, capturing mechanical responses and orientation evolution across grains. This reduces the number of nonlinear problems to solve, with cluster responses propagated throughout each group. The thermodynamics-based artificial neural network-extracted features are further processed using local material state clusters to account for history-dependent deformation and evolving microstructures. Additionally, representative volume element simulations with rate-tangent crystal plasticity fast Fourier transform provide reliable datasets for model training. The proposed model demonstrates high efficiency, accuracy, self-consistency, and enhanced generalizability in predicting strain–stress responses and orientation evolution at both individual grain and aggregate scales under complex loading conditions, such as biaxial tension and arbitrary loading scenarios.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Advances and Challenges in Low‐Temperature Upcycling of Waste Polyolefins via Tandem Catalysis

Abstract Polyolefin waste is the largest polymer waste stream that could potentially serve as an advantageous hydrocarbon feedstock. Upcycling polyolefins poses significant challenges due to their inherent kinetic and thermodynamic stability. Traditional methods, such as thermal and catalytic cracking, are straightforward but require temperatures exceeding 400 °C for complete conversion because of thermodynamic constraints. We summarize and critically compare recent advances in upgrading spent polyolefins and model reactants via kinetic (and thermodynamic) coupling of the endothermic C─C bond cleavage of polyolefins with exothermic reactions including hydrogenation, hydrogenolysis, metathesis, cyclization, oxidation, and alkylation. These approaches enable complete conversion to desired products at low temperatures (<300 °C). The goal is to identify challenges and possible pathways for catalytic conversions that minimize energy and carbon footprints.

Zhang, Wei [State Key Laboratory of Petroleum Mole↗
DOE OSTI2025

Advances and Challenges in Low‐Temperature Upcycling of Waste Polyolefins via Tandem Catalysis

Abstract Polyolefin waste is the largest polymer waste stream that could potentially serve as an advantageous hydrocarbon feedstock. Upcycling polyolefins poses significant challenges due to their inherent kinetic and thermodynamic stability. Traditional methods, such as thermal and catalytic cracking, are straightforward but require temperatures exceeding 400 °C for complete conversion because of thermodynamic constraints. We summarize and critically compare recent advances in upgrading spent polyolefins and model reactants via kinetic (and thermodynamic) coupling of the endothermic C─C bond cleavage of polyolefins with exothermic reactions including hydrogenation, hydrogenolysis, metathesis, cyclization, oxidation, and alkylation. These approaches enable complete conversion to desired products at low temperatures (<300 °C). The goal is to identify challenges and possible pathways for catalytic conversions that minimize energy and carbon footprints.

C-C cleavage↗