DESARROLLO DE LOS EJEMPLOS DEL LIBRO ROBER E.

TRAYBAL
CAPITULO 12

EJEMPLO 12.4

Un sólido cristalino insoluble húmedo con agua se coloca en una charola de 0.7 m
por 0.7 m, 25mm de profundidad, hecho de fierro galvanizado de 0.8mm de
espesor. La charola se coloca en una corriente de aire a 65°C con humedad de
0.01 kg de agua/ kg de aire seco, que fluye en forma paralela a la superficie
superior e inferior a una velocidad de 3 m/s.

La superficie superior del sólido esta colocada directamente hacia las tuberías
calentadas con vapor cuya temperatura superficial es de 120 °C a una distancia de
la superficie del sólido de 100mm

a) Calcule la rapidez del secado a rapidez constante.

Con las condiciones de: Y’ = 0.01 kg de agua/ kg de aire seco y tG = 338 K
Se calcula:
V’H = 0.97272 m3/Kg aire seco

m 1.01
G    1.03833 kg/ m3
VH 0.97272

GS  v *  G  3*(1.03833)= 3.11497 kg/ m2*s

4(sec c.transv. flujo) 4 * (0.7)(0.1)

d  = 0.175 m
perimetro (0.7  0.1) * 2

Para calcular hc es posible usar la siguiente expresión para una temperatura del
aire de 45 a 150 ˚C y una velocidad de masa G de 2450-29300 Kg/h*m2 o una
velocidad de 0,61 a 7,6 m/s.[1]

hc  0,02004 * G 0,8

hc  0,0204 * (11232) 0,8

hc  35,5 W/m2*K

Se toma la emisividad del sólido como  =0.94, como Tr=393 K

  (5.73 *10 8 )(TR 4  TS 4 ) 0.94(5.73 *10 8 )(393 4  (TS  311) 4 )
hR   = 9.52407
(TR  TS ) (393  (TS  311))

(5.39 *10 8 )( 2.38 *1010  (TS  273.15) 4 ) 1282.82  5.39 *10 8 (TS  273.15) 4 )
hR  hR 
(393  (TS  273.15)) (393  (TS  273.15))

1
UK 
(1 / 35.5)(0.49 / 0.56)  (0.0008 / 45)(0.49 / 0.56)  (0.025 / 3.5)(0.49 / 0.49)

U K  31.44W / m 2 k

Con lo cual se llega a:

( y s  y1 )(s *10 3 ) 31.44 1282.82  5.39 *10 8 (TS  273.15) 4 (TR  TS )

UK   (1  )(TG1  TS ) 
1023.8 35.5 (TR  TS ) * 35.5

Utilizando el programa ejemplo 12.4 que consiste en un proceso iterativo se halla

λs

TS YS y Nc.

λs = 2413.4 KJ/Kg aire seco

TS =37.2

YS =0.0415 Kg agua/Kg aire seco

Nc.=5.0056*10-4 Kg de agua evaporada

b) Cuando no hay radiación ni conducción de calor a través del solido, la

superficie que se esta secando toma la temperatura de bulbo húmedo del
aire. Para el sistema aire – agua a esta humedad, las líneas de saturación
adiabática de la carta psicométrica sirven como líneas de bulbo húmedo. La
línea que interseca la curva de humedad de saturación en la condición de
bulbo húmedo, TS = 28.5 ˚C YS = 0.025. A esta temperatura, λs = 2434
KJ/kg.
 hc (t g  t s ) 35.5(338.15  301.65)
NC    5.3235 *10 4 Kg / m 2 s
( s ) 2434 *10 3

Y la rapidez de evaporación es ( 5.3235 *10 4 )(0.49) = 2.6085*10-4 Kg/s

Ejemplo 12.5

Una torta de un precipitado cristalino se va a secar pasando aire a través de la

