Engineering PapersSearch

SEARCH · Engineering Papers

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 19 records

LASSO for CALPHAD Model Selection Enables Data-Efficient Thermodynamic Modeling: An Application in Thermochemical Hydrogen Production Materials

Phenomenological CALPHAD (CALculation of PHAse Diagrams) models, widely used for multicomponent materials, often contain a considerable number of parameters and require fitting using data from a relatively small number of experimental measurements or theoretical calculations. Sometimes these parameters are introduced for the purpose of improving model fits but without clear physical justification, which leads to overparametrized models with poor generalization performance. Automated approaches for optimal model selection based on the available data therefore become critical. Here, in this work, a least absolute shrinkage and selection operator (LASSO)-based approach is developed for model selection by leveraging the linearity of the CALPHAD model with respect to its parameters to convert the model selection and fitting to a LASSO minimization problem. We demonstrate its utility for thermodynamic modeling of thermochemical hydrogen (TCH) production materials using lanthanum strontium manganite (LSM) as an example. Various TCH-relevant properties, including oxygen stoichiometry as a function of oxygen partial pressure, enthalpy of reduction, and entropy of reduction, are successfully predicted with reasonable accuracy using a minimal set of model parameters. Importantly, the model selection and fitting involve minimal human decision; it can therefore be applied to high-throughput DFT defect calculations and yield efficient workflows for TCH material modeling and optimization.

CALPHAD

Guiding Principles for Geochemical/Thermodynamic Model Development and Validation in Nuclear Waste Disposal: A Close Examination of Recent Thermodynamic Models for H + —Nd 3+ —NO 3 - (—Oxalate) Systems

Development of a defensible source-term model (STM), usually a thermodynamical model for radionuclide solubility calculations, is critical to a performance assessment (PA) of a geologic repository for nuclear waste disposal. Such a model is generally subjected to rigorous regulatory scrutiny. In this article, we highlight key guiding principles for STM model development and validation in nuclear waste management. We illustrate these principles by closely examining three recently developed thermodynamic models with the Pitzer formulism for aqueous H + —Nd 3+ —NO 3 - (—oxalate) systems in a reverse alphabetical order of the authors: the XW model developed by Xiong and Wang, the OWC model developed by Oakes et al., and the GLC model developed by Guignot et al., among which the XW model deals with trace activity coefficients for Nd(III), while the OWC and GLC models are for concentrated Nd(NO 3 ) 3 electrolyte solutions. The principles highlighted include the following: (1) Principle 1. Validation against independent experimental data: A model should be validated against experimental data or field observations that have not been used in the original model parameterization. We tested the XW model against multiple independent experimental data sets including electromotive force (EMF), solubility, water vapor, and water activity measurements. The results show that the XW model is accurate and valid for its intended use for predicting trace activity coefficients and therefore Nd solubility in repository environments. (2) Principle 2. Testing for relevant and sensitive variables: Solution pH is such a variable for an STM and easily acquirable. All three models are checked for their ability to predict pH conditions in Nd(NO 3 ) 3 electrolyte solutions. The OWC model fails to provide a reasonable estimate for solution pH conditions, thus casting serious doubt on its validity for a source-term calculation. In contrast, both the XW and GLC models predict close-to-neutral pH values, in agreement with experimental measurements. (3) Principle 3. Honoring physical constraints: Upon close examination, it is found that the Nd(III)-NO 3 association schema in the OWC model suffers from two shortcomings. Firstly, its second stepwise stability constant for Nd(NO 3 ) 2+ (log K 2 ) is much higher than the first stepwise stability constant for NdNO 3 2+ (log K 1 ), thus violating the general rule of (log K 2 –log K 1 ) < 0, or $\frac{K1}{K2}$>1. Secondly, the OWC model predicts abnormally high activity coefficients for Nd(NO 3 ) 2 + (up to ~900) as the concentration increases. (4) Principle 4. Minimizing degrees of freedom for model fitting: The OWC model with nine fitted parameters is compared with the GLC model with five fitted parameters, as both models apply to the concentrated region for Nd(NO 3 ) 3 electrolyte solutions. The latter appears superior to the former because the latter can fit osmotic coefficient data equally well with fewer model parameters. The work presented here thus illustrates the salient points of geochemical model development, selection, and validation in nuclear waste management.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W

