Calculo Del Equilibrio Gas-Liquido

PROPIEDADES DE LOS FLUIDOS
INSTITUTO TECNOLOGICO SUPERIOR

DE VILLA LA VENTA
Ing. Cristino Santiago Ruiz Huerta
w w w . p e p .p e m e x .c o m
Agosto-2014

UNIDAD 4
TEMA: CALCULO DEL EQUILIBRIO GAS-LIQUIDO
CON ECUACIONES DE ESTADO
4.1.SUBTEMAS
Sustancias Puras
4.1.1.- Calculo del potencial qumico

4.1.2.- Calculo de la fugacidad
4.1.3.- Calculo del coeficiente de fugacidad
4.2
Mezclas
4.2.1.- Calculo del potencial qumico de mezcla

4.2.2.- Calculo del coeficiente de fugacidad
en una mezcla

UNIDAD 3
TEMA: COMPORTAMIENTO DE FASES
4.3.- Ecuacin de Van Der Waals
4.4.- Ecuaciones de estado en el punto critico
4.5.- Ecuacin de estado cubica de Van Der Waals
SUBTEMAS
de dos parmetros
4.6.- Ecuacin de estado de Redlich - Kwong
4.7.- Reglas de mezclado para la ecuacin de estado
de Redlich Kwong.
4.8.- Ecuacin de estado de Soave- Redlich Kwong.

UNIDAD 3
TEMA: COMPORTAMIENTO DE FASES
4.9.- Ecuacin de estado de Peng-Robinson.
4.10.- Reglas de mezclado para las ecuaciones
de estado de Soave- Redlich Kwong y
SUBTEMAS
Peng-Robinson

Sustancias Puras
Las Ecuaciones de Estado de acuerdo a van der Waals produce Isotermas
Como se muestra en la siguiente fig. para una sustancia pura .
La isoterma calculada de la fig. para temperaturas en de la temperatura
critica mostradas por van der Waals para ciertas temperaturas las ---Presiones calc. en el trayectoria son negativas como se muestran. Esta
Trayectoria no aparece experimentalmente. Para una sust. pura una lnea

Horizontal une o conecta el EQUILIBRIO gas liquido.
En LOOP de van der Waals es usado para determinar el volumen molar
cuando el gas y el liquido estn en equilibrio y entonces es reemplazado
por una la lnea de unin entre esos dos volmenes. Estas son las bases
para calcular el potencial Qumico para gases iguales en equilibrio y el
potencial qumico del Liquido en equilibrio.

Calculo del potencial qumico
El potencial qumico G para un fluido puro a temperatura constante ----Puede ser calculado como:
dG = VM dp
El potencial qumico tambin es llamado Energa Molar libre de GIBBS
Para las isotermas calculadas debajo de la temperatura critica.
El punto a es seleccionado arbitrariamente a lo largo de la parte liquida
De la isoterma. El punto d es seleccionado arbitrariamente a lo largo de
La parte de gas de la isoterma. Los puntos b y c son los puntos maxiMos y minimos del nudo de van der Waals. El punto e representa los pun
Tos a los largo del nudo para cada potencial quimico de la ecuacin. El
Punto f es el punto a lo largo de la lnea bc que tienen la misma presin
Que
w welw .punto
p e p .p e me e x .c o m

Calculo de la fugacidad
Para un gas ideal puro:
dG = (RT/p) dp = RT d (ln p)
El potencial qumico de un fluido real puede expresarse reemplazando
La presin p de la ecuacin anterior con una propiedad llamada fugaci
Dad f .
dG = RT d ( ln f )
Esta es la ecuacin que define la fugacidad en forma diferencial en refeRencia a valores que son requeridos.
Para una sustancia pura :
d ( ln f ) = ( VM/RT ) dp
La fugacidad tiene unidades de presin.
Calculo de la fugacidad
Para una sustancia pura la relacin entre la fugacidad y la presin
f/p se llama Coeficiente deFugacidad.
--
ln f/p = z 1 ln z + 1/RT (RT/VM p ) d VM

Ecuacin de Van Der Waals
Durante el desarrollo de la ecuacin de estado para gases ideales, por
ejemplo la ecuacin 2.23, se consideraron dos suposiciones, la primer
suposicin considera que el volumen de las molculas de gas es
insignificante en comparacin con el volumen del contenedor
(recipiente que lo contiene) y la distancia entre las molculas,
y la segunda suposicin considera la inexistencia de fuerzas de
atraccin o de repulsin entre las molculas del gas o entre el gas
y las paredes del recipiente que las contiene.
pVM RT
----------------- 4.3

