Ejercicio 20 3 Expansor

Descargar como pdf o txt
Descargar como pdf o txt
Está en la página 1de 28

1-. Un expansor de gas con metano inicialmente a 600 Lpca y -100 F.

El valor P1/P2 es de 5:1


El proceso es Adiabatico, pero las irreversibilidades hacen que le trabajo se reduzca a 80%
del valor isentropico. Determine:
a) El trabajo real.
b) El trabajo maximo.
c) La temperatura real final
Utilice el metodo termodinamico como el analitico sabiendo que circula 25 MMpcnd de Metano
Ts =
Ps =

-100
600

F
Lpca

Tc =
Pc =

Pd =
Tasa =
nis =

120
2744
80%

Lpca
Lbmol/hr

PM =

-116.66 F
667 Lpca
16.043 Lbm/Lbmol

Utilizando la grafica de Metano (Mollier) podemos obtener los siguiente datos


Para un proceso adiabatico, tenemos:
Con P1 = 600 Lpca y T1 = -100 F obtenemos
S1 = S2 = 1.00 Btu/LbmF (vapor sobrecalentado)
h1 =

380

Btu/Lbm

Condiciones de salida de la turbina


Con P2 = 120 Lpca y S2 = S1 = 1.00 Btu/LbmF obtenemos
h2 =

340

Btu/Lbm

T2 =

-200

(vapor sobrecalentado)

Calculo del Trabajo Maximo

Wm ax =
W max =

m(h2 h1 )
2544
-43.1

hp

Calculo del Trabajo Real

h2

real

= ( h2 h1 ) is + h1

h2 real =

Wreal =
W real =

348

Btu/Lbm

m(h2 real h1 )
2544
-34.5

hp

Preparado por Ing. Edwards Nava

14/06/2007

Calculo del Temperatura real de salida


Con P2 = 120 Lpca y h2real = 348 Btu/Lbm obtenemos
T2 real =

-200

Metodo Analitoco
Ts =

-100

Ps =

600

Lpca

Pd =

120

Lpca

r=

0.200

K@Tprom=

1.3371

Tomado Tecnical Data Book - Petroleum Refining. Pag 7-7


K 1

)(rreal _ 2 ) K

T d 2 = (T s 2 + 460
Td =

-220

Tprom =
F

460

-160

Metodo de Brill, J.P. Y Begg, H.D,


Zs =
0.644
Zd =
0.873

H is =

1545 Z aver T s
PM ( K 1) K

His =

K 1

r K 1

-34788.77 Lbf.Pie/Lbm

Calculo del Trabajo Maximo

Wm ax =
W max =

m H is
33000
-48.21

HP

Calculo del Trabajo Real

Wreal =
W real =

m H is is
33000
-38.57

HP

Preparado por Ing. Edwards Nava

14/06/2007

2-. Un turbo-expansor trabaja con gas natural de composicion conocida inicialmente a 800 Lpca y 20 F.
El valor P1/P2 es de 2:1. El proceso tiene una eficiencia adiabatica del 82%. Determine:
a) El trabajo real.
b) El trabajo maximo.
c) La temperatura real final

Comp
C1

% yi
80.8

C2
C3
iC4
nC4
iC5
nC5
nC6
C7+

12.2
4.2
0.7
1.1
0.3
0.4
0.1
0.2
100

PMC7+ =

122

C7+ =

0.73

Ts =
Ps =
Pd =
Tasa =
Moles =

20
800
400
500
52219

Lbm/lbmol

F
Lpca
Lpca
MMpcnd
Lbmol/hr

Determinacin de las Propiedades del C7+


Correlacin de Winn-Sim Daubert

PM 0.9371

Tb =
5
1 .4350476 10
Tb C7+ =

725.5

2.3776

Tc = exp 3 .9934718 Tb 0.08615 0.04614


Tc C7+ =

1034.9

Pc = 3 .48242 10 9 Tb 2.3177 2.4853


Pc C7+ =

wC 7 + =

w C7+ =

373.4

Lpca

3 Log sPc
7
14 .7
sTc

1
Tbn

0.4117

Propiedades de los compuestos


Comp
C1
C2
C3
iC4
nC4
iC5
nC5
nC6
C7+

yi
0.808
0.122
0.042
0.007
0.011
0.003
0.004
0.001
0.002
1.00

w
0.0108
0.0972
0.1515
0.1852
0.1981
0.2286
0.2510
0.2990
0.4117

Tci
343.3
550.1
665.9
734.4
765.5
829.0
845.7
911.8
1034.9

Pci
667.0
707.8
615.0
527.9
548.8
490.4
488.1
439.5
373.4

Tbn
201.3
332.5
416.3
470.8
491.1
540.1
556.9
615.7
725.5

Preparado por Ing. Edwards Nava

14/06/2007

Calculo de las fases a la entrada del Turbo-expansor


Inicialmente se estimo un Pk = 1500 Lpca, pero el proceso iterativo nos llevo hasta un Pk = 2000
Lpca. Con la temperatura 20F y presion 800 Lpca obtenemos los Ki pero no disponemos de KC7+,
por lo que debemos realizar un ajuste (lnea recta) con Ki vs Tci2

Pk

fw (xi)

Ki

sTcC2+

Pk (Hadden)

y = 1
x = 1

Comparar

donde:

x
sTc C 2 + =

PM i Tc i

xi =

2
n

x PM
i

zi
( K i 1) V + 1

y i = K i xi

A la temperatura de 20 F y Presion de 800 Lpca


Comp
C1
C2
C3
iC4
nC4
iC5
nC5
nC6
C7+

zi
0.8080
0.1220
0.0420
0.0070
0.0110
0.0030
0.0040
0.0010
0.0020

sTcC2+ =

670.0
210.0

Ki
2.8000
0.5400
0.1900
0.0890
0.0650
0.0315
0.0265
0.0097
0.0015

V=
xi
0.298
0.217
0.183
0.053
0.100
0.038
0.054
0.017
0.040
1.000

