MODELAMIEN'TO MATEMATICO y SIMULACION

DE LAS OPERACIONES DE HUMIDIFICAClON

CESAR GOftlZALEZ S.ABOGAL·

Ingeniero QuiD1ico. lUSo
M.Sc.
-,
UnJversid"d Industrial de Snntander
CRISOSTOMO BARAJAS FERREIRA.· .;~

In,enJero Quirnico. Universidad Industrial'de Sant,ander

H'UMBERTO ESCALAN'fE HERNA.NDEZ" ~~

Ingeniero Qutmico. Universidad Industrial de Santander
JAIRO ADOLFO FLOREZ ALONSO·

Inpnlero Quimlco. Univer81dad Industrial de S&ntandor

RESUMEN

Este trabajo trata de la elaboración de un programa de computador para ma.

Dej~ las ecuaciones generales de disefto en torres de enfriamiento de llquido8 y
gases para cualquier valor del número de I.lewls. Una vez definido el programa
se probó su funcionamiento con sendos casos para liquido y g88 consiguiéndose
buenos resultados.

INTRODUCCION

En la8 operaciones de humidificación cuando el número de I~ewiB es 1,0 los cálculos pueden
simplificarse de tal forma que no se requieren métodos más allá de la simple integración gráfica.

Sin embargo, cuando difiere de 1,0 dado que los datos de operación no ptledel1 representarse por
medio de una linea recta, es necesario recurrir a procedimientos iterativos muy extensos para
los cuales es preciso disponer de la ayuda de U? computador.

• Departamento de InpnJerla Quimlea. Universidad Indu8trial de Santander. Bucaramanga, Colombia.

Rev. ION. Buearamanra (Colombia), 6(1):121-144, Mayo 1982

121
Este trabajo tuvo por objeto estructurar un programa de computador que permitiera ejecutar de
una manera rllpid.a, exacta y eficiente, cálculos para torres de enfriamiento tanto de l1quido como
de gases, cualquiera que fuese el valor del número de Lewis.

Se buscó que el programa permitiera establecer los valores del·area real de la torre carga mínima
y má.xima -de liquido, caudal de inundaci6Jl, diáInetro de la columna, entalpía a lo largo de la torre,
humedad relativa, perfil de temperatura, caída de presión, restricciones para evitar la niebla y
algunas otras que durante su desarrollo se irán explicando.

Al final se h.ará una prueba de su operabilidad mediante un ejemplo cuyos valores son previamente
conocidos. Por último es necesario indicar, que para todo el trabajo ,se utilizará el sistema inter­
nacional de Unidades.

LAS ECUACIONES DEL MODELO MATEMATICO DE LA HUMIDIFICACION

La Figura 1, muestra un esquema de la torre con las principales variables de cada una de las dife­
rentes corrientes que intervienen en una operación de humidificación.

Para o"per~~iones de contacto continuo adiabático en contracorriente de gas líquido y para

aquellas en las cuales el número de Lewis es diferente de 1,0 se pueden aplicar las siguientes
ecuaciones (Treybal, 1973 P. 243), las cuales se obtienen realizando balances de masa y de calor a
la intefase en ~un elemento diferencial de la torre, mostrado en la Figura 2, que conllevan a las
siguientes expresiones:

[1]

El gradiente de humedad en esta expresión se obtiene de la ecu,ación:

dY' = MA, FO'~~ 'In PT - P.A.I

[2]
dZ G'. Pr - PAO

El gradiente de temperatura del gas empleado en la ecuación [1] se toma de la ecuación:

~ =- h'G~ • (T a -T,) [3]

¡dZ G'.C'H
A menos que se conozcan separadamente aM y Se (lo CUa1110 es usual) habrá dificultad en calcular
h'o exactamente.

Por lo tanto aquí se supone que a e = a m = 8,,11'

b'o ae: = G'.CPA (dY'/dZ)/[1- EXP G'.CPA(dY'/dZ) ha Se] [4]

122 Rev. ION. Bucaramanga (Colombia), 6(1):121-144, Mayo 1982

 L2 I<Q/h-m 2 -.. G~
r ­_ _- - . . . - - . .__ lG2
t L2
Hu H'z
Vi.

L', G~ Kg. Seco / h~ m2

te1
t L1
H;
HL 1 yl,

FIGURA J Torre de humidificacion de contacto continuo ·gas - liquido

en contracorriente.
 l t dl GiS

tl + dt L te ..... dte
y' .. dyl
I
H' + dH

--­ ..... _-­ ----------­

r---­ ----­ --------r r--------­ - -------1
-----~---- ---~*---

I I I
I I F. P I\P I
1 1 1
yn
I
. ~

°IL ti
dz

m~ I
II I
I~ _
I I
I I
I I
I I
I
L
Liquido
JI Gas
I
-_o . _

t.:

FIGURA 2 Elemento diferencial de una torre de relleno para

humidicación
El gradiente de temperatura del liquido sé toma de la ecuaci6n:

dTL = G'.C PL (dY'/dZ) - H L 8 c (T,* TL)/L'C PL [5]

dZ

SOLUCION DE LAS ECUACIONES DEL MODELO MATEMATICO

El procedimiento seguido en el desarrollo del cálculo es el propuesto en el texto de «Operaciones

de Transferencia de masa») de Robert E. Treybal detallado 8- continuación para, el enfriamiento
de gases.