Una de las primeras ecuaciones de estado para representar el comporta
miento de los gases reales fue desarrollada por van der Waals en 1873.
van der Waals se enfoc en tratar de eliminar las dos suposiciones mencio
nadas anteriormente mediante su ecuacin de estado, explicando que las
molculas de gas ocupan una fraccin significante de volumen a
presiones altas.
Para contrarrestar la primera suposicin van der Waals propuso que el
volumen de las molculas, representado por el parmetro b debe de
restarse del volumen molar real VM en la ecuacin 2.23; arreglando
esta ltima ecuacin se tiene,
pVM RT
RT
p
VM

Aplicando el parmetro b
RT
p
VM b
------------- 4.4
en donde el parmetro b representa el volumen de las molculas que

se restan y VM es la variable que representa el volumen molar real,
ambos en ft3 por una unidad de lbm-mol.
Para eliminar la segunda suposicin, van der Waals rest un trmino
de correccin, denominado a/V2M a la ecuacin 4.4 para tomar en cuenta
las fuerzas de atraccin entre las molculas. van der Waals propuso
la expresin matemtica siguiente:
RT
a
p
2
VM b VM
---4.5
O bien
a
p 2
VM
--- 4.6
VM b RT

en donde la ecuacin 4.6 representa la ecuacin de van der Waals,
la cul es diferente de la ecuacin de los gases ideales (ec. 2.23)
por la adicin del trmino a/V2M a la presin y la sustraccin del
parmetro b del volumen molar.
El trmino a/V2M representa un intento para corregir la presin debido
a las fuerzas de atraccin entre las molculas. Es decir, la presin real
ejercida sobre las paredes del recipiente que contiene al gas real es
menor por la cantidad a/V2M que la presin ejercida en el mismo
recipiente por un gas ideal.

El parmetro b representa un intento por corregir el volumen
molar
debido al volumen ocupado por las molculas. Las constantes
a y b dependen del tipo de gas. De las ecuaciones 4.4 a 4.6 R
representa la constante universal de los gases e igual a 10.73
(lb/pg2abs-ft3)/(lbm-mol -oR), p es la presin del sistema en
lb/pg2abs, T es la temperatura del sistema en oR y V es el
volumen
molar en ft3/ mol.
Los parmetros a y b representan constantes y caracterizan las
propiedades moleculares de cada componente de la mezcla.

A la ecuacin de van der Waals (ecuacin 4.5) se le denomina ecuacin
de estado de dos constantes (aunque es realidad contiene tres
constantes: a, b y R) ecuacin de estado cbica. Cualesquiera ecuacin
de estado se puede representar en una forma general como
RT
a
p
2
VM b VM
p p repulsin p atraccin
En donde
p repulsin
RT
(VM b)
p atraccin
2
VM
----------- 4.5

Ecuaciones de estado en el punto critico
van der Waals al tratar de determinar experimentalmente los valores
de las constantes a y b para cualquier sustancia pura, en un diagrama
de fase de presin contra volumen molar observ que la isoterma
crtica, representa una pendiente horizontal y un punto de inflexin
en el punto crtico (Fig. 4.1). Esta observacin para una sustancia pura
se representa matemticamente como:
dp
dVM
0
TC
------ 4.10
d2p
2
dVM
0
TC
------- 4.11
dp
dVM
P r e s i pn ,( lb /p2agb s )
pc
0,
Tc
d2p
2
dVM
0
Tc
Tc
VMc
Volumen molar, VM (ft3)
Fig. 4.1 - Diagrama de presin contra volumen para una sustancia

pura
w w mostrando
w . p e p .p e m e x las
.c o mcondiciones en el punto crtico

Obteniendo la primera y segunda derivada de la ecuacin de van der
Waals (ecuacin 4.5) con respecto al volumen al punto crtico, e
igualando ambas ecuaciones a cero se obtiene respectivamente
dp
dVM
d 2P
2
dVM
TC
TC
RTc
MC
2a
3 0
2
VM C
b
2 RTc
VM C b
---------- 4.12
6a
4 0 ------------- 4.13
VM C

Resolviendo las ecuaciones 4.12 y 4.13 simultneamente para los
parmetros a y b, se obtiene
9
a VM C Tc R
8
1
b
VM C
3
--- 4.14
--4.15
Ahora bien, expresando las ecuacin de estado de van der Waals

(ecuacin 4.6) en el punto crtico, se tiene:
a
pc
VM b RT c
2
C
V
w w w . p e p .p e m e xM.c oCm
----------- 4.16
combinando las ecuaciones 4.12, 4.13 y 4.16 se obtiene,
2
c
27 R T
a
64 pc
---- 4.17
RT c
b
8 pc
----- 4.18
Ejemplo 4.1Constantes de la ecuacin de estado de van der Waals.

