UNIVERSIDAD COMPLUTENSE DE MADRID

K
i
=
∏(
m jγ
j )
j
m γ
c c
i i
Apéndice 1. Cálculo de la especiación en soluciones acuosas.
q
ij
(A1.5)
donde Ki es la constante de equilibrio termodinámica de la asociación, γ ic es el
coeficiente de actividad del complejo i, y γj es el coeficiente de actividad del ion j.
El numerador representa la actividad del complejo i, mientras que el denominador
representa el producto de las actividades de todos los iones j1, j2, j3...) que
intervienen en la reacción (A1.4) elevadas a su correspondiente coeficiente
estequiométrico.
Para realizar este cálculo tan complejo es necesario realizar un tratamiento
informático. En este trabajo se ha optado por utilizar programas específicos para
cada uno de los sistemas SS-SA estudiados. Estos programas son similares al
ideado por Lundager Madsen (1984) para sistemas de precipitación de oxalatos y
fosfatos. El programa en cuestión esta basado en un algoritmo de resolución de
ecuaciones por medio de un proceso iterativo, en el que se hace uso del método
de aproximaciones sucesivas de Newton-Ramphson y es una modificación del
diseñado por M. Prieto (Fernández-González, 1996) para el cálculo de
especiación del sistema (Ca,Cd)CO3-H2O. La descripción de cómo funciona el
programa es la siguiente: Se parte de la suposición de que las constantes de los
diferentes equilibrios se obtienen a partir de la ley de acción de masas, en la que
se consideran las actividades iguales a las concentraciones, siendo por tanto los
coeficientes de actividad igual a la unidad. Seguidamente se resuelve un sistema
de ecuaciones que proporciona una fuerza iónica provisional, la cual se empleará
para calcular, mediante la ecuación de Debye-Hückel, unos coeficientes de
actividad también provisionales. El proceso se repetirá hasta que la diferencia
entre los coeficientes de actividad correspondientes a dos iteraciones sucesivas
sea menor que un número que se considere el error máximo admisible. Una vez
llegado a este punto el programa facilitará las actividades de los iones
procediéndose al cálculo de la sobresaturación.
261