Search NASA⌕ Search

SEARCH · Search NASA

Results for “Thermodynamic Model”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Degradation mechanism of blended cement pastes in sulfate-bearing environments under applied electric fields: Sulfate attack vs. decalcification

Applied electric fields, the reason behind stray currents, accelerate the ingress of sulfate into cementitious materials. To identify mitigation approaches, this study investigates the effects of slag and fly ash on sodium sulfate attacks of cement pastes under a constant electric current. The mineralogical alterations induced by the attacks were analyzed using X-ray diffractometry, thermogravimetric analysis, scanning electron microscopy, and thermodynamic modeling. The dissolution of aluminates in the slag and fly ash induced a monosulfate-rich area (>~20 mm from the cathode surface) to form next to the ettringite-rich area on the sample surface (<~20 mm). This effect reduced the availability of SO 4 2– in the pore solution, thereby hindering the penetration of sulfate. Meanwhile, the consumption of portlandite by the pozzolanic reaction lowered the decalcification resistance of the materials. This produced a wide area that endured the decomposition of portlandite and carbonate-AFm. Altogether, blending 30% fly ash did not improve the resistance of the material to either sulfate attacks or decalcification; blending 50% slag can effectively mitigate sulfate ingress, though decalcification may become a governing mechanism in the degradation process.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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 OSTI2020

Performance of sweeping gas membrane distillation for treating produced water: Modeling and experiments

Membrane distillation (MD) has shown promise for desalination of produced water (PW), but there are very few studies reported on the performance of MD with real PW and the characteristics of the fouling and scaling in MD. Among different MD configurations, sweeping gas MD (SGMD) is the least studied configuration. Here we report the application of SGMD for treating high-salinity PW (123 g L -1 ). The SGMD operating conditions were first carefully optimized through pure water MD experiments by varying operating conditions. The flux evaluation was accompanied by a detailed Aspen Custom Modeler simulation based on the electrolyte-NRTL thermodynamic model to better understand the heat and mass transfer in the feed and permeate channel. Supported by experiments, modeling results confirmed the partial condensation of the vapor in the permeate channel. Furthermore, it was concluded that the saturation of the sweeping gas is controlling the overall mass transfer through the membrane. Selected membranes and optimized conditions were then used to conduct short-term and long-term performance tests with PW. In long-term MD test, the concentration of the PW increased by a factor of two (50% water recovery), while the membrane exhibited small flux decline. Elemental analysis showed that strontium and sodium were the major scalant compounds.

42 ENGINEERING↗
DOE OSTI2021

Modeling Phase Equilibrium of Common Sugars Glucose, Fructose, and Sucrose in Mixed Solvents

The industrial processing of sugars and sugar-containing mixtures is gaining widespread use as a form of renewable manufacturing from biomass. To aid in process modeling for design and optimization, a commonly available thermodynamic model is needed that describes the phase equilibrium of these compounds. This work compiles and compares models for solid–liquid and vapor–liquid phase equilibrium from the available data for the representative sugars glucose, fructose, and sucrose in the representative solvents water, methanol, and ethanol, including data for multisugar, multisolvent systems. Additionally, the nonrandom two-liquid (NRTL) model was chosen for these systems because of its widespread use in industry and the availability of parameters for many solvents and cosolutes. The association-NRTL (aNRTL) model was investigated as an improvement for modeling sugars, which may experience a high degree of association because of their many hydroxy groups. Both models accurately capture the data and are able to predict the behavior of multisugar systems with only solute–solvent interaction parameters. The aNRTL model shows an improvement over the baseline NRTL model that is most significant for sucrose, the component with the highest association strength, and least significant for fructose, which has the lowest association strength.

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

First-Principles Thermodynamic Assessments of Sr-Containing Secondary Phase Formation in La1-xSrxMnO3±δ Perovskites for Solid Oxide Cell Applications

Sr-secondary phase formation is a potentially significant degradation mode threatening solid-oxide cell (SOC) commercial viability. A first-principles thermodynamic study was performed for rhombohedral perovskite (La1-xSrx) MnO3±δ (LSM) to assess its stability against Sr secondary phase formation in SOC applications. In this work, the Sr secondary phase formation reaction free energies were determined by combining ab initio lattice dynamics calculations for the solid phases and an ab initio thermodynamics approach for the gas phases. Furthermore, this approach goes beyond previous thermodynamic modeling studies by integrating first-principles based point-defect equilibria into the analyses. The modeling results indicate an increased tendency to form SrO oxide from LSM upon decreasing the oxygen partial pressure. Additionally, enhancing factors to form the Sr-related secondary phase from the associated SrO activity in LSM are further quantified by considering the equilibrium of SrO reacting with contaminant gas species as a function of temperature and gas pressure.

