Search NASA⌕ Search

SEARCH · Search NASA

Results for “Activity coefficient”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Impact of phase state and non-ideal mixing on equilibration timescales of secondary organic aerosol partitioning

Abstract. Evidence has accumulated that secondary organic aerosols (SOAs) exhibit complex morphologies with multiple phases that can adopt amorphous semisolid or glassy phase states. However, experimental analysis and numerical modeling on the formation and evolution of SOA still often employ equilibrium partitioning with an ideal mixing assumption in the particle phase. Here we apply the kinetic multilayer model of gas–particle partitioning (KM-GAP) to simulate condensation of semi-volatile species into a core–shell phase-separated particle to evaluate equilibration timescales of SOA partitioning. By varying bulk diffusivity and the activity coefficient of the condensing species in the shell, we probe the complex interplay of mass transfer kinetics and the thermodynamics of partitioning. We found that the interplay of non-ideality and phase state can impact SOA partitioning kinetics significantly. The effect of non-ideality on SOA partitioning is slight for liquid particles but becomes prominent in semisolid or solid particles. If the condensing species is miscible with a low activity coefficient in the viscous shell phase, the particle can reach equilibrium with the gas phase long before the dissolution of concentration gradients in the particle bulk. For the condensation of immiscible species with a high activity coefficient in the semisolid shell, the mass concentration in the shell may become higher or overshoot its equilibrium concentration due to slow bulk diffusion through the viscous shell for excess mass to be transferred to the core phase. Equilibration timescales are shorter for the condensation of lower-volatility species into semisolid shell; as the volatility increases, re-evaporation becomes significant as desorption is faster for volatile species than bulk diffusion in a semisolid matrix, leading to an increase in equilibration timescale. We also show that the equilibration timescale is longer in an open system relative to a closed system especially for partitioning of miscible species; hence, caution should be exercised when interpreting and extrapolating closed-system chamber experimental results to atmosphere conditions. Our results provide a possible explanation for discrepancies between experimental observations of fast particle–particle mixing and predictions of long mixing timescales in viscous particles and provide useful insights into description and treatment of SOA in aerosol models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

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↗
DOE OSTI2025

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↗
DOE OSTI2023

On the Thermodynamic Condition for Adsorption Azeotropes

Adsorption azeotropy is a common phenomenon in mixed-gas adsorption equilibria. Since the first literature report of binary adsorption azeotropes in 1933, numerous studies investigated the thermodynamic conditions for adsorption azeotropes but failed to reach definitive conclusions. Based on the generalized Langmuir isotherm model for multicomponent adsorption equilibria which takes into account the vacant site as part of the adsorbed phase, this study presents the thermodynamic condition for adsorption azeotropes as derived from the generalized Langmuir isotherm to be the equality of the ratios of adsorbed phase activity coefficient γi and adsorption equilibrium constant $K$$^{o}_{i}$ for the two adsorbates in the binary 1–2 adsorption system, i.e., γ 1 /$K$$^{o}_{1}$ = γ 2 /$K$$^{o}_{2}$. This adsorption azeotropic condition is analogous to the vapor–liquid equilibrium azeotropic condition, i.e., the equality of the products of liquid phase activity coefficient and saturation vapor pressure Pisat for the two components, i.e., γ 1 $P$$^{sat}_{1}$ = γ 2 $P$$^{sat}_{2}$. We validated the thermodynamic condition for adsorption azeotropes with 14 azeotrope-forming adsorption systems in this study. As a result, we investigated the effects of the pressure, temperature, and adsorbed phase nonideality on azeotrope formation.

42 ENGINEERING↗
DOE OSTI2023

Verification Testing of OLI Systems Mixed Solvent Electrolyte Model for the Na-K-Mg-Ca-H-Cl-SO 4 -OH-HCO 3 -CO 3 -CO 2 -H 2 ) System to High Ionic Strength at 25°C

