0 °C Determinar ???

T
273.15 K Numero de etapas
1.3 atm
P
988 mmHg

% a Eliminar 0.98 𝑃 𝑁𝐻3=𝑦∗𝑃𝑡𝑜𝑡 --> y= P

G 40 kmol/h NH3/Ptot
L/G=1.5 (Lmin/G) Fracciones molares Fracciones molares Libres d
MM NH3 17 g yN+1 0.259 YN+1
MM H20 18 g y1 0.0052 Y1

yA=pA/PT
Solubilidad del NH3 Y X
P NH3
en agua g / 100 g de y NH3 x NH3 Libre de Libre de
(mmHg)
agua sto sto
100 947 0.96 0.51 23.10 1.06
90 785 0.79 0.49 3.87 0.95
80 636 0.64 0.46 1.81 0.85
70 500 0.51 0.43 1.02 0.74
60 380 0.38 0.39 0.63 0.64
50 275 0.28 0.35 0.39 0.53
40 190 0.19 0.30 0.24 0.42
30 119 0.12 0.24 0.14 0.32
25 89.5 0.09 0.21 0.10 0.26
20 64 0.06 0.17 0.07 0.21
15 42.7 0.04 0.14 0.05 0.16
10 25.1 0.03 0.10 0.03 0.11
7.5 17.7 0.02 0.07 0.02 0.08
5 11.2 0.01 0.05 0.01 0.05
0 0 0.00 0.00 0.00 0.00

Recta de Operación

Top
X0 0
Y1 0.0052

Botton
XN 0.5115 5115
YN+1 0.35

Linea Auxiliar de Botton con la botonera busco el valor de

0.35 0 XN en el que la recta de
operacion se encuentra con la
0.35 0.6 curva de eq

Recta de Operación Caudal Minimo de Agua

y = 0.674x + 0.005

L´/G´ 0.67
G´ 40 kmol/h constante
L´ 26.93 kmol/h

Recta de Operación Utilizada

L=1.5 Lmin= 1.5 L´

Factor 1.5
L 40.39 kmol/h

1.00
0.90
G 40 constante

𝐿/𝐺=(𝑌𝑛+1
−𝑌1)/(𝑋𝑁 −𝑋0)
L/G 1.01 1.00
0.90
0.80
Top
0.70
X0 0
Y1 0.0052 0.60
0.50

𝑋𝑁=(𝑌𝑛+1

Y
Botton

−𝑌1)/(𝐿/𝐺)+𝑋0
0.40
XN 0.34
YN+1 0.35 0.30 f(x) = 1.00973810

0.20
0.10
0.00
0.00 0.10 0

Metodo Grafico Numero de Etapas

Numero de Etapas 3 Etapas

0.60

0.50

0.40

f(x) = 1.00973810323044 x + 0.0052069721155

auxiliar 1 0.30
x1 y1
Y

0.341 0.350 0.20

0.341 0.155
0.10
155
1
0.00
auxiliar 2 0.00 2
0.10 0.20 0.30 0.40 0
3
x1 y1 -0.10
0.34 0.155
X
0.148 0.155

148

auxiliar 3
x1 y1
0.148 0.155
0.148 0.042
 42

auxiliar 4
x1 y1
0.148 0.042
0.043 0.042

43

auxiliar 5
x1 y1
0.043 0.042
0.043 0.008

auxiliar 4
x1 y1
0.043 0.008
Este es el objetivo 0.0052 0.008
debe coincidir con dato
52
 yAn+1=0,0052 xA1=0

G L

𝑋𝑁𝐻3=(𝑚𝑁𝐻3/𝑚𝑚𝑁𝐻3)/(𝑚𝑁𝐻3/
𝑚𝑚𝑁𝐻3+100𝑔/𝑚𝑚𝐻20)=(𝑀𝑜𝑙𝑒𝑠 𝑑𝑒
𝑁𝐻3)/(𝑀𝑜𝑙𝑒𝑠 𝑑𝑒 𝑠𝑜𝑙𝑢𝑐𝑖𝑜𝑛)

Fracciones molares Libres de Sto

0.35 X0 0
0.0052 XN ?
L
yAn=0,259 xAn= ??

Grafico x vs y
1.20

1.00

0.80

0.60
y

0.40

0.20

0.00
0.00 0.10 0.20 0.30 0.40
x

Grafico X vs Y ( Fracciones Molares Libres de Sto)

1.00
0.90
0.80
0.70
0.60 Linea de
 Grafico X vs Y ( Fracciones Molares Libres de Sto)
1.00
0.90
0.80
0.70
0.60 Linea de
Recta d
0.50 Linear (
Y

0.40 Linea A

0.30 f(x) = 0.673158735486961 x + 0.0052069721155586

0.20
0.10
0.00
0.00 0.10 0.20 0.30 0.40 0.50 0.60 0.70 0.80 0.90 1.00
X

Grafico X vs Y
1.00
0.90
 Grafico X vs Y
1.00
0.90
0.80
0.70
0.60 Linea de Equilibrio
Recta de Operación
0.50 Linear (Recta de Oper
Y

0.40 Linea Auxiliar de Bott

0.30 f(x) = 1.00973810323044 x + 0.00520697211555862

0.20
0.10
0.00
0.00 0.10 0.20 0.30 0.40 0.50 0.60 0.70 0.80 0.90 1.00
X

Grafico X vs Y

Linea de Equilibrio
Recta de Operación
Linear (Recta de Operación )
810323044 x + 0.00520697211555862
Linea Auxiliar de Botton
aux 1
aux 2
aux 3
aux 4
aux 5
aux 6

0.20 0.30 0.40 0.50 0.60 0.70 0.80 0.90 1.00

X
 xA1=0

n+1

xAn= ??

0.40 0.50 0.60

es de Sto)

Linea de Equilibrio
es de Sto)

Linea de Equilibrio
Recta de Operación
Linear (Recta de Operación)
Linea Auxiliar de Botton

1.00
Linea de Equilibrio
Recta de Operación
Linear (Recta de Operación)
Linea Auxiliar de Botton
Caudal Volumetrico GR
560 m3/h Temperatura 20
y NH3 0.08 Presión Total 2
Caudal Vol. NH3 44.8 m3/h Constante R G.I. 0.082057
y Aire 0.92 MM NH3 17
Caudal Vol. Aire 515.2 m3/h MM medio Aire 29

Caudal másico NH3 928.14 kg/h Caudal molar NH3 54.60

Caudal másico Aire 18207.83 kg/h Caudal molar Aire 627.86
Caudal másico GR 19135.97 kg/h Caudal molar GR 682.45

MM Agua 18 kg/kmol

X (kmol Y (kmol y (kmol

x (kmol NH3/kmol
kg NH3 / 100 kg H2O p NH3 (mmHg) NH3/kmol NH3/kmol NH3/kmol
tot)
H2O) H2O) tot)
60 945.0 0.6353 1.6435 0.3885 0.6217
50 686.0 0.5294 0.8225 0.3462 0.4513
40 470.0 0.4235 0.4476 0.2975 0.3092
30 298.0 0.3176 0.2439 0.2411 0.1961
25 227.0 0.2647 0.1756 0.2093 0.1493
20 166.0 0.2118 0.1226 0.1748 0.1092
15 114.0 0.1588 0.0811 0.1371 0.0750
10 69.6 0.1059 0.0480 0.0957 0.0458
7.5 50.0 0.0794 0.0340 0.0736 0.0329
5 31.7 0.0529 0.0213 0.0503 0.0209
4 24.9 0.0424 0.0167 0.0406 0.0164
3 18.2 0.0318 0.0121 0.0308 0.0120
2 12.0 0.0212 0.0080 0.0207 0.0079

Bottom Top
yN+1 NH3 0.08 y1 0.008
YN+1 NH3 0.0870 Y1 0.0087
XN NH3 0.16 X0 0
Valor Base 0.09
Incremento 70
xN 0.138
PRENDIENTE DE
LA RECTA
 PRENDIENTE DE
(L'/G')min 0.489 LA RECTA
G' (molar) 627.86 kmol/h
L' min (molar) 307.02 kmol/h
Caudal másico AMS min. 5526.39 kg/h
°C
atm
atm m3/ kmol °C
kg/kmol
kg/kmol

kmol/h
kmol/h
kmol/h

Curva de Equilibrio con fracs. molares libres de sto.

1.8

1.6

1.4

1.2

1.0
Y

0.8

0.6

0.4

0.2

0.0 f(x) = 0.489130434782609 x + 0.00869565217391305

