UNIVERSIDAD MAYOR REAL Y PONTIFICIA DE

SAN FRANCISCO XAVIER DE CHUQUISACA

FACULTAD DE TECNOLOGÍA

INGENIERIA DE LA REACCION QUIMICA II (PRQ-207)

CARRERA: INGENIERÍA QUÍMICA

PRACTICA N.º 3: DIAGRAMAS DE EQUILIBRIO LIQUIDO-VAPOR

NOMBRE DEL ESTUDIANTE: TELLEZ JALDIN MARIA FERNANDA
FECHA DE REALIZACION DE LA PRACTICA:25/09/2023
FECHA DE ENTREGA DE LA PRACTICA:10/10/2023
DOCENTE: ING. FRANCISCO JAVIER CAMACHO CALDERÓN

sucre – Bolivia
2023
 DIAGRMAS DE EQUILIBRIO LIQUIDO-VAPOR
1. OBJETIVOS
 Determinar experimental y analíticamente la temperatura de burbuja del
sistema etanol - agua a distintas composiciones.
 Construir los siguientes Diagramas:
1. T burbuja, experimental – x Etanol.
2. T burbuja, teórico – x Etanol.
3. T rocío – x Etanol.
4. T burbuja, experimental – x – y para etanol.
5. T burbuja, experimental – x – y para etanol.
6. x – y para etanol.
7. ºBrix – x Etanol.
8. ºGL – x Etanol.
9. ºGL – ºBrix.
10. IR – X Etanol.
11. Combinar los dos diagramas 1) y 2).
12. Combinar los dos diagramas 1) y 3).
13. Combinar los dos diagramas 2) y 3).
2. MARCO TEORICO
a) Equilibrio de fases:
El equilibrio de fases es el estudio del equilibrio que existe entre o en los
diferentes estados de materia; es decir sólido, líquido y gas. el equilibrio se define
como una etapa cuando el potencial químico de cualquier componente presente
en el sistema permanece estable con el tiempo. la fase es una región en la que la
interacción intermolecular es espacialmente uniforme o en otras palabras las
propiedades físicas y químicas del sistema son iguales en toda la región. dentro
del mismo estado, un componente puede existir en dos fases diferentes, como los
alótropos de un elemento. además, dos compuestos inmiscibles en el mismo
estado líquido pueden coexistir en dos fases.
b) Condiciones de equilibrio.
Se trata de un equilibrio sólido-líquido, por lo tanto, consideraremos sus
potenciales químicos a T y p constantes.
dG=−SdT +Vdp+ ∑ μfase d nfase
fases

En equilibrio dG=0 a T y p ctes

La energía libre de Gibbs como condición de equilibrio.
Por lo tanto : ∑ μ fase d n fase =0
fases

Para el equlibrio solido−liquido : μ sol dn sol + μliq dn liq=0