Defect and phase stability↗
DOE OSTI2022

Generalization of thermodynamic Langmuir isotherm for mixed‐gas adsorption equilibria

Abstract This work presents a comprehensive thermodynamic model for both pure component isotherms and mixed‐gas adsorption equilibria. A generalization of thermodynamic Langmuir isotherm, the proposed model assumes competitive adsorption of multiple adsorbates on adsorbent surface for mixed‐gas adsorption equilibria, and it applies an area‐based adsorption nonrandom two‐liquid activity coefficient model in the activity coefficient calculations for the adsorbate phase. The resulting generalized Langmuir isotherm properly captures both surface loading dependence and adsorbate phase composition dependence for mixed‐gas adsorption equilibria. The model is validated with accurate representations of gas adsorption equilibrium data for varieties of unary, binary, and ternary gas systems. The model results are further compared with those calculated from extended Langmuir isotherm and Ideal Adsorbed Solution Theory.

Hamid, Usman↗
DOE OSTI2024

Lambda-PFLOTRAN 1.0: a workflow for incorporating organic matter chemistry informed by ultra high resolution mass spectrometry into biogeochemical modeling

Abstract. Organic matter (OM) composition plays a central role in microbial respiration of dissolved organic matter and subsequent biogeochemical reactions. Here, a direct connection of organic matter chemistry and thermodynamics to reactive transport simulators has been achieved through the newly developed Lambda-PFLOTRAN workflow tool that succinctly incorporates carbon chemistry data generated from Fourier transform ion cyclotron resonance mass spectrometry (FTICR-MS) into reaction networks to simulate organic matter degradation and the resulting biogeochemistry. Lambda-PFLOTRAN is a Python-based workflow, executed through a Jupyter notebook interface, that digests raw FTICR-MS data, develops a representative reaction network based on substrate-explicit thermodynamic modeling (also termed lambda modeling due to its key thermodynamic parameter λ used therein), and completes a biogeochemical simulation with the open source, reactive flow and transport code PFLOTRAN. The workflow consists of the following five steps: configuration, thermodynamic (lambda) analysis, sensitivity analysis, parameter estimation, and simulation output and visualization. Two test cases are provided to demonstrate the functionality of the Lambda-PFLOTRAN workflow. The first test case uses laboratory incubation data of temporal oxygen depletion to fit lambda parameters (i.e., maximum utilization rate and microbial carrying capacity). A slightly more complex second test case fits multiple lambda formulation and soil organic matter release parameters to temporal greenhouse gas generation measured during a soil incubation. Overall, the Lambda-PFLOTRAN workflow facilitates upscaling by using molecular-scale characterization to inform biogeochemical processes occurring at larger scales.

58 GEOSCIENCES↗
DOE OSTI2020

Nonrandom two‐liquid activity coefficient model with association theory

Abstract Considered one of the most versatile and widely used classical thermodynamic models to correlate phase equilibria behavior of nonideal systems, the nonrandom two‐liquid (NRTL) theory does not explicitly account for specific chemical associations such as hydrogen bonding. This deficiency has been recognized as the cause for unsatisfactory representation of association systems like methanol–alkanes binaries. This work presents a practical approach to integrate Wertheim's perturbation theory for association contribution with the classical NRTL model. Specifically, the association contribution is calculated from pre‐determined molecule‐specific association strengths while the physical interaction contribution is captured with NRTL binary interaction parameters. The resulting association NRTL model correlates fluid phase equilibria of association systems with few adjustable parameters and offers improved predictive capability for higher order systems.

Hao, Yifan↗
DOE OSTI2025

An eigenvalue-based method for computing the relaxed pressure in compressible multiphase flow with N phases