Tci2 R2
117882.4
302577
443449.4
539358
586005.6
687174.7
715208.5
831379.2
1070967

0.951
yi
0.834
0.117
0.035
0.005
0.006
0.001
0.001
0.000
0.000
1.000

xiPMi

xiPMiTci

1.957
1.263
0.309
0.639
0.174
0.289
0.072
0.172
4.876

1076.6
841.0
226.7
489.4
144.5
244.1
65.8
178.4
3266.5

R
F

x
sTc

C 2+

PM

Tc i

2
n

PM

En la Fig. 23-2, Cap.23 del GPSA, el valor de la (Tc)C2+ se encuentra entre la (Tc)C3 = 205.92 F
y la (Tc)iC4 = 274.41F. En la Fig. 25-11 del GPSA (graf de Hadden), se construye una curva la
cual estar entre los componentes C3 y CiC4, y con la (Tc)C2+ se corta la curva de unin entre el
C3 y el iC4, luego este punto se une con el punto del componente C1; se traza una curva y con
la T = 20 F se corta la curva y se obtiene una Pk = 1500 lpca.

Pk

Ki

fw (xi)

sTcC2+

Pk (Hadden)

y = 1
x = 1

Comparar

Preparado por Ing. Edwards Nava

14/06/2007

Calculo de la Entalpa a la entrada del Turbo-expansor

Fase gaseosa
Temp =
presin =
C7+ =

20.0
800
0.73

Comp
C1
C2
C3
iC4
nC4
iC5
nC5
nC6
C7+

yi
0.8340
0.1172
0.0348
0.0047
0.0065
0.0012
0.0014
0.0002
0.0001

sTr =
sPr =

1.247
1.198

F
Lpca

PMi
16.04
30.07
44.10
58.12
58.12
72.15
72.15
86.18
122.00

hi
239.2
146.9
124.5
118.0
124.5
118.0
120.5
107.1
102.0

yiPMihi
3199.9
517.6
191.0
32.1
46.9
10.2
12.5
1.6
0.8
4012.5

(0)

H H P ,T
R sTc

(I )

H H P ,T

R sTc

yiTci
286.4
64.4
23.2
3.4
5.0
1.0
1.2
0.2
0.1
384.8

1.10

0.28

H H P ,T
H H P ,T

+ w
= R sTc

R sTc
R sTc

=
846.73 Btu/Lbmol

Hp , t = Hm H H P ,T

(0 )

H H P ,T

H H P ,T

Hp,t =

3165.8

(I )

yiPci
556.3
82.9
21.4
2.5
3.6
0.6
0.7
0.1
0.0
668.0

yiwi
0.0090
0.0114
0.0053
0.0009
0.0013
0.0003
0.0004
0.0001
0.0000
0.0285

PMiewi
16.217
33.140
51.310
69.949
70.857
90.681
92.735

184.144

Grafico
TciHi0.5/PMi
331.0
221.7
168.5
137.2
146.9
124.8
128.6
85.7

Btu/Lbmol

Calculo de la entropa a la entrada del Turbo-expansor

Fase gaseosa
Temp =
Presin =

20
800

Comp
C1
C2
C3
iC4
nC4
iC5
nC5
nC6
C7+

yi
0.8340
0.1172
0.0348
0.0047
0.0065
0.0012
0.0014
0.0002
0.0001

F
Lpca
PMi
16.043
30.070
44.097
58.123
58.123
72.150
72.150
86.177
122.00

si
3.049
1.956
1.545
1.316
1.255
1.185
1.166
1.125
1.099

yiPMisi
40.79
6.89
2.37
0.36
0.47
0.10
0.12
0.02
0.01
51.1

yilnyi
-0.151
-0.251
-0.117
-0.025
-0.033
-0.008
-0.009
-0.001
-0.001
-0.597

yiTci
286.4
64.4
23.2
3.4
5.0
1.0
1.2
0.2
0.1
384.8

Preparado por Ing. Edwards Nava

yiPci
556.3
82.9
21.4
2.5
3.6
0.6
0.7
0.1
0.0
668.0

xiwi
0.0090
0.0114
0.0053
0.0009
0.0013
0.0003
0.0004
0.0001
0.0000
0.0285

14/06/2007

sTr =
sPr =

S S P ,T

1.247
1.198

S S P ,T

(0)

0.50

Grafico
Tc.ln(PM)/ln(Tbn)

180
322
418
485
502
564
572

(I )

0.24

Sm = Yi Si PMi R Yi Ln (Yi )
Sm =

S S P ,T

S S P ,T

52.32

43.374

38.311
15.521
11.043
9.012
8.257
7.486
6.896

Btu/Lbmol R

S S P ,T ( 0 )
S S P ,T

+ w
= R
R
R

= 8.94 Btu/Lbmol R

Sp,t =

Ln(PM).ln(1/w).S

(I )

+ Ln ( P / 14 .7 )

Btu/Lbmol R

Calculo de la Temperatura de rocio a la salida del Turbo-expansor


Inicialmente se estimo un Pk = 2000 Lpca, pero el proceso iterativo nos llevo hasta un Pk = 1500
Lpca. Supondremos una Temp como temperatura de roci hasta cumplir zi/Ki = 1; pero como no
disponemos de KC7+ debemos realizar un ajuste (lnea recta) con Ki vs Tci2

Pk

xi = zi/Ki

fw (xi)

sTcC2+

Pk (Hadden)

zi/Ki = 1

Comparar
donde:

x
sTc C 2 + =

PM i Tc i

2
n

x PM
i

Suponemos la primera temperatura de 20 F


Comp
C1
C2
C3
iC4
nC4
iC5
nC5
nC6
C7+

zi
0.8340
0.1172
0.0348
0.0047
0.0065
0.0012
0.0014
0.0002
0.0001

Ki
4.7000
0.8000
0.2300
0.0960
0.0700
0.0280
0.0240
0.0077
0.0009

