Equilibrio agua-aire

Para encontrar las concentraciones de equilibrio del sistema debemos igualar los potenciales qu¡micos.

La ecuaci¢n de potenciales del gas en ambas fases la reemplazamos por la ecuaci¢n de la solubilidad del aire en agua.

$$\left(1-y_w\right)=K P x_a$$ (26)

Igualando los potenciales qu¡micos del agua en ambas fases y reemplazando $y_w$ llegamos a:

$$RT \ln fs = P\left(1-\frac{P_x}{P}\right)^2\left(2 A_{aw}-A_{aa}\right)$$ (27)

$$+A_{ww}P \left(2\left(\frac{P_x}{P}-\frac{P_v}{P}\right)+\frac{P_v}{P}-\left(\frac{P_x}{P}\right)^2\right)$$ (28)

$$+A_{www} P^2 \left(\frac{P_v}{P}\right)^2 \left(\frac32 \left(fs^2-1\right) - \left(1-\frac{P_v}{P}fs^3\right)\right)$$ (29)

$$+V_l P \left(1-\frac{P_v}{P}\right)$$ (30)

$$-RT \ln \left(1-k \left(1-fs \frac{P_v}{P}\right)P\right)$$ (31)

donde

$$P_x=x_w P \quad fs=\frac{P_x}{P_v}$$ (32)

$V_l$ es el volumen del l¡quido y $k$ es la solubilidad del aire en agua.

Esta ecuaci¢n debe ser iterada, pero si hacemos $fs=1$ obtenemos de entrada una primera aproximaci¢n bastante buena, as¡ procede el [7] que escribe la ecuaci¢n con esta simplificaci¢n efectuada y reemplaza el logaritmo por una aproximaci¢n (la edici¢n que tengo tiene algunos errores de tipeo) El Psicro.bas itera pues con muy pocos ciclos converge.