Search NASA⌕ Search

SEARCH · Search NASA

Results for “predictive thermodynamics”

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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Henry’s Solubility and Diffusion Coefficients for 29 Volatile Organic Compounds in Polydimethylsiloxane Sylgard 184 at 293 K

Two-dimensional (2D) inverse gas chromatography (IGC) enables simultaneous determination of Henry’s solubility and Fickian diffusion coefficients for volatile organic compounds (VOCs) in polymer films. This technique offers a significant advantage over traditional cylindrical column IGC by providing precise control and measurement of the film thickness (here, 0.064 ± 0.002 mm), which is the critical length scale for accurate diffusivity determination. We apply this methodology to characterize VOC transport in Sylgard 184, a widely used polydimethylsiloxane (PDMS)-based polymer containing substantial silica filler content. At room temperature (20 °C), we measured solubility and diffusion coefficients for 29 common VOCs spanning diverse chemical functionalities, including alkanes, aromatics, chlorinated solvents, ketones, esters, and alcohols. Comparison with literature data for pure PDMS reveals that VOC solubility in Sylgard 184 is generally higher; for most non-hydrogen-bonding compounds it remains within a factor of 2 of pure PDMS, whereas alcohols are enhanced by roughly 1.8 to 3.7 times, which we attribute to favorable interactions with residual silanol groups on the silanized silica filler. Diffusion coefficients range from 1.0 × 10 –6 cm 2 /s (n-undecane) to 8.9 × 10–5 cm 2 /s (acetonitrile) and align well with extrapolated literature values for PDMS. This study provides essential thermodynamic and transport data for predicting VOC permeation in Sylgard 184 while demonstrating the utility of 2D IGC as a robust technique for characterizing rubbery polymer membranes across diverse industrial applications.

organic↗
DOE OSTI2022

Prediction of Above-Room-Temperature Superconductivity in Lanthanide/Actinide Extreme Superhydrides

Achieving superconductivity at or above room temperature has been a century long held dream for physicists since the discovery of superconductivity in mercury in 1911. Following the recent predictions and ensuing synthesis of clathrate superhydride LaH 10 under pressure exhibiting extraordinary superconducting critical temperatures (T c ) of 250 260 K, we predict via advanced crystal structure search methods a new class of extremely hydrogen rich clathrate superhydrides. These MH 18 (M: rare earth/actinide metal atom) stoichiometric compounds consisting of H36 cage networks are predicted to host T c values above room temperature up to 330 K at pressures of 350 GPa. The bonding and electronic properties of these MH 18 clathrate superhydrides parallel those of atomic metallic hydrogen, giving rise to the highest superconducting temperatures predicted thus far for a thermodynamically stable hydride compound. In depth examination of these extreme superhydrides offers key insights for elucidating and further exploring phonon mediated superconductivity above room temperature in hydrogen rich and other low Z materials.

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

Thermodynamic insight into AnO 2 + bonding from ThH + /UH + reactions studied by inductively coupled plasma tandem mass spectrometry

The bonding mechanisms of actinides have been a focus of fundamental research over the past few decades. In the present study, reactions of the simplest actinide-containing species, ThH + and UH + , with O 2 and CO 2 are investigated by inductively coupled plasma tandem mass spectrometry. The reactions of ThH + and UH + with O 2 are efficient, and the reactions of ThH + and UH + with CO 2 display reduced reaction efficiencies. For both reactions involving CO 2 , ThO 2 + and UO 2 + are observed; however, there is a clear barrier to ThO 2 + formation whereas UO 2 + forms through an exothermic, barrierless process. The experimental observations and available thermodynamic information are used to predict the outcomes of reactions involving the later AnH + . The anticipated reaction enthalpies for Pa–Am display a clear correlation with the promotion energy of An + to a 6d 2 electronic configuration, E p (6d 2 ), although a shift in the slope of the correlation of reaction enthalpies and E p (6d 2 ) suggests that there is likely a change in bonding mechanism that starts with Np + . Similar shifts have also been noted in previous studies. Beginning with Np + , the 6d orbitals become less accessible than they are for the earlier An + as measured by E p (6d 2 ), and this accessibility of the 6d orbitals may drive actinide bonding.