0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7

X
Curva de Equilibrio Línea de Operación
Linear (Línea de Operación)
 0.8625

0.7

de Operación
Caudal Másic 5264.2 10000 20000 30000 50000 60000
N° de etapas Fracción molar
2 0.07012 0.06377 0.05266 0.04338 0.02944 0.02444
3 0.06959 0.06290 0.05106 0.04043 0.02355 0.01767
4 0.06930 0.06246 0.05053 0.03926 0.02019 0.01365
5 0.06911 0.06218 0.05034 0.03879 0.01809 0.01098
6 0.06898 0.06199 0.05026 0.03860 0.01671 0.00909
7 0.06889 0.06185 0.05023 0.03852 0.01576 0.00767
8 0.06883 0.06174 0.05021 0.03849 0.01511 0.00658
9 0.06878 0.06165 0.05021 0.03848 0.01466 0.00571
10 0.06875 0.06158 0.05020 0.03847 0.01434 0.00500
11 0.01412 0.00442
12 0.01397 0.00393
13 0.01386 0.00352
14 0.01379 0.00316
15 0.01374 0.00286

Caudal Másic 10000 20000 30000 50000 52000

N° de etapas Fraccion Molar
20 0.06129 0.05020 0.03847 0.01365 0.01110
30 0.06120 0.05020 0.03847 0.01364 0.01106
40 0.06116 0.05020 0.03847 0.01364 0.01106
50 0.06114 0.05020 0.03847 0.01364 0.01106
60
70

y1 0.008 y2 0.008 Puntos Recta Auxiliar (Valor Objetiv

x1 0 x2 16

Fracs. molares vs n° de etapas para distintos caudale

0.03

0.025

0.02

0.015

0.01
0.015

0.01

0.005

0
0 2 4 6 8 10 12 14 16
 55000 57500 56500 56000
molar
0.02680 0.02559 0.02606 0.02631
0.02040 0.01899 0.01954 0.01982
0.01665 0.01508 0.01569 0.01600
0.01420 0.01250 0.01316 0.01350
0.01249 0.01068 0.01138 0.01174
0.01124 0.00933 0.01007 0.01045
0.01031 0.00830 0.00907 0.00947
0.00960 0.00749 0.00830 0.00872
0.00906 0.00684 0.00768 0.00813
0.00863 0.00631 0.00719 0.00766
0.00829 0.00588 0.00680 0.00728
0.00803 0.00552 0.00648 0.00698
0.00783 0.00523 0.00622 0.00674
0.00767 0.00499 0.00601 0.00655

Se trata de buscar un caudal de AMS no muy elevado junto con

54000 un razonable número de etapas. Ambas variables claro que
deben cumplir la condición impuesta de fraccion molar de
amoniaco en el gas de salida con un valor igual o menor a
0.00855 0,008.
0.00847 En los cuadros se resaltan tales valores, correspondientes a
0.00847 cierto caudal másico de líquido de entrada, y para cierto
0.00847 número de etapas.
Se puede concluir a partir de los resultados que una opción a
0.00847 elegir es la de un caudal de AMS de 56 mil kg/h, que permite
0.00847 llegar a valores de y1 menores al impuesto apartir de la etapa
número 11. Esto también se puede visualizar en el gráfico que
Recta Auxiliar (Valor Objetivo) está debajo de este texto.

para distintos caudales másicos

55000
56000
56500
57500
Valor objetivo
60000
 55000
56000
56500
57500
Valor objetivo
60000

12 14 16 18
 Amoniaco - Aire
Entra 3 %Vol
Sale 0.05 %Vol

G (en
condiciones 500 Nm3/h
normales) G2 = (P1*G1*T2)/(T1*P2)
G 554.91 m3/h

30 °C 0 °C
T2 T1
303.15 K 273.15 K
1 atm 1 atm
P2 P1
760 mmHg 760 mmHg

Relacion de Equilibrio para soluciones acuosas de Amoniaco

y=0.746x
H 0.746

x y X Y
0 0 0 0
0.05 0.04 0.05 0.04
0.1 0.07 0.11 0.08
0.15 0.11 0.18 0.13
0.2 0.15 0.25 0.18
0.25 0.19 0.33 0.23
0.3 0.22 0.43 0.29
0.35 0.26 0.54 0.35
0.4 0.30 0.67 0.43
0.45 0.34 0.82 0.51
0.5 0.37 1.00 0.59
0.55 0.41 1.22 0.70
 LUPA

Curva Auxiliar 1
Y X
0.030 0.00000
0.030 0.05000

Curva Auxiliar 2 5
Y X 125
0.000 0.040
0.050 0.040

b) 60 % más del caudal mínimo de agua

L'/G' min 0.74

L'/G' 1.18
LUPA

Curva Auxiliar 1
Y X
0.030 0.00000
0.030 0.05000

Curva Auxiliar 2
Y X
0.000 0.025
0.050 0.025

Metodo Grafico

Numero de Etapas 5 Etapas

 aprox 5

Metodo Analitico

L'/G' 1.18
H 0.574 Para encontrar H
A 2.056 realizo un ajuste
linea en la ec de
equilibrio
Y 0.017
n teóricas 4.779
n red 5

Por el gráfico

Y 0.017
n teóricas 4.9
n red 5
 Y vs X
0.8
0.7 f(x) = 0.574494964668233 x + 0.024764942692788
0.6 R² = 0.994226727346802
0.5
0.4
Y

0.3
0.2
0.1
0
0 0.2 0.4 0.6 0.8 1 1.2 1.4
X
 Y vs X
0.06 Top
0.05 X0
Y1
0.04
0.03
Bottom
Y

0.02 XN
0.01 YN+1
0
0 0.01 0.02 0.03 0.04 0.05 0.06
L'/G' min
X L' min

Curva de equilibrio Recta X

Recta Y Recta de operación

Y vs X Top
0.06 X0
0.05 Y1
0.04
0.03
Bottom
Y

0.02
XN
0.01 YN+1
0
0 0.01 0.02 0.03 0.04 0.05 0.06
L'/G' min
X
L' min

Curva de equilibrio Recta X Recta Y

Recta de operación Recta de op Real

Y vs X
0.06
0.05
0.04
0.03
Y

0.02
0.01
 0.06
0.05
0.04
0.03

Y
0.02
0.01
0
0 0.01 0.02 0.03 0.04 0.05 0.06
X

Curva de equilibrio Recta X Recta Y

Recta de operación Recta de op Real
 Top
0
0.0005

Bottom
0.04 2
0.03 50

0.7375
409.249725 m3/h

Top
0
0.0005

Bottom
0.025
0.03

1.18
0.352112 m3/h
 Datos y Primeros Cálculos

Caudal Volumétrico GR 0.8 m3/s Temperatura 30

800 L/s 303.15
Presión Total 1
y N+1 0.07 750.06
Se supondrá que el resto del gas es Aire Constante R de G.I. 0.08314
Caudal volumétrico SO2 0.056 m3/s MM SO2 64
Caudal volumétrico Aire 0.744 m3/s MM Aire. 29

Caudal Másico SO2 0.14 kg/s Caudal Másico GR 1.00

Caudal Másico Aire 0.86 kg/s 3594
Se supondrá absorcion de 95% de SO2 Caudal Molar GR 0.03
Fracción de Absorción 0.95 114.27
y1 0.00350

Caudal Másico liquido lavador 3.8 kg/s Densidad LP 1235

13680 kg/h Viscosidad LP 2.50E-03
MM H2O 18 kg/kmol

Línea de Equilibrio y de Caudal Mínimo

g SO2/100 g H2O p parcial (mmHg) X Y

0.01 0.12 2.81E-05 1.60E-04
0.05 1.68 1.41E-04 0.0022
0.1 4.62 2.81E-04 0.0062
0.15 8.07 4.22E-04 0.0109
0.2 11.8 5.63E-04 0.0160
0.25 15.7 7.03E-04 0.0214
0.3 19.8 8.44E-04 0.0271
0.4 28.3 0.0011 0.0392
0.5 37.1 0.0014 0.0520
1 83.7 0.0028 0.1256
 Pendiente recta Ley de Henry
X Y
0rigen 0 0 1.75E-03 Valor base
Top 2.08E-03 0.07527 33 Valor variable

36.1869
Constante de Henry aproximada

Línea Auxiliar
X Y
0 0.07527
2.08E-03 0.07527

Recta de Operación con Lmin/G

X Y
Top 0 0.00351
Bottom 0.00208 0.07527