Thermodynamic modeling of aqueous acetic acid, butyric acid, and lactic acid solutions

Based on the activity coefficient – fugacity coefficient approach, a rigorous thermodynamic modeling study is presented for accurate correlation of vapor-liquid equilibrium data of aqueous solutions of acetic acid (293 to 391 K), butyric acid (325 to 436 K), lactic acid (378 to 409 K), and acetic acid + butyric acid binary mixture (358 to 421 K). In addition, the pH data of the three aqueous, single carboxylic acid solutions were measured at 298 to 328 K and successfully correlated. Given that these aqueous carboxylic acid solutions exhibit various degrees of association behavior in both vapor and liquid phases, the thermodynamic models considered for this study include the Redlich-Kwong equation of state (RK-EoS) and the Hayden-O’Connell equation of state (HOC-EoS) for the vapor phase fugacity coefficients and the electrolyte non-random two-liquid model (eNRTL) and the association electrolyte non-random two-liquid model (AeNRTL) for the liquid phase activity coefficients. The combination of the HOC-EoS for the vapor phase and the AeNRTL model for the liquid phase is found to provide the best correlation results, consistent with the fact that the HOC-EoS and the AeNRTL model explicitly account for association behaviors in the vapor phase and the liquid phase, respectively.

09 BIOMASS FUELS

REDOTHERM: a thermodynamic modeling framework for redox-based thermochemical processes

Two-step thermochemical redox cycles are being developed as a potential pathway for the production of hydrogen and syngas. While there are many possible reactor and system configurations, moving oxide systems are considered promising in terms of the redox thermodynamics, due to the potential implementation of a countercurrent system that can achieve higher performance compared to other configurations. There is a lack of a robust thermodynamic modeling framework in the field, with multiple models incorporating incorrect thermodynamic assumptions that violate the second law of thermodynamics. We present in this work REDOTHERM, an open-source system model for moving oxides that incorporates the correct thermodynamic limits, as well as various options for the system auxiliary units including product separation, heat recovery, and oxygen removal. The model is agnostic to the energy source, and could be used for solar thermal or other configurations. We highlight the uses of this model, presenting some of the tradeoffs and challenges in redox-active material selection and how they affect the entire thermochemical hydrogen production process. This model could be easily adapted and used for material exploration, system/reactor design, and technoeconomic analysis.

08 HYDROGEN

Thermodynamic Modeling of Point Defects in Triple Conducting Perovskite Ba 0.95 La 0.05 FeO 3- $_δ$ with Incorporation of the Hydride Defect Formation Reaction for Solid Oxide Cells

Distinct from the proton defect, the hydride defect species may be present in certain perovskite materials in reducing environments such as in fuel electrodes of solid oxide cells (SOCs) or on the reducing side of ceramic membranes. A generalized defect thermodynamic model was developed for the triple-conducting perovskites (La,Ba)Fe 1-x M x O 3-δ (M = Y and Zr) to allow inclusion of the hydride defect formation reaction in addition to the other three main defect reactions, namely, the oxygen vacancy formation, hydration, and charge disproportionation reactions. This comprehensive defect model also allows the incorporation of polynomial functional forms of oxygen nonstoichiometry δ to describe the defect reaction energies and entropies and to enable refinements of the defect reaction equilibrium constants in the defect thermodynamic analysis. As a first step, the developed model is applied to the Ba 0.95 La 0.05 FeO 3-δ material as an illustrative system to obtain its Brouwer diagrams with both the proton and hydride defects in relevant SOC conditions, particularly for more reducing environments. In conclusion, the results provide direct guidance on the influence of electronic and ionic defect concentrations upon thermodynamic properties and ultimately on the performance of Ba 0.95 La 0.05 FeO 3-δ and potentially other (La,Ba)Fe 1-x M x O 3-δ perovskite materials involved in SOC applications.