Para un sistema de un solo componente, el cambio en una de las fases es igual al
cambio en la otra. Es decir, que los potenciales químicos de las dos fases son
iguales. Por lo tanto, en equilibrio, los potenciales químicos de las múltiples fases
del mismo componente son iguales
c) La constante de equilibrio
Para el equilibrio liquido-vapor se define para cada especie i un llamado valor K (o
relación de equilibrio vapor-liquido), que viene dado por:
yi
K i=
xi
Siendo “y” la fracción molar en la fase vapor y “x” la fracción en la fase liquida.
Existen dos valores de K, uno para cada fase liquida en equilibrio con la única fase
de vapor. Por tanto:
I yi
Ki = I
xi
II yi
Ki = II
xi
d) Simplificaciones a la Constante de Equilibrio.
Sistemas Ideales: Se define como solución ideal a aquella que se comporta según
la ley de Raoult para cada uno de los componentes, en todo el rango de
concentración. Los sistemas reales pueden presentar mayor o menor alejamiento
del comportamiento ideal, según el grado de interacciones intermoleculares que
existe en la fase líquida.
s
Pi
k i=
P
Ley de Dalton. - La ley de Dalton o ley de las proporciones múltiples formulada en
1803 por John Dalton, es una de las leyes estequiometrias más básicas. Fue
demostrada por el químico y físico francés Joseph Gay-Lussac. Esta ley afirma
que cuando dos elementos se combinan para originar diferentes compuestos,
dada una cantidad fija de uno de ellos, las diferentes cantidades del otro se
combinan con dicha cantidad fija para dar como producto los compuestos, están
en relación de números enteros sencillos.
Pi= y i P
Ley de Raoult. - La Ley de Raoult establece que la relación entre la presión de
vapor de cada componente en una solución ideal es dependiente de la presión de
vapor de cada componente individual y de la fracción molar de cada componente
en la solución. La ley debe su nombre al químico francés François Marie Raoult
(1830-1901).
S
Pi=x i P i
Su aplicación se cumple sólo para disoluciones ideales no obstante es una buena
aproximación cualitativa, se emplea en teoría de la destilación. Presión menor o
igual a 3 atmosferas y Componentes Polares. - El coeficiente de actividad es
diferente de uno mientras que el factor z es igual a uno.
Sistemas No Ideales.
Coeficientes de Actividad. - El coeficiente de actividad es un número que expresa
el factor de actividad química de una substancia en su concentración molar. La
concentración medida de una sustancia no puede ser un indicador preciso de su
eficacia química como se representa por la ecuación de una reacción particular; en
dados casos, la actividad es calculada al multiplicar la concentración por el
coeficiente de actividad. En los solutos, el coeficiente de actividad es una medida
de cuanto la solución difiere de una disolución ideal. Presión mayor a 3 atmosferas
y Componentes No Polares. - En este caso el cálculo de la constante de equilibrio
se realiza a partir de la división del coeficiente de fugacidad parcial del
componente i en fase liquida entre el coeficiente de fugacidad parcial del
componente i en fase vapor.
ϕ iL
k i=
ϕi V
Presión mayor a 3 atmosferas y Componentes Polares. - Se utiliza
adicionalmente el coeficiente de actividad del componente i del líquido.
γ i , l∗ϕ i ,l
k i=
ϕ i ,l∗V
e) Ecuaciones para el cálculo de las presiones de vapor (Antoine,
Hardlacher, Riedel- Planck-Miller, Thodos-Gomes-Nieto),
coeficientes de fugacidad (Ecuaciones de estados VDW, RK, SRK,
PR) y de actividad (NRTL, WILSON, UNIQUAC, UNIFAC).
f) Ecuaciones para el cálculo de las presiones de vapor
Antoine
B
S B A−
ln P = A− → P S=e T +C
T +C

Ecuación de Riedel-Planck-Miller

−G
∗[ 1−T r +k ( 3+T r )∗( 1−T r ) ]
2 3
ln Pr =
Tr

Donde

G=0,4835+04605 h

h=T ebr∗
( ln Pc
1−T ebr )
 h
−( 1+ T ebr )
G
k= 2
( 3+T ebr ) ¿ ( 1−T ebr )

T eb =temperatura de ebullicion normal reducida a 0 msnm

r

Thodos-Gomes-Nieto
T b∗ln ( Pc )
S=
T c−T b
0.29 0.72
m=0.0052∗M ∗T c
2.464 9.8∗10 −6
∗M ∗T c
γ= e
M
1
−1
T br
a=
1−T 7br
1
m
−1
T br
b= 7
1−T br
γ −a∗s
β=
b

ln ⁡(Pr)=β
[ 1
m
T br ]
−1 +γ [ T br −1 ]
7

P=Pr∗Pc
g) Ecuaciones para el cálculo de Coeficientes de Fugacidad.
Ecuaciones de estado:
Van Der Waals.
Redlich-Kwong.
Soave-Redlich-Kwong.
Peng Robinson.
 h) Coeficiente de actividad (NRTL, WILSON, UNIQUAC, UNIFAC).