torta. Las partículas de la torta no son porosas; su diámetro promedio es de 0.20
mm; puesto que son insolubles en agua, tienen un contenido despreciable de
humedad en el equilibrio. La torta tiene 18 mm de espesor y la densidad aparente
es 1350 kg sólido eco/m3 (85 lb,/ft’). Se va a secar desde 2.5 hasta 0.1% de
humedad. El aire va a entrar a la torta a 0.24 kg aire seco/s(m2 sección transversal
lecho)( = 177 lbm/h ft2), a una temperatura de bulbo seco de 32 °C y 50% de
humedad. Calcular el tiempo para el secado.

Calculo Humedad absoluta

𝑀𝐴 𝜑𝑃 𝑠
𝑌′ =
̅𝐵 𝑃 − 𝜑𝑃𝑠
𝑀

𝑘𝑔 𝑎𝑔𝑢𝑎
𝑌 ′ = 0,0149 ⁄𝑘𝑔 𝑎𝑖𝑟𝑒 𝑠𝑒𝑐𝑜

Calculo de la T saturación adiabática y Y’

SA
TSA  TG1   YSA  Y1 
CH 1
𝐾𝑔𝑎𝑔𝑢𝑎
Y’sa= 0.018430 𝐾𝑔𝑎𝑖𝑟𝑒𝑠𝑒𝑐𝑜

Tsa= 23.629℃

En base seca calculamos el 𝑋1 𝑦 𝑋2

𝑥
𝑋=
1−𝑥
𝑋1 = 0,0256𝐾𝑔𝑎𝑔𝑢𝑎/𝑘𝑔𝑠𝑜𝑙𝑖𝑑𝑜𝑠𝑒𝑐𝑜

𝑋2 = 0,001𝐾𝑔𝑎𝑔𝑢𝑎/𝑘𝑔𝑠𝑜𝑙𝑖𝑑𝑜𝑠𝑒𝑐

Calculo de 𝐺𝑎𝑝𝑟𝑜𝑥𝑖𝑚𝑎𝑑𝑜 𝑝𝑟

0,24(0,0149 + 0,018430) 𝑘𝑔
𝐺𝑎𝑝𝑟𝑜𝑥𝑖𝑚𝑎𝑑𝑜 𝑝𝑟 = 0,24 + = 0,24799 ⁄𝑚2 ∙ 𝑠
2

Calculo de 𝑁𝑚𝑎𝑥

𝑁𝑚𝑎𝑥 = 𝐺𝑠 (𝑌𝑠𝑎 − 𝑌1 ) = 8.472 ∗ 10−4𝐾𝑔𝑒𝑣𝑎𝑝 / 𝑚2 𝑠

Dado el espesor del solido 𝑧𝑠 =0.018m y su densidad 𝜌𝑠=1350𝑘𝑔𝑠𝑜𝑙𝑖𝑑𝑜𝑠𝑒𝑐𝑜/𝑚3

Con el promedio de las temperaturas calculamos la viscosidad del aire a T=28℃

𝜇 = 1.8 ∗ 10−5 𝑘𝑔/𝑚𝑠

0,215
0,273 𝑑𝑝 𝐺
𝑁𝑡𝐺 = 0,35 ( ) (𝑋𝜌𝑠 𝑧𝑠 )0,64
𝑑𝑝 𝜇

𝑑𝑝 𝐺 (2.0 ∗ 10−4 ∗ 0,24799)

= = 2,7554
𝜇 1.8 ∗ 10−5
0,215
0,273 𝑑𝑝 𝐺
𝑁𝑡𝐺 = 0,35 ( ) (𝑋𝜌𝑠 𝑧𝑠 )0,64 = 51.55𝑋 0.64
𝑑𝑝 𝜇

𝑁
= 1 − 𝑒 −𝑁𝑡𝐺
𝑁𝑚𝑎𝑥