25 ENERGY STORAGE

Thermodynamic modeling of countercurrent chemical looping reverse water gas shift process for redox material screening

The reverse water gas shift (RWGS) reaction is a key pathway for CO 2 utilization, particularly within Power-to-X process chains aimed at sustainable fuel and chemical production. Countercurrent chemical looping (CL-RWGS) using non-stoichiometric oxides can overcome equilibrium limitations of conventional RWGS reactors, enabling significantly higher CO 2 conversions. However, modeling the limiting performance of such systems is challenging due to their multiphase nature and coupled spatial and temporal variation in chemical composition. In this work, we present a discretized batch equilibrium model that simulates CL-RWGS reactors as a series of localized equilibrium exchanges between gas and solid elements. The model is numerically stable, computationally efficient, and free of kinetic source terms, making it well-suited for parametric studies and system-level integration. It is validated against established convection–diffusion models and shown to predict reasonable upper bounds on experimental results. Application of the model to a range of oxygen carrier materials identifies cerium–zirconium solid solutions, particularly Ce 0.80 Zr 0.20 O 2 , as a promising class offering superior oxygen storage characteristics compared to state-of-the-art La 0.6 Sr 0.4 FeO 3 . This framework provides a robust platform for materials screening, reactor sizing, and performance optimization in chemical looping systems. The model implementation is available as open-source software to support further research and development.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH

Thermodynamic Modeling of Complex Solid Solutions in the Lu-H-N System via Graph Neural Network Accelerated Monte Carlo Simulations

Metal hydrides are important across diverse applications, such as hydrogen storage, batteries, gas sensors, nuclear reactions, and high-temperature superconductivity. Previous computational studies of metal hydrides under extreme pressures, e.g., 𝑂⁡(10 2 ) ⁢GPa, usually treat them as stoichiometric compounds without considering interstitial lattice disorder. As pressures become more moderate in the 𝑂⁡(10 0 ) ⁢GPa and below range, hydrogen disorder at interstitial lattice sites becomes prominent, e.g., in the N-doped Lu hydride that was recently claimed superconducting near 1 GPa. Further adding compositional complexity from alloying and/or multielement interstitial occupation makes elucidating pressure- and temperature-dependent observables intractable by first-principles calculations alone. We therefore propose a lattice graph neural-network surrogate modeling approach to predict configuration- and pressure-dependent equation-of-state properties. Their efficiency permits Monte Carlo simulations to calculate Gibbs energies and pressure-dependent phase diagrams, thereby revealing insights into the synthesis conditions required for achieving desired phase equilibria. We demonstrate this concept for the compositionally complex cubic Lu(H,N,Va) 3 system where three constituents (hydrogen, nitrogen and vacancy) have disordered multielement interstitial occupancies and insights into pressure-dependent phase equilibria are critically needed, e.g., N-doping levels can significantly lower dehydrogenation temperatures and provide a new strategy to optimize hydrogen-storage alloys. This work can improve the thermodynamic understanding of the Lu-H-N system and help rational synthesis of N-doped Lu hydrides, but more generally demonstrates an efficient approach to model pressure-dependent thermodynamics of multicomponent solid solutions.

Monte Carlo methods

Solvent-mediated contaminant removal from plastic waste using thermodynamic modeling

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