Tci2 R2
117882
302577
443449
539358
586006
687175
715208
831379
1070967

zi/Ki
0.1775
0.1464
0.1513
0.0487
0.0925
0.0430
0.0599
0.0218
0.0657
0.807

dado que zi/Ki < 1 debemos asumir temperaturas inferiores a la supuesta (20 F)

Preparado por Ing. Edwards Nava

14/06/2007

Suponemos la segunda temperatura de 0 F


Comp
C1
C2
C3
iC4
nC4
iC5
nC5
nC6
C7+

zi
0.8340
0.1172
0.0348
0.0047
0.0065
0.0012
0.0014
0.0002
0.0001

Ki
4.4000
0.6700
0.1750
0.0700
0.0520
0.0190
0.0146
0.0047
0.0005

Tci2 R2
117882
302577
443449
539358
586006
687175
715208
831379
1070967

zi/Ki
0.190
0.175
0.199
0.067
0.125
0.063
0.098
0.036
0.126
1.078

dado que zi/Ki > 1 debemos asumir temperaturas superiores a la supuesta (0 F)


Se observa muy cercana a 1.00 haremos una interpolacin entre los Ki @ 20 F y los Ki @ 0 F
Plantearemos la ecuacin matemtica en excel con el fin de simplificar los calculos
La interpolacin la haremos con Ln(Ki) vs T y luego las ln(Ki) obtenidas se le deber aplicar el
exponencial base 2.71. La ecuacin de interpolacin es la siguiente

T T100
K 3 = exp ln( K 100 ) + (ln( K 80 ) ln( K 100 )) *
T80 T100

Obtenindose la siguiente temperatura de roco de


Comp
C1
C2
C3
iC4
nC4
iC5
nC5
nC6
C7+

zi
0.8340
0.1172
0.0348
0.0047
0.0065
0.0012
0.0014
0.0002
0.0001

sTcC2+ =

613.1
153.1

Ki
4.4729
0.7003
0.1873
0.0757
0.0560
0.0209
0.0165
0.0053
0.0006

R
F

xi=yi/Ki
0.1865
0.1673
0.1858
0.0618
0.1157
0.0575
0.0870
0.0316
0.1069
1.0000

4.98

Tci2 R2
117882
302577
443449
539358
586006
687175
715208
831379
1070967

xiPMi

xiPMiTci

3.523
1.535
0.272
0.377
0.087
0.104
0.014
0.007
5.918

1937.9
1021.9
199.6
288.2
72.0
87.7
13.2
7.7
3628

x
sTc

C 2+

PM

Tc i

2
n

PM

En la Fig. 23-2, Cap.23 del GPSA, el valor de la (Tc)C2+ se encuentra entre la (Tc)C2 = 90.07 F
y la (Tc)C3 = 205.92F. En la Fig. 25-11 del GPSA (graf de Hadden), se construye una curva la
cual estar entre los componentes C2 y C3, y con la (Tc)C2+ se corta la curva de unin entre el
C2 y el C3, luego este punto se une con el punto del componente C1; se traza una curva y con
la T = 5 F se corta la curva y se obtiene una Pk = 1500 lpca.

Pk

xi = zi/Ki

fw (xi)

sTc2+

Pk (Hadden)

zi/Ki = 1

Comparar

Preparado por Ing. Edwards Nava

14/06/2007

Temperatura de salida del Turbo-expansor (-20F)


Inicialmente se estimo un Pk = 2000 Lpca, pero el proceso iterativo nos llevo hasta un Pk = 1500
Lpca. Con la temperatura -20F y presion 400 Lpca obtenemos los Ki pero no disponemos de KC7+,
por lo que debemos realizar un ajuste (lnea recta) con Ki vs Tci2

Pk

fw (xi)

Ki

sTcC2+

Pk (Hadden)

y = 1
x = 1

Comparar

donde:

x
sTc C 2 + =

PM i Tc i

zi
( K i 1) V + 1

xi =

x i PM i

y i = K i xi

Suponemos la primera temperatura de -20 F


Comp
C1
C2
C3
iC4
nC4
iC5
nC5
nC6
C7+

zi
0.8340
0.1172
0.0348
0.0047
0.0065
0.0012
0.0014
0.0002
0.0001

sTcC2+ =

649.6
189.6

Ki
4.1000
0.5100
0.1200
0.0490
0.0340
0.0120
0.0091
0.0027
0.0015

V=
xi
0.206
0.226
0.258
0.072
0.129
0.042
0.056
0.009
0.003
1.000

Tci2 R2
117882
302577
443449
539358
586006
687175
715208
831379
1070967

0.983
yi
0.845
0.115
0.031
0.004
0.004
0.001
0.001
0.000
0.000
1.000

xiPMi

xiPMiTci

3.626
7.758
3.168
7.478
2.436
4.009
0.620
0.286
29.38

1245.1
4267.3
2109.8
5492.2
1865.0
3323.3
524.7
260.8
19088

R
F

x
sTc

C 2+

PM

Tc i

2
n

PM

En la Fig. 23-2, Cap.23 del GPSA, el valor de la (Tc)C2+ se encuentra entre la (Tc)C2 = 90.07 F
y la (Tc)C3 = 205.92F. En la Fig. 25-11 del GPSA (graf de Hadden), se construye una curva la
cual estar entre los componentes C2 y C3, y con la (Tc)C2+ se corta la curva de unin entre el
C2 y el C3, luego este punto se une con el punto del componente C1; se traza una curva y con
la T = -20 F se corta la curva y se obtiene una Pk = 1500 lpca.

Pk

Ki

fw (xi)

sTcC2+

Pk (Hadden)

y = 1
x = 1

Comparar

Preparado por Ing. Edwards Nava

14/06/2007

Calculo de la Entalpa a la salida del Turbo-expansor (-20 F)

Fase gaseosa
Temp =
presin =
C7+ =

-20.0
400
0.73