Calcular las constantes de van der Waals para el 3-metil-hexano.
Solucin.
Para el 3-metil-hexano la Tc y pc se obtienen de la Tabla A-1 del Apndice
A como 963.8 oR y 408.1 lb/pg2abs, respectivamente. Sustituyendo
valores en la ecuacin 4.17,
a
( 27 )( 10 .732 lb/pg 2 abs ft 3 / lbm mol R )2 ( 963.8 oR )2

( 64 )( 408 .1 lb/pg 2 abs )
( lb / pg 2 abs ) ( ft 3 )2
110 ,599.29
2
(
lbm
mol
)
(10.732 lb/pg 2 abs ft 3 /lbm mol R )( 963.8 R)
w w w . p e p .p e m e (8)(408.1
x .c o m
lb/pg 2 abs )
ft 3
3.168
lbm
mol
Ecuaciones de Estado para Gases Reales

EOS de Van der Waals
EOS Redlich-Kwong
EOS Soave-Redlich-Kwong
EOS Peng-Robinson
Reglas de Mezclado
Ecuaciones de Estado
Se defini una ecuacin de estado como una relacin matemtica entre

las variables Temperatura, Presin y Volumen (ms la concentracin
en el caso de mezclas). Las ecuaciones de estado constituyen una
herramienta usual en Termodinmica para el clculo de propiedades
fsicas y termodinmicas de fluidos. Si se conoce la relacin
f(T,P,V)=0 para fluidos puros o bien g(T,P;V,X)=0 para mezclas,
entonces es posible determinar el valor de propiedades tales como
densidad, entalpa, entropa o energa interna.
Ecuacin de Van der Waals
La ecuacin de Van der Waals corrige

drsticas del modelo de Gas Ideal:
dos
de
las
aproximaciones
mas
1.- El volumen de las molculas no es despreciable con respecto al volumen

del sistema
2.- Las molculas se atraen y producen una presin de atraccin, que van
der Waals model como a / V~2
2
M
M
RT
a
p
V b V
Los valores de a y b son caractersticos de un gas en particular:
TC
d 2P
2
dVM
TC
RTc
MC
2a
0
3
VM C
2 RTc
6a
0
3
4
VM C
VM C b
27 R 2Tc2
a
64 pc
RTc
8 pc
dp
dVM
pc
Presin, p (lb/pg2abs)
dp
dVM
0,
Tc
d2 p
2
dVM
0
Tc
Tc
VMc
Volumen molar, VM (ft3)
a EOS fue planteada considerando substancias puras; se deben usar reglas de mezclad
a el caso de gases reales.
Calcular las constantes de van der Waals para el 3-metil-hexano.

Solucin.
Para el 3-metil-hexano la Tc y pc se obtienen de la Tabla A-1
963.8 oR y 408.1 lb/pg2abs, respectivamente.
a
( 27 )( 10 .732 lb/pg 2 abs ft 3 / lbm mol R )2 ( 963.8 oR )2

2
( 64 )( 408 .1 lb/pg abs )

b
( lb / pg 2 abs ) ( ft 3 )2
110 ,599.29
( lbm mol )2
(10.732 lb/pg 2 abs ft 3 /lbm mol R )( 963.8 R)

(8)(408.1 lb/pg 2 abs )
ft 3
lbm mol
3.168
ecuacin de van der Waals puede ser escrita en forma polinomial:
ab
RT 2 a
VM VM
0
VM3 b
p
p
p
VM
zRT
p
z 3 ( 1 B )z 2 Az AB 0
A
ap
R 2T 2
bp
RT
Clculo de las densidades de fluidos con la ecuacin de estado van der Waals.
En un cilindro cerrado se encuentra propano puro a 100 F teniendo las fases
lquido y vapor presentes. Empleando la ecuacin cbica de estado de van der
Waals de dos parmetros, calcular la densidad de las fases lquido y vapor
(gas).
Temperatura (F)
Temperatura (F)
Solucin.
Etapa 1. Estimando
la pv del propano
puro a partir de
la cartas de Cox
Presin (lb/pg abs.)
Presin (lb/pg2 abs.)
pv = 185 lb/pg2abs
Temperatura (F)
Temperatura (F)
Etapa 2. Calculando los parmetros a y b
27 R 2Tc2
a
64 pc
a (0.421875)
Tc
Pc
(10.732 lb / pg 2 abs ft 3 / lbm mol R ) 2 (666 R ) 2