1) Datos principales conocidos: Número de componentes, flujo de liquido a la entrada,

flujo de gas seco, temperatura de gas a la entrada, humedad de entrada, temperatura
del gas a la salida. o --- I

2) Los cálculos se inician por la parte inferior de la columna.

3) Se supone una humedad de salida del gas, acostumbrándose como primera suposición

gas saturado. Esta humedad será verificada al final del proceso.

4) Con la humedad se calcula L 1 y T Lt haciendo balances de masa y calor.

5) En el extremo inferior de la torre, Z = 0, se supone T, y se calcula po Al' Sustituyendo

en la ecuaci6n [2] se consigue dY' IdZ.,¡

6) Debe evaluarse h' o· 8 e

7) De la ecu~ci6n [3] se calcula dTaldZ.

8) Con los resultados de 5,6~7, y de la ecuación [1] se obtiene T,; comprobándose con la
supuesta y si no se cumple, se repite, 5, 6, 7 Y8.

9) Se calcula

dY' = dY'/dZ
dTa dToldZ

10) Se elige un incremento de temperatura del gas; con este incremento y con el gradiente
calculado, supuesto constante en un pequefto intervalo, se evalúan las condiciones en
el extremo .del incremento. Para ello se sigue:'

a) Se calcula AZ =- AToldToldZ

b) Se obtiene To a la salida d.él nuevo incremento T o = T + ATo

OI

e) Se calcula Y' a la salida delincre~entocon Y'=Y't + dY' *ATo

dTo

Rev. ION. Bucaramanga (Colombia), 6(1):121-144, Mayo 1982 125

 d) Se calcula la presión del gas en el increnlento con la ecuación

PT * Y't

'y' 1 + MAllVIs

e) Se hace un balance de liquido ~T se obtiene L'. para el incremento:

f) Por un·balance de calor se obtiene T L en el incremento.

11) J?e igual'manera, se calculan nuevos gradientes dY'/dZ y dTaldZ a este 'nuevo nivel de
la torr~ y para otro intervalo.

El proceso se repite hasta que la temperatura del gas alcance la de las condiciones del gas
de salida. .

12) La humedad obtenida se compara con la supuesta en el numeral [2]. El proceso termina
cuando la humed.ad calculada sea igual a la supuesta, como también las condiciones de
liquido a la entrada sean iguales a las calculadas.

En caso de enfriamiento dell1quido el procedimiento es el siguiente:

1) Datos principales conocidos: Número de componentes, flujo dellfquido a la entrada,

temperatura de entrada del líquido, flujo de gas a la entrada, temperatura del gas a la
entrada, humedad de entrada, temperatura de salida dell1q~ido.

2) Los cálculos se ilúcian en la parte inferior de la columna.

3) Se supone una temperatura de salida del gas. Como primera aproximaci6n igual a la
temperatura pe e~trada del liquido (puede empezarse con cualquier temperatura, ya
sea mayor o menor que la de entrada del liquido). Con ésto se calcula la humedad de sa­
lida del gas como saturado.

4) Con esta humedad s~ hace un balance de masa y así se halla La­

6) ,En el extremo inferior de la torre, Z =.; O se supone T, y se calcula po Al. Sustituyendo

en la ecuación [2] se determina (dY' IdZ).

6) Se encuentra h'o* a e

7) De la ecuaci6n [3] se calcula dToldZ'

8) De la ecuación [6] se obtiene dTLldZ.

_.
. --_.­

9) Con los resultados de 5,6,7 Y de la.ecuación [1] se obtiene T, se comprueba con la su­
puesta, si no se cUlnple se repite 6, 6, 7 Y 9

126 Rev. ION. Bucaramanga (Colombia), 6(1):121-144, Mayo 1982

 dY' :d),'/dZ dY'
10) Se calculan: ----_ ... Y'
dTa dToldZ d'fl,

11) Se elige un incremento de temperatura dell1quido, con este increrrlento y con el gra­
(iiente (~a.lculado, sllpuesto COllstante en un pequefto intervalo, se evalúan. las condicio­
nes en el extremo del incremento. Para ello se sigue de la siguiente manera:

a) Se ~~;alc111a: I!Z =

b) Se obtiene TI. a la entrada del ~uevo incremento. TI. :: 1'1.1 + ATL

c) Se calc\lla Y' a la salida del incremento con

d) Se calcula la presión del gas en el incremento con la ecuación:

Pr • Y'.

e) Se hace un balance de liquido y se obtiene L' para el incremento.

f) Se calcula Toen el increnlento por medio de las siguientes ecuaciones:

ATo = AZ· dTo/dZ T o = TOI + ATo

12) De igual manera se calculan nuevos gradientes dY'/dZ, .4'fQ!~_ y dTL a un nuevo nivel