Lmin/G 34.4983 kmol/kmol

Número de Etapas (suponiendo cumplimiento de Ley de Henry) y Línea de Operación

y1/y2 (yN+1/y1) 20
G' 106.27 kmol/h
Entro con Nog y
corto en la linea
NOG A = H/ L/G L' (kmol/h) vertical para leer a
3 0.10 38456 la derecha el valor
de A
4 0.39 9860
5 0.55 6992
6 0.69 5573
8 0.81 4748
10 0.89 4321
15 0.96 4006
19 1.00 3846
28 1.0489 3666

Caudal AMS vs N° de etapas

45000

40000

35000

30000

25000

20000

15000

10000
 35000

30000

25000

20000

15000

10000

5000

0
0 2 4 6 8 10 12 14 16 18 20 22 24 26 28 30

Caudal AMS vs N° de etapas

Se toman 8 etapas Ya que a partir de ese número la variación de caudal co

L' 4748 kmol/h

Recta de Operación con NOG = 8

X Y
Top 0 0.00351
Bottom 0.00161 0.07527

Líneas de Equilibrio y de CaudaL Mínimo

0.08

0.07 f(x) = 38.1068598706971 x − 0.00387363943558793

R² = 0.990511878969297
0.06

0.05

0.04

LINEA DE OP REAL
0.03

0.02

0.01

0.00
0.0E+00 5.0E-04 1.0E-03 1.5E-03 2.0E-03 2.5E-03

Estimación de Diámetro de Torre

Masa Molar ponderada 31.45 kg/kmol

Densidad GR 1.25 kg/m3
Densidad LP 1235 kg/m3
FLV 1.42
 DP Relleno Intalox Cerámica 25 1.22 mm H2O/ft
mm (aprox 1 in)
4.00 mm H2O/m

K4 Diseño 0.101
K4 Inundación 0.375
% inundación 52%

Fp (Intalox Cerámica 25 mm) 300 1/m

V*W 0.38 kg/s m2
g' (másico) 0.93 kg/s
Área de Flujo en zona rellena (AF) 2.42 m2

Diámetro de Torre Calculado 1.76 m

Diámetro Escogido 1.8 m
Área de flujo recalculada 2.54 m2
% inundación corregido 49%

Estimación de Altura de Torre (HOG) con Método de Cornell

Volúmen de difusión Aire 20.1

Volúmen de difusión SO2 41.1
Difusividad GR 1.31E-05 m2/s

f 2.6
Volumen molar Agua a Teb 0.0189
Difusividad L 1.05E-06 m2/s

Viscosidad GR 1.81E-05 Pa s

Sc G 1.10
Sc L 1.92

K3 0.97

Y 60
L* = L / Ac 1.49 kg/s m2

F n 3.65E-02

Viscosidad del Agua a 30°C 8.20E-04 Pa s

Densidad del Agua a 30°C 992.81 kg/m3

f1 1.20
f2 1.31
f3 (se considerará 1) 1

Altura por etapa 0.5 m

Altura Inicial (HOG * NOG) 4m

Z iterativo 74.8389571875 m

HL 2.42E-02 m

HG 9.33526634271 m

HOG Calculado 9.35E+00 m

Z calculado 7.48E+01 m
Error Relativo Porcentual 0.00%

Resultados a informar

a) Altura Redondeada 75 m

b) Diámetro Escogido 1.8 m

Puede que el gran valor de altura obtenido se deba a que no se consideró completamente que
el líquido AMS contiene una sal de amonio. Esto se ve por ejemplo cuando se traza la línea de
equilibrio, relacionando el SO2 sólo con agua. También se ve en que no se tuvo del todo en cuenta en
el cálculo de la Difusividad de líquido, ni en el f3 (relación de tensiones superficiales).
 G1; y1

°C
K
bar
mmHg
L0; x0
m3 bar/kmol K
kg/kmol
kg/kmol

kg/s
kg/h
kmol/s
G N+1; Y N+1
kmol/h

kg/m3
kg/m s

LN; xN

Líneas de Equilibrio y d
0.08
Fracción molar libre de soluto
en fase líquida (X)
0.07 f(x) = 38.1068598706971 x − 0.00387363943558793
R² = 0.990511878969297
Fracción molar libre de soluto
en fase gaseosa (Y) 0.06

0.05
Ley de Dalton

0.04

0.03

0.02

0.01
 0.04

0.03

0.02

0.01

0.00
0.0E+00 5.0E-04 1.0E-03 1.5E-03

Lupa

0.07

0.05
1.0E-03 1.2E-03 1.4E-03 1.6E-03 1.8E-03 2.0E-03 2.2E-03 2.4E-03
Operación

con Nog y
en la linea
al para leer a
echa el valor
 26 28 30

e número la variación de caudal con respecto a la de etapas es insignificante.

LINEA DE OP REAL
v
s de Equilibrio y de CaudaL Mínimo

71 x − 0.00387363943558793
97
1.0E-03 1.5E-03 2.0E-03 2.5E-03
Eliminar 99.9 % de Ac sulfidrico
G1 Composicion
Molar %(y)
899.60 kg/h 4%
46.31 kmol/h 80%
16%
0.006%
100%

L0

42997.2
1993

GN+1

1000
49.26

LN Composicion
Molar x
43097.55 kg/h 0.00
1995.72 kmol/h 0.04
0.96
1
 Datos de equilibrio Aire/Sulfuro de hidrógeno SRK

Curva de equilibrio
x y X Y
0 0 0 0
0.0003 0.0103 0.0003 0.0104
0.0007 0.0203 0.0007 0.0208
0.0010 0.0301 0.0010 0.0311
0.0013 0.0397 0.0013 0.0414
0.0017 0.0491 0.0017 0.0516
0.0020 0.0582 0.0020 0.0618
0.0023 0.0671 0.0023 0.0719
0.0027 0.0758 0.0027 0.0820
0.0030 0.0843 0.0030 0.0921
0.0033 0.0926 0.0033 0.1021
0.0037 0.1008 0.0037 0.1121
0.0040 0.1087 0.0040 0.1220
0.0043 0.1165 0.0044 0.1319
0.0047 0.1241 0.0047 0.1417
0.0050 0.1316 0.0050 0.1515
0.0053 0.1389 0.0054 0.1613
0.0057 0.1460 0.0057 0.1710
0.0060 0.1530 0.0060 0.1807
0.0063 0.1599 0.0064 0.1903
0.0067 0.1666 0.0067 0.1999
0.0070 0.1732 0.0070 0.2094
0.0073 0.1796 0.0074 0.2189
0.0077 0.1859 0.0077 0.2284
0.0080 0.1921 0.0081 0.2378
0.0083 0.1982 0.0084 0.2472
0.0087 0.2042 0.0087 0.2566
0.0090 0.2100 0.0091 0.2659
0.0093 0.2158 0.0094 0.2751
0.0097 0.2214 0.0098 0.2843
0.0100 0.2269 0.0101 0.2935

Pendiente recta Ley de Henry

X Y
0rigen 0 0 4.00E-12 Valor base
Bottom 2.06E-03 6.38E-02 206 Valor variable

Constante de Henry aproximada 30.9853

Línea Auxiliar
X Y
0 0.06383
2.06E-03 0.06383

Número de Etapas (suponiendo cumplimiento de Ley de Henry) y Línea de Operación

y1/y2 (yN+1/y1) 1000.00

G 49.26 kmol/h
G' 46.31 Kmol/h

NOG A = H/ (L'/G') L' (kmol/h)

7 0.10 14348
8 0.29 4948
10 0.38 3776
12 0.49 2928
14 0.56 2562
15 0.60 2391
19 0.70 2050
20 0.72 1993
22 0.75 1913
26 0.79 1816

Caudal AMS vs N° de etapas

16000

14000

12000
Caudal AMS (Kmol/h)

10000

8000

6000

4000

