INTRODUCTION
The prediction of the solubility in aqueous electrolyte solutions is important for a variety of applications in the chemical and geochemical processes, seawater systems, and evaporation as well as desalination. Salt solubility data are important as a tool for the design and simulation of unit operations such as drowning-out crystallization or liquidliquid extraction. The investigation of the thermodynamics and phase diagram of the system is of theoretical and practical importance.1,2 In the salt lakes of western China, these alkali metal salts coexist with other minerals containing boron, calcium, magnesium, and chloride.3 To scientifically exploit these natural resources preliminary investigation of sodium, potassium and calcium salt solution chemistry is necessary. Thermodynamic properties of the ternary systems (NaCl + CaCl2 + H2O and KCl + CaCl2 + H2O) are of essential importance in the extraction of sodium or potassium from natural salt brine mainly containing sodium, potassium, lithium, magnesium and calcium.
Pitzer’s ion-interaction model4 and its extended Harvie and Weare model5-7 are very reliable for predicting the mineral solubility in natural water system with high ionic strengths over the wide temperature range from 0 to 300 ℃.8-11 The solubilities of calcium chloride and magnesium chloride were determined by prutton,12 and calculating the activity coefficient and osmotic coefficients of electrolytes in seawater and synthetic salt lake brines were calculated by Song, P. S.13,14 Li Ya-hong15-17 reported. The solubility in the ternary HCl-LiCl-H2O, HCl-MgCl2-H2O and LiCl-MgCl2-H2O systems at 273 and 293 K, using the ion-interaction Pitzer model and predicting the solubility isotherm of the NaCl-RbCl-H2O, KCl-CsCl-H2O and KBr-CsBr-H2O systems at 298 K. Some other systems were also reported, such as NaCH3COO + NaCl + H2O,18 NaCl and KNO3,19 LiCl-MgCl2-H2O system.20 However, the application of the Pitzer model for predicting the component solubility of salt lake brine systems of NaCl + CaCl2 + H2O and KCl+CaCl2+H2O at 50 ℃, has never been reported. Therefore, the solubility evaluation for the NaCl + CaCl2 + H2O and KCl + CaCl2 + H2O systems at 50 ℃ was performed in our lab.
In this paper, the solubility of the ternary systems NaCl + CaCl2 + H2O and KCl + CaCl2 + H2O was elaborately measured at 50 ℃ and an empirically Pitzer ion-interaction model was established based on those solubility data. The physicochemical properties (refractive index) of the equilibrium solutions were determined and a study on the prediction of the solubility was also performed.
EXPERIMENTAL SECTION
Apparatus and Reagents
A thermostatic shaker (model HS-4) whose temperature could be controlled to 0.02 K was used for the measurement of phase equilibrium. The chemicals used were of analytical grade and obtained from either the Tianjin Chemical Reagent Manufactory or the Shanghai Chemical Plant: calcium chloride (CaCl2·6H2O, 99.5 mass %), sodium chloride NaCl, g 99.8 mass %) and potassium chloride (KCl, g 99.5 mass %). Doubly deionized water was used to prepare the series of saturated solutions.
Experimental method
Various mixtures of salt and water were made by starting with a ground 250 cm3 Erlenmeyer flask containing only one salt and water and in each subsequent run more of the second salt was added to the solution and solid left from the previous run. The flask was immersed in the thermostat and the solution and solid in the flask were stirred with a magnetic stirrer. Each sample was stirred at a specific constant temperature for 72 h, and then kept static for about 6 h. A sample of the saturated solution was then taken with a pipette. The sample was transferred to a weighed 30 cm3 ground quartz beaker with cover. The salt concentration in pure saturated salt aqueous solution and the total salt concentration in the three component solutions were determined by evaporation to dryness, fusing, and weighing. The wet residuals were analyzed in the same way as for the solution. The composition of the solid phase in the wet residues was identified by the method of Schreinemaker.
Analytical method
The Cl‒ ion concentration in the liquids and their corresponding wet residues of the solid phases were analyzed by titration with a standard solution of AgNO3 in the presence of three drops of 0.1% (w/v) KCrO4 as an indicator (precision within ± 0. 2 mass %. The Ca2+ ion concentration was determined by titration with EDTA standard solution in the presence of indicator of Eriochrome Black-T21 (uncertainty of ± 0.2%). An Abbe refractometer (model WZS-1) was used to measure the refractive index (nD) with an accuracy of ± 0.0001.
RUSULT AND DISCUSSION
The experimental data on the solubilities and refractive index of the ternary systems NaCl + CaCl2 + H2O and KCl + CaCl2 + H2O at 50 ℃ are presented in Tables 1 and 2, respectively. The electrolyte concentration values of the liquid phase in the equilibrium solution are expressed in mass fraction. According to the experimental data in Tables 1 and 2, the equilibrium phase diagrams of the systems NaCl + CaCl2 + H2O and KCl + CaCl2 + H2O are plotted, as shown in Figs. 1 and 2 at 50 ℃, respectively. The solid phases in equilibrium with saturated solution are CaCl2·2H2O and MCl in the MCl (M = Na and K) + CaCl2 + H2O system at 50 ℃.
In Figs. 1 and 2 with solid lines, points A and B are the solubilities of the single-salts of sodium chloride (or potassium chloride) and calcium chloride dihydrate. Point C is a eutectic point of sodium chloride(or potassium chloride) and calcium chloride dihydrate (NaCl + CaCl2·2H2O or KCl + CaCl2·2H2O). There are two saturated curves corresponding to curves AC and BC, indicating the saturation of single salts. The phase diagram consists of two crystallization regions corresponding to the large area of NaCl (or KCl) and the relative small area of CaCl2·2H2O. Obviously, the system belongs to the simple eutectic type, and neither double salts nor solid solutions are found.
On the basis of experimental data in Tables 1 and 2, relationship of the solution physicochemical property with the concentration of calcium chloride is shown in Figs. 3 and 4. It can be found that the refractive indexes of the aqueous solutions in the ternary system, changed gradually and regularly with the content change of calcium chloride.
Table 1.Solubility data of the NaCl + CaCl2+ H2O system at 50 ℃
Table 2.Solubility data of the KCl+ CaCl2+ H2O system at 50 ℃
In Figs. 1 and 2 with solid lines, points A and B are the solubilities of the single-salts of sodium chloride (or potassium chloride) and calcium chloride dihydrate. Point C is a eutectic point of sodium chloride (or potassium chloride) and calcium chloride dihydrate (NaCl + CaCl2·2H2O or KCl + CaCl2·2H2O). There are two saturated curves corresponding to curves AC and BC, indicating the saturation of single salts. The phase diagram consists of two crystallization regions corresponding to the large area of NaCl (or KCl) and the relative small area of CaCl2·2H2O. Obviously, the system belongs to the simple eutectic type, and neither double salts nor solid solutions are found.
On the basis of experimental data in Tables 1 and 2, relationship of the solution physicochemical property with the concentration of calcium chloride is shown in Figs. 3 and 4. It can be found that the refractive indexes of the aqueous solutions in the ternary system, changed gradually and regularly with the content change of calcium chloride.
Fig. 1.The phase diagram of the ternary system NaCl + CaCl2 + H2O at 50 ℃.
Fig. 2.The phase diagram of the ternary system KCl + CaCl2 + H2O at 50 ℃.
Fig. 3.Physicochemical property vs composition in the ternary system (NaCl-CaCl2-H2O) at 50 ℃.
Fig. 4.Physicochemical property vs composition in the ternary system (KCl-CaCl2-H2O) at 50 ℃.
Calculation of the standard solubility product KMX
The solubility of hydrated salt in concentrated electrolyte solutions can be calculated from thermodynamic considerations provided that equilibrium constants are known and activity coefficients can be obtained. For a hydrated salt MνM XνX·ν0H2O, the solubility equilibrium constant Ksp, at a de finite temperature and pressure for the dissolution reaction
is expressed by
where mi and γi represent the concentration expressed as molality and activity coefficient of the ions, respectively. The water activity aW is related to the osmotic coefficient
Here, Mw = 0.018015 kg․mol-1 is the molecular mass of water; ф is the osmotic coefficient.
Parameterization
Pitzer’s binary parameters (β(0), β(1), β(2) and CФ) for pure electrolytes of CaCl2 + NaCl + H2O and CaCl2 + KCl + H2O systems at 50 ℃ were available in the literature. The parameters of single-salt CaCl2 were fitted from osmotic coefficients by least-square method. The parameters of single-salts NaCl and KCl were calculated by using the expressions of the temperature dependency of the ion interaction parameters. These values are listed in Table 3.
For the CaCl2 + NaCl + H2O and CaCl2 + KCl + H2O systems, the expression for the osmotic coefficient is
Table 3.Pitzer’s binary parameters for single electrolytes at 50 ℃
where M is Na and K. and
The constants a1 and a2 are normally 1.4 and 2.0 mol1/2·kg-1/2 respectively, θM, Ca and ΨCa, M,Cl are the mixing parameters.
The activity coefficients of MCl (M = Na or K) and CaCl2 in the ternary mixture are
The function F for the mixtures is given by:
Here
Table 4.Calculated logarithm of the solubility equilibrium constant
Table 5.Pitzer mixing parameters for systems investigated at 50 ℃
These equations contain the additional functions, which were described:24
where b = 1.2 mol.Kg-1. Value of the Debye-Hückel coefficient Aφ is 0.408 kg1/2.mol-1/2 at 50 ℃.25
For θM,Ca and ψCa,M,Cl measurements of the activities of the complex salts in the system containing high concentrations of MCl (M = Na or K) and CaCl2 are not available, hence the evaluation of the mixing parameters θM,Ca and ψCa,M,Cl relied on solubility data and solubility product of hydrated solid at 50 ℃ was obtained by a least-squares optimization procedure on the solubility data of corresponding ternary system. All of the parameters used in our calculations are listed in Tables 4 and 5.
Calculated Solubilities
The solubility data and the relevant physicochemical property data of the MCl (M = Na or K) + CaCl2 + H2O system at 50 ℃ were measured, and the results are shown in Tables 1 and 2, respectively. Using the chemical equilibrium model and the above parameters, the calculated results of the solubility are shown in Figs. 5 and 6 at 50 ℃, respectively. It is shown that the predicted values using the chemical equilibrium model agree well with experimental values. For the MCl (M = Na or K) + CaCl2 + H2O system with high ionic strengths, this agreement indicates that the parameters obtained in this work are reliable and that the chemical equilibrium model of Harvie is capable of predicting equilibiria in the system studied.
Fig. 5.Calculated and experimental solution isotherms of NaCl + CaCl2 + H2O at 50 ℃.
Fig. 6.Calculated and experimental solution isotherms of KCl + CaCl2 + H2O at 50 ℃.
CONCLUSION
The experimental solubility data and the relevant physicochemical property data of the aqueous systems NaCl + CaCl2 + H2O and KCl + CaCl2 + H2O at 50 ℃ were determined. Based on the solubility data measured, the isothermal phase diagram and diagram of physicochemical property vs composition are constructed. The single-salt Pitzer parameters of sodium chloride, potassium chloride and calcium chloride were calculated, and the mixing ion-interaction parameters of the NaCl + CaCl2 + H2O and KCl + CaCl2 + H2O systems could be fitted satisfactorily. The solubility product of NaCl, KCl and CaCl2·2H2O have been calculated. The calculated results agree well with the experimental values. This study could be useful for solubility prediction for more complicated systems and supply a theoretical basis for the exaction of sodium, potassium and calcium from salt lake brine.
References
- Raatikainen, T.; Laaksonen, A. Atmos. Chem. Phys. 2005, 5,2475. https://doi.org/10.5194/acp-5-2475-2005
- Yuanbing, L.; George, P. J. Chem. Eng. Data 2004, 49, 1263. https://doi.org/10.1021/je034238p
- Mining Association of Qinghai, Survey of Mining in Qinghai,people’s Press of Qinghai, 1998; p 85.
- Pitzer, K. S. J. Phys. Chem. 1973, 77, 268. https://doi.org/10.1021/j100621a026
- Harvie, C. E.; Weare, J. H. Geochim. Cosmochim. Acta 1980,44, 981. https://doi.org/10.1016/0016-7037(80)90287-2
- Harvie, C. E.; Eugster, H. P.; Weare, J. H. Geochim. Cosmochim. Acta 1982, 46, 1603. https://doi.org/10.1016/0016-7037(82)90317-9
- Harvie, C. E.; Moller, N.; Weare, J. H. Geochim. Cosmochim. Acta 1984, 48, 723. https://doi.org/10.1016/0016-7037(84)90098-X
- Felmy, A. R.; Weare, J. H. Geochim. Cosmochim. Acta 1986,50, 2771. https://doi.org/10.1016/0016-7037(86)90226-7
- Moller, N. Geochim. Cosmochim. Acta 1988, 52, 821. https://doi.org/10.1016/0016-7037(88)90354-7
- Pabalan, R. T.; Pitzer, K. S. Geochim. Cosmochim. Acta 1987,51, 2429. https://doi.org/10.1016/0016-7037(87)90295-X
- Spencer, R. T.; Moller, N. Geochim. Cosmochim. Acta 1990,54, 575. https://doi.org/10.1016/0016-7037(90)90354-N
- Purutton C. F.; Tower, O. F. J. Am. Chem. Soc. 1932, 54, 3040. https://doi.org/10.1021/ja01347a003
- Song, P. S.; Yao, Y. CALPHAD 2001, 25, 329. https://doi.org/10.1016/S0364-5916(01)00053-0
- Song, P. S.; Yao, Y. CALPHAD 2003, 27, 342.
- Li, Y.-H.; Song, P.-S.; Xia, S.-P.; Gao, S.-Y. CALPHAD 2000,24, 295. https://doi.org/10.1016/S0364-5916(01)00006-2
- Li, Y.-H.; Song, P.-S.; Xia, S.-P.; Li, W.; Gao, S.-Y. Chin. J. Chem. 2005, 23, 953. https://doi.org/10.1002/cjoc.200590953
- Li, Y.-H.; Song, P.-S.; Xia, S.-P.; Li, W.; Gao, S.-Y. CALPHAD2006, 30, 61. https://doi.org/10.1016/j.calphad.2005.10.008
- Joachim, D.; Michael S. J. Chem. Eng. Data 2007, 52, 1784. https://doi.org/10.1021/je7001495
- Park, T. Y. J. Kor. Chem. Soc. 2009, 53(4), 453. https://doi.org/10.5012/jkcs.2009.53.4.453
- Yang, J. M.; Ji, J. Russian J. Phys. Chem. 2010, 84(7), (inPress).
- Qinghai Institute of Salt Lakes of CAS, Analytical Methods of Brines and Salts, 2nd, Science Press: Beijing, 1988; p 36.
- Ananthasway, J.; Atkinson, G. J. Chem. Eng. Data. 1985, 30,120. https://doi.org/10.1021/je00039a035
- Holmes, H. M. ; Mesmer, R. E. J. Phys. Chem. 1983, 87, 1242. https://doi.org/10.1021/j100230a030
- Rard, J. A.; Platford, R. F.; Pitzer, K. S. Activity Coefficients in Electrolyte Solutions, 2nd, CRC Press: Boca Raton, FL, 1991;p 75.
- Balarew, C.; Tepavitcharova, S.; Voigt, W. J. Solution Chem.2001, 30, 815. https://doi.org/10.1023/A:1012288316884
Cited by
- The phase equilibriums in the NH4Cl-CaCl2-H2O system at 50 and 75°C and their Pitzer model representations vol.88, pp.13, 2014, https://doi.org/10.1134/S0036024414130214
- THC simulation of halite scaling in deep geothermal single well production vol.65, 2017, https://doi.org/10.1016/j.geothermics.2016.09.009
- Thermodynamic modeling of aqueous Ca2+– Na+ – K+ – Cl− quaternary system vol.409, 2016, https://doi.org/10.1016/j.fluid.2015.09.048
- Solid-liquid phase equilibria at 50 and 75°C in the NaCl + MgCl2 + H2O system and the pitzer model representations vol.87, pp.13, 2013, https://doi.org/10.1134/S003602441313013X
- The phase diagrams and Pitzer model representations for the system KCl + MgCl2 + H2O at 50 and 75°C vol.86, pp.13, 2012, https://doi.org/10.1134/S0036024412130146