Comp
C1
C2
C3
iC4
nC4
iC5
nC5
nC6
C7+

yi
0.8448
0.1153
0.0310
0.0035
0.0044
0.0005
0.0005
0.0000
0.0000

sTr =
sPr =

1.155
0.598

F
Lpca

PMi
16.04
30.07
44.10
58.12
58.12
72.15
72.15
86.18
122.00

hi
218.2
131.6
110.1
103.6
110.1
103.6
106.1
93.5
88.1

yiPMihi
2957.2
456.2
150.3
21.2
28.0
3.8
3.9
0.2
0.1
3620.8

(0)

H H P ,T
R sTc

(I )

H H P ,T

R sTc

yiTci
290.1
63.4
20.6
2.6
3.3
0.4
0.4
0.0
0.0
380.9

0.50

0.23

H H P ,T
H H P ,T

+ w
= R sTc

R sTc
R sTc

=
382.89 Btu/Lbmol

Hp , t = Hm H H P ,T

(0 )

H H P ,T

H H P ,T

Hp,t =

3237.9

(I )

yiPci
563.5
81.6
19.0
1.9
2.4
0.2
0.2
0.0
0.0
668.9

yiwi
0.00912
0.01121
0.00469
0.00065
0.00087
0.00011
0.00013
0.00001
0.00000
0.02679

PMiewi
16.217
33.140
51.310
69.949
70.857
90.681
92.735

184.144

Grafico
TciHi0.5/PMi
316.1
209.9
158.5
128.6
138.2
116.9
120.7
79.6

Btu/Lbmol

Fase Liquida
Temp =
presin =
C7+ =
Comp
C1
C2
C3
iC4
nC4
iC5
nC5
nC6
C7+

-20.0
400
0.73
xi
0.2061
0.2260
0.2580
0.0718
0.1287
0.0419
0.0556
0.0086
0.0033

F
Lpca

PMi
16.04
30.07
44.10
58.12
58.12
72.15
72.15
86.18
122.00

hi
218.2
131.6
110.1
103.6
110.1
103.6
106.1
93.5
88.1

xiPMihi
721.3
894.6
1252.6
432.7
823.4
313.3
425.4
69.3
35.7
4968.2

xiTci
70.7
124.3
171.8
52.8
98.5
34.7
47.0
7.8
3.4
611.2

xiPci
137.4
160.0
158.7
37.9
70.6
20.6
27.1
3.8
1.2
617.3

Preparado por Ing. Edwards Nava

xiwi
0.00223
0.02197
0.03909
0.01331
0.02549
0.00958
0.01395
0.00257
0.00137
0.12954

14/06/2007

sTr =
sPr =

0.720
0.648

H H P ,T

R sTc

(0)

H H P ,T

R sTc

(I )

4.70

6.60

H H P ,T
H H P ,T

+ w
= R sTc

R sTc
R sTc

=
6742.44 Btu/Lbmol

Hp , t = Hm H H P ,T

(0 )

H H P ,T

H H P ,T

Hp,t =

(I )

-1774.2 Btu/Lbmol

H Mezcla = H V V + H L L
Hmezcla =

3153.1

Btu/Lbmol

Calculo de la entropa a la salida del Turbo-expansor (-20 F)

Fase gaseosa
Temp =
Presin =

-20
400

Comp
C1
C2
C3
iC4
nC4
iC5
nC5
nC6
C7+

yi
0.8448
0.1153
0.0310
0.0035
0.0044
0.0005
0.0005
0.0000
0.0000

sTr =
sPr =

F
Lpca
PMi
16.043
30.070
44.097
58.123
58.123
72.150
72.150
86.177
122.00

si
3.004
1.923
1.514
1.284
1.224
1.154
1.136
1.095
1.068

S S P ,T

1.1552
0.5980

S S P ,T

yiPMisi
40.71
6.67
2.07
0.26
0.31
0.04
0.04
0.00
0.00
50.1

yilnyi
-0.142
-0.249
-0.108
-0.020
-0.024
-0.004
-0.004
0.000
0.000
-0.551

yiTci
290.1
63.4
20.6
2.6
3.3
0.4
0.4
0.0
0.0
380.9

S S P ,T

S S P ,T

Sp,t =

51.20

0.28

Grafico
Tc.ln(PM)/ln(Tbn)

(I )

0.21

180
322
418
485
502
564
572

Ln(PM).ln(1/w).S

37.746
15.258
10.821
8.796
8.053
7.289
6.718

Btu/Lbmol R

S S P ,T
S S P ,T

+ w
= R
R
R

= 7.13 Btu/Lbmol R

(0 )

44.069

yiwi
0.00912
0.01121
0.00469
0.00065
0.00087
0.00011
0.00013
0.00001
0.00000
0.02679

(0)

Sm = Yi Si PMi R Yi Ln (Yi )
Sm =

yiPci
563.5
81.6
19.0
1.9
2.4
0.2
0.2
0.0
0.0
668.9

(I )

+ Ln ( P / 14 .7 )

Btu/Lbmol R

Preparado por Ing. Edwards Nava

14/06/2007

Fase Liquida
Temp =
Presin =

-20
400

Comp
C1
C2
C3
iC4
nC4
iC5
nC5
nC6
C7+

xi
0.2061
0.2260
0.2580
0.0718
0.1287
0.0419
0.0556
0.0086
0.0033

sTr =
sPr =

F
Lpca
PMi
16.043
30.070
44.097
58.123
58.123
72.150
72.150
86.177
122.00

0.7199
0.6480

si
3.004
1.923
1.514
1.284
1.224
1.154
1.136
1.095
1.068

xiPMisi
9.93
13.07
17.23
5.36
9.16
3.49
4.55
0.81
0.43
64.0

S S P ,T

(0)

S S P ,T

(I )

xilnxi
-0.325
-0.336
-0.350
-0.189
-0.264
-0.133
-0.161
-0.041
-0.019
-1.818

4.70

7.10