𝑁 = 1 − 𝑒 −𝑁𝑡𝐺 ∗ 𝑁𝑚𝑎𝑥 = 8,472 ∗ 10−4 ∗ (1 − 𝑒𝑥𝑝(51,55𝑋 0,64 )

𝑋1
𝑆𝑠 𝑑𝑋
𝑡 = ( )∫
𝐴 𝑋2 𝑁
Esta integral se puede resolver numéricamente por el método de trapecio cuya
expresión es:

 Yo  Yn 
b

 f ( x)dx  h
a
2
 Y1  Y2  Y3  Y4  ...  Yn 


Xn  Xo
En donde h  y n=particiones
n
𝑥𝑛 = 0,0256𝑦 𝑥0 = 0,001 𝑛 = 10 ℎ = 0,00246

X N 1/N
0,0256 0,00084112 1188,8866
0,02314 0,00083892 1192,00638
0,02068 0,00083579 1196,4791
0,01822 0,00083124 1203,02481
0,01576 0,0008245 1212,85421
0,0133 0,00081425 1228,11709
0,01084 0,00079813 1252,92752
0,00838 0,00077157 1296,06374
0,00592 0,00072475 1379,78037
0,00346 0,00063227 1581,59911
0,001 0,00039135 2555,27432

3000

2500
m^2.S/Kg de agua evap

2000

1500
1/N
1000

500

0
0 0.01 0.02 0.03
x kg agua/kg de solido seco
 𝑠
Puesto que ( 𝐴𝑠 ) = 𝜌𝑆 𝑍𝑆 = 1350 ∗ 0,018 = 24,3𝐾𝑔𝑠𝑜𝑙𝑖𝑑𝑜𝑠𝑒𝑐𝑜/𝑚2

El tiempo de secado es:

𝑡 =33,0007346*24,3 = 801.92 𝑆 = 13.36𝑚𝑖𝑛𝑢𝑡𝑜𝑠

EJEMPLO 12.6

Se van a secar unas lentejas de catalizador porosas y húmedas, cuya forma es de

pequeños cilindros; tienen 13.5 mm de diámetro, 13.0 mm longitud. Se secaran
mediante un secador de circulación transversal. Las lentejas se van a colocar en
lechos de 50 mm de profundidad, sobre mallas y se van a secar con aire que fluye
a 1.1 kg aire seco/s (m2 sección transversal lecho); el aire entra a una temperatura
de bulbo seco de 82 ºC, humedad 0.01 kg agua/kg aire seco.
La densidad aparente del lecho es 600 kg sólidos secos/m3; la superficie de la
partícula es 280 m2/m3 de lecho. Calcular la rapidez de secado, y la humedad y
temperatura del aire que sale del lecho durante el periodo de rapidez constante.

Dado: Y’ = 0.01 kg de agua/ kg de aire secoy tG = 82 ºC

Calculando
SA
TSA  TG1   YSA  Y1 
CH 1

𝑌𝑆𝑎 = 0.0310 kg de agua/ kg de aire seco

𝑇𝑆𝑎 = 32.2102 ºC

𝐺𝑆 = 1.1𝐾𝑔 𝑑𝑒 𝑎𝑖𝑟𝑒 𝑠𝑒𝑐𝑜/𝑚2 𝑠

El 𝐺 promedio aproximado
1.1(0.0310 + 0,01) 𝑘𝑔
𝐺 = 1.1 + = 1.1224 ⁄𝑚2 ∙ 𝑠
2

Viscosidad promedio del aire aproximada 𝜇 = 1.9*10-5 Kg/m.s

Superficie de cada partícula:

La superficie de cada partícula es la suma de las áreas del cilindro y de la círculo

𝑠𝑝 = 2𝜋R2 + 𝜋*R*h= 2*𝜋*(0.01352/4) + 𝜋*(0.0135*0.0130) = 8.376*10-4 m2

El diámetro de una esfera de la misma area se halla despejando esta ecuación
𝑠𝑝 = 𝜋𝑑𝑝2 donde 𝑑𝑝 = 0.01633 mm

a= 280 m2/m3
zs= 0.050 m

𝑑𝑝 𝐺
𝑅𝑒 =
𝜇

𝑑𝑝 𝐺 (0.01633 ∗ 1.1225)
𝑅𝑒 = = = 964.759
𝜇 1.9 ∗ 10 − 5

𝐴𝑑𝐴 0.01633 ∗ 280

𝜀 =1− =1− = 0.2379
6 6
2.06
𝐽𝐷 = ( ) ∗ 𝑅𝑒 −0.575 = 0.1665
𝜀
𝑗𝑑 ∗ 𝐺𝑠 0.2574𝑘𝑔𝑎𝑔𝑢𝑎
𝑘𝑦 = 2 =
𝑚2 𝑠∆𝑦
𝑠𝑐 3
𝐾𝑌 ∗ 𝐴 ∗ 𝑍𝑆
𝑁𝑡𝐺 = = 3.276
𝐺𝑆

𝑁𝑚𝑎𝑥 = 𝐺𝑆 (𝑌𝑠𝑎 − 𝑌1 ) =1.1*(0.0310-0.01)= 0.0231 Kg/m2*s

𝑁 = 1 − 𝑒 −𝑁𝑡𝐺 ∗ 𝑁𝑚𝑎𝑥

𝑁 = 1 − 𝑒 −𝑁𝑡𝐺 ∗0.0231 = 0.0223 Kg/m2*s en el periodo de rapidez constante.

N (Y2 − Y1)
= = 1 − 𝑒 −𝑁𝑡𝐺
Nmax (Ysa − Y1)

Despejando obtenemos el valor de Y2 = 0.03027 kg de agua/ kg de aire seco para

el aire saliente.

Se calcula la temperatura del gas a la salida, con la ecuación teniendo en cuenta

que la la curva de saturación adiabatica:
 (Y2  Y1 )(Tsa  TG1 )
TG 2   TG1 = 33.94 ºC
(Ysa  Y1 )

EJEMPLO 12.9

Después de ser centrifugadas, se van a secar al aire madejas de rayón húmedo,

de 46 a 8.5% de agua, en un secador de túnel, continuo, a contracorriente. Las
madejas se colgaran de perchas que van a viajar a través del secador. El aire
entrara a 82 ºC, humedad 0.03 kg agua/kg aire seco; se va a descargar con una
humedad de 0.08. La temperatura del aire se mantendrá constante a 82 ºC
mediante espirales de calentamiento dentro del secador. El flujo del aire va a ser
de 1.36 kg/m2*s.
El contenido de humedad crítica de las madejas de rayón es del 50%; su
porcentaje de humedad en el equilibrio a 82 ºC puede tomarse como un cuarto del
porcentaje de humedad relativa del aire. La rapidez de secado es entonces:

𝑑𝑥
− = 0.0137. 𝐺 1.47 (𝑋 − 𝑋 ∗ )(𝑌𝑊 − 𝑌)
𝑑𝑡

en donde Yw, es la humedad de saturación del aire en la temperatura de bulbo

húmedo correspondiente a Y. (Nota: aquí G esta expresada como kg/m 2*s y en
segundos,
Calcúlese el tiempo que deberá permanecer el rayón en el secador.

La humedad en base seca del rayón a la entrada y salida del secador es:

0.46 0.085
X1   0.852 X2   0.093
1  0.46 1  0.085

Sabiendo que Y1=0.08 kg agua/kg aire seco y Y2=0.03 kg agua/kg aire seco

Realizando un balance de agua: Ss( X 1  X 2 )  Gs(Y1  Y2 )

Ss 0.08  0.03
  0.0660 Kg sólido seco/Kg aire
Gs 0.852  0.093

Debido a que el contenido inicial de humedad del rayón es menor a la crítica, el

secado se realiza por completo, dentro de la zona III, es decir, cuando ocurre el
secado de la superficie no saturada y cuando la rapidez de secado depende
estrictamente de las condiciones predominantes en cualquier momento,
independientemente de cómo haya sido hasta entonces el secado.
 dX KyA( X  X  )(Ys  Y ) kyA
  donde  0.0137G 1.47
dt Ss( X c  X  ) 
Ss ( X c  X )

Al sustituir 𝐺= 1.36 kg/m2*s, se obtiene:

𝜃 𝑋1
1 𝑑𝑥
𝜃𝐼𝐼𝐼 = ∫ 𝑑𝜃 = 1.47
∫ ∗
0 0.0137𝐺 𝑋2 (𝑋 − 𝑋 )(𝑌𝑊 − 𝑌)

𝜃 0.852
𝑑𝑥
𝜃𝐼𝐼𝐼 = ∫ 𝑑𝜃 = 46.4 ∫
0 0.093 (𝑋 − 𝑋 ∗ )(𝑌𝑊 − 𝑌)

El calculo de la integral se desarrolla realizando el siguiente procedimiento

1. Suponemos valores de X dentro del intervalo de las humedades dadas

𝑆
2. Para cada valor de X calculamos y’ , 𝑌 ′ = 𝑌2′ + (𝑋 − 𝑋2 ) ∗ 𝐺𝑆 remplazando
𝑆
los valores tenemos
𝑌 ′ = 0.03 + (𝑋 − 0.093𝑋2 ) ∗ 0.066

3. Con T=82℃ y para cada valor de Y’ leemos en la carta psicrometrica, la

humedad correspondiente a saturación Yw
4. Calculamos la humedad relativa
5. Calculamos la humedad en el equilibrio

X Y' Y'w %Hr X* 1/(X-X*)(Y'w-Y')

0,852 0,08 0,0945 22,53 0,05968 87,04
0,7761 0,075 0,0895 21,28 0,05619 95,79
0,7002 0,0701 0,086 20,03 0,0527 97,13
0,6243 0,0651 0,081 18,74 0,04915 109,35
0,5484 0,06 0,076 17,424 0,04554 125,07
0,4725 0,055 0,072 16,0659 0,04185 136,59
0,3966 0,05 0,067 14,708 0,03817 164,11
0,3207 0,045 0,062 13,342 0,0345 205,53
0,2448 0,04 0,058 11,95 0,0308 259,6
0,1689 0,035 0,054 10,54 0,0271 371,17
0,093 0,03 0,048 9,1 0,0233 797,07
  Yo  Yn 
b

 f ( x)dx  h
a
2
 Y1  Y2  Y3  Y4  ...  Yn 


Xn  Xo
En donde h  y n=particiones
n

Donde el tiempo de secado es

𝜃𝐼𝐼𝐼 =46.4*152,285=70.66 s = 1.96 h

900
800
700
dX/(X-X*)(Yw-Y)

600
500
400
300 Series1
200
100
0
0 0.2 0.4 0.6 0.8 1
X
DESARROLLO DE LOS EJEMPLOS DEL LIBRO ROBER E. TRAYBAL
CAPITULO 12

LAURA SMITH MORENO ARCINIEGAS

INGRID JOHANA SUAREZ P

LILIANA SOFIA GOMEZ RODRIGUEZ

FERNADO CORZO ALVAREZ

JUAN CARLOS PATIÑO

ANDRES ESTEBAN AMAYA

RICARDO ALONSO BARRERA

JUAN CARLOS DIAZ

SERGIO ALBERTO SOLANO

HARVEY STEVE ROJAS

JESUS DAVID BADILLO

UNIVERSIDAD INDUSTTRIAL DE SANTANDER

OPERACIONES UNITARIAS I

INGENIERIA QUIMICA

BUCARAMANGA

2008
 PROGRAMAS

 Programa para calcular Ys y Ts:

tg1=input('temperatura al a entrada del secador ºC= ');

Y1=input('HUEMDAD a la entrada del secador Kg agua/Kg gas seco = ');

tr=input('temperatura de radiación en el secador ºC= ');

A1=23.7093;

B1=4111;

C1=237.7;

A2= 23.1863;

B2=3809.4;

C2=226.7;

a=0;

for i=1:0.1:100

t=i-1;

if (t>=0)

if (t<57)

Ps=exp(A1-B1/(C1+t));

Else

Ps=exp(A2-B2/(C2+t));

end

end

Ys=0.622*Ps/(101325-Ps);
for k=1:1:101

T(k)=k-1;

end

CC=[2501.4,2494.0,2494.3,2496.7,2491.9,2489.6,2487.2,2484.8,2482.5,2480.1,24
77.7,2475.4,2473,2470.7,2468.3,2465.9,2463.5,2461.2,2458.8,2456.5,2454.1,245
1.8,2449.4,2447.0,2444.7,2442.3,2439.9,2437.6,2435.2,2432.8,2430.5,2428.1,242
5.7,2423.4,2421,2418.6,2416.2,2413.9,2411.5,2409.1,2406.7,2404.3,2401.9,2399.
5,2397.2,2394.8,2392.4,2390,2387.6,2385.2,2382.7,2380.3,2377.9,2375.5,2373.1,
2370.7,2368.2,2365.8,2363.4,2360.9,2358.5,2356.0,2353.6,2351.1,2348.7,2346.2,
2343.7,2341.3,2338.8,2336.3,2333.8,2331.4,2338.9,2326.4,2323.9,2321.4,2318.9,
2316.3,2312,2311.3,2308.8,2306.2,2303.7,2301.1,2298.6,2296,2293.5,2290.9,228
8.3,2285.8,2283.2,2280.6,2278,2275.4,2272.8,2270.2,2267.6,2264.9,2262.3,2259.
7,2257];

for j=1:101

CC(j);

end

t1=fix(t);

t2=t1+1;

CC1=CC(t1+1);

CC2=CC(t2+1);

lambda=CC1-((t1-t)/(t1-t2))*(CC1-CC2);

lambda;
F=abs((1+(31.44/35.5))*(tg1-t)+((1282.82-5.39*10^(-8)*t)*(tr-t)/((393-t)*35.5))-((Ys-
Y1)*(lambda*1000)/(1023.8)));

if F<=0.5

a=1;

break

else

end

k=0;

j=0;

end

if a==1

Nc=(35.5+31.44*(tg1-t)+((1282.82-5.39*10^(-8)*t)/(393-t))*(tr-t))/(lambda*1000)

Ys

else

disp ('EL SECADOR NO ES OPERABLE A LAS CONDICIONES DADAS')

end
  Programa para calcular Ts

clc

clear

t=input('temperatura de saturación en ºC: ');

tg1=input ('temperatura a la entrada del gas en °C: ')

while t<0

disp('la temperatura debe estar entre 0 y 100')

t=input('temperatura: ');

end

for i=1:101

T(i)=i-1;

end

CC=[2501.4,2494.0,2494.3,2496.7,2491.9,2489.6,2487.2,2484.8,2482.5,2480.1,24
77.7,2475.4,2473,2470.7,2468.3,2465.9,2463.5,2461.2,2458.8,2456.5,2454.1,245
1.8,2449.4,2447.0,2444.7,2442.3,2439.9,2437.6,2435.2,2432.8,2430.5,2428.1,242
5.7,2423.4,2421,2418.6,2416.2,2413.9,2411.5,2409.1,2406.7,2404.3,2401.9,2399.
5,2397.2,2394.8,2392.4,2390,2387.6,2385.2,2382.7,2380.3,2377.9,2375.5,2373.1,
2370.7,2368.2,2365.8,2363.4,2360.9,2358.5,2356.0,2353.6,2351.1,2348.7,2346.2,
2343.7,2341.3,2338.8,2336.3,2333.8,2331.4,2338.9,2326.4,2323.9,2321.4,2318.9,
2316.3,2312,2311.3,2308.8,2306.2,2303.7,2301.1,2298.6,2296,2293.5,2290.9,228
8.3,2285.8,2283.2,2280.6,2278,2275.4,2272.8,2270.2,2267.6,2264.9,2262.3,2259.
7,2257];

for j=1:101

CC(j);

end
t1=fix(t)

t2=t1+1

CC1=CC(t1+1);

CC2=CC(t2+1);

lambda=CC1-((t1-t)/(t1-t2))*(CC1-CC2);

disp('el calor de condensacion es [KJ/Kg]')

lambda

disp('la pérdida de humedad es [Kg agua evaporada/(m^2*s)]')

Nc=35.5*(tg1-t)/(lambda*1000)
 BIBLIOGRAFIA

GEANKOPIIS C J. Procesos de transporte y operaciones unitarias Tercera edición. Editorial

CECSA. México, 1998. Pág. 604