5f orbitals↗
DOE OSTI2020

The role of chemical potential in compensation control in Si:AlGaN

Reduction in compensation in Si-doped Al-rich AlGaN is demonstrated via chemical potential control (CPC). The chemical potentials and the resulting formation energies of carbon on the nitrogen site (C N ) and cation vacancy complex with Si (V III + nSi III ) were related to growth variables through a thermodynamic supersaturation model, which quantitatively predicted the incorporation of C N and the generation of the V III + nSi III complex. The compensation “knee” behavior, i.e., decreasing conductivity with increasing Si incorporation beyond a certain concentration, was successfully controlled. The maximum free carrier concentration was improved by impeding the formation of V III + nSIi III complexes under III-richer conditions, while the impurity compensation by C N was reduced by making the growth environment N-richer. The results of Hall effect measurement and photoluminescence agreed well with quantitative theoretical predictions of the CPC model. Based on the developed model, the highest conductivity of 160 Ω -1 cm -1 with free carrier concentration of 3 × 10 19 cm -3 in Al 0.7 Ga 0.3 N ever reported was achieved on single crystal AlN substrates. The demonstrated predictive power of the CPC model should greatly reduce the empirical analysis or iterative experimentation that would otherwise be necessary.

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

Electrical control of two-dimensional electron-hole fluids in the quantum Hall regime

We study the influence of quantizing perpendicular magnetic fields on the ground state of a bilayer with electron and hole fluids separated by an opaque tunnel barrier. In the absence of a field, the ground state at low carrier densities is a condensate of s-wave excitons that has spontaneous interlayer phase coherence. We find that a series of phase transitions emerge at strong perpendicular fields between condensed states and incompressible incoherent states with full electron and hole Landau levels. When the electron and hole densities are unequal, condensation can occur in higher angular momentum electron-hole pair states and, at weak fields, break rotational symmetry. Here, we explain how this physics is expressed in dual-gate phase diagrams, and predict transport and capacitively probed thermodynamic signatures that distinguish different states.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Property Measurements of NaCl-UCl 3 and NaCl-KCl-UCl 3 Molten Salts (Rev.1)

Thermochemical and thermophysical property values of several salt compositions of interest are needed by molten salt reactor (MSR) developers to design, license, and operate their reactors. Thermochemical and thermophysical properties being measured at Argonne include thermal transitions, phase behavior, heat capacity, density, volumetric thermal expansion of the liquid phase, thermal diffusivity, thermal conductivity, and viscosity. Several properties of eutectic compositions in the ternary NaCl-KCl-UCl 3 and binary NaCl-UCl 3 systems that may be used by MSR developers as fuel bearing salts are being measured. A 65.8 mol % NaCl–34.2 mol % UCl 3 mixture and a near-eutectic mixture of 50.9 mol % NaCl–24.4 mol % KCl–24.7 mol % UCl 3 were synthesized and the thermochemical properties of the mixtures were measured by using differential scanning calorimetry (DSC). The measured transition temperatures were compared to transition temperatures predicted using two models. A thermodynamic model of the binary NaCl-UCl 3 system was constructed using data in the Molten Salt Thermal Properties Database–Thermochemical Version 2.0 (MSTDB-TC V2.0). A ternary NaClKCl-UCl 3 system model constructed at Argonne and was described in a previous report. These comparisons can be used to validate the models. Thermophysical property values of molten salts are needed to model how salt retains and transfers heat in an MSR system. These property values are essential to the entire MSR design because molten salt is used as both the fuel and the coolant material in a salt fueled reactor. Heat capacity of the synthesized NaCl-UCl 3 and NaCl-KCl-UCl 3 salt mixtures was measured by using DSC and thermal diffusivity was measured by using laser flash analysis (LFA) at temperatures spanning the typical operating range of an MSR.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

