Ejercicio 7.

14

Se van a enfriar, a 27 ºC, 5 m3/s (10600 ft3/min) de un gas pobre, 65 % N2,

35% CO, inicialmente seco a 1 atm a 315 ºC, mediante contacto a
contracorriente con agua que entra a 18 ºC. Se va a utilizar una torre
empacada con anillos de Rasching de cerámica de 50 mm, con L’2/G’2 = 2.0.
Especificar el diámetro de una torre adecuada para este propósito y la altura
empacada.

G2

TL2 Y2

L2 Tg2

L2

TL1
G2

L1 Y2

L2 Tg2

Definimos para la torre

A=como agua

B= como mezcla de monóxido de carbono y nitrógeno

Para la corriente de gas entrante

Y’1 = 0

M B = PMN2*0.65+PMCO*0.35 = 28.012 kg/kmol

MG1 = M B = 28.012 kg/kmol

PT = 1 atm = 101325 Pa

TG1 = 315 ºC = 588.15 K

M G1 * PT
DG1= = 0.5805 kg/m3
R * TG1

Q = 5 m3/s

G1 = Q*G1 = 5*0.5805 = 2.9025 kg/s

Gs = G1 = 2.9025 kg/s

Para la corriente gas saliente

TG2 = 27 ºC =300.15 K

Y’s = 0.0234526 kg A /kg B

Se supone que el gas a la salida esta casi saturado entonces

Y’2= 0.021213 kg A /kg B (valor a verificar)

G2 = G1*(1+ Y’2) = 2.9641 kg/s

MG2 = 27.6928 kg/kmol

M G 2 * PT
G2 = = 1.1245 kg/m3
R * TG 2
A la hora de obtener el valor G’2 debemos buscar en la grafica de Eckert, en la
abcisa y con condiciones en el tope de la torre ( a T L2 = 18 ºC , DL2 =998.595
kg/m3 tomada del Perry ) :

1/ 2
L´   G 
   0.0672
G´   L   G 

De la figura 6.34 del Treybal para una caída de presión de 400 (N/m2)/m. Se
lee en la ordenada (de la tabla 6.3 para anillos Rasching Cf = 65 ; µL = 1.056E-
3 Pa*s tomada de Hysys ) :

G´2 Cf L J
0. 1
= 0.073 J=1, gc = 1 para unidades S.I.
G (  L  G ) gc

De esta ecuación obtenemos el valor de 1.5811 kg/m2s

Luego encontramos la sección transversal de la torre será G2/G’2 = 1.8747 m2

El diámetro correspondiente es 1.545 m, que se puede aproximar a 1.5 m para

el cual el área de sección transversal es = 1.7671 m2 , realizando la corrección
para Gs se tiene:

G’s = 1.64235 kg /m2s

G’2 = 1.6772 kg /m2s

L’2 = 3.3544 kg /m2s

Haciendo un balance global de agua: L’2 = L’1 - G’2* Y’2

L’1 = 3.3195 kg /m2s

Las capacidades caloríficas se obtuvieron como un valor promedio en el

intervalo de temperaturas de trabajo:
CpN2 =1058.35 J/(Kg*K)

Cp CO=1070.27 J/(Kg*K)

CpA = 1930.27 J/(Kg*K)

CAL = 4187 J/kg K.

De donde:

C’H=1062.52+1930.27*Y’

Utilizando como temperatura base, To = 0 °C, para la cual λo = 2501400 J/kg.

H’=(1062.52+1930.27*Y’)*T+2501.4*Y’

A partir de un balance global de entalpía:

L’2*H L2 +G’s*H1’ = L’1*H L1 + G’s*H’2

se obtiene : TL1 = 47.761 °C,

temperatura para la cual se tienen las siguientes propiedades del liquido:

μL1 = 5.6705E-4 Pa*s

L1= 989.022 kg/m3

Conductividad térmica (kth) = 0.64077 W/mºC

PrL = 3.7026

Para el gas entrante:

μG1 = 2.9945E-5 kg/ms

Difusividad del vapor de agua (DAB) = 8.3684E-5 m2/s

ScG = 0.61648
Se procede a realizar el cálculo de la superficie interfacial, de la tabla 6.5 para
anillos Rasching:

ds = 0.0725 m

β = 1.508ds0.376 = 0.562

(2.09 *10 6 )(737.5 * L´) 

φLtw =  0.03196 m3/ m3
ds 2

2.47 *10 8
φLsw =  0.005911 m3/ m3
ds1.21

φLow = φLtw - φLsw = 0.02605 m3/ m3

0.1737 0.262LOGL'
975.7 L'0.57  L
  
0.13

H  0.84    0.93062
 L (2.024 L  1)  0.073 
'0.43

φLo = φLow * H =0.024245

0.0486 L  0.55
0.02

 LS   0.00549258 m 3 / m 3
dS
1.21
L 0.37

LT = LS + LO = 0.0297085 m3/m3.

De la tabla 6.4 se tiene m = 34.03; n = 0; p = 0.362

n
 808G´  P
aAw = m 
  0.5  L´ = 52,5402751 m / m
2 3

 ´G 

avw = 0.85 avw*( φLtw / φLow ) = 54,7923116 m2/ m3

De la tabla 6.3 se lee:  = 0.73

 Lo =  - φLt

Lo = 0,71029148

Se calcula FG a través de la ecuación y luego se multiplica por el área avw :

0.36
 
2/3  
FG S CG  d s G' 
 1.195
G   G (1   Lo ) 
 
 

FG = 0,00313353 kmol/m2s

FGa = 0,17169312 kmol/m3s

