Se desarrolla un procedimiento analítico para predecir el equilibrio sólido-líquido de soluciones acuosas de dos electrólitos 1:1 con unión común, usando la ecuación de Reilly, Wood y Robinson en el cálculo de coeficientes de actividad de los electrólitos mezclados y el modelo de Pitzer para evaluar las propiedades de los sistemas binarios involucrados. Se dan valores de los parámetros de los modelos mencionados para los sistemas NaNO3 - NaCl -&nbsp;H2O. KCl - NaCl - H2O, KCI- KNO3&nbsp;- H2O&nbsp;y NaNO3 - KNO3 - H2O,&nbsp;a 25°C. Los resultados obtenidos son satisfactorios respecto a la información experimental disponible, excepto para el último sistema considerado

Correa H., Horacio

Vega V., Ricardo

Contribuciones Científicas y Tecnológicas

Contribuciones Científicas y Tecnológicas NO 52- Area Ingeniería Química V - Enero de 1982 equilibrio sólido -líquido en soluciones acuosas de electrólitos fuertes.! l. uso de los modelos de reilly- wood- robinson y de pitzer* Horacio Correa H. y Ricardo Vega V.** RtSUJl1é-JV: Se desarrolla un procedimiento analítico para pre-decir el equilibrio sólido-líquido de soluciones acuosas de dos electrólitos 1:1 con w1 ion comzín, usando la ecuación de Reilly, Wood y Robinson en el cálculo de coeficientes de actividad de los electrólitos mezclados y el modelo de Pitzer para e1aluar las propiedades de los sistemas binarios involucrados. Se dan valores de los parámetros de los modelos mencio-nados para los sistemas NaN0 3 - NaCI- H2 0. KCI - NaCI-H2 0, KCI- KN03 - H 2 0 y NaN0 3 - KN03 -H2 0. a 25°C. Los resultados obtenidos son satisfactorios respecto a la información experimental disponible, excepto para el úl timo sistema considerado. SUMJ\JAR Y: An analyt ical method has be en developed for predicting the solid-liquid equilibrio of aqueous solutions con-taining two 1:1 electro~l'tes with a common ion , using the equation of Reil(v, Wood and Robinson to determine the acti-vity coefficients o[ the mixed electrolytes. Bina1y s_vstems ín-formation required for tllis equation is obtained from the mo-del of Pitzer. Parameters of bot/1 equations are given for the systems NaN03 - NaC/- H2 0. KCI - NaCI- H 2 O, KCI- KN0 3 - Hz D and NaN03 - KN03 - H 2 O, at 25°C. Predictions are in good agreement with available experi-mental data, except for the last system above. * Manuscrito revisado y aprobado en forma definitiva en Agosto de 198 J. ** Departamento de Ingeniería Química, Facu ltad de lngenierÍ<t, Universidad de Santiago de Chile, (ex-Universidad Técnica del Estado). 17 INTRODUCCION En el estudio de separación de sales por cristalización fraccionada, es indispensable conocer las curvas de saturación de los electrólitos presentes. Con el objeto de disponer de un método analítico para dicho estudio, se han empleado algunos modelos de predicción, con resultados relativamente satisfactorios, como es el caso del modelo de Correa. Sin embargo, por tratarse de una correlación empírica cuya aplicación se garantiza para sistemas simples, se han .probado otros modelos semi-empíricos que permiten estudiar sistemas más complejos. En una publicación reciente1 se analizó el equilibrio sólido-líquido de soluciones acuo-sas de dos electrólitos 1:1 con un ion común, a 25°C utilizando el modelo de Reilly, Wood y Robinson2 , el cual requiere de información binaria experimental o predicha por algún mo-delo dado. Para tal efecto, se recurrió a la ecuación de Bromley3 , debido a su sencillez y precisión. Los resultados no mejoraron ostensiblemente respecto a los publicados por Correa, Vega y Flores4 , lo cual se atribuyó al uso de la ecuación de Bromley, a pesar que su único parámetro se reevaluó para hacerlo extensivo a la región de sobresaturación. En el presénte trabajo se emplea el modelo de Reilly, Wood y Robinson en la predic-ción del equilibrio sólido-líquido de sistemas ternarios. Para la información binaria requerida, se usa el modelo de Pitzer5 , porque su intervalo de aplicación es más amplio que el de Brom.ley y además considera interacciones iónicas de rango largo así como aquellas de rango corto. MODELOS USADOS Modelo de R~illy-Wood-Robinson (R.W.R.) Para soluciones con dos electrólitos 1: 1, con un ion común, los coeficientes de activi-dad iónico medio de los electrólitos presentes pueden calcularse~a partir del modelo de R.W.R. mediante la ecuación 1 : 1n /' · = 1n /'~ +Y. (if>~ - if>0 ) + 0.5 Y. m g .. + Y. m - -[ o g¡j] 1 1 1 1 1 J IJ 1 0 m (1) donde, /'¡, /'r = Coeficientes de actividad iónico medio del electrólito i en la solución ternaria a su molalidad m¡ y en la solución binaria a la molalidad total m, respectivamente. Y., Y. 1 J Coeficientes osmótico medio de las soluciones binarias de i y j , a la molalidad total m. Fracción molar de i y j en la solución ternaria, defmida como Y. = m . 1m, Y.= m. 1m 1 1 J J en que m = m. + m. 1 J gij Parámetro que considera las interacciones de igual carga. Se ha demostrado experimentalmente5 que el parámetro g¡j varía linealmente con la molalidad total en soluciones de dos electrólitos 1:1 , con un ion común (para molalidades hasta 6): g¡_¡ = 2T + F m (2) Reemplazando g¡j y su derivada en la ecuación (1), se llega a: ln /'· = ln /'0 +Y. (if>0 - if>1?) +Y. m [T + 0.5 (1 +Y.) F m] 1 1 J J J 1 (3) 18 siendo T y F parámetros empíricos. Modelo de Pitzer Los valores de 'Yf, cp f y ct>j de la ecuación (3) se obtienen del modelo de Pitzer, a la molalidad total m: y en que, In 'Y~ = f. + m B. + m2 C. 1 1 1 1 1 - cjJ0 = f. '+ m B. ' + m2 C. 1 1 1 1 f.'= 1 A 3 A 3 + - 1 n (1 + b m0 · 5 ) [ mo.s 2 J 1 + b m0 · 5 b [ mo.s J 1 + b m0 ' 5 (4) (5) (6) (7) B1 = 2 B0 ¡ + (2 B¡¡ / a2 m) [1- (1 +a m0 · 5 - 0.5 a2 m) exp (-a mo.s)] (8) donde , A = 1.17202 para electrólitos en agua, a 25°C. b 1.2 para electrólitos en agua, a 25°C. a 2 para electrólitos 1:1 y 1:2. B0 i' B¡¡, C¡ = Parámetros empíricos de electrólito i. Pitzer entrega valores de los parámetros empíricos para diversos electrólitos a 25°C6 , utilizables hasta 6 molal. Puesto que las ecuaciones (4) y (5) deben usarse a la molalidad total y ella es, en la mayoría de los casos, superior al límite de aplicación, se hace necesario reevaluar los parámetros para asegurar su empleo a concentraciones mayores. Para ellos, se recurrió a la información de Sangster y Lenzi7 sobre actividad de agua a molalidades de sobresaturación. EQUILIBRIO SOLIDO-LIQUIDO La condición de equilibrio para sistemas de dos electrólitos 1: 1 con un ion común, saturados en el componente i, toma la forma: 1 (m.* 'Y-*)2 = m.(m.+m.)-y.2 1 1 1 1 J 1 (1 O) en que el asterisco corresponde al sistema binario saturado en i y los ténninos restantes, al ternario saturado en el mismo componente. Una ecuación análoga puede plantearse para el sistema saturado en j. Para predecir la curva de equilibrio del ternario saturado en i debe resolverse simultá-neamente las ecuaciones (3) y (10). El cálculo de las composiciones de equilibrio en el temario se realiza mediante un proceso iterativo. Dado un valor de m¡, se itera con m_¡ en las ecuaciones ya citadas. Los parámetros T y F de la ecuación (3) pueden determiriarse usando diferentes alternativas. 1 Una de ellas se basa en los valores dados por Pitzer. 8 La otra alternativa utiliza 19 la ecuación (3) reordenada: (11) Esta relación corresponde a la ecuación de una recta de parámetros t y F. En ella, 'Y¡ se calcula con la ecuación (1 O) usando datos experimentales de saturación; -yf, 4> i y 4> j se obtienen del modelo de Pitzer (ecuaciones 4 y 5). En esta segunda alternativa, T y F pueden determinarse por un ajuste de mínimos cuadrados de los valores de las variables de la ecuación (11) (Variante 1). Otro método posible es suponer despreciable el valor de F, dado que es muy pequeño, y calcular T como el promedio aritmético de los valores del primer miembro de la relación (11) (Variante 11) o usar datos del punto de doble saturación (Variante III). En este último caso, se obtienen dos valores de T (T 1 y T z), uno de ellos mediante la ecuación (11) para el componente i y el otro, aplicando la expresión equivalente para j. RESULTADOS Y DISCUSION El método de predicción propuesto se aplicó a los siguientes sistemas, a 25°C: l. NaN03 (2) - NaCl (3) - Hz O (1) 2. KCl (2) - NaCl (3) - Hz O (1) 3. KCl (2) - KN03 (3) - H2 O (1) 4. NaN03 (2) - KN03 (3) - Hz O (1) La información experimental de estos sistemas aparece resumida en trabajos anterio-res.9 , 1 o Valores de los parámetros del modelo binario de Pitzer, a 25°C, aparecen en la Tabla l. Los parámetros T y F dados por Pitzer8 y los valores evaluados por las diferentes variantes se dan en las Tablas 2 y 3. TABLA 1 Parámetros del modelo de Pitz.er, a 25°C. SISTEMA Bo B¡ e NaN03 - Hz O 0.0068 0.1783 -0.00108 NaCl- Hz O 0.0765 0.2664 0.00191 KC1-H2 0 0.0484 0.2122 -0.00125 KN03 - Hz O -0.0816 0.0494 0.00990 TABLA2 Parámetros del modelo de RWR publicados por Pitzer a 25°C. SISTEMA T F NaN03 - NaCl - Hz O 0.016 -0.006 KCl - NaCl- H2 0 -0.012 -0.0018 KCl- KN03 - Hz O 0.016 -0.006 NaN03 - KN03 - Hz O -0.012 -0.0012 20 TABLA 3 Parámetros del modelo de RWR obtenidos por las variantes 1, 11 y III. Sistema Variante 1 Variante li Variante m T F T T¡ T2 NaN0 3 (2)- NaCl (3)- H2 O (1) -0.084 0.0068 -0.027 -0.030 -0.030 KCl (2)- NaCI (3)- H2 O (1) -0.025 0.0009 -0.021 -0.020 -0.019 KCl (2)- KN03 (3) - H2 0 (1) 0.036 -0.0103 -0.012 -0.022 -0.012 NaN03 (2)- KN0 3 (3)- H2 O (1) -0.032 0.0009 -0.023 -0.029 -0.035 Las predicciones de las curvas de equilibrio, a 25°C, de los sistemas temarios conside-rados se pueden apreciar en las Figuras 1 a 5. De la observación de las Figuras 1 y 2, correspondiente al sistema 1, se concluye que al usar los parámetros T y F dados por Pitzer, la predicción es bastante deficiente. En cambio, las tres variantes describen las curvas de solubilidad en forma satisfactoria. La Figura 3 (sistema 2) muestra que todas las alternativas propuestas dan resultados altamente califica-dos. Del análisis de las Figuras 4 y 5 (sistema 3) se observa que el mejor resultado se logró con los parámetros dados por Pitzer. Las tres variantes dan predicciones satisfactorias con sólo pequeñas diferencias entre ellas. Respecto al sistema 4, ninguna alternativa fue capaz de dar una predicción, pues no se consiguió convergencia en la resolución simultánea de las ecuaciones (3) y (1 0). A través del análisis recién realizado, se puede afirmar que para los primeros tres sistemas, el método propuesto pernüte hacer predicciones aceptables, sin tener gran influen-cia las alternativas descritas. El hecho que no se obtengan resultados para el sistema 4, posiblemente se debe a dos causas. Primero, que las expresiones para calcular las propiedades binarias no sean apropiadas cuando se usan en la región de sobresaturación; y segundo, que la ecuación (2) para la función g¡j no sea lineal en m para el sistema 4. En general, las predicciones obtenidas son mejores que las del trabajo publicado ante-riormente1 , excepto para el sistema 4 cuyas curvas de solubilidad no fue posible calcular. 21 e EJC.ptl 2SOC - Cllc:d. con param P1tz.er • • 002 000 ~~--~~~~~~--~~--~-4·--~~ 0.18 0.14 0.,12 o 10 x2 0 .08 0.06 0.04 0.02 0.00 000 000 22 0 .02 0.02 0.04 0 .06 XJ 008 o 10 0.12 • Exptl, 25 OC - · •• Variante 1 - --Vlfilnte 11 - - . Variante 111 0.08 Figura 1. Predicción de las curvas de satu-ración del sistema Na N O 3 (2)-NaCI(J)-H2 0(1), a 25°C, usando parámetros dados por Pitzer. Figura 2. Predicción de las curvas de satu-ración del sistema NaN0 3 (2)-NaCI(J)- H2 0(1), a 25°C, usan-o.12 do las variantes 1,11 y 111. 0.08 [\ 0.07 0.06 0.05 x2 0.04 0.03 0.02 0.01 0.00 0.00 0.02 0.04 0.08 ---.. ------. .. 0.07 006 0.05 0.04 0.03 0.02 • Exptl. 250C _ Calcd con parám. de Piuer 0.01 • Exptl . 25°C _ ParWn. Pitzer v Variantes 1.11 y lll 0.06 0.08 o.oo L_¡___J_..J....._¡__L__¡___J__L__¡____!._.L.-...J....~.....J 0.00 0.01 0 .02 0.03 o 04 0.05 0.06 x, Figura 3. Predicción de las curvas de satu-ración del sistema KCI(2)-NaCI-(3)- H2 0(1), a 25°C. Figura 4. Predicción de las curvas de satu-ración del sistema KCI (2)-KN03 (3)-H2 0(1), a 25°C, usando pará-metros dados por Pitzer. 23 0.07 0.06 1--o.oo -0.04 -0.03 1-0.02 0.01 r-0.00 0.00 • bptl. 25 oc ---- - Veriante 1 _- - Variante 11 - · - · - Variante 111 1 0.01 BIBLIOGRAFIA 1 0.02 1 1 1 1 _l 0.03 0.04 XJ 1 1 -----Figura 5. Predicción de las curvas de satu-ración del sistema KCI(2)-KN0 3 (3)-H 2 0(1). a 25°C, usando las variantes, 1, 11 y 11 l. l. CORREA, H., VEGA, R. y FLORES, F., "Equilibrio sólido-líquido en soluciones acuosas de electró-litos fuertes. l. Uso de modelos de Rcilly-Wood-Robinson y de Bromlcy", Contribuciones, 44, 5- 13 (1980). 2. RElLLY, P.J ., WOOD, R.H. y ROBI SO , R.A., " Prcdiction of osmolic and act ivity cocfficients in mixcd-elcctrolytc solutions", J . Phys. Chem., 76, 9, 1305-1 314 (1971). 3. BROMLEY, LA. , "Thermodynamics properties of strong electro! y tes in aqueous solutions", Al CHE J., 19, 2, 313-320 (1973). 4. 'LORREA , H., VEGA, R. y f LORES, f. , "Métodos para predecir propiedades termodinámicas de soluciones de electrólitos fuertes", Contribuciones, 32, 5-16 (1978). 5. PITZER, K.S. , "Thermodynamics of eleetro lytes. l. TheoreticaJ basis and general equations", J. Phys. Chem., 77, 2, 268-277 (1973) . 6. PITZER, K. S. y MA YORGA, G., "Thermody namics of elcctrolytes. 11. Aetivi ty and osmotic cocl~ ficien ts for electrolytes with one or both ions univalent", J. Phys. Che m. 77, 19, 2300-2308 ( 1973). 7. SANGSTER, J ., TENG, T.T. y LENZI, f., "A general method of calculating thc water activity of super saturated aqueous solu tions from ternary data", Can. J . Chem. , 51, 2, 2626-263 1 (1973). 8. PITZER, K. S. y KIM, J.J. , "Thermodynamics of electrolytes. IV. Activity and osmotie coefficien ts fo r mixed electrolytes", J. Am. Chem. Soc. , 96, 18, 5701 -5707 (1974). 9. CORREA, l-1. , y VERA, J.l-1. , "Coeficientes de actividad de electrólitos fuertes en soluciones ..:onccn-tradas acuosas", Contribuciones, 19, 5-2 1 (1975). LO. CORREA, H., VERA, J.H . y VEGA, R., "Equilibrio sólido-l íquido de soluciones acuosas de dc<.:tró-litos fuertes". Contribudones 25, 5-.13 ( 1977). 24 

Equilibrio sólido -líquido en soluciones acuosas de electrólitos fuertes. II. uso de los modelos de reilly- wood- robinson y de pitzer

http://www.revistas.usach.cl/ojs/index.php/contribuciones/article/download/4190/26003253/

Equilibrio sólido -líquido en soluciones acuosas de electrólitos fuertes. II. uso de los modelos de reilly- wood- robinson y de pitzer

Abstract

Similar works

Full text

Available Versions

Contribuciones Científicas y Tecnológicas