Research Article |
Corresponding author: Christomir Christov ( ch.christov@shu.bg ) Academic editor: Michaela Beltcheva
© 2022 Stanislav Donchev, Tsvetan Tsenov, Christomir Christov.
This is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY 4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Citation:
Donchev S, Tsenov T, Christov C (2022) Development of accurate chemical thermodynamic database for geochemical storage of nuclear waste. Part II: Models for predicting solution properties and solid-liquid equilibrium in binary nitrate systems. In: Chankova S, Peneva V, Metcheva R, Beltcheva M, Vassilev K, Radeva G, Danova K (Eds) Current trends of ecology. BioRisk 17: 389-406. https://doi.org/10.3897/biorisk.17.77487
|
The main purpose of this study is to develop new thermodynamic models for solution behavior and solid-liquid equilibrium in 10 nitrate binary systems of the type 2–1 (Mg(NO3)2-H2O, Ca(NO3)2-H2O, Ba(NO3)2-H2O, Sr(NO3)2-H2O, and UO2(NO3)2-H2O), 3–1 (Cr(NO3)3-H2O, Al(NO3)3-H2O, La(NO3)3-H2O, Lu(NO3)3-H2O), and 4–1 (Th(NO3)4-H2O) from low to very high concentration at 25 °C. To construct models, we used different versions of standard molality-based Pitzer approach. To parameterize models, we used all available raw experimental osmotic coefficients data (φ) for whole concentration range of solutions, and up to supersaturation zone. The predictions of developed models are in excellent agreement with φ-data, and with recommendations on activity coefficients (γ±) in binary solutions from low to very high concentration. The Deliquescence Relative Humidity (DRH), and thermodynamic solubility product (as ln K°sp) of 12 nitrate solid phases, precipitating from saturated binary solutions have been calculated. The concentration-independent models for nitrate systems described in this study are of high importance for development of strategies and programs for nuclear waste geochemical storage.
Nuclear waste sequestration, Chemical modelling, Pitzer approach, DRH and K°sp of Mg(NO3)2.6H2O(s), Ca(NO3)2.4H2O(s), Ca(NO3)2.3H2O(s), Ba(NO3)2(s), Sr(NO3)2(s), UO2(NO3)2.6H2O(s), Al(NO3)3.9H2O(s), Cr(NO3)3(s), La(NO3)3.6H2O (s), La(NO3)3(s), Lu(NO3)3.5H2O(s) and Th(NO3)4.6H2O(s)
Computer models that predict solution behavior and solid-liquid-gas equilibria close to experimental accuracy have wide applicability. They can simulate the complex changes that occur in nature and can replicate conditions that are difficult or expensive to duplicate in the laboratory. Such models can be powerful predictive and interpretive tools to study the geochemistry of natural waters and mineral deposits, solve environmental problems and optimize industrial processes. However, development of comprehensive models for natural systems, with their complexity and sensitivity, is a very difficult, time consuming and challenging task. The specific interaction approach for describing electrolyte solutions to high concentration introduced by
A long term safety assessment of a repository for radioactive waste requires evidence, that all relevant processes are known and understood, which might have a significant positive or negative impact on its safety. It has to be demonstrated, that the initiated chemical reactions don’t lead to an un-due release of radionuclides into the environmental geo-, hydro-, and bio-sphere. One key parameter to assess the propagation of a radionuclide is its solubility in solutions interacting with the waste. Solubility estimations can either be based on experimental data determined at conditions close to those in the repository or on thermodynamic calculations. A so called “thermodynamic database” created from experimental data is the basis for thermodynamic model calculations. Since the disposal of radioactive waste is a task encompassing decades, the database is projected to operate on a long-term basis. Chemical models that predict equilibrium involving mineral, gas and aqueous phases over a broad range of solution compositions and temperatures are useful for studying the interactions between used nuclear fuel waste and its surroundings. The reliability of such predictions depends largely on the thermodynamic database. Waters of high salinity are not a typical of many geochemical environments which may be chosen as future nuclear waste repository sites. This suggests that an accurate description of highly saline waters should be required for modeling of chemical interactions in and around nuclear repositories. Currently, the most accurate description of saline waters uses the Pitzer ion interaction model. Extensive thermodynamic databases, which are based on the Pitzer ion interaction model was developed within the Yucca Mountain Project (YMTDB: data0.ypf.r2) (
Nitrates are expected to play a significant role in the context of the underground geochemical repository of nuclear waste (
In our previous study (
The models for nitrate binary systems have been developed on the basis of Pitzer’s semi-empirical equations (
According to Pitzer theory electrolytes are completely dissociated and in the solution there are only ions interacting with one to another (
(1)
Equation (1) is symmetric for anions. The subscripts c and a in eqn 1 refer to cations and anions, and m is their molality; z is the charge of the M+ ion. B and Φ represent measurable combinations of the second virial coefficients; C and ψ represent measurable combinations of third virial coefficients. B and C are parameterized from single electrolyte data, and Φ and ψ are parameterized from mixed solution data. The function F is the sum of the Debye-Hueckel term,
-Aφ [√ I / (1 + b√ I) + (2/b) (ln(1 + b√ I)] , (2)
and terms with the derivatives of the second virial coefficients with respect to ionic strength (see
For the interaction of any cation M and any anion X in a binary system MX-H2O, Pitzer assumes that in Eq. (1) B has the ionic strength dependent form:
BMX = β(°)MX + β(1)MX g(α1√I) + (3)
+β(2)MX g(α2√I), (3A)
where g(x) = 2[1 - (1 + x)e-x] / x2 with × = α1√ I or α2√ I. α terms are function of electrolyte type and does not vary with concentration or temperature.
In Eq. 1, the Φ terms account for interactions between two ions i and j of like charges. In the expression for Φ,
Φij = θij + Eθij (I), (4)
θij is the only adjustable parameter. The Eθij (I) term accounts for electrostatic unsymmetric mixing effects that depend only on the charges of ions i and j and the total ionic strength. The ψijk parameters are used for each triple ion interaction where the ions are not all of the same sign. Their inclusion is generally important for describing solubilities in concentrated multicomponent systems. Therefore, according to the basic Pitzer equations, at constant temperature and pressure, the solution model parameters to be evaluated are: 1) pure electrolyte β(0), β(1), and Cφ for each cation-anion pair; 2) mixing θ for each unlike cation-cation or anion-anion pair; 3) mixing ψ for each triple ion interaction where the ions are all not of the same sign.
Fluids commonly encountered in natural systems include dissolved neutral species (such as carbon dioxide (CO2(aq), SiO2(aq), and Al(OH)3°(aq)). To account neutral specie interactions in aqueous solutions the UCSD Chemical Modelling Group included in their models additional terms to Pitzer equations, denoted as λN,X or λN,A, and ζ N, A,X. (Eq. (1)) (
Some authors found that there are some restrictions limited the potential of the model to describe correctly activity and solubility properties in some binary electrolyte systems with minimum one univalent ion (see
To describe the high concentration solution behaviour of systems showing a “smooth” maximum on γ± vs. m dependence, and to account strong association reactions at high molality,
In this study we developed new thermodynamic models for solution behavior and solid-liquid equilibrium in 10 nitrate binary systems of the type 2–1 (Mg(NO3)2-H2O, Ca(NO3)2-H2O, Ba(NO3)2-H2O, Sr(NO3)2-H2O, and UO2(NO3)2-H2O), 3–1 (Cr(NO3)3-H2O, Al(NO3)3-H2O, La(NO3)3-H2O, Lu(NO3)3-H2O), and 4–1 (Th(NO3)4-H2O) from low to very high concentration at 298.15 K. New sets of Pitzer ion interaction binary parameters are evaluated using available raw experimental osmotic coefficients (φ) data for whole molality range of solutions. Rard and co-authors (1977, 1981) reported an extensive experimental activity database for rare earth nitrate systems. Data of
In parameterization we used the value of Debye-Hückel term (Aφ) equals to 0.39147 (
On next Figure
Comparison of model calculated (lines) osmotic coefficients (φ) of Mg(NO3)2 Ca(NO3)2, Ba(NO3)2, Sr(NO3)2, UO2(NO3)2, Cr(NO3)3, Al(NO3)3, La(NO3)3, Lu(NO3)3, and Th(NO3)4 in binary solutions 2–1 (Mg(NO3)2-H2O, Ca(NO3)2-H2O, Ba(NO3)2-H2O, Sr(NO3)2-H2O, and UO2(NO3)2-H2O), 3–1 (Cr(NO3)3-H2O, Al(NO3)3-H2O, La(NO3)3-H2O, Lu(NO3)3-H2O), and 4–1 (Th(NO3)4-H2O) against molality at T = 298.15 K, with recommendations in literature (symbols). For Mg(NO3)2-H2O and Ca(NO3)2-H2O systems an enlargement of the low molality corner is also given. Heavy solid lines represent the predictions of the developed in this study and accepted models. Dashed-dotted, dashed and light solid lines represent the predictions of the reference models of
The models for all nitrate binary systems under study are also validated by comparison with recommendations given in literature (
Deliquescence of single inorganic salt or their mixture is a process of spontaneous solid-liquid phase change. It is a process in which a soluble solid substance sorbs water vapor from the air to form a thermodynamically stable saturated aqueous solution on the surface of the particle. It is occurring when relative humidity (RH) in the gas-phase environment is at, or above deliquescence relative humidity (DRH) of the salt, or mutual deliquescence relative humidity (MDRH) of a salt mixture. Within the solid-liquid equilibrium model, relative humidity is related to water activity(aw) (
aw = Pw/ P°w = RH/100, (5)
where Pw and P°w are the vapor pressure of the saturation solution and pure water, respectively, at given temperature. As a result, both DRH and MDRH of saturated surface solutions depend of temperature, the salt stoichiometry, and the solution composition. This process is of interest in many areas, such as heterogeneous chemistry of inorganic salts, corrosion of metals in wet atmosphere, in studies of chemistry of sea-type aerosol atmospheric system (
Comparison between model calculated and recommended values of the Deliquescence Relative Humidity [DRH (%) = aw (sat) . 100; where aw(sat) is activity of water at saturation] and of the logarithm of the thermodynamic solubility product (as lnK°sp) of nitrate solid phases crystallizing from saturated binary solutions at T = 25° C.
Salt composition | m (sat) (exp) (mol.kg-1) | lnK°sp | DRH(%) | ||
---|---|---|---|---|---|
This work calculated | Reference data | This work calculated | Reference data a | ||
Mg(NO3)2.6H2O(cr) | 5.06 a | 7.0098 | 7.02b | 52.32 | 52.90 |
Ca(NO3) 2 .4H2O(cr) (stable solid) | 8.41 a | 4.4362 | 4.53b | 49.07 | 49.10 |
Ca(NO3) 2 .3H2O(cr) (metastable solid) | 14.77 a | 6.6449 | 5.34b (m(sat) = 15.0 m) | 22.52 | - |
Ba(NO3)2(cr) | 0.39 a | -5.125 | - | 98.61 | 98.60 |
Sr(NO3)2 .4H2O (cr) | 3.76 a | 0.0327 | - | 84.83 | 84.80 |
UO2(NO3)2 .6H2O(cr) | 3.21 a | 5.3022 | 5.251c | 73.44 | 73.60 |
Al(NO3)3 .9H2O (cr) | 3.16 a | 4.3081 | - | 59.88 | 60.20 |
Cr(NO3)3 (cr) e | 1.4 e | 1.2097 | - | 86.38 | - |
La(NO3)3 .6H2O (cr) d | 4.615 d | 2.1599 | 2.97d | 62.38 | - |
La(NO3)3 (cr) a | 2.94 a | 1.4704 | - | 77.73 | 77.60 |
Lu(NO3)3 5H2O (cr) | 6.815 d | 10.7681 | 10.67 | 31.27 | - |
Th(NO3)4 .6H2O(cr) | 4.00 a | 4.4886 | 4.71c (as Th(NO3)4.5H2O) | 54.46 | 55.0 |
In this study we determine the thermodynamic solubility products (as K°sp) of solid phases, precipitating from saturated nitrate binary solutions, s.a. anhydrous Ba(NO3)2(s) and hydrate Ca(NO3)2.3H2O(s), precipitating in Ba(NO3)2-H2O and Ca(NO3)2-H2O. The K°sp have been determined on the basis of evaluated binary parameters and using experimental m(sat) solubility data, and using the following relationships (
K°sp (Ba(NO3)2) = 4 . γ(±)(sat)3 . m(sat)3
K°sp(Ca(NO3)2.3H2O) = 4 . γ(±)(sat) 3 . m(sat) 3 . aw(sat) 3 (6)
As a next step, using the accepted new developed parameterizations, and experimentally determined molalities (m(sat) of the saturated binary solutions (
In this study we developed new thermodynamic models for solution behavior and solid-liquid equilibrium in 10 nitrate binary systems of the type 2–1 (Mg(NO3)2-H2O, Ca(NO3)2-H2O, Ba(NO3)2-H2O, Sr(NO3)2-H2O, and UO2(NO3)2-H2O), 3–1 (Cr(NO3)3-H2O, Al(NO3)3-H2O, La(NO3)3-H2O, Lu(NO3)3-H2O), and 4–1 (Th(NO3)4-H2O) from low to very high concentration at 25 °C. To parameterize models for binary systems we used all available raw experimental osmotic coefficients data (φ) for whole concentration range of solutions, and up to saturation point. Data for supersaturation zone, available for Ca(NO3)2-H2O, UO2(NO3)2-H2O, and Th(NO3)4-H2O systems, are also included in parameterization. To construct models, we used different versions of standard molality-based Pitzer approach. It was established that with only 2 exceptions (Ba(NO3)2-H2O, and UO2(NO3)2-H2O) application of extended approach with 4 parameters (β(0), β(1), β(2),and Cφ) and variation of α2 term in fundamental Pitzer equations leads to the lowest values of standard model-experiment deviation. The predictions of new developed here models are in excellent agreement with experimental osmotic coefficients data (see Fig.
We wish to thank the reviewers (Dr. Krasimir Kostov and anonymous reviewer) for their constructive suggestions and helpful comments. The manuscript was improved considerably through their comments. The work was supported by the European Regional Development Fund, Project BG05M2OP001-1.001-0004, and by Shumen University Research Program, Project RD-08-131/04.02.2021