de la torre y para otro intervalo.

El proceso se repite hasta que la 1temperatura I dell1quido alcance la de las condiciones

del liquido de entrada. .

13) I.J8 h\lmedad y T o obtenidas se comparan con las supuestas en el numeral [3]. El pro­
ceso termina cuando la temperatura 1 y la humedad del gas sean iguales a las supuestas,
como también las .condiciones del liquido a la entrada sean iguales a las calculadas.

Los datos conocidos, mencionados antes, tanto para enfriamiento de liquido como para enfria­
miento de gas, son solo los que comúnmente se conocen para diseftar determinado tipo de torre de
enfriamiento; pero como el programa está en condiciones de calcular otras variables de importancia
que tienen incidencia en el disefio de la torte, entonces se requieren otros datos, los cuales se
pueden ver en el diagrama de flujo del programa principal.

REv.ION. Bucaramanga (Colombia). 6(1):121"144, Mayo 1982 127

 ENCABEZAMIENTO

ENFLIQ: r.

ENFGAS = 2.
Gs , TL •T'2.TGI , v: .
A J' L J AIRE

d t di, e , op

e PL , k , 1< L, .J{ A t ~ L • -A( V,

..AVT Pt e'

FIGURA 3 Diagrama de flujo

del programa de

CARGAM so~ucion de las ecuaclones del

(ARGUMENTOS)

modelo.
 DIAGRAMA DEL FLUJO D'EL I)ROG!t.i~MA DE SOLtJCION D}tJ

L.AS ECUACIONES DEL MODELO

El progranla consta de una linea principal, la cual tonl8. la .decisión sobre elllam.adú de laa rlife..
rerites rutinas que componen el resto del programa. En este arti.culo solo se presellta esta linea
principal, indlc~da en la Figura 3e

En total el programa consta de 8 rutinas-, las cuales en su detalle se t.ratan a continuación:

I~a Rutina GAS que calcula, para. un determinado delta de temperatura de gas; la humedad
m,olar, el flujo de líquido, la temperatura de liquido, la variac~6n d~ la humed~d.con la tempera­
tura de gas, entalpfa de la mezcla gaseosa. la temperatura de Jntertase, el coefiCiente de trallsfe­
rencias de masa y de calor para el gas y ell1quido t el área interfacial, la altura de relleno, el área y
diá.metro de la columna, la humedad relativa.

I-Aa Rutina FlJNCTION PV funciona cada vez que se necesita una presión para una temperatura
dada. Solo requiere 'de la tempernttlra como argumentoe Es llamada por las rutinas GAS y LI­
QUID. '. .

CAIPRE es una rutina que calcula la caída de presi6n seca y la caid'a de presión total 8 lo largo
de la torre. Está realizada en base a los métodos de Chilton-Colburns y ~Y9~~QY(Sawisto\\rskit
1967, p.41-46).

El flujo mínimo y máximo de liquido lo calcula la rutina CARGAM basado en el tipo de empaque,
flujo de gas y flujo de liquido. Esta rutina fue programada en base a los datos de Monis y Jackson
(Sawistowski, 1967, p. 15).

Para la Rutina LIQUID se tiene un procedimiento de cálculo similar a la rutina GAS, pero con
una diferencia fun.damentsl, la cual consiste en que la rutina de enfriamiento de un gas no tierle
rutinas internas. La rutina de enfriamiento de un liquido consta de las rutinas internas: VISDD,
ARINTF y FG HG HL. Estas rutinas en ·su orden calculan: la viscosidad de la mezcla gaseosa
como un promedio entre el fondo y el tope; área interfacial y fracción de vacío; y el coeficiente
de transferencia de masa, la transferencia de calor de la fase gaseosa y la transferencia de calor
de la fase liquida.. ~a rutina LIQUID, pa~a cada delta de temperatura de liquido, calcula: la hume­
dad molar, el flujo dell1quido, la temperatura de gas, el gradiente de humedad, la entalpía la altura
de relleno, los coeficiente de transferencia de masa y de calor, el área interfacial y la temperatura
de interfase. Además el área, el diámetro de la columna y la humedad relativa.

PECUAJ..IIARIDADES DEL PROGRAMA DESARROLLADO

Este apartado presenta los elementos necesarios que se deben conocer para manejar este pro­
grama. Se empieza por dar las' especificaciones de entrada, o sea la forma de suministrar los datos.

Rev.ION. Bucaramanp(Colombia). 6(1):121-144. Mayo 1982

129
En primera instancia para enfriamiento de liquido:

FORMATO DATO UNIDADES FORMATO DATO UNIDADES

PRIMERA TARJETA

12 N Adimensional F8.3 LDO kJ/kg

F8.2 1~2 kg/h F7.2 RODLIQ kg/m J
F8.2 G. kg/h F7.4 ROAIRE kg/m 3
':f:" TLJ oC Fó.2 PMLlQ kg/kmol
F6.2 TL.I oC F6.2 TO.I oC
F6.4 YMll kg ,vapor/kg
gas seco

F8.6 DABO m 2/h

SEGUNDA TARJETA

F6.4 CPL kJlkg oc

F7.4 K kJ/kg oc

F7.4 KPL kJ/h.m. oc

F8.4 UA kg/m.h:,

F.8.4 UL kg/m.h

F8.4 UVF kg/m.h

F8.4 UVT kg/m.h

F7.3 PT kPa

F6.4 CHLIQ kJ/h.m. OC

TERCERA TARJETA

F7.4 PM(~) kg/kmol F8.4 UGS(2) kg/m.h

F7.4 CP(2) kJ/kgOC F7.4 PM(N) kg/kmol
F7.4 W(2) kg/kg tot gas F7.4 CP(N) kJ/kgOC
F8.4 UGE(2) kg/m.h F7.4 W(N) kg/kg' tot gas
F8.4 UGE(N) kg/m.h
F8.4 UGS(N) kg/m.h

ULTIMA TARJETA

F6.4 DS m
F6.4 DR m
F4.1
DRiN pulg*
F5.3 EPSILO m 3 /m.1

F6.2 APRIMA m 2/mJ

130 Rev.ION. Bucaramanp(Colombla), 6(U:121·144, Mayo 1982

 p o o
o 2 o
ii o N
UI

."
@ Hum.dad ICG vapor I KG QOs seco
:o
J> ~
~

rn
~
~.
Q

~.
CD
.....
::J
o
Q.
CD

s::
::J
_.
.a
e
o.:
~

~ ~
3 (JI
'O

..a
ct

e
a
o
o

r
¡ .
t
 O
0.0360-1 U
1&1
(1)

U» 172
44..J o· 034O i ~ j/ O
u
'"
en

/0 //Í
<:)

42..1
0.0320
1 ~ca: 3
o m
O.03~ ~
~
/00# y": f (z) _ ....______
c:J

~
49~ t-t= f (z) ~
o
~ TL=f(z) _ _. _ ""i2
.'l'~#'~
38J~ 0C280i Q
TG= fez) ___ I
8w ~
..;,J
~ 0.0260 /..
~ t!
36 J ~ :;)
. /"~ ,/" z
~ % lIJ

34J ~
a: 0.0240 /' /" .
"
//"
/

2
tIJ
f­ 0.0220
/ ",¿oo //
32 -1 / o" ,,//
,. //
/
/
0.020°1 ,/
30J ~~
,~~

O.0180J ~/

/~/'
27 J 0.0165 J i <" i i r

o I 2 3 4 5 6 7 8
'z' M ALTURA DE RELLENO
FIGURA 5 Variacion de las temperaturas del gas y del liquido 1 entalpia y h.umedad del gas con la
alfura de la columna en el enfriamiento de un liquido.
 L2 = 5eooo.00 ka/HORA
Tl2 = 44.00 CENTGRA ~ 881: 18000.00 ,.11 HORA
TG2 = 42.87 CENTGRA
VM2 e 0.033884 1<8 (V) IICI (1 SECO)

Z = 4.12 NETAOI.
AREA = 5. 180 NE7ROS . . . 2
OC = 2.6844752 METROS.

L I • 84 s.e. 54 ICG HORA

TLI • 30.00 CENTGRA
L VN.aO.OIUOO ICG {VII K8
TGI = 50.00 CENTGRA.
leSeCOI

ca PRIMA • 1180.27 ICG HORA *

N •• 2
HUMEDAD RELATIVA 1: eO.80 POR CIENTO.

DELTAP : O. 361 KPA DElTPS : 0.267 KPA

L MIN = 7305.85 KG/HORA * M ... 2
L MAX= 12947.65 KG/HORA .. M ... 2
L PRIMA =9681.23 KG/HORA I M** 2
CARGA DE LIQUIDO DENTRO Oe: lOS LIMITES

FIGURA 6 Resultado del primer caso

I Enfriamiento de un liquido'
 l48 rutina de enfriamiento de gas 8010 difiere en la primera tarjeta d.e la rutina IJIQUID~ La forma
de d,ar loa datos es la siguiente: . .

F"ORl'..lATO ' DA'l'O UNIDADES FORMATO.,. DAfrO UNIDADES

12 N Adimensional F6.2 YMI1 kg vapor/kg gas

F8.2 L2 kg/h F6.4 DI"\BO nl~/h

F8.2 G.r kg/h F8.3 LDO k~/~~..

F6.2 TL2 oC F7.3 RODLIQ kg/m J

F6.2 TG 2 oC F8.4 R()AIRE kg/m 3

F6.2 TGII oC

i~l programa principal está est~cturadode tal manera que puede trabajar s6lo para enfriamiento
de líquido, o s610 enfriamiento de gas, o los dos casos simultárleos. Quien desee trabajar con las
dos rutinas lo deja tal y como está en ést.e trabajo, pero quien quiera dise:ftar una torre de enfria­
miellto de líqui.do solamente, entonces cambia la instrucción ENFGAS = 2 por ENFGAS = O
lo miSlno sucede para el otro caso.- s610 que en lugar de cambial· la tarjeta de GAS ca~bia
ENFLIQ = 1 por ENFLIQ = o.

PRUEBA DEL PROGRAMA DESARROLLADO

El programa fue probado con los siguientes casos:

1) Una planta requiere 56000 kg de agua/hora para circular a través de los condensadores a 440C
se plantea enfriarla para volverla a usar, por contacto con aire en una torre de enfriamiento. '

Las condiciones elegidas para el proyecto son: la temperatura de bulbo seco 30°C, la humedad
de entrada YMl = 0.0165 kg vapor de agua/kg de aire seco, flujo de gas seco 18000 kg/h,· l
el agua se requiere a 30°C, el tipo de relleno es anillos Raschig de 0,0508 m (2 pulgadas).

Para el desarrollo ~del programa se presenta el siguiente listado con las variables indispensables
para su funcionmniento , de acuerdo al lenguaje utilizado en el programa. Además de las variables
que aparecen en el enunciado, se reportan los demás datos de entrada para el problema, las cuales
han sido tomadas de las referencias dadas;

N 2 UVT 0,0396
L2 55000 PT 101,.325
GS 18000 CHLIQ 1,8723
TL1 44,00 PM(l) 28,000 N 1 =1 0 1 =2
TL¡ 30,00 CP(l) 1,0383
YM11 0.0165 W(l) 0,790
DABO 0,001347 UGE(1) 0,0655
LDO 2500,776 U08(1) 0,06·88

134 REv. ION. Bucafamanga (Colornblat. 6(1):121·144, Mayo 1982

_ llIIPl'"7.nll_msrtttu' "11" I'liTUrt;'

RODLIQ 993,55 PM(2) 32,00

ItOAIRE 1,2062 CP(2) 0,8424
PMI.lIQ 18,02 W(2) 0,2100
T(}}I 30,00 UOE(2) 0,0720
CPL 1,1868 U08(2) 0,0756
K 0,0963 ns 0,0726
KPL 2,2609 DR 0,0508
UA O,033(j DRIN 2,0
U1.1 ~,5290 EPSILO 0,740
l.JVF 0,Q3(;)0 AP'RIMA 91,84

Como resultado de este caso se presenta en la Figura 4 la variación de la humeda.d de la mezcla

gaseosa con la temperatura del gas. Se observa que la humedad relativa está bastante lejos de
la saturación lo que corrobora el hecho que para lograr lIn buen grado de separación -seria el caso
de una torre de humidificación- se necesita. una torre de bastante altura, por lo tanto al diseftar
wm torre lo más importante es el delta de temperatura útil capaz de hacer el enfriamiento ·deseado.
La Figura 6 indica los perfiles de temperatura tanto del liquido como del gas a través de toda
la torre. De igual manera se presenta la variación, en toda la sección de la torre, de la entalpía
y la humedad del gas. En la Fi.gura 6 se presenta un resumen de los resultados conseguido~ para
el disefto de esta torre de enfriamiento de agua.

Un generador de gas produce 65% de Nitrógeno y 350/0 de monóxido de Carbono a 101,325 kPa
. y 315°C. El gas inicialmente seco será enfriado hasta 27°C por contaco en contracorriente con
agua a 18°~~yeevaa usar una torre con relleno' irregular usando anillos Raschig de 0,0508 m
(2 =
pulgadas) con G', = 230 kg/hora y del Uql1ido L2 "1000 kg/h. Calculm las condiciones de
la torre.

Para el desarrollo del programa se presenta el siguiente li~tadQ I con las variables indispensables
para su funcionamiento, de acuerdD allen~8ie utilizado en el programa. Además de las variables
que aparecen en el enunciado, se reportan los demás datos en entrada para el problema, las cuales han
sido tomadas de las referencias dadas. - ,

N 2 p'r 101,325
L 1 1000 CHLIQ 1,8723
OS 230 PM(l) 28N]=1 CO=2
TL 1 18 CP(l) 1,0383
~G2 27 W(l) 0.65
TGII 315 UGE(l) 0,1051
YMII 0,00 UGS(1) 0,0641
DABO 0,0188 PM(2) 28
LDO 2500,776 CP(2) 1,0666
RODLIQ 1000 W(2) 0,35
ROAIRE' 1.2062 UGE(2) 0,1044
CPI~ 4,1868 UGS(2) 0,0648
K 0,0427 ns 0,0726
KPL 2,2176 DR 0,0508
UA 0,0336 DRIN 2,0
UL 0.9 EPSILO 0,74
UVF 0,00 APRIMA 91,84
UVT 0,0356

Rev. ION. Bucaramang. {Colombia), aU):121-144, M.yo 1182 135