Quantifying Uncertainties in Earth's Energy Budget by Cloud Feedback and Ocean Heat Uptake Using E3SM-Slab Ocean Configurations

In order to improve predictive skills of Earth System Model, we need to better understand processes that control Earth's energy budget via ocean, atmosphere, and cryosphere interactions. Simulated energy budget in comprehensive Earth System Models shows a wide range, leading to large uncertainties in predicting Earth system dynamic and thermodynamic variations and associated social-economic impacts. Uncertainties in cloud feedbacks have been identified as the main cause of the large inter-model spread, but oceanic adjustments, especially those associated with ocean heat uptake (OHU) and the Atlantic Meridional Overturning Circulation (AMOC), also play an important role. In this proposed work, we focus on understanding the individual and combined roles of cloud feedbacks and ocean adjustments on modulating Earth's energy balance. This research is motivated by our overarching hypothesis that oceanic adjustment is a key source of uncertainty, in addition to those associated with the cloud feedbacks; further, the ocean adjustment and associated OHU work through the cloud feedbacks to modulate Earth's energy budget and temperature variations. We test this hypothesis using numerical experiments where we systematically enable and disable cloud feedbacks in conjunction with perturbations to OHU.

58 GEOSCIENCES↗
DOE OSTI2026

Organic Rankine Cycle Integration and Optimization for High Efficiency CHP Genset Systems (Final Technical Report)

This project successfully advanced the integration of Organic Rankine Cycle (ORC) technology with reciprocating engine–based combined heat and power (CHP) systems to improve electrical efficiency, total CHP efficiency, and grid-responsive operation. Over three budget periods, the work progressed from high-temperature ORC component development and thermodynamic model validation to next-generation system design, working fluid transition, and techno-economic analysis. Key technical accomplishments include development and validation of a thermodynamic model capable of accurately predicting ORC performance across an expanded temperature and pressure envelope; successful identification and validation of low-global-warming-potential (GWP) working fluids—most notably R1233zd(E)—as viable replacements for R245fa; and demonstration of scalable ORC architectures suitable for integration with 1–20 MW class reciprocating engines. These advances enable flexible CHP configurations that can increase electrical output while maintaining high overall utilization of available thermal energy. The project also produced a clean-sheet design for a next-generation ORC system targeting substantially higher power output per unit, supported by detailed component selection, heat exchanger evaluation, and system-level modeling. Techno-economic analyses indicate that ORC-enabled flexible CHP systems can meet or exceed Department of Energy (DOE) efficiency targets while providing value to both facility operators and the electric grid.. Late-stage testing of the largest next-generation ORC prototype identified limitations related to pump net positive suction head (NPSH) requirements and condenser flooding under certain operating conditions. Although these issues constrained full validation of that configuration within the project timeframe, they provided clear and actionable design guidance for future system refinements. Importantly, validated modeling, smaller-scale testing, and working fluid evaluations confirmed the technical viability of the overall approach. In aggregate, this project met its core objectives by establishing validated design tools, de-risking key ORC technologies for CHP applications, and defining a credible pathway toward commercialization of flexible, high-efficiency CHP systems. The results form a strong foundation for continued development and deployment beyond the conclusion of the DOE-funded effort.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2021

Fine Aerosol Acidity and Water during Summer in the Eastern North Atlantic