This technical report summarizes model verification results and summary statistics for 41 evaporite mineral solubility cases evaluated by Savannah River National Laboratory using OLI Systems’ aqueous electrolyte thermodynamic modeling software. The 41 verification cases containing a total of 60 solubility curves comprise mineral solubility data from low to high ionic strength at 25°C for the eight-component system Na-K-Mg-Ca-H-Cl-SO 4 -OH-HCO 3 -CO 3 -CO 2 -H 2 O as reported by Harvie et al. (1984). Thermodynamic calculations were executed using OLI Systems’ Stream Analyzer computation module within the OLI Studio software platform (Ver. 11.0, Rev. 11.0.1.9). The Mixed Solvent Electrolyte (MSE) thermodynamic framework was chosen for this investigation because of its superiority in modeling high ionic-strength inorganic salt solutions and actinide redox chemistry and solubility, both of which are relevant to the geological repository conditions at the Waste Isolation Pilot Plant in Carlsbad, New Mexico. Mineral solubility data in various inorganic salt solutions were digitized and extracted from figures generated by Harvie et al. (1984). For each of the 60 solubility curves, a case-specific chemistry model and input file were generated in OLI Studio using OLI Stream Analyzer and the MSE (H 3 O + ion) public databank provided by OLI Systems. Model simulation results were exported to Microsoft Excel to calculate summary statistics and to generate graphs comparing the OLI model predictions to the solubility data. Summary statistics include residuals (model – data) and concordance (accuracy × precision, where precision is indicated by the Pearson correlation coefficient and accuracy accounts for bias and scale differential). Private databanks were not developed, and activity coefficient model regressions were not performed to improve OLI model fits to the data. Of the 41 model verification plots, 83% have a mean of the percent residuals less than or equal to 25%. Similarly, 75% display a concordance greater than or equal to 0.75. Only seven of the 41 verification plots fail to show good agreement between the model and data. Of these seven, three are relevant to the WIPP repository because they involve the Mg-OH-Cl-SO 4 -CO 3 aqueous system. The remaining four address salt solubilities at the pH extremes (strong acid and strong base). It should be noted that in two of the three Mg-OH-Cl-SO 4 -CO 3 system cases, the regressed Harvie et al. (1984) solubility curve also deviated from the data. Lack of agreement between the OLI model-predicted solubility curves and the data is attributable to one or more of the following: specific solid species are not included in the OLI MSE databank; there is significant variation among the different solubility datasets chosen by Harvie et al. (1984); the OLI MSE model’s thermodynamic parameters were determined using different solubility datasets; and the activity coefficient parameters for certain relevant ion-ion and ion-molecule pairs have not been optimized via data regression. Two recommendations for future work are to (1) evaluate solubility data for the Mg-OH-Cl-SO 4 -CO 3 system at high ionic strength and, if necessary, develop a private OLI MSE database that includes missing species and, where necessary, regressed standard state properties and interaction parameters; (2) perform similar verification testing of the OLI model for actinide solubility data.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2024

Development of a rate-based ENRTL-RK process model for a water-lean solvent

Advanced water-lean solvents (WLS) for post-combustion CO2 capture offer several advantages over the aqueous amine solvents . WLS have lower parasitic energy penalty, lower corrosion, lower temperature and high-pressure CO2 regeneration leading to lower cost of CO2 capture. RTI International, with funding from the US Department of Energy, has been developing its novel water-lean solvent, that has shown specific reboiler duty of 2.3 GJ/t-CO2 at the 60-kWe pilot testing unit (Tiller Plant, SINTEF, Norway) and 2.6 GJ/t-CO2 at the engineering scale testing system (12 MWe) at the Technology Centre Mongstad (TCM) in Norway. All heat duties, including the one from TCM testing, were consistent with Aspen Plus modeling of the specific configuration of each test plant. This work focuses on the development of a detailed process model using in-house laboratory measurements and process data at pilot scale. The eNTRL-RK model used in this work is based on an unsymmetric activity coefficient model with the reference states chosen to be pure liquids for solvents and ideal dilute solution at unit solute molality (resulting in activity coefficient of unity at infinite dilution) for electrolytes. It uses the Redlich-Kwong equation of state for vapor phase properties and Henry’s law for solubility of supercritical gases. The model was validated using process data from the pilot-scale campaign at the Tiller plant, and the engineering scale test campaign at TCM. Data on CO2 capture rate, absorber, and regenerator temperature profiles and specific reboiler duties from two different test campaigns at Tiller and TCM, were used to further refine and validate the model and the model compares favorably to experimental data. The validation results against TCM campaign will be presented in this work.

CO2 capture↗
DOE OSTI2023

Addressing spreading pressure dependence of real adsorbed solution theory with generalized Langmuir isotherm

