UNIVERSIDAD NACIONAL AUTONOMA DE MÉXICO

FACULTAD DE ESTUDIOS SUPERIORES CUAUTITLAN

CAMPO 1

INGENIERÍA QUÍMCA

ASIGNATURA:
TRANSFERENCIA DE MASA II

GRUPO: 2701

SERIE 4 – Destilación multicomponentes

INTEGRANTES:
Cortés Cabañas Raúl Jair
Miranda Carmona Issac
Trejo Galván José Raúl
Rubio Abelino Eduardo
Yépez Torres Víctor Erubey

FECHA DE ENTREGA: 03/05/2023

12.7. Para la operación de destilación que se indica seguidamente, calcúlese el
número mínimo de etapas de equilibrio y la distribución de los componentes no
clave utilizando la ecuación de Fenske y los valores K de la Fig. 7.5 y
determine el número de etapas y el plato de alimentación si R=1.2 Rmín.

Σ 3930 kmol/h
De inicio tenemos como datos que en el destilado hay 15 kmol/h de HK, por lo tanto, en
los fondos tendríamos 85 kmol/h de iC5. De igual forma que cuando tenemos 6 kmol/h de
LK en los fondos, realizando una resta, obtenemos 594 kmol/h de nC4.
Suponemos que los componentes por debajo de HK se van a los fondos y que los
componentes por arriba de LK se van al destilado obtendríamos que:
kmol
D=2,500+400+ (600−6 )+ 15=3509
h
kmol
B=40+50+ 40+200+ ( 100−15 ) +6=421
h

Composición
Compuestos Zi di=DXi,D bi=BXi,B Xi,D Xi,B
(kmol/h)
C3 2500 0.63613232 2500 0 0.71245369 0
iC4 400 0.10178117 400 0 0.11399259 0
LK nC4 600 0.15267176 594 6 0.169279 0.01425178
HK iC5 100 0.02544529 15 85 0.00427472 0.20190024
nC5 200 0.05089059 0 200 0 0.47505938
nC6 40 0.01017812 0 40 0 0.09501188
nC7 50 0.01272265 0 50 0 0.11876485
nC8 40 0.01017812 0 40 0 0.09501188
Σ 3930 1 3509 421 1 1
Las fracciones molares se calcularon de la siguiente manera:

di bi
X Di= X Bi=
D B
15 85
XD = =0.00427 X B = =0.2019
nC 4
3509 421
nC 4

Suponemos unas temperaturas de 123°F a 340 °F. Mediante la figura 7.5 obtenemos los
valores de K de los componentes.

Posteriormente con la definición de volatilidad relativa:

 ki
α i , HK =
k HK

Donde:

k HK =k i C =0.36 @ 123 ° F
5

Destilado a 123 °F Fondos 340 °F

Destilado Ki Fondos Ki
C3 2.3 6.38888889 C3 8.2 3.15384615
iC4 1.1 3.05555556 iC4 5.2 2
nC4 LK 0.76 2.11111111 nC4 4.4 1.69230769
iC5 HK 0.36 1 iC5 2.6 1
nC5 2.8 7.77777778 nC5 2.3 0.88461538
nC6 0.11 0.30555556 nC6 1.3 0.5
nC7 0.045 0.125 nC7 0.74 0.28461538
nC8 0.19 0.52777778 nC8 0.45 0.17307692

Ahora calculamos una volatilidad relativa media, multiplicando la volatilidad relativa

del destilado por la del fondo para cada componente:
1 /2
α i , HK =[ ( α i , HK )D ( α i , HK ) B ]

Obtenemos:
Componen Volat.
te media
C3 4.488827
54
iC4 2.472066
16
nC4 1.890145
38
iC5 1
nC5 2.623040
58
nC6 0.390867
98
nC7 0.188618
46
nC8 0.302235
26
Posteriormente, con la ecuación de Fenske calculamos el número mínimo de etapas:

N min + 1=
ln
[( )( )]
x D , LK
x D , HK
x B , LK
x B , HK
ln ⁡(α LK , HK )