(616.3 lb / pg 2 abs )
lb / pg 2 abs ( ft 3 ) 2
34,970.465
2
(lbm mol )
RTc
8 Pc
b (0.125)
(10.732 lb / pg 2 abs ft 3 / lbm mol R )( 666 R )
(616.3 lb / pg abs )
1.4496
ft 3
lbm mol
Etapa 3. Calcular los coeficientes A y B aplicando las ecuaciones 4.25 y 4.26
ap
R 2T 2
34 ,970.465 lb / pg
A
abs ft 6 /( lbm mol)2 ( 185 lb / pg 2 abs )

0.17912
2
3
o
2
o
2
( 10.732 lb / pg abs ft / lbm mol R ) ( 560 R )
bp
RT
(1.4496 ft 3 / lbm mol )(185 lb / pg 2 abs )

B
0.04462
2
3
o
o
(10.732 lb / pg abs ft / lbm mol R )( 560 R)
Etapa 4. Sustituyendo los valores de A y B dentro de la ecuacin:
z ( 1 B )z Az AB 0
3
2
z
1
.
04462
z
0polinomio
.17912de
z tercer
0.007993
0
Etapa 5. Resolviendo el
grado,obteniendo
3
el mayor
y menor valores (races) del polinomio. Para resolver este polinomio
se utilizan mtodos directos o iterativos. Para la fase vapor se
obtiene,
z v = 0.84350
y para la fase lquida,
z L = 0.07534
Etapa 6. Calculando las densidades de las fases lquido y vapor. Para la fase vapor se tiene,
pM C3
z V RT
lbm
( 185 lb / pg 2 abs )( 44 lbm / lbm mol )

1.60572 3
2
3
o
o
( 0.84350 )( 10.732 lb / pg abs ft / lbm mol R )( 560 R )
ft
v
pM C3
z L RT
lbm

17.97753 3
2
3
o
o
( 0.07534 )( 10.732 lb / pg abs ft / lbm mol R )( 560 R )
ft
L
Ecuacin de Redlich-Kwong(RK)
En 1948 Redlich
y
Kwong modificaron el trmino de fuerzas de atraccin
(fuerzas de presin a/V2 ) de la ecuacin de van der Waals, lo cul mejora en
forma considerable la prediccin de las propiedades fsicas y volumtricas de
la fase gas. Redlich-Kwong sustituyeron el trmino de fuerzas de atraccin de
presin, con un trmino general de dependencia de la temperatura. La ecuacin
de Redlich-Kwong se expresa como,
RT
a
p
V b V V b T 0.5
R 2Tc2.5
a a
pc
a 0.42747
RT c
b b
pc
b 0.08664
La EOS de RK expresada en forma polinomial:
z 3 z 2 ( A B B 2 ) z AB 0
ap
A 2 2.5
R T
bp
B
RT
En un cilindro cerrado se encuentra propano puro a 100 F teniendo las fases

lquido y vapor presentes. Empleando la ecuacin cbica de estado de van der
Waals de dos parmetros, calcular la densidad de las fases lquido y vapor
(gas).
Etapa 1. Calculando los parmetros a, b, A y B.
R 2Tc2.5
a a
pc
b b
RTc
pc
ap
R 2 T 2. 5
bp
RT
a (0.42747)
(10.732 lb / pg 2 abs ft 3 / lbm mol R ) 2 ( 666 R ) 2.5

(616.3 lb / pg 2 abs )
( lb / pg 2 abs ) ( ft 3 ) 2 R 0.5
a 914 ,450.86
( lbm mol ) 2
b (0.08664)
(10.732 lb / pg 2 abs ft 3 / lbm mol R )( 666 R )

2
(616.3 lb / pg abs )
lbm
3
ft
1.0048
( lb / pg 2 abs ) ( ft 3 ) 2 R 0.5
914 ,450.863
( 185 lb / pg 2 abs )
2
( lbm mol )
A
0.19793
2
3
o
( 10.732 lb / pg abs ft / lbm mol R ) 2 ( 560 o R ) 2.5
1.0048
lbm
2
(
185
lb
/
pg
abs )
3
ft
0.03093
2
3
o
o
( 10.732 lb / pg abs ft / lbm mol R )( 560 R )
, ..........................................................................................................................................(4.42)
sustituyendo valores,
Etapa 2. Sustituyendo los valores calculados de A y B en:
z 3 z 2 ( A B B 2 )z AB 0
z z 0.16604 z 0.00612 0
3
Etapa 3. Resolviendo el polinomio de tercer grado de la ecuacin 4.40, obteniendo los

valores mayor y menor se tiene,
zv = 0.80263 y
zL = 0.05271
Etapa 4. Calculando las densidades de las fases:

,
Para la fase vapor se tiene,
pM C3
lbm
1
.
68749
2
3
o
o
3
z RT ( 0.80263 )( 10.732 lb / pg abs ft / lbm mol R )( 560 R )
ft
v
y para la fase lquida,
pM C 3
185 lb / pg 2 abs 44 lbm / lbm mol
L
25.69583
z RT ( 0.05271 )( 10.732 lb / pg 2 abs ft 3 / lbm mol o R )( 560 o R )
L
lbm
3
ft
Ecuacin de Soave-Redlich-Kwong(RK)
En 1972, G. Soave modific el trmino relacionado a la atraccin molecular de

la EDE de RK, quedando,
a
RT
p
V b V V b
1 m 1 Tr1 / 2
m 0.480 1.574 0.176 2
(log p vr 1 )
Factor acntrico de Pitzer

Medida de la acentricidad
(no-esfericidad) de las molculas
un parmetro importante porque la fuerza intermolecular en molculas complejas es
suma de las interacciones entre las diferentes partes de las molculas no solo en
tros
de estados correspondientes: todos los gases que presenten el mismo factor acentri
portan de igual manera en condiciones reducidas.

El factor acntrico de un componente qumico puro se define con referencia
a su presin de saturacin reducida:
f prsat
dp
dVM
0
Tc
d2p
2
dVM
0
Tc
R 2Tc2
a 0.427480
pc
b 0.086640
La EOS de SRK expresada en forma polinomial:
z 3 z 2 ( A B B 2 ) z AB 0
A
ap
( RT ) 2
bp
B
RT
RTc
pc
Clculo de las densidades de las dos fases de una sustancia pura con la
ecuacin de Soave-Redlich-Kwong, SRK .
Empleando la ecuacin de SoaveRedlich-Kwong, SRK , calcular la densidad de las fases vapor y lquido para
propano puro a 100 F que se encuentra en un cilindro cerrado teniendo las
fases1.lquido
y vapor
Etapa
Calculando
la pc, Tpresentes.
y a partir de la Tabla A-1.
c
pc=616 lb/pg2abs, Tc=666.06 R y =0.1522.

Etapa 2. Calculando la Tr.
Tr
T
560 R
0.8407
Tc 666.06 R
Etapa 3. Calculando el parmetro m,
m 0.480 1.574 0.176 2
m 0.480 1.574( 0.1522 ) 0.176( 0.1522 ) 2 0.71549
Etapa 4. Calculando el parmetro ,
1 m1 T
1/ 2
r
1 0.71549 1 0.8407 1 / 2 1.12245

2
Etapa 5. Calculo de los coeficientes a y b,

2
a a
2
c
RT
0.42748
pc
lb / pg 2 abs ft 3
2
10.732
666
.
06
lb / pg 2 abs ( ft 3 ) 2
lbm mol o R
35,441.465
2
2
616.3lb / pg abs
lbm mol
lb / pg 2 abs ( ft 3 ) 2
lb / pg 2 abs ( ft 3 ) 2
1.12245 39781.27239
a a 35441.465
2
2
lbm mol
lbm mol
RTc
b b
0.08664
pc
lb / pg 2 abs ft 3
666.06 R
10.732
lbm mol o R
ft 3
1.00471
lbm - mol
616.3lb / pg 2 abs
Etapa 6. Calculo de los coeficientes A y B,
lb / pg 2 abs ( ft 3 ) 2
2
39781.27239
185
lb
/
pg
abs
2
a p
lbm mol
0.20376
2
2
2
3
( RT )
lb / pg abs ft
10.732
o
560.0 R
lbm
mol
ft 3
185lb / pg 2 abs
1.00471
lbm - mol
bp
0.03331
2
3
RT
lb / pg abs ft
10.732
560.0 R
o
lbm
mol
Etapa 7. Resolviendo la ecuacin 4.40 para zv y zL,
z 3 z 2 ( A B B 2 )z AB 0
z 3 z 2 ( 0.20376 0.03331 0.033312 )z ( 0.20376 x.03331 ) 0
z 3 z 2 0.16934 z 0.00679
zv = 0.79859 y zL = 0.06018
Etapa 8. Calculando las densidades de las fases vapor y lquido, respectivamente:
pM C3
lbm
1
.
69602
2
3
o
o
3
ft
v
pM C3
lbm
22
.
506
2
3
o
o
3
ft
L
Ecuacin de Peng-Robinson(PR)
En 1975, Peng y Robinson propusieron una ecuacin de dos constantes, que