2000
 Caudal AMS (Kmol/
10000

8000

6000

4000

2000

0
6 8 10 12 14 16 18 20 22 24 26 28

N° de etapas

Caudal AMS vs N° de etapas

Se toman 20 etapas Ya que a partir de ese número la variación de cau
L' 1993 kmol/h
L 1993 Kmol/h

Recta de Operación con NOG = 20

X Y
Top 0 0.00006
Bottom 0.00158 0.06383

Estimación de Diámetro de Torre

Masa Molar ponderada 20.3 kg/kmol

Densidad GR 8.11 kg/m3
Densidad LP 1013 kg/m3
Viscosidad LP 1.07E-03 Kg/(m*s)
FLV 3.85
DP Relleno Intalox Cerámica 38 mm 1.09 mm H2O/ft
(1.5 in) 3.58 mm H2O/m

K4 Diseño 0.070
K4 Inundación 0.130
% inundación 73%

Fp (Intalox Cerámica 38 mm) 170 1/m

V*W 1.01 kg/s m2
g' (másico) 0.25 kg/s
Área de Flujo en zona rellena (AF) 0.25 m2

Diámetro de Torre Calculado 0.56 m

Diámetro Escogido 1.1 m
Área de flujo recalculada 0.95 m2
% inundación corregido 19%
 Estimación de Altura de Torre (HOG) con Método de Cornell

H2S CH4 C2H6 CO2

σ (A) 3.62 3.758 4.443 3.941
ε/κ (K) 301.10 148.6 215.7 195.2
y 0.06 0.75 0.15 0.04
σAB 3.6905 - 4.033 3.782
εAB/κ (K) 211.53 - 254.85 242.43
T* 1.43 - 1.19 1.25
ΩD 1.2163 - 1.3212 1.2971
PM 34 16 30 44
DAB (cm2/s) 1.45E-01 - 1.12E-01 1.30E-01
y' - 0.79787234 0.159574468085 0.042553191489

1.38E-01 cm2/s
Difusividad de H2S en la GR
1.38E-05 m2/s

Estimación de difusividad en el líquido: Método de WILKE-CHANG

f 2.6
Volumen molar H2S a Teb 0.0329
Difusividad L 1.93E-06 m2/s
Viscosidad GR 1.12E-05 Pa s

Sc G 0.0995
Sc L 0.54

K3 0.8

Y 81

Lw* = Lw/ Ac 12.57 kg/s m2

F n 1.20E-01

Viscosidad del Agua a 30°C 8.20E-04 Pa s

Densidad del Agua a 30°C 992.81 kg/m3
Tensión superficial del agua a 30 °C 7.12E-02 N/m
Tensión superficial LP 6.66E-02 N/m
f1 1.04
f2 1.03
f3 0.948

Altura por etapa 0.5 m

Altura Inicial (HOG * NOG) 10 m

Z iterativo 1.83 m

HL 0.020 m

HG 0.2761800925 m

HOG Calculado 0.29 m

Z calculado 5.8118980785 m
Error Relativo Porcentual -2.175900589

Resultados a informar

a) Altura Redondeada 5.81 m

b) Diámetro Escogido 1.1 m

Relacion de esbeltez 5.2836363636

 Composicion
Molar (y) Gw(kg/h) G(kmol/h) Pesos Moleculares (Kg/Kmol)
0.0425 86.70 1.97 44 CO2
0.7978 591.13 36.95 16 CH4
0.1596 221.67 7.39 30 C2H6
6E-05 0.09 0.0028 34 H2S
1 899.60 46.31 19.43 Peso Molecular Promedio

Composicion Composición Pesos Moleculares

L0
Molar x másica x (w) Lw(kg/h) L(kmol/h) (Kg/Kmol)
kg/h 0 0 0 0 34
kmol/h 0.04 0.20 8599.71 81.79 105.14
0.96 0.80 34397.45 1910.97 18
1 1 42997.2 1993 21.58

GN+1 Composicion Composicion Pesos Moleculares

Molar %(y) Molar (y) Gw(kg/h) G(kmol/h) (Kg/Kmol)
kg/h 4% 0.04 86.70 1.97 44
kmol/h 75% 0.75 591.13 36.95 16
15% 0.15 221.67 7.39 30
6% 0.06 100.49 2.96 34
100% 1 1000 49.26 20.3

Lw(kg/h) L(kmol/h) Pesos Moleculares (Kg/Kmol)

100.40 2.95 34 H2S
8599.71 81.79 105.14 Dietanolamina
34397.45 1910.97 18 H2O
43097.55 1995.72 21.60 Peso molecular promedio

Constante R de G 0.08314 m3 bar/kmol K

P 9.8 bar
30 °C
T
 T
303.15 K

Sulfuro de hidrógeno SRK

Y vs X
0.35
0.3
f(x) = 29.0522128540612 x + 0.00356940601992928
0.25 R² = 0.999603530402166
0.2
0.15
Y

0.1
0.05
0
0 0.002 0.004 0.006 0.008 0.01 0.012
X

Curva de equilibrio Linear (Curva de equilibrio)

Linear (Curva de equilibrio)

Y vs X
0.07
0.06
0.05 f(x) = 30.9121065411852 x
R² = 0.999990849906106
0.04
0.03
Y

0.02
 Y vs X
0.07
0.06
0.05 f(x) = 30.9121065411852 x
R² = 0.999990849906106
0.04
0.03

Y
0.02
0.01
0
0 0.0005 0.001 0.0015 0.002 0.0025
X

Column E Linear (Column E)

Pendiente recta Ley de Henry Línea Auxiliar

Operación
 24 26 28

de ese número la variación de caudal con respecto a la de etapas es insignificante.

v
 El componente mayoritario es el metano. Se calcula la difusividad del H2S en el metano, etano
30 °C
T
303.15 K
Pesos Moleculares
(Kg/Kmol)
H2S
Di Etanol Amina (DNA)
H2O
Peso Molecular Promedio

Pesos Moleculares
(Kg/Kmol)
CO2
CH4
C2H6
H2S
Peso Molecular promedio
0.01 0.012

brio)
ad del H2S en el metano, etano y CO2
 DATOS
30 °C
T
303 K
1 atm
P
760 mmHg
ya 8%
A a eliminar 75.0%
Gs 𝑉 ̇
50 m3/min
3000 m3/h
𝜌_𝑙𝑖𝑞𝑢𝑖𝑑𝑜 996 kg/m3
996 g/L Evaluadas a
𝜇_𝑙𝑖𝑞𝑢𝑖𝑑𝑜 0.000798 kg/sm 30°C

𝑀𝑀 𝐶𝑙_2
0.798 mN/sm

𝑀𝑀 𝐻_2 𝑂
35.45 g/mol

𝑀𝑀 𝑎𝑖𝑟𝑒
18 g/mol
29 g/mol
R 0.082057 m3atm/kmolK

𝑃_𝐴 (mmHg) ya Solubilidad

g H2O mol H2O g Cl2
Solubilidad (g Cl2 / L (g Cl2 / g
H2O) H2O)
 Solubilidad
Solubilidad (g Cl2 / L (mmHg) (g Cl2 / g g H2O mol H2O g Cl2
H2O) H2O)
0.424 5 0.007 0.0004257 1.000 0.0556 0.000425703
0.553 10 0.013 0.00055522 1.000 0.0556 0.000555221
0.873 30 0.039 0.00087651 1.000 0.0556 0.000876506
1.106 50 0.066 0.00111044 1.000 0.0556 0.001110442
1.573 100 0.132 0.00157932 1.000 0.0556 0.001579317
1.966 150 0.197 0.0019739 1.000 0.0556 0.001973896
2.34 200 0.263 0.0023494 1.000 0.0556 0.002349398
2.69 250 0.329 0.0027008 1.000 0.0556 0.002700803
3.03 300 0.395 0.00304217 1.000 0.0556 0.003042169

ya1 8%

𝑦_𝐴2=𝑦_𝐴1 (1−𝑅)
YA1 0.08695652
ya2 0.0200 0.16
YA2 0.0204 f(x) = 81359226.973951 x³ +
0.14
R² = 0.999937366536379
0.12
Realizamos un balance de masa
∆𝑊 ̇=𝑊 ̇^𝑖
0.10
0.08

−𝑊 ̇_𝐴^𝑖𝐺=𝑊 ̇_𝐴
𝑊 ̇_𝐴^𝑖𝐿−𝑊 ̇_𝐴^𝑖𝐺=−𝑊 ̇_𝐴2^𝐿 〖 +𝑊 ̇ 〗 _𝐴1^𝐿+𝑊 ̇_𝐴2^𝐺−𝑊 ̇_𝐴1^𝐺=0

Y
2^𝐺−𝑊 ̇_𝐴1^𝐺
0.06

(𝑊 ̇_𝐴2^𝐺)/𝐴_𝑡 −(𝑊 ̇_𝐴1^𝐺)/𝐴_𝑡 =(𝑊 ̇_𝐴2^𝐿)/𝐴_𝑡 −(𝑊 ̇_𝐴1^𝐿)/𝐴_𝑡

0.04