Modelo NRTL (Non-Random-Two-Liquid)
G
E
G21 τ 21 G12 τ 12
= +
x1 x 2 RT x 1+ x 2 G21 x 2 + x 1 G12
2
lnγ 1=X 2 ¿
2
lnγ 2=X 1 ¿
Gij =exp (−α τ ij )
bij
τ ij =
RT

Metodo UNIFAC (UNIQUAC Functional-group Activity Coefficients)

ln T k =Qk ¿ ¿
G

X m Qm
∑ V (mj) x j
j−1
θm = G
Xm= G G

∑ XnQn ∑ ∑ V (mj ) x j
n−1 j−1 m −1

T mk =exp ( −aT ) mk

Modelo UNIQUAQ (Universal Quasi-Chemical)

( ) ( )
n
G
E
Ψi z n θi n
=∑ x i ln + ∑ q i x i ln −∑ qi xi ln ¿ ¿ ¿ ¿ ¿ ¿ ¿
RT i=1 x i 2 i=1 Ψ i i =1
xi r i cion
Ψ i= n
segmento
de
∑ xi ri
i=1

xi qi
θi = n
fraccion de area
∑ x i qi
i=1

z=numero de coordinacion de reticulo=10

τ ji =exp ( −μ ji −μii
RT )
Modelo de Wilson
G
E
− x1 x2
= ln ( x 1 +∧12 x2 ) − ln ( x2 +∧21 x 1 )
x1 x 2 RT x 1 x 2 x 1 x2
 ∧12 ∧21
ln γ 1=−ln ( x1 +∧12 x 2 ) + x 2 ( ¿ ¿− )¿¿
x 1 +∧12 x2 x 2 +∧21 x1
∧12 ∧21
ln γ 2=−ln ( x2 +∧21 x 1 ) + x 1 ( ¿ ¿− )¿ ¿
x 1 +∧12 x 2 x 2 +∧21 x1
Puntos de burbuja y de rocío.
La temperatura de burbuja se refiere a la temperatura a la cual la primera burbuja
de tamaño diferencial egresa de una solución líquida con n componentes de
composición conocida a una presión constante.
El punto de rocío o temperatura de rocío es la temperatura a la que empieza a
condensarse el vapor de agua contenido en el aire, produciendo rocío, la
temperatura del punto de roció depende de la presión del gas, y de la
composición.
Medidas de concentración (ºGL, ºBrix, IR).
 GL: Los grados Gay Lussac es la medida de alcohol contenida en volumen,
es decir, la concentración de alcohol contenida en una bebida
 Brix: sirven para determinar el cociente total de materia seca (generalmente
azúcares) disuelta en un líquido. Una solución de 25 °Bx contiene 25 g de
sólido disuelto por 100 g de disolución total.
 IR: El índice de refracción de un medio homogéneo es una medida que
determina la reducción de la velocidad de la luz al propagarse por un medio.
De forma más precisa, el índice de refracción es el cambio de la fase por
unidad de longitud, esto es, el número de onda en el medio (k) será n veces
más grande que el número de onda en el vacío (k 0).

3. DESCRIPCIÓN DE LA PRÁCTICA
a) Materiales y Equipos
 1 litro de alcohol comercial.
 2 litros de agua (vital)
 Cinta Masquin
 Un paquete de clínex
 1 paquete de vasos desechables.
 14 botes de 120ml
 1 paquete de bombillas
 1 estilete
 2 repasadores o servilletas de tela
 2 botes de 2L (plásticos)
 1 marcador de CD permanente
 Caja de cartón para todo el material
 b) Verificar la concentración el alcohol etílico comercial, para ello:
 Vaciar alcohol etílico comercial en una probeta de 500 ml.
 Medir la concentración en ºGL del alcohol etílico comercial utilizando el
alcoholímetro y el densímetro digital, para verificar.
c) Preparar varias soluciones a distintas concentraciones de etanol y agua,
según la siguiente tabla:
VALCOHOL VAGUA
(ml) (ml)
100 0
95 5
90 10
85 15
80 20
70 30
60 40
50 50
40 60
30 70
20 80
10 90
5 95
0 100