Zhou, Panzheng [University of, Wisconsin, Madison,

Thermodynamic Modeling of Intrinsic Defects in MnBi₂Te₄

This repository contains the computational data supporting the manuscript titled “The critical role of intrinsic defects and many-body interactions on the stability of MnBi₂Te₄.” It includes: 1. DFT data generated using VASP, used for training and benchmarking electronic structure models. 2. Quantum Monte Carlo (QMC) data produced with QMCPACK, used to apply many-body corrections and validate the electronic and magnetic properties of MnBi₂Te₄. 3. Relevant scripts used to run, analyze, and process the calculations, enabling reproducibility and transparency of the workflows.

36 MATERIALS SCIENCE

Development of a kinetic-thermodynamic model for lime-stabilization of Na-bentonite

This study presents the first kinetic model to predict the solid and pore solution composition of Na-bentonite clay reacting with slaked lime over a period of 720 days. The model successfully accounts for most experimental data using a single kinetic rate constant. The following sequence of reactions was predicted by the model: initial rapid dissolution of portlandite within the first 7 days, leading to a decrease in pH and dissolved calcium, and concurrent formation of calcium silicate hydrates (C-S-H: jennite), calcium aluminate hydrate (C-A-H: C₄AH₁₃), calcium aluminosilicate hydrates (stratlingite) and hydrotalcite. After 7 days, jennite and stratlingite are predicted to transform into tobermorite-II, contributing to strength development up to 28 days. From 28 to 90 days, continued montmorillonite dissolution is predicted, along with minor formation of ettringite, partial tobermorite-II dissolution, and precipitation of secondary phases such as albite and talc. Experimentally, portlandite dissolution was confirmed by TGA and XRD and found to be complete within 7 days, in agreement with model predictions. However, other predicted solid-phase transformations (e.g., tobermorite-II formation and dissolution, ettringite, albite, and talc formation) could not be conclusively verified through experimental techniques. Aqueous phase measurements confirmed that the pH and Ca trends in solution, and that equilibrium was reached by 90 days.

Chemical kinetics

PyJMAK: An Open-Source Python Toolkit for Modeling Solid-State Metallurgical Phase Transformations

Accurate prediction of metallurgical phase transformations is an essential basis for autonomous optimization and rapid part qualification. Several methods can be used to estimate the evolution of phase fractions such as JMAK kinetics-based models, phase-field models, thermodynamic models, and data-driven machine learning models. Thermodynamic and phase-field-based methodologies solve multiphysics equations requiring numerous calibration parameters and significant computational resources. As a result, the computation domain is limited to a point or on order of micron-meters. The data-driven models rely on large datasets from experiments and simulations. While the JMAK model only provides information about phase fraction evolution, it can predict this evolution in near real-time using thermal history and thermodynamic data without restriction on the domain. JMAK models have been popularly used by researchers to model phase transformations occuring during additive manufacturing or over arbitrary temperature profiles. Commercial proprietary software such as Abaqus and Ansys or closed-source in-house implementations offer the ability to model JMAK based kinetics to predict phase transformation. However, these software packages are not open-source or freely available for use and development in conjunction with manufacturing machines, sensors, and machine learning algorithms. In addition, the use of the model is restricted by a license token. In contrast, given temperature profiles at multiple points in the domain, this Python-based PyJMAK model can compute phase evolution in parallel due to its stand-alone modular, voxel-based structure, and it can be executed on high-performance computing resources without any license restrictions.

Prabhune, Bhagya [Oak Ridge National Laboratory (O

REDOTHERM (Redox Countercurrent Thermodynamic Limits Model) [SWR-24-88]

REDOTHERM is an open-source, MATLAB-based thermodynamic modeling framework developed to evaluate the performance of redox-active materials for water (H2O) and carbon dioxide (CO2) splitting. It includes models of all major unit operations and supports comparative analysis of different redox-active material candidates. The model is tailored for systems of moving oxide under a parallel/cocurrent flow (PF) and countercurrent flow (CF) configurations. Unvalidated mixed flow reactor (MFR, also known as CSTR) model is also included as an optional addition.

Lidor, Alon [National Renewable Energy Laboratory

Modeling Pb(II) Adsorption on Mineral Surfaces: Bridging Density Functional Theory and Experiment with Thermodynamic Insights

Despite decades of work on aqueous lead (Pb) adsorption on a-Fe2O3 (hematite) and a-Al2O3 (alumina), gaps between measurements and modeling obscure molecular-level understanding. Achieving well-matched geometries between theory and experimental for mineral-water interfaces is a hurdle, as surface functional group type and distribution must be accounted for in determining mechanisms. Additionally, computational methods that can describe the substrate are often not appropriate to capture aqueous effects. Progress requires focusing on well-studied and relevant systems, such as key facets (001),(012), and (110) of hematite and alumina, and ubiquitous contaminants such as aqueous Pb. In the past, bulk-parametrized bond-valence principles were used to rationalize Pb(II) adsorption trends. These approaches can break down at surfaces, where flexible bonding environments and adsorption-induced surface relaxations play a critical role. Here, we adapt and apply a density functional theory (DFT) and thermodynamics framework, integrating DFT-calculated energies with experimental data and electrochemical principles, to predict Pb(II) adsorption. Our model results capture trends across the full set of surfaces and predict that inner-sphere Pb(II) sorption on (001) alumina varies from unfavorable to weakly favorable across a range of pH conditions. This aligns with experiment insights that Pb(II) interacts at that surface through outer-sphere interactions. Extending to Fe(II) adsorption, we demonstrate a coverage-dependent site preference, potentially explaining disorder in overlayers grown by the oxidative adsorption of Fe(II) on hematite (001).

lead contamination

Toward an improved nonlocal thermodynamic equilibrium model for more predictive simulations of ignition scale hohlraums

Recently, nonlocal thermodynamic equilibrium (NLTE) modeling has been identified as the primary reason for discrepant predictions of the peak neutron production time in indirectly driven inertial confinement fusion (ICF) platforms. It has also been observed that predictions of collisional excitation rates differ by as much as 50% from measurements. Theoretical uncertainties in dielectronic recombination rates have also been posited as possibly contributing to errors in NLTE predictions. This work examines the impact of multipliers on collisional excitation and dielectronic recombination rates on simulations of a directly driven gold sphere and an indirect drive ICF implosion. It is found that multipliers on the collisional excitation rates have a strong impact on radiant intensity and electron temperature and a weaker impact on ionization state, whereas multipliers on dielectronic recombinations rates strongly impact ionization state with a smaller impact on radiant intensity and electron temperature. A self-consistent NLTE model which places multipliers on differing transitions, as motivated by experimental measurements and more detailed atomic physics predictions, improves agreement but does not completely eliminate discrepancies with measurements of the radiant intensity within the 2–4 keV spectral range.

Farmer, W. A. [Lawrence Livermore National Laborat

Adsorption Thermodynamics for Process Simulation

Adsorption has rapidly evolved in recent decades and is an established separation technology extensively practiced in gas separation industries and others. However, rigorous thermodynamic modeling of multicomponent adsorption equilibrium remains elusive, and industrial practitioners rely heavily on expensive and time-consuming trial-and-error pilot studies to develop adsorption units. Here, this article highlights the need for rigorous adsorption thermodynamic models and the limitations and deficiencies of existing models such as the extended Langmuir isotherm, dual-process Langmuir isotherm, and adsorbed solution theory. It further presents a series of recent advances in the generalization of the classical Langmuir isotherm of single-component adsorption by deriving an activity coefficient model to account for the adsorbed phase adsorbate–adsorbent interactions, substituting adsorbed phase adsorbate and vacant site concentrations with activities, and extending to multicomponent competitive adsorption equilibrium, both monolayer and multilayer. Requiring a minimum set of physically meaningful model parameters, the generalized Langmuir isotherm for monolayer adsorption and the generalized Brunauer–Emmett–Teller isotherm for multilayer adsorption address various thermodynamic modeling challenges including adsorbent surface heterogeneity, isosteric enthalpies of adsorption, BET surface areas, adsorbed phase nonideality, adsorption azeotrope formation, and multilayer adsorption. Also discussed is the importance of quality adsorption data that cover sufficient temperature, pressure, and composition ranges for reliable determination of the model parameters to support adsorption process simulation, design, and optimization.

09 BIOMASS FUELS