𝑊 ̇_𝐴^𝑖𝐿=𝑊 ̇_𝐴1^𝐿
0.02

−𝑊 ̇_𝐴2^𝐿
0.00

𝐺_2 𝑦_𝐴2−𝐺_1 𝑦_𝐴1=𝐿_2

-0.0002 0.0000

𝑥_𝐴2−𝐿_1 𝑥_𝐴1
-0.02

𝐺𝑦_𝐴=𝐺_𝑠 𝑌_𝐴 𝐺_𝑠 𝑌_𝐴2−𝐺_𝑠 𝑌_𝐴1=𝐿_𝑠 𝑋_𝐴2−𝐿_𝑠 𝑋_𝐴1

𝐿𝑥_𝐴=𝐿_𝑠 𝑋_𝐴 𝐺_𝑠 (𝑌_𝐴2−𝑌_𝐴1 )=𝐿_𝑠 (𝑋_𝐴2−𝑋_𝐴1 )
Curva de Equilibrio
Recta Mínima de Operaci

(𝑌_𝐴2−𝑌_𝐴1)/(𝑋_𝐴2−𝑋_𝐴1 )=𝐿_𝑠/𝐺_𝑠 =𝑝𝑒𝑛𝑑𝑖𝑒𝑛𝑡𝑒 Del grafico:

A
B
C
Calculo la pendiente de la recta mínima de operación: D
Pendiente 107.17 (L'/G')min
Recta de Top
Calculo el flujo molar: XA
𝑃𝑉 ̇=𝑛 ̇𝑅𝑇 ⟹ 𝑛 ̇=(𝑃𝑉 ̇)/𝑅𝑇 0
0.0006

Gs 120.66 kmol/h
Ls min 12931.47 kmol/h Recta Ley de Henry
XA YA
0 0
0.000618 0.0870
618
 Calculo la pendiente (H) de la recta de Ley de Henry:
H 140.71

Método Analítico para determinar el número de etapas:

𝐴=(𝑚𝐺_𝑚)/𝐿_𝑚 ⟹ 𝐹𝑎𝑐𝑡𝑜𝑟
𝑑𝑒 𝑎𝑏𝑠𝑜𝑟𝑐𝑖ó𝑛
ya1/ya2 4

Del gráfico:
NOG A Ls' (kmol/h) m vendria siendo H
2 0.59 28775.55 Constante de henry
3 1 16977.58
4 1.15 14763.11
5 1.2 14147.98
6 1.3 13059.67 es el AMS
7 1.31 12959.98
8 1.32 12861.80
9 1.33 12765.09
10 1.33 12765.09
11 1.33 12765.09

Caudal de AMS vs Número de Etapas

35000

30000

25000
Caudal Ls (kmol/h)

20000

15000

10000

5000

0
0 2 4 6 8 10 12
Número de Etapas

Etapas Teóricas 4
Este valor es razonable debido a que
14763.11 kmol/h para mayores etapas, no varía
Ls'
4.10 kmol/s notablemente el caudal de AMS.

Eficiencia 75% Elijo este valor, segun mi criterio.

 Elijo este valor, segun mi criterio.
5.33
Etapas Reales 6.00
13059.67 kmol/h
Ls' nuevo El Ls que corresponde para ese numero de etapas
3.63 kmol/s

DIÁMETRO DE LA COLUMNA
Relleno Pall 1''
𝑀𝑀=𝑚 ̇/𝑛 ̇ ⟹ 𝑚 ̇=𝑛 ̇𝑀𝑀_𝑝𝑟𝑜𝑚𝑒𝑑𝑖𝑜
120.66 kmol/h
3561 kg/h
Gs
0.0335 kmol/s
0.99 kg/s
𝑀𝑀=𝑚 ̇/𝑛 ̇ ⟹ 𝑚 ̇=𝑛 ̇𝑀𝑀
13059.67 kmol/h
235074.13 kg/h
Ls'
3.63 kmol/s
65.30 kg/s

De Tabla 11.2:
Factor de relleno Fp 160 1/m

𝜌_(𝐶𝑙_
Sacar las densidades de cada uno

2)
evaluados a la T° correspondiente
𝜌_𝑎𝑖𝑟𝑒
2.85 kg/m3 Evaluadas a

𝜌_𝑔𝑎𝑠
1.164 kg/m3 30°C

𝜌𝑔𝑎𝑠=𝜌𝐶𝑙2∗𝑦𝐴1+𝜌𝑎𝑖𝑟𝑒∗(1−𝑦𝐴1)
1.30 kg/m3