Aerosol pH governs many important atmospheric processes that occur in the marine boundary layer, including regulating halogen and sulfur chemistries, and nutrient fertilization of surface ocean waters. In this study, we investigated the acidity of PM1 over the eastern North Atlantic during the Aerosol and Cloud Experiments in Eastern North Atlantic (ACE-ENA) aircraft campaign. The ISORROPIA-II thermodynamic model was used to predict PM1 pH and water. We first investigated the sensitivities of PM1 pH and water predictions to gas-phase NH3 and HNO3 concentrations. Our sensitivity analysis indicated that even though NH3 and HNO3 were present at very low concentrations in the eastern North Atlantic during the campaign, PM1 pH calculations can still be sensitive to NH3 concentrations. Specifically, NH3 was needed to constrain the pH of populations of PM1 that had low mass concentrations of NH4+ and non-volatile cations (NVCs). We next assumed that gas-phase NH3 and HNO3 concentrations during the campaign were 0.15 and 0.09 µg m−3, respectively, based on previous measurements conducted in the eastern North Atlantic. Using the assumption that PM1 were internally mixed (i.e., bulk PM1), we determined that PM1 pH ranged from 0.3–8.6, with a mean pH of 5.0 ± 2.3. The pH depended on both Hair+ and Wi. Hair+ was controlled primarily by the NVCs/SO42− molar ratio, while Wi was controlled by the SO42− mass concentration and RH. Changes in pH with altitude were driven primarily by changes in SO42−. Since aerosols in marine atmospheres are rarely internally mixed, the scenario where non-sea salt species and sea-salt species were present in two separate aerosol modes in the PM1 (i.e., completely externally mixed) was also considered. Smaller pH values were predicted for the aerosol mode comprised only of non-sea salt species compared to the bulk PM1 (difference of around 1 unit on average). This was due to the exclusion of sea-salt species (especially hygroscopic alkaline NVCs) in this aerosol mode, which led to increases in Hair+ values and decreases in Wi values. This result demonstrated that assumptions of aerosol mixing states can impact aerosol pH predictions substantially, which will have important implications for evaluating the nature and magnitude of pH-dependent atmospheric processes that occur in the marine boundary layer.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Primary biomolecular adsorption energetics of core–shell nanocomplexes: Implications for biological interactions

Molecular organization at the nano-bio interface governing the colloidal stability, reactivity, immune recognition, and drug delivery performance of nanoparticles remains difficult to predict. Quantifying the primary hydration energetics of biomolecule-coated nanomaterials can determine those interactions and provide a basis for engineered nanocarriers with tailored behavior in biological systems. Here, we measured the thermodynamics of water adsorption on patchy dry magnetite (Fe 3 O 4 ) nanoparticles coated with three model biomolecules, bovine serum albumin, potato starch, and lauric acid and compared these properties to the hydration energetics of the corresponding free dry biomolecules. The results demonstrate how the surface functionalization alters the hydrophilicity, the accessible hydrophilic surface, and the interaction potential of the nanocomplex surface with biological media. The protein coating increases the interaction potential of the surface of the nanocomplex. The weaker interaction potential of the polysaccharide coating and the relatively large hydrophilic surface area allow dynamic and reversible binding, while the fatty acid rearranges into a partial bilayer with very strong hydrophilicity. The findings establish the hydration enthalpy as a quantitative basis to determine and interpret nanoparticle interactions with proteins, membranes, and biological fluids, and provide a thermodynamic foundation for designing nanocarriers with predictable biological reactivity.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Uncertainty analysis of correlated parameters in automated reaction mechanism generation

Abstract Uncertainty analysis is a useful tool for inspecting and improving detailed kinetic mechanisms because it can identify the greatest sources of model output error. Owing to the very nonlinear relationship between kinetic and thermodynamic parameters and computed concentrations, model predictions can be extremely sensitive to uncertainties in some parameters while uncertainties in other parameters can be irrelevant. Error propagation becomes even more convoluted in automatically generated kinetic models, where input uncertainties are correlated through kinetic rate rules and thermodynamic group values. Local and global uncertainty analyses were implemented and used to analyze error propagation in Reaction Mechanism Generator (RMG), an open‐source software for generating kinetic models. A framework for automatically assigning parameter uncertainties to estimated thermodynamics and kinetics was created, enabling tracking of correlated uncertainties. Local first‐order uncertainty propagation was implemented using sensitivities computed natively within RMG. Global uncertainty analysis was implemented using adaptive Smolyak pseudospectral approximations as implemented in the MIT Uncertainty Quantification Library to efficiently compute and construct polynomial chaos expansions to approximate the dependence of outputs on a subset of uncertain inputs. Cantera was used as a backend for simulating the reactor system in the global analysis. Analyses were performed for a phenyldodecane pyrolysis model. Local and global methods demonstrated similar trends; however, many uncertainties were significantly overestimated by the local analysis. Both local and global analyses show that correlated uncertainties based on kinetic rate rules and thermochemical groups drastically reduce a model's degrees of freedom and have a large impact on the determination of the most influential input parameters. These results highlight the necessity of incorporating uncertainty analysis in the mechanism generation workflow.