The modeling of compressible multiphase flows is a decades-old area of study with many applications across various fields. Many of these application areas use stiff pressure relaxation. This process involves the solution of a nonlinear system with N + 1 equations and N + 1 unknowns, where N is the number of phases. The resolution of this system with general equations of state (EOSs) is difficult. Furthermore, nonlinear systems can admit multiple solutions, and current solution methods do not address this possibility. Very recently, a thermodynamic relaxation method was introduced, which effectively maps a relatively simple predictor equation of state onto a more complex target equation of state. In this context, the target EOSs are the chosen EOSs for the thermodynamic model. Furthermore, this thermodynamic relaxation has the benefit of simplifying the stiff pressure relaxation system of equations. In this article, we show this system reduces to a polynomial of degree N, which can be recast as an eigenvalue problem through the use of the associated companion matrix. We show that although this eigenvalue method is generally less efficient than Newton–Raphson iteration, it does not suffer from convergence issues and finds all N roots of the polynomial. Hence, the method provides a fail-safe for root-finding iterative methods and a way to address the issue of multiple solutions to the nonlinear system of equations in stiff pressure relaxation.

Eigenvalue algorithm↗
DOE OSTI2021

Low global warming potential (GWP) refrigerant supermarket refrigeration system modeling and its application

As an environmentally friendly low global warming potential (GWP) refrigerant, Carbon dioxide (CO 2 ) has continuously gained popularity and research attention as alternative refrigerant for supermarket refrigeration system. In this paper, to fulfill the increasing need of accurate Low-GWP supermarket refrigeration models for development of supervisory level control and optimization strategies, a high fidelity model is developed for CO 2 transcritical supermarket refrigeration system which includes compressor rack of low temperature (LT) compressors and medium temperature (MT) compressors, air-cooled gas cooler, evaporator, expansion valves and other auxiliary equipment. A resistance-capacity model structure is proposed to simulate the display cases. Semi-thermodynamic models are proposed to estimate reciprocating compressors volumetric efficiency and power consumption. Furthermore the zone modeling approach is used for evaporator simulation, and air-cooled gas cooler is modeled with distributed modeling method. The expansion valve simulation is based on orifice flow model. To calibrate these models, both manufacture performance data and experimental data are used. The experiments are conducted with a full instrumental CO 2 supermarket refrigeration system installed in Oak Ridge National Lab Flexible Research Platform (FRP). The simulation model can predict the system performance, including power consumption, cooling capacity, mass flow rate, temperature, and pressure, with high accuracy (within ±4%) compared to experimental data. In addition, this developed model has been used to create the system optimum high side pressure for high side expansion valve control, and to generate wide operating range simulation data for developing the batter-equivalent commercial refrigeration system model which can be used in grid interactive control development.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2021

Developing Practical Models of Complex Salts for Molten Salt Reactors

Molten salt reactors (MSRs) utilize salts as coolant or as the fuel and coolant together with fissile isotopes dissolved in the salt. It is necessary to therefore understand the behavior of the salts to effectively design, operate, and regulate such reactors, and thus there is a need for thermodynamic models for the salt systems. Molten salts, however, are difficult to represent as they exhibit short-range order that is dependent on both composition and temperature. A widely useful approach is the modified quasichemical model in the quadruplet approximation that provides for consideration of first- and second-nearest-neighbor coordination and interactions. Its use in the CALPHAD approach to system modeling requires fitting parameters using standard thermodynamic data such as phase equilibria, heat capacity, and others. A shortcoming of the model is its inability to directly vary coordination numbers with composition or temperature. Another issue is the difficulty in fitting model parameters using regression methods without already having very good initial values. The proposed paper will discuss these issues and note some practical methods for the effective generation of useful models.

Besmann, Theodore M. (ORCID:0000000155980550)↗
DOE OSTI2021

A method to apply Zavitsas' aqueous electrolyte model to multicomponent solutions and its equivalence to Zdanovskii's rule

Abstract Zavitsas developed an aqueous electrolyte thermodynamic model and applied it to single electrolyte solutions previously. Here, a method to apply the model to multicomponent solutions is derived by setting the water activity to the mole fraction of free water in the mixture. The method is shown to work effectively for the CsBr‐KNO 3 ‐H 2 O and CsCl‐KNO 3 ‐H 2 O systems at 100.3°C even at concentrations greater than 6 molal. The largest error in the prediction of water activities was just −0.00365 for the CsBr‐KNO 3 ‐H 2 O system and −0.00631 for the CsCl‐KNO 3 ‐H 2 O system. Zavitsas' model naturally obeys Zdanovskii's rule because when two electrolyte solutions with the same fraction of free water are mixed together, water is obviously going to have that same fraction in the mixture. This means that Zavitsas' model will likely work well for the many electrolytes that obey Zdanovskii's rule.

