Mezclas líquidas TERMODINÁMICA AVANZADA ! Funciones de exceso en mezclas binarias Expansión de Wohl y ecuaciones der
Views 604 Downloads 53 File size 2MB
Mezclas líquidas
TERMODINÁMICA AVANZADA
!
Funciones de exceso en mezclas binarias
Expansión de Wohl y ecuaciones derivadas
Unidad III: Termodinámica del Equilibrio
• La expansión de Wohl es un método general para representar gE aplicable a mezclas multicomponentes. • Tiene como vantaja el hecho de poder atribuir algún significado físico a cada uno de sus parámetros.
Ecuaciones para el coeficiente de actividad ! Funciones de exceso para mezclas multicomponentes !
9/27/10
• La expansión de Wohl presenta gE de una solución binaria como una serie de potencias en z1 y z2, las fracciones efectivas volumétricas de los componentes:
Rafael Gamero
1
9/27/10
Mezclas líquidas !
!
Expansión de Wohl y ecuaciones derivadas
La expansión de Wohl contiene dos tipos de parámetros: q y a: • q: volumen efectivo de las moléculas (medida del tamaño de la molécula).
+ 4a1222 z1z23 + 6a1122 z12 z22 + ...
9/27/10
z1 "
x1q1 x1q1 + x 2q2
Rafael Gamero
!
z2 "
• a: parámetros de interacción cuyo significado físico es de alguna forma similar al de los coeficientes viriales.
x 2q2 x1q1 + x 2q2 3
!
Funciones de exceso en mezclas binarias
Expansión de Wohl y ecuaciones derivadas
gE = 2a12 z1z2 + 3a112 z12 z2 + 3a122 z1z22 + 4a1112 z13 RT(x1q1 + x 2q2 )
Fracciones efectivas volumétricas de los componentes
2
Mezclas líquidas
Funciones de exceso en mezclas binarias
!
Rafael Gamero
9/27/10
Rafael Gamero
4
Mezclas líquidas !
Mezclas líquidas
Expansión de Wohl y ecuaciones derivadas
!
Ecuaciones de Van Laar
Ecuaciones de Van Laar
Sea una mezcla binaria con componentes no muy diferentes químicamente, salvo por el tamaño de sus moléculas. Si la expansión da Wohl se trunca hasta a12: E
• La razón entre A' y B' es la misma que entre los volúmenes efectivos q1 y q2, y también es igual a la relación entre ln!1" y ln!2".
g 2a x x q q = 12 1 2 1 2 RT x1q1 + x 2q2
ln "1 =
A' # A' x1 & 2 %1+ ( $ B' x 2 '
ln "1 = ln " 2 =
9/27/10
Expansión de Wohl y ecuaciones derivadas
B' # B' x 2 & 2 %1+ ( ! $ A' x1 '
• Las ecuaciones de ln!1 y ln!2 contienen dos parámetros, mientra la ecuación de gE/RT parece ser una ecuación de tres parámetros. • Los valores de A' y B', pueden ser estimados empíricamente.
Donde: A' = 2q1a12 B' = 2q2a12
Rafael Gamero
5
9/27/10
Rafael Gamero
6
!
!
Mezclas líquidas !
Mezclas líquidas
Expansión de Wohl y ecuaciones derivadas
!
Ecuaciones de Van Laar
Expansión de Wohl y ecuaciones derivadas
Ecuaciones de Van Laar
• No es posible encontrar un valor del coeficiente de interacción a12, salvo que se considere un valor arbitrario de q1 o de q2.
• Las suposiciones aplicadas en lsa ecuaciones de Van Laar sugieren que éstas deben ser válidas para soluciones de fluidos relativamente simples, preferiblemente no polares.
• Para aplicaciones prácticas, no es necesario conocecer los valores independientes de q1 e q2, debido a que sólo la relación entre ellos es importante.
• Sin embargo, se ha demostrado empíricamente que estas ecuaciones son capaces de representar los coeficientes de actividad de mezclas complejas. • En estas mezclas, A' y B' dejan de tener un significado físico específico, volviéndose términos empíricos.
9/27/10
Rafael Gamero
7
9/27/10
Rafael Gamero
8
Mezclas líquidas !
Mezclas líquidas
Expansión de Wohl y ecuaciones derivadas
Ecuaciones de Van Laar
Ecuaciones de Margules
• Un caso especial donde las constantes de Van Laar A' e B' son iguales, las ecuaciones de Van Laar son idénticas a las ecuaciones de Margules de dos sufijos.
ln "1 =
A' # A' x & 2 1 %1+ ( $ B' x 2 '
A 2 ln "1 = x2 RT
ln " 2 =
Margules de dos sufijos
! !
B' # B' x & 2 2 %1+ ( $ A' x1 '
Sea una solución binaria cuyos componentes no son muy diferentes en tamaño molecular. La expansión de Wohl con q1 = q2:
ln "1 = A'x22 + B'x23 + C'x24
A 2 ln " 2 = x1 RT Margules de dos sufijos
Rafael Gamero !
3 8 ln " 2 = (A'+ B'+2C')x12 # (B'+ C')x13 + C' x14 2 3
!
!
9/27/10
Expansión de Wohl y ecuaciones derivadas
!
9
9/27/10
Rafael Gamero
10
!
Mezclas líquidas !
Mezclas líquidas
Expansión de Wohl y ecuaciones derivadas
!
Ecuaciones de Margules
Expansión de Wohl y ecuaciones derivadas
Ecuaciones de Margules Las ecuaciones de Margules truncada después de dos términos:
Los coeficientes A’, B’ y C’ correspondientes se definen entonces de la manera siguiente:
Truncadas hasta el término cúbico (Si C’ = 0): Ecuaciones de Margules de tres sufijos (con dos parámetros)
A'= q(2a12 + 6a112 " 3a122 + 12a1112 " 6a1122 ) B'= q(6a122 " 6a112 " 24a1112 " 8a1222 + 24a1122 )
Truncadas hasta el término a la cuarta potencia: Ecuaciones de Margules de tres sufijos (con tres parámetros)
! C'= q(12a1112 + 12a1222 "18a1122 ) ! 9/27/10 !
Rafael Gamero
11
9/27/10
Rafael Gamero
12
Mezclas líquidas !
Mezclas líquidas
Expansión de Wohl y ecuaciones derivadas
!
Ecuaciones de Scatchard-Hamer
Expansión de Wohl y ecuaciones derivadas
Ecuaciones de Scatchard-Hamer
Una forma más de la expansión de Wohl, si se considera la siguiente relación:
A'= v1 (2a12 + 6a112 " 3a122 ) Con:
Es:
q1 v1 = q2 v 2
B'= v1 (6a122 " 6a112 )
ln "1 = A' z22 + B' z23
El comportamiento de las ecuaciones de de Scatchard-Hamer pueden ser consideradas como un intermedio ! entre el comportamiento de las ecuaciones de Van Laar y Margules de tres sufijos.
#v & # 3 v & ln " 2 = % A'+ B'(% 2 (z12 ) B'% 2 (z13 $ 2 '$ v1 ' $ v1 '
Donde: v1 y v2 son los volumes molares de los líquidos puros!1 y 2 a temperatura de la solución.
!9/27/10
!
Rafael Gamero
13
9/27/10
Rafael Gamero
14
!
Mezclas líquidas !
Mezclas líquidas
Ecuación de Wilson
!
Ecuación de Wilson Fracción molar global : x1 = x 2 = 1/2
Las fuerzas intermoleculares
"Moléculas de 2 alrededor de% # & $una molécula central 1 ' x 21 = "Moléculas totales alrededor de% # & $una molécula central 1 '
!
!
Basicamente, el concepto de composión local establece que la composición del sistema en las vecindades de una molécula dada no es igual a la composición 'bulk', debido a las fuerzas intermoleculares.
! 9/27/10
Rafael Gamero
15 de tipo 2
Fracciones molares locales :
La ecuación de Wilson (1964) es un modelo que relaciona a gE con las fracciones molares basado parcialmente en consideraciones moleculares mediante un concepto de composición local.
15
! !
15 de tipo 1
9/27/10
x 21 + x11 = 1
x12 + x 22 = 1
x11 " 3/8
x 21 " 5 /8
! !
Rafael Gamero
16
Mezclas líquidas !
Mezclas líquidas
Ecuación de Wilson
Los términos del modelo:
!
Ecuación de Wilson
Los términos del modelo:
La energía libre de Gibbs de exceso:
gE = "x1 ln(x1 + #12 x 2 ) " x 2 ln(x 2 + # 21 x1 ) RT
Los coeficientes de actividad:
% $12 $ 21 ( ln "1 = #ln(x1 + $12 x 2 ) + x 2 ' # * & x1 + $12 x 2 $ 21 x1 + x 2 )
El ! valor de gE está definido en relación a una solución ideal de LewisRandall, por lo tanto, esta ecuación obedece a la condición límite:
!
% $12 $ 21 ( ln " 2 = #ln(x 2 + $ 21 x1 ) # x1' # * & x1 + $12 x 2 $ 21 x1 + x 2 )
gE = 0 cuando x1 ó x2 =0 9/27/10
Rafael Gamero
17
!
9/27/10
Mezclas líquidas !
!
& % $% ) v1 exp($ 12 11 + ' v2 RT *
9/27/10
!
" 21 #
Rafael Gamero
!
Ecuación de Wilson
Características de la ecuación de Wilson:
Los parámetros de las ecuaciones::
• El modelo de la ecuación de Wilson es ampliamente utilizado como una buena representación de la energía libre de Gibbs de exceso para una gran variedad de mezclas, especialmente para soluciones de compuestos polares o con tendencia a asociación en solventes no polares, para las cuales las ecuaciones de Van Laar y Margules no son suficientes.
Los parámetros de las ecuaciones de Wilson están relacionados con los volúmenes molares de los líquidos y la diferencia de energías características (#ij).
"12 #
18
Mezclas líquidas
Ecuación de Wilson
Los términos del modelo:
Rafael Gamero
& % $ %22 ) v2 exp($ 21 + ' v1 RT *
• Puede ser fácilmente extendido a mezclas multicomponentes. • Desventaja: No predice miscibilidades parciales, por tanto no aplica para equilibrio líquido-líquido. 19
9/27/10
Rafael Gamero
20
Mezclas líquidas !
Mezclas líquidas
Ecuación de Wilson
!
Ecuación NTRL
Características de la ecuación de Wilson:
NTRL (Non-random two liquids):
• Los parámetros de Wilson se reportan para diversos sistemas binarios, ternarios y cuaternarios en:
• El modelo NRTL es, al igual que la ecuación de Wilson, basado en el concepto de composición local, pero si es aplicable a sistemas de miscibilidad parcial.
Gmehling J., U. Onken, W. Arlt (1981). Vapour-Liquid Equilibrium Data Collection, Chemistry Data Series, DECHEMA, Frankfurt/Main, Germany.
9/27/10
Rafael Gamero
# " G gE "12G12 & = x1 x 2 % 21 21 + ( RT $ x1 + x 2G21 x 2 + x1G12 '
G12 = exp("#12$12 )
!
21
9/27/10
g21 # g11 RT 22
Mezclas líquidas
Ecuación NTRL
!
• El significado de los elementos gij es similar al de los elementos #ij de la ecuación de Wilson: parámetros de energía característicos de interacciones i-j.
Ecuación NTRL
Los coeficientes de actividad se puedn calcular de acuerdo a; 2 * $ G21 ' #12G12 -/ , ln "1 = x # 21& ) + 2 ,+ % x1 + x 2G21 ( (x 2 + x1G12 ) /. 2 2
• El parámetro $12 está relacionado con la no-aleatoriedad de la mezcla, o sea que los componentes de la mezcla no se distribuyen aleatoriamente (es decir, uniformemente), sino que siguen un patrón definido por la composición local.
!
• Cuando $12 = 0, la mezcla es completamente aleatoria y al ecuación se reduce a la ecuación de Margules de dos sufijos. Rafael Gamero
" 21 =
!
!
Mezclas líquidas
9/27/10
g12 # g22 RT
Rafael Gamero
!
!
! G21 = exp("#12$ 21 )
"12 =
23
9/27/10
2 * $ G12 ' # 21G21 -/ , ln " 2 = x #12 & ) + 2 ,+ % x 2 + x1G12 ( (x1 + x 2G21 ) /. 2 1
Rafael Gamero
!
24
Mezclas líquidas !
Mezclas líquidas
Ecuación NTRL
!
Ecuación UNIQUAC
• Para sistemas ideales o moderadamente ideales, NRTL no ofrece muchas ventajas sobre Van Laar o Margules de tres sufijos, pero para sistemas fuertemente no ideales, esta ecuación puede representar una buena simulación de datos experimentales.
• Debido a que no es sencillo obtener datos de buena calidad para estimar los tres parámetros de NRTL, se han efectuado varios intentos por desarrollar un modelo para gE con dos parámetros y con las ventajas de la ecuación de Wilson, pese a estar limitada a la miscibilidad total.
• Se requieren datos experimentales de buena calidad para la estimación de los tres parámetros.
• Abrams y Prausnitz (1975) desarrollaron una ecuación conocida como Teoria Quasi-química Universal, o UNIQUAC, por las siglas en inglés.
• NRTL también puede ser fácilmente extendido para mezclas multicomponentes.
• Este modelo es una especie de extensión de la teoría quasi-química de Guggenheim (1952) para moléculas no-aleatorias en mezclas que concomponentes de diferentes tamaños.
9/27/10
Rafael Gamero
25
9/27/10
Mezclas líquidas !
!
• La ecuación UNIQUAC consiste en dos partes:
Ecuación UNIQUAC
• Para una mezcla binaria:
• Parte combinatoria: Contribución entrópica de los componentes
E gcomb "* "* z $ # # ' = x1 ln 1 + x 2 ln 2 + &q1 x1 ln 1* + q2 x 2 ln 2* ) RT x1 x2 2 % "1 "2 (
• Parte residual: Fuerzas intermoleculares responsables de la entalpía de la mezcla
!
E gres = "q'1 x1 ln(# '1 +# '2 $ 21 ) " q'2 x 2 ln(# '2 +# '1 $12 ) RT
!
g E = g E (combinatoria ) + g E (residual ) Rafael Gamero
26
Mezclas líquidas
Ecuación UNIQUAC
9/27/10
Rafael Gamero
27
!
9/27/10
Rafael Gamero
28
Mezclas líquidas !
Mezclas líquidas
Ecuación UNIQUAC
!
Para cada mezcla existen dos parámetros ajustables !12 y !21, que se definen como:
Donde el número de coordinación z = 10 y además: Fracciones de segmento
!
Fracciones de área
"1* =
x1r1 x1r1 + x 2 r2
"1 =
x1q1 x1q1 + x 2q2
" '1 =
x1q'1 x1q'1 + x 2q'2
"*2 =
x 2 r2 x1r1 + x 2 r2 !
"2 =
x 2q2 x1q1 + x 2q2 !
" '2 =
x 2q'2 x1q'1 + x 2q'2
9/27/10
!
!
Ecuación UNIQUAC
Rafael Gamero
% $u ( % a ( "12 = exp'# 12 * + exp' # 12 * & RT ) & T)
! 29
Donde %u12 y %u21 son las energías características de interacción, las ! cuales son dependientes da temperatura. 9/27/10
Ecuación UNIQUAC ln "1 = ln
!
9/27/10
!
30
Mezclas líquidas
Las ecuaciones para el coeficiente de actividad son:
!
& #1* z $ r ) + q1 ln 1* + #*2 (l1 % 1 l2 + x1 2 #1 r2 * '
Ecuación UNIQUAC
Los términos l1 y l2 de las cuaciones para el coeficiente de actividad del modelo UNIQUAC están definidos como:
& ) , 21 ,12 % q'1 ln($ '1 +$ '2 , 21 ) + $ '2 q'1 ( % + ' $ '1 +$ '2 , 21 $ '2 +$ '1 ,12 *
ln " 2 = ln
Rafael Gamero
!
Mezclas líquidas !
% $u ( % a ( " 21 = exp'# 21 * + exp' # 21 * & RT ) & T )
z l1 = (r1 " q1 ) " (r1 "1) 2
UNIQUAC es aplicable a una gran variedad de soluciones líquidad no
& #*2 z $ r ) + q2 ln 2* + #*2 (l 2 % 2 l1 + x2 2 #2 r1 * '
& ) ,12 , 21 % q'2 ln($ '2 +$ '2 ,12 ) + $ '1 q'2 ( % + ' $ '2 +$ '1 ,12 $ '1 +$ '2 , 21 * Rafael Gamero
z l2 = (r2 " q2 ) " (r2 "1) 2
polares o no polares incluyendo ! electrolíticas que contengan componentes !
sistemas de miscibilidad paracial. Gmehling et al (1981) reporta parámetros UNIQUAC para una variedad de sistemas.
31
9/27/10
Rafael Gamero
32
Mezclas líquidas !
Mezclas multicomponentes
Funciones de exceso para mezclas multicomponentes
!
Ecuación de Wilson
La ecuación de Wilson extendida a mezclas multicomponentes está definidos como:
• Para los casos prácticos más comunes, las funciones de exceso son aplicables a sistemas multicomponentes.
" ij =
% $ # $ii ( vj exp' # ij * vi RT ) &
" ji =
% $ # $ jj ( vi exp' # ji * vj RT ) &
La energía libre de Gibbs de exceso:
% ( g = "# x i ln'# x j $ ij * RT '& j *) i E
• Las ecuaciones resultantes de las funciones de exceso para mezclas líquidas multicomponentes son extensiones de aquellas para sistemas binarios. • En algunos casos, existen ciertos coeficientes adicionales agregados a las ecuaciones para su extensión a mezclas multicomponentes. 9/27/10
Rafael Gamero
! & ) x i $ ik ln " k = #ln(% x j $ kj + + 1# % (' j +* i % x j $ ij j !
! 33
9/27/10
Coeficientes de actividad para cualquier componente k
Rafael Gamero
34
!
Mezclas multicomponentes !
Mezclas multicomponentes
Ecuación NRTL
!
La ecuación NRTL extendida a mezclas multicomponentes está definidos como:
Ecuación UNIQUAC
Fracciones de segmento
La ecuación UNIQUAC extendida a mezclas multicomponentes es la suma de dos partes:
"*i =
ri x i # rj x j j
gE = " xi RT i
"#
jiG ji x j
j
"G
ki x k
G ji = exp("# ji$ ji ) con : (# ji = # ij )
" ji =
g ji # gii RT
E gcomb #* z $ = " x i ln i + " qi x i ln i* ! RT x 2 # i i i i
k
$# ln " i =
!
!
G!ji x j
$G k
9/27/10
ji
j ki
xk
& ! $x # G ) r rj rj + x j Gij ( r +$ ( # ij % $ Gkj x k ++ j $ Gkj x k ( ' * k k Rafael Gamero
Fracciones de área
"i =
qi x i #q j x j j
& ) g = "# q'i x i ln(( #$ ' j % ji ++ RT i ' j * E res
! 35
9/27/10
!
Rafael Gamero
" 'i = !
!
q'i x i # q' j x j j 36
Mezclas multicomponentes !
Mezclas multicomponentes
Ecuación UNIQUAC
!
El número de coordinación z =10 y el coeficiente de actividad para cualquier componente i se define:
ln " i = ln
j
%# '
k
$ kj
lj =
z (rj " qj )" (rj " 1) 2
Rafael Gamero
!
Grupos contenidos en una mezcla propanol-hexano.
Donde:
#' j k
9/27/10
El método de contribución de grupos (Langmuir, 1925) se aplica una mezcla contiene componentes con moléculas complejas, en las cuales los grupos funcionales tiene diferentes contribuciones.
( + #*i z $ #* + qi ln i* + l i % i & x j l j % q'i ln**&$ ' j ' ji -- + xi 2 #i xi j ) j ,
+q'i "q'i %
!
Contribución de grupos
37
9/27/10
Mezclas multicomponentes
Contribución de grupos •
De acuerdo a la literatura, existen muchas propuestas para definir los coeficientes de actividad para mezclas según los grupos funcionales.
•
El método más frecuente es la combinación del concepto de solución por grupos para la parte combinatoria, de la forma en que se aplica en UNIQUAC para la parte residual. Este método se conoce como UNIFAC (UNIQUAC Functional-group Activity Coefficient).
•
!
Ecuación UNIFAC Un modelo consistente en la suma de contribuciones ya fue descrito antes como la suma de un término combinatorio y un termino residual:
ln " i = ln " icom + ln " ires
!
UNIFAC es un método difícil de aplicar a casos de mezclas simples. Requiere una extensa base de datos de equilibrio líquido-vapor para diversos grupos funcionales. Rafael Gamero
li =
Donde la parte combinatoria es:
ln "
com i
# z $ # = ln i + qi ln i + l i % i & x j l! j xi 2 #i xi j ri = #" k(i ) Rk k
9/27/10
38
!
Mezclas multicomponentes !
Rafael Gamero
39
!
k
9/27/10
!
qi = #" k(i )Qk
"i =
ri xi # rj xj j
"i = !
Rafael Gamero
!
z ( ri " qi ) " ri +1 2
qi xi
#q x
j j
j 40
!
Mezclas multicomponentes !
Mezclas multicomponentes
Ecuación UNIFAC
!
Los parámetros ri y qi son calculados como la suma de los parámetros de área y volumen del grupo: ri = #" k(i ) Rk
!
Rk =
Qk =
Awk 2.5x109
"m =
grupos
ln " ires =
$# [ln% (i) k
k
k
Vwk ! 15.17
Área superficial y volumen de los grupos de Van der Waals
!
9/27/10
!
Rafael Gamero
41
!
9/27/10
!
k
& ln%k(i) ]
Qm X m
#Q X n
n
n
0 ' * / $ m %km! 2 ln"k = Qk /1# ln) & $ m %mk , # & 2 (m + m & $ n %nm 2 /. 1 n $ a ' "mn = exp mn ) % RT (
!
Fracción de área del grupo m:
La parte residual:
qi = #" k(i )Qk
k
Ecuación UNIFAC
Parámetro de interacción energética entre los grupos m y n: Rafael Gamero
Fracción molar del grupo m: M
#" Xm =
(i ) m i
x
i M
N
##" i
(i ) j i
x
j 42
!