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.
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.
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.
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.
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.
Computational Design of Eutectic Molten Salt Mixtures: What Can Thermodynamic Models Do?
In the search for efficient energy storage battery technologies, designing stable electrolytes has been a long-standing challenge. Electrolytes based on molten salt eutectics are known for their stability with minimum parasitic reactions when compared to their widely used organic counterparts. However, the operating temperatures of these molten salt electrolyte-based batteries are dictated by the melting point of the eutectic mixtures. Design and high throughput screening of low melting temperature eutectic molten salt mixtures have been hindered by the lack of computational models. In this work, we develop thermodynamic models to predict the eutectic points of several molten salt mixtures. The framework of the COSMO-SAC model is used for the predictions and is compared with experimental data and other thermodynamic approaches. Rapid thermodynamics-based approaches, as shown in this study, can accelerate the discovery of new materials, complementing experimental techniques.
Thermodynamic models of the chemistry of lunar volcanic gases
Thermodynamic models and mass-balance arguments are used to constrain the chemistry of lunar volcanic gases. The results predict that lunar gases were dominated by reduced C and S gases such as CO, COS, CS2, S2. The more oxidized gases CO2 and SO2 were also important, but only in limited temperature ranges. Gases such as Cl2, CCl4, and CF4 were more abundant than HF and HCl, which were the two major H compounds in the lunar gases. Chlorides and fluorides were important species for transporting many volatile and ore-forming metals, and the implications for fractionating and concentrating metals into lunar ore-deposits merit further study.
Thermodynamic models for bounding pressurant mass requirements of cryogenic tanks
Thermodynamic models have been formulated to predict lower and upper bounds for the mass of pressurant gas required to pressurize a cryogenic tank and then expel liquid from the tank. Limiting conditions are based on either thermal equilibrium or zero energy exchange between the pressurant gas and initial tank contents. The models are independent of gravity level and allow specification of autogenous or non-condensible pressurants. Partial liquid fill levels may be specified for initial and final conditions. Model predictions are shown to successfully bound results from limited normal-gravity tests with condensable and non-condensable pressurant gases. Representative maximum collapse factor maps are presented for liquid hydrogen to show the effects of initial and final fill level on the range of pressurant gas requirements. Maximum collapse factors occur for partial expulsions with large final liquid fill fractions.
CFD and Thermodynamic Model Predictions of No-Vent Tank Filling in Microgravity
Cryogenic Fluid Management (CFM) of propellant fluids under low gravity conditions present an important challenge in the pursuit of NASA’s future long-duration space missions. Success of these missions relies to a large extent on the ability to store and transfer the cryogenic liquid fuels efficiently and safely. To this end, the Zero-Boil-Off Tank Filling and Transfer Experiment (ZBOT-FT) undertaken by University of Bremen in collaboration with Case Western Reserve University aims at studying two phase flow tank-to-tank transfer operations in microgravity using a simulant fluid, Perfluoro-n-pentane (PNP). As part of this effort, in this work, we present a system level homogeneous thermodynamic model and a two-phase Computational Fluid Dynamics (CFD) model of a no-vent cylindrical tank filling in microgravity. Both 2D axisymmetric and 3D parametric simulations case studies are performed to show the detailed physics of the filling process and the evolution of tank pressure for filling flow rates between 1.0 and 1.65 mL/s. A zero-dimensional homogeneous thermodynamic model is also developed to predict the final steady state tank pressure analytically. Comparison of the 2D and 3D CFD predictions of final pressures to the steady state thermodynamic results displays very good agreement. However, it is shown that the 2D axisymmetric model predicts a maxima in tank pressure before reaching steady state. It is indicted that this pressure maxima is caused by the constraints inherent in the axisymmetric model and is absent in the 3D and thermodynamic pressure predictions.
Designing Molten Salt Eutectics: A Combined Thermodynamic Modeling and Machine Learning Approach
Designing stable electrolytes with target properties is an important challenge in realizing next generation energy storage devices. Molten salt eutectics-based electrolytes are known for their stability with minimal parasitic reactions when compared to traditional organic electrolytes and are an attractive option for different battery chemistries. The operating temperature of the molten salt batteries depends on the melting temperature of the eutectic and hence there is a necessity to discover novel low melting temperature molten salt eutectic mixtures for energy storage applications. In this work we develop a high throughput computational screening approach for molten salt mixtures using thermodynamic modeling and machine learning (ML). COSMO-SAC model and ML approaches were independently developed based on the existing experimental data and these models were further used to predict the eutectic melting temperature and composition of several new binary, ternary, and quaternary mixtures. We show that combining ML and thermodynamic modeling strategies is effective in exploring the vast design space of molten salt mixtures.
A nonequilibrium thermodynamic model of ion transport in a three-compartment system
Nonequilibrium thermodynamic model of ion transport in three-compartment system
Analysis of Experimental Ice Accretion Data and Assessment of a Thermodynamic Model During Ice Crystal Icing
This paper evaluates a thermodynamic ice crystal icing model that has been previously presented to describe the possible mechanisms of icing within the core of a turbofan jet engine. The model functions between two distinct ice accretions based on a surface energy balance: freeze-dominated icing and melt-dominated icing. Freeze-dominated icing occurs when liquid water (from melted ice crystals) freezes and accretes on a surface along with the existing ice of the impinging water and ice mass. This freeze-dominated icing is characterized as having strong adhesion to the surface. The amount of ice accretion is partially dictated by a freeze fraction, which is the fraction of impinging liquid water that freezes. Melt-dominated icing occurs as unmelted ice on a surface accumulates. This melt-dominated icing is characterized by weakly bonded surface adhesion. The amount of ice accumulation is partially dictated by a melt fraction, which is the fraction of impinging ice crystals that melts. Experimentally observed ice growth rates suggest that only a small fraction of the impinging ice remains on the surface, implying a mass loss mechanism such as splash, runback, bounce, or erosion. The fraction of mass loss must be determined in conjunction with the fraction of freezing liquid water or fraction of melting ice on an icing surface for a given ice growth rate. This mass loss parameter, however, along with the freeze fraction and melt fraction, are the only experimental parameters that are currently not measured directly. Using icing growth rates from ice crystal icing experiments, a methodology that has been previously proposed is used to determine these unknown parameters. This work takes ice accretion data from tests conducted by the National Aeronautics and Space Administration (NASA) at the Glenn Research Center in 2018 that examined the fundamental physics of ice crystal icing. This paper continues evaluation of the thermodynamic model from a previous effort, with additions to the model that account for sub-freezing temperatures that have been observed at the leading edge of the airfoil during icing. The predicted temperatures were generally in good agreement with measured temperatures. Other key findings include the total wet-bulb temperature being a good first order indicator of whether icing is freeze-dominated (sub-freezing values) or melt-dominated (above freezing). Maximum sticking efficiency values, the fraction of impinging mass that adheres to a surface, was calculated to be about 0.2, and retained this maximum value for a range of melt ratios (0.3 to 0.65 and possibly higher), which is defined as the ratio of liquid water content to total water content. Higher air velocities reduced the maximum sticking efficiency and shifted the icing regime to higher melt ratio values. Finally, the leading edge ice accretion angle was found to be related to ice growth (lower growth rates for smaller angles) and melt ratio (smaller melt ratios resulted in smaller angles, likely due to erosion effects).
Proton production in relativistic heavy ion collisions; comparison with a thermodynamical model
Experimental results concerning proton production in nuclear collisions, obtained at Saturne with the Diogene 4 pi facility, are compared with the predictions of a thermodynamical model, using collective velocity distributions combined with a statistical thermodynamics in local rest frames. Experimental differential cross sections for alpha + nucleus and Neon + nucleus central collisions at incident energies between 200 and 800 MeV per nucleon are well reproduced by the model, for an angular range 30-110 degrees in the laboratory system. Extracted values of the temperatures are compared with those given by other authors.
Improved thermodynamic modeling of the no-vent fill process and correlation with experimental data
The United States' plans to establish a permanent manned presence in space and to explore the Solar System created the need to efficiently handle large quantities of subcritical cryogenic fluids, particularly propellants such as liquid hydrogen and liquid oxygen, in low- to zero-gravity environments. One of the key technologies to be developed for fluid handling is the ability to transfer the cryogens between storage and spacecraft tanks. The no-vent fill method was identified as one way to perform this transfer. In order to understand how to apply this method, a model of the no-vent fill process is being developed and correlated with experimental data. The verified models then can be used to design and analyze configurations for tankage and subcritical fluid depots. The development of an improved macroscopic thermodynamic model is discussed of the no-vent fill process and the analytical results from the computer program implementation of the model are correlated with experimental results for two different test tanks.
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.
Improved thermodynamic modelling of the no-vent fill process and correlation with experimental data
The United States plans to establish a permanent manned presence in space and to explore the Solar System have created the need to efficiently handle large quantities of subcritical cryogenic fluids, particularly propellants such as liquid hydrogen and liquid oxygen, in low- to zero-gravity environments. One of the key technologies to be developed for fluid handling is the ability to transfer the cryogens between storage and spacecraft tanks. The no-vent fill method has been identified as one way to perform this transfer. In order to understand how to apply this method, a model of the no-vent fill process is being developed and correlated with experimental data. The verified models then can be used to design and analyze configurations for tankage and subcritical fluid depots. This paper discusses the development of an improved macroscopic thermodynamic model of the no-vent fill process and correlates the analytical results from the computer program implementation of the model with experimental results for two different test tanks at NASA Lewis Research Center.
A Case Study of AI-assisted Creation of a Thermodynamics Model of Precipitation Formation During Rapid Depressurization of a Vented Container
Precipitation may form in humid containers undergoing rapid depressurization. This precipitation may be liquid, i.e. fog, if the dewpoint is crossed above the freezing point of water, or direct snow crystallization if the dewpoint is crossed below the freezing point. Accurate modeling of this effect is potentially important for rapidly ascending vented containers in aircraft, spacecraft, and launch vehicles, as well as rapidly depressurizing vacuum chambers. A transient thermodynamics model of precipitation formation during the rapid depressurization of a container was developed in python. The model is written for a generic container and includes an optional water pool and water vapor source. Details of the model and results from several example cases spanning the full capabilities of the model, including a validation case, will be presented. Although the model is not novel, in contrast to prior works, this one was treated as a case study of the assistance of AI Large Language Models (LLMs) to create physical models. Impressions, performance, time, and cost of using AI for this task will be discussed.
Thermodynamic Modeling of Calcium and Magnesium Partitioning Between Sulfide and Silicate and Implications for Mercury’s Differentiation
Mercury is the least oxidized terrestrial planet of the solar system and is enriched in sulfur, which may be saturated in its mantle, possibly forming Ca- and Mg-rich sulfides. Using literature experimental data, we developed a thermodynamic model for calcium and magnesium partitioning between sulfide and silicate melts to enhance our understanding of the chemical equilibria involved in Mercury’s mantle-crust differentiation. Our models show that temperature and carbon abundance in sulfide have a positive effect on Ca partitioning into sulfides, while pressure decreases it. In addition, carbon abundance and sulfur concentration in the silicate melt are the primary controls on Mg partitioning, both favoring Mg partitioning into sulfides. Our results will be used to better understand the abundance as well as presence of calcium and magnesium sulfides within Mercury’s mantle and crust.