xiTci
70.7
124.3
171.8
52.8
98.5
34.7
47.0
7.8
3.4
611.2

xiPci
137.4
160.0
158.7
37.9
70.6
20.6
27.1
3.8
1.2
617.3

xiwi
0.00223
0.02197
0.03909
0.01331
0.02549
0.00958
0.01395
0.00257
0.00137
0.12954

Sm = Yi Si PMi R Yi Ln (Yi )
Sm =

S S P ,T

S S P ,T

Sp,t =

67.65

Btu/Lbmol R

S S P ,T ( 0 )
S S P ,T

=
R
+ w

R
R

= 17.72 Btu/Lbmol R

49.924

(I )

+ Ln ( P / 14 .7 )

Btu/Lbmol R

S Mezcla = S V V + S L L
Smezcla =

44.168

Btu/Lbmol R

Temperatura de salida del Turbo-expansor (-40F)


Inicialmente se estimo un Pk = 2000 Lpca, pero el proceso iterativo nos llevo hasta un Pk = 1500
Lpca. Con la temperatura -40F y presion 400 Lpca obtenemos los Ki pero no disponemos de KC7+,
por lo que debemos realizar un ajuste (lnea recta) con Ki vs Tci2

Pk

fw (xi)

Ki

sTcC2+

Pk (Hadden)

y = 1
x = 1

Comparar

donde:

x
sTc C 2 + =

PM i Tc i

xi =

2
n

x i PM i

zi
( K i 1) V + 1

y i = K i xi

Preparado por Ing. Edwards Nava

14/06/2007

Suponemos la primera temperatura de -40 F


Comp
C1
C2
C3
iC4
nC4
iC5
nC5
nC6
C7+

zi
0.8340
0.1172
0.0348
0.0047
0.0065
0.0012
0.0014
0.0002
0.0001

sTcC2+ =

682.7
222.7

Ki
3.8000
0.3900
0.0800
0.0350
0.0220
0.0075
0.0051
0.0012
0.0001

V=
xi
0.227
0.280
0.283
0.058
0.096
0.022
0.028
0.004
0.001
1.000

Tci2 R2
117882.4
302577
443449.4
539358
586005.6
687174.7
715208.5
831379.2
1070967

0.953
yi
0.864
0.109
0.023
0.002
0.002
0.000
0.000
0.000
0.000
1.000

xiPMi

xiPMiTci

8.419
12.489
3.399
5.572
1.616
2.015
0.303
0.160
33.974

4631.3
8316.8
2496.4
4265.3
1339.5
1704.3
276.4
165.4
23195

R
F

x
sTc

C 2+

PM

Tc i

PM

En la Fig. 23-2, Cap.23 del GPSA, el valor de la (Tc)C2+ se encuentra entre la (Tc)C3 = 205.92 F
y la (Tc)iC4 = 274.41F. En la Fig. 25-11 del GPSA (graf de Hadden), se construye una curva la
cual estar entre los componentes C2 y C3, y con la (Tc)C2+ se corta la curva de unin entre el
C3 y el CiC4 luego este punto se une con el punto del componente C1; se traza una curva y con
la T = -40 F se corta la curva y se obtiene una Pk = 1500 lpca.

Pk

fw (xi)

Ki

sTcC2+

Pk (Hadden)

y = 1
x = 1

Comparar

Calculo de la Entalpa a la salida del Turbo-expansor (-40 F)

Fase gaseosa
Temp =
presin =
C7+ =
Comp
C1
C2
C3
iC4
nC4
iC5
nC5
nC6
C7+

-40.0
400
0.73
yi
0.8637
0.1092
0.0227
0.0020
0.0021
0.0002
0.0001
0.0000
0.0000

F
Lpca

PMi
16.04
30.07
44.10
58.12
58.12
72.15
72.15
86.18
122.00

hi
207.8
124.2
103.2
96.7
103.2
96.7
99.2
87.0
81.5

yiPMihi
2879.8
407.9
103.2
11.5
12.7
1.2
1.0
0.0
0.0
3417.3

yiTci
296.5
60.1
15.1
1.5
1.6
0.1
0.1
0.0
0.0
375.1

yiPci
576.1
77.3
13.9
1.1
1.2
0.1
0.1
0.0
0.0
669.7

Preparado por Ing. Edwards Nava

yiwi
0.00933
0.01061
0.00343
0.00038
0.00042
0.00004
0.00004
0.00000
0.00000
0.02425

14/06/2007

sTr =
sPr =

(0)

H H P ,T
R sTc

(I )

H H P ,T

R sTc

1.120
0.597

0.55

0.26

H H P ,T ( 0 )
H H P ,T

+ w
= R sTc

R
sTc

R sTc

=
414.38 Btu/Lbmol

Hp , t = Hm H H P ,T

H H P ,T

H H P ,T

Hp,t =

3002.9

(I )

PMiewi
16.217
33.140
51.310
69.949
70.857
90.681
92.735

184.144

Grafico
TciHi0.5/PMi
308.5
203.9
153.4
124.3
133.8
113.0
116.8
76.6

Btu/Lbmol

Fase Liquida
Temp =
presin =
C7+ =
Comp
C1
C2
C3
iC4
nC4
iC5
nC5
nC6
C7+

sTr =
sPr =

-40.0
400
0.73
xi
0.2273
0.2800
0.2832
0.0585
0.0959
0.0224
0.0279
0.0035
0.0013

F
Lpca

PMi
16.04
30.07
44.10
58.12
58.12
72.15
72.15
86.18
122.00

0.7195
0.6308

hi
207.8
124.2
103.2
96.7
103.2
96.7
99.2
87.0
81.5

xiPMihi
757.8
1046.0
1289.4
328.9
575.3
156.3
200.0
26.4
13.0
4393.1

H H P ,T

R sTc

(0)

H H P ,T

R sTc

(I )

xiTci
78.0
154.0
188.6
43.0
73.4
18.6
23.6
3.2
1.4
583.7