Reynolds, Jacob G.↗
DOE OSTI2019

Dalton’s and Amagat’s laws fail in gas mixtures with shock propagation

A shock propagating through a gas mixture leads to pressure, temperature, and density increases across the shock front. Rankine-Hugoniot relations correlating pre- and post-shock quantities describe a calorically perfect gas but deliver a good approximation for real gases, provided the pre-shock conditions are well characterized with a thermodynamic mixing model. Two classic thermodynamic models of gas mixtures are Dalton’s law of partial pressures and Amagat’s law of partial volumes. We measure post-shock temperature and pressure in experiments with nonreacting binary mixtures of sulfur hexafluoride and helium (two dramatically disparate gases) and show that neither model can accurately predict the observed values, on time scales much longer than that of the shock front passage, due to the models’ implicit assumptions about mixture behavior on the molecular level. However, kinetic molecular theory can help account for the discrepancy. Our results provide starting points for future theoretical work, experiments, and code validation.

42 ENGINEERING↗
DOE OSTI2025

Experiments to validate Thermodynamics and Transport models of Strongly Coupled Dusty Plasma Matter (Final Technical Report for DE-SC0023416)

The goal of this two-year grant is to provide access to the PI to dusty plasma experimental facilities at the DOE-funded Collaborative Research Facility Magnetized Plasma Research Laboratory, Auburn University to become a user of that facility, to obtain experimental data to support another ongoing grant DE-SC0021146 (an Early Career Award to the PI that is focused on modeling of dusty plasma thermodynamics and transport processes), to generate experimental data for funding proposals, and to provide exposure to University of Memphis students to advanced experimental techniques. The following technical accomplishments were made: 1. Development of a novel Bidirectional Electrode Control Arms Assembly (BECAA) for producing perfect 2D grain layers for complex plasma experimentation. BECAA uses movable electrode arms to tilt or move the electrode in a RF discharge from outside the chamber, allowing for the manipulation of grain clouds without needing to change the plasma parameters or gas pressure. This work addresses a longstanding gap in the literature for a method to produce clusters of selectable number of grains and that are perfectly two dimensional as opposed to being only quasi-2D. 2. Experimental investigation of the structural properties of finite-N clusters with N=2 to 50. Individual particle behavior in clusters could vary from grain to grain and this study measured systematically produced clusters for two different grain sizes. Analysis (funded by another grant DE-SC0021146) is currently underway to quantify the differences between grains that are found on the surface vs. the interior of clusters, the shell structure, and the decay of correlations in position, velocity, and kinetic energy. 3. An experimental method to measure the structural entropy of clusters was developed by observing the self-induced structural transitions between various possible arrangements. In a series of heating and cooling cycles, the number of times each possible arrangement was attained was experimentally observed and used to compute the probability of existence of that arrangement, and subsequently the configurational entropy of the cluster. Analysis (funded by another grant DE-SC0021146) is currently underway to produce the entropy of clusters as a function of the number of grains and use the same to compute thermodynamic state variables for 2D complex plasma/grain clusters. 4. A preliminary experimental study of multibody collisions between N grains (N=2 – 10) was conducted. The clusters were produced using the BECAA technique and velocities were imparted to the grains using manipulation laser pulses. Analysis (funded by another grant DE-SC0021146) is currently underway to develop a theoretical framework to describe multibody collisions analogous to classical two-body interactions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Topical Report on Water Balance Model Development and Initial Application