This work addresses spreading pressure dependence of Real Adsorbed Solution Theory (RAST) for mixed-gas adsorption equilibria using generalized Langmuir (gL) isotherm. Considering vacant sites as an integral part of competitive multicomponent adsorption on a constant adsorbent surface area, the gL isotherm properly accounts for surface heterogeneity and loading, adsorbate composition, and temperature dependence. We show the spreading pressure dependence of adsorbate activity coefficient expression in the RAST framework can be generated from the gL isotherm. The procedure is illustrated with a spreading pressure dependent adsorption Nonrandom Furthermore, two-Liquid activity coefficient model, and the results are validated for ten binary mixed-gas adsorption equilibria systems including two highly nonideal azeotropic systems.

42 ENGINEERING↗
DOE OSTI2020

Influence of Proton Activity Gaps between Electrodes on Open-Circuit Potential of H 2 /H 2 and H 2 /Air Cells

Polymer -electrolyte-fuel-cell open-circuit voltages (OCVs) are exactly defined by equation (1), where cathode and anode proton activities [(aH+)cathode and (aH+)anode, respectively] usually are identical, so the third term in the right-hand side of equation (2) is ignored. OCV=E0+RT/2F*ln(a1/2 O2*(a2 H+)cathode/aH2O)-RT/2F*ln((a2 H+)anode/aH2) (1) =E0+RT/2F*ln(a1/2 O2*aH2/aH2O)+RT/2F*ln((a2 H+)cathode/(a2 H+)anode) (2) Water vapor pressure is a colligative property that fundamentally correlates to electrolyte concentrations in aqueous solutions. Proton activity is a function of acid concentration, such as pH, when electrolytes are acids. In polymer-electrolyte membranes, water vapor pressure and acid concentration are understood as relative humidity (RH) and water uptake (λ), respectively, where λ represents number of water molecules per sulfonic acid molecule. Several investigations have reported the relationship between RH and λ, meaning that proton activities and associated water uptakes are intimately related to RH. In actual fuel cell operation, cathode RH is determined by ambient-atmosphere and/or humidifier RH(s), and anode RH depends on hydrogen-circulator RH. Therefore, RH is not always identical at both electrodes, and the difference between electrode RHs is considerable during dry operation of polymer electrolyte fuel cells. Therefore, the third term in the right-hand side of equation (2) may be significant for dry operation. We measured OCVs when hydrogen was supplied to both electrodes at 80°C. One electrode (A) was fixed at 30% RH, while RH at the other electrode (B) was varied (0, 5, 10, 20, and 30%). Measured OCVs varied from 0 to 75 mV. For fuel cell tests, electrode A was supplied with hydrogen at 30% RH; electrode B, oxygen at 0, 5, 10, 20, and 30% RH. OCVs deviated from that measured when RH at electrode B was 30%, increasing from 0 to 60 mV with decreasing RH at electrode B. Results are also shown in Figure 1. Proton activities of both electrodes were thermodynamically calculated. The Gibbs–Duhem relation was applied to obtain molar Gibbs free energies of water and sulfonic acid, and proton activity coefficient was calculated using the Gibbs free energy of sulfonic acid and the relationship between RH and λ1–4, assuming that protons and sulfonic anions show identical ionic-activity coefficients. OCVs were estimated using the third term in the right-hand side of equation (2). Results are shown in Figure 1. Fuel-cell current–voltage performance was poor when RHs at the anode and cathode were 30 and 20%, respectively. To determine kinetic current, we measured the oxygen-reduction reaction (ORR) using a rotating-disk electrode (RDE) in concentrated-acid aqueous solutions, which modeled catalyst-layer ionomers. Kinetic currents decreased with acid concentrations. References T. A. Zawodzinski, Jr., C. Derouin, S. Radzinski, R. J. Sherman, V. T. Smith, T. E. Springer and S. Gottesfeld , J. Electrochem. Soc., 140,1041 (1993) P. K. Das and A. Z. Weber, Proceedings of the ASME 2013 11th Fuel Cell Science, Engineering and Technology Conference, Fuel Cell 18010 (2013) V. A. Sethuraman, J. W. Weidner, A. T. Haug, S. Motupally,b and L. V. Protsailo, J. Electrochem. Soc., 155, B50 (2008) A. Kusoglu and A. Z. Weber, Chem. Rev., 117, 987 (2017) Figure 1