1. debo calc
2. Luego a p
. Los valores
desorción 1
Elijo diseñar
(correspond
3.Se busca la
4. Realizo lo
 𝐾_4=(13.1𝐺_𝑠^(′ 2) 𝐹_𝑝 (𝜇_𝑙𝑖𝑞𝑢𝑖𝑑𝑜/𝜌_𝑙𝑖𝑞𝑢𝑖𝑑𝑜 )^0.1)/(𝜌_𝑔𝑎𝑠 (𝜌_𝑙𝑖𝑞𝑢𝑖𝑑𝑜−𝜌_𝑔

𝐹_𝐿𝑉=(𝐿_𝑠^′)/𝐺_𝑠
FLV 2.38 √(𝜌_𝑔𝑎𝑠/𝜌_𝑙𝑖𝑞𝑢𝑖𝑑𝑜 )

Caída de presión por metro de relleno 21 mmH2O/m

Elijo diseñar para una caída de presión de 21 mm de H2O por metro de relle
Del gráfico:
K4 Diseño 0.14

% 𝑖𝑛𝑢𝑛𝑑𝑎𝑐𝑖ó𝑛=√((𝐾_4 𝑑)/(𝐾_4 𝑖))∗100

K4 Inundación 0.24

% Inundación 76%

Área columna requerida 𝐴𝑟𝑒𝑎=𝐺_𝑠/(𝐺_𝑠^′

1.67 m2 )=(𝐶𝑎𝑢𝑑𝑎𝑙 𝑚𝑜𝑙𝑎𝑟 𝑑𝑒 𝑔𝑎𝑠)/(𝐹𝑙𝑢𝑗𝑜 𝑚𝑜𝑙𝑎𝑟 𝑑𝑒 𝑔

𝐷_𝑐=√(4𝐴𝑟𝑒𝑎/𝜋)
Dc 1.46 m
Dc redond. 1.5 m

Área de columna nueva 1.77 m2

% 𝑖𝑛𝑢𝑛𝑑𝑎𝑐𝑖ó𝑛_𝑛𝑢𝑒𝑣𝑜=%𝑖𝑛𝑢𝑛𝑑𝑎𝑐𝑖ó𝑛∗𝐴_𝑟𝑒𝑞𝑢𝑒𝑟𝑖𝑑𝑎/𝐴_𝑛𝑢𝑒𝑣𝑎
% Inundación nuevo 72%

ESTIMACIÓN DE LA ALTURA - Método de Cornell

Calculamos los coeficientes de difusión
 De Tabla 8.5:
Cl2 𝐷_𝑔𝑎𝑠=(1.013∗ 〖37.7
10 〗 ^(−7) 𝑇^1.75 (1/(𝑀𝑀_𝑎 )+1/(𝑀𝑀_𝑏 ))^0.5)/(𝑃(𝜈_𝑎^(1/3)+𝜈
Aire 20.1

Dgas 1.51E-05 m2/s

De Tabla 8.6:
Vm Cl2 0.0484 m3/kmol

φ 𝐷_𝑙𝑖𝑞𝑢𝑖𝑑𝑜=(1.173∗
2.6 〖 10 〗 ^(−13) (𝜙𝑀𝑀)^0.5 𝑇)/(𝜇𝑉_𝑚^0.6 )
Dliquido 1.87E-09 m2/s

𝜇_(𝐶𝑙_2 )
𝜇_𝑎𝑖𝑟𝑒 1.86E-05 kg/sm𝜇_𝑔𝑎𝑠=(𝑦_(
30°C〖𝐶𝑙〗 _2 )/𝜇_( 〖𝐶𝑙〗 _2 ) +𝑦_𝑎𝑖𝑟𝑒/𝜇_𝑎𝑖𝑟𝑒 )^(−1)
1.36E-05 kg/sm Evaluadas a

𝜇_𝑔𝑎𝑠 1.81E-05 kg/sm

𝑆𝑐_𝐺 𝑆𝑐_𝐿=𝜇_𝑙𝑖𝑞𝑢𝑖𝑑𝑜/(𝜌_𝑙𝑖𝑞𝑢𝑖𝑑𝑜
𝑆𝑐_𝐺=𝜇_𝑔𝑎𝑠/(𝜌_𝑔𝑎𝑠 𝐷_𝑔𝑎𝑠 ) 𝐷_𝑙𝑖𝑞𝑢𝑖𝑑𝑜 )
𝑆𝑐_𝐿
0.92
427.36
𝐿_𝑠^′′=(𝐿_𝑠^′)/(𝐴𝑟𝑒𝑎 𝑑𝑒 𝑐𝑜𝑙𝑢𝑚𝑛𝑎 𝑛𝑢𝑒𝑣𝑎)
𝐿_𝑠^′′ 36.95 kg/sm2

De la Figura 11.51:
K3 0.71

De la Figura 11.52:
ψh 57

De la Figura 11.53:
φh 0.15
 Suponiendo:

𝑍=𝑎𝑙𝑡𝑢𝑟𝑎 𝑑𝑒 𝑙𝑎 𝑐𝑜𝑙𝑢𝑚𝑛𝑎=𝑛𝑟𝑜 𝑒𝑡𝑎𝑝𝑎𝑠∗𝐻_𝑂𝐺

HOG' 1.332 m 1332
Z 7.992 m
𝐻_𝐿=0.305𝜙_ℎ 𝑆𝑐_𝐿^0.5 𝐾_3 (𝑍/3.05)^0.15

𝑓_1=(𝜇_𝑙𝑖𝑞𝑢𝑖𝑑𝑜
HL 0.7759 m

20 °C f1 0.9642
𝑓_2=(𝜌_𝑙𝑖𝑞𝑢𝑖𝑑𝑜
Tw
𝜇_𝑤
293 K

𝜌_𝑤
0.001002 kg/sm f2 0.9975

𝑓_3=(𝜎_𝑤/𝜎_𝑙𝑖
Evaluadas a
𝜎_𝑤
998 kg/m3 20°C
0.0728 kg/s2 f3 1.0179

𝜎_𝑙𝑖𝑞𝑢𝑖𝑑𝑜 0.0712 kg/s2 Ev. a 30°C

𝐻_𝐺=(0.011Ψ_ℎ 𝑆𝑐_𝑔^0.5
(𝐷_𝑐/0.305)^1.11
(𝑍/3.05)^0.33)/(𝐿_𝑠^′ 𝑓_1 𝑓_2
HG 0.3158 m
𝑓_3 )^0.5
HOG 1.3245 m 𝐻_𝑂𝐺=𝐻_𝐺+𝐴𝐻_𝐿

Z calculado 7.9470 m
Error 0.56%

RESUMEN DE DATOS
Etapas 6.00
Altura 8.0 m
Diametro 1.5 m
𝑥𝐴=𝑥𝐶𝑙2=(𝑚𝐶𝑙2/𝑚𝑚𝐶𝑙2)/(𝑚𝐶𝑙2/
𝑚𝑚𝐶𝑙2+1𝑔/𝑚𝑚𝐻20)=(𝑀𝑜𝑙𝑒𝑠 𝑑𝑒
𝐶𝑙2)/(𝑀𝑜𝑙𝑒𝑠 𝑑𝑒 𝑠𝑜𝑙𝑢𝑐𝑖𝑜𝑛)

Curva de Equilibrio
Solubilidad
mol Cl2 (mol Cl2/ xa
mol H2O) XA YA
 Solubilidad
mol Cl2 (mol Cl2/ xa
mol H2O) XA YA
1.20E-05 0.0002 0.0002 0.0002 0.0066
1.57E-05 0.0003 0.0003 0.0003 0.0133
2.47E-05 0.0004 0.0004 0.0004 0.0411
3.13E-05 0.0006 0.0006 0.0006 0.0704
4.46E-05 0.0008 0.0008 0.0008 0.1515
5.57E-05 0.0010 0.0010 0.0010 0.2459
6.63E-05 0.0012 0.0012 0.0012 0.3571
7.62E-05 0.0014 0.0014 0.0014 0.4902
8.58E-05 0.0015 0.0015 0.0015 0.6522

Grafico X vs Y
0.16
f(x) = 81359226.973951 x³ + 112443.018268824 x² + 62.7240607726128 x − 0.0144804781050036
0.14
R² = 0.999937366536379
0.12
0.10
0.08

𝑊 ̇_𝐴1^𝐺=0
0.06
0.04
0.02
0.00 Top = 2
0002 0.0000 0.0002 0.0004 Bottom0.0006
=1 0.0008 0.0010
-0.02

X
Curva de Equilibrio Polynomial (Curva de Equilibrio) Top
Recta Mínima de Operación Recta Ley de Henry

Con solver encuentro XA1

80000000 correspondiente a YA1 previamente seria XA2 el q encuentro es para el punto bottom
112443 calculado:
62.724 XA1 0.0006
-0.0145 YA1 0.08695632

Recta de Top Recta Mín. de Operación

YA XA YA
0.0870 0 0.0204 Top (2)
0.0870 0.0006 0.0870 Bottom (1)
ar el número de etapas:

Ingreso con yA1/yA2= 4 y trazo vertical. Lueg

vertical con distintos valores de NOG (arbrita
extrayendo los distintos valores de A (curvas
leo el valor en la derecha) para despues desp
AMS)

iterio.
iterio.

umero de etapas

1. debo calcular el FLV para ingresar por abcisas como indica la formula. Trazo vertical
2. Luego a partir de lo siguiente:
. Los valores de diseño recomendados, milímetros de agua por metro de relleno, son: Absorbedores y columnas de
desorción 15 a 50 (pag 803 Sinnot Towler)
Elijo diseñar para una caída de presión de 21 mm de H2O por metro de relleno.
(corresponde a la (21) la cual indica líneas de caída de presión constante (parametro))
3.Se busca la interseccion de la curva (21) con la vertical de FLV y en eje y leo K4
4. Realizo lo mismo pero esta vez con la curva de inundacion.
𝑔𝑎𝑠 (𝜌_𝑙𝑖𝑞𝑢𝑖𝑑𝑜−𝜌_𝑔𝑎𝑠)) ⟹ 𝐺_𝑠^′=√((𝐾_4 𝜌_𝑔𝑎𝑠 (𝜌_𝑙𝑖𝑞𝑢𝑖𝑑𝑜−𝜌_𝑔𝑎𝑠))/(13.1𝐹_𝑝 (𝜇_𝑙𝑖𝑞𝑢𝑖𝑑𝑜/𝜌_𝑙𝑖𝑞𝑢𝑖𝑑𝑜 )^0
K4= factor que relaciona los
*caudales del líquido y el
vapor
*las propiedades físicas del
sistema
*las características del relleno
*caudal másico del gas por

𝐺_𝑠^′
unidad de área transversal
0.5927 kg/sm2 A su vez en el grafico tengo
líneas de caída de presión
constante como un
parámetro
de H2O por metro de relleno.

Seguridad operativa y eficiencia

(𝐹𝑙𝑢𝑗𝑜 𝑚𝑜𝑙𝑎𝑟 𝑑𝑒 𝑔𝑎𝑠)=𝜋𝑟^2=(𝜋𝐷_𝑐^2)/4
Las torres de relleno se diseñan para operar por
debajo del punto de inundación, típicamente entre
60% y 80% del flujo de gas correspondiente a la
inundación.

En este caso el 76% es satisfactorio

𝑞𝑢𝑒𝑟𝑖𝑑𝑎/𝐴_𝑛𝑢𝑒𝑣𝑎
.5)/(𝑃(𝜈_𝑎^(1/3)+𝜈_𝑏^(1/3) )^2 )

𝜇𝑉_𝑚^0.6 )

𝑒/𝜇_𝑎𝑖𝑟𝑒 )^(−1)
𝑓_1=(𝜇_𝑙𝑖𝑞𝑢𝑖𝑑𝑜/𝜇_𝑤 )^0.16

𝑓_2=(𝜌_𝑙𝑖𝑞𝑢𝑖𝑑𝑜/𝜌_𝑤 )^1.25Podria haber considerado f1=f2=f3 por tratarse

de liquido como el agua a 20°C

𝑓_3=(𝜎_𝑤/𝜎_𝑙𝑖𝑞𝑢𝑖𝑑𝑜 )^0.8

(𝐷_𝑐/0.305)^1.11
Corrección 2.3
para el punto bottom
yA1/yA2= 4 y trazo vertical. Luego intersecto a esta
distintos valores de NOG (arbritarios) y voy
os distintos valores de A (curvas que voy cortando
n la derecha) para despues despejar LS (Caudal
y columnas de
(𝜇_𝑙𝑖𝑞𝑢𝑖𝑑𝑜/𝜌_𝑙𝑖𝑞𝑢𝑖𝑑𝑜 )^0.1 ))
aciona los
ido y el

físicas del

s del relleno
el gas por
nsversal
fico tengo
e presión
un
 CURVA DE EQUILIBRIO

% m/m SO2 Pi SO2 (mmHg) moles SO2

0.05 1.200 0.00078
0.1 3.200 0.00156
0.15 5.800 0.00234
0.2 8.500 0.00313
0.3 14.100 0.00469
0.5 26.000 0.00781
0.7 39.000 0.01094
1 59.000 0.01563
1.5 92.000 0.02344

303 K
T=
30 ºc
1 atm
Ptotal= 760 mmHg
101325 Pa

H= 32.81

CARACTERIZACIÓN DE LA CORRIENTE GASEO

GN+1 PM(kg/kmol) ρ(kg/m3)

SO2 64
Aire 29
Total 31.45 1.248

0.8 m3/s
Gvol 2880 m3/h
115914.030427433 mol/h
Gmolar 115.914030427433 kmol/h
DETERMINACIÓN DEL NÚMERO DE ETAPAS

y1= 0.07
y2= 7.00E-04
y1/y2= 100.00

NOG A
5 0.1
6 0.3
 7 0.42
8 0.53
9 0.6
10 0.65
11 0.68
12 0.73
13 0.75
14 0.77
15 0.79
16 0.82
17 0.84
18 0.85
19 0.855
20 0.87
21 0.88
22 0.89
23 0.895
24 0.9

Se toman
Nº etapas (gráfico) 10
L'= 5441.08 kmol/h

Se construye el gráfico de caudal de AMS L' vs Nº de etapas NOG. Se elige el número de

etapas hasta lograr el que caudal de AMS no varíe significativamente al seguir aumentando
el NOG

RECTA DE OPERACIÓN -(L'/G') Límite

CURVA DE EQUILIBRIO
X Y
0.0001 0.0016
0.0003 0.0042
0.0004 0.0077
0.0006 0.0113
0.0008 0.0189
0.0014 0.0354
0.0020 0.0541
0.0028 0.0842
0.0043 0.1377
 LO
G1

𝐺^′ 𝑌1+1+𝐿^′ 𝑋𝑜=𝐺^′

Del balance de materia:

𝐺^′ (𝑌1−𝑌𝑛+1)=𝐿^′ (𝑋𝑜−

Despejando

𝑋𝑛=𝑋𝑜−(𝑌1−𝑌𝑛+1)/(𝐿^′
Por lo tanto:

TORRE DE ABSORCIÓN

LN GN+1

BALANCE DE MASA
 El relleno Intalox de cerámica de 25 mm se selecciona en este caso porque tiene
características que lo hacen adecuado para la absorción de gases en torres empacadas,
especialmente en la eliminación de dióxido de azufre (SO₂) mediante lavado con una
solución acuosa. Algunas razones clave para su elección incluyen:
Alta eficiencia de contacto gas-líquido:
Su diseño estructurado maximiza la superficie de contacto entre el gas y el líquido, lo
que mejora la transferencia de masa y aumenta la eficiencia de absorción del SO₂.
Baja caída de presión:
La forma del Intalox permite un flujo de gas eficiente con menor resistencia al paso del
gas, lo que reduce el consumo de energía.
Resistencia química y térmica:
Al ser de cerámica, el material es altamente resistente a la corrosión química causada
por el gas ácido (SO₂) y el lavado con soluciones acuosas.
Buena capacidad de mojado:
La cerámica tiene una alta afinidad por los líquidos, asegurando una distribución
uniforme de la solución de lavado sobre la superficie del relleno.
Tamaño de 25 mm adecuado para este tipo de operación:
Un tamaño intermedio como 25 mm equilibra eficiencia de contacto y resistencia
hidráulica, siendo ideal para torres de absorción de gases a caudales moderados.
En resumen, este tipo de relleno es una excelente opción para procesos de absorción de
gases debido a su eficiencia, durabilidad y bajo costo operativo.

VERIFICACION DEL BALANCE DE MASA

Entrada
Gn+1 3645.50 kg/h
Lo 97939.495 kg/h
Total 101584.99 kg/h
Salida
G1 3131.034 kg/h
Ln 98453.96 kg/h
Total 101584.991 kg/h

Error 0.000

ESTIMACIÓN DEL DIÁMETRO DE LA TORRE

𝑃𝑎𝑟𝑎 𝑔𝑎𝑠𝑒𝑠 𝑖𝑑𝑒𝑎𝑙𝑒𝑠:
𝑝𝑉=𝑛𝑅𝑇
𝑝𝑉=𝑚/𝑃𝑀 𝑅𝑇
Masa Molar ponderada 31.45 kg/kmol

𝑝.𝑃𝑀= 𝜌 𝑅𝑇
 𝑃𝑎𝑟𝑎 𝑔𝑎𝑠𝑒𝑠 𝑖𝑑𝑒𝑎𝑙𝑒𝑠:
𝑝𝑉=𝑛𝑅𝑇
𝑝𝑉=𝑚/𝑃𝑀 𝑅𝑇
𝑝.𝑃𝑀= 𝜌 𝑅𝑇
Densidad GR 1.25 kg/m3
Densidad LP (AMS) 1235 kg/m3
FLV 0.996
DP Relleno Intalox Cerámica 1-1/2 in 1.09 mm H2O/ft
TABLA 9.28 3.58 mm H2O/m

K4 Diseño 0.13
GRAFICO
K4 Inundación 0.7
% inundación 43%

μ AMS 2.50E-03 Pa.s

Fp (Intalox Cerámica 38 mm) - TABLA 11.2 300 1/m
V*W 0.43 kg/s m2
G' (másico) 0.8684 kg/s
Área de Flujo en zona rellena (AF) 1.9974 m2

Diámetro de Torre Calculado 1.59 m

Diámetro Escogido 1.59471280872134 m
Área de flujo recalculada 2.00 m2
% inundación corregido 43%

Estimación de Altura de Torre (HOG) con Método de Cornell - CALCULO DE DIFUSIVIDADES Y Nº DE SC

Volumen de disfusion aire 20.1

Volumen de difusion SO2 41.1
Dv (difusividad del gas rico) 1.31E-05 m2/s

φ 1
Vm Agua (Teb) 0.0189 m3/kmol
DL (difusividad del liquido) 6.52E-07 m2/s

Viscosidad GR 1.71E-05 Pa s
Viscosidad L 2.50E-03 Pa.s

Sc G = μv/(ρv.Dv) 1.04
Sc L=μL/(ρL.DL) 3.1025
 Estimación de Altura de Torre (HOG) con Método de Cornell - CALCULO DE ALTURA DE LA TORRE POR ELEMENTOS

K3 (Figura 11.43) 1

Y (Figura 11.44) 48

L 27.205 kg/ s
L* = L / Ac 13.62 kg/s m2

F n (Figura 11.45) 1.50E-01

Viscosidad del Agua a 30°C 7.98E-04 Pa s

Densidad del Agua a 30°C 995.71 kg/m3

f1 (μL/μw)^0.16 1.20
f2 (ρL/ρw)^1.25 1.31
f3 (σL/ σw)^0.8 1

Altura por etapa 0.5 m

Altura Inicial (HOG * NOG) 5m

Z iterativo 12.19 m

HL 9.92E-02 m

HG 1.15 m

A 0.65
HOG Calculado 1.22 m

Z calculado 12.19 m
Error Relativo Porcentual 0.00%

RESULTADOS FINALES

Altura Redondeada 12.2 m

Diámetro Escogido 1.59471280872134 m

ILIBRIO

moles H2O x y X Y
5.553 0.0001 0.0016 0.0001 0.0016
5.550 0.0003 0.0042 0.0003 0.0042
5.547 0.0004 0.0076 0.0004 0.0077
5.544 0.0006 0.0112 0.0006 0.0113
5.539 0.0008 0.0186 0.0008 0.0189
5.528 0.0014 0.0342 0.0014 0.0354
5.517 0.0020 0.0513 0.0020 0.0541
5.500 0.0028 0.0776 0.0028 0.0842
5.472 0.0043 0.1211 0.0043 0.1377

Base(masa solucion) 100.00 g solucion

PM SO2 64
PM H2O 18

CARACTERIZACIÓN DE LA CORRIENTE GASEOSA

μ (Pa.s) T (K) Gw (kg/h) G(kmol/h) yi

1.28E-05 303 519.295 8.114 0.07
1.81E-05 303 3126.201 107.800 0.93
1.71E-05 303 3645.496 115.914 1

G'= 107.80004829751

MERO DE ETAPAS

Con la relacion y1/y2 se ingresa al gráfico. Se toman valores

arbitrarios de numero de etapas NOG y se obtiene del grafico
el factor de absorcion A y para cada uno se calcula el L' (caudal
de liquido absorbente AMS necesario)
L' (kmol/h)
3.54E+04
1.18E+04
 8.42E+03
6.67E+03
5.89E+03
5.44E+03
5.20E+03
4.84E+03
4.72E+03
4.59E+03
4.48E+03
4.31E+03
4.21E+03
4.16E+03
4.14E+03
4.07E+03
4.02E+03
3.97E+03
3.95E+03
3.93E+03
L' vs NOG
4.00E+04
3.50E+04
3.00E+04
2.50E+04
2.00E+04
Ls

1.50E+04
1.00E+04
5.00E+03
0.00E+00
0 5 10 15 20 25 30
Etapas

-(L'/G') Límite

Puntos de operación - RECTA DE OPERACIÓN

x y X Y
0 0.0007 0 0.0007 Top
? 0.07 0.0753 Bottom

L' 5441.083 kmol/h

G 115.914 kmol/h
G' 107.800 kmol/h
(L'/G') 50.474
 Recta de operación real
XA YA
0 0.0007
0.001 0.075

𝐺^′ 𝑌1+1+𝐿^′ 𝑋𝑜=𝐺^′ 𝑌𝑛+𝐿^′ 𝑋𝑛

Del balance de materia:

𝐺^′ (𝑌1−𝑌𝑛+1)=𝐿^′ (𝑋𝑜−𝑋𝑛)

Despejando

𝑋𝑛=𝑋𝑜−(𝑌1−𝑌𝑛+1)/(𝐿^′/𝐺′)
Por lo tanto:

MASA

G (kmol/h) Gw (kg/h) PM (kg/kmol)

G1 107.876 3131.034
SO2 0.076 4.833 64
Aire 107.800 3126.201 29
y1 0.001
Y1 0.001
G' 107.800
so porque tiene
en torres empacadas,
ante lavado con una

el gas y el líquido, lo
bsorción del SO₂.
resistencia al paso del

ión química causada

una distribución

cto y resistencia
ales moderados.
ocesos de absorción de L (kmol/h) Lw (kg/h) PM (kg/kmol)
LN 5449.12 98453.96
AMS 5441.08 97939.49 18
SO2 8.04 514.46 64
Aire 0.00 0.00 29
xn 0.0015
Xn 0.0015
L' 5441.08

IMPORTANTE: El error tiene que ser cero o muy cercano

TRO DE LA TORRE
𝑃𝑎𝑟𝑎 𝑔𝑎𝑠𝑒𝑠 𝑖𝑑𝑒𝑎𝑙𝑒𝑠:
𝑝𝑉=𝑛𝑅𝑇
𝑝𝑉=𝑚/𝑃𝑀 𝑅𝑇
𝑝.𝑃𝑀= 𝜌 𝑅𝑇
 𝑃𝑎𝑟𝑎 𝑔𝑎𝑠𝑒𝑠 𝑖𝑑𝑒𝑎𝑙𝑒𝑠:
𝑝𝑉=𝑛𝑅𝑇
𝑝𝑉=𝑚/𝑃𝑀 𝑅𝑇
𝑝.𝑃𝑀= 𝜌 𝑅𝑇

ell - CALCULO DE DIFUSIVIDADES Y Nº DE SCHMIDT

O DE ALTURA DE LA TORRE POR ELEMENTOS DE TRANSMISION

PROCEDIMIENTO ITERATIVO
1) Se supone una altura inicial por etapa y se multiplica por el núme
etapas obteniendo una altura inicial.
2) Luego se calcula con ese valor de Z, el HL y HG con las expresio
correspondientes.
3) A partir de dicho valor, se obtiene el HOG (altura de la etapa)
multiplica por el numero de etapas obteniendo un nuevo Z (Z calculado)
4) Se calcula el error porcentual entre el Z supue
Z calculado. Si es bajo, el Z es correcto, si no, se comienza una nueva iterac
tomando como Z inicial el Z calculado en la etapa anterior
 2) Luego se calcula con ese valor de Z, el HL y HG con las expresio
correspondientes.
3) A partir de dicho valor, se obtiene el HOG (altura de la etapa)
multiplica por el numero de etapas obteniendo un nuevo Z (Z calculado)
4) Se calcula el error porcentual entre el Z supue
Z calculado. Si es bajo, el Z es correcto, si no, se comienza una nueva iterac
tomando como Z inicial el Z calculado en la etapa anterior