N min + 1=
ln
[( 0.169279
0.00427472 0.01425178 )]
)( 0.20190024
ln ⁡(1.89014538)
N min + 1=9.94217

N min =8.94217 ≈ 9 etapas

Para el cálculo de reflujo mínimo

Con la ecuación de Underwood:

NC
α z
∑ α −θ
i i
=1−q
i=1 i

Donde:
NC: número de componentes
αi: volatilidad relativa promedio de i respecto a HK
zi: composicion
Suponiendo que se trata con líquido saturado: q=1

NC
α z
∑ α −θ
i i
=0
i=1 i

Sustituyendo:

4.488(
2500 2.472 400
3930
)
+
(
3930
+
)
1.8901(
600
3930
) 1(
+
100
3930
) 2.623(
+
200
3930
) 0.3908 (
+
40
3930
) 0.1886(
+
50
3930
) 0.3
+
4.488−θ 2.472−θ 1.8901−θ 1−θ 2.623−θ 0.3908−θ 0.1886−θ 0
Calculamos:
θ=1.01933
Con el valor calculado utilizamos la segunda ecuación de Underwood:
NC
α x
∑ α i−θD =1+ R min
i

i=1 i

Teniendo todos los valores podemos calcular Rmin:

4.488 (0.7125) 2.472(0.1139) 1.8901(0.1693) 1(0.0043) 2.623 (0) 0.3908 (0)
+ + + + + +
4.488−1.01933 2.472−1.01933 1.8901−1.01933 1−1.01933 2.623−1.01933 0.3908−1.01933 0.1

Rmin =0.2938
Finalmente calculamos el número de etapas teóricas con la relación de reflujo:

CORRELACIÓN DE GILLILAND:

N −N min
N +1
=0.75 1−
[ (
R−Rmin
R +1 )]
Donde:
N min=8.94

Rmin =0.2938

R=1.2 Rmin =0.35256

Sustituimos en la ecuación y obtenemos el número de etapas teóricas:

N −8.94
N +1
=0.75 1−
[ (
0.35256−0.2938
0.35256+1 )]
N=29.098 etapasteóricas
-Problema 15.8 del Henley
Utilizando el método BP de Wang-Henke calcúlense composición de los
productos, temperaturas de las etapas, composiciones y flujos de las corrientes de
vapor y líquido interetapas, servicio del ebullidor y servicio del condensador para
las especificaciones de la siguiente columna de destilación.
Alimentación (liquido saturado a 250 psia y 213.9°F)

Componente lb mol/h
Etano 3.0
Propano 20.0
n-Butano 37.0
n-Pentano 35.0
n-Hexano 5.0

Presión de la columna=250 psia Condensador parcial y ebullidor parcial

Flujo de destilado=23 lb mol/h
Reflujo=150 lb mol/h Número de etapas de equilibrio (excluidos el condensador y
ebullidor )=15
La alimentación se envia a la etapa intermedia

Para este sistema a 150 psia, los valores K y las entalpías pueden calcularse
dentro de un intervalo de temperaturas de 50°F a 350°F mediante las ecuaciones
polinómicas del ejemplo 12.8 utilizando las siguientes constantes (Amundson y
Pontinen).

Constantes de los valores K Constantes de entalpia de

vapor
Especies αi 4
β i x 10 γi x 10
6 8
δ i x 10 A B C2
C2 1.665 -1.50 73.5 -3.00 8310 13.2 60.2
C3 0.840 46.6 49.4 -3.033 11950 12.35 221
nC4 -0.177 49.5 -4.15 2.22 16550 9.01 378
nC5 -0.0879 17.7 0.2031 1.310 19200 29.5 72
nC6 0.0930 -15.39 10.37 -0.1590 22000 37 0

Constantes de entalpia del liquido

Especies a b c2
C2 6120 9.35 150
C3 3620 40.5 -132
nC4 6280 25 465
nC5 6840 38.5 252
 nC6 7920 43.5 302

Solución
Lo primero que debemos suponer son temperaturas y valores V, en la diapositiva
se nos proponen los siguientes valores. Se supone que se tienen flujos de vapor
internos constantes. Es claro que para el plato 1 tendremos un flujo de vapor igual
al destilado y que a partir del segundo plato será la suma del reflujo + el valor del
destilado, como se muestra a continuación:
Plato T(°F) V
1 110 23
2 120 173
3 130 173
4 140 173
5 150 173
6 160 173
7 170 173
8 180 173
9 190 173
10 200 173
11 210 173
12 220 173
13 230 173
14 240 173
15 250 173
16 260 173
17 270 173

Los polinomios que nos ayudaran a obtener los valores de la constante K, así
como la constante K, así como las entalpías de las entalpías del líquido y del
vapor son:
2 3
K i=α i + β i T +γ i T +δ i T
2
H vi = Ai + Bi T +C i T
2
H Li=ai+ bi T + ci T

Componente α β γ δ
Etano 1.665 -0.00015 0.00007350 -0.00000003
Propano 0.84 -0.00466 0.00004940 -0.00000003
Butano -0.177 0.004950 -0.00000415 0.000000022
Pentano -0.0879 0.001170 0.00000020 0.000000013
Hexano 0.093 -0.001539 0.00001037 -0.0000000015

Posteriormente, con los datos proporcionados de flujo de cada componente,

obtenemos sus fracciones molares a la entrada de la torre:
 Componente lb mol/h Zi
Etano 3.0 0.03
Propano 20.0 0.2
n-Butano 37.0 0.37
n-Pentano 35.0 0.35
n-Hexano 5.0 0.05

Con la ecuación de Antoine y la presión proporcionada por el ejercicio, podemos

obtener la temperatura para cada componente, para así, verificar si estamos en el
rango correcto de temperaturas:
La presión 250 psia=17.2369 mmHg
B
T= −C
A−ln ⁡(P)
A B C T(°F)
Etano 3.9540 663.72 256.681 9.59004
Propano 3.9282 803.997 247.04 124.9546
n-Butano 3.9326 935.773 238.789 226.9075
n-Pentano 3.9778 1064.84 232.014 313.5472
n-Hexano 4.0014 1170.875 224.317 390.4817
Si realizamos un promedio de las 5, el rango de temperaturas es correcto
MATRIZ 1
A j x i , j−1+ B j x i , j +C x i , j+1=D j
En donde: j=Platoi=Componente
Así mismo:
j−1
A j=V j+ ∑ (F ¿ ¿ m−W m −U m )−V i 2 ≤ j≤ N ¿
m=1
Lo que nos indica que empezaremos a partir del segundo plato. Del segundo al
noveno plato tendremos:V-Vinicial y a partir del décimo plato tendremos:F+V-Vinicial
j
B j =−V j +1+ ∑ (F ¿ ¿ m−W m −U m )−V 1 +U j+ ( V j −W j ) ( K i , j ) 1 ≤ j≤ N ¿
m=1
Lo que nos indica es que realizaremos el calculo para todos los platos del primero
al octavo plato tendremos: B j =−[ V j +1−D+ ( V j K i , j ) ]y a partir del noveno plato en
adelante tendremos: : B j =−[ V j +1−D+ ( V j K i , j ) ]+ F

C j=V j+1 K i , j 1 ≤ j≤ N −1
Lo que nos indica que el calculo se realizará para todos los platos excepto para el
ultimo, en donde será necesario ocupar el valor de V y K del plato posterior al
plato en donde necesitamos obtener ese dato.

D j=F j z i , j 1≤ j ≤ N
Donde solo será posible obtener el valor correspondiente al plato de alimentación
porque es el único plato en donde tenemos entrada de flujo de alimentación.
Con todos estos valores obtenidos, armamos nuestra matriz de la siguiente
manera, para así obtener las composiciones de nuestro componente (en este caso
del etano) en cada plato:

Para poder resolver la matriz podemos ocupar cualquier método, pero se resolverá
por el método de Thomas, el cual ocupa las siguientes ecuaciones:
C1 Cj
p1 = p j =
B1 B j− A j p j−1

Lo que nos quiere decir esta ecuación es que para el plato 1 obtendremos
solamente su C/D correspondiente y para cualquier plato con la segunda
ecuación.
D1 D j− A j q j−1
q 1= q j =
B1 B j −A j p j −1

Lo que nos quiere decir esta ecuación es que para el plato 1 obtendremos
solamente su D/B correspondiente y para cualquier plato con la segunda ecuación.

Finalmente, para obtener la fracción mol del componente i en cada plato, se

calculará con la siguiente ecuación:
x i , j=q j− p j x i , j+1
Es decir, calcularemos las fracciones desde el último plato hacia el primer plato,
dicho en otras palabras, de abajo hacia arriba.
Con las ecuaciones anteriores aplicadas a cada componente, podremos calcular
las composiciones de cada componente, en cada plato.
 ETANO
T (°F) V Ki Ai Bi Ci Di
1 110 23 2.49792 0 -207.45216 459.06588 0
2 120 173 2.65356 150 -609.06588 488.16102 0
3 130 173 2.82174 150 -638.16102 519.39444 0
4 140 173 3.00228 150 -669.39444 552.735 0
5 150 173 3.195 150 -702.735 588.15156 0
6 160 173 3.39972 150 -738.15156 625.61298 0
7 170 173 3.61626 150 -775.61298 665.08812 0
8 180 173 3.84444 150 -815.08812 706.54584 0
9 190 173 4.08408 150 -956.54584 749.955 -3
10 200 173 4.335 250 -999.955 795.28446 0
11 210 173 4.59702 250 -1045.28446 842.50308 0
12 220 173 4.86996 250 -1092.50308 891.57972 0
13 230 173 5.15364 250 -1141.57972 942.48324 0
14 240 173 5.44788 250 -1192.48324 995.1825 0
15 250 173 5.7525 250 -1245.1825 1049.64636 0
16 260 173 6.06732 250 -1299.64636 1105.84368 0
17 270 173 6.39216 250 -1182.84368 0 0

Una vez obtenidos los valores de A, B, C y D armamos nuestra matriz:

MATRIZ 1
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17
1 -207.45216 459.06588 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
2 150 -609.06588 488.16102 0 0 0 0 0 0 0 0 0 0 0 0 0 0
3 0 150 -638.16102 519.39444 0 0 0 0 0 0 0 0 0 0 0 0 0
4 0 0 150 -669.39444 552.735 0 0 0 0 0 0 0 0 0 0 0 0
5 0 0 0 150 -702.735 588.15156 0 0 0 0 0 0 0 0 0 0 0
6 0 0 0 0 150 -738.15156 625.61298 0 0 0 0 0 0 0 0 0 0
7 0 0 0 0 0 150 -775.61298 665.08812 0 0 0 0 0 0 0 0 0
8 0 0 0 0 0 0 150 -815.08812 706.54584 0 0 0 0 0 0 0 0
9 0 0 0 0 0 0 0 150 -956.54584 749.955 0 0 0 0 0 0 0
10 0 0 0 0 0 0 0 0 250 -999.955 795.28446 0 0 0 0 0 0
11 0 0 0 0 0 0 0 0 0 250 -1045.28446 842.50308 0 0 0 0 0
12 0 0 0 0 0 0 0 0 0 0 250 -1092.50308 891.57972 0 0 0 0
13 0 0 0 0 0 0 0 0 0 0 0 250 -1141.57972 942.48324 0 0 0
14 0 0 0 0 0 0 0 0 0 0 0 0 250 -1192.48324 995.1825 0 0
15 0 0 0 0 0 0 0 0 0 0 0 0 0 250 -1245.1825 1049.64636 0
16 0 0 0 0 0 0 0 0 0 0 0 0 0 0 250 -1299.64636 1105.84368
17 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 250 -1182.84368

Posteriormente, obtenemos nuestros valores de p y q correspondientes, para así

obtener la fracción mol del ETANO en cada plato:
Pi qi X i,j
-2.21287587 0 0.052217146
-1.76145887 0 0.023596961
-1.38896988 0 0.01339626
-1.19886401 0 0.009644745
-1.12477622 0 0.008044903
-1.09865541 0 0.007152448
-1.0888542 0 0.006510183
 -1.08405832 0 0.00597893
-0.94460255 0.00377864 0.005515321
-1.04121487 0.00123678 0.001838535
-1.07327866 0.00039389 0.000577933
-1.08177343 0.00011948 0.000171479
-1.08190093 3.4288E-05 4.80691E-05
-1.07936427 9.2971E-06 1.27378E-05
-1.0761835 2.383E-06 3.18763E-06
-1.07300908 5.7807E-07 7.47623E-07
0 1.5801E-07 1.58014E-07

Se realiza el mismo procedimiento para cada componente.

Una vez que tenemos la composición de cada componente en cada plato
realizamos la suma de todos los componentes en cada plato. Resultando:

PLATO SUM x i , j
1 1.005822
2 1.037951
3 1.032561
4 0.989979
5 0.92612
6 0.860611
7 0.811288
8 0.802607
9 0.905174
10 0.926031
11 0.943313
12 0.956386
13 0.962217
14 0.956615
15 0.940657
16 0.935164
17 1.009651
Una vez que sumamos las composiciones en todos los platos, NORMALIZAMOS
las composiciones con la siguiente ecuación:

xi , j
x i , j= c

∑ xi , j
i=1
Los resultados de la normalización de todos los componentes de muestran a
continuación:
  Composiciones normalizadas para cada plato en la torre

PLATO SUM x i , j x i ,j N ORM x i,j NORM x i,j NORM x i,j NORM x i,j NORM
ETANO PROPANO BUTANO PENTANO HEXANO
1 1.00582222 0.05191489 0.93613171 0.01195332 8.58E-08 5.31E-10
2 1.03795071 0.02273418 0.95052724 0.02673764 9.32E-07 8.19E-09
3 1.03256068 0.01297382 0.93345936 0.05355827 8.43E-06 1.10E-07
4 0.98997894 0.00974237 0.89132542 0.09886504 6.59E-05 1.31E-06
5 0.92611994 0.00868668 0.82225625 0.16859611 0.00044753 1.34E-05
6 0.86061121 0.00831089 0.72553985 0.26339047 0.00264024 0.00011855
7 0.81128796 0.0080245 0.60601388 0.37167736 0.01339958 0.00088468
8 0.80260651 0.00744939 0.46960284 0.46041655 0.05709753 0.00543369
9 0.90517363 0.00609311 0.31596218 0.46292915 0.18954973 0.02546584
10 0.92603118 0.00198539 0.26872299 0.51363976 0.19022928 0.02542258
11 0.94331326 0.00061266 0.21379225 0.56776614 0.19228193 0.02554701
12 0.95638634 0.0001793 0.15915288 0.61867433 0.19612324 0.02587025
13 0.96221691 5.00E-05 0.1110871 0.65897049 0.20334477 0.02654768
14 0.95661496 1.33E-05 0.07282037 0.68114297 0.21800784 0.02801551
15 0.94065726 3.39E-06 0.04450808 0.67441894 0.24955945 0.03151015
16 0.93516362 7.99E-07 0.02444165 0.61674471 0.31810285 0.04070999
17 1.00965149 1.57E-07 0.01079273 0.47469243 0.45020035 0.06431434

El siguiente paso es obtener la temperatura de punto con punto de burbuja. En

este punto sabemos que la suma de las fracciones mol en el vapor debe ser 1:
x

∑ yi , j=1
i=1
Así mismo sabemos que:
x

∑ k i , j x i , j=1
i=1

Con la función de SOLVER de Excel, dejando el criterio que la sumatoria de

nuestras fracciones debe ser 1 encontramos la temperatura a la cual se cumpla
esta condición (1° iteración).

Obtenemos lo siguiente:
 PLATO T (° F)
1 117.057646
2 126.64968
3 132.278758
4 137.857336
5 144.057224
6 151.50797
7 160.734539
8 174.592811
9 203.695563
10 213.842612
11 223.376897
12 231.987778
13 239.413756
14 246.175122
15 253.948437
16 266.256077
17 288.747114

Una vez que se ha encontrado esta temperatura proseguimos con el

procedimiento.

El siguiente paso es encontrar los flujos de vapor y de líquido, para lo cual

ocupamos las siguientes ecuaciones:

-Con las temperaturas calculadas anteriormente, podemos encontrar las entalpias

del líquido y del vapor para cada componente con los polinomios:

H V = A i + Bi T + C i T 2
i

2
H L =a i+ bi T +c i T
i

Obtenemos los siguientes resultados:

 Entalpía del liquido
Etano Propano Butano Pentano Hexano
7420.02638 8179.96176 9843.60705 11692.0222 13425.8229
7544.77663 8537.58217 10192.1086 12120.2242 13913.6734
7619.27144 8746.32047 10400.6106 12373.6735 14202.5556
7694.03577 8952.36079 10610.1494 12626.4245 14490.7344
7778.2223 9180.38479 10846.4211 12909.1657 14813.2143
7880.9195 9453.07122 11135.0912 13251.5144 15203.8276
8010.40182 9788.71901 11499.7185 13679.3367 15692.1873
8209.68253 10288.6379 12062.2635 14329.986 16435.3633
8646.93175 11321.9775 13301.7616 15727.8746 18033.8119
8805.35835 11677.0074 13752.4481 16225.3028 18603.1592
8957.03256 12008.1208 14184.644 16697.4209 19143.7916
9096.36066 12305.1031 14582.2468 17127.7513 19636.7819
9218.30281 12559.647 14930.6749 17501.8671 20065.5306
9330.77025 12790.1435 15252.3799 17844.9174 20458.804
9647621.88 13406573.6 165263228 18827472.3 210694134

Así mismo para el vapor:

Entalpía del vapor

Etano Propano Butano Pentano Hexano
9937.65 13698.49 18122.64 22751.86 26331.13
10078.34 13868.61 18297.43 23051.65 26686.04
10161.42 13970.34 18403.24 23228.21 26894.31
10244.12 14072.54 18510.47 23403.62 27100.72
10336.49 14187.74 18632.4 23599.11 27330.12
10448.09 14328.42 18782.77 23834.76 27605.79
10587.23 14506.04 18974.8 24127.69 27947.18
10798.13 14779.89 19275.33 24569.96 28459.93
11248.56 15382.61 19953.69 25507.76 29536.74
11408.01 15601.56 20205.27 25837.6 29912.18
11558.96 15811.43 20448.74 26148.88 30264.95
11696.23 16004.43 20674.54 26431.13 30583.55
11815.32 16173.51 20873.77 26675.4 30858.31
11924.34 16329.57 21058.8 26898.5 31108.48
12050.35 16511.49 21275.79 27155.81 31396.09
12251.35 16804.98 21628.7 27564.98 31851.47
12623.38 17358.61 22303.18 28318.34 32683.64
Para aplicar las siguientes ecuaciones necesitamos obtener una entalpía global de
plato, por lo cual, con las fracciones en la fase vapor anteriormente obtenidas,
podemos encontrar una “entalpía de mezcla” para el líquido y para el vapor:

Entalpia global
Liquido Vapor
8163.58471 13239.9529
8563.16901 13753.2829
8825.87986 14082.0458
9113.8321 14409.4971
9468.48827 14707.1336
9923.87989 15004.3865
10512.8002 15332.0466
11410.7449 15884.8921
13246.5058 17637.2742
13835.089 18457.4872
14402.3515 19254.2781
14925.4338 19990.9867
15387.8978 20626.7121
15825.7952 21178.4666

Una vez que obtuvimos la enta.lpia global de la mezcla lo que prosigue es obtener
la segunda matriz para así obtener los flujos V y L. Para esto se necesitan obtener
los siguientes valores:
Se calcula el nuevo conjunto de variables de tanteo V mediante la aplicación del
siguiente balance de energía modificado, obtenido por combinación y reordenado:

α i V i + β i V i+1=γ i
En donde:
γ i=H Lj−1−H Vj β i=H Vj+1−H Lj

[∑ ( ]
j−1
γ i= Fm −W m −U m )−V 1 ( H Lj −H Lj −1) + F j ( H Lj−H F )+ W j ( H Vj −H Lj ) +Q j
m=1

Se realiza la matriz y se resuelve para cada valor de V:

Así mismo, basándonos en Henley podemos obtener los flujos de V mediante:
La ecuación matricial se resuelve de forma inmediata, comenzado a resolver por la
ecuación superior donde V2 es conocida y operando hacia debajo de la misma
forma:
Sabemos que:
V
V1 23
V2 173

Por lo cual podemos obtener los diferentes flujos de V y L de la siguiente manera:

N Alfa Beta Gamma V L

1 V1 23 150
2 -5609.01058 5540.45372 -9173.70815 V2 173 150.484911
3 -5540.45372 5620.08217 -6004.77997 V3 173.484911 146.95843
4 -5620.08217 5665.3249 -6530.10912 V4 169.95843 144.448514
5 -5665.3249 5668.20616 -7971.52583 V5 167.448514 142.957039
6 -5668.20616 5619.79622 -10186.8727 V6 165.957039 142.573946
7 -5619.79622 5660.90404 -13202.0036 V7 165.573946 139.03946
8 -5660.90404 6407.61867 -20381.7699 V8 162.03946 116.975256
9 -6407.61867 5523.76933 -150237.277 V9 139.975256 212.174143
10 -5523.76933 5867.66795 42380.6092 V10 135.174143 211.474446
11 -5867.66795 6063.80751 42144.865 V11 134.474446 214.074974
12 -6063.80751 6108.05634 40261.2573 V12 137.074974 219.673458
13 -6108.05634 6084.49018 36753.0232 V13 142.673458 226.266498
14 -6084.49018 6115.73675 35258.1717 V14 149.266498 231.269019
15 -6115.73675 6350.57554 42192.9877 V15 154.269019 232.208247
16 -6350.57554 6987.99089 68895.0804 V16 155.208247 227.909868
17 V17 150.909868 77
Una vez obtenidos estos flujos, debemos saber que las temperaturas obtenidas
son las correctas, entonces se realiza el siguiente calculo:
N
τ =∑ [T kj ¿−T k−1 2
j ] ≤ 0.01 N ¿
j=1

Por lo cual nuestra tolerancia para el número de platos deberia ser menor a 0.17
Para la primera iteración tenemos:

Tcal-Tsup Tcal-Tsup^2
7.05764599 49.81036686
6.64967996 44.21824352
2.27875843 5.192739972
-2.14266417 4.591009758
-5.94277577 35.31658381
-8.49202959 72.11456649
-9.26546117 85.84877062
-5.40718899 29.2376928
13.6955633 187.5684533
13.8426115 191.6178924
13.3768969 178.9413698
11.9877779 143.7068184
9.41375572 88.61879672
6.17512186 38.13212993
3.94843725 15.59015671
6.25607701 39.13849961
18.7471139 351.4542792
1561.09837

Nuestra suma nos dio 1561.0983 por lo cual comenzamos a iterar ya que no está
dentro de la tolerancia permitida.

-Con los nuevos valores de flujos V y T obtenidos obtenemos las fracciones

normalizadas de x.
-Duplicamos esa hoja de cálculo y con SOLVER obtenemos las nuevas
temperaturas calculadas.
-Con estas nuevas temperaturas obtenemos las entalpias de líquido y vapor.
-Para obtener las nuevas entalpias de líquido y vapor globales realizamos el
mismo calculo que al principio, ocupando las constantes K obtenidas con la nueva
temperatura y con las fracciones “x” y “y” obtenidos con la temperatura y V
anteriores.
-Una vez obtenidas estas entalpias globales podemos realizar la segunda matriz
para obtener alfa, beta y gamma, para finalmente obtener los nuevos flujos V y L.
-Realizamos el calculo final de diferencia de temperaturas al cuadrado para así
corroborar si nuestra temperatura es la correcta.
Después de iterar aproximadamente 15 veces, obtenemos:
Temperaturas en cada etapa:
PLATO T ( °F)
1 117.057646
2 126.64968
3 132.278758
4 137.857336
5 144.057224
6 151.50797
7 160.734539
8 174.592811
9 203.695563
10 213.842612
11 223.376897
12 231.987778
13 239.413756
14 246.175122
15 253.948437
16 266.256077
17 288.747114

Composiciones de cada plato en líquido:

X i,j
Plato T (°F) Etano Propano Butano Pentano Hexano
1 119.08510 0.04951162 0.91544697 0.03504106 3.51E-07 2.28E-09
9
2 129.79705 0.02138774 0.90612561 0.07248342 3.21E-06 3.09E-08
3
3 137.31040 0.01211877 0.85362371 0.13423194 2.52E-05 3.67E-07
3
4 145.51142 0.00900176 0.76555569 0.22526695 0.0001718 3.79E-06
7 2
5 155.3887 0.00778389 0.65048193 0.34070029 0.0010006 3.32E-05
5
6 166.57903 0.0071247 0.52468305 0.46299556 0.0049531 0.0002435
7 2
7 178.73434 0.00662371 0.40579659 0.56524446 0.0208535 0.0014816
3 9 5
8 193.58486 0.00613386 0.30237543 0.61100622 0.0731968 0.0072876
6 1 7
9 216.40665 0.00555248 0.21373706 0.55506909 0.1990761 0.0265652
4 1 6
10 225.39335 0.0016765 0.16479007 0.60666747 0.2002695 0.0265963
9 7 9
11 232.92535 0.00048054 0.11948221 0.65167736 0.2017068 0.0266530
 8 6 4
12 239.24113 0.0001319 0.08218948 0.6869219 0.2039808 0.0267758
5 6
13 244.59733 3.49E-05 0.05393751 0.70951436 0.2093249 0.0271882
8 5 8
14 249.78650 8.93E-06 0.03375755 0.71473763 0.2230091 0.0284867
1 2 7
15 256.62065 2.19E-06 0.01987697 0.69162358 0.2562594 0.0322378
6 7
16 268.47055 5.02E-07 0.01057777 0.61911923 0.3281190 0.0421834
2 2 8
17 289.57740 9.92E-08 0.00465486 0.47616084 0.4542861 0.0648980
3 8 3
Composiciones de cada plato en vapor:
Y i,j
Plato T (°F) Etano Propano Butano Pentano Hexano
1 119.085109 0.13065109 0.85539107 0.01370491 2.68E-08 1.23E-10
2 129.797053 0.06027499 0.90710213 0.03219167 3.08E-07 1.99E-09
3 137.310403 0.03578087 0.89887254 0.06468835 2.79E-06 2.68E-08
4 145.511427 0.02796847 0.8531695 0.11799694 2.18E-05 3.18E-07
5 155.3887 0.0257167 0.77725221 0.19599232 0.00014806 3.27E-06
6 166.57903 0.0252276 0.67911089 0.29401392 0.00085782 2.85E-05
7 178.734343 0.02526895 0.57300473 0.3967551 0.00422298 0.00020762
8 193.584866 0.02559503 0.47446698 0.48072492 0.01765807 0.00125461
9 216.406654 0.02648888 0.39277428 0.51335656 0.06123028 0.00609625
10 225.393359 0.00841876 0.32167041 0.59578907 0.06731647 0.00677495
11 232.925358 0.00251735 0.24511065 0.67212112 0.07285426 0.00738425
12 239.24113 0.0007156 0.17566302 0.7375506 0.0781285 0.00794108
13 244.597338 0.00019499 0.11930069 0.78780226 0.08417623 0.00852929
14 249.786501 5.13E-05 0.07715391 0.81944531 0.09392828 0.00942391
15 256.620656 1.31E-05 0.04740219 0.826592 0.11457488 0.01141537
16 268.470552 3.19E-06 0.02710665 0.79395288 0.1622022 0.0167253
17 289.577403 7.00E-07 0.01348289 0.68923418 0.26622884 0.03104126

Las sumatorias de “y” en cada plato:

 Plato T (°F) Sum (y i,j)
1 119.085109 0.99974709
2 129.797053 0.9995691
3 137.310403 0.99934456
4 145.511427 0.99915704
5 155.3887 0.99911256
6 166.57903 0.99923873
7 178.734343 0.99945939
8 193.584866 0.99969961
9 216.406654 0.99994624
10 225.393359 0.99996966
11 232.925358 0.99998763
12 239.24113 0.9999988
13 244.597338 1.00000345
14 249.786501 1.00000271
15 256.620656 0.99999751
16 268.470552 0.99999022
17 289.577403 0.99998786
Los flujos de líquido y vapor:

PLATO V L
1 23 150
2 173 151.271239
3 174.271239 148.61678
4 171.61678 143.934406
5 166.934406 138.51854
6 161.51854 133.118065
7 156.118065 126.958154
8 149.958154 117.387061
9 140.387061 226.152334
10 149.152334 230.51129
11 153.51129 234.37921
12 157.37921 237.874352
13 160.874352 240.454551
14 163.454551 241.234471
15 164.234471 239.120425
16 162.120425 233.979603
17 156.979603 77

Finalmente, la diferencia de temperaturas al cuadrado y la suma de estas es:

 Tcal-Tsup Tcal-Tsup^2
0.040759518 0.001661338
0.066834039 0.004466789
0.103393407 0.010690197
0.14200423 0.020165201
0.169849778 0.028848947
0.177276648 0.03142701
0.16339555 0.026698106
0.133044447 0.017700825
0.087186622 0.007601507
0.088287905 0.007794754
0.079617303 0.006338915
0.064922933 0.004214987
0.048762328 0.002377765
0.034688784 0.001203312
0.024377471 0.000594261
0.017618563 0.000310414
0.012318783 0.000151752
0.170024608
Las entalpias de líquido y vapor (globales) para cada plato:
 PLATO T (° F) Liquido Vapor
1 119.085109 8273.19725 13301.731
2 129.797053 8750.1805 13837.5131
3 137.310403 9137.57534 14210.9556
4 145.511427 9594.39634 14630.8957
5 155.3887 10158.8475 15178.8597
6 166.57903 10807.0223 15850.1861
7 178.734343 11529.4489 16599.7451
8 193.584866 12466.6845 17449.3195
9 216.406654 14013.3547 18636.601
10 225.393359 14541.6313 19390.3211
11 232.925358 14997.6269 20030.8978
12 239.24113 15387.8035 20575.0856
13 244.597338 15727.0776 21030.8389
14 249.786501 16070.4005 21440.8601
15 256.620656 16548.497 21918.0134
16 268.470552 17411.7205 22699.4372
17 289.577403 18987.9933 24187.5133

SERVICIOS DEL EBULLIDOR Y CONDENSADOR:

-Para obtener Q1 tenemos la ecuación:

Q1=F 1 H F 1 +V 2 hV 2−( L1 +U 1 ) h L1−( V 1+W j ) hV 1

Despreciando términos:

Q1=V 2 hV 2 −( L1 h L1 ) −(V 1 hV 1 )

Sustituyendo valores tenemos:

BTU
Q 1=846,970
h

-Para obtener QN usamos la ecuación:

N
Q N =∑ ( F j h F , j −U j h L, J ) + L N h L, N
j=1
Despreciando términos:
 Q N =F 9 hF 9−( Q1−V 1 hV 1−L17 h17 )

Sustituyendo obtenemos:

BTU
Q N =2,350,996.327
h

Xi,j vs j
nC2 nC3 nC4 nC5 nC6
18
16
14
12
j (platos)

10
8
6
4
2
0
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
Xi,j