Gao, Connie W.↗
DOE OSTI2025

Erratum: Review—Materials Science Predictions of Thermal Runaway in Layered Metal-Oxide Cathodes: A Review of Thermodynamics [ J. Electrochem. Soc. , 167, 090543 (2020)]

It was found that corrections are required for the published article cited in the title of this document. A revised Table VII is shown below, in which the values for S (entropy of the liquid solvents) have been corrected for consistency with the original sources. No corrections were required for the enthalpies, but they are included for completeness. The cited references and table footnotes found in the originally published version of Table VII remain applicable.

Shurtz, Randy C. [Sandia National Laboratories (SN↗
DOE OSTI2022

Correlating Thermodynamic and Kinetic Hydricities of Rhenium Hydrides

The kinetics of hydride transfer from Re( R bpy)(CO) 3 H (bpy = 4,4'-R-2,2'-bipyridine; R = OMe, t Bu, Me, H, Br, COOMe, CF 3 ) to CO 2 and seven different cationic N-heterocycles were determined. Additionally, the thermodynamic hydricities of complexes of the type Re( R bpy)(CO) 3 H were established primarily using computational methods. Linear free-energy relationships (LFERs) derived by correlating thermodynamic and kinetic hydricities indicate that, in general, the rate of hydride transfer increases as the thermodynamic driving force for the reaction increases. Kinetic isotope effects range from inverse for hydride transfer reactions with a small driving force to normal for reactions with a large driving force. Hammett analysis indicates that hydride transfer reactions with greater thermodynamic driving force are less sensitive to changes in the electronic properties of the metal hydride, presumably because there is less buildup of charge in the increasingly early transition state. Bronsted α values were obtained for a range of hydride transfer reactions and along with DFT calculations suggest the reactions are concerted, which enables the use of Marcus theory to analyze hydride transfer reactions involving transition metal hydrides. It is notable, however, that even slight perturbations in the steric properties of the Re hydride or the hydride acceptor result in large deviations in the predicted rate of hydride transfer based on thermodynamic driving forces. This indicates that thermodynamic considerations alone cannot be used to predict the rate of hydride transfer, which has implications for catalyst design.

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

Proton-coupled electron transfer across benzimidazole bridges in bioinspired proton wires

Designing molecular platforms for controlling proton and electron movement in artificial photosynthetic systems is crucial to efficient catalysis and solar energy conversion. The transfer of both protons and electrons during a reaction is known as proton-coupled electron transfer (PCET) and is used by nature in myriad ways to provide low overpotential pathways for redox reactions and redox leveling, as well as to generate bioenergetic proton currents. Herein, we describe theoretical and electrochemical studies of a series of bioinspired benzimidazole-phenol (BIP) derivatives and a series of dibenzimidazole-phenol (BI 2 P) analogs with each series bearing the same set of terminal proton-accepting (TPA) groups. The set of TPAs spans more than 6 pK a units. These compounds have been designed to explore the role of the bridging benzimidazole(s) in a one-electron oxidation process coupled to intramolecular proton translocation across either two (the BIP series) or three (the BI 2 P series) acid/base sites. These molecular constructs feature an electrochemically active phenol connected to the TPA group through a benzimidazole-based bridge, which together with the phenol and TPA group form a covalent framework supporting a Grotthuss-type hydrogen-bonded network. Infrared spectroelectrochemistry demonstrates that upon oxidation of the phenol, protons translocate across this well-defined hydrogen-bonded network to a TPA group. The experimental data show the benzimidazole bridges are non-innocent participants in the PCET process in that the addition of each benzimidazole unit lowers the redox potential of the phenoxyl radical/phenol couple by 60 mV, regardless of the nature of the TPA group. Using a series of hypothetical thermodynamic steps, density functional theory calculations correctly predicted the dependence of the redox potential of the phenoxyl radical/phenol couple on the nature of the final protonated species and provided insight into the thermodynamic role of dibenzimidazole units in the PCET process. This information is crucial for developing molecular “dry proton wires” with these moieties, which can transfer protons via a Grotthuss-type mechanism over long distances without the intervention of water molecules.

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

Optimized thermodynamic properties of REE aqueous species (REE 3+ and REEOH 2+ ) and experimental database for modeling the solubility of REE phosphate minerals (monazite, xenotime, and rhabdophane) from 25 to 300 °C

Rare earth elements (REE) are critical elements found in monazite, xenotime, and hydrated REE phosphates which typically form in hydrothermal mineral deposits. Accurate predictions of the solubility of these REE phosphates and the speciation of REE in aqueous fluids are both key to understanding the controls on the transport, fractionation, and deposition of REE in natural systems. Previous monazite and xenotime solubility experiments indicate the presence of large discrepancies between experimentally derived solubility constants versus calculated solubilities by combining different data sources for the thermodynamic properties of minerals and aqueous species at hydrothermal conditions. In this study, these discrepancies were resolved by using the program GEMSFITS to optimize the standard partial molal Gibbs energy of formation (Δ f G° 298 ) of REE aqueous species (REE 3+ and REE hydroxyl complexes) at 298.15 K and 1 bar while keeping the thermodynamic properties fixed for the REE phosphates. A comprehensive experimental database was compiled using solubility data available between 25 and 300 °C. The latter permits conducting thermodynamic parameter optimization of Δ f G° 298 for REE aqueous species. Optimal matching of the rhabdophane solubility data between 25 and 100 °C requires modifying the Δ f G° 298 values of REE 3+ by 1–6 kJ/mol, whereas matching of the monazite solubility data between 100 and 300 °C requires modifying the Δ f G° 298 values of both REE 3+ and REEOH 2+ by ~15–31 kJ/mol and ~2–10 kJ/mol, respectively. For xenotime, adjustments of Δ f G° 298 values by 1–26 kJ/mol are only necessary for the REE 3+ species. The optimizations indicate that the solubility of monazite in acidic solutions is controlled by the light (L)REE 3+ species at <150 °C and the LREEOH 2+ species at >150 °C, whereas the solubility of xenotime is controlled by the heavy (H)REE 3+ species between 25 and 300 °C. Based on the optimization results, we conclude that the revised Helgeson-Kirkham-Flowers equation of state does not reliably predict the thermodynamic properties of REE 3+ , REEOH 2+ , and likely other REE hydroxyl species at hydrothermal conditions. We therefore provide an experimental database (ThermoExp_REE) as a basic framework for future updates, extensions with other ligands, and optimizations as new experimental REE data become available. As a result, the optimized thermodynamic properties of aqueous species and minerals are available open access to accurately predict the solubility of REE phosphates in fluid-rock systems.

58 GEOSCIENCES↗
DOE OSTI2024

Thermodynamic driving forces in contact electrification between polymeric materials

Abstract Contact electrification, or contact charging, refers to the process of static charge accumulation after rubbing, or even simple touching, of two materials. Despite its relevance in static electricity, various natural phenomena, and numerous technologies, contact charging remains poorly understood. For insulating materials, even the species of charge carrier may be unknown, and the direction of charge-transfer lacks firm molecular-level explanation. Here, we use all-atom molecular dynamics simulations to investigate whether thermodynamics can explain contact charging between insulating polymers. Based on prior work suggesting that water-ions, such as hydronium and hydroxide ions, are potential charge carriers, we predict preferred directions of charge-transfer between polymer surfaces according to the free energy of water-ions within water droplets on such surfaces. Broad agreement between our predictions and experimental triboelectric series indicate that thermodynamically driven ion-transfer likely influences contact charging of polymers. Furthermore, simulation analyses reveal how specific interactions of water and water-ions proximate to the polymer-water interface explain observed trends. This study establishes relevance of thermodynamic driving forces in contact charging of insulators with new evidence informed by molecular-level interactions. These insights have direct implications for future mechanistic studies and applications of contact charging involving polymeric materials.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Evaluation of New Additions to OLI Software in Predicting Mercuric and Mercurous Species in Liquid Waste Operations

Speciation of mercury during the pretreatment steps of tank waste processing is critical to successful mercury removal prior to vitrification during Liquid Waste Operations (LWO) at SRS. OLI software has been used to predict mercury speciation and activity throughout LWO. The OLI software operates based on a thermodynamic framework called the Mixed Solvent Electrolyte (MSE) framework. The MSE framework allows prediction in theoretically infinitely dilute to concentrated mixtures (e.g., purely solute solutions). Before modification to the MSE framework databanks, certain critical mercury species were missing in the MSE databank, and some thermodynamic data needed to be updated for the OLI software to accurately predict mercury chemical species in SRS waste tanks. To better reflect streams across LWO, new mercury species were integrated into the MSE database. To evaluate the changes to the OLI MSE framework per the Technical Task Request (TTR) and the Task Technical and Quality Assurance Plan (TTQAP), waste stream compositions from Tanks 38, 43, and Tank 50 decontaminated salt solution (DSS) were used as model inputs. Models were developed and executed using both the old and new databases. Compositional analyses from caustic Tank 50 DSS and caustic Tanks 38 and 43 were used as the input streams. These streams represent the most comprehensive chemical data sets where both mercury and tank constituents were measured together. Results for Tank 50 DSS predict HgO as the predominant species in both databases. Both methyl and dimethyl Hg species are present when the new database is ‘on’ and are not predicted with the new database turned ‘off’. The new database predicts a greater amount of HgO and a greater fraction of it in the solid phase. Pourbaix diagrams (potential vs. pH) generated for each Tank 50 DSS were identical regardless of which database was used. Elemental Hg and HgO were predicted in the water stable region under basic conditions. Tanks 38 and 43 follow similar trends as the Tank 50 DSS models. Unlike Tanks 38 and 50 DSS, the Tank 43 Pourbaix plot shows a region of stability for an aqueous HgOHCO3 - species between approximately pH 7-11. In all streams, when MeHg+ is included in the inputs, the new database predicts aqueous MeHgOH as the dominant species. If elemental or dimethyl mercury is in the waste stream, the new database model predicts they are unchanged and remain in those states and quantities. Additionally, the total mercury values are reported for both the measured input data and the OLI output data for all considered tanks. The summary indicates that the percentage error between the measured and calculated values is less than 1% in all cases The reconciliations and generation of the Pourbaix diagrams for Tank 50 DSS took approximately ten times longer with the new database ‘on’. In addition, over the course of that time, models with the new database ‘on’ were more likely to crash or display an error. Some modest performance improvements were noted when modeling with an i7 processor versus an i5. An example error is found in Appendix A. Furthermore, Appendix B provides V&V for two chemical systems analyzed with the OLI software, results were satisfactory. It is recommended to utilize the new databases (i.e., HCO.ddb and SR-Hg.ddb) in future Savannah River Mission Completion applications of OLI to represent pseudo steady-state. Furthermore, the integration and utilization of the new databases (i.e., HCO.ddb and SR-Hg.ddb) in modeling applications (e.g., Aspen) is also recommended.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