Universidad Simón Bolívar

Departamento de Termodinámica y Fenómenos de Transferencia

Termodinámica II (TF-2323)
Profesor: Freddy Figueira

Problema de Ecuaciones Cúbicas

Se tiene el sistema cilindro pistón mostrado en la figura, que contiene inicialmente una mezcla líquido
vapor de agua con una calidad de 0,6. El sistema se calienta hasta que el pistón apenas toca los topes.
Luego se introduce un 1 kg de agua, momento en el cual la temperatura alcanza los 200 ºC. Determine:
1. Masa inicial de agua.
2. Temperatura del agua al finalizar el proceso de calentamiento (Estado 2)
3. Presión final del sistema.

Utilizar Ecuación de Estado de

1m

Peng-Robinson
1m

1 m2

Propiedades para el agua:

kg kJ
Tc  647.1 K Pc  22048.32 kPa ω  0.345 M  18.015 RU  8.31451
kmol kmol·K

RU
La constante específica del agua es: R  R  0.46153
M

Las Constantes de Antoine son: A1  6.532470 A2  3985.439 A3  38.9974

Solución:
Estado 1

P1  100 kPa (El sistema se encuentra a presión atmosférica)

x 1  0.6

Estamos en saturación líquido vapor, se debe aplicar el criterio de Maxwell para determinar los volúmenes
específicos de las fases líquido y vapor, además de la temperatura de saturación. Se plantea el
procedimiento iterativo:

1. Suponer Tsat, para lo cual se emplea la ecuación de Antoine:

 ln
Psat  A2
 = A1  T  A3 Donde T va en Kelvin
 Pc

Despejando:

A2
T=  A3
A1  ln
Psat 

 Pc 
Así, el primer estimado es:

A2
T1   A3 T1  373.1139 K ó T1  273.15  99.9639
 P1 
A1  ln 
 Pc 
2. Se resuelve la ecuación cúbica (Ec. Peng-Robinson):

R T ac α
P= 
vb 2 2
v  2 b v  b

Donde:
2 2
R  Tc R  Tc 3
ac  0.45724  ac  1.84976 b  0.07780  b  1.05385  10
Pc Pc

En forma cúbica:
3 2
v  c2  v  c1  v  c0 = 0

Con:

2
3 R T b  ac α b ac α  2  b  R T 2 R T
c0 = b  c1 =  3 b c2 = b 
P P P

 
2
 2 T 
α = 1  0.37464  1.54226  ω  0.26992  ω   1  
  Tc

Para la primera iteración:

  1 
2
 2 T1 
α1  1  0.37464  1.54226  ω  0.26992  ω  α1  1.46527
  Tc

 2 
 3 R T1  b  ac α1 b  5
c0  b  c0  2.66498  10
 
 P1


 ac α1  2  b R T1 2
c1    3 b  c1  0.02347
 P1

 R T1
c2  b  c2  1.72099
P1

Las raíces de la ecuación cúbica son:

Coeficientes: Raíces

 c0 
 1.2499  10 3  Líquido
 c1   
coef    vv  polyroots( coef ) vv   0.01249 
 c2   1.70725 
1   Gas
 
3
v f1  vv v f1  1.2499  10 v g1  vv v g1  1.70725
0 2

3. Comprobamos con el criterio de Maxwell:

v
 g
1
Psat =   P dv
v g  v f v
f

En este caso (Revisar Guía Relaciones termodinámicas EOS):

  vg  b  ac α  v g   1  2  b  v f   1  2   b 

 R T ln
1
Psat =   ln
vg  vf 
  v f  b  2  2  b  v g   1  2  b  v f   1  2   b 

Comprobamos

  vg1  b  ac α1  vg1   1  2   b  vf1   1  2   b 

 R T1  ln
1
Psat1   ln
v g1  v f1 
  v f1  b  2 2 b  vg1   1  2   b  vf1   1  2   b 

Psat1  95.77063 diferente de P1  100

Suponiendo una nueva temperatura (debe ser mayor, puesto quedamos por debajo de la presión)

T2  101  273.15

 
2
 2 T2 
α2  1  0.37464  1.54226  ω  0.26992  ω   1   α2  1.46304
  Tc

 2 
 3 R T2  b  ac α2 b  5
c0  b  c0  2.6601  10
 
 P1

  ac α2  2  b R T2 2
c1    3 b  c1  0.02342
 P1


R T2
c2  b  c2  1.72577
P1

Recalculemos las raíces del polinomio cúbico:

Coeficientes: Raíces

 c0 
 1.25109  10 3  Líquido
 c1   
coef    vv  polyroots( coef ) vv   0.01242 
 c2   
1  1.7121  Gas
 
3
v f2  vv v f2  1.25109  10 v g2  vv v g2  1.7121
0 2

Comprobemos:

  vg2  b  ac α2  vg2   1  2   b  vf2   1  2   b 

 R T2  ln
1
Psat2   ln
v g2  v f2 
 v f2  b
    vg2   1  2   b  vf2   1  2   b 
 2 2 b
   

Psat2  99.56147

El error es esá por debajo del 1%, podríamos quedarnos con este valor, sin
embargo, realicemos una iteración adicional:

Extrapolemos a P=100

 Psat1   T1  100

T3  linterp     T3  374.26986 K ó T3  273.15  101.11986
 Psat2   T2  

Comprobemos nuevamente:

 
2
 2 T3 
α3  1  0.37464  1.54226  ω  0.26992  ω   1   α2  1.46304
  Tc

 2 
 3 R T3  b  ac α3 b  5