Al igual que el caso anterior para el en-friamiento de la mezcla gaseosa se presenta en las Figu­
ras 7, 8 Y 9 las mismas variaciones a través de la columna asi como los valores de disefto. Es con­
veniente hacer notar en la Figura 7, como a medida que se asciende en la torre el gas se enfrfa
y humidifica hasta llegar a un limite donde el potencial se invierte y entonces el gas comienza a
dehumidificarse.

CONCLUSIONES

Del presente trabajo se pueden obtener algunas observaciones importantes:

a) Este programa pr~senta una precisión bastan.te buena respecto a la prueba de ensayo
con la cual fue programado.

b) Emplea el método iterativo en el proceso de cálculo. Por este nlotivo su convergencia

es rápida, su tiempo corto y por ende-económico.

c) La utilización de ecuaciones eJnpfri~as para la cuantificaci6~ de la mayoria de las

variables que intervienen en este programa y el cálculo de otras Ilor ecuaciones de
balance, hacen del programa un elemento básico para la comparación de las diferentes
variables para cada intervalo y una guia fundamental para los estudiosos de esta clase
de operaciones unitarias.

NOMENCLATURA DEL PROGRAMA

LENGUAJE LENGUAJE
PROGRAMA INGENIERIL

AAW Superficie interfacial específica para absor­