FINALES
 0.1600

0.1400
f(x) = 32.8084355403039 x − 0.00711037565716675
0.1200 R² = 0.99542202697113

0.1000

0.0800
Y

0.0600

0.0400

0.0200

0.0000
0.0000 0.0005 0.0010 0.0015 0.0020 0.0025 0.0030 0.0035 0.0040 0.0045
X

0
20 25 30

0.1600
0.1400
0.1200
0.1000
0.0800
Y

0.0600
0.0400
0.0200
0.0000
-0.0010 0.0000 0.0010 0.0020 0.0030 0.0040 0.0050 0.0060
-0.0200
X

Equilibrio Linea de Operación

LINEA AUX 2 L'/G' limite
 (L'/G') limite
XA YA Base 25
0 0.0007 Incremento 0.0001
0.0025 0.075
(L'/G') limite 29.83 L' min 3215.387695
% aumento 69.22
Para encontrar la relacion (L'/G') limite, se iteran (con botonera)
Linea auxiliar 2 diferentes valores de Xn, hasta encontrar la recta que une el punto de la
XA YA cabeza de la torre (Xo, Y1) con el punto de la curva de equilibrio
0 0.0753 correspondiente a la entrada de la torre Yn+1. La pendiente de la recta
de operación indica el L'/G' limite
0.8 0.0753