c0  b  c0  2.65954  10
 
 P1

  ac α3  2  b R T3 2
c1    3 b  c1  0.02341
 P1

R T3
c2  b  c2  1.72632
P1

Recalculemos las raíces del polinomio cúbico:

Coeficientes: Raíces

 c0 
 1.25123  10 3  Líquido
 c1   
coef    vv  polyroots( coef ) vv   0.01241 
 c2   1.71266 
1   Gas
 
3
v f3  vv v f3  1.25123  10 v g3  vv v g3  1.71266
0 2

Comprobemos:

  vg3  b  ac α3  vg3   1  2   b  vf3   1  2   b 

 R T2  ln
1
Psat3   ln
v g3  v f3 
  v f3  b  2 2 b  vg3   1  2   b  vf3   1  2   b 

Psat3  99.70522 Nos quedamos con esta temperatura.

Así:

T1  T3 T1  374.26986

El volumen específico de la mezcla es:


v 1  v f3  x 1  v g3  v f3  v 1  1.0281

La masa viene dada por:

3
El volumen total en el estado 1 es: V1  1 m

V1
m1  m1  0.97267 kg
v1
Estado 2

El sistema se expande hasta tocar solo los topes:

La presión es: P2  100 kPa

La mas es la misma: m2  m1
2m

3
El nuevo volumen del sistema es: V2  2 m
V2
El volumen específico es: v 2 
m2
1 m2 3
m
v 2  2.05619
kg

Como la presión es la misma, los volumenes de líquido saturado y vapor saturado son los mismos que se
calcularon para el edo. 1. En este caso se tiene que el estado es vapor sobrecalentado ya que:

v2  vg

Debemos calcular la temperatura en este estado, lamentablemente, para la Ec. de Peng-Robinson, el

cálculo no es directo puesto que la Temperatura no es despejable. Debemos realizar un procedimiento
iterativo.

1. Suponer una temperatura:

T1  400 K

2. Calcular la presión mediante la cúbica:

 
2
 2  T1 
ac 1  0.37464  1.54226  ω  0.26992  ω   1  
P1 
R T1

  Tc

v2  b 2 2
v2  2 b v2  b

3. Verificar

P1  89.21409 diferente de: P2  100

Segunda iteración

T2  450
  
2
 2 T2 
ac 1  0.37464  1.54226  ω  0.26992  ω   1  
P2 
R T2

  Tc

v2  b 2 2
v2  2 b v2  b

P2  100.48534 El error está por debajo del 1%, sin embargo hagamos una iteración más.

Interpolando:

 P1   T1  100
T3  linterp     T3  447.84701 K ó T3  273.15  174.69701 ºC
 P2   T2  

 
2
 2  T3 
ac 1  0.37464  1.54226  ω  0.26992  ω   1  
P3 
R T3

  Tc

v2  b 2 2
v2  2 b v2  b

P3  100.00008 Nos quedamos con este valor

Así: T2  T3 T2  447.84701 K ó T2  273.15  174.69701 ºC

Estado 3
Ahora cambia la masa del sistema y por ende su volumen específico.

V2 3
m
v 3  v 3  1.01385
m2  1 kg

Se sabe:

T3  200  273.15 T3  473.15 K

Nuevamente deberíamos aplicar el criterio de Maxwell para comprobar el estado del sistema.

1. Estimado de presión de saturación por Antoine

ln
Psat  A2
 = A1  T  A3
 Pc 

Psat3  Pc exp A1 

A2 
Psat3  1561.88743
T3  A3 
kPa
 
2. Se resuelve la ecuación cúbica (Ec. Peng-Robinson):

En forma cúbica:
 3 2
v  c2  v  c1  v  c0 = 0
Con:
2
3 R T b  ac α b ac α  2  b  R T 2 R T
c0 = b  c1 =  3 b c2 = b 
P P P

Para la primera iteración:

Ahora el valor de α será constante:

2
  T3 


α3  1  0.37464  1.54226  ω  0.26992  ω
2  1 


Tc 
α3  1.26953


 2 
 3 R T3 b  ac α3 b  6
c0  b  c0  1.42803  10
 
 Psat3


 ac α3  2  b R T3 2 3
c1    3 b  c1  1.2055  10
 Psat3


R  T3
c2  b  c2  0.13876
Psat3

Las raíces de la ecuación cúbica son:

Coeficientes: Raíces

 c0 
 1.41164  10 3  Líquido
 c1   
coef    vv  polyroots( coef ) vv    3
7.80933  10
 c2    Gas
1  0.12954 
 

v f3  vv v g3  vv
0 2

3. Comprobamos con el criterio de Maxwell:

v
 g
1
Psat =   P dv
v g  v f v
f
Comprobamos

  v g3  b  ac α3  v g3   1  2  b  v f3   1  2  b 

 R T3  ln
1
Psat3calc   ln
v g3  v f3 
  vf3  b  2  2  b  v g3   1  2  b  v f3   1  2  b 
 3
Psat3calc  1.55732  10

El error viene dado por:

Psat3calc  Psat3 Menor al 2%, no es necesario realizar una nueva iteración, ya

 100  0.29358
Psat3calc tenemos los volúmenes

3 v g3  0.12954
v f3  1.41164  10

v 3  1.01385 como v 3  v g3 tenemos vapor sobrecalentado, podemos calcular la presión del

sistema evaluando directamente P de la ecuación de estado:

2
  T3 


ac 1  0.37464  1.54226  ω  0.26992  ω
2  1 


R  T3 Tc 
P3  

v3  b 2 2
v3  2 b v3  b

P3  213.33438 kPa