Yoshida, Toshihiko↗
DOE OSTI2024

Generalized Brunauer–Emmett–Teller Isotherm for Mixed-Gas Multilayer Adsorption Equilibria

Here, this work presents a rigorous thermodynamic framework to calculate mixed-gas multilayer adsorption equilibria from single-component isotherms and vapor–liquid equilibria. Named the generalized Brunauer–Emmett–Teller (gBET) isotherm, the newly formulated isotherm considers adsorbent surface heterogeneity, competitive adsorption on the monolayer, condensation–evaporation on the subsequent layers, and adsorbed phase nonideality for the monolayer and the subsequent layers. The monolayer adsorbed phase nonideality is tracked by using an area-based adsorption nonrandom two-liquid activity coefficient model. The adsorbed phase composition and corresponding nonideality in the subsequent layers are calculated at the dew point condition of the mixed gas with either an equation of state or an activity coefficient model for the vapor–liquid equilibria. The proposed model is validated with six single and three binary multilayer adsorption equilibrium systems, and the model results are compared against those from the classical BET isotherm for single-component adsorption and ideal adsorbed solution theory for mixed-gas adsorption equilibria.

09 BIOMASS FUELS↗
DOE OSTI2023

Activity model for 36 elements in Fe-Ni-Si-S-C liquids with application to terrestrial planet accretion and mantle geochemistry: New data for Ru, Re, Pt, Os, Ti, Nb, and Ta

Understanding siderophile element partitioning between metal and silicate melts under diverse conditions can be used to place important constraints on the materials and conditions of planetary accretion and core formation, as well as post core formation processes. However, the effects of Si on the partitioning and activity coefficients for these elements are not well known, despite Si likely being one of the dominant light elements in Earth’s core. To address this gap in understanding, we have undertaken a systematic study of the highly siderophile elements Re, Pt, Os, and Ru, and the refractory lithophile elements Nb, Ta and Ti at 1600 °C and 1 GPa, to derive epsilon interaction parameters for these elements in FeSi metallic liquids. Positive epsilon interaction parameters were measured for Nb, Ta, Ti, Ru, Re, Pt, and Os, indicating that dissolved Si in Fe liquids causes a decrease in their metal/silicate partition coefficients (or ‘silicophobic’ behavior). Furthermore, ε$_{Re, Os, or Ru}^{Si}$ > ε$_{Re, Os, or Ru}^{S}$ which means Si causes a larger decrease in D(metal/silicate) than S, and the chalcophile behavior expected from some elements will be completely masked by the presence of Si in a metallic liquid. The new parameters are used to update an activity model that now includes 36 siderophile elements in Fe-Ni-Si-S-C liquids (27 trace elements considered here). Systematic assessment of these 27 elements shows which have the strongest affinity for Si, C, and S, and also how activity coefficients for these elements would vary during accretion and core formation in Earth, Mars, and Mercury of widely differing fO 2 and core compositional conditions. The activity model is combined with new partitioning expressions for Mo, W, Cr, Re Ru, Pt, and Os and applied to aspects of post core formation mantle geochemistry of Earth, Mars, and Mercury. Our updated expressions show that the BSE Mo/W ratio can easily be achieved with metal/silicate partitioning during growth of the Earth, whereas Re, Os and Ru become lower than and highly fractionated compared with BSE values during core formation and accretion, and thus nearly 99% of their BSE abundances are likely contributed by late accretion. Ru isotopes should be a very good indicator of the source material for the late accretion. The high Pt/Os and Re/Os developed in a deepening magma ocean during the growth of the Earth, indicates 186 Os and 187 Os isotopes could be coupled if this ancient material remained isolated and subsequently became entrained in mantle plumes and measured in surficial lavas. The extent to which this occurred will be limited by the low Os content of this ancient material, thus requiring mixing as a major component in plume sources. Martian mantle Hf/W ratio stays low during accretion and core formation modelling, suggesting that W isotope anomalies are more likely due to solid/liquid silicate fractionation than to core formation. Finally, Ti contents measured by MESSENGER at Mercury’s surface can be explained by segregation of either a metallic core (IW-6 to -8) or metallic core + sulfide (IW-4 to -7.5) followed by mantle melting.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2020