4.77

6.80

H H P ,T ( 0 )
H H P ,T

=
R

sTc

+ w

R
sTc

R sTc

=
6445.31 Btu/Lbmol

Hp , t = Hm H H P ,T

H H P ,T

H H P ,T

Hp,t =

(I )

xiPci
151.6
198.2
174.2
30.9
52.6
11.0
13.6
1.5
0.5
634.1

xiwi
0.00245
0.02722
0.04291
0.01083
0.01899
0.00512
0.00701
0.00105
0.00054
0.11612

-2052.2 Btu/Lbmol

H Mezcla = H V V + H L L
Hmezcla =

2767.3

Btu/Lbmol

Preparado por Ing. Edwards Nava

14/06/2007

Calculo de la entropa a la salida del Turbo-expansor (-40 F)

Fase gaseosa
Temp =
Presin =

-40
400

Comp
C1
C2
C3
iC4
nC4
iC5
nC5
nC6
C7+

yi
0.8637
0.1092
0.0227
0.0020
0.0021
0.0002
0.0001
0.0000
0.0000

sTr =
sPr =

F
Lpca
PMi
16.043
30.070
44.097
58.123
58.123
72.150
72.150
86.177
122.00

si
2.980
1.906
1.498
1.268
1.209
1.139
1.121
1.080
1.052

S S P ,T

1.1198
0.5973

S S P ,T

yiPMisi
41.29
6.26
1.50
0.15
0.15
0.01
0.01
0.00
0.00
49.4

yilnyi
-0.127
-0.242
-0.086
-0.013
-0.013
-0.001
-0.001
0.000
0.000
-0.483

yiTci
296.5
60.1
15.1
1.5
1.6
0.1
0.1
0.0
0.0
375.1

S S P ,T

S S P ,T

Sp,t =

50.33

0.34

Grafico
Tc.ln(PM)/ln(Tbn)

(I )

0.27

Ln(PM).ln(1/w).S

180
322
418
485
502
564
572

37.448
15.123
10.707
8.685
7.950
7.189
6.629

xiPci
151.6
198.2
174.2
30.9
52.6
11.0
13.6
1.5
0.5
634.1

xiwi
0.00245
0.02722
0.04291
0.01083
0.01899
0.00512
0.00701
0.00105
0.00054
0.11612

Btu/Lbmol R

S S P ,T
S S P ,T

+ w
= R
R
R

= 7.24 Btu/Lbmol R

(0 )

43.089

yiwi
0.00933
0.01061
0.00343
0.00038
0.00042
0.00004
0.00004
0.00000
0.00000
0.02425

(0)

Sm = Yi Si PMi R Yi Ln (Yi )
Sm =

yiPci
576.1
77.3
13.9
1.1
1.2
0.1
0.1
0.0
0.0
669.7

(I )

+ Ln ( P / 14 .7 )

Btu/Lbmol R

Fase Liquida
Temp =
Presin =

-40
400

Comp
C1
C2
C3
iC4
nC4
iC5
nC5
nC6
C7+

xi
0.2273
0.2800
0.2832
0.0585
0.0959
0.0224
0.0279
0.0035
0.0013

F
Lpca
PMi
16.043
30.070
44.097
58.123
58.123
72.150
72.150
86.177
122.00

si
2.980
1.906
1.498
1.268
1.209
1.139
1.121
1.080
1.052

xiPMisi
10.87
16.05
18.71
4.31
6.73
1.84
2.26
0.33
0.17
61.3

xilnxi
-0.337
-0.356
-0.357
-0.166
-0.225
-0.085
-0.100
-0.020
-0.009
-1.655

xiTci
78.0
154.0
188.6
43.0
73.4
18.6
23.6
3.2
1.4
583.7

Preparado por Ing. Edwards Nava

14/06/2007

sTr =
sPr =

0.7195
0.6308

S S P ,T

(0)

S S P ,T

(I )

4.75

7.20

Sm = Yi Si PMi R Yi Ln (Yi )
Sm =

64.55

S S P ,T

S S P ,T

Sp,t =

Btu/Lbmol R

S S P ,T ( 0 )
S S P ,T

+ w
= R
R
R

= 17.65 Btu/Lbmol R

46.898

(I )

+ Ln ( P / 14 .7 )

Btu/Lbmol R

S Mezcla = S V V + S L L
Smezcla =

43.266

Btu/Lbmol R

Temperatura de salida del Turbo-expansor (-60F)


Inicialmente se estimo un Pk = 2000 Lpca, pero el proceso iterativo nos llevo hasta un Pk = 1500
Lpca. Con la temperatura -60F y presion 400 Lpca obtenemos los Ki pero no disponemos de KC7+,
por lo que debemos realizar un ajuste (lnea recta) con Ki vs Tci2

Pk

fw (xi)

Ki

sTcC2+

Pk (Hadden)

y = 1
x = 1

Comparar

donde:

x
sTc C 2 + =

PM i Tc i

xi =

2
n

x PM
i

zi
( K i 1) V + 1

y i = K i xi

Suponemos la primera temperatura de -60 F


Comp
C1
C2
C3
iC4
nC4
iC5
nC5
nC6
C7+

zi
0.8340
0.1172
0.0348
0.0047
0.0065
0.0012
0.0014
0.0002
0.0001

sTcC2+ =

681.0
221.0

Ki
3.7000
0.2850
0.0550
0.0235
0.0150
0.0042
0.0028
0.0005
0.0000

R
F

V=
xi
0.239
0.343
0.267
0.046
0.069
0.014
0.018
0.002
0.001
1.000

Tci2 R2
117882
302577
443449
539358
586006
687175
715208
831379
1070967

0.921
yi
0.885
0.098
0.015
0.001
0.001
0.000
0.000
0.000
0.000
1.000

xiPMi

xiPMiTci

5.498
8.043
2.040
4.036
0.839
1.264
0.152
0.066
21.94

