1.
- Determine el coeficiente de fugacidad y la fugacidad a partir de la expresión
G /RT para una mezcla binaria de metano y n- Hexano, a 75C y 8 atm.
E
La energía libre de Gibbs de exceso se define como lo siguiente
G E=G( solucionreal ,T ,P y x )−G( solucionideal a la mismaT , P , X )
Para llegar a esto tenemos que la energía libre de Gibbs molar parcial en función
de la fugacidad se describe por lo siguiente
Ǵ i=Γ i ( T ) + RT ln ^f i
Para un gas ideal
Ǵ iid =Γ i (T ) + RT ln x i f i
La diferencia entre estas dos es precisamente la energía libre de Gibbs de exceso
f^ i
Ǵ i−Ǵi id =RT ln
xi f i
Donde la última fracción se define como el coeficiente de la actividad
f^ i
γ i=
xi f i
Sustituyendo el coeficiente de actividad queda, que
Ǵ i E=RT ln yi
El cual se puede despejar
Ǵi E
=ln y i
RT
Para el calculo de la fugacidad se puede hacer uso de la ecuación virial
2.- Determine el coeficiente de fugacidad y la fugacidad a partir de la expresión
G /RT para una mezcla binaria de metano y n- Hexano, a 75°C y 8 atm.
R
Para la energía libre de Gibbs residual podemos empezar proponiendo lo siguiente
P
dP
ln φ=∫ ( Z−1 )
0 P
Podemos sustituir el coeficiente de fugacidad por la energía libre de Gibbs residual
P
GR dP
=∫ ( Z−1 )
RT 0 P
Se igualan entonces la fugacidad y la energía libre de Gibbs para cualquier
especie
Gi R
ln φ i=
RT
Si bien la energía libre de Gibbs residual en un principio es lo siguiente
Gi R=G real −G ideal
Relacionado con el coeficiente de fugacidad, tenemos
R fi
Gi =RT ln =RT ln φ
P
Para el cual el coeficiente de fugacidad es
fi
φ i=
P
3.- La energía libre de exceso del metil ciclohexano y el tetra hidrofurano a
303.15K sigue la siguiente ecuación:
GE = RTx (1-x) (0.4857 – 0.1077(2x-1) + 0.0191 (2x – 1)2 )
Donde x es la fracción mol del metil ciclohexano. Calcule la energía libre de Gibbs
de la disolución de 1 mol de metil ciclohexano y 3 moles de tetra hidrofurano
Entonces, primero se tiene que calcular la fracción mol del metil ciclohexano, para
esto tenemos la siguiente formula
nmetil ciclohexano
X metilciclohexano =
nT
1 mol
X metilciclohexano = =0.25
1 mol+ 3 mol
ntetra hidrofurano
X tetra hidrofurano=
nT
3 mol
X tetra hidrofurano= =0.75
1 mol+3 mol
Como la x que nos mencionan en el problema es la fracción mol del metil
ciclohexano ya tenemos el dato
x=0.25
Se sustituye x, R que es la constante de los gases ideales y la temperatura
en la ecuación de G E y
G E=( 8.314 J /molK )( 303.15 K ) (0.25)(1−0.25)(0.4857 – 0.1077( 2× 0.25−1)+ 0.0191(2× 0.25 – 1)2)
J
G E=246.2345
mol
Ahora, de G E se puede despejar la Energía libre de Gibbs, ya que :
G E=G−Gid
Entonces despejamos G
G=G E +Gid
Para esto calculamos G id , donde
Gid =RT ¿
x A= la fracción mol del metil ciclohexano
x B =¿la fracción mol del tetra hidrofurano
La a es la actividad, pero como es ideal
a i=x i
Entonces ya podemos sustituir para la G id
J
G id =(8.314 ) ( 303.15 K ) [ (0.25)ln 0.25+ ( 0.75 ) ln 0.75 ]
molK
Gid =−1417.3 J /mol
Susutituimos para calcular la G en la formula que despejamos anteriormente
G=246.2345 J /mol + (−1417.3 J /mol )
J kJ
G=−1171 ó−1.171
mol mol
4.- Goodwin reporta los siguientes datos para el tolueno (J. Phys. Chem.Ref. Data
18, 1565 (1989)) en función del factor de compresibilidad Z, a diferentes presiones
y temperaturas. A partir de estos datos determinar el coeficiente de fugacidad del
tolueno a 600K y : a) 30 bar, b) 1000 bar
a) 30 bar
Para calcular el coeficiente de fugacidad de tiene que
Pr
Z −1
ln φ=∫
P0
[ ]
P
d Pr
Para esto necesitamos la presión reducida
P
Pr =
Pc
La presión critica del tolueno es de 40.8 bar
¯¿
Pr =30 ¿
40.8 ¿ =0.7352 ¿
¯
Ahora para la integral tenemos que
Z −1
ln φ= [ ]P
( ¿ P r−Po ) ¿
Ahora sustituimos con los datos, considerando que Z para 30 bar 0.70117
ln φ=¿ ¿
ln φ=¿−7.3233× 10−3 ¿
Despejamos y que es el coeficiente de fugacidad
−3
φ=e−7.3233× 10
φ=0.9927
b)1000 bar
Ocupamos la misma formula
Pr
Z −1
ln φ=∫
P0
[ ]
P
d Pr
Obtenemos la presión reducida, cuando la presión critica del tolueno es de 40.8
bar
P
Pr =
Pc
¯¿
Pr =1000 ¿
40.8 ¯¿ =24.5 ¿
En la integral tenemos que
Z −1
ln φ= [ ]P
( ¿ P r−Po ) ¿
Y sustituimos, dado que en la tabla el valor de Z para 1000 bar es de 2.48836
ln φ=¿ ¿
ln φ=0.0364
Despejamos y que es el coeficiente de fugacidad
φ=e 0.0364
φ=1.03