Equilibrium Thermodynamic Properties of Aqueous Solutions of Ionic Liquid 1-Ethyl-3-Methylimidazolium Methanesulfonate [EMIM][MeSO3]

Abstract The ionic liquid 1-ethyl-3-methylimidazolium methanesulfonate ([EMIM][MeSO 3 ]) has been considered as a promising alternative desiccant to triethylene glycol and lithium bromide commonly used in the industry. In this paper, the water activity coefficient of this binary system was measured from 303 K to 363 K with water concentration from 18% to 92%. The interaction energies between the ionic liquid molecules ( $${g}_{22}$$ g 22 ) and between the ionic liquid and water molecules ( $${g}_{12}$$ g 12 ) for the [EMIM][MeSO 3 ]/water binary system were determined from the water activity coefficient data using the Non-Random Two-Liquid (NRTL) model. The magnitude of the interaction energy between the [EMIM][MeSO 3 ] and water molecules ( $${g}_{12}$$ g 12 ) was found to be in the range of 45~49 kJ/mol, which was about 20% larger than that between the water molecules ( $${g}_{11}$$ g 11 ) in the [EMIM][MeSO 3 ]/water system. The large ( $${g}_{12}$$ g 12 ) can explain many observed macroscopic thermodynamic properties such as strong hygroscopicity in the ionic liquid [EMIM][MeSO 3 ]. These interaction energies were used to determine the heat of desorption of the [EMIM][MeSO 3 ]/water system, and the obtained heat of desorption was in good agreement with that calculated from the conventional Clausius-Clapeyron Equation.

42 ENGINEERING↗
DOE OSTI2022

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↗
DOE OSTI2021

Experimental investigation of elemental and isotopic evaporation processes by laser heating in an aerodynamic levitation furnace

In this report we carried out evaporation experiments on a B-type calcium–aluminium-rich inclusion (CAI) melt in a gas-mixing aerodynamic levitation laser furnace, at 1873 K and an oxygen partial pressure of 10 -9.1 atm, for durations ranging from 60 to 600 s. Evaporation of SiO 2 and MgO follow the same trend as those observed in vacuum furnace experiments at the same temperature and starting composition, showing that their evaporation relative to one another from the melt is independent of pressure, oxygen fugacity, and hydrodynamical regime specific to the furnace. Isotopic ratios of Mg and Si in evaporation residues are used to derive fractionation factors of α 26/24 Mg vap-liq = 0.9906± 0.0004 and α 30/28 Si vap-lip = 0.9943±0.0003, which are both significantly closer to unity than those found for evaporation in a vacuum, which translates to less isotope fractionation. The residues are also less isotopically fractionated than expected for cases in which transport of the gas species away from the melt is diffusion-controlled at 1-atm. By analysing the flow regimes in our furnace, we find that advection by the levitating gas is the primary mode of mass transport away from the melt surface, as opposed to diffusion-limited transport in a vacuum or 1-atm tube furnace. A modified Hertz–Knudsen–Langmuir formulation accounts for this process, and shows that isotopic fractionation of both Si and Mg reflect a saturation factor (ratio of the pressure of the evaporating species to vapour saturation pressure) equal to 0.75. This is in perfect accord with recent measurements of Cu isotopic fractionation using a similar furnace. The fact that three elements (Mg, Si, Cu) with varying equilibrium vapour pressures, activity coefficients in the liquid, and diffusion coefficients in the gas have the same scaling behaviour to saturation pressure is a strong indication that the mechanism controlling evaporation is driven by the hydrodynamical regime imposed in the furnace. Therefore, this class of experiments can be used to constrain processes in which advection dominates over diffusion, such as (but not limited to) planetary ejecta, tektites, giant impacts, nebular condensation in a turbulent flow, or nuclear fallout material. Finally, the possibility to reach high temperatures (in excess of 3500 K) in this furnace allows it to be used to evaluate the activity coefficients of melt components in extreme conditions relevant to molten planetary interiors (i.e., magma oceans), with a specific focus on refractory elements.

58 GEOSCIENCES↗
DOE OSTI2020

Simulations of activities, solubilities, transport properties, and nucleation rates for aqueous electrolyte solutions