Completar tercero

L (kmol/h) Lw (kg/h) PM (kg/kmol)

Lo 5441.083
AMS 5441.083 97939.495 18 Completar segundo
SO2 0.000 0.000 64
Aire 0.000 0.000 29
xo 0.000
Xo 0.000
L' 5441.083

G (kmol/h) Gw (kg/h) PM (kg/kmol)

GN+1 115.91 3645.50
SO2 8.11 519.29 64
Aire 107.80 3126.20 29
yn+1 0.07
Yn+1 0.08
G' 107.80

Completar primero
L
 al por etapa y se multiplica por el número de

e valor de Z, el HL y HG con las expresiones

, se obtiene el HOG (altura de la etapa) y se

eniendo un nuevo Z (Z calculado)
ula el error porcentual entre el Z supuesto y el
cto, si no, se comienza una nueva iteración
en la etapa anterior
e valor de Z, el HL y HG con las expresiones

, se obtiene el HOG (altura de la etapa) y se

eniendo un nuevo Z (Z calculado)
ula el error porcentual entre el Z supuesto y el
cto, si no, se comienza una nueva iteración
en la etapa anterior
teran (con botonera)
ta que une el punto de la
la curva de equilibrio
La pendiente de la recta
limite
 DATOS EQUILIBRIO agua
x y
0 0
0.025 0.99191
0.05 0.998896
0.075 0.999694
0.1 0.99988
0.15 0.999969

0.2 0.999988
0.3 0.999996
0.4 0.999998
0.5 0.999999
 DATOS EQUILIBRIO
x y X Y

0 0 0 0
0.025 0.999988 0.025641026 83332.33333
0.05 0.999994 0.052631579 166665.6667
0.1 0.999997 0.111111111 333332.3333
0.15 0.999998 0.176470588 499999
0.2 0.999998 0.25 499999
0.3 0.999998 0.428571429 499999
0.4 0.999999 0.666666667 999999
0.5 0.999999 1 999999
0.6 0.999999 1.5 999999

Chart Title
600000

500000
f(x) = 2825112.44442848 x + 9949.95729136627
R² = 0.997973785385208

400000

300000

200000

100000

0
0 0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.16 0.18 0.2
100000

0
0 0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.16 0.18 0.2