ción, desorción y destilación, m 2 /m J
APRIMA Superficie especifica de rellello, m 2 /m 3
AREA A Area de la columna, m 2
AVW Superficie específica interfacial para el con­
tacto de un gas con un líquido puro, m 2 /m 3
BETA (J Constante utilizado para el cálculo de 'tw
B b Constante utilizado para el cálculo de +.rw en
la rutina CAIPRE.
CH C'H Calor húmedo especifico kJ Iltg o e
CHLIQ C'HL Calor especifico del vapor de liquido, kJIkg o e
CMH Caudal mínimo de humectaci611, m 3 /h.m
CP Capacidad calorífica del gas, kJ/kg oC
CPL Capacidad calorífica del líquido, kJIkg °e

136 Rev. ION. Buearamanp (Colombia), 6(1):121-144, Mayo 1982

 o o o o
.b b b
8 o
CJ'
"O N
O
o

(Jo'
o
1\

en
o
~

~J>
--..,
c.o
o

rr1
~
Q
Ñ
'3 o
éi·
::J
o
a. <:ñ
CD o

e:;:,

¡¡
e o
Q
ti)

~
o

..
I\l
o
 o
u
••
!
-o y= f( z)
u
<:)
.­g. L400
400J~
o
0°·020
'a
x
~
o
Q.
HG=f(z) _0 ___-
•a a e
Ql ! >
C)
T =1(z) -
L
____
x
.,
"-
o....

