DESTILACION FLASH

1.- CALCULO DE TEMPERATURA EN EL PUNTO DE ROCIO:

Pantalla de cálculo de propiedades de la corriente en el punto de rocío,

determinado por el simulador de ASPEN HYSYS

Composición
2.- ESPECIFICAR LAS CORRIENTES DE SALIDA DE UN FLASH Y MOSTRAR
RESULTADOS DE TODO EL PROCESO DE ITERACIONES:

Algoritmo de resolución del flash

Conocidos presión y temperatura ( P y T )

2.1.- calcular presiones de saturación para cada componente.

2.2.- resolver el algoritmo ideal para el proceso

2.3.- calcular Ki , real

2.4.- suponer β

Zi Ki Zi
2.5.- calcular y i= ; x i=
(1+ β )(K i−1) (1+ β )( K i−1)

2.6.- verificar y i , x i

- SI→ FIN

- NO → NORMALIAR

xi yi
x i= ; xi = ;
∑ xi ∑ yi
x i ,k ≠ x i , k+1
y i , k ≠ y i ,k +1

- Si cumple volver a paso 3

- No cumple suponer otra fracción vaporizada

Relación de equilibrio

yi
K i=
xi

Criterio para resolver el flash

N
Z (1−K )
∑ 1+iβ (K −1)
i
=0
i=1 i

Datos conocidos:

Temperatura es igual al intermedio entre en punto de Burbuja y puto de Rocío:

Temperaturas calculadas en simulador ASPEN HYSYS

 Calculo de temperatura en el p8unto de burbuja P y Xi conocidos.

 Calculo de temperatura en el punto de roció P y Yi conocidos.

  Calculo de temperatura flash

T(burbuja) ◦C T(rocío) ◦C T(flash)◦C T(flash) K P(kpa)

-199,7 0,7879 -99,45605 173,69395 5230,888

T burbuja +T rocio
T flash=
2

−199,7 ( ℃ )+ 0,7879(℃ )
T flash=
2

T flash=−99,45605 ℃=173,6939 K

Composición

Componente x %xi
N2 0,160 0,0016
CO2 6,600 0,066
C1 80,910 0,8091
C2 7,610 0,0761
C3 3,020 0,0302
iC4 0,330 0,0033
nC4 0,780 0,0078
iC5 0,210 0,0021
nC5 0,200 0,002
C6(+) 0,080 0,0008
C7(+) 0,100 0,001
100,000 1
2.1.- Calcular presiones de saturación para cada componente.

Constates para cálculo de presiones de saturación:

Componente Tc[K] Pc[bar] A B C D ECUACION

N2 126,2 33,9 -6,09676 1,1367 -1,04072 -1,93306 3
CO2 304,1 73,8 -6,95626 1,19695 -3,12614 2,99448 1
C1 190,4 46 -6,00435 1,1885 -0,83408 -1,22833 3
C2 305,4 48,8 -6,34307 1,0163 -1,19116 -2,03539 1
C3 369,8 42,5 -6,72219 1,33236 -2,13868 -1,38551 1
iC4 408,2 36,5 -6,95579 1,5009 -2,52717 -1,49776 1
nC4 425,2 38 -6,88709 1,15157 -1,99873 -3,13003 1
iC5 460,4 33,8 -7,56707 2,15002 -3,8903 -3,70363 1
nC5 469,7 33,7 -7,28936 1,53679 -3,08367 -1,02456 1
C6(+) 507,5 30,1 -7,46765 1,44211 -3,28222 -2,50941 1
C7(+) 540,3 27,4 -7,67468 1,37068 -3,5362 -3,20243 1
Calculo de presiones de saturación con las respectivas ecuaciones de vapor

EC. PVAP 3:

B
ln Ps =A−
( T +C)

EC. PVAP 1:

Ps
ln ( )
Pc
−1 1,5 3 6
=(1−X ) ∗[ A∗X + B∗X +C∗X + D∗X ]

T
X =1−
Tc

Componente X ln(Ps/Pc) Psat,(bar) Psat,(KPA)

N2 -0,376338748 0 0,002235381 0,22353807
CO2 0,428826208 -5,03314873 0,481047145 48,1047145
C1 0,087741859 0 0,002451083 0,24510825
C2 0,431257531 -4,49466257 0,545020291 54,5020291
C3 0,530303002 -7,23877773 0,030523034 3,05230344
iC4 0,574488119 -9,10779885 0,004044139 0,40441391
 nC4 0,591500588 -10,0307109 0,00167302 0,16730203
iC5 0,622732515 -12,7526785 9,7836E-05 0,0097836
nC5 0,630202363 -12,6039819 0,000113185 0,01131853
C6(+) 0,657745911 -15,4263242 6,01172E-06 0,00060117
C7(+) 0,678523135 -18,2238008 3,33622E-07 3,3362E-05

2.2.- Resolver el algoritmo ideal para el proceso

Primera iteración ideal

Psaturacion
K ideal =
P

N
Z i( 1−K i )
Criterio para resolver el flash f ( FVAP ) =∑ =0
i=1 1+ β ( K i −1)

Z K Z
y i= ¿i ¿ i ; x i= ¿ ¿i

Verificación ∑ ( xi − y i )=0

FVAP=0,32

Nro. de VERIFICAIO
Ki f(β) yi xi Yi,norm xi,norm
comp. N
4,273E- 0,00235 4,72E- 1,60E-
1 1,005E-07 0,002 1,467454
05 3 05 03
9,196E- 0,09575 4,17E- 6,58E-
2 8,887E-04 0,097
03 2 01 02
4,686E- 1,18977 2,62E- 8,10E-
3 5,575E-05 1,190
05 1 02 01
1,042E- 0,11020 5,44E- 7,58E-
4 1,160E-03 0,111
02 5 01 02
5,835E- 0,04437 1,22E- 3,02E-
5 2,591E-05 0,044 0,00E+00
04 4 02 02
7,731E- 0,00485 1,76E- 3,30E-
6 3,752E-07 0,005 FIN
05 2 04 03
3,198E- 0,01147 1,72E- 7,81E-
7 3,669E-07 0,011
05 0 04 03
8 1,870E- 0,00308 5,776E-09 0,003 2,71E- 2,10E-
 06 8 06 03