origin grandes expectativas para el mejoramiento de las predicciones de
densidades para el lquido, la cual se expresa matemticamente como:
a
RT
p
VM b VM (VM b) b(VM b)
R 2Tc2
a a
pc
b b
RTc
pc
a 0.457234
b 0.077796
1 / 2 1 m( 1 Tr0.5 )
m 0 .37464 1.54226 -0.26992 2
La EOS de PR expresada en forma polinomial:
z 3 ( B 1 )z 2 ( A 2 B 3 B 2 )z ( AB B 2 B 3 ) 0
A
ap
( RT ) 2
bp
RT
Ecuacin
Van der Waals (1873)
RT
p
v b
Redlich-Kwong (1949)
a
vvv2a
a
0.5
v(b
b
T
b
ab)
(v)
Fuerzas de atraccin
RT
p
v b
Parmetros
27 R 2Tc2
a
64 pc
RTc
8 pc
RTc
R 2Tc2.5
b
a a
b
pc
pc
RTc
R 2Tc2
b b
a a
pc
pc
Soave-Redlich-Kwong
(1972)
p
RT
v b
1 m 1 Tr1 / 2
a 0.42747
b 0.08664
a 0.42747
b 0.08664
(log p vr 1 )
m 0.480 1.574 0.176 2

R 2Tc2 b RTc b 0.077796
a a
b
pc
pc
0.457234
Peng-Robinson (1976)
RT
v b
1 / 2 1 m( 1 Tr0.5 )
m 0 .37464 1.54226 -0 .26992 2
ecuaciones de estado cbicas se pueden representar de forma general,
hmidt y Wenzel 1980):
a
RT
p
v b v 2 k1 bv k 2b 2
Ecuacin
k1
k2
VdW
0.421875
0.12500
R-K
0.427470
0.08664
1/Tc0.5
S-R-K
0.427470
0.08664
P-R
-1
0.457240
0.07780
1 m1 Tr1 / 2
1 m1 Tr1 / 2
m 0.480 1.574 0.176 2
m 0 .37464 1.54226 -0 .26992 2
Ecuaciones de estado para mezclas de hidrocarburos

Reglas de mezclado para determinar los parmetros a m y b m
Ecuacin
am
S-R-K
P-R
xi x j ai a j
VdW
R-K
am
i 1
am
i 1
bm
n
0.5
bm xi bi
j 1
i 1
x x a a 1 k
n
j 1
0.5
ij
bm xi bi
i 1
Ejemplo:
Para una mezcla de tres componentes, n=3, a y b:
Clculo de la densidad de lquido de una mezcla de hidrocarburos con la

ecuacin de Redlich-Kwong. Calcular la densidad de un hidrocarburo lquido
cuya composicin se proporciona en la Tabla y que se encuentra a 4,000 lb/pg2
abs y 160 F. Emplear la ecuacin de estado de Redlich-Kwong.
Composicin de los componentes lquido de la mezcla de hidrocarburos
Componente
xj
(fraccin
mol)
Metano, C1H4
0.45
Etano, C2 H6
0.05
Propano, C3 H8
0.05
Butano normal, n-C4 H10
0.03
Pentano, C5 H12
0.01
Hexano, C6 H14
0.01
Heptano,
C7 H16+ del C =215 lbm/lbm-mol
0.40
Peso
molecular
7+
pc del C7+ = 285 lb/pg2 abs.

Tc del C7+ = 825 R.
cin.
a 1. Calculando los parmetros a j
y b j para cada componente
R 2Tc2.5
a a
pc
b b
RTc
pc
Etapa 2. Calculando los parmetros de mezclado am y bm , respectivamente.

n
am
i 1
x x a a
n
j 1
0 .5
1,136,233.11
bm x j b j 1.45263
j 1
Etapa 3. Calculando los coeficientes A y B , respectivamente,
( lb / pg 2 abs ) ( ft 3 ) 2 R 0.5
1,136 ,233.11
( 4000 lb / pg 2 abs )
2
( lbm mol )
a p
A 2 m 2 .5
4.12275
2
3
o
R T
( 10.732 lb / pg abs ft / lbm mol R )2 ( 620 o R ) 2.5
lbm
(4000 lb / pg 2 abs )
3
ft
b p
B m
0.87326
2
3
o
o
RT (10.732 lb / pg abs ft / lbm mol R)( 620 R)
1.452633
Etapa 4. Resolviendo la ecuacin para encontrar la raz con el menor valor,
z 3 z 2 ( A B B 2 ) z AB 0
z 3 z 2 ( 4.12275 0.87326 0.87326 2 )z ( 4.12275 x0.87326 ) 0
z 3 z 2 2.48676 z 3.6001
Etapa 5. Calculando el peso molecular aparente de la mezcla lquida,
Ma = 100.2547 lbm/lbm-mol
zL = 1.27134
Etapa 6. Calculando la densidad de la fase lquida, considerando el peso molecular

aparente de la mezcla,
pM a
z L RT
lbm
pM a
4 ,000 lb / pg 2 abs 100.2547 lbm / lbm mol
L
47.40571
2
3
o
o
3
ft
L
Clculo de la densidad del gas de una mezcla de hidrocarburos con la ecuacin de