This topical report summarizes results from developing a generic dynamic water balance tool which can be used by coal-fired power plant operators to evaluate FGD water management strategies. EPRI is using a statistical software model created on the GoldSim platform and data from coal-fired energy plants to determine the fraction of power plant water needs that can be satisfied by electrodialysis reversal treatment and, correspondingly, the discharge volume and/or byproducts that still need to be managed. The generic model is designed to accept a number of easily obtainable process variables and generate probabilities of resultant water characteristics. The model is also constructed as a universal process flow; one that can be adjusted to accommodate most major components of an FGD blowdown handling system. The generic model is dynamic and utilizes user-inputted operational variation in water quality and quantity to best provide a probabilistic representation of potential performance outcomes. In order to account for complexities of water composition and treatment, a separate chemical modeling software is used to evaluate water chemistry and unit operation performance. The thermodynamic modeling software, OLI Studio, is employed to convert ion data from laboratory analyses into balanced molecular concentrations that are classified as dissolved and suspended solids. The dynamic water model is designed to import a total of 11 species commonly found in FGD blowdown water as molecular concentrations.Water management unit operations are used in the model and categorized by either physical unit operation blocks or chemical unit operations blocks. As the name implies, the physical water treatment blocks are unit operations that are exclusively physical changes to the system (e.g., separations and volume changes). The chemical water treatment blocks are unit operations where chemical changes occur and require the use of OLI Studio to characterize subsequent streams. Water and material balances from three unique coal-fired power plants were obtained and preliminary modeling on the dynamic water balance tool was performed using the extracted data from these power plants. These three power plants were chosen due to their unique configurations and water profiles that represented a variety of operational changes. Results from using the dynamic tool with the three plant configurations are presented and discussed in the context of the existing tool and planned future capabilities.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2021

Modeling dissociation of ionic liquids with electrolyte NRTL model

Dissociation of ionic liquids (ILs) is critical to electrochemical process performance by affecting ionic conductivity. Prior studies of thermodynamic modeling either assumed no dissociation or complete dissociation that were not consistent with experimental observations. This work presented an eNRTL model for three aqueous imidazolium-based IL systems: [Emim][EtSO 4 ], [Emim][TFA], and [Emim][TfO]. The model accurately correlates the IL dissociation and VLE data in the entire concentration range from pure IL to infinite dilution aqueous solution. The system non-ideality is dominated by the dissociation chemistry. Dissociation behavior is driven by strong short-range interaction of the molecular IL-water pair relative to the dissociated IL-water pair, consistent with prior studies indicating that molecular IL and water can form an energetically favorable complex. The interaction energy parameters calculated in this work qualitatively agrees with the hydrogen bonding ability of the anions previously reported. This research not only provides a model that well describes the IL dissociation but also demonstrates the physical significance of eNRTL model.

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

Optimization of Desalination Systems with Detailed Water Chemistry through Integration of Reaktoro in WaterTAP

Chemistry predictions are critical for an accurate estimation of performance and costs in desalination process models, which allows for the estimation of the value of new technologies and the viability of treating new water sources. Herein, we present how an implicit function formulation can be used to integrate the chemical modeling package, Reaktoro, into the techno-economic assessment and modeling platform, WaterTAP. This approach resolves the critical issues of integrating large-scale thermodynamic models and databases into equation-oriented process models while allowing more flexibility relative to previously presented surrogate-based methods. We describe how this integration into Pyomo and WaterTAP models is implemented and used through the open-source package Reaktoro-PSE . We first validate this integration approach by performing optimization on a previously presented desalination treatment train with softening and acid addition as the pretreatment steps. Then, to demonstrate the value of this approach, we extend the cost-optimization problem to include the simultaneous addition of lime and soda ash for softening, and HCl and H 2 SO 4 in the acidification steps. Finally, we were able to confirm the previously established results that were obtained by using surrogate models and demonstrate that the implicit function approach enables exploration of different feedwater compositions and a larger number of chemicals and their combinations.

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

Uncertainty propagation in pore water chemical composition calculation using surrogate models

Performance assessment in deep geological nuclear waste repository systems necessitates an extended knowledge of the pore water chemical conditions prevailing in host-rock formations. In the last two decades, important progress has been made in the experimental characterization and thermodynamic modeling of pore water speciation, but the influence of experimental artifacts and uncertainties of thermodynamic input parameters are seldom evaluated. In this respect, we conducted an uncertainty propagation study in a reference geochemical model describing the pore water chemistry of the Callovian-Oxfordian clay formation. Nineteen model input parameters were perturbed, including those associated to experimental characterization (leached anions, exchanged cations, cation exchange selectivity coefficients) and those associated to generic thermodynamic databases (solubilities). A set of 13 quantities of interest were studied by the use of polynomial chaos expansions built non-intrusively with a least-squares forward stepwise regression approach. Training and validation sets of simulations were carried out using the geochemical speciation code PHREEQC. The statistical results explored the marginal distribution of each quantity of interest, their bivariate correlations as well as their global sensitivity indices. The influence of the assumed distributions for input parameters uncertainties was evaluated by considering two parametric domain sizes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