d) Con cada mezcla del inciso c), realizar el siguiente procedimiento:

1. Sacar una muestra de aproximadamente 5 ml.
2. Medir la concentración de la muestra en: ºBrix, IR, densidad y ºGL, con los
respectivos equipos: Refractómetro, densímetro digital.
3. Calentar la solución hasta alcanzar el punto de burbuja, cuando empieza a
hervir y se forman las primeras burbujas de vapor.
4. Estar atentos al momento de alcanzar el punto de burbuja, pues si no leen
correctamente la temperatura, deberán volver a preparar la mezcla y repetir
la lectura.
e) Determinar la presión de operación de las pruebas experimentales
utilizando:
 Tablas
 Ecuaciones de presión de vapor
  Recurrir a Senami
f) Determinar el Modelo Matemático a utilizar.
g) El procedimiento de la práctica, a realizar con cada mezcla del inciso b),
es el siguiente:
1) Sacar una muestra de aproximadamente 5 ml, y medir su
concentración en ºBrix, IR, densidad y ºGL con los respectivos equipos
(Refractómetro, densímetro digital).
2) Calentar la solución hasta alcanzar el punto de burbuja, cuando
empieza a hervir y se forman las primeras burbujas de vapor, y leer la
temperatura con diferentes equipos (termómetro de vidrio, termocupla y
termómetro digital)
3) Calcular la fracción molar, el % en volumen.
4) Calcular la temperatura de burbuja y comparar con la lectura
experimental.

4. CÁLCULOS Y GRAFICAS
TABULACIÓN DE DATOS
Tabla 1

No °Brix °GL Tb Tb
Termometro Termometro
de Hg (°C) digital (°C)
1 18,5 97,0 69,0 68,3
2 20,0 83,0 70,0 70,6
3 21,0 89,0 68,0 67,8
4 21,3 84,0 69,0 69,1
5 19,0 80,0 74,0 72,7
6 18,0 71,0 71,0 71,0
7 17,5 61,5 73,0 72,5
8 16,0 51,3 68,0 69,0
9 14,5 41,0 69,0 68,0
10 10,0 31,8 75,0 73,9
11 6,2 19,5 75,1 76,1
12 2,5 9,8 83,0 81,9
13 1,2 0,5 85,3 86,0
14 0,0 0,0 88,6 89,0
 5. CÁLCULOS Y TABULACIÓN DE RESULTADOS
a) Cálculo de la fracción molar.
'
V al =V al∗° GL
' '
V agua =V al∗0 , 04+V agua
ml∗gr
∗mol
1 [ ] ml
n al =V al∗ρ al ¿
P M al gr
ml∗gr
∗mol
1 ml
n agua=V agua∗ρ agua [ ¿ ]
P M agua gr
nT =n al + nagua
nal
x al =
nT
n agua
x agua =
nT

ρagua 0,998 PM H2O 18,01528

gr/ml

ρal 0,786 PM 46,06904

gr/ml al

Donde
Tabla 2

Tabla 3- Tabulación de resultados de numero de moles calculados

V alcohol V agua n alcohol n agua n total x al x agua

97,0 3,0 1,655 0,166 1,821 0,909 0,091

78,9 25,0 1,345 1,387 2,733 0,492 0,508
80,1 27,9 1,367 1,546 2,913 0,469 0,531
71,4 39,3 1,218 2,178 3,396 0,359 0,641
64,0 48,8 1,092 2,703 3,795 0,288 0,712
 49,7 65,2 0,848 3,612 4,460 0,190 0,810
36,9 77,9 0,630 4,313 4,943 0,127 0,873
25,7 87,2 0,438 4,828 5,266 0,083 0,917
16,4 93,4 0,280 5,176 5,456 0,051 0,949
9,5 97,2 0,163 5,382 5,545 0,029 0,971
3,9 99,2 0,067 5,497 5,563 0,012 0,988
1,0 99,9 0,017 5,534 5,551 0,003 0,997
0,3 100,0 0,004 5,539 5,543 0,001 0,999
0,0 100,0 0,000 5,540 5,540 0 1