This article reviews recent molecular simulation studies of "collective" properties of aqueous electrolyte solutions, specifically free energies and activity coefficients, solubilities, nucleation rates of crystals, and transport coefficients. These are important fundamental properties for biology and geoscience, but also relevant for many technological applications. Their determination from molecular-scale calculations requires large systems and long sampling times, as well as specialized sampling algorithms. As a result, such properties have not typically been taken into account during optimization of force field parameters; thus, they provide stringent tests for the transferability and range of applicability of proposed molecular models. There has been significant progress on simulation algorithms to enable the determination of these properties with good statistical uncertainties. Comparisons of simulation results to experimental data reveal deficiencies shared by many commonly used models. Moreover, there appear to exist specific tradeoffs within existing modeling frameworks, so that good prediction of some properties is linked to poor prediction for specific other properties. For example, non-polarizable models that utilize full charges on the ions generally fail to predict accurately both activity coefficients and solubilities; the concentration dependence of viscosity and diffusivity for these models is also incorrect. Scaled-charge models improve the dynamic properties and could also perform well for solubilities, but fail in the prediction of nucleation rates. Even models that do well at room temperature for some properties generally fail to capture their experimentally observed temperature dependence. Finally, the main conclusion from the present review is that qualitatively new physics will need to be incorporated in future models of electrolyte solutions to allow description of collective properties for broad ranges of concentrations, temperatures, and solvent conditions.

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

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↗
DOE OSTI2023

Effect of Na Concentration on Cs Distribution with Crystalline Silicotitanate in Tank Waste Simulants

Crystalline silicotitanate (CST) ion exchange media is currently utilized in the tank side cesium removal (TSCR) system on the Hanford site to remove Cs-137 from the tank waste supernate. As the main dose contributor to the liquid supernate in Hanford tank waste, it is necessary to remove the Cs-137 to expedite low activity waste processing and immobilization. Ongoing batch contact experiments with CST have been done to tease out the impact of group II metals, anions, potassium, and temperature on Cs removal to better understand bounding conditions for operations as well as aid in development of a wholesome isotherm model to predict Cs behavior in the tank waste. This paper discusses the impact of Na concentration on Cs exchange and provides insight into the challenges associated with activity coefficient estimations for modeling Cs distribution in tank waste matrices.

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

Modeling swelling behavior of hydrogels in aqueous organic solvents

Hydrogel, a three-dimensional elastic structure of cross-linked polymer chains, swells or deswells when immersed in aqueous organic solutions. The swelling of hydrogels depends on temperature, solvent composition, and hydrogel structure as the thermodynamic driving forces for hydrogel swelling are due to solvent – hydrogel interactions and expansion of elastic network. Here in this work, polymer Nonrandom Two-Liquid activity coefficient model is used to account for the solvent – hydrogel interactions and a semi-empirical expression for rubber elasticity is used to describe the expansion of elastic network. The model successfully correlates data for both swelling and equilibrium phase compositions for n-isopropyl acrylamide hydrogels in water – ethanol, water – acetone, water – n-butanol, and water – methyl isobutyl ketone mixed solvents.

42 ENGINEERING↗
DOE OSTI2025

Thermodynamic Water Activity Explains the Unusual Electrochemical Stability of Aqueous Deep Eutectic Solvents

The presence of water in nonaqueous deep eutectic solvent (DES) electrolytes has been debated in recent years, with efforts ranging from its complete removal to willful addition. It was shown that controlled amounts of water can be beneficial, as it not only enhances the physicochemical properties of these electrolytes but also has no significant detrimental effects on their electrochemical stability. Despite these advantages, there is still limited understanding of how water interacts with DES systems at the molecular level. This study examines the water activity in ethylene glycol and glycerol, as well as their binary mixtures with choline chloride to form the DESs ethaline and glyceline, respectively. In this work, we show that the high electrochemical stability of glyceline is related to its lower water activity compared to ethaline and can be attributed to the robust H-bonding network formed by the three hydroxyl groups of glycerol. Its 3D H-bond network effectively integrates water molecules within its solvent structure, reducing degradation and maintaining stability at higher water contents. The deviations from the ideal Raoult's law behavior are reflected in the water activity and activity coefficients, which highlight the intricate H-bond interactions within DES-water mixtures. Water acts like a lubricant within the more viscous DES mixtures without being detrimental to their electrochemical performance. The presented results emphasize the necessity of customizing DES-water compositions to enhance their performance as electrolytes, especially in flow battery applications where electrochemical stability, ionic conductivity, and fluidity are of utmost importance.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