UNIVERSIDAD MAYOR REAL Y PONTIFICIA DE
SAN FRANCISCO XAVIER DE CHUQUISACA
FACULTAD DE CIENCIAS Y TECNOLOGIA
INGENIERIA QUIMICA
DISEÑO Y SIMULACION DE PROCESOS INDUSTRIALES I
(PRQ 350)
Tema : Practica N°1
Fecha de la realización de la practica :
Grupo :1
Estudiante: Mamani Gonzales Mayra Dayana
Docente : Ing. Maximiliano Churata
Fecha de la realización de la practica : 06-05-22
Fecha de presentación: 15-07-2022
SUCRE-BOLIVIA
2022
1. OBJETIVOS
Caracterizar un modelo termodinámico para la selección en diseño de procesos.
Evaluar un modelo termodinámico en Excel.
2. FUNDAMENTO
2.1. Principios Termodinámicos
El equilibrio de fases como de liquido y vapor existe cuando no se presenta cambios con
respecto a la temperatura y presión total de la mezcla y las fugacidades parciales de los
componentes puros son igual en todo el sistema se expresa de la siguiente forma
^f = ^f
iV iV
Y términos de coeficiente de fugacidad queda representado por :
^ iV P=x i Φ
yi Φ ^ iL P
Para esta ecuación es necesario saber el coeficiente de fugacidad entre las mas conocidas y
utilizadas son Ley de Gases Ideales, la ecuación de Peng-Robinson y la de Redlich-Kwong-
Soave.
Para la fase liquida, es mas compleja encontrar una ecuación con el coeficiente de
fugacidad que represente el comportamiento de las mezclas por lo que frecuentemente se
recurre por medio del coeficiente de la actividad , y en este caso la ecuación de equilibrio
es:
Ecuación (3)
^ iV P=x i γ i f iL
yi Φ
Es importante mencionar que también existe otra medida de análisis del ELV, y ésta es por
medio del coeficiente de distribución vapor-líquido, también llamado razón de eq|uilibrio
líquido-vapor.
Ecuación (4)
^
yi Φ
K i= = iL
^ iV
xi Φ
Mezcla liquido no ideal -vapor ideal.
Se toma en cuenta que un vapor sea gas ideal cuando la presión del sistema sea menor o a
tres atmosferas por lo que el coeficiente de fugacidad en la saturación y el coeficiente de
fugacidad parcial en la fase vapor serán iguales a uno. Pero determinando el punto de
burbuja para la mezcla alcohol-agua siendo componente polar no se los puede considerar
como ideal por lo que el coeficiente de actividad (γ i) no será igual a 1. (Henley, 1990)
Para evaluar los sistemas no ideales se tiene:
UNIFAC
El método UNIFAC (De sus siglas en inglés UNIQUAC Functional-group Activity
Coefficients, no confundir con UNIQUAC)1 es un sistema semi-empírico para la
predicción de la actividad de no-electrolitos en mezclas no ideales. La UNIFAC utiliza los
grupos funcionales presentes en las moléculas que componen la mezcla líquida para
calcular los coeficientes de actividad. Utilizando interacciones para cada uno de los grupos
funcionales presentes en las moléculas, así como algunos coeficientes de interacción
binaria, se puede calcular la actividad de cada una de las soluciones. Esta información
puede ser usada para obtener información sobre los equilibrios de líquidos, lo cual es útil en
muchos cálculos termodinámicos, tales como el diseño de reactores químicos y en la
destilación.
combinatorio residual
ln ( γ i ) =ln ( γ i ) +ln ( γ i )
combinatorio
ln ( γ i ) =1+ J i+ ln J i−5 qi 1−
( Ji
Li
+ln
Ji
Li )
ln ( γ i )
residual
[ (
=q i∗ 1−∑ θ k
β ik
sk
β
−eik ln ik
sk )]
θ k=
∑ x i qi e ik ; e = v ik ∗Qk ; β =∑ e ∗τ
∑ x j q j ik qi ik mi mk
−amk
sk =∑ Qm∗τ mk ; τ mk =−e T
q i=∑ v k∗¿ Qk ; r i=∑ v k∗¿ R k ¿ ¿
i i
3. DESCRIPCION DE LA PRÁCTICA
3.1. Componentes
Agua
Formula desarrollada Formula
estructural formula tridimensional
El agua es una sustancia cuya molécula está compuesta por dos átomos de hidrógeno y uno
de oxígeno (H2O) unidos por un enlace covalente.2 El término agua, generalmente, se
refiere a la sustancia en su estado líquido, aunque esta puede hallarse en su forma sólida,
llamada hielo, y en su forma gaseosa, denominada vapor.2 Es una sustancia bastante común
en la Tierra y el sistema solar, donde se encuentra principalmente en forma de vapor o de
hielo. Es indispensable para el origen y sustento de la vida.
Metanol
Formula desarrollada Formula estructural formula tridimensional
El compuesto químico metanol, también conocido como alcohol de madera, alcohol
metílico o raramente alcohol de quemar, es el alcohol más sencillo. A temperatura ambiente
se presenta como un líquido ligero (de baja densidad), incoloro, inflamable y tóxico que se
emplea como anticongelante, disolvente y combustible.
3.3. Datos de la práctica
Datos a condiciones normales
CN
Z 10
T (°C) 25
P (atm) 1
ɥ (Kg/m*s) Metanol 0.544
ɥ (Kg/m*s) agua 0.000891
p (Kg/m^3) agua 997
p (Kg/m^3) Metanol 792
Datos de las constantes de metanol y agua
D C Di Ci
Agua 0.581672975 1.534285714 0.379116464 0.65176909
Metanol 0.379116464 0.651769088 0.581672975 1.53428571
Rk agua 0.92 Rk 1.4311
metanol
Qk agua 1.4 Qk 1.432
metanol
3.4. Procedimiento de la Práctica
3.4.1. Planteamiento de la práctica:
Hallar el coeficiente de actividad en una mezcla de H2O y CH3OH en condiciones
normales con un 40% mol de H2O
Con la ecuación UNIFAC que es la siguiente.
C R
ln γ i=ln γ i + ln γ i
Obteniendo los siguientes datos
ln T k 0.007561174
ln T k i 1.31814274
R
ln γ I -0.00090137
Ѱ i, k 0.96078944
ln γ i -2.1870643
γi 0.11224579
Diagrama de proceso realizado
Se realizo la practica de manera presencial y se utilizo las herramientas del Excel como se
ve en el grafico
3.4.2. Descripción de la práctica
Se trabajo con la ecuación de estado UNIFAC crean un programa en Excel que en el
momento de meter las constantes y datos de los componentes nos saque el resultado del
coeficiente de actividad
Se empezó con el análisis de la ecuación UNIFAC viendo que constantes se necesita en
este caso se trabajo con una mezcla de agua y metanol con referencias bibliográficas
hallamos las constantes para poder realizar los cálculos en el programa de Excel
Con datos de las propiedades de los componentes de la mezcla pudimos hallar algunas
variables que necesitábamos para poner en practica la ecuación de UNIFAC
4. CALCULOS
4.1. Datos
θ φ
r q l
A 8.9 0.0 0.0 0.9 0. 0.6
g 368 006 00 97 60 99
u 1E- 327 96 71 54 99
a 07 41 28 66 11 83
7 3 5 7
m 0.0 0.0 0.0 0.9 0. 0.3
et 006 009 00 99 39 00
a 868 842 98 01 45 00
4.2. Cálculos Matemáticos
n 69 57 48 26 88 16
o 8 5 5 3
Cálculo del coeficiente de actividad
l
C
ln γ i
-2.18616292
-2.1870643
Ѱ i, k
ln γ i
0.96078944
Se separa por grupos funcionales a los
componentes y guiándose por las tablas
del libro Eliott- Lira ln T k 0.007561174 pag 540
ln T k i 1.31814274
AGUA ETANOL Rk Qk
H2O vk= 1 vk = 0 0.92 1.4
CH3 vk = 0 vk= 1 0.9011 0.848
CH2 vk = 0 vk= 1 0.6744 0.54
OH vk = 0 vk= 1 1 1.2
q i=∑ v k∗¿ Rk =1∗1.4 +1∗0.848+ 1∗0.54 +1∗1.2=3.988 ¿
i
r i=∑ v k ∗¿ R k =1∗0.92+1∗0.9011+1∗0.6744+1∗1=3.4955 ¿
i
Combinatorio
combinatorio
ln ( γ i )
(
=1+ J i+ ln J i−5 qi 1−
Ji
Li
+ln
Ji
Li )
ri
J i=
∑ r i∗xi
Para el etanol
3.4955
=1.22549
3.4955∗0.816
Para el agua
3.4955
=5.43478
3.4955∗0.184
qi
Li = =¿
∑ q i∗x i
Para el etanol
3.98
=1.22549
3.98∗0.816
Para el agua
3.98
=5.43478
3.98∗0.184
Por lo que:
Li=J i
Remplazando datos en: ln ( γ i )
combinatorio
(
=1+ J i+ ln J i−5 qi 1−
Ji
Li
+ln
Ji
Li )
combinatorio
ln ( γ i ) =−1.52177E-07
Residual
residual
ln ( γ i )
[ (
=q i∗ 1−∑ θ k
β ik
sk
β
−eik ln ik
sk )]
Para el grupo funcional H2O:
i
v k∗Q k 1∗1.4
e ik = = =0.351
qi 3.988
θ k=
∑ x i qi e ik
∑ x jq j
Para el grupo funcional H2O:
0.184∗3.988∗0.351
θ k= =0.0654
3.988∗(0.184+ 0.816)
−amk
T
τ mk =−e
Interacción Binaria
Por lo que: H2O (8) CH3 (1) CH2 (1) OH (5)
H2O (8) 0 342.4a mk 342.4 -122.4
CH3 (1) 1300.1H2O 0 -122.4 0 644.6
CH2 (1) 1300.1CH3 0 644.6 0 644.6
OH (5) 28.73 328.2 328.2 0
CH2 644.6
OH 328.2
Entonces a Temperatura
−amk − (−122.4 k )
T 342.15k
τ mk =−e =−e =−0.699256372
β ik =∑ e mi∗τ mk =0.351053159∗(−0.699256372 )=−0.245476159
sk =∑ Qm∗τ mk =1.4∗(−0.699256372 )=−2.788634412
Por lo que para todos los componentes:
a mk e ki θk τ mk β ik sk
H2O -122.4 0.3511 0.0646 -0.6993 -0.2455 -2.7886
CH3 644.6 0.2126 0.1735 -6.5796 -1.3991 -26.2393
CH2 644.6 0.1354 0.1105 -6.5796 -0.8909 -26.2393
OH 328.2 0.3009 0.2455 -2.6097 -0.7853 -10.4074
Remplazando en:
residual
ln ( γ i )
[ (
=q i∗ 1−∑ θ k
β ik
sk
β
−eik ln ik
sk )]
residual
ln ( γ i ) =−3.55079E-07
combinatorio residual
ln ( γ i ) =ln ( γ i ) +ln ( γ i )
ln ( γ i ) =−1.52177E-07−3.55079E-07
ln ( γ i ) =−5.07255E-07
γ i=0.999999493 ≅ 1
4.3. Resultados
Observando los resultados obtenidos, se puede observar que la temperatura de burbuja más
alta es de la mezcla Que depende mucho las propiedades de T para el resultado final
También se puede observar que los errores de las presiones son de un alto grado ya que es
una ecuación netamente experimiental
5. CONCLUSIONES.
Se logro usar adecuadamente la ecuación de estado UNIFAC, aplicándola a un ejemplo
para encontrar el coeficiente de actividad de metanol– agua, así aprendiendo la importancia
como ingeniería química del uso de estas ecuaciones dadas, para futuro uso profesional.
6. BIBLIOGRAFÍA
henley, 1990. s.l.:s.n.
Henley, E. J. &. S. J. D., 1990. Operaciones de separación por etapas de equilibrio en
ingeniería Química. Barcelona: Reverte S.A..