2,164E- 0,00294 2,99E- 2,00E-
9 6,364E-09 0,003
06 1 06 03
1,149E- 0,00117 6,34E- 8,01E-
10 1,352E-10 0,001
07 6 08 04
6,378E- 0,00147 9,37929E- 4,40E- 1,00E-
11 0,001
09 1 12 09 03
1,4695851 1,00E+0 1,00E+0
1,467 0,0021
3 0 0
f(β)

Verificación ∑ ( xi − y i )=0

No cumple

Segunda iteración ideal:

FVAP=0,80

Nro. de VERIFICAI
Ki f(β) yi xi Yi,norm xi,norm
comp. ON
4,273E- 0,00200 8,547E- 4,71E- 1,60E-
1 0,002 0, 0002
05 0 08 05 03
9,196E- 0,08155 7,570E- 4,17E- 6,59E-
2 0,082
03 4 04 01 02
4,686E- 1,01131 4,739E- 2,61E- 8,09E-
3 1,011
05 6 05 02 01
1,042E- 0,09388 9,886E- 5,44E- 7,59E-
4 0,095
02 9 04 01 02
5,835E- 0,03772 2,202E- 1,21E- 3,02E-
5 0,038 0,00E+00
04 2 05 02 02
7,731E- 0,00412 3,189E- 1,76E- 3,30E-
6 0,004 FIN
05 5 07 04 03
3,198E- 0,00975 3,118E- 1,72E- 7,80E-
7 0,010
05 0 07 04 03
1,870E- 0,00262 4,910E- 2,70E- 2,10E-
8 0,003
06 5 09 06 03
2,164E- 0,00250 5,409E- 2,98E- 2,00E-
9 0,002
06 0 09 06 03
10 1,149E- 0,00100 1,149E- 0,001 6,33E- 8,00E-
 07 0 10 08 04
6,378E- 0,00125 7,9724E- 4,39E- 1,00E-
11 0,001
09 0 12 09 03
1,0249546 1,00E+0 1,00E+0
1,248 1,0018
09 0 0
f(β)

Verificación ∑ ( xi − y i )=0

CUMPLE DONDE ∑ ( xi − y i )=0,0002

2.3.- Calcular Ki , real

 Ejemplo de cálculo a detalle de componente 1 en ultima iteración para

factor de vaporización (B) = 0,81

CÁLCULO DE COEFICIENTE DE FUGACIDAD DE SATURACIÓN PARA CADA

COMPONENTE USANDO ECUACIÓN DE ESTADO SOAVE-R-W

Calculo de volumen específico para sustancia pura con Ecuación de Soave-R-W;

Tabla de constantes y propiedades críticas de cada componente del gas:

Componente Xi Tc[K] Pc[bar] Pc[Kpa] w R[Kpa*m3/Kmol*K]

N2 0,0016 126,2 33,9 3390 0,039
CO2 0,0659 304,1 73,8 7380 0,239
C1 0,8094 190,4 46 4600 0,011
C2 0,0759 305,4 48,8 4880 0,099
C3 0,0302 369,8 42,5 4250 0,153
iC4 0,0033 408,2 36,5 3650 0,183 8,31
nC4 0,0078 425,2 38 3800 0,199
iC5 0,0021 460,4 33,8 3380 0,227
nC5 0,0020 469,7 33,7 3370 0,251
C6(+) 0,0008 507,5 30,1 3010 0,299
C7(+) 0,0010 540,3 27,4 2740 0,349
 RT aα
P= −
V −b V (V +b)

RT 2 aα −bRT −P b2 aαb
V 3− ( ) (
P
V +
P
V−
P) =0

Dónde:

R 2 T 2c
a=0,42748
Pc

RTc
b=0,0867
Pc

s=0,4858+1,55171 w−0,15613 w2

α =¿ ¿

Cal. de volumen especifico con ec. de est.soave-r.k

a B s α V[m3/Kmol]
138,6873532 0,026821297 0,546078826 0,8198089 0,0013
369,9085135 0,029687948 0,847737998 1,45697256 0,087
232,6450973 0,029821483 0,502849808 1,045643 0,316
564,2039102 0,045088868 0,63788807 1,33824334 0,386
949,8668514 0,062689975 0,719555253 1,50408655 0,006
1347,630975 0,080574989 0,764532462 1,60229575 0,0755
1404,49683 0,080617584 0,788405396 1,64995129 0,0802
1851,278587 0,098138346 0,829990677 1,7429094 0,4023
1932,542563 0,100417818 0,865440354 1,79334285 0,4111
2525,942417 0,121475773 0,935800122 1,9274707 0,0118
3145,120002 0,142070702 1,00832651 2,06386602 0,10021

Calcular factor Z, como sustancia pura para cada componente:

Z3 −Z 2 + ( A−B−B2 ) Z− AB=0

Dónde:

aαP
A=
( RT )2
 bP
B=
RT

Ecuación de Estado Soave-R-W para sustancia pura

A B
ln ∅si =Z −1−ln ( Z−B )−
B ( )
ln 1+
Z

A B Z ln∅s ∅s
0,1253184
1,2199E-05 0,8 0,19350 1,2135
3
0,0124440 29,850491
0,68 0,00000 1,0000
5 7
0,1527815
2,862E-05 0,692 0,30960 1,3629
2
0,0197520 51,364808
0,9 0,00000 1,0000
6 3
0,0020931 3,9995347
0,9 0,01000 1,0101
1 3
0,0004191 0,6810986
0,72 2,96632 0,1290
5 3
0,0001860 0,2819127
0,6 0,74518 2,1068
9 1
0,0200687
1,5152E-05 0,8 0,04853 1,0497
9
0,0237566
1,8828E-05 0,9 0,03209 1,0326
2
0,0015264
1,4049E-06 0,9 0,00706 1,0071
1
1,0394E-07 9,907E-05 0,71 0,05263 1,0540

Cálculo de coeficiente de fugacidad parcial molar como mezcla:

CALCULO DE COEFICIENTE DE FUGACIDAD PARCIAL DE VAPOR COMO
MEZCLA (∅ V ¿, CON ECUACIÓN DE ESTADO SOAVE-R-W:

Calculo de factor Z de la mezcla:

Z3 −Z 2 + ( A−B−B2 ) Z− AB=0

Ecuación de Estado Soave para mezclas

Bi A Bi 2 B
ln Ǿ vi =
B [
( z−1 ) −ln ( z−B ) + ]( )
− ∑ yi ( aα ) ij ln 1+
B B aα z

R2 T 2c
ai=0,42748
Pc

R Tc
bi=0,0867
Pc

si=0,4858+1,55171 w−0,15613 w2

αi=¿ ¿

Componen Tc[K Pc[Kp R[Kpa*m3/Kmol*

xi W ai bi si α1 (aα)i
te ] a] K]
0,001 126, 0,03 8,31 0,0 0,5 0,8
N2 3390 138,69 113,7
6 2 9 3 5 2
0,065 304, 0,23 0,0 0,8 1,4
CO2 7380 369,91 538,9
9 1 9 3 5 6
0,809 190, 0,01 0,0 0,5 1,0
C1 4600 232,65 243,3
4 4 1 3 0 5
0,075 305, 0,09 0,0 0,6 1,3
C2 4880 564,20 755,0
9 4 9 5 4 4
0,030 369, 0,15 0,0 0,7 1,5 1428,
C3 4250 949,87
2 8 3 6 2 0 7
0,003 408, 0,18 1347,6 0,0 0,7 1,6 2159,
iC4 3650
3 2 3 3 8 6 0 3
0,007 425, 0,19 1404,5 0,0 0,7 1,6 2317,
nC4 3800
8 2 9 0 8 9 5 4
iC5 0,002 460, 3380 0,22 1851,2 0,1 0,8 1,7 3226,
1 4 7 8 0 3 4 6
 0,002 469, 0,25 1932,5 0,1 0,8 1,7 3465,
nC5 3370
0 7 1 4 0 7 9 7
0,000 507, 0,29 2525,9 0,1 0,9 1,9 4868,
C6(+) 3010
8 5 9 4 2 4 3 7
0,001 540, 0,34 3145,1 0,1 1,0 2,0 6491,
C7(+) 2740
0 3 9 2 4 1 6 1

( aα ) ij=( 1−kij )∗√( aα ) i ( aα ) j

kii=0

Desarrollo de ecuación para componente (1)

( aα )11=( 1−k 11 )∗√ ( aα )1 ( aα )1=¿

( aα )12=( 1−k 12)∗√ ( aα )1 ( aα )2 =¿

( aα )13=( 1−k 13 )∗√ ( aα )1 ( aα )3=¿

( aα )14=( 1−k 14 )∗√ ( aα )1 ( aα )4 =¿

( aα )15=( 1−k 15 )∗√ ( aα )1 ( aα )5=¿

( aα )16=( 1−k 16 )∗ √ ( aα )1 ( aα )6=¿

( aα )17=( 1−k 17 )∗√ ( aα )1 ( aα )7=¿

( aα )18=( 1−k 18 )∗√ ( aα )1 ( aα )8=¿

( aα )19=( 1−k 19 )∗√ ( aα )1 ( aα )9=¿

( aα )110=( 1−k 110 )∗√ ( aα )1 ( aα )10 =¿

( aα )111=( 1−k 111)∗√ ( aα )1 ( aα )11=¿

(aα)ij
k
ij 1 2 3 4 5 6 7 8 9 10 11
 1,1E+ 2,5E+ 1,7E+ 2,9E+ 4,0E+ 5,0E+ 5,1E+ 6,1E+ 6,3E+ 7,4E+ 8,6E+
0 02 02 02 02 02 02 02 02 02 02 02
2,5E+ 5,4E+ 3,6E+ 6,4E+ 8,8E+ 1,1E+ 1,1E+ 1,3E+ 1,4E+ 1,6E+ 1,9E+
0 02 02 02 02 02 03 03 03 03 03 03
1,7E+ 3,6E+ 2,4E+ 4,3E+ 5,9E+ 7,2E+ 7,5E+ 8,9E+ 9,2E+ 1,1E+ 1,3E+
0 02 02 02 02 02 02 02 02 02 03 03
2,9E+ 6,4E+ 4,3E+ 7,6E+ 1,0E+ 1,3E+ 1,3E+ 1,6E+ 1,6E+ 1,9E+ 2,2E+
0 02 02 02 02 03 03 03 03 03 03 03
4,0E+ 8,8E+ 5,9E+ 1,0E+ 1,4E+ 1,8E+ 1,8E+ 2,1E+ 2,2E+ 2,6E+ 3,0E+
0 02 02 02 03 03 03 03 03 03 03 03
5,0E+ 1,1E+ 7,2E+ 1,3E+ 1,8E+ 2,2E+ 2,2E+ 2,6E+ 2,7E+ 3,2E+ 3,7E+
0 02 03 02 03 03 03 03 03 03 03 03
5,1E+ 1,1E+ 7,5E+ 1,3E+ 1,8E+ 2,2E+ 2,3E+ 2,7E+ 2,8E+ 3,4E+ 3,9E+
0 02 03 02 03 03 03 03 03 03 03 03
6,1E+ 1,3E+ 8,9E+ 1,6E+ 2,1E+ 2,6E+ 2,7E+ 3,2E+ 3,3E+ 4,0E+ 4,6E+
0 02 03 02 03 03 03 03 03 03 03 03
6,3E+ 1,4E+ 9,2E+ 1,6E+ 2,2E+ 2,7E+ 2,8E+ 3,3E+ 3,5E+ 4,1E+ 4,7E+
0 02 03 02 03 03 03 03 03 03 03 03
7,4E+ 1,6E+ 1,1E+ 1,9E+ 2,6E+ 3,2E+ 3,4E+ 4,0E+ 4,1E+ 4,9E+ 5,6E+
0 02 03 03 03 03 03 03 03 03 03 03
8,6E+ 1,9E+ 1,3E+ 2,2E+ 3,0E+ 3,7E+ 3,9E+ 4,6E+ 4,7E+ 5,6E+ 6,5E+
0 02 03 03 03 03 03 03 03 03 03 03

R2 T 2c
a i=0,42748
Pc

aα =∑ ¿ ∑ yiyj ( aα ) ij

aα = y 21 ( aα )11 + y 22 ( aα )22+ y 23 ( aα )33+ y 24 ( aα )44 + y 25 ( aα )55 + y 26 ( aα )66 + y 27 ( aα )77 + y 28 ( aα )88+ y 29 ( aα )99+ y 210 ( aα )1010 + y 211 ( a

Aα

y^2i*(aα)ii yi*yj*(aα)i9
1 2 3 4 5 6 7 8 9 10 11 aα
0,00 0,03 0,22 0,04 0,02 0,00 3,29 0,00 0,00 0,00 0,00
14648,78841

0,03 2,34 19,31 3,19 1,75 0,23 641,94 0,18 0,18 0,09 0,12
0,22 19,31 159,36 26,34 14,41 1,94 3560,19 1,51 1,49 0,70 1,02
167,56
0,04 3,19 26,34 4,35 2,38 0,32 1036,63 0,25 0,25 0,12 0,17
0,02 1,75 14,41 2,38 1,30 0,18 780,33 0,14 0,13 0,06 0,09
0,00 0,23 1,94 0,32 0,18 0,02 128,89 0,02 0,02 0,01 0,01
 0,01 0,57 4,74 0,78 0,43 0,06 326,95 0,04 0,04 0,02 0,03
0,00 0,18 1,51 0,25 0,14 0,02 122,56 0,01 0,01 0,01 0,01
0,00 0,18 1,49 0,25 0,13 0,02 125,38 0,01 0,01 0,01 0,01
0,00 0,09 0,70 0,12 0,06 0,01 70,45 0,01 0,01 0,00 0,00
0,00 0,12 1,02 0,17 0,09 0,01 117,41 0,01 0,01 0,00 0,01
0,31 27,99 231,03 38,18 20,89 2,81 6914,03 2,18 2,16 1,02 0,01

R Tc
bi=0,0867
Pc

b=∑ yibi

b= y 1 b1+ y 2 b 2+ y 3 b3 + y 4 b4 + y 5 b5 + y 6 b6 + y 7 b7 + y 8 b 8+ y 9 b 9+ y 10 b10 + y 11 b 11

( aα ) i∗P
Ai=
(RT )2

Component
b Ai
e
4,29292E- 0,28546536
N2
05 1
1,35316298
CO2 0,00195562
8
0,02413704 0,61077511
C1
4 2
0,00342359 1,89572507
C2
1 3
0,00189364 3,58707031
C3
9 8
0,00026598 5,42148290
iC4
9 7
0,00062904 5,81830337
nC4
1 2
0,00020616 8,10123363
iC5
5 5
0,00020090 8,70155807
nC5
8 2
C6(+) 9,72159E- 12,2240747
 05 6
0,00014212
C7(+) 16,2975989
2
0,03299427
5

Aij =(1−k ij )( A i A j )0,5

k ii=0

Desarrollo de ecuación para componente (1)

A11 =(1−k 11 )( A1 A 1 )0,5

A12=(1−k 12)(A 1 A 2)0,5

A13=(1−k 13)( A1 A 3 )0,5

A14 =(1−k 14)( A 1 A 4 )0,5

A15=(1−k 15)( A1 A 5 )0,5

A16=(1−k 16)( A 1 A6 )0,5

A17=(1−k 17)( A 1 A7 )0,5

A18=(1−k 18)( A1 A 8 )0,5

A19=(1−k 19)( A1 A 9 )0,5

A110 =(1−k 110 )( A1 A 10)0,5

A111=(1−k 111)( A1 A 11 )0,5

Aij
1 2 3 4 5 6 7 8 9 10 11
 0,285 0,622 0,418 0,736 1,012 1,244 1,289 1,521 1,576 1,868 2,157
0,622 1,353 0,909 1,602 2,203 2,709 2,806 3,311 3,431 4,067 4,696
0,418 0,909 0,611 1,076 1,480 1,820 1,885 2,224 2,305 2,732 3,155
0,736 1,602 1,076 1,896 2,608 3,206 3,321 3,919 4,061 4,814 5,558
1,012 2,203 1,480 2,608 3,587 4,410 4,568 5,391 5,587 6,622 7,646
1,244 2,709 1,820 3,206 4,410 5,421 5,616 6,627 6,868 8,141 9,400
1,289 2,806 1,885 3,321 4,568 5,616 5,818 6,866 7,115 8,433 9,738
1,521 3,311 2,224 3,919 5,391 6,627 6,866 8,101 8,396 9,951 11,490
1,576 3,431 2,305 4,061 5,587 6,868 7,115 8,396 8,702 10,314 11,909
1,868 4,067 2,732 4,814 6,622 8,141 8,433 9,951 10,314 12,224 14,115
2,157 4,696 3,155 5,558 7,646 9,400 9,738 11,490 11,909 14,115 16,298

A=∑ ¿ ∑ yiyjAij

A= y 21 A11 + y 22 A22 + y 23 A 33+ y24 A 44 + y 25 A 55+ y 26 A66 + y 27 A 77 + y 28 A 88+ y29 A 99+ y 210 A 1010 + y 211 A111 +2[ y 1 y 2 A 12+ y 1 y 3

yiyjAij
y^2i*A
ii 1 2 3 4 5 6 7 8 9 10 11 A
7E- 7E- 5E- 9E- 5E- 7E- 2E- 5E- 5E- 2E- 3,45358E-
07 05 04 05 05 06 05 06 06 06 06
7E- 6E- 5E- 8E- 4E- 6E- 1E- 5E- 5E- 2E- 0,00030945
05 03 02 03 03 04 03 04 04 04 6
5E- 5E- 4E- 7E- 4E- 5E- 1E- 4E- 4E- 2E- 0,00255455
04 02 01 02 02 03 02 03 03 03 4
9E- 8E- 7E- 1E- 6E- 8E- 2E- 6E- 6E- 3E- 0,00042220
05 03 02 02 03 04 03 04 04 04 2
5E- 4E- 4E- 6E- 3E- 4E- 1E- 3E- 3E- 2E- 0,00023104
05 03 02 03 03 04 03 04 04 04 2 2,095152091
0,420 7E- 6E- 5E- 8E- 4E- 6E- 1E- 5E- 5E- 2E- 3,10414E-
06 04 03 04 04 05 04 05 05 05 05
2E- 1E- 1E- 2E- 1E- 1E- 4E- 1E- 1E- 5E- 7,60093E-
05 03 02 03 03 04 04 04 04 05 05
5E- 5E- 4E- 6E- 3E- 5E- 1E- 4E- 4E- 2E- 2,41475E-
06 04 03 04 04 05 04 05 05 05 05
5E- 5E- 4E- 6E- 3E- 5E- 1E- 4E- 3E- 2E- 2,38345E-
06 04 03 04 04 05 04 05 05 05 05
2E- 2E- 2E- 3E- 2E- 2E- 5E- 2E- 2E- 8E- 1,12999E-
06 04 03 04 04 05 05 05 05 06 05
3E- 3E- 3E- 4E- 2E- 3E- 8E- 2E- 2E- 1E- 1,63094E-
 06 04 03 04 04 05 05 05 05 05 05
8E- 7E- 6E- 1E- 5E- 7E- 2E- 5E- 5E- 3E- 1,63094E-
04 02 01 01 02 03 02 03 03 03 05

biP
Bi=
RT

B=∑ yi Bi

B= y 1 B 1+ y 2 B 2+ y 3 B3 + y 4 B4 + y 5 B 5+ y 6 B6 + y 7 B7 + y 8 B8 + y 9 B 9+ y10 B 10 + y 11 B11

Component
Bi B
e
0,00015557
N2 2,49009E-07
6
0,00708719 0,00046685
CO2
2 1
0,08747295 0,07079924
C1
3 6
0,01240713 0,00094207
C2
9 2
0,00686260 0,00020729
C3
7 6
0,00096394
iC4 3,18212E-06
7
0,00227965
nC4 1,77876E-05
1
0,00074714
iC5 1,56958E-06
6
0,00072809
nC5 1,45672E-06
5
0,00035231
C6(+) 2,81952E-07
2
0,00051505
C7(+) 5,1524E-07
3
 0,07244050
8

Calculo de factor Z de la mezcla:

Z3 −Z 2 + ( A−B−B2 ) Z− AB=0

Zmez cla=0,95

Ecuación de Estado Soave para mezclas

Bi A Bi 2 B
ln Ǿ vi =
B [
( z−1 ) −ln ( z−B ) + ]( )
− ∑ yi ( aα ) ij ln 1+
B B aα z

Componente Bi B lnǾ^vi Ǿ^vi

N2 0,000155576 2,49009E-07 0,190263719 1,210
CO2 0,007087192 0,000466851 0,0122 1,012
C1 0,087472953 0,070799246 0,13337 1,143
C2 0,012407139 0,000942072 0,001 1,001
C3 0,006862607 0,000207296 2,735014164 0,990
iC4 0,000963947 3,18212E-06 0,49364757 1,638
nC4 0,002279651 1,77876E-05 0,987389755 2,684
iC5 0,000747146 1,56958E-06 0,411937134 1,510
nC5 0,000728095 1,45672E-06 0,404548281 1,499
C6(+) 0,000352312 2,81952E-07 0,263327363 1,301
C7(+) 0,000515053 5,1524E-07 0,335893584 1,399

Calculo de FP

Datos necesarios:

component
xi Psat,(KPA)
e
0,22353807
N2 0,000
3
48,1047144
CO2 0,417
8
0,24510825
C1 0,026
1
 54,5020291
C2 0,544
4
3,05230344
C3 0,012
5
0,40441390
iC4 0,000
6
0,16730202
nC4 0,000
7
0,00978360
iC5 0,000
2
0,01131853
nC5 0,000
2
0,00060117
C6(+) 0,000
2
3,33622E-
C7(+) 0,000
05

POR METODO DE SIMPSON Y TRAPECIO

Psist
1
FP i=exp ∫ V liq dP
RT Psat

P−Psat
∆ P=
n

Dónde:

n=numero de intervalos

Integración numérica por método de Simpson:

∫ f ( x ) dx= h3 [ ( y 0 + y n ) + 4 ( y 1+ y 3+ y n−1) + 2 ( y 2 + y 4+ y n−2 ) ]

a

N2
V(m3/kmo
ΔP P(kpa) I 1/RT FP
l)
 0,22353807
0,0012
3
523,289984
0,0213
3
1046,35643 0,00167
1569,42287
0,01666
7
2092,48932
0,0021
3
523,06644 2615,55576 274,365786 0,0006928
0,00104 1,209350478
6 9 6 1
3138,62221
0,3021
5
3661,68866
0,0299
1
4184,75510
0,00843
8
4707,82155
0,0088
4
5230,888 0,031

DETERMINACIÓN DE COEFICIENTE DE ACTIVIDAD POR CORRELACIÓN DE

UNIFAC:

ln γ i=ln γ ic + ln γ Gi

C φi Z θi φi N
ln γ =ln + q i ln −l i − ∑ x i l i
i
xi 2 φi x i j=1

Z=35,2−0,1272 T ( K ) +0,00014 T 2 (K )

Z
l i= ( r i−q i )−( r i−1 )
2
Desarrollo para componente (1) y (2)

Z
l 1= ( r 1−q 1) −( r 1−1 )
2

10,0000
l 1= ( 1,06−3,882 )−( 1,06−1 )
2

l 1=−1,42E+01

Z
l 2= ( r 2−q 2 )−( r 2 −1 )
2

10,0000
l 2= ( 1,0021−3,002 )− (1,0021−1 )
2

l 2=−1,00E+01

Componente Z ri Qi li
N2 1,06 3,882 -1,42E+01
CO2 1,0021 3,002 -1,00E+01
C1 1,3454 3,1176 -9,21E+00
C2 1,1167 1,0252 3,41E-01
C3 1,923 1,023 3,58E+00
iC4 10,0000 1,542 0,892 2,71E+00
nC4 1,7233 1,002 2,88E+00
iC5 1,022 1,003 7,30E-02
nC5 1,0023 0,832 8,49E-01
C6(+) 1,182 0,998 7,38E-01
C7(+) 0,0909 0,0823 9,52E-01

x k qi
θi = N

∑ xjqj
j=1

Desarrollo para componente (1) y (2)

 x 1 q1
θ1 =
x 1 q 1+ x2 q 2+ x 3 q3 + x 4 q 4 + x 5 q5 + x 6 q 6 + x 7 q7 + x 8 q 8 + x 9 q9 + x 10 q 10+ x 11 q 11

xk ri
φ i= N

∑ x jr j
j=1

Desarrollo para componente (1) y (2)

x1 r1
φ 1=
x 1 r 1+ x 2 r 2+ x3 r 3 + x 4 r 4 + x 5 r 5 + x 6 r 6 + x 7 r 7 + x 8 r 8 + x 9 r 9 + x 10 r 10 + x 11 r 11

φ 1=(4,99E-05)/(4,99E-05+ 4,18E-01+3,51E-02+6,08E-01+2,33E-02+2,71E-04+2,96E-04 +2,76E-06+2,99E

φ 1=4,60E-05

Componente xi*qi θi xi*ri φ1

N2 1,83E-04 9,60E-05 4,99E-05 4,60E-05
CO2 1,25E+00 6,57E-01 4,18E-01 3,85E-01
C1 8,14E-02 4,27E-02 3,51E-02 3,24E-02
C2 5,58E-01 2,93E-01 6,08E-01 5,60E-01
C3 1,24E-02 6,52E-03 2,33E-02 2,15E-02
iC4 1,57E-04 8,23E-05 2,71E-04 2,50E-04
nC4 1,72E-04 9,04E-05 2,96E-04 2,73E-04
iC5 2,71E-06 1,42E-06 2,76E-06 2,55E-06
nC5 2,48E-06 1,30E-06 2,99E-06 2,75E-06
C6(+) 6,32E-08 3,32E-08 7,48E-08 6,90E-08
C7(+) 3,61E-10 1,90E-10 3,99E-10 3,68E-10
1,90E+00 1,08E+00

φi Z θ φ N
ln γ Ci =ln + q i ln i −l i − i ∑ x i l i
xi 2 φi x i j=1

4,60E-05 10,0000 9,60E-05 4,60E-05

ln γ Ci =ln + 4,60E-05∗ln −1,42E+01− ∗(−4,18E+00)
0,000 2 4,60E-05 xi

ln γ Ci =3,25E+01

Component
Z li xi Θi φ1 xi*li lnγc
e
 0,00 9,60E- 4,60E-
N2 -1,42E+01 -6,67E-04 3,25E+01
0 05 05
0,41 6,57E- 3,85E- -
CO2 -1,00E+01 2,18E+01
7 01 01 4,17E+00
0,02 4,27E- 3,24E-
C1 -9,21E+00 -2,40E-01 1,89E+01
6 02 02
0,54 2,93E- 5,60E-
C2 3,41E-01 1,86E-01 6,70E-01
4 01 01
0,01 6,52E- 2,15E-
C3 3,58E+00 4,34E-02 -1,70E+00
2 03 02
10,000 0,00 8,23E- 2,50E-
iC4 2,71E+00 4,76E-04 -1,36E+00
0 0 05 04
0,00 9,04E- 2,73E-
nC4 2,88E+00 4,95E-04 -1,31E+00
0 05 04
0,00 1,42E- 2,55E-
iC5 7,30E-02 1,97E-07 8,91E-01
0 06 06
0,00 1,30E- 2,75E-
nC5 8,49E-01 2,53E-06 -1,81E-01
0 06 06
0,00 3,32E- 6,90E-
C6(+) 7,38E-01 4,67E-08 2,52E-01
0 08 08
0,00 1,90E- 3,68E-
C7(+) 9,52E-01 4,18E-09 0,00E+00
0 10 10
-
4,18E+00
N
ln γ Gi =∑ v ki (ln T k −ln T (ki) )
k=1

ln T k =1−C k −ln Dk

ln T (ki) =1−C (ki)−ln D (ki)=0

ω ki=exp ( −aT )
ki

a ii=1

Aki
1 2 3 4 5 6 7 8 9 10 11
 1,0 109,1 8,6 21,6 2075,5 160,6 4,1 0,0 48,6 60,6 5,001
109,1 1,0 37,6 94,3 9055,8 700,6 17,7 0,0 212,1 264,5 21,82
8,6 37,6 1,0 7,4 715,1 55,3 1,4 0,0 16,7 20,9 1,723
21,6 94,3 7,4 1,0 1792,9 138,7 3,5 0,0 42,0 52,4 4,32
2075,5 9055,8 715,1 1792,9 1,0 13325,9 337,0 0,4 4034,0 5030,5 415,021
160,6 700,6 55,3 138,7 13325,9 1,0 26,1 0,0 312,1 389,2 32,109
4,1 17,7 1,4 3,5 337,0 26,1 1,0 0,0 7,9 9,8 0,812
0,0 0,0 0,0 0,0 0,4 0,0 0,0 1,0 0,0 0,0 0,001
48,6 212,1 16,7 42,0 4034,0 312,1 7,9 0,0 1,0 117,8 9,72
60,6 264,5 20,9 52,4 5030,5 389,2 9,8 0,0 117,8 1,0 12,121
5,0 21,8 1,7 4,3 415,0 32,1 0,8 0,0 9,7 12,1 1

Desarrollo para componente (1)

ω 11=exp ( −aT )
11

−1,00
ω 11=exp ( 173,6939 K)

ω 11=0,994

−a12
ω 12=exp ( )
T

−109,1
ω 12=exp ( 173,6939 K)

ω 12=0,534

ω 13=exp ( −aT )
13

−8,6
ω 13=exp ( 173,6939 K)

ω 13=0,952

−a14
ω 14=exp ( )
T
 −21,6
ω 14=exp ( 173,6939 K)

ω 14=0,883

−a 15
ω 15=exp ( )
T

−2075,5
ω 15=exp ( 173,6939 K)

ω 15=0,000

ω 16=exp ( −aT )
16

−160,6
ω 16=exp ( 173,6939 K)

ω 16=0,397

ω 17=exp ( −aT )
17

−4,1
ω 17=exp ( 173,6939 K)

ω 17=0,977

−a 18
ω 18=exp ( )
T

−0,0
ω 18=exp ( 173,6939 K)

ω 18=1,000

ω 19=exp ( −aT )
19
 −48,6
ω 19=exp ( 173,6939 K)

ω 19=0,756

−a110
ω 110=exp ( )
T

−60,6
ω 110=exp ( 173,6939 K)

ω 110=0,705

ω 111=exp ( −aT )
111

−5,0
ω 111=exp ( 173,6939 K)

ω 111=0,972

wki
1 2 3 4 5 6 7 8 9 10 11
0,994 0,534 0,952 0,883 0,000 0,397 0,977 1,000 0,756 0,705 0,972
0,53
0,064 0,805 0,581 0,000 0,018 0,903 1,000 0,295 0,218 0,882
4
0,95
0,805 0,983 0,958 0,016 0,727 0,992 1,000 0,908 0,887 0,990
2
0,88
0,581 0,958 0,898 0,000 0,450 0,980 1,000 0,785 0,740 0,975
3
0,00
0,000 0,016 0,000 0,000 0,000 0,144 0,998 0,000 0,000 0,092
0
0,39
0,018 0,727 0,450 0,000 0,003 0,861 1,000 0,166 0,106 0,831
7
0,97
0,903 0,992 0,980 0,144 0,861 0,996 1,000 0,956 0,945 0,995
7
1,00
1,000 1,000 1,000 0,998 1,000 1,000 1,000 1,000 1,000 1,000
0
 0,75
0,295 0,908 0,785 0,000 0,166 0,956 1,000 0,580 0,507 0,946
6
0,70
0,218 0,887 0,740 0,000 0,106 0,945 1,000 0,507 0,429 0,933
5
0,97
0,882 0,990 0,975 0,092 0,831 0,995 1,000 0,946 0,933 0,994
2

N
Dk =∑ θi ωik
i=1

D1=θ 1 ω 11 +θ2 ω12+ θ3 ω13+ θ4 ω 14+θ 5 ω 15+θ 6 ω16 +θ7 ω17 +θ8 ω 18+θ 9 ω 19+θ 10 ω 110 +θ11 ω111

D 1=9,60E-05∗0,994 +6,57E-01∗0,534 +4,27E-02∗0,952+2,93E-01∗0,883+ 6,52E-03∗0,000+ 8,23E-05∗0,3

D1=8,3E-05

Di
1 2 3 4 5 6 7 8 9 10 11 Di
8,3E-05 5,1E-05 9,1E-05 8,5E-05 6,2E-10 3,8E-05 9,4E-05 9,6E-05 7,3E-05 6,8E-05 9,3E-05
3,5E-01 4,2E-02 5,3E-01 3,8E-01 1,5E-23 1,2E-02 5,9E-01 6,6E-01 1,9E-01 1,4E-01 5,8E-01
4,1E-02 3,4E-02 4,2E-02 4,1E-02 7,0E-04 3,1E-02 4,2E-02 4,3E-02 3,9E-02 3,8E-02 4,2E-02
2,6E-01 1,7E-01 2,8E-01 2,6E-01 9,6E-06 1,3E-01 2,9E-01 2,9E-01 2,3E-01 2,2E-01 2,9E-01
0,0E+0

72,30685545
4,2E-08 1,5E-25 1,1E-04 2,1E-07 0 3,1E-36 9,4E-04 6,5E-03 5,3E-13 1,7E-15 6,0E-04
3,3E-05 1,5E-06 6,0E-05 3,7E-05 3,9E-38 2,2E-07 7,1E-05 8,2E-05 1,4E-05 8,8E-06 6,8E-05
8,8E-05 8,2E-05 9,0E-05 8,9E-05 1,3E-05 7,8E-05 9,0E-05 9,0E-05 8,6E-05 8,5E-05 9,0E-05
1,4E-06 1,4E-06 1,4E-06 1,4E-06 1,4E-06 1,4E-06 1,4E-06 1,4E-06 1,4E-06 1,4E-06 1,4E-06
9,8E-07 3,8E-07 1,2E-06 1,0E-06 1,1E-16 2,2E-07 1,2E-06 1,3E-06 7,6E-07 6,6E-07 1,2E-06
2,3E-08 7,2E-09 2,9E-08 2,5E-08 8,8E-21 3,5E-09 3,1E-08 3,3E-08 1,7E-08 1,4E-08 3,1E-08
1,8E-10 1,7E-10 1,9E-10 1,9E-10 1,7E-11 1,6E-10 1,9E-10 1,9E-10 1,8E-10 1,8E-10 1,9E-10
1,7E+0 2,2E+0 3,9E+0 4,7E+0 5,0E+0 6,2E+0 7,9E+0 9,0E+0 9,5E+0 1,0E+0 1,2E+0
0 0 0 0 0 0 0 0 0 1 1

N
θi ωik
C k =∑
i=1 Di

θ 1 ω 11 +θ2 ω12 +θ3 ω13 +θ4 ω 14+ θ5 ω15+ θ6 ω16 +θ7 ω17 +θ8 ω 18+θ 9 ω 19+θ 10 ω 110 +θ11 ω 111 θ 1 ω 21+θ2 ω 22+θ3 ω 23+
C 1= +
D1
 9,5E-05+θ2 ω12+ θ3 ω13 +θ4 ω 14+θ 5 ω15+ θ6 ω16 +θ7 ω17 +θ8 ω 18+θ 9 ω 19+θ 10 ω 110 +θ11 ω111 θ1 ω 21+θ2 ω22 +θ3 ω
C 1= +
72,30685545

θi wij Ck
1,650444 43,81054
03 67
2,247330 32,17455
69 08
3,852608 18,76828
45 56
4,686470 15,42885
76 02
5,000720 14,45928
45 77
6,174748 11,71008
78 86
7,924457 9,124518
23 35
8,999894 8,034189
14 55
9,463041 7,640973
87 85
10,39829 6,953720
71 85
11,90884 6,071694
19 95

ln γ i=ln γ ic + ln γ Gi

Calculo de coeficiente de actividad

Componente yi lnγc lnγG lnγ γ

N2 4,7073E-05 32,5088365 -179,112046 -10,8027546 0
CO2 0,41690409 21,8091962 -25,5874351 -10,1750977 0
 C1 0,02610099 18,9374793 -15,992147 -0,1795567 0,83564057
C2 0,54446435 0,66987753 -12,7937176 -15,3036525 2,2579E-07
C3 0,01213035 -1,69942717 -14,3929323 -16,7682968 5,2194E-08
iC4 0,00017564 -1,36401114 -9,59528818 -13,894568 9,2399E-07
nC4 0,00017175 -1,31412503 -7,99607348 -11,5085972 1,0043E-05
iC5 2,7041E-06 0,89066169 -6,39685879 -8,34074067 0,0002386
nC5 2,9794E-06 -0,18052128 -4,79764409 -9,06888901 0,00011519
C6(+) 6,3298E-08 0,25152816 -3,19842939 -8,04383474 0,00032108
C7(+) 4,3909E-09 0 -0,31984294 0 0

Calculo de constante de equilibrio:

γ 1 ∅1s Ps1 ( FP)1

K 1= v
∅1 P

Zi Ki
y i=
(1+ β )(K i−1)

Zi
x i=
(1+ β )( K i−1)

Verificación∑ ( xi − y i )=0

CALCULO DE CONDICION DE EQUILIBRIO

N
Z i (K i−1)
f (FVAP)=∑ =0
i=1 1+ β( K i−1)

Primera iteración del proceso real

B,SUP 0,5
Componente KI Zi f(B) Xi yi Verificacion
0,0002802
N2 9,42 0,0016 -0,00263958 0,00291979 -0,762643471
1
0,1310694
CO2 0,06771 0,066 0,130138814 0,000930593
1
0,2439252
C1 4,634 0,8091 -1,130349533 1,374274766
3
0,1501771
C2 0,01347 0,0761 0,148154229 0,002022886
1
0,0597280
C3 0,01125 0,0302 0,059056119 0,000671941
6
0,0065986
iC4 0,00021 0,0033 0,006597229 1,38571E-06
1
0,0155986
nC4 0,000084 0,0078 0,015597379 1,31029E-06
9
0,0041999
iC5 0,00002 0,0021 0,004199832 8,39983E-08
2
0,0039999
nC5 0,0000078 0,002 0,003999938 3,11998E-08
7
C6(+) 0,00000074 0,0008 0,001599998 0,0016 1,184E-09
C7(+) 1 0,001 0,002 0,0010021 0
0,6181793
1 -0,761645576 1,380822788
2

Verificación∑ ( xi − y i )=0,76264

No cumple

Segunda iteración real con suposición de factor de vaporización (B)

 B,SUP 0,7677

Componente KI Zi f(B) Xi yi Verificacion

N2 10,42 0,0016 -0,001830963 0,00019437 0,002025333
CO2 0,0071 0,066 0,275630769 0,27760174 0,001970972
C1 4,634 0,8091 -0,775833154 0,21349289 0,989326042
C2 0,01347 0,0761 0,309407553 0,31363218 0,004224625
C3 0,01125 0,0302 0,123934043 0,12534416 0,001410122 0,001160755
iC4 0,00021 0,0033 0,014192935 0,01419592 2,98114E-06
nC4 0,000084 0,0078 0,033565133 0,03356795 2,81971E-06
iC5 0,00002 0,0021 0,009039256 0,00903944 1,80789E-07
nC5 0,0000078 0,002 0,008609268 0,00860933 6,71528E-08
C6(+) 0,00000074 0,0008 0,003443812 0,00344381 2,54842E-09
C7(+) 0 0,001 0,004304778 0,0010021 0
1 0,004463429 1,0001239 0,998963145

Verificación∑ ( xi − y i )=0,0011= APRO X IMADO A 0

POR LO TANTO, CUMPLE Y NORMALIZAR

Normalización
Componente Xi Yi Verificación
N2 0,000194346 0,002027435
CO2 0,277567351 0,001973018
C1 0,213466439 0,990352894
C2 0,313593324 0,00422901
C3 0,125328636 0,001411585
iC4 0,014194158 2,98424E-06 0
nC4 0,033563794 2,82263E-06
iC5 0,009038317 1,80976E-07
nC5 0,008608268 6,72225E-08
C6(+) 0,003443388 2,55107E-09
C7(+) 0,00100198 0
1 1

 Tabla final de resultados para la iteración final y para factor de

vaporización igual a 0,7677:

T de FLASH
 R K C
T
173,69395 -99,45605
P DEL FLASH
Kpa psi bar
P
5230,888 493,05 52,30888

B,SUP 0,7677
Compone Psat, V[m3/Km Ǿ^v (FP)
Z ∅s γ KI Zi f(B) Xi Yi
nte (KPA) ol] i i
-
1,2 1,21 1,20 0,00 10,42 0,00 0,00 0,00
N2 0,224 0,001 0,8 0,00
1 0 9 0 0 2 0 2
2
1,0 1,01 1,23 0,00 0,06 0,27 0,27 0,00
CO2 48,105 0,087 0,68 0,007
0 2 2 0 6 6 8 2
-
0,69 1,3 1,14 1,03 0,70 0,80 0,21 0,99
C1 0,245 0,316 4,634 0,77
2 6 3 0 8 9 3 0
6
1,0 1,00 1,52 0,00 0,07 0,30 0,31 0,00
C2 54,502 0,386 0,9 0,013
0 1 4 0 6 9 4 4
1,0 0,99 1,06 0,00 0,03 0,12 0,12 0,00
C3 3,052 0,006 0,9 0,011
1 0 3 0 0 4 5 1
0,1 1,63 1,01 0,00 0,00 0,01 0,01 0,00
iC4 0,404 0,076 0,72 0,000
3 8 3 0 3 4 4 0
2,1 2,68 1,01 0,00 0,00 0,03 0,03 0,00
nC4 0,167 0,080 0,6 0,000
1 4 1 0 8 4 4 0
1,0 1,51 1,00 0,00 0,00 0,00 0,00 0,00
iC5 0,010 0,402 0,8 0,000
5 0 8 0 2 9 9 0
1,0 1,49 1,01 0,00 0,00 0,00 0,00 0,00
nC5 0,011 0,411 0,9 0,000
3 9 2 0 2 9 9 0
1,0 1,30 1,01 0,00 0,00 0,00 0,00 0,00
C6(+) 0,001 0,012 0,9 0,000
1 1 8 0 1 3 3 0
1,0 1,39 1,04 0,00 0,00 0,00 0,00
C7(+) 0,000 0,100 0,71 0 0
5 9 6 1 4 1 0
0,00
1 1 1
4
3.- ESPECIFICACIÓN DESTILACIÓN FLASH MEDIANTE SIMULADOR ASPEN
HYSYS

Elección de paquete de fluidos para resolución del proceso de destilación flash:

Introducción de composición en el simulador

Calculo de propiedades del gas y factor de vaporización (B)

Pantalla de muestreo de k de la fase vapor

Pantalla de muestreo de k la fase liquida

Composición de salida de la fase liquida y vapor

4.- COMPARACION DE RESULTADOS Y DETERMINACION DE PORCENTAJE
DE ERROR

FVAP,cal. FVAP,hysys e %
0,7677 0,838 0,08389021 8,38902148

((FVAP)calculado−( FVAP)hysys)
%e= ∗100
(FVAP)hysys

(0,7677−0,838)
%e= ∗100
0,838

e=8,3890 %