3024.5
5355.7
1497.9
3089.5
695.9
1068.8
138.2
68.6
14939

x
sTc

C 2+

PM

Tc i

2
n

PM

Preparado por Ing. Edwards Nava

14/06/2007

En la Fig. 23-2, Cap.23 del GPSA, el valor de la (Tc)C2+ se encuentra entre la (Tc)C3 = 205.92 F
y la (Tc)iC4 = 274.41F. En la Fig. 25-11 del GPSA (graf de Hadden), se construye una curva la
cual estar entre los componentes C2 y C3, y con la (Tc)C2+ se corta la curva de unin entre el
C3 y el CiC4 luego este punto se une con el punto del componente C1; se traza una curva y con
la T = -60 F se corta la curva y se obtiene una Pk = 1500 lpca.

Pk

fw (xi)

Ki

sTcC2+

Pk (Hadden)

y = 1
x = 1

Comparar

Calculo de la Entalpa a la salida del Turbo-expansor (-60 F)

Fase gaseosa
Temp =
presin =
C7+ =

-60.0
400
0.73

Comp
C1
C2
C3
iC4
nC4
iC5
nC5
nC6
C7+

yi
0.8854
0.0977
0.0147
0.0011
0.0010
0.0001
0.0000
0.0000
0.0000

sTr =
sPr =

1.083
0.597

F
Lpca

PMi
16.04
30.07
44.10
58.12
58.12
72.15
72.15
86.18
122.00

hi
197.6
117.0
96.6
90.1
96.6
90.1
92.6
80.8
75.2

yiPMihi
2806.4
343.8
62.7
5.7
5.8
0.4
0.3
0.0
0.0
3225.1

H H P ,T

R sTc

(0)

H H P ,T

R sTc

(I )

yiTci
304.0
53.7
9.8
0.8
0.8
0.1
0.0
0.0
0.0
369.2

0.60

0.34

H H P ,T ( 0 )
H H P ,T

=
R

sTc

+ w

sTc

R sTc

=
445.35 Btu/Lbmol

Hp , t = Hm H H P ,T

H H P ,T

H H P ,T

Hp,t =

2779.8

(I )

yiPci
590.5
69.1
9.0
0.6
0.6
0.0
0.0
0.0
0.0
669.9

yiwi
0.00956
0.00949
0.00223
0.00020
0.00021
0.00001
0.00001
0.00000
0.00000
0.02172

PMiewi
16.217
33.140
51.310
69.949
70.857
90.681
92.735
184.144

Grafico
TciHi0.5/PMi
300.8
197.9
148.4
119.9
129.4
109.1
112.8
73.5

Btu/Lbmol

Preparado por Ing. Edwards Nava

14/06/2007

Fase Liquida
Temp =
presin =
C7+ =

-60.0
400
0.73

Comp
C1
C2
C3
iC4
nC4
iC5
nC5
nC6
C7+

xi
0.2393
0.3427
0.2675
0.0463
0.0694
0.0144
0.0175
0.0021
0.0008

sTr =
sPr =

0.707
0.619

F
Lpca

PMi
16.04
30.07
44.10
58.12
58.12
72.15
72.15
86.18
122.00

hi
197.6
117.0
96.6
90.1
96.6
90.1
92.6
80.8
75.2

xiPMihi
758.5
1206.3
1139.2
242.2
389.8
93.9
117.0
14.6
7.1
3968.5

H H P ,T

R sTc

(0)

H H P ,T

R sTc

(I )

xiTci
82.2
188.5
178.1
34.0
53.2
12.0
14.8
1.9
0.8
565.4

4.77

6.80

H H P ,T
H H P ,T

+ w
= R sTc

R sTc
R sTc

=
6176.21 Btu/Lbmol

Hp , t = Hm H H P ,T

(0 )

H H P ,T

H H P ,T

Hp,t =

(I )

xiPci
159.6
242.6
164.5
24.4
38.1
7.1
8.5
0.9
0.3
646.0

xiwi
0.00258
0.03331
0.04052
0.00857
0.01376
0.00330
0.00440
0.00063
0.00032
0.10738

-2207.7 Btu/Lbmol

H Mezcla = H V V + H L L
Hmezcla =

2383.3

Btu/Lbmol

Calculo de la entropa a la salida del Turbo-expansor (-60 F)

Fase gaseosa
Temp =
Presin =

-60
400

Comp
C1
C2
C3
iC4
nC4
iC5
nC5
nC6
C7+

yi
0.8854
0.0977
0.0147
0.0011
0.0010
0.0001
0.0000
0.0000
0.0000

F
Lpca
PMi
16.043
30.070
44.097
58.123
58.123
72.150
72.150
86.177
122.00

si
2.955
1.889
1.482
1.251
1.193
1.122
1.106
1.065
1.036

yiPMisi
41.98
5.55
0.96
0.08
0.07
0.00
0.00
0.00
0.00
48.6

yilnyi
-0.108
-0.227
-0.062
-0.007
-0.007
-0.001
0.000
0.000
0.000
-0.413

yiTci
304.0
53.7
9.8
0.8
0.8
0.1
0.0
0.0
0.0
369.2

Preparado por Ing. Edwards Nava

yiPci
590.5
69.1
9.0
0.6
0.6
0.0
0.0
0.0
0.0
669.9

yiwi
0.00956
0.00949
0.00223
0.00020
0.00021
0.00001
0.00001
0.00000
0.00000
0.02172

14/06/2007

sTr =
sPr =

S S P ,T

1.0834
0.5971

S S P ,T

(0)

0.38

Grafico
Tc.ln(PM)/ln(Tbn)

(I )

0.32

Sm = Yi Si PMi R Yi Ln (Yi )
Sm =

S S P ,T

S S P ,T

Sp,t =

49.46

37.138
14.984
10.591
8.573
7.846
7.088
6.539

xiPci
159.6
242.6
164.5
24.4
38.1
7.1
8.5
0.9
0.3
646.0