Redlich-Kwong. Empleando la ecuacin de estado de Redlich-Kwong calcular la
densidad de una mezcla de gases a 4,000 lb/pg2abs y 160 F . La composicin de
la mezcla de gases se proporciona en la Tabla :
Composicin de los componentes gaseosas de la mezcla de
hidrocarburos
Componente
yj
(fraccin mol)
Metano, C 1 H 4
Etano, C 2 H 6
Propano, C 3 H 8
Butano normal, n-C 4 H 10
Pentano, C 5 H 12
Hexano, C 6 H 14
Heptano, C 7 H 16+
0.86
0.05
0.05
0.02
0.01
0.005
0.005
Peso molecular del C7+=215 lbm/lbm-mol

pc del C7+ = 285 lb/pg2 abs.
Tc del C7+ = 825 R.
Solucin.
Etapa 1. Obteniendo el peso molecular, la presin y temperatura crtica y los parmetros a y
para cada componente de la mezcla.
R 2Tc2.5
a a
pc
b b
RTc
pc
Etapa 2. Calculando los parmetros de mezclado am y bm, respectivamente,

n
am
i 1
y y a a
n
j 1
0 .5
234,506
bm y j b j 0.563204
j 1
Etapa 3. Calculando los coeficientes A y B, respectivamente:
( lb / pg 2 abs ) ( ft 3 ) 2 R 0.5
234 ,506.00
( lbm mol ) 2
am p
A 2 2.5
R T
( 10.732 lb / pg 2 abs ft 3 / lbm mol
0.563204
( 4000 lb / pg 2 abs )
R ) ( 620 R )
lbm
2
(
4000
lb
/
pg
abs )
3
ft
2.5
0.85089
bm p
0.33857
2
3
o
o
RT ( 10.732 lb / pg abs ft / lbm mol R )( 620 R )
Etapa 4. Resolviendo la ecuacin para encontrar la raz con el mayor valor,
z 3 z 2 ( A B B 2 ) z AB 0
z 3 z 2 (0.85089 0.33857 0.33857 2 ) z (0.85089 x0.33857) 0
z 3 z 2 0.39769 z 0.28809
zv = 0.91068
Etapa 5. Calculando la densidad de la mezcla de gases,

considerando el peso molecular de la mezcla,
pM a
v
z RT
lbm
pM a
( 4 ,000 lb / pg 2 abs )( 20.89 lbm / lbm mol )
13
.
789
2
3
o
o
3
ft
v
Clculo de las densidades de las dos fases de una mezcla de hidrocarburos con la
ecuacin de Soave-Redlich-Kwong, SRK . Una mezcla de hidrocarburos cuya
composicin se muestra en la tabla inferior existe en equilibrio en dos fases a
4,000 lb/pg2abs y 160 F . Calcular las densidades de las fases vapor y lquido
empleando la ecuacin de Soave-Redlich-Kwong, SRK . Considerar que no existe
interaccin binaria entre los componentes de la mezcla.
Peso molecular del C7+ = 215 lbm/lbm-mol
pc del C7+ = 285 lb/pg2abs

Tc del C7+ = 700 R
= 0.52
Solucin.
Etapa 1. Calculando los parmetros , a y b
Etapa 2.
Calculando los parmetros de mezclado: para la fase lquida
a m
=104,373.9 y b=1.8901.
Para la fase vapor se calcula

a m
=9,244.11 y b=0.56897.
Etapa 3. Calculando los coeficientes A y B para cada fase. Para la fase lquida,
a mp 104,373.9 lb / pg 2 abs ft 6 / lbm mol 4,000lb / pg 2 abs

A
2
2
2
( RT )
10,732lb / pg 2 abs ft 3 / lbm mol R 620 R
9.4298
( 1.8901 ft 3 / lbm mol )( 4000 lb / pg 2 abs )

B
1.136
2
3
o
o
Para la fase vapor,
2
a mp
9244.11lb / pg 2 abs ft 6 / lbm mol 4,000lb / pg 2 abs
A
2
2
( RT ) 2
10,732lb / pg 2 abs ft 3 / lbm mol R 620 R
0.8352
( 0.56897 ft 3 / lbm mol )( 4000 lb / pg 2 abs )