b) Calcular el % en volumen.

V al
% V al = ∗100
V al +V agua
V agua
% V agua= ∗100
V al +V agua
Tabla 4 - % en volumen.

%V al %V agua

97 3
75,9 24,1
74,2 25,8
64,5 35,5
56,7 43,3
43,3 56,7
32,2 67,8
22,7 77,3
14,9 85,1
8,9 91,1
3,8 96,2
1,0 99,0
0,2 99,8
0 100
 c) Determinar la presión de operación de las pruebas experimentales
utilizando
Tablas
T mercurio T digital
90,5 92,8
Tabla 5

T mercurio P
90 70,14
90,5 71,58
95 84,55
T digital P
90 70,14
92,8 78,21
95 84,55
Ecuación de Antoine
sat B mmhg∗0 , 133322 kpa
ln ⁡(P ¿ ¿ i )=A− [¿ ] [ ¿ ] kpa ¿
T i+ C 1mmhg
Donde
Tabla 6

T digital
A B C T agua
mercurio °C
°C

8.0713 1730.6300 233.4260 90,5 92,8

Cálculo de presión atmosférica con ecuación de Antoine

s
P =ⅇ
( A−
T +C )
B
i
i

Tabulación de resultados de presión calculados para ecuación de Antoine

Tabla 7

Ps Ps digital
mercurio

71.37435341 77.84124914
8 6

Ecuación de Riedel-Planck-Miller

−G
∗[ 1−T r +k ( 3+T r )∗( 1−T r ) ]
2 3
ln Pr =
Tr

Donde

( )
G=0,4835+04605 H ln P c
H=T eb r∗
1−T eb r
 H P
−( 1+ T eb r ) Pr = → P=P r∗Pc
G Pc
k= 2
( 3+T ebr ) ¿ ( 1−T ebr )
T eb normal
T eb r =
Tc

T
T r=
Tc

Tabla 8
Pc Tc T eb normal T eb normal
atm k Mercurio Digital
k k
218,307 647,3 363,65 365,95

Cálculo de presión atmosférica con ecuación de Riedel-Planck-Miller

P =ⅇ
( −G
T
∗[ 1−T +k (3+ T )∗( 1−T ) ]
r
2
r r
) r
3

Tabulación de resultados de presión calculados para ecuación de Riedel-Planck-Miller

Tabla 9

Tr H G K Pr P (atm) P (kpa)
T 0,562 6,905 3,663 0,472 0,00461 1,006 102,0
mercurio
T digital 0,565 7,005 3,709 0,480 0,00454 0,991 100,4

Ecuación Thodos-Gomes-Nieto.

S=T e b
ln ( Pc
1,01325 )
r
1−T e b r

m=0,0052 ( PM 0 ,29∗T 0c , 72)

2,464
exp ( 9 ,8 x 10 ∗PM ¿ T c )
6
γ=
m
 1
−1
T eb r
a= 7
1−T ⅇb r
1
m
−1
T ⅇb r
b= 7
1−T ⅇb r
γ −a ⋅ s
β=
b

ln ( Pr )=β
[ ]
1
m
Tr
−1 + γ [ T r −1 ]
7

Donde
Tabla 10
Pc Tc T eb normal T eb normal
atm k Mercurio Digital
k k
218,307 647,3 363,65 365,95

Tabulación de resultados de presión calculados para ecuación de Thodos-Gomes-Nieto

Tabla 11

Tr S M Y A b B Pr P P(kpa)
(atm)
T 0,562 6,888 1,271 1,367 0,7940 1,100 - 3,7 0,0046 1,01 102,8
mercuri 3
o
T digital 0,565 6,988 1,271 1,367 0,7833 1,084 - 3,79 0,0049 1,07 108,4

Grafica 1Grafica de dispersión

 Dispersion
80
72 72
70 72 72

60

50