xiwi
0.00258
0.03331
0.04052
0.00857
0.01376
0.00330
0.00440
0.00063
0.00032
0.10738

Btu/Lbmol R

S S P ,T ( 0 )
S S P ,T

+ w
= R
R
R

= 7.33 Btu/Lbmol R

42.134

Ln(PM).ln(1/w).S

180
322
418
485
502
564
572

(I )

+ Ln ( P / 14 .7 )

Btu/Lbmol R

Fase Liquida
Temp =
Presin =

-60
400

Comp
C1
C2
C3
iC4
nC4
iC5
nC5
nC6
C7+

xi
0.2393
0.3427
0.2675
0.0463
0.0694
0.0144
0.0175
0.0021
0.0008

sTr =
sPr =

F
Lpca
PMi
16.043
30.070
44.097
58.123
58.123
72.150
72.150
86.177
122.00

#REF!
0.6192

si
2.955
1.889
1.482
1.251
1.193
1.122
1.106
1.065
1.036

xiPMisi
11.34
19.46
17.48
3.36
4.81
1.17
1.40
0.19
0.10
59.3

S S P ,T

(0)

S S P ,T

(I )

xilnxi
-0.342
-0.367
-0.353
-0.142
-0.185
-0.061
-0.071
-0.013
-0.006
-1.540

4.55

5.30

xiTci
#REF!
188.5
178.1
34.0
53.2
12.0
14.8
1.9
0.8
#REF!

Sm = Yi Si PMi R Yi Ln (Yi )
Sm =

S S P ,T

S S P ,T

Sp,t =

62.38

Btu/Lbmol R

(I )
S S P ,T ( 0 )

S S P ,T

+ w
+ Ln ( P / 14 .7 )
= R

R
R

= 16.73 Btu/Lbmol R

45.656

Btu/Lbmol R

S Mezcla = S V V + S L L
Smezcla =

42.414

Btu/Lbmol R

Preparado por Ing. Edwards Nava

14/06/2007

Con las entropas obtenidas a las distintas temperaturas sern utilizadas para graficar
T vs S donde entraremos con la entropa de entrada (S 2 = S 1 ) y obtendremos la
temperatura de salida en el turbo-expansor (isoentropico)
T
F
-20
-40
-60

Stotal
Btu/LbmolF
44.168
43.266
42.414

Htotal
Btu/Lbmol
3153.1
2767.3
2383.3

De la grafica antes mencionada se obtiene el siguiente ajuste

S = 6.237305E- 05 T 2 + 4.883568E- 02 T + 4.511996E+ 01


donde para una S = 43.374 Btu/LbmolF
se obtiene una Tsalida = -37.55 F del turbo expansor; ideal a entropia constante.
Con la temperatura obtenida entramos al grafico T vs H donde obtendremos H 2
H2 =
2814.451 Btu/Lbmol
H1 =
3165.79 Btu/Lbmol

h2

real

= ( h2 h1 ) is + h1

H2 real =

2877.692 Btu/Lbmol

Calculo del trabajo Maximo

W
W

max

max

= m (h2 h1 )
-7212

hp

Calculo del trabajo Real


W real = m (h2
W

real

real

-5914

h1 )
hp

Calculo de la temperatura real a la salida del Turbo-expansor


Con las entalpias obtenidas a las distintas temperaturas sern utilizadas para graficar
T vs H donde entraremos con la entalpia de salida (H 2 real) y obtendremos la
temperatura de salida en el turbo-expansor (real)

H = 2.254074E- 03 T 2 + 1.942341E+ 01 T + 3.540622E+ 03


donde para una H = 2877.69 Btu/Lbmol
se obtiene una Tsalida = -34.26 F del turbo expansor

Preparado por Ing. Edwards Nava

14/06/2007

Graficas para Determinar KC7+ , SC6 y SC7+


y HC6 a la entrada del Turbo-expansor

Entrada Turbo-expansor
T = 20F P = 800 Lpca
2

Ki vs Tc
10.0000

Ki

1.0000

0.1000

y = 6.262348E+00e-7.766982E-06x
R2 = 9.982822E-01
0.0100

0.0010
0

100000

200000

300000

400000

500000

Tc2

600000

700000

800000

900000

Entrada Turbo-expansor
T = 20F P = 800 Lpca
Tci.Ln(PMi)/Ln(Tbni) vs Si.Ln(PMi).ln(1/wi)

Si.Ln(PMi).ln(1/wi))

100

10

y = 7.028338E+04x-1.451188E+00
R2 = 9.978140E-01

1
100

1000
Tci.Ln(PMi)/Ln(Tbni)

Entrada Turbo-expansor
T = 20F P = 800 Lpca
PMi.ewi vs Tci.Hi0.5/PMi
1000

100

Tci.Hi0.5/PMi

y = 1.534315E+03x-5.551037E-01
R2 = 9.950624E-01

10

1
10

100

PMi.ewi

1000

Graficas para Determinar KC7+ , SC6 y SC7+


y HC6 a la salida del Turbo-expansor
a temperatura de -20F

Graficas para Determinar KC7+ , SC6 y SC7+


y HC6 a la salida del Turbo-expansor
a temperatura de -40F

Graficas para Determinar KC7+ , SC6 y SC7+


y HC6 a la salida del Turbo-expansor
a temperatura de -60F

Salida Turbo-expansor

44.4
44.2
44
43.8
43.6

43.4
43.2
43
42.8

y = 6.237305E-05x2 + 4.883568E-02x + 4.511996E+01


R2 = 1.000000E+00

42.6
42.4
42.2
-70

-60

-50

-40

-30
Tsalida, F

-20

-10

Salida Turbo-expansor

3500

3000 H2 real

2500

H2 ideal

2000

1500

y = 2.254074E-03x2 + 1.942341E+01x + 3.540622E+03


R2 = 1.000000E+00

1000

500
Tideal

Treal

0
-70

-60

-50

-40

-30
Tsalida, F

-20

-10

También podría gustarte