B
0.3420
2
3
o
o
Etapa 4. Calculando el valor de zL,
z 3 z 2 ( 9.4298 1.136 1.136 2 )z ( 9.4298 x1.136 ) 0
z 3 z 2 7.0033 z 10.71225
zL =1.41221
Etapa 5. Calculando el valor de zv,
z 3 z 2 ( 0.8352 0.3420 0.3420 2 )z ( 0.8352 x0.3420 ) 0
z 3 z 2 0.37624 z 0.28564
zv =0.92663
Etapa 6. Calculando el peso molecular aparente del lquido y del vapor a partir de la
composicin de la mezcla. Para la fase lquida, se tiene:
Ma= 100.25 lbm/lbm-mol
Para la fase vapor,
Ma = 20.89 lbm/lbm-mol
Etapa 7. Calculando la densidad de cada fase, para la fase lquida y para la fase vapor:
lbm
pM a
42
.
67492
2
3
o
o
3
ft
L
lbm
pM a
13
.
55251
2
3
o
o
3
ft
v

Calculo Del Equilibrio Gas-Liquido

Cargado por

Copyright:

Formatos disponibles

Calculo Del Equilibrio Gas-Liquido

Cargado por

Información del documento

Descripción original:

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Calculo Del Equilibrio Gas-Liquido

Cargado por

Copyright:

Formatos disponibles

PROPIEDADES DE LOS FLUIDOS

INSTITUTO TECNOLOGICO SUPERIOR

PROPIEDADES DE LOS FLUIDOS

4.1.1.- Calculo del potencial qumico

4.2.1.- Calculo del potencial qumico de mezcla

PROPIEDADES DE LOS FLUIDOS

PROPIEDADES DE LOS FLUIDOS

PROPIEDADES DE LOS FLUIDOS

Trayectoria no aparece experimentalmente. Para una sust. pura una lnea

PROPIEDADES DE LOS FLUIDOS

PROPIEDADES DE LOS FLUIDOS

PROPIEDADES DE LOS FLUIDOS

PROPIEDADES DE LOS FLUIDOS

La fugacidad tiene unidades de presin.

PROPIEDADES DE LOS FLUIDOS

ln f/p = z 1 ln z + 1/RT (RT/VM p ) d VM

PROPIEDADES DE LOS FLUIDOS

PROPIEDADES DE LOS FLUIDOS

PROPIEDADES DE LOS FLUIDOS

en donde el parmetro b representa el volumen de las molculas que

PROPIEDADES DE LOS FLUIDOS

PROPIEDADES DE LOS FLUIDOS

PROPIEDADES DE LOS FLUIDOS

PROPIEDADES DE LOS FLUIDOS

PROPIEDADES DE LOS FLUIDOS

Fig. 4.1 - Diagrama de presin contra volumen para una sustancia

PROPIEDADES DE LOS FLUIDOS

PROPIEDADES DE LOS FLUIDOS

Ahora bien, expresando las ecuacin de estado de van der Waals

PROPIEDADES DE LOS FLUIDOS

combinando las ecuaciones 4.12, 4.13 y 4.16 se obtiene,

Ejemplo 4.1Constantes de la ecuacin de estado de van der Waals.

( 27 )( 10 .732 lb/pg 2 abs ft 3 / lbm mol R )2 ( 963.8 oR )2

(10.732 lb/pg 2 abs ft 3 /lbm mol R )( 963.8 R)

Ecuaciones de Estado para Gases Reales

Se defini una ecuacin de estado como una relacin matemtica entre

Ecuacin de Van der Waals

La ecuacin de Van der Waals corrige

1.- El volumen de las molculas no es despreciable con respecto al volumen

Los valores de a y b son caractersticos de un gas en particular:

a el caso de gases reales.

Calcular las constantes de van der Waals para el 3-metil-hexano.

( 27 )( 10 .732 lb/pg 2 abs ft 3 / lbm mol R )2 ( 963.8 oR )2

( 64 )( 408 .1 lb/pg abs )

(10.732 lb/pg 2 abs ft 3 /lbm mol R )( 963.8 R)

ecuacin de van der Waals puede ser escrita en forma polinomial:

Presin (lb/pg abs.)

Presin (lb/pg2 abs.)

Etapa 2. Calculando los parmetros a y b

(10.732 lb / pg 2 abs ft 3 / lbm mol R ) 2 (666 R ) 2

(10.732 lb / pg 2 abs ft 3 / lbm mol R )( 666 R )

Etapa 3. Calcular los coeficientes A y B aplicando las ecuaciones 4.25 y 4.26

abs ft 6 /( lbm mol)2 ( 185 lb / pg 2 abs )

(1.4496 ft 3 / lbm mol )(185 lb / pg 2 abs )

Etapa 4. Sustituyendo los valores de A y B dentro de la ecuacin:

La EOS de RK expresada en forma polinomial:

En un cilindro cerrado se encuentra propano puro a 100 F teniendo las fases

Etapa 1. Calculando los parmetros a, b, A y B.