...
o ~ 3JC

. ...
:::1
'6
=
e;

"o T =f{z) _ .. _ _ 0_ ea. .

.,
• ~ G
..,¡
4)

300J E" ~O.015 ~

I 1-300
~ ~

,~.,,\.
%

30J I .
,'\ ~.

200J I ------­
j 0.010
"
'. ,"",
,",
00

25
.. ~~. '"
100
" ". ,
",~ "­ 100

­ - ­ -----­
0.0 1

-~
'""
.~
20 ------

0.0
18
0.00
.~.

'---, _ _ _._
'--"-- '-- '-­
­ --­
......... ~ ..
-....;;;;
.. .. ..

0.0 0.5 1.0 1.5 2.0 2.5

Aatura 'z' M
de retleno,

FIGURA 8 Variocion de las temperaturas del gas y ~eI liquido, entalpia y humedad "del

gas con la altura de Ia forre en el enfriamiento de gas.

.......". ~._ .. .....,~-..,..:.:-~ ..'..,..._n.•_:"~~ ........~.. c-~.:::-~",'::;;;==----,,:~_,_~.~~:..~,~~~-, ~~-.. ~.~ - ~~:c;. .;':':~',,:;;::;_...-~', ..,,,,,~~;.-.>-¡;:;:':-.~' -F ~':'-:"'';-~''~~-H -..L~-¡:~3f':::::.z::;.~,:.~:··,'"-~.;:>~·~::;·,·,,--;··~ .,::""""<0" ~~~;-~~- ,..' :..:..'- )::-];,.?;;."'"-:"'1- ._'- ";"""~~iñiñ'-;--~~ ,-.._- ::;.=-~

L2 • toOO.OO kG / HORA es: 230.00 KG/HORA

