Parámetros críticos y factor de comp.

para el agua

Pc := 220.55bar ω := 0.345

Tc := 647.1K

J
R := 8.314
K mol

Calculo de parámetros de la ecuacion de Peng Robinson EOS

T := 773K Temperatura de trabajo

2
2 Tc 2
ac := 0.45723553 R  κ := 0.37464 + 1.54226  ω - 0.26993 ω
Pc

2
 
α( T) := 1 + κ  1 -
T 

  Tc 

Tc
a( T) := ac  α( T) b := 0.07779607 R 
Pc

α( 773K) = 0.844 κ = 0.875

2
 m3  -5 m
3
a( 773K) = 0.506 Pa   b = 1.898  10
 mol  mol

R T a ( T) Presión, según Peng-Robinson

P_PR( V , T) := -
( V - b) V ( V + b ) + b  ( V - b ) en función de V y T
¿Cálculo de Presión teniendo V y T?

3
cm
Vol := 8514 Temp := 773K
mol

 3 
cm 5 Evaluar la función en V y T
P_PR 8514 , 773K = 7.496  10 Pa
 mol  (Chequear en Tabla de Vapor
el valor de P obtenido)

Gráficamente

Vmin := 1.2 b
Se define una variable de rango
para el volumen
V := Vmin , 1.1Vmin .. 1000Vmin

6
2.510

6
210

6
1.510
P_PR( V , 773K)
6
110
(Pa)
5
510

0
3 3 3 3 3 4
2 10 3.610 5.2 10 6.8 10 8.4 10 1 10
V 1000000

(cm3/mol)
 ¿Cómo calcular V si tengo P y T?
El comando Given-Find no funciona con unidades!
Se definen nuevamente las constantes, sin unidades
T := 773 en K

R := 8.314 en J / K mol
-5
a := 0.506 b := 1.898 10 en unidades S.I
(calculados arriba)

κ = 0.875

5
P := 7.5 10 en Pa

-3
V := 8  10 Valor Semilla (m 3 mol -1)

Given

R T a
P =  - 

 (V - b) V ( V + b ) + b  ( V - b ) 

-3
Vol := Find( V) = 8.509  10 En m 3 mol -1
3
Vol  1000000 = 8.509  10 En cm 3 mol -1
 A T=773 K un único volumen
para todas las presiones (fluido supercrítico):

¿Qué pasa cuando a la T de trabajo puedo tener

líquido y/o vapor? P. ej a 600K (327 ºC)

Pc := 220.55bar

Tc := 647.1K

J
R := 8.314
K mol

Calculo de parámetros de la ecuacion de Peng Robinson EOS

2
2 Tc 2
ac := 0.45723553 R  κ := 0.37464 + 1.54226  ω - 0.26993 ω
Pc
 2
 
α( T) := 1 + κ  1 -
T 

  Tc 

Tc
a( T) := ac  α( T) b := 0.07779607 R 
Pc

3
-5 m
κ = 0.875 b = 1.898  10 Vmin := 1.2 b
mol

V := Vmin , 1.1Vmin .. 1000Vmin

R T a ( T)
P_PR( V , T) := -
( V - b) V ( V + b ) + b  ( V - b )

7
310
7
2.710
7
2.410
7
2.110
7
1.810
7
P_PR( V , 600K) 1.510
7
(Pa) 1.210
6
910
6
610
6
310

0
3
0 100 200 300 400 500 600 700 800 900 1 10
V 1000000

(cm3/mol)
A 600 K Peq=1.25x107 Pa

7
310
6
3010
7
2.710

7
2.410

7
2.110

Líquido
7
1.810

7
P_PR ( V , 600K) 1.510

7
1.210
(Pa)
6
910

6
610

Vapor
6
310

0
0
3
0 100 200 300 400 500 600 700 800 900 110
0 V 1000000 1000

Vliq sat = 27 cm3/mol

Vvap sat = 240 cm3/mol

Para P < Peq hay una única fase vapor.

Para P> Peq hay una única fase líquido

¿Cómo calculo V teniendo P y T?

5
Temp := 600 en Pc := 220.10 ω := 0.345
K
Tc := 647.1

R := 8.314

-5
a := 0.64 b := 1.898  10 en Unidades
S.I
6
P := 6  10 en Pa

Esta presión es menor a la Peq a esta T con lo cual hay fase vapor
 -3
Valor Semilla V := 1  10 en m -3 mol -1

Given

R Temp a
P =  - 

 (V - b) V ( V + b ) + b  ( V - b ) 

Vol := Find( V)

-4
Vol = 7.118  10
en m -3 /mol

Vol  1000000 = 711.815 en cm -3 /mol

7
310
6
30 10
7
2.710

7
2.410

7
2.110

Líquido
7
1.810

7
P_PR ( V , 600K) 1.510

7
1.210
(Pa)
6
910

3
6
610
2 Vapor 1
6
310

0
0
3
0 100 200 300 400 500 600 700 800 900 110
0 V  1000000 1000

La única raíz que tiene sentido en estas condiciones de P y T es la (1) que

corresponde al volumen del vapor. Las otras dos carecen de sentido físico.

La raíz que se obtenga como resultado de la resolución de la ecuación de

Peng Robinson depende el valor semilla utilizado para V.
Robinson EOS