Se calcula hG a través de las ecuaciones:

2/3
FG SCG J H C pG '
JH  hG  2/3
G PrG

hG = 82,5742275 W/m2K

hGa = 4524,4328 W/m3K

hL se calcula a partir de:

0.45
k  d L' 
hL  25.1 th  s 2 
0 .5
PrL
ds  L 

hL = 6498,35056 W/m2K

hLa = 356059,649 W/m2mK

el cálculo de la temperatura en la interfase en el tope de la torre requiere un
proceso iterativo que permite encontrar simultáneamente Ti , dY’/dZ, y dTg/dZ
.

El orden de cálculo es el siguiente :

1. Se supone un valor de Ti . Al inicio Ti=TL
2. Se calcula la presión en la interfase ( Pi) a través de la siguiente
ecuación tomada del Perry.
Pi = Exp(81.768 - 6876 / (Ti + 273.15) - 8.7078 * Log(Ti + 273.15)
+ 0.0000071926 * (Ti + 273.15) ^ 2)

3. Se resuelven consecutivamente las expresiones dY’/dZ, h’Ga, dTG /dZ

dadas por las ecuaciones :

dY ' M A FG a 1  PAi / PT
 ln
dZ G' s 1  PAG / PT

G ' sC A ( dY ' / dZ )
h'G aH  
1  eG ' sC A ( dY ' / dZ ) / hG a

dTG h' a (T  Ti)

 G H G
dZ G' sCs

4. Se calcula un valor de Ti mediante la siguiente ecuación:


Gs' * C s *  G
dT   C * T  C * T  C  C * T   * dY ' 
  dZ  A G A, L L A, L A o o dZ 
Ti cal  TL  '
Gs' * C A, L * (dY )  hL a
dZ

5. Hacemos Ti sup = Ti cal.

El proceso iterativo termina cuando Ti supuesta - Ti cal <1E-4.

Al final de este procedimiento se obtuvieron los siguientes valores:

 Valores en el tope

Ti (ºC) 47,9579421

dY’/dz 2,1927E-01

dTg/dz -640,546802

Posteriormente se determina un T con el cual va a variar la temperatura del

gas a lo largo de la torre, el cual toma un valor cada vez menor a medida que
nos acercamos al tope.

Se calcula un Z con la siguiente ecuación:

T
Z 
dt G dZ 

Con el T se calcula la nueva TG

TG K+1= TG K + T

En el fondo se tomó T= -30 ºC, de donde TG =285ºC y Z=0,04683499 m

Después se halla la nueva Y’

dY '
Y ' K 1  Y ' K  * Z
dZ

Y’ = 0,010269565 KgA/KgB
Finalmente se calcula, mediante un balance de masa, el L’ K+1 y mediante un
balance de entalpía se obtiene T L K+1 .

L’ K+1 = 3,336410286 Kg/m2s

T L K+1 = 47,85963829 ºC

Se repite el proceso hasta que la Tg alcance el valor del tope, momento en el

cual se debe verificar el Y’2 supuesto. Por último se realiza la sumatoria de los
Z para obtener la altura de la zona empacada en la torre.

Las propiedades de los fluidos se han calculado en cada sección de acuerdo a

las siguientes expresiones, obtenidas correlacionando datos tomados de
Haysys (H) o del Perry (P):

µ L (Pa*s) = -3,65744E-09*T L
3
+ 6,44446E-07*T L
2
- 4,61041E-05*T L +
1,69744E-03 (H)

R2 = 9,99952E-01

 L (Kg/m3) = 2,22315E-05*T L
3
- 6,46541E-03*T L
2
+ 2,69773E-02*T L +
1,00006E+03 (P)

R2 = 9,99998E-01

 L (N/m) = -1,74355E-04* T L + 7,64495E-02 (H)

R2 = 9,99967E-01
K th (W/mK) = -7,81352E-06*T L
2
+ 1,87756E-03*T L + 5,68926E-01 (H)

R2 = 9,99998E-01

Se consideró la variación de la viscosidad del gas con la temperatura y la

composición por lo que se realiza una interpolación a partir de los siguientes
valores obtenidos de Hysys
La difusividad en el gas se obtuvo considerando la difusividad del agua en el N2
y en el CO por separado, de acuerdo a los métodos presentados en el capítulo
2 del Treybal:

DA-N2 = (0,000000002900874*(Tg+273,15)^1,5))/ (-
0,0006*(Tg+273,15)+0,8482)

DA-CO = (0,000000003000728*(Tg+273,15)^1,5))/ (-
0,0006*(Tg+273,15)+0,8627)

D A-G = 1/(0.65/ DA-N2 +0.35/ DA-CO )

Los resultados de todos los cálculos se muestran a continuación habiéndose

trabajado todo tal como se presentaron los cálculos para el tope.

Estos resultados corresponden a la última iteración del cálculo de la torre en la

cual efectivamente se comprueba que Y’2= 0.021213 kg A /kg B
DIAGRAMA DE H´ Vs Tg
 Temperaturas y humedades para ejemplo

0.050
0.045
Y'= Kg de agua/ Kg de gas seco

0.040
0.035
0.030
0.025
0.020
0.015
0.010
0.005
0.000
0 50 100 150 200 250 300 350
Temperatura ºC
 TRABAJO DE OPERACIONES UNITARIAS

(Ejercicios 7.11 y 7.14)

JAVIER DARIO MALDONADO RINCON

Presentado a:

DR. CRISOSTOMO BARAJAS FERREIRA

UNIVERSIDAD INDUSTRIAL DE SANTANDER

INGENIERIA QUIMICA

BUCARAMANGA

2008