'rLa. ".00 CENTeRA 'Te2 = 2'7.00 CENTEGRA
YM2 e 0.014514 k8 (V) I 1('8 (CI :SECO)

Z;: 2. 18 N ETROS
AREA 11 O. 104 METROS • • 2

DC • O. 1810397 NETROS

LI •
TLI
Itl.98
12.•1
K8 I HORA
CEHT8RA
LVNI ·
TG'
0.0000
a15.00
Ke (V J/ K8(e SECO
CENTGRA
I

GPRIMA 2221. 93 Ke/ HORA M... 2

HUMEDAD RELAtIVA 62.28 POR CIENTO.

OELTAP 0.112 J<PA

DELTAPS 0.083 kPA
L MIN 7337.87 KI/HORA M•• 2
L MAX 15955.6S KG/HORA M**2
L PRIMA 9645.89 kG/ HORA M *. 2
CARGA DE LIQUIDO DENTRO DE LOS LIMITES'

FIGURA 9 Resultado del segundo caso·

1 Enfriamiento de un gas
CPROM CP Capacidad calorlfica prornedio del gas,

kJ/kg oC

DAB Difusividad nlolecular, m 1(he

DABO Difusividad nlo1ecular a O°C y 101,326

kPam:l/h

De de Diámetro de la columna" m

DELTP 4Pia Caída de' presión húmeda., 'kPa

DELTAZ AZ Diferencial de altura que corresponde a cada

delta de temperatura, m~~

DELTPS Caída de presión seca, kPa

DELTL Incremento de tem.peratura de líquido tomado _

como Op5 en la rutina de Uquido y en la de

gas como igual a dTLldZ*AZ, oC ~
DELTG ATo Incremento de temperatura de ga.s tomando

como ·2 ell la rutina de gas y en la liquido

como igu.al a dTGldZ*AZ, oc.

DR da Diámetro del anillo, m

DS d. Diámetro de una esfera de igual superficie

que un solo elemento de relleno, m,\

DTLENZ dTL/dZ Gradiente de la variación de temperatura de

líquido con respecto a la altura, OC/m

DYENDZ dY'dZ Gradiente de la variaci6n de la humedad del

gas con respecto a la altura, kg . vapor/,kg l

gas seco/m

DYENTL dY'/dTL Gradiente de la va~iación de la hUlnedad

del gas con respecto a la temperatura del

liquido, kg vapor/kg gas sec%C

DTGENZ dY'/dZ Gradiente de la variación de la temperatura

de gas con respecto a la variación de altura,

oC/m,·

ENTALP H Entalpía de la mezcla gaseosa, kJ /kg

EPSIL tO Volumen de vacío, fraccionario, en· rellenos

irrigados, m 3 /m 3

EPSILO Volumen de vacios, fraccionario, en un lecho

relleno seco, m 3 /m 3
FA Fa Factor de corrección
FG Fo CoeJiciente de transferencia de masa,
kmol/h.m 1
FIOW Fracción de volunlen de líquido en operación
para agua, m J /m 3 ,
FISW Fracción de volumen de l1quido estático
para agu.a, m 3 /m 3

FITE Retención total de líquido, rn J /m 3

140 Rev. ION. Bucarl'.manga fColombia), 6(1):121-144, Mayo 1982

~.-.

FITEW +,. l¡'faetor de corrección

F, Factor de corrección

FS
f F8ct.~r de rozamiento u1odificado.

FPRIMA
G Velocidad superficial molar del gas,

G
kmol/h.nl1

G'iNU," Caudal de inundación para el gas, kg/h.m1

GINUM
G'. Velocidad superficial de la masa de gas seco,

OPRIMA
kg/h.m 2

GPROM G Velocidad de la mezcla gas seco-vapor pro­

medio entre el fondo y el tope, kg/h.m1
OS
G. Flujo de gas seco kg/h,,}
HG
hG Coefi.ciente de transferencia de calor para el
. gas, kJ/h.m 2 oC
HL hL Coeficiente de transferencia de calor para el
liquido, kJ/h.m 2 oC I

HPGAM h'oaM Coeficiente corregido de transferencia calóri­

ca para el caso de transferencia simultánea
de masa por umdad de volumen kJ/h.m 3 0C.
HR
H. Porcentaje de :hum~dad l relativa
K
k Conductividad térmica para la fase gaseosa,
kJ/h.m.oC
KPL k'L Conductividad térnuca para la fase liquida,
kJJh.m.oC
L L Flujo promedio de líquido entre el fondo y el
tope, kg/h.

LDO lo Calor latente del liquido a tempertura de

referencia, kJl.kg~

LH L"
Caudal de humectación, m 3 /h.m

LMAX L".,...
Carga máxima de líquido, kg/h.m 2

LMIN L""N Carga mínima de liquido, kg Ih.m 1

LPRIMA L' Flujo promedio de líquido por unidad de

área, kg Ih.m l

M
m Constante empírica""

N
N N~ de componentes en el gas seco

NZ
n Constante empfrica
p
p Constancia empfrica
PA pO
A Presión parcial del vapor en la mezcla ga..
seosa, kPa
PI n Nllmero adimensional usado para el cálculo
de AP".
PM M Peso molecular del gas, kg fkmoll
,PMLIQ Mil. Peso molecular dell1quido, kgrkmol

REv. ION. Buearamanp (Colombia), e(1):111·1..... Mayo 1982 141

PMPROM M lleso nlolecular promedio para el gas seco,

Kg/kmol¡

PRANDL Prl, Número de Prandtl, lJara el liquido, adilnen­

sionaJt

Número'd.e'"Prandtl para el ~as, adimensional

PS pO Presión de saturaciÓn, kPa

PT PT Presi6n total igual. a 101,325 kPa

-- --'

PV pO Presión de vapor en la rutina PV(TG), kPa

RAGLIQ Fa(Vo/VL ) Rela.ción de velocidad de gas a velocidad de
liquido por un factol- de correccióDt adimen­
8ional~

RE Re Núnlero de Reynolds para el gas, adimensio­

nal.
liEL ReL Núnlero de Reynolds para.ellíquido, adimen­
'J sional
ROAIRE Qaire Denªidad del aire a 20°C y 101,325 kPa
~k.g/mJ .
ROATI Qti Densidad de gas a temperatura de interfase,
kg'lm 3
RODGÁL\.S 80 Densidad del gas a temperatura de opera­
ción, kg'/m 3
ROGPRO Qa Densidad del gas promedio entre el fondo y
el tope, kg 1m 3
SCMIDT Se NÚlnero de Schmidt, adimensiona1.
TG To Temperatura del gas, oC
TO To Temperatura de referencia, tomada como O°C
TL TL Temperatura del liquido, oC
TI TI Telnperatura de interfase, oC
DA IlÁ Viscosidad del agua a 20°C, .kglm.h·
UG (lo Viscosidad de la mezcla gaseosa, calculada
como un promedio entre el fondo y el tope,
kg:/m.h",
UGE #loE viscosidad del gas seco a la entrada de la
torre, :kg Im.h'Q
uas IJos Viscosidad del gas seco a la salida de la
torre ltg.'/m .h ,.
UL JAL Viscosidad del liquido a la temperatura de
proceso, ,kg/h.m,~
VG vo Caudal volUInétrico, m 3 /h
VGAVL VO/V L Relación de caudales de gas a liquido, adi­
mensional.

VL
VL Caudal volumétrico, m3/ht~

XM
Xm Fracción molar de la mezcla gases secos

Rev. ION. 8ucaramo.nga (Colombia), 6(1):121-144, Mayo 1982

142
w w Fracción en peso de los componentes de
la mezcla gaseosa
YM Y' Humedad másica, kg vapor/kg gas seco.
Z z Altura total de la torre, m

OTRAS VARIABLES UTII..IZADAS EN EL PROGRAMA

CAMBTL o. Esto significa que cuando el programa haya

funcionado y cumplido el requisito de hunle­
dad y no cumpla por temperatura de liquido,
entonces suponga la misma humedad pero
I

TLl 'sea igual a rrLII+PARETL y haga

todo el proceso ,~tra vez.
ITER Número de iteraciones.
PARETG Criterio de pare por temperatura de gas.
PAREYM Criterio de pare p~r humedad
PAREL2 Criterio de pare por liquido
PARETL Criterio de pare por temperatura de 1iquido~
CONST Se llama const al término'
d,O' ) = 0,36
( 0(1- +)
ABSTI Criterio de comparaci6n de la T, supuesta,
con la calculada.

SUBINDICES

F Final de cada incremento, para T, Y', L. Fondo para el resto de cálculos.

I Inicial. Al entrar a la torre solamente.

1 Entrada a la torre o a cada incremento.

2 Salida de la torre.

T Tope.

AB8TRACT

This work sbows the elaboratioD of a computer program for banling tbe general
desigD equatioDs of cooling towers for any vahie 01 Lewia Dumber. Tbe system
W88 tested wit~ case problems and good results were obtained.

Rev.ION. Bucaramanga (Colombia), 6(1):121-144, Mayo 1982

143
 BIBLIOGRAFli\

1) AIIUMADA, N."y GOMEZ, G. Trabajo de Dlsefto 11. UIS. 1979.

2) BARAHONA, A. y BARAHONA, F. Metodologia de trabajos cientificos. 28. ed., Bogotá IPLER, 1979.

3) BARAJAS, C., ESCALANTE, H. y FLOREZ. J. Operaciones de humidificación con Lewis diferente de 1,0. Trabajo

Disefio 11. UIS, 1981.

4) BROWN, 0.0. y Colaboradores. Operaciones Básicas de la Ingenierla Quimica. 3a. OO., Barcelona, Mmn, 1966.

6) FOUST, A.S. y Colaboradores. Principios de Operaciones Unitarias. 11a. OO., Mé~co. C.E.C.S.A., 1978.

6) McCABE. W.L. y SMITH, J. Operaciones básic8s de la Ingenierla Química. 58" oo., Barcelona, Reverte, 1975.

7) ORGANICK, E.I. Fortran IV. Mass8chusetts, ~~o Educativo Interamericano, 1972.

8) SAWISTOWSKI, H. y SMITH, W. Métodos de célculo en 108 procesos de transferencia de materia. Madrid, Alhambra,
1967.

9) TREYBAL, R.E. Operaciones con Transferencia de Masa. 2a. oo.,: BUeD08 Aires, Hasft, 1973.

Firmado digitalmente por Crisóstomo Barajas Ferreira

Nombre de reconocimiento (DN): cn=Crisóstomo Barajas Ferreira,
o=Universidad Industrial de Santander, ou=Escuela de Ingeniería Química,
[email protected], c=CO
Motivo: Soy el autor de este documento
Fecha: 2008.12.11 16:17:18 -05'00'

144 Rev.ION. Bucaramanga (Colombia). 6(1):121-144. Mayo 1982