Bayesian Model Selection for Thermodynamic Models of Redox Active Metal Oxides.
Abstract not provided.
Johnston, Katherine↗
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.
Chemical degradation of fly ash blended concrete with the seasonal variation of rainwater in a radioactive waste repository: A thermodynamic modeling approach
Highlights: • The thermodynamic model for the cementitious system was successfully developed. • The chemical degradation of concrete was greatly affected by rainwater types. • Inorganic carbon species in rainwater reduced the chemical resistance of concrete. • The durability of concrete to rainwater decreased with increasing fly ash content. This study presents a long-term degradation behavior of fly ash blended concrete in a vault type low- and intermediate-level radioactive waste (LILW) repository by thermodynamic equilibrium calculations using PHREEQC combined with CEMDATA18 database. Since rainwater plays a predominant role as leachate for the concrete waste container in the repository, the effect of chemical properties of seasonal rainwater on the chemical degradation of concrete was analyzed. In addition, the impact of the blending with fly ash on the chemical degradation was evaluated through a gradual replacement of ordinary Portland cement (OPC) by fly ash. Regardless of the replacement level, autumn rainwater showed the greatest influence on the concrete degradation owing to the presence of inorganic carbon species, which induce the carbonation. The resistance of concrete to chemical degradation by rainwater was persistently reduced due to the decreasing formation of portlandite and calcium silicate hydrate (C-S-H) with increasing replacement level.
36 MATERIALS SCIENCE↗
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↗
Extended thermodynamic model for high salinity produced waters
To support the development of desalination processes for high salinity produced waters, we developed a comprehensive thermodynamic model that reliably predicts the thermodynamic properties and phase behavior of produced water. Here this study extends a previously developed thermodynamic model for the aqueous hexary oceanic salt system by incorporating the Sr 2+ and Ba 2+ ions. Based on the electrolyte nonrandom two-liquid theory, the model requires two adjustable binary interaction parameters for each water-electrolyte and electrolyte-electrolyte pair that has a common ion. The binary interaction parameters for the electrolyte-electrolyte pairs involving Sr 2+ and Ba 2+ ions were identified using thermodynamic and salt solubility data. The model was validated for temperatures from 273.15 K to 473.15 K and electrolyte concentrations from infinite dilution to salt saturation. Extension of the model to include HCO 3 − , CO 3 2− , and CO 2 is currently in progress.
42 ENGINEERING↗
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 Model and Initial Experimental Investigation of Air Dehumidification through Electrically Charged Vapor Capturing Electrostatic Droplets
The removal of water vapor from the air to reduce relative humidity is a well known indoor environmental comfort requirement. Common dehumidification approaches require a substantial amount of energy and usually involve the cooling of atmospheric humid air below its dew point or the use of absorbent/adsorbent materials to extract water vapor out of the air. More recently, researchers investigated the effect of electrostatic forces for enhancing water vapor condensation with the goal to reduce the energy consumption associated with dehumidification. However, the studies are limited, and there is a lack of correlations that can predict the dehumidification rate. In this thesis, a broad theoretical investigation to study the electrostatically enhanced condensation processes was carried out. These processes consist of the use of highly charged particles, preferably highly charged water droplets to attract polar water vapor molecules to their surfaces and promote condensation, a phenomenon known as dielectrophoresis. The electrical charge promotes the reduction of the vapor pressure on the droplets' surface with respect to the saturated pressure predicted by the Kelvin equation on curved surfaces and, consequently, the equilibrium between evaporation and condensation is shifted towards condensation. This investigation resulted in the development of a thermodynamics model, which was able to predict an effective size range of the charged droplets for optimal dehumidification, under ideal conditions. The range resulted in about 2 to 4 m in diameter, while the electrical charge was kept to the maximum limit predicted by the Rayleigh model. A sensitivity analysis on the variation of the size and the charge of the charged droplets was also performed. In terms of dehumidification rates, when six electrospray heads, i.e. the maximum number of heads tested in the experimental campaign, were considered, the model predicted very limited rates, way lower than the target of this research, i.e. a dehumidification rate in terms of relative humidity of 5 % with an air flow rate of 5 cfm. While ways to increase this rate existed, their implementation proved difficult for this initial investigation. According to the requirements predicted by the model, the use of electrosprays appeared the most suitable solution for the production of small but highly charged droplets. The electrospray features and operational modes were studied in detail and an electrospray assembly was designed to be tested in this initial experimental investigation. For this initial investigation, an experimental system was designed and built to validate the thermodynamics model and with the goal to achieve a dehumidification rate of 5 % in terms of relative humidity with an air flow rate of 5 cfm. The experimental system consisted of a wind tunnel test apparatus, dew point sensors, thermocouples for temperature measurements, pressure sensors and other devices. The experimental system was equipped with a computerized data acquisition and storage system. The system was suitable to evaluate the air water content differential before and after the test section where the electrospray heads were installed. The air water content differential was evaluated in terms of the ∆ω/ω_1 ratio, to eliminate the dependence on the dry bulb temperature along the test apparatus. This work presented initial experimental data for electrostatically enhanced dehumidification processes. Working parameters, such as air and water flow rates, high voltage potentials and polarity, deionized water types with different electrical conductivity, were varied to find the best combination for an improved dehumidification. With the set up used and the conditions considered, it was never possible to achieve a 5 % dehumidification rate with 5 cfm. In general, the dehumidification rate was always limited and lower than 1 % for all the air flow rates considered and mainly within the uncertainty of the dew point sensors. These results confirmed the already limited predictions of the thermodynamics model, for ideal conditions and the interesting and promising results obtained in some experimental investigations available in the open literature were not achievable in this thesis. The more probable reasons were the higher air flow rates considered and the complete absence of the use of cooling power, which was mainly implemented to facilitate the vapor condensation.
Morcelli, Stefano↗
Thermodynamic modeling of MDEA(aq)-NH 3 (aq)-K 2 CO 3 (aq)-CO 2 (aq) using the Extended UNIQUAC model
Here, the development of advanced processes for carbon dioxide capture requires reliable thermodynamic models to ensure accurate process design. This work presents an optimized parameter set of the Extended UNIQUAC model for advanced mixed solvent solutions containing CO 2 -MDEA-NH 3 -K 2 CO 3 in water. The model was parameterized using vapor–liquid equilibrium, pure component saturation pressure, excess enthalpy, heat of absorption, molar heat capacity, apparent molar heat capacity, and solid–liquid equilibrium. These thermal and equilibrium properties cover temperatures between 253.6 K to 600 K and pressures up to 75.6 bar. The current model represents an extension of previous versions that were suitable for modeling aqueous NH 3 and K 2 CO 3 solutions. This parametrization accounts for molecular interactions between MDEA and these two components, which allows for the design of advanced capture processes. This paper demonstrates that the new model parameters are valid within the range of operational conditions typically found in CO 2 capture processes.
42 ENGINEERING↗
Thermodynamic Modeling of Aqueous LiCl, LiBr, LiI, and LiNO 3 Solutions
Thermodynamic models are essential to facilitate the advancement of process design, optimization, and operation of electrolyte systems. In this work, a comprehensive thermodynamic framework based on the Electrolyte Nonrandom Two-Liquid model is developed to calculate phase equilibria behavior and salt solubility of aqueous LiCl, LiBr, LiI, and LiNO 3 solutions. The model describes the non-ideality of the electrolyte solutions by using two binary interaction parameters for each electrolyte-molecule pair in the system. Each binary interaction parameter is further expressed with up to three temperature coefficients which are regressed from experimental data. To take into account the hydration of lithium ion, two separate chemistries for the dissociation of lithium salts are investigated. In the first case, the lithium ion is considered as a bare ion, Li + , while in the second case hydration of the lithium ion from Li + to Li(H 2 O) + is considered. Here, the calculated thermodynamic properties compare adequately with the experimental data for both sets of chemistries for concentrations up to saturation and temperatures from 273.15 K up to 623.15 K. Moderate to significant improvements are observed with the incorporation of the hydration chemistry for aqueous LiCl, LiBr, and LiI solutions when compared to the non-hydrated lithium ion model results.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Liquid and Gamma-BCC Phase Equilibria in the Uranium-Niobium System: A Simple Thermodynamic Model for Phase Diagram and Other Predictions
A simple thermodynamic model for the uranium-niobium (U-Nb) system was created using a Calphad-type approach, circa 2002. Only the liquid and gamma-BCC phases were modeled, with their subregular solution mixing terms optimized by trial and error. The custom (user-defined) database file (.tdb) and experimental data file (.exp) are listed, both of which are compatible with Thermo-Calc software. This model is used to calculate the solidus, liquidus, miscibility gap on the U-Nb phase diagram and also the enthalpy of melting in the U-6Nb alloy.
73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
Bayesian automated weighting of aggregated DFT, MD, and experimental data for candidate thermodynamic models of aluminum with uncertainty quantification
Atomic-scale modeling methods such as density functional theory (DFT) and molecular dynamics (MD) can predict the thermodynamic properties of materials at a lower cost than experimental measurements. However, their regular usage in thermodynamic model construction is hampered by the lack of quantitative agreement with experimental measurements and the lack of uncertainty estimates on the data. To make regular usage of this atomistic simulation data, it is important to assess whether the atomistic simulation datasets, by themselves or in combination with experimental measurements, result in the same physics-informed models best supported by experimental measurements alone. Here, models of aluminum thermodynamic properties are discussed using three data sources: atomistic calculations (DFT and MD), experiments, and a combination of atomistic calculations and experiments. The study shows that, after ensuring self-consistency in predicting key invariant points, both experimental measurements and atomistic calculations can significantly contribute to an optimal model.
36 MATERIALS SCIENCE↗
Thermodynamic Modeling of Heat Engines Including Heat Transfer and Compression–Expansion Irreversibilities
In this paper, a thermodynamic model based on an endoreversible engine approach is developed to analyze the performance of heat engines operating under different thermodynamic cycles. The model considers finite heat transfer rate, variable heat source and sink temperatures, and irreversibilities associated with the expansion and compression. Expressions for the maximum power and efficiency at maximum power output are obtained as a function of hot and cold reservoir temperatures, the equivalent isentropic efficiency of compression and expansion components, and the effective conductance ratio between heat exchangers. In all cases, the Curzon–Ahlborn efficiency is retrieved at constant reservoir temperatures and neglected compression–expansion irreversibilities. The proposed model allows assessing the effect of isentropic efficiencies and heat exchanger design and operation characteristics for different thermodynamic cycles.
33 ADVANCED PROPULSION SYSTEMS↗
Sustainable bioleaching of lithium-ion batteries for critical metal recovery: Process optimization through design of experiments and thermodynamic modeling
Recycling spent lithium-ion batteries (LIBs) could alleviate supply risks for critical metals and be less harmful to the environment compared to new production of metals from mining. Developing a cost-effective LIB bioleaching process could be a promising alternative to traditional energy-intensive recycling technologies. Here, this study aimed to optimize bioleaching conditions for maximum economic competitiveness through design of experiments using iterative response surface methodology (RSM), assisted by thermodynamic modeling. The optimal condition was identified as 2.5% pulp density in 75 mM gluconic acid biolixiviant at 55°C for 30 h which could recover 57%–84% of nickel, 71%–86% of cobalt, and 100% of lithium and manganese, yielding a 17%–26% net profit margin. The recommended pulp density and acid concentrations, together with the observed metal solubilization, were supported by thermodynamic modeling predictions. Our study demonstrated that combining RSM with thermodynamic simulations could be a powerful tool for optimizing bioleaching conditions.
60 APPLIED LIFE SCIENCES↗
Boron nitride nanotube precursor formation during high-temperature synthesis: kinetic and thermodynamic modelling
In this work, we performed integrated modelling of the chemical pathways of formation for boron nitride nanotube (BNNT) precursors during high-temperature synthesis in a B/N 2 mixture. Integrated modelling includes quantum chemistry, Quantum–classical molecular dynamics, thermodynamic modelling, and kinetic approaches. We demonstrate that BN compounds are formed via the interaction of molecular nitrogen with small boron clusters, rather than through interactions with less reactive liquid boron. (This process can also be described as N 2 molecule fixation.) Liquid boron evaporates to produce these boron clusters (B m with m ≤ 5), which are subsequently converted into B m N n chains. The production of such chains is crucial to the growth of BNNTs because these chains form the building blocks of bigger and longer BN chains and rings, which are in turn the building blocks of fullborenes and BNNTs. Additionally, kinetic modelling revealed that B 4 N 4 and B 5 N 4 species in particular play a major role in the N 2 molecule fixation process. Furthermore, the formation of these species via reactions with B 4 and B 5 clusters is not adequately described under the assumption of thermodynamic equilibrium, as is demonstrated in our kinetic modelling. Thus, the accumulation of both B 4 N 4 and B 5 N 4 depends on the background gas pressure and the gas cooling rate. Long BN chains and rings, which are precursors of the fullborene and BNNT growth, form via self-assembly of components B 4 N 4 and B 5 N 4 . Our modelling results—particularly the increased densities of B 4 N 4 and B 5 N 4 species at higher gas pressures—explain the experimentally observed effect of gas pressure on the yield of high-quality BNNTs. The catalytic role of hydrogen was also studied; it is shown that HBNH molecules can be the main precursor of BNNT synthesis in the presence of hydrogen.
77 NANOSCIENCE AND NANOTECHNOLOGY↗
Thermodynamic modeling of aqueous polyelectrolyte solutions with mixed-valent counterions
Thermodynamic modeling of aqueous polyelectrolyte solutions with salts is of significant interest for many industrial applications. Here this study applies the polyelectrolyte Nonrandom Two-liquid activity coefficient model to aqueous polyelectrolyte solutions with mixed-valent counterions. A modified Delocalized Binding Theory was proposed to determine the polyion condensation fractions of the mixed counterions. This modified theory accounts for the electrostatic binding of the counterions on the polyion, the dissociation entropy of the counterions, and the electrostatic interactions between the uncondensed ionic species. Given the polyion condensation fractions, the critical value of Manning’s parameter ξ and the amounts of uncondensed polyions and counterions can be calculated along with the activity coefficients of mobile ions. The model successfully correlates experimental data for various aqueous polyelectrolyte systems with mixed-valent counterions.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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↗