40

30

20

10

0
65.00 70.00 75.00 80.00 85.00 90.00 95.00 100.00 105.00
Series2
 12

10

0
0 2 4 6 8 10 12

Grafica 2 Grafica de comparación de errores

comparacion de error

102.00 102.80

110
71.58 71.37
90 72 72 72 72

70

50

30
Series2
10 Series1
1 2 3 4
Series1 72 72 72 72
Series2 71.58 71.37 102.00 102.80

Series1 Series2
ALGORITMO
Datos
P=0,721 bar
X alcohol = conocido
1. Suponer Temperatura de burbuja
2. Calcular PSi (Peng-Robinson)
3. Calcular coeficiente de actividad γ i (UNIQUAC )
sat
γ i∗P i
4. Calcular la constante de equilibrio k i=
P
5. Calcular y i=k i∗x i
6. ¿ ∑ yi =1?
7. NO → volver a paso 1
8. Si → FIN
Leer T, y i
Tabla 12 Paso 1 temperatura supuesta

T sup T sup K

68,77454373 341,9245437
68,56622055 341,7162206
68,76441548 341,9144155
69,1092207 342,2592207
69,49192996 342,64193
70,29063125 343,4406312
71,13751049 344,2875105
72,12048043 345,2704804
73,4506966 346,6006966
75,32919751 348,4791975
78,23654647 351,3865465
82,9821622 356,1321622
86,39861383 359,5486138
Tabla 13 Paso 2 cálculo de presión parcial de cada componente
Ps alcohol Ps gua Kpa
Kpa

Tabla 13 Paso 3 cálculo de coeficiente de actividad

Λ21 ln ⁡(γ ¿¿ alco h ol)¿ ln ( γ agua ) γ alco h ol γ agua

Λ12
Tabla 14 Paso 4 Calculo de la constante de equilibrio
Kalcohol kagua

Tabla 15 Calcular y i=k i∗x i

Kalcohol kagua
Tabla 16 Calculo de Yi

Y alcohol y agua

Tabla 17 ∑ y i

∑ yi
 ANALISIS DE RESULTADOS
En la comparación de los resultados de presión teórica con la presión atmosférica
obtenido de senamhi oscilan de 30%. Estos errores pueden ser debido a una mala
manipulación de los instrumentos o que los instrumentos ya no están en condiciones para
dar una lectura fiable.
Para el cálculo de temperatura del punto burbuja si bien se asignó la ecuación de estado
de pen-Robinson el cálculo se realizó con la ecuación
CONCLUSIONES Y RECOMENDACIONES
1.1. Conclusiones:
 Se determino experimental y analíticamente la temperatura de burbuja del sistema
etanol - agua a distintas composiciones.
 Se logró construir todos los diagramas como estaba previsto en los objetivos
T burbuja, experimental – x Etanol.
T burbuja, teórico – x Etanol.
T rocío – x Etanol.
T burbuja, experimental – x – y para etanol.
T burbuja, experimental – x – y para etanol.
x – y para etanol.
ºBrix – x Etanol.
ºGL – x Etanol.
ºGL – ºBrix.
IR – X Etanol.
Combinar los dos diagramas 1) y 2).
Combinar los dos diagramas 1) y 3).
Combinar los dos diagramas 2) y 3).

1.2. Recomendaciones:
 Se debe prestar mucha atención en la selección de ecuación de estado para
no tener muchos errores, yo que no todas las ecuaciones son para cualquier
sustancia.
 A la hora de medir el volumen del alcohol, se tiene que hacer lo más rápido
posible, ya que el alcohol es más volátil que el agua.

BIBLIOGRAFIA
Bibliografía
Henao, J. C. (03 de 06 de 2022). youtube. Obtenido de [Link]
john m. prausnitz, r. n. (s.f.). molecular thermodynamics of fluid-phase equilibria. new jersey: joan l
mcNamara.

Walas, S. M. (1985). Phase Equilibria in Chemical Engineering. Butterworth Heinemann.

[Link]