Tesis Carlos Durán

ESCUELA SUPERIOR POLITCNICA DEL LITORAL
Facultad de Ingeniera en Mecnica y Ciencias de la

Produccin
Desarrollo de una Tcnica Lagrangiana para el Dimensionamiento

y Anlisis Paramtrico de Secadores Rotativos.
TESIS DE GRADO
Previa a la obtencin del Ttulo de:
INGENIERO MECNICO
Presentada por:
Carlos Geovanny Durn Oscuez
GUAYAQUIL ECUADOR
2014
AGRADECIMIENTO
A Dios y a mi familia.
Al Ing. Vicente Adum Gilbert,
por su colaboracin en la
realizacin de este trabajo.

DEDICATORIA
A MIS PADRES Y
HERMANOS
TRIBUNAL DE GRADUACIN
________________________ ________________________
Dr. Kleber Barcia V., Ph.D. Ing. Vicente Adum G., M.Sc.
DECANO DE LA FIMCP DIRECTOR
PRESIDENTE
_________________________
Dr. Guillermo Soriano I., Ph.D.
VOCAL
DECLARACIN EXPRESA
La responsabilidad del contenido de esta
Tesis de Grado, me corresponde
exclusivamente; y el patrimonio intelectual
de la misma a la ESCUELA SUPERIOR
DEL LITORAL
(Reglamento de Graduacin de la ESPOL)
________________________
Carlos Geovanny Durn Oscuez
II
RESUMEN
En este trabajo se pretende probar la hiptesis de que es factible analizar
secadores rotativos desde una perspectiva lagrangiana, la cual combine la
transferencia de calor y masa con el movimiento de una partcula individual a
travs del secador.
Un secador rotativo es un cilindro metlico giratorio dispuesto de forma
horizontal, con o sin inclinacin, por el cual pasa aire caliente para extraer
humedad del producto. En el interior del cilindro se encuentran montadas
unas paletas que recogen el material granular desde el fondo y lo descargan
hacia la corriente de aire para formar una cascada de partculas.
En la primera parte del estudio, se desarrolla un anlisis termodinmico del
secador para obtener el flujo msico y las propiedades del aire requeridas
para secar cierta cantidad de producto hmedo. Luego se establece un
ngulo de cada representativo de toda la cortina de slidos, con el cual se
estima la altura promedio de cada. A partir de all, el anlisis se concentra en
simular el movimiento de una sola partcula que cae desde dicha altura hacia
el fondo del tambor, donde es nuevamente recogida por una paleta para
repetir el ciclo de levantamiento y cada a lo largo del secador.
Durante su recorrido, se determina la tasa de transferencia de calor y masa
en pasos de tiempo definidos, registrando las variaciones en las propiedades

III
de la partcula y del aire, y sumando los desplazamientos sucesivos de la
partcula hasta obtener la longitud de tambor correspondiente al nivel de
humedad deseado en el producto.
Para medir la validez de este planteamiento, se dise un software que
permiti aplicar la metodologa propuesta a dos secadores rotativos
previamente analizados con otros modelos existentes, uno para partculas de
madera y otro para concentrados de zinc. De este modo, se compararon los
datos tericos obtenidos con el modelo propuesto contra aquellos de los
modelos existentes, y contra la informacin experimental disponible.
Los resultados indican una diferencia relativa de 6.36% entre las estimaciones
del modelo propuesto y las del modelo existente para el contenido de
humedad final en la madera. Asimismo, se encontr un buen ajuste entre el
perfil de humedad terico generado por el modelo propuesto y los datos
experimentales de humedad para concentrados de zinc, en funcin de la
posicin axial en el secador, obtenindose un error relativo de 4.62% en el
contenido de humedad final del producto.

IV
NDICE GENERAL
Pg.
RESUMEN ...................................................................................................... II
ABREVIATURAS .......................................................................................... VII
SIMBOLOGA .............................................................................................. VIII
NDICE DE FIGURAS .................................................................................... XI
NCIDE DE TABLAS ....................................................................................XIV
INTRODUCCIN ..........................................................................................XV
CAPTULO 1
1. GENERALIDADES.............................................................................. 18
1.1 Caracterizacin de los slidos granulares .................................... 18
1.2 Fundamentos en el secado de slidos ......................................... 22
1.3 Transporte de lquido en los slidos ............................................. 27
1.4 Equipos de secado ....................................................................... 30
1.5 Caractersticas del secador rotativo ............................................. 32
CAPTULO 2
2. MTODOLOGAS DE DISEO DE SECADORES ROTATIVOS ........ 46
2.1 Diseo mediante correlaciones disponibles.................................. 47
2.2 Diseo basado en pruebas de planta ........................................... 50
2.3 Avances en el estudio de secadores rotativos .............................. 51
2.4 Metodologa propuesta ................................................................. 57
CAPTULO 3
3. ANLISIS TERMODINMICO ............................................................ 60
3.1 Calentador .................................................................................... 60
3.2 Secador ........................................................................................ 64
3.3 Determinacin del flujo msico del aire de secado ....................... 67
V
CAPTULO 4
4. MOVIMIENTO DE LOS SLIDOS EN EL SECADOR ......................... 72
4.1 Cada de las partculas desde las paletas .................................... 74
4.2 ngulo representativo y altura promedio de cada ....................... 88
4.3 Suposiciones sobre la partcula rastreada .................................... 92
4.4 Movimiento longitudinal ................................................................ 97
4.4.1 Flujo paralelo ........................................................................... 100
4.4.2 Contraflujo ............................................................................... 105
CAPTULO 5
5. TRANSFERENCIA DE CALOR Y DE MASA ..................................... 109
5.1 Suposiciones generales ............................................................. 110
5.2 Capa lmite y grupos adimensionales de importancia ................ 113
5.3 Transferencia de calor ................................................................ 114
5.4 Transferencia de masa ............................................................... 116
5.5 Analoga de la transferencia de calor y masa ............................. 122
5.6 Enfriamiento evaporativo ............................................................ 124
CAPTULO 6
6. IMPLEMENTACIN DE SOFTWARE DE ANLISIS DE SECADORES
ROTATIVOS ............................................................................................ 127
6.1 EES (Engineering Equations Solver) .......................................... 127
6.2 Anlisis termodinmico ............................................................... 130
6.3 Masa de producto en las paletas en funcin del tiempo ............. 133
6.4 ngulo representativo y altura promedio de cada ..................... 142
6.5 Secado en flujo paralelo ............................................................. 145
6.6 Secado en contraflujo ................................................................. 160
CAPTULO 7
7. RESULTADOS Y ANLISIS .............................................................. 163
7.1 Parmetros de operacin y dimensiones de secadores tipo ...... 163
7.2 Resultados obtenidos en EES .................................................... 177
7.3 Comparaciones .......................................................................... 188
7.4 Anlisis sobre el modelo lagrangiano ......................................... 196
VI
CAPTULO 8
8. CONCLUSIONES Y RECOMENDACIONES..................................... 208
APNDICES
REFERENCIAS BIBLIOGRFICAS
VII
ABREVIATURAS
cte. Constante
CFD Dinmica de Fluidos Computacional
EES Engineering Equations Solver
FCA Factor de Correccin de rea
PI Informacin de Proceso
SCLR Sistema de Coordenadas Local Rotatorio
SCG Sistema de Coordenadas Global No Rotatorio
gPROMS Software de Modelado de Procesos
cm Centmetro
m Metro
mm Milmetro
ft Pie
hr Hora
min Minuto
s Segundo
C Grado Centgrado
F Grado Fahrenheit
K Kelvin
J Joule
BTU Unidad Trmica Britnica
Grado
CFM Pies cbicos por minuto
rpm Revoluciones por minuto
atm Atmsferas
VIII
SIMBOLOGA
Aceleracin de la gravedad m/s2

Agua evaporada desde una partcula kg
Altura promedio de cada m
aciado ngulo de vaciado , rad
ngulo dinmico de reposo , rad
A ngulo entre segmentos de paleta , rad
ngulo representativo de la cada de , rad
slidos
A rea m2
Calor agregado W
cp Calor especfico J/kg-K
Calor latente de vaporizacin J/kg
E Cambio de energa cintica J
E Cambio de energa potencial J
2
DA Coeficiente de difusin binario m /s
Coeficiente de friccin dinmico
Coeficiente global de transferencia de W/m3-K
a
calor volumtrico
Coeficiente promedio de transferencia de W/m2-K
calor
Coeficiente promedio de transferencia de m/s
masa
D Coeficiente de arrastre
Conductividad trmica W/m-K
Contenido de humedad en el equilibrio
bs (Base Seca)
Contenido de humedad final (Base
Hmeda)
Contenido de humedad inicial (Base
Hmeda)
Desalineamiento angular con respecto a la , rad
lnea radial que va desde el centro del
tambor a la raz de la paleta.
b l Densidad al granel o aparente kg/m3
s Densidad de grano seco kg/m3
Desplazamiento de la partcula en el i- m
i
simo ciclo
D Dimetro del tambor m
D rano Dimetro promedio de grano o partcula m
IX
RO Distancia desde el centro del tambor hasta m

la punta del labio de la paleta
Entalpa J/kg
e Espaciamiento entre paletas , rad
Flujo msico de agua kg/s
Flujo msico de aire (mezcla aire-vapor de kg/s
aire
agua)
d Flujo msico de aire seco kg/s
p Flujo msico de producto hmedo kg/s
s Flujo msico de producto seco kg/s
Flujo msico de vapor de agua kg/s
part c las Flujo total de partculas UND/s
Fraccin de producto que ha cado desde
i
la paleta
Fuerzas de arrastre en las direcciones Z y N
D D
Y
Humedad absoluta del aire kg agua/kg
aire seco
, HR Humedad relativa del aire
Inclinacin del tambor
l Longitud de la base de la paleta m
l Longitud del labio m
L Longitud del tambor m
Masa kg
s Masa de grano seco kg
rano Masa del grano hmedo kg
n Nmero de ciclos de secado
Le Nmero de Lewis
Nmero de Nusselt
r Nmero de Prandtl
Re Nmero de Reynolds
Posicin angular de la paleta , rad
d Profundidad o ancho de paleta m
R Radio del tambor m
Tasa de evaporacin de agua desde una kg/s
nA
partcula
Tasa total de agua evaporada kg/s
T Temperatura C
t Tiempo s
tres Tiempo de residencia de la partcula min
tsecado Tiempo de secado s
Tiempo de vuelo de la partcula en el i- s
ti
simo ciclo
X
Trabajo mecnico J
Velocidad m/s
rpm,
Velocidad angular del tambor rad/s
Velocidad msica del gas kg/s- m2
Viscosidad dinmica del aire kg/m-s
V Volumen m3
Volumen de grano m3
Volumen especfico m3/kg
XI
NDICE DE FIGURAS
Pg.
Figura 1.1 Flujo de energa y masa en el secado ..................................... 23
Figura 1.2 Contenido de agua en el equilibrio de algunos materiales
Adaptado de (3) ....................................................................... 26
Figura 1.3 Secador rotatorio de contraflujo ............................................... 32
Figura 1.4 Secador rotatorio de flujo paralelo ........................................... 33
Figura 1.5 Sistema de secado .................................................................. 33
Figura 1.6 Modelos de paletas elevadoras ............................................... 44
Figura 1.7 Sellos de aire Tomado de (4)................................................ 45
Figura 2.1 Modelo de compartimientos Tomado de (12) ....................... 56
Figura 3.1 Diagrama termodinmico del calentador ................................. 63
Figura 3.2 Proceso de calentamiento del aire ........................................... 63
Figura 3.3 Idealizacin del proceso de secado ......................................... 65
Figura 3.4 Proceso isoentlpico ................................................................ 66
Figura 3.5 Proceso de evaporacin .......................................................... 68
Figura 4.1 Ubicacin del sistema de coordenadas ................................... 72
Figura 4.2 Primer cuadrante de la seccin transversal del tambor ........... 75
Figura 4.3 Determinacin de y .......................................................... 77
Figura 4.4 Recorrido de la paleta .............................................................. 79
Figura 4.5 Etapas del vaciado .................................................................. 81
Figura 4.6 Relacin entre y ................................................. 83
Figura 4.7 Coordenadas del punto ........................................................ 84
Figura 4.8 rea de material en paleta - Caso 1 ........................................ 85
Figura 4.9 rea de material en paleta - Caso 2 ........................................ 86
Figura 4.10 Teorema del valor medio ......................................................... 89
Figura 4.11 Aparato para medicin del ngulo dinmico de reposo
Tomado de (18) ....................................................................... 95
Figura 4.12 Recorrido de la partcula en el secador ................................... 97
Figura 4.13 DCL de la partcula y diagrama de velocidades ..................... 100
Figura 4.14 Coeficientes de resistencia de cuerpos lisos Tomado de (19)
............................................................................................... 102
Figura 4.15 Curva de ajuste para el coeficiente de resistencia del aire .... 104
Figura 4.16 DCL de la partcula y diagrama de velocidades (contraflujo) . 106
Figura 5.1 Transferencia de masa .......................................................... 117
XII
Figura 5.2 Intercambio de calor latente y sensible en una interfaz gas-

lquido Copiado de (21)....................................................... 125
Figura 6.1 Diagrama de flujo para el anlisis de secadores rotativos en
EES........................................................................................ 129
Figura 7.1 Partculas de madera Douglas-fir Tomado de (21) ............. 167
Figura 7.2 Secador de partculas de madera Tomado de (11) ............ 168
Figura 7.3 Posibles configuraciones de paletas para el modelo de secado
A Tomado de (11) ............................................................... 169
Figura 7.4 Comparacin entre los resultados experimentales y simulados
para un secador con paletas de relleno central Tomado de
(11) ........................................................................................ 170
Figura 7.5 Comparacin entre las simulaciones para un secador con y sin
paletas de relleno central Tomado de (11) ......................... 170
Figura 7.6 Representacin esquemtica de la cmara de combustin y del
secador rotativo Tomado de (12) ........................................ 171
Figura 7.7 Resultados tericos y experimentales para el contenido de
humedad en los slidos Tomado de (12) ............................ 173
Figura 7.8 Ajuste del perfil de humedad en los slidos con FCA variable
Tomado de (12) ..................................................................... 174
Figura 7.9 Simulacin del perfil de temperatura de los slidos con FCA
variable Tomado de (12) ..................................................... 175
Figura 7.10 Simulacin del perfil de temperatura del aire con FCA variable
Tomado de (12) ..................................................................... 176
Figura 7.11 Trayectoria de la partcula ..................................................... 182
Figura 7.12 Trayectoria de la partcula (acercamiento)............................. 182
Figura 7.13 Terfiles de humedad (madera) ............................................... 183
Figura 7.14 Perfiles de temperatura (madera) .......................................... 183
Figura 7.15 Trayectoria de la partcula (Regin B) ................................... 185
Figura 7.16 Trayectoria de la partcula (Regin C) ................................... 185
Figura 7.17 Trayectoria de la partcula (Regin D) ................................... 186
Figura 7.18 Perfiles de humedad (zinc) .................................................... 186
Figura 7.19 Perfiles de temperatura (zinc) ................................................ 187
Figura 7.20 Humedad de partcula vs Posicin axial ................................ 189
Figura 7.21 Temperatura del aire vs Posicin axial .................................. 189
Figura 7.22 Temperatura de partcula vs Posicin.................................... 190
Figura 7.23 Humedad de partcula vs Posicin axial ................................ 191
Figura 7.24 Temperatura del aire vs Posicin .......................................... 192
Figura 7.25 Temperatura de partcula vs Posicin.................................... 193
XIII
Figura 7.26 Humedad de partcula vs Posicin (Dimetro de Partcula

Constante) ............................................................................. 194
Figura 7.27 Temperatura del aire vs Posicin (Dimetro de Partcula
Constante) ............................................................................. 195
Figura 7.28 Temperatura de partcula vs Posicin (Dimetro de Partcula
Constante) ............................................................................. 195
Figura 7.29 Efecto de la variacin del tamao de partcula ...................... 205
Figura 7.30 Efecto de la variacin del flujo de producto hmedo ............. 205
Figura 7.31 Efecto de la variacin del flujo msico de aire de secado...... 206
Figura 7.32 Efecto de la variacin de la temperatura inicial del aire ......... 206
Figura 7.33 Efecto de la variacin del contenido de humedad inicial del
producto ................................................................................. 207
XIV
NCIDE DE TABLAS
Pg.
Tabla 1 Datos del secador analizado con el Modelo A ............................... 165
Tabla 2 Datos del secador analizado con el Modelo B ............................... 166
Tabla 3 Resumen de resultados para madera ............................................ 181
Tabla 4 Resumen de resultados para zinc .................................................. 184
XV
INTRODUCCIN
Los secadores continuos de tipo rotativo son ampliamente utilizados para
secar productos granulados en la industria alimenticia, farmacutica, mineral
y de fertilizantes. Su desempeo est gobernado por tres factores:
movimiento, transferencia de masa y transferencia de calor. Estos factores
estn acoplados.
La literatura presenta modelos matemticos basados en la determinacin
experimental de transferencia de masa y coeficientes globales de
transferencia de calor para obtener los perfiles de humedad y temperatura
del producto y del aire de secado. Sin embargo, estos modelos son
aplicables solamente a ciertos tipos de secadores, con configuraciones y
parmetros de operacin particulares.
Un modelo simple, rpidamente ajustable a diferentes condiciones de trabajo,
le permitira al diseador contar con una herramienta para analizar el
comportamiento de secado en funcin de los parmetros de operacin. Con
ello, eventualmente se podra optimizar el diseo para disminuir los costos
operativos y de fabricacin.
Con ese horizonte, el presente trabajo plantea una estrategia de anlisis
basada en principios bsicos de ingeniera, en la que se escoge a una
partcula y se la rastrea durante su paso por el secador, llevando el registro
de la variacin en sus propiedades debido a la transferencia de calor y masa

XVI
ocasionada por la interaccin con el aire de secado. A esta metodologa se la
denomina aqu como lagrangiana.
En mecnica de fluidos existen dos representaciones del flujo, la euleriana
que es la que generalmente se utiliza, y la lagrangiana. La primera dirige su
atencin a un volumen de control fijo en el espacio, mientras que la segunda
persigue a un volumen de material especfico e identificable que es
arrastrado con el flujo. De este modo, la trayectoria de una partcula en un
campo lagrangiano se define como ( ), donde es la posicin inicial
de la partcula y es el tiempo.
En el caso del secador rotativo, la posicin inicial de la partcula est dada
por la altura promedio de cada de los slidos desde las paletas elevadoras.
Luego, estudiando la interaccin de la partcula con la corriente de aire,
ambos de propiedades variables en el tiempo , se dibuja su trayectoria a
travs del secador.
El captulo 1 de este documento es una introduccin a las propiedades de los
slidos granulares, el secado, y el funcionamiento de los secadores rotativos.
En el captulo 2 se presentan algunas de las teoras desarrolladas y avances
con respecto al modelado de secadores rotativos. En el captulo 3 se realiza
un anlisis termodinmico del secador, del cual se obtiene el flujo msico de
aire y sus condiciones en los extremos para condiciones de operacin
definidas. En los captulos 4 y 5 se presentan todas las ecuaciones y

XVII
correlaciones bsicas que permitirn resolver, de manera simultnea, los
problemas de movimiento, transferencia de calor y transferencia de masa.
Posteriormente, el captulo 6 describe paso a paso la estructura del software
desarrollado en EES para el anlisis de secadores rotativos.
Este software es luego ejecutado en dos secadores previamente analizados
por otros investigadores. En cada caso, se comparan los perfiles de
humedad y de temperatura del producto y del aire obtenidos de la simulacin
en EES contra los datos de los modelos planteados por cada investigador.
En algunos casos se cuenta con informacin experimental, en otros no.
Los resultados y las conclusiones con respecto a la metodologa de anlisis
planteada en este trabajo se exponen en el captulo 7 y 8, respectivamente.

18
CAPTULO 1
1. GENERALIDADES
1.1 Caracterizacin de los Slidos Granulares.
Una de las primeras caractersticas a ser evaluadas en la
manipulacin de slidos al granel es la l idez o pe ajosidad, que
se relaciona principalmente con las fuerzas de cohesin entre las
partculas constituyentes. As, cuando estas fuerzas de atraccin
son bajas, el material granular puede fluir libremente bajo la
influencia de la gravedad con partculas movindose de forma
individual una con respecto a otra.
El aumento de las fuerzas de cohesin, que puede ser provocado
por la humedad o por cargas electrostticas, resulta en la tendencia
de generar aglomerados que hacen que el material se mueva
errticamente en forma de r os. Eje plos de solidos al ranel
que normalmente exhiben este comportamiento son harina de trigo,
cacao en polvo y azcar en polvo (1).

19
En secadores rotativos, el parmetro que mide la fluidez de un
material dado es el ngulo dinmico de reposo, que es el ngulo
formado entre la superficie libre de material en una paleta y la
horizontal. Existen materiales para los cuales el ngulo dinmico de
reposo es fuertemente dependiente del contenido de humedad en el
producto, as como tambin hay otros materiales, generalmente los
que presentan buena fluidez en cierto rango de humedad, en los
que el ngulo dinmico de reposo se puede considerar
independiente de la humedad en el producto.
Otras propiedades importantes en el mbito de los materiales de
tipo granular se presentan a continuacin.
Pureza: se refiere a la cantidad de partculas extraas presentes en
un grano individual.
Color: propiedad relacionada con la composicin mineral y, por
ende, es una medida aproximada del ndice de pureza.
Tamao de grano: puede ser determinado mediante varios
mtodos clasificados de acuerdo al mecanismo de operacin, el
cual puede incluir principios mecnicos, pticos y electrnicos,
dinmicos, fsicos y qumicos. El tamizado, el cual est basado en
principios mecnicos, es la tcnica ms simple y la ms utilizada
para determinar distribuciones de tamaos de partcula (2). Los

20
tamices, cuyas medidas estn estandarizadas, son mallas de
alambre tejido que tienen aperturas de forma cuadrada.
Dimetro equivalente: es el dimetro que hace que la partcula no
esfrica sea equivalente a una partcula esfrica; su definicin
depende del criterio de equivalencia empleado. As, se puede tener
dimetro de tamiz, dimetro de Martin, dimetro de Feret, dimetro
del rea proyectada, dimetro de superficie equivalente, dimetro
de volumen equivalente, dimetro de Sauter, o dimetro dinmico.
El dimetro de tamiz, por ejemplo, se define como el ancho mnimo
de la apertura cuadrada a travs de la cual pasar la partcula (2).
Distribucin de tamaos: se puede definir ya sea por el nmero
de partculas o por la masa de las partculas en un determinado
rango de tamaos; ambos tipos de distribuciones son
intercambiables. La primera se puede obtener usando microscopa
o algn otro mtodo ptico; la segunda se puede obtener por
tamizado.
Tamao promedio: se establece en funcin de la distribucin de
tamaos empleada. Para cada distribucin, se utiliza un factor de
ponderacin segn la aplicacin. As, se puede tener un dimetro
medio aritmtico, un dimetro medio superficial, un dimetro
volumtrico medio, un dimetro medio de Sauter, o un dimetro

21
medio de DeBroucker. El dimetro medio aritmtico, por ejemplo, es
el dimetro promedio basado en la funcin de densidad de la
muestra.
Refractariedad: habilidad de un material para soportar
temperaturas elevadas sin fundirse o desmoronarse.
Forma: se refiere a la esfericidad del grano (o angularidad). Toma
el valor de 1 para un grano perfectamente esfrico.
Granulometra: es la determinacin de la distribucin de tamaos
de un conjunto de granos o partculas.
Tenacidad: resistencia del grano a desintegrarse (o partirse)
durante su manipulacin y procesamiento.
Porosidad: valor entre 0 y 1 que define la fraccin del volumen de
grano que es ocupada por poro o espacio vaco.
Gravedad especfica: relacin entre la densidad del grano seco y
la densidad del agua.
ndice de vaco: cantidad relativa de intersticios en muestras al
granel. El porcentaje de espacios es funcin de la granulometra y
forma. Tiene que ver con el empaquetamiento de los granos, mas
no con la porosidad de un grano individual.

22
Peso: es funcin del ndice de vaco y de la gravedad especfica.
Usualmente se expresa en trminos de la densidad al granel o
densidad aparente. La densidad aparente depende del contenido
de humedad en la muestra.
1.2 Fundamentos en el Secado de Slidos.
El secado habitualmente se refiere a la remocin de lquido desde
un slido hmedo mediante evaporacin. Los mtodos mecnicos
utilizados para este efecto, como el exprimido o el centrifugado, no
se consideran operaciones de secado (3).
Denotando al grano hmedo como un sistema termodinmico, ste
estara conformado por el material slido a secar, el agua lquida y
el vapor de agua generado a partir de la evaporacin. En este caso
particular, el entorno ser una corriente de aire a una determinada
velocidad, temperatura y humedad relativa. Este es un sistema
termodinmico abierto porque existe intercambio de masa y energa
entre el sistema y el entorno. Es preciso notar, adems, que en el
interior del sistema habr una transformacin de fase lquido-
gaseosa debido a la evaporacin del agua, as como movimiento de
humedad en el interior del slido. En el diagrama de la figura 1.1 se
indica el flujo de energa y masa del sistema aqu descrito.

23
FIGURA 1.1 FLUJO DE ENERGA Y MASA EN EL SECADO
Visto de este modo, el secado comprendera dos procesos
simultneos:
1. se transfiere calor hacia el sistema para evaporar el lquido;
2. se transfiere masa en forma de lquido o vapor dentro del
sistema, y en forma de vapor desde la frontera del sistema
hacia el entorno.
Los factores que gobiernan la rapidez de estos procesos
determinan la tasa o rapidez de secado. Estos factores pueden ser

24
externos o internos a la partcula. Los factores externos se vinculan
principalmente con las condiciones del aire de secado (presin,
temperatura, humedad relativa) y su velocidad; los factores internos
estn relacionados con el mecanismo de transporte de lquido o
vapor en el interior de la partcula (difusin, capilaridad, gradientes
de presin, etc.).
A medida que el sistema de la figura 1.1 recibe energa del aire
caliente, cierta masa de lquido se evapora en el interior del mismo.
As como la transferencia de calor, la transferencia de masa ocurre
cuando existe una diferencia de concentraciones entre dos regiones
contiguas. Dado que la mayor parte del vapor de agua se encuentra
en una zona cercana a la interfaz gas-lquido, donde inicialmente
ocurre la evaporacin, la concentracin de vapor en el aire en esta
regin es ms alta que la concentracin de vapor en el entorno,
generando la transferencia de masa hacia la regin de menor
concentracin. El vapor que se desplaza hacia fuera del sistema es
entonces arrastrado por la corriente de aire, cuya humedad relativa
aumenta paulatinamente con su avance en el secador.
La diferencia de concentraciones es lo que promueve la
transferencia de masa. Mientras ms amplia es esta diferencia,
mayor es la velocidad de secado. Por eso es importante que el aire
entre al secador con una humedad relativa baja; la humedad

25
relativa es una medida de la concentracin de vapor de agua en el
mismo. El perfil de la capa lmite de concentracin (vase figura 1.1)
es una forma de representar la variacin de la concentracin de
cierta especie en funcin de la posicin con respecto a la interfaz
gas-lquido, tenindose la concentracin ms alta en dicha interfaz y
la ms baja en regiones alejadas de la partcula.
En etapas avanzadas del secado, generalmente cuando el
contenido de humedad en el slido es bajo, puede ser que la
concentracin de vapor en la regin cercana a la interfaz gas-lquido
se acerque a la concentracin de vapor en el entorno, ya sea por la
disminucin de lquido en el sistema o por el aumento de la
humedad relativa del aire. Esto achatara el perfil de
concentraciones, disminuyendo la tasa de transferencia de masa
entre el sistema y el entorno.
Cuando las concentraciones de vapor son iguales, la transferencia
de masa se cancela; la cantidad de agua presente en el sistema en
ese instante, por unidad de masa de slido seco, es lo que se
conoce como contenido de humedad en el equilibrio, Xbs*. Por
ejemplo, si un trozo de madera hmeda es expuesto a una corriente
de aire con una humedad relativa del 60%, punto A en la figura 1.2,
la madera perder agua hasta que el contenido de humedad de la
misma sea de mas o menos el 12%. Una exposicin de la madera a

26
esta misma corriente de aire durante tiempos ms prolongados no
lograr efecto alguno en su contenido de humedad (3).
Aqu se considerar al contenido de humedad en el equilibrio como
un indicador, referencial solamente, de cunta agua se puede
extraer bajo determinadas condiciones de humedad y temperatura
del aire, y condiciones de humedad en los slidos. Es referencial
porque, en ocasiones, es posible extraer humedad ms all del
contenido de humedad en el equilibrio.
FIGURA 1.2 CONTENIDO DE AGUA EN EL EQUILIBRIO DE
ALGUNOS MATERIALES ADAPTADO DE (3)

27
El anterior es un tema importante en el diseo de secadores
convectivos de tipo continuo. En estos secadores, la humedad
relativa del aire aumenta y el contenido de humedad en los slidos
disminuye a medida que se produce el secado. Si se quiere llegar a
niveles de humedad muy bajos en el producto, es necesario que el
contenido de humedad en el equilibrio del material manipulado sea
cercano a cero para cualquier valor de humedad relativa del aire de
secado; esto equivale a tener una curva mas o menos vertical en el
diagrama de la figura 1.2. De otro modo, existir un punto en el
secador donde el aire ya no ser capaz de extraer agua de los
slidos, tornndose en una operacin ineficiente.
1.3 Transporte de Lquido en los Slidos.
El mecanismo de transporte de lquido en el interior de los slidos
depende de su estructura. Estos mecanismos pueden incluir
difusin, flujo capilar, flujo causado por gradientes de presin, flujo
causado por gravedad, y flujo causado por una secuencia de
vaporizacin-condensacin (4).
Es difcil saber con certeza cul de estos mecanismos de transporte
predomina en alguna etapa del secado. ando se abla de etapa
de secado se ace re erencia al contenido de edad en los

28
slidos en cierto momento del secado. Existen experimentos que
han sido diseados para, de alguna manera, tener una idea de qu
mecanismo de transporte se acoplara mejor a la etapa de secado
dada.
Cuando un slido es secado bajo la accin de un entorno de
propiedades constantes, por ejemplo, se suele medir su contenido
de humedad en base seca durante intervalos de tiempo
controlados. Luego se grafica la velocidad de evaporacin contra el
contenido de humedad en el slido.
Bajo estas condiciones experimentales, en el secado de algunos
materiales se pueden identificar dos periodos. En el primero, la
velocidad de evaporacin es constante debido a que la humedad en
los slidos puede migrar fcilmente hacia el exterior. En el segundo,
la rapidez de migracin de lquido decrece provocando que en la
superficie se generen islas secas; esto hace que la velocidad de
evaporacin tambin disminuya.
As, quien ejecuta el experimento podra asociar a la primera etapa
de secado con los mecanismos de flujo capilar y gravedad, y a la
segunda etapa con el mecanismo de difusin, por ejemplo. Esto se
hace en base a criterios previamente establecidos, los cuales estn

29
relacionados con la estructura interna de los materiales. Si no se
hace un anlisis profundo, esto puede llegar a ser muy subjetivo.
En un secador convectivo continuo se tiene un escenario diferente
al del experimento descrito anteriormente. Los slidos en este caso
estn expuestos a corrientes de aire con propiedades
psicromtricas variables, lo que hara an ms complicado
identificar los mecanismos de transporte implicados a diferentes
niveles de humedad.
En este trabajo se supondr que la migracin de lquido desde el
interior de los slidos es suficientemente rpida como para
mantener una superficie expuesta ojada. En otras palabras, se
trabajar bajo la suposicin de que siempre habr lquido
uniformemente distribuido en la superficie de los slidos, y en
contacto directo con la corriente de aire. Esto independiza el
problema de la estructura interna de los slidos y da lugar a que el
secado sea controlado nica y exclusivamente por la tasa de
transferencia de calor convectivo desde el entorno hacia la
superficie evaporativa, y no por el mecanismo de transporte interno
de lquido.
30
1.4 Equipos de Secado.
Los equipos de secado comnmente se clasifican de acuerdo al
mtodo de transferencia de calor empleado para este propsito. Los
ms importantes son los secadores de tipo directo y los secadores
de tipo indirecto. En los secadores directos o convectivos, la
transferencia de calor se logra por contacto directo entre el slido
hmedo y un gas a mayor temperatura, el cual arrastra consigo al
lquido vaporizado. En los secadores de tipo indirecto, el medio de
calentamiento es separado fsicamente de los slidos mediante una
pared o tubo metlico.
Los secadores convectivos se pueden subdividir en secadores
continuos y secadores por lotes (tipo batch). En aquellos del primer
tipo, la operacin es continua y sin interrupcin, siempre y cuando la
materia prima sea suministrada de manera constante. Por otro lado,
en el secado por lotes los equipos son diseados para operar con
cierto tamao definido de lote de producto para ciclos de tiempo
dados.
Este estudio se fundamenta en el anlisis de secadores rotativos,
los cuales se incluyen en el grupo de secadores directos de tipo
continuo. Estos secadores estn constituidos por un cilindro
giratorio que puede tener una ligera inclinacin con respecto a la

31
horizontal. El producto a secar es suministrado en un extremo del
cilindro, movindose a travs del mismo debido a la rotacin, efecto
del arrastre de la corriente de gas, inclinacin del cilindro si
existiese, o a la combinacin de estos factores. Finalmente, el
producto es descargado en el otro extremo para su procesamiento
posterior.
Los gases que circulan a travs del cilindro pueden retardar o
incrementar la velocidad de avance de los slidos, dependiendo si
el gas se mueve en la misma direccin (flujo paralelo) o en direccin
opuesta al avance del producto (contraflujo). La direccin del flujo
de gas relativa a los slidos depende principalmente de las
propiedades del material procesado. El flujo paralelo es utilizado
para materiales sensibles al calor, inclusive para altas temperaturas
de gas a la entrada, debido al enfriamiento rpido del aire durante la
evaporacin inicial de la humedad superficial. Para otro tipo de
materiales, donde el control de temperatura en el producto no es tan
importante, se usa el secado en contraflujo ya que es posible
alcanzar una mayor eficiencia trmica de esta forma (5).
Comercialmente, los secadores rotativos son muy utilizados para el
secado de slidos granulares no muy pegajosos que pueden ser
arrojados sin temor de romperlos (materiales tenaces).

32
1.5 Caractersticas del Secador Rotativo.
Las figuras 1.3 y 1.4 son esquemas de secadores rotatorios de
contraflujo y de flujo paralelo, respectivamente. Aqu se pueden
observar las partes ms importantes que constituyen al secador. La
figura 1.5 es un diagrama general tpico del sistema de secado,
donde se pueden observar equipos auxiliares y accesorios como el
calentador, ventilador, cicln, transportador y chimenea.
FIGURA 1.3 SECADOR ROTATORIO DE CONTRAFLUJO

33
FIGURA 1.4 SECADOR ROTATORIO DE FLUJO PARALELO
FIGURA 1.5 SISTEMA DE SECADO

34
Tambor
El secador rotatorio de tipo directo es un cilindro de metal que
puede estar aislado trmicamente o no, dependiendo de la
temperatura de operacin; en cilindros aislados se puede utilizar
lana de vidrio como aislamiento. La temperatura de operacin, que
es la temperatura del aire a la salida del calentador, est limitada
por varios factores:
1. la resistencia del metal empleado para la fabricacin casco,
que generalmente es de acero al carbono o acero
inoxidable;
2. la resistencia del material granular al desmoronamiento o al
cambio de fase (fusin);
3. aspectos relacionados con el procesamiento del material,
principalmente en productos alimenticios.
En la prctica comercial, se ha encontrado que la longitud de los
secadores ms eficientes vara entre 4 y 10 veces su dimetro, que
normalmente est entre 0.3 y 3 m (4) (5). La pendiente del tambor
con respecto a la horizontal vara de 0 a 8 cm/m, o sea que es de, a
lo mucho, 5.
35
Transmisin
Se suelen emplear dos mtodos de transmisin de potencia para
brindar el carcter giratorio al cilindro del secador. El primero
consiste en fijar un engrane exterior en el tambor, el cual se acopla
a un pin que se encarga de transmitir el movimiento proveniente
de un moto reductor. El segundo mtodo remplaza al engrane por
una catarina sobre la cual se monta una cadena, tambin movida
por un motor.
La potencia del motor estar en funcin del tamao del secador y
de la velocidad de giro en la que se trabaje. Estos equipos operan a
velocidades angulares muy bajas, generalmente entre 3 y 8
revoluciones por minuto.
Adicionalmente, el cilindro tiene instalado, dependiendo de su
longitud, uno o ms anillos metlicos, llamados cintas de apoyo,
que descansan sobre rodillos giratorios. Estos sirven como soportes
simples de rodadura para brindar estabilidad y evitar la eventual
flexin de la estructura del tambor.

36
Calentador
Para calentar el aire se puede utilizar un quemador, un
intercambiador de calor, o alguna otra forma de energa disponible.
En el primer caso el calentamiento es directo. Este mtodo se utiliza
cuando los slidos pueden someterse a elevadas temperaturas sin
temor a ser contaminados por los productos de la combustin. Para
ello, se coloca una cmara de combustin en un extremo del
secador, donde se quema la mezcla aire-combustible y los
productos generados a partir de la combustin son arrastrados a
travs del tambor por medio de un ventilador. Se puede utilizar
diesel o gas natural para este fin.
En el segundo caso el calentamiento es indirecto. Este mtodo se
utiliza cuando existe riesgo de contaminacin de los slidos debido
a los productos de la combustin, o cuando los slidos no pueden
calentarse excesivamente. De un lado del secador se instala un
intercambiador de calor, el cual puede estar conformado por una
carcasa metlica y unos serpentines aleteados. El fluido de
calentamiento puede ser gas caliente de un quemador o vapor de
una caldera.
37
Ventilador
Los gases son forzados a pasar por el tambor mediante un
ventilador. La velocidad del aire en el interior del tambor rara vez es
menor a 0.5 m/s en un diseo econmico (4).
En principio, el ventilador puede ser colocado antes (tiro forzado) o
despus del secador (tiro inducido). Normalmente, se utilizan
ventiladores centrfugos movidos por una banda acoplada a la polea
de un motor elctrico exterior a la carcasa; lo anterior es para evitar
la rpida contaminacin por polvo de las partes mviles del motor.
El uso de un ventilador en tiro forzado (o blower) por s solo es una
prctica inusual debido a que el aire de secado, junto con el
material particulado, tal vez sea soplado hacia los sellos y luego
hacia el exterior, en caso de tener un sellado deficiente. Esto se
puede dar porque la presin en el interior del tambor estara por
encima de la presin atmosfrica (cilindro en presin positiva).
Si se usa un ventilador en tiro inducido, por otro lado, la presin en
el interior del tambor sera ligeramente menor que la presin del
entorno (cilindro en presin negativa), provocando que la eventual
filtracin de aire sea hacia adentro; esto ayuda a controlar la
generacin de polvo.
38
Si la cada de presin en el intercambiador de calor es grande, se
puede colocar un blower aguas arriba para compensar aquello,
adicional al ventilador en tiro inducido ubicado despus del cicln o
colector de polvos. Esta combinacin blower-extractor ayuda a tener
un mejor control de la presin en el secador.
La presin en el interior del tambor debe estar entre 0.25 y 2.5 mm
de agua por debajo de la presin atmosfrica (4). Las variaciones
de presin por cadas de presin en el tambor son despreciables.
En este estudio se supondr que la presin en el secador es
constante e igual a la presin baromtrica.
La presin baromtrica tiene mucha importancia si se desea operar
en ciudades de altura como Quito o La Paz. Las variaciones en el
desempeo del equipo son significativas para estos casos, si se
compara con ciudades del llano como Guayaquil.
Cicln
Cuando el producto a secar contiene partculas muy finas, el aire
probablemente arrastre consigo una cantidad considerable de
material particulado. El material acarreado por la corriente de aire
saliente es comnmente removido en colectores o separadores de
cicln.
39
Un cicln es como un embudo metlico con un sello en la parte
inferior y un tubo de escape en la parte superior. El aire entra por la
parte de arriba tangencialmente a la circunferencia del embudo, y
va formando una espiral a medida que desciende por el mismo.
Este movimiento, o efecto cicln, hace que el material particulado
se dispare hacia las paredes metlicas del separador por accin de
la fuerza centrfuga, y luego se deposite en el fondo por gravedad.
Como el aire parcialmente limpio no tiene otra salida ms que el
tubo en la parte superior, ste asciende por el centro del embudo y
es llevado por el ventilador hacia la chimenea.
El sello en la parte inferior del cicln puede ser un roto-feeder o
alimentador rotatorio. El roto-feeder es un cilindro metlico giratorio
con pedazos de plancha montados en su interior, de forma radial,
que dividen al cilindro en varios compartimientos. En cada
compartimiento se hace un agujero por la parte exterior del cilindro
de tal modo que el dimetro proyectado coincida con el del orificio
en la parte baja del cicln. El roto-feeder se instala con su eje
longitudinal perpendicular al eje longitudinal del cicln. El orifico
inferior del cicln se fabrica de tal modo que la chapa sea siempre
tangencial al movimiento giratorio del roto-feeder; se hace una
boca de pescado. De ese modo, cuando un agujero del roto-
feeder coincide con el agujero del cicln, una cantidad justa de

40
material cae hacia el compartimiento del roto-feeder donde se
encuentra el agujero; cuando no, la superficie cilndrica del roto-
feeder se encarga de sellar el agujero del cicln impidiendo la salida
de aire y producto.
Este accesorio tiene dos propsitos. El primero es impedir que se
escape el aire por la parte inferior del cicln y el segundo es
suministrar (o evacuar del cicln) cantidades precisas de material
en intervalos de tiempo definidos por su velocidad de giro. Los
granos recogidos tal vez sean enviados de vuelta al secador para
un reproceso o simplemente son depositados en un transportador
mecnico, que puede ser un tornillo sinfn, para ser llevados, junto
con el producto seco proveniente del secador, a una siguiente etapa
del proceso. Se puede instalar un tamiz antes del roto-feeder,
adentro del cicln, para separar el producto de granulometra
inadecuada. Este afrecho es luego removido del cicln mediante
una compuerta instalada en el cicln.
Existen tambin otros sistemas donde en lugar de utilizar un cicln,
se instala un lavador de aire. Este equipo es una torre cilndrica alta
que tiene instalado en el interior unas duchas de agua. Cuando el
aire asciende por la torre, el agua proveniente de las duchas cae y
arrastra consigo al material particulado. Abajo del lavador se puede

41
colocar una bomba que lleva el agua con producto a una etapa de
reproceso.
Suministro de material
El mtodo de suministro de material en el secador rotatorio depende
de las caractersticas de los slidos y de la ubicacin y tipo del
equipo de proceso ubicado antes del secador.
Cuando el material viene desde arriba, se utiliza una tolva que se
extiende hacia el interior del cilindro, tal como se muestra en las
figuras 1.3 y 1.4.
Por propsitos de sellado, o si el suministro por gravedad resulta
inconveniente, normalmente se utiliza un tornillo sinfn. Con el
tornillo sinfn se puede controlar mejor el aprovisionamiento de la
materia prima al secador.
En equipos de flujo paralelo, se recomienda el uso de chaquetas de
agua fra para la tolva o transportador si stos estn en contacto
con la corriente de gas caliente. Esto impedira que la temperatura
de la plancha metlica suba demasiado, evitando el
sobrecalentamiento de los slidos sensibles al calor (4).

42
Paletas elevadoras
El secador rotativo est usualmente equipado con paletas montadas
en el interior para el levantamiento y duchado de los slidos hacia la
corriente de gas durante su paso a travs del cilindro. Estas paletas
pueden estar soldadas o empernadas al cilindro. A medida que la
materia prima es depositada en un extremo del tambor, ya sea
mediante una tolva o con un tornillo sinfn, las paletas recogen este
material y lo llevan consigo mediante el giro del cilindro. Existe un
punto en el cual el material empieza a caer desde las paletas,
formando una cascada de producto dentro de la corriente de aire.
Este patrn se repite a lo largo de todo el secador.
Con el uso de paletas se promueve el contacto ntimo entre los
slidos y los gases calientes, siempre y cuando la geometra de la
paleta sea la adecuada.
Las paletas usualmente se ubican de manera seguida, en el sentido
longitudinal del tambor, para asegurar que las cortinas de slidos en
el gas sean lo ms continuas y uniformes como sea posible.
La longitud radial de las paletas suele variar desde un doceavo a un
octavo del dimetro del tambor (4).

43
La forma de las paletas depende de cun pegajosos son los slidos.
Para los materiales que fluyen fcilmente, por ejemplo, se puede
utilizar una paleta de labio recto o a 45. Para materiales pegajosos,
se suele utilizar una paleta radial plana sin labio.
Cuando las caractersticas del material cambian durante el secado,
el diseo de la paleta puede ser modificado de acuerdo a la
evolucin del producto a travs del secador. Algunos secadores de
diseo estndar emplean paletas sin labio en el primer tercio del
secador, paletas con labio de 45 en la parte media y paletas con
labio de 90 en el ltimo tercio del cilindro. Tambin se suelen
colocar paletas espirales justo despus de la descarga de la tolva
de alimentacin, o transportador mecnico, para que el material se
acelere hacia adelante y para prevenir filtraciones de slidos hacia
los sellos de aire (5).
En la figura 1.6 se presentan algunos ejemplos de configuraciones
comunes de paletas comerciales.

44
FIGURA 1.6 MODELOS DE PALETAS ELEVADORAS
Cuando se utilizan secadores de flujo paralelo, es recomendable
dejar el ltimo metro del secador libre de paletas elevadoras para
reducir el material granular que es arrastrado por la corriente de aire
hacia el ducto de escape. Esto ayuda de cierto modo a disminuir la
cantidad de material recogido en el cicln.
Sellos de aire
Existen varios mtodos para sellar el cilindro rotatorio y prevenir la
fuga de aire a travs de la abertura anular entre el tambor y las
partes estacionarias (cabezales de carga y descarga). En la figura
1.7 se muestran tres ejemplos comunes de sellos de aire.

45
En los secadores rotatorios directos no se espera que los sellos
impidan por completo la filtracin de aire, pero con un control
cuidadoso de la presin interna se puede prevenir la produccin de
polvo hacia el exterior, si el cilindro est en presin positiva, o la
filtracin del aire exterior hacia el interior, si el cilindro est en
presin negativa.
FIGURA 1.7 SELLOS DE AIRE TOMADO DE (4)

46
CAPTULO 2
2. MTODOLOGAS DE DISEO DE SECADORES
ROTATIVOS
El problema principal en el anlisis terico y diseo de secadores
rotativos est en desarrollar un modelo que combine el movimiento en
cascada de las partculas con la transferencia de calor y de masa a lo
largo del tambor. Hasta el momento, el diseo de secadores industriales
se ha basado en la experiencia de los fabricantes y proveedores, quienes
han usado informacin obtenida de secadores rotativos en plantas piloto
y de las caractersticas operativas de unidades ya instaladas (5).
En este captulo se citan algunas de las teoras desarrolladas, avances y
resultados obtenidos por diferentes investigadores con respecto a este
tema. Posteriormente, se propone un enfoque alternativo, o
complementario, para el diseo y estudio paramtrico de secadores
rotativos.
47
2.1 Diseo mediante Correlaciones Disponibles.
Un mtodo comn de diseo suele expresar al secado directo como
un mecanismo de transferencia de calor, solamente, de la siguiente
manera:
( )
Donde
= calor total transferido,
= coeficiente global de transferencia de calor volumtrico,
( )
= volumen del secador,
( ) = diferencia media real entre el gas caliente y el material,
Segn Moyers, cuando se desconoce la temperatura de los slidos,
una buena aproximacin para ( ) es la media logartmica entre
las depresiones de bulbo hmedo a la entrada y a la salida del
secador (4). La depresin de bulbo hmedo es la diferencia entre la
temperatura de bulbo seco y la temperatura de bulbo hmedo del
aire en una seccin.
Sin embargo, Friedman y Marshall afirman que el uso de la
diferencia de temperatura media logartmica es incorrecto, debido a

48
que la temperatura de los slidos no vara linealmente con la
temperatura del gas (6).
Para evaluar el coeficiente de transferencia de calor volumtrico, se
han desarrollado correlaciones de la forma , donde
es una constante de proporcionalidad, la velocidad msica del
gas, el dimetro del secador, y una constante.
Por ejemplo, para la fabricacin de secadores comerciales
modernos con una cuenta de paleta por ciclo de 2.4 a 3.0 veces el
dimetro, y operando a velocidades perifricas de tambor entre 0.3
y 0.4 m/s, se recomienda el uso de la ecuacin 2.1 (4) (5).
( ) ( ) ( )
(2.1)
Donde
= calor total transferido,
= longitud del secador,
= dimetro,
= velocidad msica del aire, ( )
( ) = temperatura media logartmica entre la depresin de bulbo
hmedo del gas a la entrada y la depresin de bulbo hmedo del
gas a la salida del tambor,

49
Si se efecta un anlisis termodinmico al secador (referirse al
captulo 3 de este documento), es posible determinar el calor total
transferido durante la evaporacin. Con esta informacin, y
haciendo uso de la ecuacin 2.1, se puede determinar la longitud
del secador para un correspondiente dimetro de tambor, o
viceversa. No obstante, debe notarse que los resultados que se
obtengan son aplicables solamente para secadores con una
configuracin similar a la que se us para derivar la ecuacin
anterior.
En ningn momento se especifica, por ejemplo, el tipo o las
caractersticas del material a secar, las dimensiones de las paletas,
inclinacin del tambor, su velocidad de giro, o la direccin de flujo
del aire con respecto al movimiento del producto. Tampoco se toma
en cuenta el fenmeno de transferencia de masa.
El punto es que, si se desea obtener un resultado medianamente
confiable con este mtodo para alguna aplicacin en particular, es
necesario contar con varias correlaciones de para
configuraciones de secadores rotativos tipo. As, quien construye
estos equipos podra asociar la aplicacin de secado requerida a
alguno de estos modelos y seguir adelante en el proceso de diseo.

50
2.2 Diseo Basado en Pruebas de Planta.
Los modelos matemticos existentes solo son aplicables a diseos
especficos de secadores rotativos con modos de operacin y
configuraciones particulares. Es ms, en la literatura se encuentra
muy poca informacin con respecto al coeficiente global de
transferencia de calor volumtrico, , aplicable al diseo de
secadores rotativos industriales.
Otro mtodo consiste en dimensionar los secadores rotativos
mediante procedimientos de prueba y error en plantas piloto. Como
una primera aproximacin, se puede utilizar el procedimiento
planteado en la seccin anterior, haciendo uso de la informacin
que se encuentre disponible.
Desde este punto de vista, la experiencia del fabricante desempea
un rol importarte porque se necesita seleccionar un dimetro de
secador adecuado, escoger la geometra de paleta, definir el modo
de secado (contraflujo o flujo paralelo) y jugar con las condiciones
del aire caliente (temperatura, humedad, y velocidad), para
optimizar la operacin y economa del secado. La habilidad de
quien disea tendr afectacin directa sobre el consumo de tiempo
y recursos durante la fabricacin.

51
En lugar de seguir un proceso de optimizacin en el diseo, aun
cuando los costos capitales y operativos son mayores, es una
prctica comn construir secadores sobredimensionados e
ineficientes, pero mecnicamente funcionales y bien probados en
operacin (5).
2.3 Avances en el Estudio de Secadores Rotativos.
En 1993 se propone un modelo matemtico que predice el tiempo
de residencia del producto en un secador rotatorio. El modelo se
basa en las ecuaciones diferenciales que describen el movimiento
de las partculas (7). Esta propuesta considera la interaccin entre
los granos individuales y la corriente de aire, as como distintas
configuraciones de paletas e inclinaciones de tambor. Se ignoran
los cambios en las propiedades tanto del aire como de los slidos a
travs de su paso por el secador.
En el ao 2007 se presenta un estudio sobre el desempeo de un
secador rotativo en funcin del nmero de paletas elevadoras (8).
Una parte importante de este trabajo es que se logra determinar la
cantidad de material presente en la paleta y la altura de cada para
diferentes posiciones angulares de tambor y geometras de paleta.
Esto fue posible gracias a una relacin derivada por Schofield y

52
Glikin en 1962, la cual evala el ngulo de reposo de un polvo en
funcin de su coeficiente de friccin dinmico (9).
En el 2009 se presenta una metodologa para el modelado y control
de secadores rotativos directos continuos en flujo paralelo (10). La
estrategia consiste en dividir el secador en un nmero finito de
elementos volumtricos dispuestos en serie y aplicar a cada uno de
ellos las ecuaciones diferenciales de conservacin de masa y
energa. El modelo hace uso del concepto de coeficiente
volumtrico de transferencia de calor, el cual se determin a partir
de pruebas efectuadas en una planta piloto de secado de arena fina
ubicada en el Departamento de Ingeniera de Sistemas y
Automtica de la Universidad de Sevilla.
A continuacin se exponen dos trabajos que se tomarn como
referencia para establecer comparaciones con el modelo a
plantearse en este estudio.
Secado de partculas de madera
En su tesis doctoral presentada en 1983 en la Universidad Estatal
de Oregon Engineering Analysis of a Rotary Dryer: Drying of Wood
Particles rederic A. Ka e, desarrolla un modelo de simulacin
computacional para un secador rotativo de partculas de madera,
con y sin paletas de relleno central (11).

53
Para la ejecucin del programa, se necesita conocer el dimetro y la
longitud del secador, su velocidad de rotacin, la geometra de las
paletas, y las propiedades del material. Con ello, se efecta un
anlisis del tiempo de residencia de las partculas en el tambor, y
luego se incorporan los procesos de transferencia de calor y de
masa.
La estrategia consiste primero estimar la altura y el tiempo de cada
promedio de una partcula para poder determinar el avance
horizontal de la misma. Con esta informacin, y considerando el
tiempo que la partcula permanece en la paleta, se estima el
nmero de ciclos de secado y se divide al secador en segmentos de
longitud igual al desplazamiento longitudinal de la partcula que
desciende desde la altura promedio de cada. Se considera a cada
segmento como un volumen de control, en donde se resuelven los
problemas de transferencia de calor y de masa, asumiendo que las
propiedades de las partculas y del gas son constantes, y
considerando que el tiempo de vuelo y los desplazamientos
horizontales de la partcula son siempre los mismos.
Los problemas de movimiento, transferencia de calor y transferencia
de masa se resuelven de forma separada.

54
No se involucran constantes empricas o ecuaciones de
transferencia de calor desarrolladas para un sistema de secado
particular. El estudio est basado en principios elementales y
relaciones empricas que normalmente se pueden encontrar en los
libros de transferencia de calor. No se utilizaron correlaciones
obtenidas a partir de estudios previos en secadores rotativos.
Se considera, entre otras cosas, el arrastre del aire, la prdida de
calor desde el cilindro, la transferencia de calor en el interior de los
slidos y varios modelos de secado de la madera.
En el trabajo de Kamke se ejecutaron experimentos en un secador
rotativo industrial para verificar el desempeo del modelo de
simulacin como un pronosticador del comportamiento de secado.
Las variables de control fueron la rapidez de giro del tambor, flujo
msico de aire, y temperatura de entrada del gas. Las mediciones
de temperatura del gas, temperatura de las partculas, y contenido
de humedad de las partculas, se obtuvieron a lo largo de la longitud
del tambor.
La comparacin entre los resultados predichos por el modelo de
Kamke y los resultados experimentales fue buena, indicando un
porcentaje de error cuadrtico medio de 22.2 en la prediccin de
contenido de humedad de las partculas al final del secado.

55
No se ejecutaron ensayos experimentales en un secador sin paletas
de relleno central, pero s hay datos de simulacin computacional.
Secado de concentrados de zinc
En el 2011, Olukayode Ajayi, de la Universidad James Cook en
Australia, presenta una tesis doctoral que tiene por tema Multiscale
modelling of industrial flighted rotary dryers, donde se desarrolla un
modelo que simula los mecanismos de transporte, transferencia de
masa y transferencia de calor para un secador rotativo industrial de
concentrados de zinc (12).
El secador utilizado para el estudio est constituido por secciones
con paletas (B, C y D) y secciones sin paletas (A y E). El anlisis del
transporte de slidos en la zona con paletas se basa en una tcnica
computacional de modelado de compartimientos que divide al flujo
de slidos en una fase pasiva y una activa, con varios
co parti ientos o tanq es bien ezclados dispuestos en paralelo
y en serie. La fase activa es la que est en contacto con la corriente
de aire, mientras que la fase pasiva es la que se desplaza en la
parte baja del tambor debido a su inclinacin y rotacin (kilning).
Existe intercambio de masa entre fases.

56
La proporcin de partculas en contacto con el aire se determina
utilizando una tcnica de simulacin en CFD (Computational Fluid
Dynamics) basada en un enfoque euleriano.
Entre compartimiento y compartimiento, se resuelven las
ecuaciones diferenciales de balance de masa y energa. El calor es
transferido por conveccin y radiacin en la fase activa. Se
considera la prdida de calor desde el tambor hacia el entorno. La
figura 2.1 es una representacin del modelo planteado.
FIGURA 2.1 MODELO DE COMPARTIMIENTOS

TOMADO DE (12)
57
Se llevaron a cabo varias simulaciones utilizando el programa
gPROMS (software de modelado de procesos). El modelo de
transporte de slidos fue validado mediante comparaciones contra
distribuciones de tiempo de residencia bajo diferentes condiciones
experimentales. La validacin de los balances de energa y masa se
bas en Informacin de Proceso (PI), distribucin experimental de
tiempo de residencia y datos del perfil de humedad de los slidos.
Se manipularon parmetros del modelo relacionados con el rea
superficial en contacto con el aire de secado para ajustar el perfil de
humedad terico con el experimental. Tambin se estimaron los
perfiles de temperatura del aire y de los slidos, pero no existe
informacin experimental que los corrobore.
2.4 Metodologa Propuesta.
En el presente trabajo se renen algunas de las ideas presentadas
en publicaciones anteriores y se genera una estrategia de anlisis y
dimensionamiento de secadores rotativos para materiales de tipo
granular.
Primero, se hace un anlisis termodinmico del secador para
obtener el flujo msico de aire necesario para secar cierta cantidad

58
de producto hmedo bajo condiciones dadas, en un tiempo
determinado.
Partiendo de un dimetro de secador, se estudia el movimiento de
los granos en la seccin transversal y se calcula un ngulo de cada
representativo de toda la cortina de slidos. Con dicho ngulo, se
procede a determinar una altura de cada promedio. Esta altura de
cada, junto con la velocidad de la corriente de aire (obtenida a
partir del flujo msico y de la seccin transversal del tambor), sirven
como punto de partida para simular el movimiento de una partcula
individual a travs del secador, utilizando los modelos matemticos
pertinentes.
Simultneamente, es decir, a medida que la partcula avanza en la
direccin axial, se calcula la tasa de transferencia de calor y de
masa en intervalos de tiempo o pasos previamente definidos. Con
ello, y bajo la suposicin de que la partcula es representativa de
todas las dems, se extienden los resultados hacia todo el flujo de
partculas y se corrigen las propiedades tanto del aire como de los
slidos en cada paso consecutivo. Esto se hace simplemente
multiplicando la cantidad de agua evaporada desde la partcula, en
un intervalo de tiempo dado, por el flujo total de partculas; las
propiedades del aire se actualizan en base este ltimo resultado.

59
En el proceso, se van sumando los desplazamientos sucesivos de
la partcula hasta que sta adopte el nivel de humedad final
requerido, obtenindose as la longitud de tambor correspondiente.
Con este enfoque lagrangiano, donde se rastrea a una partcula a
travs del secador, se espera obtener la evolucin histrica de
humedades, temperaturas, flujos de masa y energticos de ambos,
los slidos y el aire de secado. En los captulos que siguen se
explica ms detenidamente el procedimiento seguido y la
herramienta computacional empleada para el desarrollo de esta
propuesta terica.
60
CAPTULO 3
3. ANLISIS TERMODINMICO
En el presente captulo se desarrollarn balances de masa y energa en
el calentador y en el secador con el objetivo de encontrar las condiciones
del aire y su flujo msico, requeridos para determinados parmetros de
operacin. No se incluyen anlisis de exerga.
3.1 Calentador.
Previo al ingreso al secador, el aire debe ser calentado para
incrementar su potencial de secado. Ya sea que el aire se caliente
directamente mediante un quemador, o indirectamente en un
intercambiador de calor, lo que se produce es un aumento de
temperatura que conlleva a una disminucin sustancial de su
humedad relativa.
61
Hay que saber diferenciar entre la humedad absoluta del aire, ,y
su humedad relativa, . La humedad absoluta constituye la cantidad
real de agua presente en el aire, mientras que la humedad relativa
es un indicador de cunta agua adicional se podra introducir en el
aire antes de que ste alcance el estado de saturacin. En la
saturacin, el aire no puede recibir ms agua y su potencial de
secado disminuye a cero.
La humedad absoluta se suele expresar en libras de agua por libra
de aire seco, mientras que la humedad relativa es una cantidad
porcentual. Por ejemplo, una humedad relativa del 40% a 25C
indica que, a esa temperatura, el aire an puede recibir un 60% de
agua antes de llegar a la saturacin.
El manejo de humedades y temperaturas para mezclas de gases se
conoce como psicrometra. En este documento, cuando se hable de
mezcla de gases siempre se har referencia al aire y al vapor de
agua. El tratamiento ser el mismo ya sea que el calentamiento se
de en un intercambiador de calor o en una cmara de combustin.
Si bien es cierto que en el ltimo caso el aire se mezcla con los
productos de la combustin del quemador, la mezcla sigue siendo
sustancialmente aire y agua, y se lo idealizar como tal. Adicional a
esto rec ente ente se tilizar el tr ino aire de manera
indistinta para hacer referencia a la mezcla aire/vapor de agua.

62
Cuando se desee hacer referencia solamente al flujo msico de aire
co o tal sin a a se tilizar el tr ino aire seco si bolizado
por . El vapor de agua arrastrado por el aire ser representado
por .
La figura 3.1 es un diagrama termodinmico del calentador, el cual
recibe aire a presin atmosfrica, con una temperatura de bulbo
seco y humedad relativa dada por las condiciones ambientales.
Estas tres propiedades (presin, temperatura y humedad) definen
por completo el estado termodinmico de la mezcla a la entrada del
calentador. La humedad absoluta puede ser determinada en base a
las tres propiedades anteriores.
El calentador es un medio cerrado que no permite el ingreso de
humedad exterior. Tampoco existe condensacin en su interior
porque el aire est siendo calentado y no enfriado. Por lo tanto, el
contenido real de humedad o humedad absoluta permanece
constante en este volumen de control. El resultado del
calentamiento es un aire con mayor temperatura y menor humedad
relativa, es decir, con mayor capacidad de secado. La figura 3.2
muestra como lucira este proceso en un diagrama pricromtrico.

63
FIGURA 3.1 DIAGRAMA TERMODINMICO DEL CALENTADOR
FIGURA 3.2 PROCESO DE CALENTAMIENTO DEL AIRE
El calor agregado generalmente proviene de un combustible fsil,
que puede ser diesel o gas natural, o del vapor generado en una
caldera. Considerando los supuestos listados en la parte inferior, la
energa requerida para calentar el aire de una temperatura a una
temperatura est dada por la ecuacin 3.1. Si se divide este calor
para la eficiencia del calentador, se obtiene la energa total que

64
debe ser entregada por la combustin en el caso de que se utilice
un quemador. Dependiendo del poder calorfico del combustible, y
de su dosificacin, se puede determinar el consumo. El mismo
principio se aplica cuando se utiliza vapor como fluido de
calentamiento.
Supuestos:
Estado estable
( ) (3.1)

En la ecuacin anterior, es el calor especfico promedio de la
mezcla aire/agua entre las temperaturas de entrada y salida en el
calentador. Este valor puede aproximarse mediante el calor
especfico del aire solamente sin provocar errores apreciables en
los clculos.
3.2 Secador.
Una vez en el tambor giratorio, la mezcla atraviesa un proceso
diferente al que ocurre en el calentador. Durante su paso por el
secador, el aire cede energa en forma de calor hacia los slidos,
que estn a menor temperatura, para producir la evaporacin de la
humedad superficial. La masa de agua evaporada sale de las

65
partculas debido a una diferencia de presiones parciales de vapor
(ver seccin 1.2) y es arrastrada por la corriente de aire,
incrementando su humedad absoluta y relativa. Aqu el intercambio
neto de energa es cero, ya que la misma energa que cede el aire
para la evaporacin es devuelta hacia l en forma de vapor de
agua. Como resultado, se obtiene un proceso donde la entalpa de
entrada del aire es igual a su entalpa de salida. Las prdidas de
calor hacia el entorno se suponen despreciables a pesar de que el
casco del secador no necesariamente est aislado. Este proceso se
representa en las figuras 3.3 y 3.4. Por fines didcticos, se ha
simbolizado a la humedad de los slidos como una masa de agua
depositada en el secador.
FIGURA 3.3 IDEALIZACIN DEL PROCESO DE SECADO

66
FIGURA 3.4 PROCESO ISOENTLPICO
El prrafo anterior es la descripcin de lo que ocurrira en una
situacin ideal; sin embargo, existen algunos detalles que merecen
ser citados.
Al entrar en contacto con el aire, los slidos hmedos a menor
temperatura adquieren calor sensible hasta que su superficie
alcanza la temperatura de evaporacin bajo las condiciones del
medio en ese momento. Esta etapa de ajuste de temperaturas
ocurre solamente en un pequeo tramo del secador (a la entrada de
los slidos), donde el calor sensible cedido por el aire es
relativamente pequeo en comparacin al calor total que cede para
la evaporacin del agua (5). Por ese motivo, este calor no ser
tomado en cuenta para los clculos, sino que se asumir que todo
el proceso en el interior del secador es netamente evaporativo. Para

67
ello, se supondr que los slidos siempre estarn abrigados por una
pelcula de agua.
3.3 Determinacin del Flujo Msico del Aire de Secado.
A parte de modelar el fenmeno psicromtrico, el objetivo principal
de realizar un anlisis termodinmico global del secador es
encontrar el flujo msico de aire necesario, y sus propiedades, para
las condiciones de secado dadas. Adicionalmente, se pueden
derivar otros parmetros como la cantidad de agua a evaporar, el
calor total de vaporizacin, la energa requerida para calentar el aire
y el consumo de combustible en el calentador si es el caso.
Para ello, se debe partir de los siguientes datos, necesarios para el
diseo o dimensionamiento de secadores rotativos:
condiciones ambientales: presin, temperatura y humedad
relativa del aire;
flujo msico de producto a la entrada o a la salida del secador;
humedad inicial y final del producto;
temperatura del aire despus del calentador;
humedad relativa del aire a la salida del secador

68
Los dos ltimos son parmetros que estn ligados directamente a la
economa del equipo. A mayor flujo y temperatura del aire, mayor
ser el costo operativo por tonelada de producto seco. Aqu entran
en juego factores como el tipo de calentamiento, la eficiencia del
calentador, y la clase de combustible u alguna otra fuente de
energa empleada.
En la figura 3.5 se puede observar que, durante el secado, la masa
slida del producto permanece constante. Lo que vara es la masa
del sistema agua-slido debido a la evaporacin (vase tambin la
figura 1.1). Lo anterior se cumple siempre y cuando el slido posea
una tenacidad y refractariedad suficiente como para no
desmoronarse debido al impacto o a la temperatura.
FIGURA 3.5 PROCESO DE EVAPORACIN

69
Los valores y representan el porcentaje de agua en base
hmeda antes y despus del secado. Normalmente, son datos
conocidos. Se dice que es de base hmeda porque se mide con
respecto a la masa total del sistema, slido ms agua.
Si es la cantidad de producto hmedo suministrada por hora,
el flujo msico de agua y la masa slida, se puede afirmar lo
siguiente:
( )
Entonces, con el dato inicial de producto hmedo a secar y su
contenido de humedad respectivo, se puede calcular fcilmente el
flujo msico de slidos secos:
( )( ) (3.2)
La cantidad de agua a evaporar, , se puede determinar como
sigue:
( ) ( ) ( ) ( )
70

(3.3)
El agua evaporada pasa al aire de secado aumentando su
humedad absoluta. El calor total de vaporizacin se puede
determinar multiplicando la cantidad de agua evaporada por el calor
latente de vaporizacin, el cual se evala a la temperatura
superficial de los slidos, tal como se expresa en la ecuacin 3.4.
Este es el calor al cual se hace referencia en la seccin 2.1.
(3.4)
Partiendo de la definicin de humedad absoluta, y considerando
que ( ) representa el contenido inicial de vapor de agua en el
aire, se obtiene que
( ) ( )

( ) ( )

71

(3.5)
A manera de resumen, se puede seguir el siguiente procedimiento
para calcular el flujo de aire y el calor suministrado en el quemador
o intercambiador de calor:
1. Datos: , ,( )
2. Calcular mediante la ecuacin 3.2
4. Datos: , , (propiedades del aire atmosfrico)
5. Calcular
6.
7. Dato: Temperatura a la salida del horno o intercambiador
de calor (la escoge el diseador)
8. Con y , determinar
9.
10. Dato: humedad relativa a la salida del secador (la escoge
el diseador)
11. Con y , calcular

72
CAPTULO 4
4. MOVIMIENTO DE LOS SLIDOS EN EL SECADOR
Para abordar este captulo, se considerar al secador como un cilindro
circular horizontal, con o sin inclinacin, y se emplear un sistema de
coordenadas alineado con sus ejes principales, tal como se muestra en
la figura 4.1.
FIGURA 4.1 UBICACIN DEL SISTEMA DE COORDENADAS

73
Primero se estudiar el movimiento en el plano XY, donde los granos son
recogidos por las paletas y luego arrojados para formar una cascada de
slidos. Para ello, se utilizar una metodologa geomtrica que fue
inicialmente desarrollada por Glikin (13), y que ha sido publicada en
obras de autores ms recientes (8) (11) (14) (15).
Luego se determinar un ngulo de cada representativo de toda la
cortina de slidos, en base al teorema del valor medio para integrales, y
se calcular la altura promedio de cada correspondiente a dicho ngulo.
La contribucin de este trabajo empieza con el anlisis del avance de
una partcula individual en el secador, la misma que es recogida por una
paleta en la parte inferior del tambor y arrojada hacia la corriente de aire
desde la altura calculada previamente. Este movimiento toma lugar en el
plano YZ y se repite a lo largo de todo el secador. Aqu es importante el
efecto que pueda tener la inclinacin del tambor sobre el movimiento de
la partcula, ya sea para secadores en flujo paralelo o en contraflujo.
En el anlisis se considerar, entre otras cosas, la influencia de la
gravedad, el arrastre del aire, y las variaciones de las propiedades tanto
de los slidos como las del medio de secado a medida que la partcula
avanza en el secador.
74
4.1 Cada de las Partculas desde las Paletas.
Si se tiene material granular en un recipiente y este ltimo es
volteado lentamente para depositarlo sobre una superficie plana, se
puede observar la formacin un pilo cnico a medida que los granos
caen paulatinamente. El ngulo interno formado entre la superficie
del pilo cnico y la superficie horizontal se conoce como ngulo de
reposo (16).
En un tambor rotatorio ocurre un efecto similar. Desde la parte baja
del tambor, las paletas elevadoras recogen cierta cantidad de
material que en un principio es mxima, y luego va disminuyendo a
medida que el material cae desde las paletas por gravedad. La capa
superior de material forma con la horizontal un ngulo que depende
del volumen de material presente en la paleta. Este ngulo es
anlogo al ngulo de reposo mencionado anteriormente. Sin
embargo, se lo denomina ngulo dinmico de reposo debido a
que el marco de referencia (tambor) est en movimiento (9).

75
FIGURA 4.2 PRIMER CUADRANTE DE LA SECCIN
TRANSVERSAL DEL TAMBOR
La figura 4.2 es un diagrama donde se definen algunas de las
variables a utilizarse en el mtodo geomtrico a exponerse a
continuacin. En este caso, se trabajar con una paleta de dos
segmentos paralela al eje longitudinal del tambor. Sin embargo,
existen investigadores que han desarrollado tratamientos
matemticos para otras geometras de paletas (17).
En la figura 4.2 se representa una paleta con labio de dimensiones
arbitrarias, dadas por la longitud del labio, , y la longitud de la
paleta, . El ngulo formado entre estos dos segmentos es .
Ntese que el segmento no necesariamente est radialmente

76
alineado con el centro del tambor. Se supondr que las paletas son
paralelas al eje longitudinal del tambor.
El segmento representa la capa superior de material siendo
transportado por la paleta en un instante dado. Los slidos estn
representados por la regin punteada. El ngulo formado entre esta
lnea y la horizontal (eje X) es el ngulo dinmico de reposo, .
La posicin angular de la paleta est dada por , que es el ngulo
formado entre la lnea que va desde el centro del tambor hasta el
filo del labio de la paleta, , y el eje horizontal. Al rotar, el punto
dibuja una circunferencia de radio . El radio del tambor es .
Otro detalle a considerar es la utilizacin de dos sistemas de
coordenadas, X-Y y x-y. El primero es un sistema de coordenadas
global no rotatorio (SCG) con origen en el centro del tambor,
mientras que el segundo es un sistema de coordenadas local
rotatorio (SCLR) con origen en el filo del labio. El eje x de este
ltimo sistema est alineado con el labio. El ngulo representa el
desfase instantneo entre los sistemas de coordenadas global y
local, medido con respecto a la horizontal; en la figura 4.2 es
positivo. De ahora en adelante, las letras en maysculas denotarn
puntos del plano referenciados al SCG, mientras que las letras en
minsculas denotarn puntos referenciados al SCLR.

77
La figura 4.3 es un esquema que permitir deducir las cantidades
y . El ngulo es un valor dado por el diseador y define el
grado de desalineamiento angular (offset) de la paleta con respecto
a la lnea radial del tambor.
FIGURA 4.3 DETERMINACIN DE Y
De la figura anterior, se obtiene lo siguiente:
Donde , , son las longitudes de labio y paleta
respectivamente.
78
Entonces,
( )

( )
Por otro lado,

( )
( )
( )
Como ya se mencion, los ngulos y deben ser
proporcionados por el diseador. Vale la pena indicar que puede
ser positivo o negativo; en la figura 4.3 ste se asumi como
positivo.
La figura 4.4 es un intento de simular lo que ocurre durante el
recorrido ascendente de la paleta. Es preciso indicar que siempre
es mayor que cero, desde que empieza el recorrido de vaciado

79
( ) hasta que se ha descargado por completo el material
( ).
FIGURA 4.4 RECORRIDO DE LA PALETA
En 1962, Schofield y Glikin desarrollaron una relacin que permite
estimar la variacin del ngulo dinmico de reposo en funcin de la
posicin angular de la paleta (9). Dicha relacin se presenta a
continuacin.
( )
(4.1)
( )
En la ecuacin anterior, es la velocidad angular del tambor y es
el coeficiente de friccin dinmico entre capas adyacentes de

80
material granular. Este coeficiente es una medida de la resistencia
que presentan los granos al deslizamiento durante el movimiento
giratorio del secador. La relacin de Schofield y Glikin es vlida
solamente para ciertos materiales, generalmente los que presentan
buena fluidez. El ngulo dinmico de reposo es la nica propiedad
del producto involucrada en este anlisis.
Volviendo a la figura 4.4, se puede apreciar que en la primera parte
del recorrido, el extremo derecho de la capa superior de material
(superficie libre) est en contacto con el cilindro del secador.
Mientras se produce el vaciado, existe un momento en que esta
capa de slidos deja de hacer contacto con la pared del tambor, y
solamente toca a la paleta.
En este punto, se deben diferenciar dos etapas del vaciado:
1. Material en contacto con paleta y tambor
2. Material en contacto solo con paleta
Para distinguir el primer caso del segundo, se utilizan los subndices
y , respectivamente, en las variables que se presentan a
continuacin.
81
ETAPA 1 ETAPA 2
FIGURA 4.5 ETAPAS DEL VACIADO
El objetivo consiste en determinar las coordenadas locales de los
puntos , , , y , para luego realizar el clculo del rea
delimitada por estos puntos, segn sea el caso.
A y B son puntos fijos en el sistema de coordenadas local, que
estn dados por las coordenadas ( )y( ). As,
( )
( )
En el sistema de coordenadas global, el punto est dado por
;
82
La recta que contiene a los puntos y , referida al SCG, tiene por
ecuacin
( )
Por otro lado, la circunferencia del tambor est definida por
La interseccin de estos dos lugares geomtricos, la recta y la
circunferencia, da como resultado las coordenadas del punto
referido al SCG. El siguiente paso consiste en representar este
punto en el sistema de coordenadas local.
Del esquema presentado en la figura 4.6, se puede deducir lo
siguiente:
( )
Dado que , se obtiene

83
Tambin se puede verificar que
FIGURA 4.6 RELACIN ENTRE ( )y( )
Resolviendo estas dos ltimas ecuaciones de manera simultnea,
se pueden obtener las coordenadas del punto referidas al SCLR,
o sea, ( ). Con esto se soluciona parcialmente el primer caso;
est pendiente determinar el rea.
El segundo caso se ilustra en la figura 4.7, de la cual se deducen
las ecuaciones de las rectas cuya interseccin genera el punto .

84
FIGURA 4.7 COORDENADAS DEL PUNTO
Recta 1:
( )
Recta 2:
( )( )
( )
La interseccin entre las rectas 1 y 2 arroja las coordenadas
( ).
Hasta el momento, se ha expuesto un mtodo para determinar las
coordenadas de los puntos , , , y . Ahora se procede a

85
calcular el rea correspondiente a la seccin transversal de la cama
de material en la paleta.
Caso 1:
FIGURA 4.8 REA DE MATERIAL EN PALETA - CASO 1
( )
Aqu la coordenada representa la posicin vertical del extremo
derecho de la capa superior de material.
Para se hace una aproximacin mediante la cuerda .
( )
86
El rea total es .
Caso 2:
FIGURA 4.9 REA DE MATERIAL EN PALETA - CASO 2
( )
( ) ( ) ( )

87
Cuando , que es equivalente a , la paleta est
completamente vaca y, por lo tanto, . Como siempre es
mayor o igual que cero, entonces debe existir un instante del
recorrido en el que el ngulo cambia de signo positivo a negativo.
Esto ocurre cuando el labio est completamente horizontal, alineado
con el eje X del SCG.
Si es la profundidad de la paleta, medida a lo largo del eje Z, y
asumiendo una seccin transversal constante, el volumen de
material en la paleta estara dado por ( ) , donde depende
de la posicin de la paleta, .
La masa de material est relacionada con el volumen mediante la
densidad al granel o densidad aparente, .
( ) ( ) (4.2)
Ahora ya es posible conocer la cantidad de material granular
presente en una paleta en funcin de la posicin angular de la
misma, o en funcin del tiempo, lo cual representa un gran paso en
el desarrollo de este estudio.

88
4.2 ngulo Representativo y Altura Promedio de Cada.
Las partculas presentes en la cascada de slidos experimentan
trayectorias diferentes, dependiendo de la posicin inicial en la que
dejaron la paleta elevadora. Por lo tanto, es imposible que los
slidos compartan exactamente las mismas propiedades al final de
la operacin de secado. Entonces, cuando se hable de contenido
de humedad, tamao de grano, posicin, velocidad, o algn otro
parmetro que haga referencia a los slidos, el lector deber pensar
en seguida en promedios o distribuciones.
Habiendo dicho esto, se propone escoger una de las tantas
partculas presentes en el secador, seguirla a lo largo de su
trayectoria, y llevar el registro de la variacin de sus propiedades.
Esta tcnica lagrangiana permitir, mediante el uso de
extrapolaciones, tener una visin general sobre la variacin de
humedad en el producto durante su paso por el secador.
La primera etapa de este mtodo consiste en hallar un ngulo que
sea representativo de toda la cada de material durante un ciclo o
revolucin de paleta. Para ello, se utilizar el teorema del valor
medio para integrales.
Este teorema establece que, si una funcin es continua en un
intervalo cerrado, entonces existe un punto en el interior del

89
intervalo tal que la integral (rea bajo la curva) es igual a la
diferencia entre los valores extremos del intervalo multiplicada por la
funcin evaluada en el valor interior.
En este caso, la funcin ser la cantidad de masa en la paleta, ,y
la variable independiente ser la posicin angular, .
Por el momento, suponga que la masa en la paleta est relacionada
con la posicin angular del tambor de una manera similar a la que
se muestra en la figura 4.10.
FIGURA 4.10 TEOREMA DEL VALOR MEDIO

90
Aplicando el teorema del valor medio, se deduce lo siguiente:
( )

(4.2)
La ecuacin anterior puede ser expresada de forma numrica como
sigue:

( ) (4.3)
Donde es la fraccin de producto que ha cado desde la paleta en
el intervalo .
La altura de cada correspondiente a este ngulo medio estara
dada por
( ) (4.4)
Esta altura promedio de cada es importante porque representa el
punto de partida para el estudio del movimiento en el plano YZ. El
tiempo de vuelo de la partcula, tomando como referencia esta
altura, es el tiempo de secado por revolucin de tambor.

91
Dado que las condiciones en el secador no son constantes de una
seccin a otra, es inevitable pensar que, tanto la altura promedio de
cada, como el tiempo de vuelo, son parmetros variables en la
direccin Z. Si se observa con cuidado la deduccin presentada en
la seccin 4.1, el lector podr notar que la nica propiedad de los
slidos que tiene influencia directa sobre la altura promedio de
cada es el ngulo dinmico de reposo.
En 1969, Kelly investig experimentalmente los factores que
contribuyen a la variacin del ngulo dinmico de reposo de polvos
de azcar en las paletas de un tambor rotativo (18). Los resultados
mostraron que el tamao de partcula, la velocidad de giro y el tipo
de material s tienen un efecto sobre el ngulo dinmico de reposo.
No obstante, cuando se tiene un control de la granulometra del
producto a la entrada del secador, y si se opera con una velocidad
angular constante, esta dependencia deja de ser relevante para el
anlisis aqu presentado.
Ms importante es la afectacin que pueda tener el contenido de
humedad en los slidos sobre el ngulo dinmico de reposo. Este
parmetro es una medida de la capacidad de flujo del producto en
la paleta. Si aumenta el contenido de humedad de la muestra al
granel, es lgico pensar que habr una mayor resistencia al

92
movimiento. No es lo mismo voltear un vaso con arena hmeda que
un vaso con arena seca.
4.3 Suposiciones sobre la Partcula Rastreada.
La propuesta de este trabajo consiste en registrar la evolucin en
las propiedades de una partcula a medida que ste avanza en el
secador, y luego extender los resultados para todo el flujo de
producto. Por lo tanto, se supondr que las propiedades de la
partcula a ser analizada son tpicas de todo el material que ingresa
al secador.
Al hablar de las propiedades, se hace referencia a la pureza,
tamao de grano, refractariedad, forma, tenacidad, porosidad,
gravedad especfica, textura, ndice de vaco, peso, entre otras.
Inicialmente, el analista debe conocer, por lo menos, la siguiente
informacin:
1. tipo de material granular (azcar, sal, madera, arena, etc.);
2. densidad de la partcula seca;
3. densidad al granel;
4. esfericidad (se asumir una esfericidad de 1);
5. distribucin de tamaos de partcula (granulometra);

93
6. variacin del ngulo dinmico de reposo en funcin del
contenido de humedad.
Tipo de material
El tipo de material granular se refiere principalmente a la resistencia
del material al impacto y a la temperatura, pegajosidad o fluidez
dentro de los niveles de humedad especificados, y facilidad de los
slidos para evacuar agua.
El conocer estas caractersticas del producto dar una pauta para
saber hasta qu punto es aplicable la teora desarrollada en este
trabajo.
Densidad real y densidad aparente
La densidad de la partcula seca es la densidad real del grano, y se
mide cuando el porcentaje de humedad en el mismo es
aproximadamente igual a cero. Esta propiedad ser necesaria para
estimar la masa slida del grano.
Por otro lado, la densidad al granel o densidad aparente es la masa
de un conjunto de granos apilados, y est directamente relacionada
con el nivel de compactacin en el pilo. Esta densidad se utilizar
principalmente para determinar la cantidad de masa presente en
una paleta en base al volumen de material en la misma.

94
Esfericidad
La esfericidad se refiere a la forma del grano o angularidad del
mismo. La esfericidad es igual a 1 para un grano perfectamente
esfrico.
En los concentrados de slice, por ejemplo, pueden existir granos
completamente irregulares en forma, as como tambin granos casi
esfricos. La partcula de madera (astilla), por otro lado, es
completamente irregular en forma.
Si bien es cierto que estos slidos no poseen forma redonda, puede
resultar vlido introducir una aproximacin a grano esfrico cuando
se analiza el problema macroscpicamente. La suposicin de grano
esfrico simplifica mucho el anlisis de movimiento. En este trabajo
se lo har de esa forma.
Distribucin de tamaos
El dimetro del grano, que pasa a ser representativo de todo el
producto en el secador, se determina como sigue:
(4.5)
95
En la ecuacin anterior, es el nmero de clases en la distribucin
de tamaos de grano, es el valor medio (marca de clase) de la
clase y es el porcentaje de las partculas (frecuencia relativa)
presentes en la clase . Esta ecuacin es similar a la ecuacin 4.4,
derivada para obtener el ngulo representativo de la cada.
A pesar de que las partculas pierden masa lquida durante su paso
por el secador, se asumir que la variacin en su dimetro
promedio es insignificante. No hay encogimiento.
Coeficiente de friccin dinmico
El aparato utilizado para medir el coeficiente de friccin dinmico es
similar al que aparece en la figura 4.11.
FIGURA 4.11 APARATO PARA MEDICIN DEL NGULO
DINMICO DE REPOSO TOMADO DE (18)

96
El aparato consiste de un tambor rotatorio horizontal que tiene
instalado un motor con ajuste variable de velocidad. En este caso,
se tienen ocho paletas distribuidas a intervalos de 45 en el interior
del tambor; estas paletas se llenan parcialmente con el material
granular a ser analizado. Durante la prueba, se toman varias
fotografas del tambor a diferentes velocidades de giro, y luego se
mide el ngulo dinmico de reposo, , y la posicin angular, , para
cada paleta en las fotografas. Para cualquier set de condiciones,
estas lecturas son analizadas en un programa de computadora que
genera el valor del coeficiente de friccin, , con la menor
desviacin estndar a partir de la ecuacin 4.1.
Se debe tener cuidado al utilizar la ecuacin de Schofield y Glikin,
puesto que pierde validez cuando se trata con materiales cuyo
ngulo dinmico de reposo vara mucho con el contenido de
humedad del producto. Deben efectuarse pruebas experimentales
para descifrar esta dependencia.

97
4.4 Movimiento Longitudinal.
Una vez encontrado el ngulo representativo y la altura de cada
correspondiente, se concentrar el anlisis en un solo grano o
partcula que cae desde dicha posicin, partiendo del reposo. Este
grano es rastreado durante todo su recorrido, tal como se muestra
en la figura 4.12, para el caso de un secador en contraflujo.
FIGURA 4.12 RECORRIDO DE LA PARTCULA EN EL SECADOR
En la figura anterior se representan los desplazamientos sucesivos
de un grano a medida que ste avanza en el secador. Cuando el
grano cae existe un pequeo rebote cuya direccin es difcil de
predecir debido a la aleatoriedad caracterstica del choque entre
partculas. En este estudio se supondr que el choque entre los
granos y la chapa metlica es perfectamente inelstico.
Consiguientemente, el nico desplazamiento que contribuye a la
longitud del secador es aquel provocado por el arrastre del aire, la
accin de la gravedad, o ambos. En los secadores donde el aire se

98
mueve en direccin opuesta al desplazamiento de los slidos, la
gravedad es el nico parmetro que contribuye al avance de las
partculas; en este caso, mientras ms inclinado est el tambor,
mayor ser la componente de la fuerza de gravedad a favor del
movimiento.
A diferencia de los secadores en flujo paralelo, la inclinacin del
tambor es un factor crtico en el diseo de los secadores en
contraflujo; no se puede prescindir de ella cuando las paletas son
paralelas al eje longitudinal del tambor.
Ya sea para un secador en flujo paralelo o en contraflujo, el clculo
de la lon it d tr ica del secador se realizar mediante la
siguiente ecuacin:
(4.6)
Donde es el nmero de ciclos de secado necesarios para que el
producto alcance cierto nivel de humedad bajo determinadas
condiciones de diseo trmico y mecnico del secador.

99
Se puede derivar una ecuacin similar para el tiempo mnimo de
secado:
(4.7)
Donde es el tiempo de vuelo de la partcula en el intervalo . Este
tiempo es el tiempo en el que los slidos se encuentran en contacto
directo con la corriente de aire. No debe confundirse con el tiempo
de residencia del producto. Para determinar el tiempo de residencia,
se debe sumar el tiempo que le toma a la partcula subir desde la
zona inferior del tambor hasta la posicin de cada. Este
desplazamiento es provocado nicamente por el ascenso de la
paleta, donde se asume un secado nulo. El tiempo de residencia
aproximado, en la longitud trmica del secador, se puede aproximar
mediante la siguiente ecuacin:

(4.8)
En la ecuacin anterior, es la separacin angular entre paletas
consecutivas.
100
4.4.1 Flujo Paralelo.
El diagrama de cuerpo libre y el diagrama de velocidades de
la figura 4.13 muestran la interaccin entre la partcula y el
aire de secado cuando este ltimo se mueve en la direccin
de avance del producto.
FIGURA 4.13 DCL DE LA PARTCULA Y DIAGRAMA DE
VELOCIDADES (FLUJO PARALELO)
El grano o partcula de masa est bajo la influencia
del campo gravitacional, , y de la fuerza de arrastre del aire,
. Se considera un sistema de coordenadas alineado con el
eje longitudinal del tambor, el cual puede tener una
inclinacin . La fuerza de arrastre depende de la magnitud
de la velocidad relativa del aire con respecto al grano, .

101
Haciendo una sumatoria de fuerzas en Z y en Y, se obtiene
( ) (4.9)
( ) (4.10)
Las componentes de la fuerza de arrastre se pueden
determinar como sigue:
donde
( ) ; ( )
En las ecuaciones anteriores, es la densidad del aire
evaluada a la temperatura flmica, es el ngulo formado
entre los vectores velocidad relativa y velocidad del aire, y
son las componentes Z y Y de la velocidad absoluta de la
partcula, es la velocidad absoluta promedio del aire en
una seccin del tambor, es el rea proyectada de la

102
partcula esfrica, y es el coeficiente de resistencia para
una esfera lisa.
Los coeficientes de resistencia dependen del nmero de
Reynolds, el cual se evala a la temperatura flmica. En la
figura 4.14 se muestran los coeficientes de resistencia de
cuerpos lisos a bajos nmeros de Mach (flujos
incompresibles).
FIGURA 4.14 COEFICIENTES DE RESISTENCIA DE

CUERPOS LISOS TOMADO DE (19)
Existen dos dificultades con respecto al manejo de las
ecuaciones 4.10 y 4.11.

103
Primero, los coeficientes de resistencia del aire dependen de
la velocidad instantnea de la partcula, la cual cambia
paulatinamente durante su cada desde la paleta hacia la
parte baja del tambor. En cada punto de la trayectoria
existir, por ende, un coeficiente de resistencia diferente que
depender del valor de la velocidad relativa entre el aire y la
partcula.
Segundo, en cada seccin del tambor o plano , las
propiedades del aire son diferentes debido a que las
condiciones de temperatura y humedad estn cambiando;
por lo tanto , y son variables en la direccin Z.
La primera dificultad se resuelve expresando el coeficiente
de resistencia en trminos del nmero de Reynolds,
haciendo uso de la curva presentada en la figura 4.14 para el
caso de una esfera lisa. Se toman varios puntos de dicha
curva y se selecciona la ecuacin que presente un mejor
ajuste. En la figura 4.15 se presenta el resultado obtenido.

104
FIGURA 4.15 CURVA DE AJUSTE PARA EL

COEFICIENTE DE RESISTENCIA DEL AIRE
De este modo, se tiene que puede ser aproximado
mediante la ecuacin 4.12. La ventaja de tener esta ecuacin
es que queda en funcin de la velocidad de la partcula.
[ ] (4.11)
donde es la viscosidad dinmica del aire a la temperatura
flmica correspondiente.
Para tratar el segundo problema, se dividir el tiempo de
vuelo de la partcula en pequeos intervalos y se resolvern
las ecuaciones de movimiento para cada uno de ellos,
incorporando los cambios en las propiedades de ambos, aire

105
y partcula, ocasionados por la transferencia de calor y masa.
El paso en el tiempo ser constante.
El problema de movimiento est totalmente resuelto cuando
se encuentra la trayectoria de la partcula. Para ello, se debe
hallar la velocidad en cada punto de la trayectoria,
resolviendo las ecuaciones diferenciales para y , e
integrar los resultados con respecto al tiempo. En el captulo
6 se expone la solucin simultnea de estas ecuaciones
mediante el uso del software EES, el cual emplea el mtodo
numrico de Runge Kutta para este efecto.
4.4.2 Contraflujo.
En la figura 4.16 se muestra el DCL y el diagrama de
velocidades la partcula en interaccin con el aire de secado
cuando este ltimo se mueve en direccin contraria al
avance del producto.

106
FIGURA 4.16 DCL DE LA PARTCULA Y DIAGRAMA DE
VELOCIDADES (CONTRAFLUJO)
El tratamiento es similar al caso del secador en flujo paralelo;
la nica diferencia radica en la direccin horizontal del vector
fuerza de arrastre. Ntese que, para que exista avance de
producto, la componente Z del peso del grano debe ser
mayor que la fuerza de arrastre provocada por la corriente de
aire.
De acuerdo a la figura anterior, la sumatoria de fuerzas en Z
y en Y quedara como sigue:
( ) (4.12)
107
( ) (4.13)
Las fuerzas de arrastre y sustentacin se determinan de la
misma forma que en el secador en flujo paralelo, con la
diferencia de que la magnitud de la velocidad relativa del aire
con respecto a la partcula es
( )
En este punto, deber notarse que la seleccin del dimetro
e inclinacin del secador debe ser tal que se evite el
estancamiento de producto. Como el flujo de aire necesario
para el secado est dado netamente por un clculo trmico,
la nica manera de jugar con la velocidad del aire, sin
interferir con el ventilador, es variando el dimetro del
tambor. Si en el diseo se selecciona un dimetro muy
pequeo, la velocidad del aire puede crecer a tal punto de
que la fuerza de arrastre en sentido contrario al avance del
producto genere el estancamiento del mismo.
Bajo las mismas circunstancias o condiciones preliminares de
diseo, el avance paulatino de una partcula en un secador
en contraflujo puede ser menor que en uno de flujo paralelo.
Este hecho, sumado a que la transferencia de calor global es

108
ms efectiva en el secado en contraflujo, conlleva al
pensamiento de que este tipo de secadores pueden llegar a
ser ms cortos que aquellos en flujo paralelo. Esto debe ser
verificado para cada aplicacin.
Cabe resaltar que la seleccin de uno u otro tipo de secador
no solo depende de cun efectiva es la transferencia de
calor, sino tambin del control de temperatura en los slidos.
La temperatura del producto se vuelve un factor importante
dependiendo de sus caractersticas o el tipo de proceso.
Con un secador en flujo paralelo, se puede controlar la
temperatura del producto de una mejor manera. No es as en
un secador en contraflujo, ya que el aire ms caliente entra
en contacto directo con los slidos que tienen el menor
porcentaje de humedad, aumentando la probabilidad de que
stos adquieran una cantidad significativa de calor sensible
que derive en el aumento de su temperatura.
En el captulo siguiente se abordar el problema de
transferencia de calor y masa en la partcula, el cual debe
resolverse simultneamente con el problema de movimiento.

109
CAPTULO 5
5. TRANSFERENCIA DE CALOR Y DE MASA
En este captulo se plantean las ecuaciones necesarias para resolver los
problemas de transferencia de calor y masa relacionados al secador
rotativo. Bsicamente, la estrategia consiste en determinar el coeficiente
de transferencia de masa a partir del conocimiento del coeficiente de
transferencia de calor, encontrado mediante las correlaciones
correspondientes al caso de flujo externo sobre una esfera. Para ello, se
utilizar la analoga de transferencia de calor y masa.
El secado de la partcula se modelar como un caso de enfriamiento
evaporativo. Solamente se considerar la transferencia de calor por
conveccin. Se ignorarn las prdidas de calor hacia el entorno, as
como tambin la transferencia de calor por radiacin desde el aire de
secado hacia la partcula.

110
5.1 Suposiciones Generales.
Una vez establecido el modelo matemtico del movimiento de los
granos en el secador, se procede a estimar el intercambio de calor y
masa experimentado durante su recorrido. Para esta parte del
estudio, se supondr lo siguiente:
1. el calor transferido desde el secador hacia el ambiente no
tiene una repercusin apreciable en el secado;
2. el calor radiante emitido por el aire en el interior del secador
es despreciable;
3. la transferencia de calor y masa ocurre solamente durante la
cada de los slidos;
4. la interaccin entre partculas adyacentes es nula.
En el secador, el aire adquiere la humedad del producto en forma
de vapor de agua. El mecanismo ms importante de vaporizacin
es la transferencia de calor convectivo desde el aire hacia los
slidos; el movimiento relativo entre el fluido y los slidos provoca
que el aire ceda calor al producto, haciendo que el agua contenida
en el mismo cambie de fase y se libere. Si se establece una
comparacin entre el calor que cede el aire hacia los slidos y el
calor transferido desde el aire hacia la chapa metlica del tambor, y
posteriormente hacia el ambiente, se encontrar que la diferencia

111
entre uno y otro es muy grande. La transferencia de calor ms
importante ocurre entre el aire y el producto hmedo. Por ende, se
considerar al secador como un medio adiabtico.
Los secadores rotativos muy pocas veces se aslan en la parte
exterior. Si bien es cierto que el aislamiento contribuye a que se
aproveche mejor el recurso energtico, el tener un tambor desnudo
ayuda a que la temperatura de la chapa metlica no crezca
demasiado, sobretodo a la entrada de aire caliente. Esto evita el
sobrecalentamiento, a veces indeseado, del producto.
El aislamiento se vuelve un elemento importante en los hornos
rotativos (que constructivamente son muy similares a los secadores
rotativos), donde la temperatura es ms elevada y la prdida de
calor por radiacin hacia el exterior tiene mayor influencia sobre el
desempeo y la economa del equipo. Aqu s se busca el
sobrecalentamiento de los slidos.
La temperatura del aire de secado en un secador rotativo
generalmente no va ms all de los 800 K. En un calentador de
llama directa que trabaje con gas natural, por ejemplo, para calentar
1 kg/s de aire de 300 a 800 K, la cantidad aproximada de
combustible requerida es de 0,012 kg/s. Suponiendo una
combustin estequiomtrica, sin exceso de aire, la concentracin de

112
vapor de agua y dixido de carbono en el aire de secado bajo estas
condiciones sera menor al 1.2%. A esta concentracin, la
emisividad de los productos de la combustin antes mencionados
es muy baja, volvindose elementos no participativos en la
transferencia de calor por radiacin.
Finalmente, los granos no solamente ganan calor mientras
interactan con la corriente de aire. Tambin existe transferencia de
calor por contacto entre granos adyacentes, o entre granos y metal,
mientras ascienden en la paleta o mientras reposan en el fondo del
tambor; y transferencia de calor por conveccin entre la capa
superior de producto en las paletas y el aire que pasa por encima
(caso convectivo de placa plana). No obstante, estos modos
secundarios de transferencia de calor no son comparables con
aquella que ocurre mientras los granos son enfrentados por el aire
en la ducha generada durante su cada. Por ese motivo, se
considerar que el nico aporte al secado de los slidos es la
transferencia de calor convectivo que ocurre durante su trayectoria
de vuelo.
113
5.2 Capa Lmite y Grupos Adimensionales de Importancia.
El concepto de capa lmite surge a partir de la necesidad de medir
la influencia que ejerce la presencia de una superficie slida sobre
un fluido, cuando este ltimo se desplaza sobre esta superficie. La
capa lmite se caracteriza por tener un espesor, el cual est dado
por la distancia entre la superficie slida y la primera pelcula fluida
que deja de ser afectada por la presencia de la superficie slida.
Existen tres tipos de capas lmite: hidrodinmica, trmica y de
concentracin. Las dos ltimas, que son las que interesan en este
captulo, son una manifestacin de la transferencia de calor y de
masa, respectivamente. La primera est directamente relacionada
con la fuerza de arrastre y sustentacin que ejerce el fluido sobre el
slido cuando existe movimiento relativo entre ellos; esto ya se
analiz en el captulo anterior.
Para el flujo sobre cualquier superficie, siempre existir una capa
lmite de velocidad y, por ello, friccin superficial. Sin embargo, una
capa lmite trmica y, de aqu la transferencia de calor por
conveccin, existe solo si difieren las temperaturas de la superficie y
del flujo libre (flujo no afectado por la superficie). De manera similar,
una capa lmite de concentracin y la transferencia de masa por
conveccin existen solo si la concentracin superficial de una
superficie difiere de su concentracin en el flujo libre. En el caso

114
particular del secador rotativo, estn presentes las tres capas lmite
cuando los slidos interactan con la corriente de aire de secado.
El estudio de la capa lmite se simplifica significativamente mediante
el uso de grupos adimensionales tales como Nusselt (Nu), Prandtl
(Pr), Reynolds (Re), Lewis (Le), entre otros. En Apndices se
presenta un resumen de los grupos adimensionales ms
importantes relacionados con la transferencia de calor y de masa.
Estos grupos adimensionales se relacionan entre ellos mediante
correlaciones que han sido desarrolladas a travs de estudios
experimentales efectuados por varios investigadores con el pasar
de los aos. Su uso est dirigido principalmente a la determinacin
de los coeficientes promedio de transferencia de calor y de masa,
y respectivamente.
5.3 Transferencia de Calor.
La transferencia de calor entre el aire y los slidos en el secador
ocurre principalmente debido a la conveccin forzada en flujo
externo alrededor de partculas con dimetro igual al dimetro
promedio de grano. Existen numerosas correlaciones de

115
transferencia de calor propuestas para este caso. Whitaker (20)
recomienda una expresin de la forma
( ) ( )
(5.1)
Donde
Solamente la viscosidad dinmica del aire, , se evala en la
temperatura superficial del grano. Las dems propiedades se
evalan en , que es la temperatura del aire hmedo fuera de la
capa lmite trmica.
El coeficiente de transferencia de calor por conveccin puede
determinarse a partir del nmero de Nusselt promedio mediante la
siguiente expresin:

(5.2)
En la ecuacin anterior, el coeficiente de conduccin para el aire
hmedo, , tambin se evala en .

116
Para poder hallar , es necesario resolver simultneamente las
ecuaciones de transferencia de calor y las ecuaciones de
movimiento. Lo anterior se debe a que el nmero de Nusselt
depende del nmero de Reynolds, el que a su vez est en funcin
de la velocidad relativa entre la partcula y la corriente de aire. En el
captulo 6 se presenta una descripcin detallada de la solucin a
este problema mediante EES.
5.4 Transferencia de Masa.
Existen dos modos de transferencia de masa: conveccin y difusin.
El primero es anlogo a la transferencia de calor por conveccin y el
segundo a la transferencia de calor por conduccin.
As como la transferencia de calor se debe a una diferencia de
temperaturas, la transferencia de masa es ocasionada por un
desequilibrio de concentraciones. En este texto, cuando se hable de
concentraciones, solamente se har referencia a dos especies, y
, donde es la especie en mayor proporcin. La especie
normalmente es un vapor que se transfiere en un flujo de gas
debido a la evaporacin o sublimacin de una superficie lquida o
slida, respectivamente. As, en el caso del secado, el aire

117
correspondera a la especie gaseosa mientras que el vapor de
agua a la especie .
En el secador rotativo, la masa se transfiere primordialmente por
conveccin. La figura 5.1 muestra un slido de superficie arbitraria
con una concentracin superficial , que est siendo
enfrentado por una corriente de aire que se mueve a una velocidad
y que tiene una concentracin ; la concentracin de agua es
un parmetro que est relacionado directamente con la humedad
absoluta del aire.
FIGURA 5.1 TRANSFERENCIA DE MASA
Asumiendo que la concentracin de vapor de agua en la superficie
lquida es mayor que la concentracin de vapor de agua en la
corriente de aire, la transferencia de masa por unidad de rea
diferencial est dada por

118
( )
(5.3)
Donde es el coeficiente puntual de transferencia de masa por
conveccin. Las unidades de las variables presentadas
anteriormente son
{ } { }
{ } { }
{ } { } { }
Si se integra el coeficiente puntual de transferencia de masa en
toda la superficie , se obtiene lo siguiente:
( ) (5.4)
Donde
Al multiplicar la ecuacin 5.4 por el peso molecular de la especie ,
se obtiene
( ) (5.5)
119
Esta ltima expresin tiene mayor valor para el lenguaje de uso
comn, ya que al trabajar con las densidades de , tanto para la
superficie lquida como para la regin fuera de la capa lmite,
est dada directamente en .
La cantidad representa la tasa de evaporacin de agua desde
una sola partcula. Para conocer la tasa total de evaporacin, es
necesario establecer el flujo de partculas que pasa por el secador.
Si la densidad y masa del grano seco son y
respectivamente, y la humedad inicial del producto es , la masa
del grano a la entrada del secador se puede estimar como sigue:
(5.6)
Donde
( )
Si el flujo msico de producto hmedo antes del secado es , el
flujo de partculas atravesando cualquier seccin del secador est
dado por
(5.7)
[ ]
120
La ecuacin anterior es equivalente a
( )

( )
Parte de la estrategia consiste en discretizar el anlisis de
movimiento, transferencia de calor y masa en intervalos de tiempo
fijos. As, la cantidad de agua transferida por una partcula en cada
intervalo de tiempo sera
( ) (5.8)
Donde es una cantidad constante. Haciendo uso de
la ecuacin 5.7, se puede extender el resultado anterior al flujo total
de partculas:
(5.9)
La cantidad representa la tasa total de evaporacin en un paso
de tiempo.
De acuerdo a la ecuacin 5.5, para determinar es necesario
conocer , y . El primer parmetro puede ser deducido a
partir de la siguiente expresin:
(5.10)
121
En la ecuacin anterior, es la humedad relativa del aire fuera de
la capa lmite y es la densidad del vapor de agua saturada
evaluada a la temperatura del aire. Este valor de densidad se puede
obtener fcilmente mediante el recproco del volumen especfico en
la curva de vapor saturado, a la temperatura de saturacin
correspondiente.
Para determinar la densidad del vapor de agua en la superficie del
slido, se considerarn las siguientes suposiciones:
1. Existe equilibrio termodinmico en la interfaz entre el gas y el
lquido. Por ende, la temperatura del vapor en la interfaz es
igual a la temperatura superficial, .
2. El vapor se encuentra en estado saturado; por lo tanto, se
pueden utilizar las tablas termodinmicas del agua para obtener
la densidad a partir de . Por ejemplo, de la tabla se leer el
volumen especfico , lo que lleva a .
Finalmente, el coeficiente promedio de transferencia de masa por
conveccin puede obtenerse de dos maneras. La primera consiste
en utilizar la correlacin de Whitaker, remplazando y con los
nmeros de Sherwood y Schmidt, y respectivamente (20).
Esta alternativa es viable siempre y cuando se cuente con tablas

122
que contengan informacin sobre los parmetros correspondientes
a la transferencia de masa.
La segunda alternativa, que puede resultar ms conveniente,
consiste en aplicar la analoga de la transferencia de calor y masa.
Este mtodo se discute a continuacin.
5.5 Analoga de la Transferencia de Calor y Masa.
Si dos o ms procesos estn gobernados por ecuaciones
adimensionales de la misma forma, se dice que los procesos son
anlogos (21). Las transferencias de calor y de masa por
conveccin son anlogas. Por lo tanto, las relaciones
adimensionales que gobiernan el comportamiento de la capa lmite
trmica deben ser las mismas que las que gobiernan la capa lmite
de concentracin. Dicho de otro modo, las relaciones de
transferencia de calor y masa para una geometra particular son
intercambiables.
Lo anterior se refleja en la siguiente relacin:
(5.11)

123
En la literatura muchas veces es difcil encontrar parmetros de los
fluidos relacionados con la transferencia de masa; no ocurre lo
mismo con aquellas propiedades asociadas a la transferencia de
calor. El resultado anterior sirve para determinar uno de los
coeficientes de conveccin, en este caso, a partir del
conocimiento del otro coeficiente, y es aplicable tanto para flujo
turbulento como para flujo laminar. Para la mayor parte de las
aplicaciones es razonable suponer un valor (21).
El coeficiente de difusin binario, , es un parmetro que
relaciona a la transferencia de masa por conduccin con el
gradiente de concentracin de A en B sobre una superficie lquida o
slida (Ley de Fick); este coeficiente depende de la temperatura. En
este caso, para vapor de agua en aire, se obtiene mediante la
ecuacin 5.12, proporcionada por Nellis y Klein (22).
[ ] [ ]
(5.12)
[ ]
Donde la temperatura flmica, , es el promedio entre la
temperatura superficial del slido (o envolvente lquida de agua) y la
temperatura del aire de secado.

124
El nmero de Lewis, , se determina de la siguiente forma:
(5.13)
En casi todos los problemas prcticos, la concentracin de la
especie es mucho mayor que la concentracin de la especie
( ). Por lo tanto, es razonable suponer que las propiedades
de la capa lmite ( , , , , , etctera) son las propiedades de la
especie .
Una aplicacin importante de la analoga de transferencia de calor y
masa tiene que ver con el proceso de enfriamiento evaporativo, que
ocurre cada vez que un gas fluye sobre un lquido.
5.6 Enfriamiento Evaporativo.
En este estudio, el proceso de secado se va a modelar como un
caso de enfriamiento evaporativo. Si se analiza a todo el secador
como un volumen de control, lo anterior sera equivalente a un
proceso de saturacin adiabtica, considerando prdidas de calor
despreciables hacia el entorno.
La evaporacin debe ocurrir a partir de la superficie del lquido, y la
energa asociada con el cambio de fase es el calor latente de

125
vaporizacin del lquido. La energa que se requiere para mantener
la evaporacin debe provenir de la energa interna del lquido, que
entonces experimenta una disminucin de temperatura. Sin
embargo, en condiciones de estado estable, la energa latente
perdida por el lquido debido a la evaporacin debe recuperarse
mediante la transferencia de energa al lquido desde los
alrededores, tal como se muestra en figura 5.2 (21). Considerando
que los efectos de radiacin son insignificantes, esta transferencia
se debe a la conveccin de energa sensible del aire.
FIGURA 5.2 INTERCAMBIO DE CALOR LATENTE Y SENSIBLE
EN UNA INTERFAZ GAS-LQUIDO COPIADO DE (21)
Si no se agrega calor por otros medios, el balance entre la
transferencia de calor por conveccin desde el gas y la prdida de
calor evaporativo desde el lquido se expresa como

126
( ) [ ( ) ]
Donde la densidad de vapor en la superficie es la que se asocia con
las condiciones saturadas en . Por ello, la magnitud del efecto de
enfriamiento se expresa como

( )[ ( ) ]
(5.14)
El factor se obtiene a partir de la ecuacin 5.14, que est
basada en la analoga de la transferencia de calor y de masa.
Las ecuaciones presentadas en este apartado funcionan siempre
que exista una pelcula de agua abrigando a los slidos. Se supone
que esto es as en el caso del secador rotativo.

127
CAPTULO 6
6. IMPLEMENTACIN DE SOFTWARE DE ANLISIS
DE SECADORES ROTATIVOS
En el presente captulo se exponen cada una de las etapas de este
trabajo desarrolladas en EES, dando una breve explicacin de las
funciones, sentencias y algoritmos utilizados.
6.1 EES (Engineering Equations Solver).
En el rea de termo-fluidos, EES es una herramienta muy usada
porque, aparte de su cualidad de solver, cuenta con funciones
predefinidas que permiten acceder a las propiedades
termodinmicas y de transporte de muchas sustancias durante la
solucin de ecuaciones. Adems, EES toma capacidades de
programacin estructurada, permitiendo el uso de sentencias

128
lgicas como ifthenelse o repeatuntil, que son tiles cuando
se requiere ejecutar operaciones repetitivas.
Para un mejor manejo y organizacin en la programacin, se ha
dividido el problema del secador rotativo en libreras, programas,
sub-programas y funciones. La figura 6.1 es un diagrama de flujo
que muestra la secuencia e anlisis utilizada.

129
FIGURA 6.1 DIAGRAMA DE FLUJO PARA EL
ANLISIS DE SECADORES ROTATIVOS EN ESS

130
6.2 Anlisis Termodinmico.
Los datos de entrada o inputs para esta parte del anlisis son las
humedades de entrada y salida del producto, que pueden ser
expresadas en base seca o en base hmeda, el flujo msico de
producto a secar o de producto seco, la humedad relativa y
temperatura del aire ambiental, la temperatura del aire a la salida
del calentador, y la humedad relativa del aire despus del secado.
Si el anlisis trata sobre un secador existente, y no se conoce
alguno de estos parmetros, hay que medirlo o estimarlo de algn
modo. En EES, la escritura de las variables de entrada se hace de
la siguiente manera:
X_1=xxx {Humedad inicial del producto en base hmeda}

X_2=xxx {Humedad final del producto en base hmeda}
m_DOT_p=xxx*CONVERT(kg/h;kg/s) {Flujo msico del producto a secar}
phi_1=xxx {Humedad relativa del aire ambiental}
T_1=xxx {Temperatura aire ambiental}
T_2=xxx {Temperatura del aire a la salida del calentador}
phi_3=xxx {Humedad relativa del aire despus del secado}
P1=xxx*CONVERT(psi;kPa) {Presin de trabajo}
Con esta informacin, se ingresan las ecuaciones expuestas en el
resumen presentado al final de la seccin 3.3. No es necesario
seguir el orden que all se especifica; sin embargo, es aconsejable
hacerlo para un mejor entendimiento.

131
m_DOT_s=(1-X_1)*m_DOT_p
DELTAm_DOT_w=m_DOT_s*(1/(1-X_1)-1/(1-X_2))
w_1=HUMRAT(AirH2O;T=T_1;P=P1;R=phi_1)
w_2=w_1
h_2=ENTHALPY(AirH2O;T=T_2;P=P1;w=w_2)
phi_2=RELHUM(AirH2O;T=T_2;P=P1;w=w_2)
h_3=h_2
T_3=TEMPERATURE(AirH2O;R=phi_3;P=P1;h=h_3)
w_3=HUMRAT(AirH2O;h=h_3;P=P1;T=T_3)
DELTAw=w_3-w_2
m_DOT_aire=DELTAm_DOT_w/DELTAw
rho_aire_3=DENSITY(AirH2O;R=phi_3;P=P1;w=w_3)
rho_aire_2=DENSITY(AirH2O;T=T_2;P=P1;w=w_2)
V_DOT_aire_out=(m_DOT_aire/rho_aire_3)*CONVERT(m^3/s;cfm)
V_DOT_aire_in=(m_DOT_aire/rho_aire_2)*CONVERT(m^3/s;cfm)
cp_aire1=SPECHEAT(AirH2O;T=T_1;P=P1;R=phi_1)
cp_aire2=SPECHEAT(AirH2O;T=T_2;P=P1;R=phi_2)
cp_aire=(cp_aire1+cp_aire2)/2
Q_DOT=m_DOT_aire*cp_aire*(T_2-T_1)
Los caracteres entre llaves son comentarios que ayudan a
identificar cada variable o lnea de programacin. Las letras o
nmeros escritos despus de un subguin representan el subndice
de la variable que lo precede. La palabra DOT seguida de un
subguin indica una variacin en funcin del tiempo. La funcin
CONVERT sirve para convertir las unidades de la variable
antecesora. Por ejemplo, en el caso de la variable , la funcin
CONVERT arroja un factor multiplicador que sirve para transformar
sus unidades de kilogramos por hora a kilogramos por segundo.
As, hay muchas otras funciones predefinidas como HUMRAT,
ENTHALPY, TEMPERATURE, DENSITY, SPECHEAT, cuyos

132
valores de salida son obvios. Cada una de ellas, por supuesto,
requiere de ciertos parmetros de entrada, colocados entre
parntesis, que hacen referencia a las condiciones o estado de la
sustancia.
En total, se han planteado veintisiete ecuaciones para veintisiete
variables existentes en este programa, el cual se almacena con el
nombre de Anlisis Termodinmico.EES.
Adems del flujo msico de aire, los resultados o outputs tambin
incluyen las condiciones del aire antes y despus del secado, la
cantidad de agua evaporada, el flujo volumtrico de aire a la
entrada y salida del secador, el calor total agregado al aire en el
calentador y el flujo msico de producto seco.
Con esta breve pero necesaria descripcin del uso de EES, el lector
podr captar de forma ms rpida la codificacin presentada en las
secciones posteriores, dirigiendo su concentracin ms hacia la
lgica y estrategia de programacin empleada.

133
6.3 Masa de Producto en las Paletas en Funcin del Tiempo.
Este anlisis toma lugar en la seccin transversal del tambor. De
acuerdo a lo desarrollado en el captulo 4, el rea transversal de la
cama de partculas transportada por una paleta est en funcin del
ngulo dinmico de reposo, de la geometra de la paleta y de su
posicin angular.
El objetivo de cuantificar la cantidad de material en las paletas en
funcin del tiempo es poder estimar el ngulo representativo y la
altura promedio de cada de las partculas.
Para ello, se ha creado una funcin en EES cuyo nombre es
( ). Esta funcin puede ser almacenada en la librera de EES
como un archivo de tipo .LIB. La ventaja de hacerlo de esta forma
es que, as como cualquier otra funcin predefinida (ENTHALPY por
ejemplo), sta puede ser llamada o utilizada desde cualquier otro
programa con extensin .EES.
Las funciones son generadas en la pantalla de trabajo habitual de
EES, en cuya primera lnea debe necesariamente aparecer la
palabra FUNCTION seguida del nombre de la funcin y de la
variable independiente entre parntesis; la estructura finaliza con la
palabra END:
134
FUNCTION masa(theta)
.
.
.
END
A diferencia del formato utilizado en Anlisis Termodinmico.EES,
aqu es obligatorio que las lneas de cdigo posean un orden lgico
y, adems, las ecuaciones ingresadas deben ser expresadas de
forma explcita, es decir, deben empezar con una variable seguida
del signo igual.
A continuacin se detallan los inputs de la funcin ( ):
FUNCTION masa(theta)
"DATOS GEOMTRICOS Y CINEMTICOS"
l_1=xxx {Longitud del labio}
l_2=xxx {Longitud de la base de la paleta}
alpha_A=xxx {ngulo entre segmentos de paleta}
d=1 {Ancho de la paleta}
gamma=0 {OFFSET}
R=xxx {Radio del tambor}
omega=xxx*CONVERT(rpm;rad/s) {Velocidad de giro del tambor}
g=9,81 {Aceleracin de la gravedad}
"DATOS DEL MATERIAL"
mu=xxx {Coeficiente de friccin dinmico}
rho_BULK=xxx {Densidad al granel promedio}
{CONTINA}
Los datos geomtricos incluyen las dimensiones de la paleta y el
radio o dimetro interior del tambor. En este caso se ha tomado el
ejemplo de una paleta conformada por dos segmentos, base y labio.
Algunos investigadores han desarrollado procedimientos
geomtricos para trabajar con paletas de tres o ms segmentos, y

135
paletas de seccin circular inclusive (17). Lo nico que cambia es la
manera de determinar el rea transversal de la cama de material.
El coeficiente de friccin dinmico y la densidad al granel son las
nicas propiedades del producto que intervienen en esta parte del
anlisis.
El coeficiente de friccin dinmico es una medida de la fluidez entre
capas adyacentes de producto. Este coeficiente se introduce en la
ecuacin 4.1 para hallar el ngulo dinmico de reposo. No obstante,
la ecuacin deducida por Schofield y Glikin es vlida solamente
para materiales no pegajosos que poseen buena fluidez (9). Aqu
surge la interrogante de cmo saber si un material determinado
cumple con estas caractersticas, o hasta qu punto las satisface.
Una forma es ejecutando ensayos similares al de la figura 4.11.
Bsicamente, la prueba consiste en hacer girar un tambor con
producto en su interior. El ensayo se realiza generalmente a
diferentes niveles de humedad, velocidades angulares, dimetros
de tambor, y tamaos de partcula. Este ltimo parmetro de control
sirve para ver el efecto que tiene la compactacin del producto
sobre el ngulo dinmico de reposo.
A medida que el tambor gira, se toman fotografas consecutivas de
la seccin transversal y se las enva a una computadora para medir

136
el ngulo dinmico de reposo en cada posicin angular de la paleta.
Luego, se introducen estos valores en la ecuacin 4.1 y se
determina el coeficiente de friccin dinmico en cada caso. Si este
coeficiente resulta ser aproximadamente constante, entonces el
material dado satisface la relacin de Schofield y Glikin. Por otro
lado, si existe una alta variabilidad en los resultados, simplemente
se concluye que el material no posee la fluidez suficiente. En ese
caso, la ecuacin 4.1 y el coeficiente de friccin dinmico son
descartados del anlisis y, en su lugar, se ingresa como valor de
entrada el ngulo dinmico de reposo.
El ngulo dinmico de reposo depende fuertemente del contenido
de humedad en el producto ms que de cualquier otro
parmetro (11) (12). Tambin se ha identificado que este ngulo es
aproximadamente constante para cualquier posicin angular de
paleta (7). Para materiales que fluyen de manera errtica y en forma
de grumos, siempre es importante conocer la variacin del ngulo
de reposo en funcin del contenido de humedad, previo al anlisis
de desempeo del secador rotativo. Aqu se recomienda graficar el
ngulo dinmico de reposo con respecto al contenido de humedad,
en base a pruebas experimentales, y a partir de all establecer una
relacin o regla de correspondencia que pueda ser utilizada en la
programacin.
137
Si la diferencia entre los niveles de humedad a la entrada y salida
del secador es pequea, un promedio sera suficiente. Tambin
puede ocurrir que la aleatoriedad de los resultados, debido al
carcter pegajoso del producto, no permita establecer una relacin
precisa entre el ngulo dinmico de reposo y su contenido de
humedad. Es muy importante que cada material sea evaluado en
cuanto a su calidad de fluidez antes de empezar el anlisis de
secado.
Otra propiedad que depende de la humedad en los slidos es la
densidad aparente. Si se desea conocer la masa de material en la
paleta, entonces la densidad al granel es importante. Sin embargo,
no debe perderse de vista que lo que se est determinando
realmente es un rea transversal de material que vara de acuerdo
a la posicin angular. El nico parmetro que interviene en esta
variacin es el ngulo dinmico de reposo, mas no la densidad
aparente.
Para estimar la masa en la paleta, lo que se hace es multiplicar el
rea transversal de la cama de material, asumida uniforme, por la
profundidad o ancho de paleta para obtener el volumen. Luego, se
multiplica dicho volumen por la densidad al granel y se obtiene un
estimado de la masa instantnea de producto.

138
Sea que se calcule la masa de material o no, el ngulo
representativo y la altura promedio de cada dependen nica y
exclusivamente de la tasa de variacin de rea de la cama de
material en la paleta. Dicho de otra manera, para este anlisis no
importa si se conoce o no el ancho de paleta ni la densidad al
granel, siempre y cuando se tenga claro cmo vara el ngulo
dinmico de reposo en las diferentes secciones del tambor.
En este trabajo se utilizar un ancho de paleta unitario. Es decir, lo
que la funcin ( ) arrojar realmente es una masa por unidad
de longitud de paleta. Se pudo haber trabajado solamente con el
rea, pero se prefiri hacerlo de esta manera para no perder de
vista el sentido fsico del problema. De una u otra manera, los
resultados referentes a la altura promedio de cada seran los
mismos.
Una vez aclarados los puntos relacionados con las variables de
entrada de la funcin ( ), se puede presentar ahora el cdigo
correspondiente a las ecuaciones de la seccin 4.1. Para el caso de
un material que obedece la relacin de Schofield y Glikin, se tiene

139
{CONTINUACIN}
phi=arctan((mu+(R*omega^2/g)*(cos(theta)-mu*sin(theta)))/(1-
(R*omega^2/g)*(mu*cos(theta)+sin(theta))))
x_A=l_1
y_A=0
x_B=l_1-l_2*cos(alpha_A)
y_B=l_2*sin(alpha_A)
IF (phi+delta>arctan(y_B/x_B)) THEN
X_O_g=R_O*cos(theta)
Y_O_g=R_O*sin(theta)
a=(tan(phi))^2+1
b=2*tan(phi)*(Y_O_g-X_O_g*tan(phi))
c=(X_O_g*tan(phi))^2-2*X_O_g*Y_O_g*tan(phi)+Y_O_g^2-R^2
X_W_g=(-b+SQRT(b^2-4*a*c))/(2*a)
Y_W_g=SQRT(R^2-X_W_g^2)
y_W=Y_W_g*cos(delta)+X_W_g*sin(delta)-R_O*sin(theta+delta)
x_W=(R_O*sin(theta)+y_W*cos(delta)-Y_W_g)/sin(delta)
beta=phi+delta-arctan(y_B/x_B)
OW=SQRT(x_W^2+y_W^2)
A_1=0,5*l_1*l_2*sin(alpha_A)
A_2=0,5*OB*OW*sin(beta)
AREA=A_1+A_2
ENDIF
IF (phi+delta<=arctan(y_B/x_B)) THEN
x_P=(l_1*tan(alpha_A))/(tan(phi+delta)+tan(alpha_A))
y_P=-tan(alpha_A)*(x_P-l_1)
AP=SQRT((x_P-l_1)^2+y_P^2)
AREA=0,5*l_1*AP*sin(alpha_A)
IF (phi+delta<=0) THEN
AREA=0
ENDIF
ENDIF
V=AREA*d {Volumen de material en la paleta}
masa=rho_BULK*V {Masa de material en la paleta}
END
140
Si se tratase de un material que no obedece la relacin de Schofield
y Glikin, simplemente se remplaza la primera lnea de cdigo por
phi=phi_0
donde es una constante igual al ngulo dinmico de reposo
promedio en la seccin. Este valor debe ser hallado
experimentalmente para cada producto.
Como se mencion previamente en el captulo 4, existen dos etapas
en el vaciado de las paletas: material en contacto con paleta y
tambor, y material en contacto con paleta solamente. Otra etapa
que est implcita es aquella en la que la paleta est completamente
vaca. Para diferenciar, o ms bien detectar, cada una de estas
etapas en el programa, se han utilizado condicionales IF que
evalan el valor del ngulo . Estas condiciones fueron
oportunamente deducidas en el captulo 4.
Para conocer el ngulo de vaciado de la paleta, se ha diseado otra
funcin muy simple en un programa llamado theta_vaciado.EES.
Esta funcin consta de un contador i que vara grado a grado el
valor de , evaluando la funcin ( ) hasta que su output sea
igual a cero. En ese momento, el ngulo de vaciado, , toma
el valor de la variable contadora i. Para ejecutar este procedimiento

141
se ha utilizado el lazo de repeticin REPEAT UNTIL, tal como se
detalla a continuacin.
FUNCTION theta_vaciado(theta_0)
i=theta_0
REPEAT
m_paleta=masa(i)
i=i+1
UNTIL m_paleta=0
theta_vaciado=i-1
END
theta_0=0
theta_vaciado=theta_vaciado(theta_0)
La variable de entrada, , es el ngulo donde inicia la cada de
material desde la paleta y coincide con el ngulo que optimiza a la
funcin rea, es decir, es el ngulo en el cual la cantidad de
material en la paleta es mxima. Al derivar la funcin de rea con
respecto a e igualar el resultado a cero, se encontrar que .
Es muy importante actualizar los datos de la funcin ( ) para
cada tipo de secador y material con que se trabaje. Dado que esta
funcin se almacena en la librera de funciones de EES, cuya
ubicacin no est a simple vista del usuario, puede ser un error
comn actualizar los datos de entrada de todos los programas y
olvidar actualizar los de esta funcin en particular.

142
6.4 ngulo Representativo y Altura Promedio de Cada.
Esta seccin se enfoca principalmente en la evaluacin de las
ecuaciones 4.4 y 4.5. La primera estima el ngulo representativo y
la segunda la altura promedio de cada.
Como de costumbre, el programa Altura promedio de cada.EES
inicia con la declaracin de las variables de entrada:
"DATOS DE ENTRADA"
l_1=xxx {Longitud de labio}
l_2=xxx {Longitud de paleta}
alpha_A=xxx {ngulo entre segmentos de paleta}
gamma=xxx {OFFSET}
R=xxx {Radio del tambor}
omega=xxx*CONVERT(rpm;rad/s) {Velocidad angular del tambor}
theta_vaciado=xxx {ngulo de vaciado}
DELTAtheta=1/xxx {Step angular en grados}
{CONTINA}
El ngulo de vaciado se calcula en el programa theta_vaciado.EES.
La variable representa el paso angular para la integracin
numrica que toma lugar en la ecuacin 4.5; por supuesto, mientras
ms pequeo sea el valor de , ms preciso ser el resultado
arrojado por esta ecuacin.
En este caso, para determinar el ngulo de cada representativo de
un ciclo o revolucin de tambor, se utiliz un conjunto de arreglos
mediante el comando DUPLICATE, formulado como sigue:

143
{CONTINUACIN}
div=1/DELTAtheta
N=theta_vaciado*div
theta[0]=0
f[0]=0
producto[0]=0
theta_BAR[0]=0
DUPLICATE i=0;N
theta[i+1]=theta[i]+DELTAtheta
f[i+1]=(masa(i/div)-masa((i+1)/div))/masa(0)
producto[i+1]=f[i+1]*theta[i+1]
theta_BAR[i+1]=theta_BAR[i]+producto[i+1]
END
theta_BAR=theta_BAR[N]
{CONTINA}
En el cdigo anterior, representa el nmero de filas en cada
arreglo y es el nmero de divisiones que existen en un grado
(es el inverso del paso angular). En total se han definido cuatro
arreglos:
[ ] = ngulo correspondiente a la posicin ;

[ ] = fraccin de partculas que han abandonado la paleta
entre [ ] y [ ];
[ ] = producto entre las variables [ ] y [ ];
[ ] = variable acumuladora o sumatoria de productos.
El valor final de , o sea [ ], es el ngulo representativo de toda la
cada de slidos. Ntese que el comando DUPLICATE es una
abreviacin de REPEAT UNTIL.

144
Una vez estimado el ngulo representativo de la cada, se procede
a utilizar la ecuacin 4.6 para determinar la altura promedio de
cada, que no es ms que la altura de cada correspondiente al
ngulo previamente calculado. Para ello, se sigue el siguiente
procedimiento:
{CONTINUACIN}
OB=SQRT(l_1^2+l_2^2-2*l_1*l_2*cos(alpha_A))
alpha_B=arcsin(l_1*sin(alpha_A)/OB)
R_O=SQRT(OB^2+R^2-2*OB*R*cos(gamma+alpha_B))
X_O_g=R_O*cos(theta_BAR)
Y_O_g=R_O*sin(theta_BAR)
Y_BAR=(Y_O_g+SQRT(R^2-X_O_g^2))
La variable es la distancia desde la punta del labio de la paleta
hacia el centro del tambor (u origen del SCG). Por otro lado, y
representan las coordenadas globales no rotatorias de la punta
del labio. Se aadi el subndice a las variables anteriores para
diferenciarlas de y , que son las coordenadas del mismo punto
en el SCL; EES no es case sensitive (no diferencia entre
maysculas y minsculas).
Aqu finaliza el programa Altura promedio de cada.EES. Hasta el
momento, todos los esfuerzos perpetrados en este captulo han sido
dirigidos exclusivamente a la determinacin de , que es uno de los
datos de entrada para los programas de secado a desarrollarse en
las secciones siguientes.

145
6.5 Secado en Flujo Paralelo.
En esta seccin se resuelven de forma simultnea los problemas de
transporte, transferencia de calor y transferencia de masa para una
partcula que es representativa de todo el flujo de producto.
La estrategia consiste en suponer que una partcula cae desde una
altura igual a con una velocidad inicial nula. A medida que sta
desciende por gravedad, se calcula la transferencia de calor y de
masa generada por la interaccin con la corriente de aire caliente,
cuya velocidad es determinada a partir de sus propiedades, flujo
msico y dimetro del tambor. Simultneamente, se resuelven las
ecuaciones de movimiento que definen la trayectoria de la partcula,
influenciada a su vez por el arrastre del aire.
Se asume que, al final de la cada, la partcula choca de forma
inelstica con el fondo del tambor, donde es inmediatamente
recogida por otra paleta para repetir el ciclo de ascenso y cada.
El programa principal, denominado Secado en Flujo Paralelo.EES,
luce tal como se muestra en la parte inferior. Aqu se declaran las
variables de entrada y se ejecuta la funcin .

146
X_1=xxx {Humedad inicial de la partcula en base

hmeda}
X_2=xxx {Humedad final de la partcula en base hmeda}
T_s=xxx {Temperatura inicial de la partcula, C}
rho_grano_seco=xxx {Densidad del grano seco, kg/m3}
D_grano=xxx {Dimetro promedio de grano, m}
m_DOT_p=xxx*CONVERT(kg/h;kg/s){Flujo msico del producto a secar}
m_DOT_aire=xxx {Flujo msico de aire, kg/s}
T_2=xxx {Temperatura del aire antes del secado, C}
phi_2=xxx {Humedad relativa del aire antes del secado}
P1=xxx*CONVERT(atm;kPa) {Presin de trabajo, kPa}
R=xxx {Radio del tambor, m}
alpha=xxx {Inclinacin del tambor, }
Y_BAR=xxx {Altura promedio de cada, m}
step=xxx {Paso en el tiempo, s}
L_secador=L_secador(X_1;X_2;T_s;rho_grano_seco;D_grano;m_DOT_p;
m_DOT_aire;T_2;phi_2;R;alpha;Y_BAR;step)
Los parmetros y corresponden al contenido de humedad del
producto antes y despus del secado. En este caso, estos valores
se expresan en base hmeda y son los mismos que se utilizaron
para el anlisis termodinmico efectuado en el apartado 6.2. Se
asume que el contenido de humedad de la partcula es igual al del
producto que se desplaza en el secador en el mismo instante.
El siguiente valor de entrada es la temperatura superficial inicial de
la partcula, . En realidad, este es un valor referencial solamente.
Puede omitirse del programa sin generar consecuencias sobre los
resultados finales. En este trabajo particularmente, se supuso que la
transferencia de calor sensible que ganan los slidos, antes de que
su superficie llegue a la temperatura inicial de evaporacin, ocurre

147
en un pequeo tramo de cilindro que se asume despreciable. En
otras palabras, el anlisis se efectuar considerando que hay un
salto escaln de temperatura superficial en el producto a la entrada
del secador, el cual va desde su temperatura antes del secado, que
no necesariamente es la temperatura ambiente, a su temperatura
inicial de evaporacin. De este modo, se supone que todo el
proceso en el interior del secador es netamente evaporativo.
La densidad del grano seco es importante porque permitir estimar
la masa de la partcula cuando el contenido de humedad en la
misma es cero. Esta densidad es la densidad de un grano
individual, mas no la densidad al granel o aparente.
El dimetro promedio de grano se calcula de acuerdo a la
distribucin de tamaos determinada a partir del cribado de las
partculas. Como se mencion en su momento, se trabajar bajo la
suposicin de que la partcula es esfrica y que su dimetro es
constante durante el secado. El dimetro promedio de grano no
solamente sirve para determinar el volumen de la partcula y con
ello su masa, sino tambin para estimar el nmero de Reynolds y
posteriormente el coeficiente de transferencia de calor a partir del
nmero de Nusselt.
148
El flujo msico de producto a secar es el flujo de producto hmedo
antes de su ingreso al secador. No debe confundirse ni con el flujo
de producto seco ni con el flujo msico que se obtiene a la salida
del secador. En este caso, el flujo de producto hmedo se utiliza
para estimar el flujo de partculas, en unidades/segundo, que pasan
por el secador bajo circunstancias estables; este clculo se realiza
mediante las ecuaciones 5.6 y 5.7.
El flujo msico de aire, , y sus propiedades a la salida del
calentador, y , se determinan mediante el programa Anlisis
Termodinmico.EES. Se considera que el flujo msico de aire es
constante durante su paso por el secador, y se utiliza para
determinar su velocidad en cada seccin a partir del rea
transversal y de su densidad, la que a su vez depende de la
temperatura, humedad relativa (o absoluta), y de la presin de
trabajo, que en este caso se considera constante e igual a la
presin baromtrica. Las propiedades del aire se actualizan a
medida que ste va extrayendo la humedad de los slidos.
Los nicos parmetros geomtricos del secador que intervienen en
esta parte del anlisis son el radio interior del tambor, , y su
inclinacin con respecto a la horizontal, . El segundo tiene

149
influencia directa sobre el avance del producto en el interior del
cilindro.
La altura de cada promedio, , se obtiene del programa Altura
promedio de cada.EES. En este caso particular, se ha considerado
una altura de cada constante a lo largo del secador. Sin embargo,
hay ocasiones en que no se tienen paletas de igual geometra en
toda la longitud del secador, o que el ngulo dinmico de reposo del
producto presenta variaciones considerables entre los niveles de
humedad de entrada y salida. En esos casos, habra que considerar
una altura promedio de cada variable durante la ejecucin del
programa; lo ms simple es establecer una regla de
correspondencia que relacione la altura promedio de cada con el
contenido de humedad del producto.
Finalmente, el step es el paso temporal en segundos. El step se
utiliza en la discretizacin de los clculos; mientras ms pequeo es
su valor, la precisin en los resultados incrementa.
El programa principal est constituido por una funcin llamada
, la que a su vez contiene a un sub-programa denominado
SECADO. Ambos poseen con sus propios valores de entrada y
salida.
150
En el caso de la funcin , sus valores de entrada son los
mismos que se describieron anteriormente para Secado en Flujo
Paralelo.EES; su nico valor de salida es la longitud del secador en
metros. Esta funcin debe escribirse a partir de la primera lnea en
la pantalla principal de EES, detallando entre parntesis todas las
variables de entrada necesarias para su ejecucin.
La funcin consta bsicamente de dos lazos repetitivos
anidados. El lazo principal se encarga de encerar el estado de la
partcula, es decir, llevarla a sus condiciones iniciales de cada, y
preguntar si sta ha alcanzado el contenido de humedad deseado.
El lazo sec ndario el c al est etido dentro del principal tiene
como funcin evaluar punto a punto el movimiento, transferencia de
masa y de calor, y propiedades del aire y de la partcula, cuando
esta ltima se encuentra en su trayectoria de vuelo. Este lazo
finaliza cuando la partcula se ha desplazado verticalmente una
cantidad igual a . Por otro lado, el lazo de repeticin principal
termina cuando la partcula ha alcanzado el nivel de humedad
requerido. En ese instante, la funcin es dada por
concluida, arrojando como resultado la longitud de secador que
resulta de la suma de los desplazamientos consecutivos de la
partcula en la direccin del eje longitudinal del tambor. El cdigo
que corresponde a la funcin se presenta a continuacin:

151
FUNCTION L_secador(X_1;X_2;T_s;rho_grano_seco;D_grano;m_DOT_p;
m_DOT_aire;T_2;phi_2;P1;R;alpha;Y_BAR;step)
Vol_grano=(4/3)*pi*(D_grano/2)^3 {Volumen del grano, m3}
A_s=4*pi*(D_grano/2)^2 {Superficie del grano, m2}
m_grano_seco=rho_grano_seco*Vol_grano {Masa del grano seco, kg}
m_grano=m_grano_seco/(1-X_1) {Masa inicial del grano hmedo, kg}
cont_agua=m_grano-m_grano_seco {Contenido inicial de agua, kg}
N_DOT_partculas=m_DOT_p/m_grano {Flujo de partculas, UND/S}
g=9,81 {Aceleracin de la gravedad, m/s2}
A_trans=pi*R^2 {rea transversal del tambor, m2}
h_2=ENTHALPY(AirH2O;T=T_2;P=P1;R=phi_2)
w_2=HUMRAT(AirH2O;T=T_2;P=P1;R=phi_2)
rho_2=DENSITY(AirH2O;T=T_2;P=P1;R=phi_2)
V_DOT_2=m_DOT_aire/rho_2 {Flujo volumtrico de aire, m3/s}
L_secador=0 {Variable acumuladora "Longitud de secador", m}
i=1
REPEAT
t_1=0
t2=step
V_z_0=0 {Velocidad inicial de la partcula en Z}
V_y_0=0 {Velocidad inicial de la partcula en Y}
h_BAR=0 {Coef. inicial de T.C., W/m2-K}
h_BAR_m=0 {Coef. inicial de T.M., m/s}
z_0=0
y_0=0
lookup('Lookup 1';i;1)=L_secador
lookup('Lookup 1';i;2)=Y_ac
lookup('Lookup 1';i;3)=X_1
lookup('Lookup 1';i;4)=T_s
lookup('Lookup 1';i;5)=phi_2
lookup('Lookup 1';i;6)=T_2
lookup('Lookup 1';i;7)=V_DOT_2
lookup('Lookup 1';i;8)=h_BAR
lookup('Lookup 1';i;9)=h_BAR_m
i=i+1
REPEAT
CALL SECADO(phi_2;w_2;h_2;T_2;D_grano;m_grano;
m_DOT_aire; N_DOT_partculas;R;alpha;P1;g;V_y_0;
V_z_0;z_0;y_0;A_trans;A_s;t_1;t2 : m_DOT_w;m_w;T_s;
h_BAR_m;h_BAR;V_z;V_y;z;y)
152
w_2=w_2+m_DOT_w/m_DOT_aire
T_2=TEMPERATURE(AirH2O;h=h_2;P=P1;w=w_2)
phi_2=RELHUM(AirH2O;T=T_2;P=P1;w=w_2)
rho_aire=DENSITY(AirH2O;T=T_aire;P=P1;w=w_aire)
V_DOT_2=m_DOT_aire/rho_2
cont_agua=cont_agua-m_w
m_grano=m_grano-m_w
X_1=cont_agua/m_grano
L_secador=L_secador+(z-z_0)
Y_ac=Y_ac-(y-y_0)
lookup('Lookup 1';i;1)=L_secador
lookup('Lookup 1';i;2)=Y_ac
lookup('Lookup 1';i;3)=X_1
lookup('Lookup 1';i;4)=T_s
lookup('Lookup 1';i;5)=phi_2
lookup('Lookup 1';i;6)=T_2
lookup('Lookup 1';i;7)=V_DOT_2
lookup('Lookup 1';i;8)=h_BAR
lookup('Lookup 1';i;9)=h_BAR_m
t_1=t2
t2=t2+step
V_y_0=V_y
V_z_0=V_z
z_0=z
y_0=y
i=i+1
UNTIL Y_ac<=0
UNTIL X_2>=X_1
END
Una vez declarada la funcin, y antes de entrar a los lazos de
repeticin, se ponen a punto todos aquellos parmetros necesarios
para empezar la ejecucin de la misma. En este caso, se desea
conocer la masa inicial del grano o partcula, su contenido inicial de
humedad, el rea superficial del grano (este parmetro se asume
constante durante el anlisis), flujo de partculas, la seccin

153
transversal del tambor, y las condiciones iniciales del aire. De all se
escriben las propiedades iniciales del aire y se inicializan las
variables , que es un acumulador de los desplazamientos de
la partcula en la direccin Z, e , que es un ndice o contador que
permitir llevar un registro y orden de cada paso consecutivo en el
tiempo.
El primer lazo de repeticin define el estado de la partcula cuando
sta se encuentra a punto de iniciar la cada, o el ciclo de secado.
Al hablar del estado inicial de la partcula, se hace referencia
principalmente a su posicin vertical, , velocidad inicial,
, y coeficientes de transferencia de calor y de masa,
. Aqu se incluyen todas las variables que deben ser
reseteadas al inicio de cada ciclo. No se incluye, por ejemplo, el
contenido de humedad de la partcula o su posicin horizontal
relativa al secador, dada por .
Las variables y definen el intervalo de tiempo en el que va a
tomar lugar el anlisis correspondiente al lazo repetitivo secundario.
Por otro lado, y son variables que servirn para calcular los
desplazamientos horizontal y vertical de la partcula en cada
intervalo de tiempo.
154
Para llevar un registro o historia de la variacin de parmetros, los
resultados arrojados por cada variable son almacenados en una
tabla de tipo LOOKUP. Esta tabla se llena automticamente durante
la ejecucin del programa principal. Al final, el usuario puede
escoger cualquier par de variables y graficarlas en EES para tener
una mejor visin de lo que est ocurriendo con el modelado del
secador. En este caso, las variables que se seleccionaron para
mostrar son las que se detallan en el cdigo anterior, donde el
ndice hace referencia al nmero de fila y el nmero ubicado
inmediatamente a la derecha especifica el nmero de columna.
Una vez definido el estado inicial de la partcula, se da comienzo al
ciclo de secado (o cada de la partcula) que toma lugar en el
segundo lazo de repeticin. En la primera lnea de este lazo se hace
un llamado o CALL al subprograma SECADO, el cual se presenta al
final de esta seccin. Este subprograma posee variables de
entrada, en azul, y variables de salida, en rojo. La ventaja de un
subprograma sobre una funcin, es que el primero permite trabajar
con mltiples variables de resultado; la funcin puede tener varias
entradas pero solo una variable de salida.
Para poder actualizar tanto las propiedades del aire como las de la
partcula, es necesario obtener del subprograma SECADO el valor

155
resultante de los siguientes parmetros: el flujo msico total de
agua evaporada entre y , , el agua evaporada por la
partcula entre y , , y la velocidad y coordenadas de la
partcula en Z y en Y en el instante , o sea, , , y .
El resto de variables de salida que aparecen en la declaracin del
subprograma ( , y ) pueden ser omitidas sin problema; no
tienen influencia alguna sobre la correccin de propiedades ya que
el subprograma las calcula automticamente en cada corrida. En
este caso particular, se han incluido por fines informativos y de
anlisis; solo se desea conocer cul es su variacin a lo largo del
secador. Si lo desea, el usuario puede dar seguimiento a la
evolucin de cualquier variable que se haya definido ya sea en el
programa principal, en la funcin, o en el subprograma; solo es
cuestin de incluirla en la tabla lookup.
Una vez obtenidos los valores de las variables de resultado
requeridas, se corrigen las propiedades psicromtricas del aire, el
contenido de agua y la masa total de la partcula rastreada. Todos
estos valores son almacenados en una nueva fila de la tabla lookup
para s posterior is alizacin. Al inal se prepara al lazo
repetitivo para continuar con el anlisis en el siguiente intervalo de
tiempo.
156
El procedimiento anterior se repite hasta que la partcula toca el
fondo del tambor, o sea, . En ese instante, se plantea la
pregunta de si la partcula ha alcanzado o no el nivel de humedad
deseado. De ser negativa la respuesta, se repite el ciclo de secado,
es decir, la partcula es inmediatamente recogida por una paleta y
llevada a su posicin inicial de cada.
Una vez que la partcula alcanza el nivel de humedad requerido, se
da por finalizado el lazo repetitivo principal y, en seguida, la funcin
arroja como resultado la suma final de todos los
desplazamientos en Z de la partcula. Otro resultado importante, por
supuesto, es la tabla lookup que muestra la evolucin de los
parmetros durante el proceso de secado.
El subprograma SECADO es el que se encarga de realizar el
trabajo s cio p esto q e es el q e res elve de forma simultnea
los problemas de movimiento, transferencia de calor y de masa. Es
un programa tpico de EES donde el nmero total de ecuaciones
debe coincidir con el nmero de incgnitas existentes. El
subprograma puede ser llamado, ya sea desde el programa
principal o desde la funcin, cuantas veces se requiera, mediante el
uso de la sentencia CALL. A continuacin se presenta su desarrollo.

157
SUBPROGRAM SECADO(phi_2;w_2;h_2;T_2;D_grano;m_grano;m_DOT_aire;
N_DOT_partculas;R;alpha;P1;g;V_y_0;V_z_0;z_0;y_0;A_trans;A_s;t_1;
t2 : m_DOT_w;m_w;T_s;h_BAR_m;h_BAR;V_z;V_y;z;y)
T_f=(T_2+T_s)/2 {Temperatura flmica}

C_D=32,628*abs(rho_B*(V_2-V_z)*D_grano/mu_B)^(-0,646)
C_L=32,628*abs(rho_B*V_y*D_grano/mu_B)^(-0,646)
{INTEGRACIN PARA DETERMINAR LA COMPONENTE "Z" DE LA

VELOCIDAD ABS DE LA PARTCULA en t=t2}
m_grano*dV_zdt=m_grano*g*sin(alpha)+0,5*rho_B*C_D*(V_2-
V_z)^2*(pi*(D_grano/2)^2)
V_z=V_z_0+integral(dV_zdt;Time;t_1;t2)
{INTEGRACIN PARA DETERMINAR LA COMPONENTE "Y" DE LA VELOCIDAD

ABS DE LA PARTCULA en t=t2}
m_grano*dV_ydt=m_grano*g*cos(alpha)-
0,5*rho_B*C_L*V_y^2*(pi*(D_grano/2)^2)
V_y=V_y_0+integral(dV_ydt;Time;t_1;t2)
{INTEGRACIN PARA DETERMINAR LA COORDENADA "Z" DE LA PARTCULA

en t=t2}
z=z_0+integral(V_z;Time;t_1;t2)
{INTEGRACIN PARA DETERMINAR LA COORDENADA "Y" DE LA PARTCULA

en t=t2}
y=y_0+integral(V_y;Time;t_1;t2)
{DETERMINACIN DE LA VELOCIDAD RELATIVA ENTRE EL AIRE Y LA

PARTCULA}
V_DOT_2=m_DOT_aire/rho_2 {Flujo volumtrico de aire}
V_2=V_DOT_2/A_trans {Velocidad del aire}
V_rel=SQRT((V_2-V_z)^2+V_y^2) {Velocidad relativa entre el
grano y el aire}
{Balance de energa}
h_BAR*(T_2-T_s)=h_fg*h_BAR_m*(rho_A_sat-rho_A_INFINITY)
h_f=ENTHALPY(Water;T=T_s;X=0)*CONVERT(KJ/kg;J/kg)
h_g=ENTHALPY(Water;T=T_s;X=1)*CONVERT(KJ/kg;J/kg)
158
h_fg=h_g-h_f {Calor latente de vaporizacin @ T_s}

rho_A_sat=DENSITY(Water;T=T_s;X=1)
rho_A_infinity=phi_2*DENSITY(Water;T=T_2;X=1)
{Evaluacin de NUSSELT mediante la correlacin de Whitaker.

Todas las propiedades, excepto mu_s se evalan a T_2.}
NUSSELT_BAR_D=2+(0,4*Re_D^(1/2)+0,06*Re_D^(2/3))*Pr_2^0,4*
(mu_2/mu_2_s)^(1/4) {Correlacin de Whitaker}
NUSSELT_BAR_D=h_BAR*D_grano/k_2 {Determinacin del coef. de T.C}
Re_D=(rho_2*V_rel*D_grano)/mu_2 {Nmero de Reynolds}
rho_2=DENSITY(AirH2O;T=T_2;P=P1;w=w_2)
mu_2=VISCOSITY(AirH2O;T=T_2;P=P1;w=w_2)
mu_2_s=VISCOSITY(AirH2O;T=T_s;P=P1;w=w_2)
Pr_2=PRANDTL(Air;T=T_2)
k_2=CONDUCTIVITY(AirH2O;T=T_2;P=P1;w=w_2)
{Analoga de Reynolds. Las propiedades de la especie B se

evalan a la temperatura flmica.}
h_BAR_m=(h_BAR*D_AB*abs(Le)^(1/3))/k_B
rho_B=DENSITY(Air;T=T_f;P=P1)
mu_B=VISCOSITY(Air;T=T_f)
k_B=CONDUCTIVITY(Air;T=T_f)
cp_B=CP(Air;T=T_f)*1000
Le=k_B/(rho_B*cp_B*D_AB)
D_AB=-2,775E-6[m^2/s]+4,479E-8[m^2/s-K]*(T_f+273,15)+1,656E-
10[m^2/s-K^2]*(T_f+273,15)^2 {Coeficiente de Difusin Aire-
Vapor de Agua}
n_A=h_BAR_m*A_s*(rho_A_sat-rho_A_INFINITY)
{Transferencia de masa para una sola partcula, kg/s}
{Agua Evaporada}
m_w=n_A*(t2-t_1) {Agua que pierde una partcula entre
t_1 y t2, kg}
m_DOT_w=m_w*N_DOT_partculas {Flujo total de vapor de agua
entregado al aire entre t_1 y t2, kg/s}
END
159
La estructura del subprograma SECADO parece estar bastante
clara. Se han utilizado comentarios entre llaves para explicar las
distintas partes del cdigo. Primero se determina la temperatura
flmica y los coeficientes de arrastre en ambas direcciones, Z y Y.
En este caso, se han utilizado las relaciones determinadas a partir
de la figura 4.14. No obstante, se recomienda realizar pruebas
experimentales para cada material que se vaya a utilizar. Las
irregularidades en la forma de la partcula pueden influir sobre los
coeficientes de arrastre, sobretodo cuando se ha trabajado bajo una
suposicin de partcula esfrica.
Posteriormente, se hace una doble integracin para determinar la
velocidad y posicin de la partcula en el tiempo . Para la siguiente
corrida, estos resultados se convierten en las condiciones iniciales
de movimiento de la partcula, es decir, , , y .
Con el problema de movimiento resuelto, se plantea la ecuacin de
balance de energa de acuerdo al volumen de control presentado en
la figura 5.2. Conjuntamente, se calculan las propiedades
necesarias para resolver esta ecuacin: calor latente de
vaporizacin, , y densidades del vapor de agua inmediatamente
sobre la interfaz gas-lquido y en el aire de secado, y
respectivamente.
160
En seguida, se determina el nmero de Nusselt a partir de la
correlacin de Whitaker y se calcula el coeficiente promedio de
transferencia de calor. Haciendo uso de la analoga de transferencia
de calor y de masa, se estima el coeficiente promedio de
transferencia de masa.
Finalmente, se calcula la cantidad de agua evaporada por la
partcula entre y y, en base a este resultado, se estima el flujo
total de agua evaporada, por todas las partculas que atraviesan la
seccin del tambor, para el mismo intervalo de tiempo.
6.6 Secado en Contraflujo.
Al igual que en el caso anterior, la partcula ser rastreada durante
su recorrido por el secador desde su estado de humedad ms alto
hasta su estado de humedad ms majo. Sin embargo, la
inicializacin de las variables correspondientes al aire se realizar
en torno a sus propiedades posteriores al secado, es decir, con su
humedad relativa ms alta. A medida que transcurra el programa de
secado en contraflujo, tanto el aire como la partcula irn perdiendo
humedad en la misma direccin.

161
En el programa principal, llamado Secado en Contraflujo.EES, en
lugar de introducir la temperatura y humedad relativa del aire
despus del calentamiento, y respectivamente, se ingresa
y , que definen su estado a la finalizacin del proceso de secado
(referirse a la figura 3.3). Estos parmetros se obtienen tambin del
programa Anlisis Termodinmico.EES.
Una vez declaradas las variables de entrada, se ejecuta la funcin
, la cual posee una estructura similar a la que aparece en el
programa de secado en flujo paralelo, con la diferencia de que en el
lazo repetitivo secundario, en lugar de escribir ,
se escribe .
As como la humedad absoluta, todas las variables que se refieran a
las propiedades del aire, tales como densidad, humedad relativa,
te perat ra y entalp a se escriben con el s b ndice 3; esto es
solo para tener presente que el potencial de secado del aire va de
menos a ms durante la ejecucin del programa.
Finalmente, las ecuaciones referentes al movimiento de la partcula,
descritas en la seccin 4.4.2, se escriben en el subprograma
SECADO de la siguiente manera:
{COEFICIENTES DE ARRASTRE EN Z y Y}
C_D=32,628*abs(rho_B*(V_3+V_z)*D_grano/mu_B)^(-0,646)
C_L=32,628*abs(rho_B*V_y*D_grano/mu_B)^(-0,646)
162
{INTEGRACIN PARA DETERMINAR LA COMPONENTE "Z" DE LA VELOCIDAD ABS

DE LA PARTCULA en t=t_2}
m_grano*dV_zdt=m_grano*g*sin(alpha)-
0,5*rho_B*C_D*(V_3+V_z)^2*(pi*(D_grano/2)^2)
V_z=V_z_0+integral(dV_zdt;Time;t_1;t_2)
{INTEGRACIN PARA DETERMINAR LA COMPONENTE "Y" DE LA VELOCIDAD ABS

DE LA PARTCULA en t=t_2}
m_grano*dV_ydt=m_grano*g*cos(alpha)-
0,5*rho_B*C_L*V_y^2*(pi*(D_grano/2)^2)
V_y=V_y_0+integral(dV_ydt;Time;t_1;t_2)
{INTEGRACIN PARA DETERMINAR LA COORDENADA "Z" DE LA PARTCULA en

t=t_2}
z=z_0+integral(V_z;Time;t_1;t_2)
{INTEGRACIN PARA DETERMINAR LA COORDENADA "Y" DE LA PARTCULA en

t=t_2}
y=y_0+integral(V_y;Time;t_1;t_2)
{DETERMINACIN DE LA VELOCIDAD RELATIVA ENTRE EL AIRE Y LA

PARTCULA}
V_DOT_3=m_DOT_aire/rho_3
V_3=V_DOT_3/A_trans
V_rel=SQRT((V_3+V_z)^2+V_y^2)
Notar que la nica diferencia con respecto al subprograma de
secado en flujo paralelo est en los signos de la sumatoria de
fuerzas en Z y en el valor de la velocidad relativa entre el aire y la
partcula. De all en adelante, todo se mantiene igual.

163
CAPTULO 7
7. RESULTADOS Y ANLISIS
A continuacin se presenta informacin referente a los casos de estudio
de dos secadores rotativos previamente analizados con modelos
desarrollados en otros trabajos. Los datos mostrados sern utilizados
para comparaciones posteriores con el modelo terico propuesto en este
estudio.
7.1 Parmetros de Operacin y Dimensiones de Secadores Tipo.
Las tablas 1 y 2 son una sntesis de las caractersticas del producto,
geometra y parmetros de operacin de los secadores estudiados
por Kamke y Ajayi para el procesamiento de partculas de madera y
concentrados de zinc, respectivamente. Las metodologas y

164
modelos planteados por estos investigadores fueron brevemente
descritos en el apartado 2.3.
Los datos que se mostrarn a continuacin fueron extrados de las
publicaciones los autores mencionados anteriormente (11) (12); en
algunos casos fue necesario deducir ciertos parmetros a partir de
la informacin disponible. Luego de las tablas 1 y 2, se presenta un
resumen de los procedimientos experimentales y resultados
obtenidos en cada trabajo, que servirn despus para establecer
comparaciones con los resultados logrados en EES.

165
TABLA 1
DATOS DEL SECADOR ANALIZADO CON EL MODELO A
Partculas de
madera
DATOS DEL
PRODUCTO
Material a secar Douglas-Fir
Densidad del grano seco (kg/m3) 450
Densidad al granel promedio (kg/m3) 200
ngulo dinmico de reposo promedio () 82,6
Dimetro promedio de partcula (mm) 1,81
Flujo de aire con respecto al flujo de producto Paralelo
Dimetro del tambor (m) 1,2
DATOS DEL
SECADOR
Longitud del tambor (m) 5,5

Inclinacin () 0
Longitud segmento 1 de paleta (m) 0,208
ngulo entre segmentos 1 y 2 () 90
Temperatura ambiente (C) 17,5
Humedad relativa ambiental 0,35
PARMETROS DE
Presin baromtrica (atm) 1

OPERACIN
Velocidad de giro (r.p.m.) 5,5

Flujo de producto seco (TON slido seco / hr) 0,283
Contenido de humedad inicial (% BASE HMEDA) 58,39
Contenido de humedad final (% BASE HMEDA) 39,28
Temperatura del aire antes del secado (C) 155,14
Temperatura del aire despus del secado (C) 73,75
Flujo volumtrico de aire a la salida (cfm) 4200
166
TABLA 2
DATOS DEL SECADOR ANALIZADO CON EL MODELO B

Concentrados
Material a secar de zinc
Densidad del grano seco (kg/m3) 4150
1595
DATOS DEL
PRODUCTO
Densidad al granel promedio (kg/m3)

ngulo dinmico de reposo promedio REGIN B () 55,96
ngulo dinmico de reposo promedio REGIN C () 52,84
ngulo dinmico de reposo promedio REGIN D () 49,15
Dimetro promedio de partcula REGIN B (mm) 15
Dimetro promedio de partcula REGIN C (mm) 12
Dimetro promedio de partcula REGIN D (mm) 8
Flujo de aire con respecto al flujo de producto Paralelo
Dimetro del tambor (m) 3,9
Inclinacin () 4
2,4
REGIN B
Longitud (m)
DATOS DEL SECADOR

3,3
REGIN C
Longitud (m)
6,6
REGIN D
Longitud (m)
Temperatura ambiente (C) 28
Humedad relativa ambiental 0,7
PARMETROS DE
Presin baromtrica (atm) 1

OPERACIN
Velocidad de giro del tambor (r.p.m.) 3

Flujo de producto seco (TON slido seco / hr) 122,2
Humedad inicial del producto (% BASE HMEDA) 16,30
Humedad final del producto (% BASE HMEDA) 12,40
Temperatura del aire antes del secado (C) 500
Temperatura del aire despus del secado (C) 131
167
El comportamiento de secado de partculas de madera Douglas-Fir
(figura 7.1) se examin mediante un sistema similar al que se
muestra en la figura 7.2. El objetivo del experimento consista en
obtener los perfiles de humedad y temperatura para ambos, el gas y
las partculas, en toda la longitud del secador. En este caso, las
variables de control estudiadas fueron temperatura del aire en la
caja de mezcla, flujo volumtrico de aire a la salida del secador y
velocidad angular del tambor. Los rangos empleados para cada
variable fueron, respectivamente, 550 a 750 C, 3180 a 4240 cfm, y
2.8 a 5.5 rpm. El rango de temperaturas correspondiente a la caja
de mezcla genera un rango de temperaturas de 150 a 270 C en la
entrada del secador (punto I en la figura 7.2).
FIGURA 7.1 PARTCULAS DE MADERA

DOUGLAS-FIR TOMADO DE (23)
168
FIGURA 7.2 SECADOR DE PARTCULAS

DE MADERA TOMADO DE (11)
En las pruebas experimentales, se utiliz una combinacin de
paletas perifricas y paletas internas o de llenado central, tal como
se muestra en la parte derecha de la figura 7.3. Para la validacin
del modelo, se compararon los resultados experimentales con
aquellos simulados en el software desarrollado por Kamke. Las
mediciones de humedad y temperatura se realizaron en los puntos
B, C, D, E y F (referirse a la figura 7.2).

169
FIGURA 7.3 POSIBLES CONFIGURACIONES DE PALETAS
PARA EL MODELO DE SECADO A TOMADO DE (11)
Los resultados obtenidos de la comparacin entre los datos
experimentales y tericos, para los parmetros de operacin de la
tabla 1, se presentan en la figura 7.4.
Para mostrar las aplicaciones de su modelo, Kamke ejecut una
simulacin para el mismo secador, operando bajo las mismas
condiciones, pero utilizando solamente paletas perifricas.
Entonces, en la figura 7.5, l compara los resultados entre las
simulaciones de un secador con paletas internas y uno sin ellas.
Las curvas de humedad y temperatura obtenidas de la simulacin
para el secador sin paletas internas son las que se utilizarn en este
trabajo para establecer comparaciones contra el modelo
lagrangiano.
170
FIGURA 7.4 COMPARACIN ENTRE LOS RESULTADOS

EXPERIMENTALES Y SIMULADOS PARA UN SECADOR CON
PALETAS DE RELLENO CENTRAL TOMADO DE (11)
FIGURA 7.5 COMPARACIN ENTRE LAS SIMULACIONES

PARA UN SECADOR CON Y SIN PALETAS DE RELLENO
CENTRAL TOMADO DE (11)
171
La figura 7.6 muestra el proceso esquemtico de secado de zinc en
la mina MMG (Queensland, Australia). El producto es suministrado
al secador mediante un tornillo sinfn, con una humedad inicial que
vara entre 16 y 18%, alcanzando a una humedad mnima del 12%
a la salida del secador. Luego de haber pasado por la cmara de
combustin, el aire entra al secador con una temperatura de 500C.
FIGURA 7.6 REPRESENTACIN ESQUEMTICA DE LA

CMARA DE COMBUSTIN Y DEL SECADOR ROTATIVO
TOMADO DE (12)
El secador consta de cinco regiones, A, B C, D y E. La primera y
ltima regin no poseen paletas elevadoras; el avance del producto
est dado nicamente por la rotacin e inclinacin del tambor. Por
otro lado, las regiones B, C y D poseen paletas perifricas de dos

172
segmentos, cuyas dimensiones se detallan en la tabla 2. Solo estas
tres regiones sern objeto de anlisis en el presente trabajo.
La distribucin de tamaos de partcula de los concentrados de zinc
de determin mediante el cribado de muestras tomadas en
diferentes secciones del secador. En total se tomaron 23 muestras,
una por cada metro de tambor. Los dimetros promedio de partcula
para las secciones A, B, C, D y E, fueron 20, 15, 12, 8 y 7 mm,
respectivamente. La disminucin en dimetro se debe al
desmoronamiento de las partculas, mas no al encogimiento de las
mismas.
Ajayi (12), quien model y estudi el desempeo de este secador
rotativo, obtuvo el perfil de humedad experimental, , a lo largo del
tambor y tambin determin el valor del ngulo dinmico de reposo
para los concentrados de zinc, , a diferentes niveles de humedad.
Las relaciones experimentales encontradas fueron
(7.1)
(7.2)
En las ecuaciones anteriores, y .

173
Las variables de control estudiadas en este caso fueron flujo msico
de producto y velocidad angular del tambor, en rangos de 97 a 135
tonslido_seco/hr y 2 a 3 rpm, respectivamente. Los resultados de las
comparaciones entre la informacin experimental y los datos
obtenidos mediante el modelo de Ajayi, para los datos de la tabla 2,
se presentan en la figura 7.7.
FIGURA 7.7 RESULTADOS TERICOS Y EXPERIMENTALES
PARA EL CONTENIDO DE HUMEDAD EN LOS SLIDOS
TOMADO DE (12)
Para determinar el rea superficial de las partculas desnudas en
contacto con la corriente de aire, Ajayi calcul el nmero de

174
partculas slidas e introdujo un factor de correccin de rea (FCA)
para incluir la incertidumbre en el rea superficial estimada.
Tambin asumi que la mitad del rea superficial de las partculas
estaba en contacto con el aire. En la figura 7.7 se muestran los
resultados de la simulacin sin factor de correccin de rea y con
un factor correccin de rea constante. Se puede observar de la
grfica que ambos enfoques resultaron en un ajuste pobre del perfil
de humedad experimental en los slidos.
Para mejorar el ajuste, Ajayi utiliz diferentes factores de correccin
para cada seccin del tambor. Los factores de correccin para las
secciones B, C y D, fueron 1.2, 0.6 y 0.42 respectivamente. Los
resultados de la simulacin corregida se muestran en la figura 7.8.
FIGURA 7.8 AJUSTE DEL PERFIL DE HUMEDAD EN LOS

SLIDOS CON FCA VARIABLE TOMADO DE (12)
175
Los perfiles reales de temperatura del aire y de los slidos a travs
del secador no pudieron ser obtenidos debido al peligro que
involucraba ejecutar ensayos experimentales en ese entorno
industrial. Como consecuencia, en el estudio de Ajayi se asumi
que si los valores de temperatura estimados a la salida del secador
concordaban con los valores experimentales, entonces los perfiles
de temperatura simulados se consideraran apropiados. Los
resultados de estas simulaciones se presentan en las figuras 7.9 y
7.10.
FIGURA 7.9 SIMULACIN DEL PERFIL DE TEMPERATURA DE

LOS SLIDOS CON FCA VARIABLE TOMADO DE (12)
176
FIGURA 7.10 SIMULACIN DEL PERFIL DE TEMPERATURA

DEL AIRE CON FCA VARIABLE TOMADO DE (12)
En el presente estudio, se combinaron las ecuaciones 7.1 y 7.2 para
determinar el ngulo dinmico de reposo promedio correspondiente
a las regiones B, C y D. Con esta informacin, y tomando en cuenta
la geometra de las paletas, se calcul la altura promedio de cada
en cada regin.
Para estimar las condiciones del aire de secado a la entrada de la
regin B, primero se realiz un balance de masa y energa global
del secador para hallar el flujo msico de aire requerido. Para ello,
se ejecut el programa Anlisis termodinmico.EES considerando
las propiedades de entrada y salida tanto del aire como de los

177
concentrados de zinc. Luego, mediante la ecuacin 7.1, se
determin el contenido de humedad a la entrada de la regin B
( ) y se calcul el porcentaje de agua evaporada en la
regin A. Finalmente, mediante un anlisis termodinmico de la
regin A solamente, se estimaron las propiedades del aire a la
salida de la regin A (entrada regin B).
Una vez obtenida la informacin necesaria, se ejecut el programa
Secado en flujo paralelo.EES para cada una de las regiones B, C y
D, tomando en cuenta el dimetro promedio de partcula, la altura
de cada y las condiciones del aire al comienzo de cada seccin.
Los resultados de las simulaciones en EES para el secador de
concentrados de zinc, as como para el secador de partculas de
madera, se presentan en el siguiente apartado.
7.2 Resultados Obtenidos en EES.
A partir de las caractersticas y parmetros de operacin de los
secadores presentados en las tablas 1 y 2, se ejecutaron los
programas elaborados en este estudio para obtener los perfiles de
humedad y temperatura de los slidos y del aire de secado.

178
Se elaboraron tablas de resumen con algunos parmetros
importantes que pueden dar una idea general del desempeo del
secador. Entre estos parmetros se tienen el flujo msico de aire,
tiempo de secado, nmero de ciclos de secado, ngulo y altura
promedio de cada, tiempo de residencia de la partcula,
temperatura de la partcula y propiedades del aire a la entrada y
salida del secador.
El tiempo de secado se puede estimar multiplicando el nmero de
ciclos de secado por el tiempo de vuelo de la partcula. En EES, el
nmero de ciclos se puede obtener ya sea del programa principal
de secado, agregando un contador en el lazo de repeticin
principal, o simplemente cuantificando manualmente los picos en la
grfica de trayectoria de la partcula.
El mtodo anterior tal vez introduzca un pequeo error debido a que
el tiempo de vuelo no necesariamente es el mismo en todos los
ciclos. Una manera de evitar aquello, es dirigirse a la tabla lookup y
observar el ndice de la ltima fila registrada durante la ejecucin
del programa. En este caso, el tiempo de secado sera igual al
nmero total de filas multiplicado por el paso temporal.
El tiempo de residencia, que es el tiempo que le toma a la partcula
viajar de un extremo del secador al otro, se estima a partir de la

179
ecuacin 4.9. El tiempo de secado es solo una pequea fraccin del
tiempo de residencia. Hay que tener presente que, cuando se habla
del tiempo de secado, nmero de ciclos, altura de cada y tiempo de
residencia, se hace referencia a un promedio que corresponde a la
partcula representativa que se est analizando.
Cuando se desea verificar la validez del modelo de movimiento en
el secador, es comn comparar las distribuciones de tiempo de
residencia (DTR) generadas por el modelo terico con los
resultados experimentales obtenidos a partir de pruebas con
trazadores. Esta ltima consiste en inyectar una cantidad conocida
de qumico inerte a la entrada del secador mientras su
concentracin es monitoreada a la salida como una funcin del
tiempo.
En este estudio, no obstante, se ha idealizado la trayectoria de una
partcula con el objeto de predecir la evolucin de los parmetros de
secado durante su recorrido. Para ello, se han planteado una serie
de suposiciones gruesas tales como que no existe interaccin entre
partculas, o que el choque entre partcula y tambor es
perfectamente inelstico. Es ms, se ha supuesto que el secador
est idealmente cargado; es decir, que la cantidad del producto que
est ingresando en el mismo es la cantidad que garantiza el llenado

180
perfecto de las paletas en la posicin . Por consiguiente, no se
espera que el tiempo de residencia calculado con el modelo
lagrangiano sea, necesariamente, un buen estimador de la media
en la DTR que se obtendra mediante una prueba de trazadores.
Todos los cdigos, ecuaciones, tablas y resultados brutos obtenidos
de los diferentes programas elaborados en EES se adjuntan en la
seccin de Apndices. All se incluyen anlisis termodinmicos,
clculos de ngulo representativo y altura promedio de cada, y
resultados de los programas de secado. En el caso del secador de
concentrados de zinc, solo se mostrar una de las tres regiones
analizadas a manera de ejemplo. Las tablas de tipo lookup que
aparecen en Apndices fueron elaboradas con un paso temporal de
0,05 segundos. No obstante, los resultados que aparecen en esta
seccin fueron obtenidos con un paso temporal de 0,01 segundos.

181
TABLA 3
RESUMEN DE RESULTADOS PARA
MADERA
( ) [ ] 0,1891
[ ] 58,39
[ ] 41,82
[ ] 1,892
[ ] 20,36
44
[ ] 47,56
[ ] 0,8283
[ ] 0,576
[ ] 3,01
Entrada 5,67
38,06
Salida 38,83
Entrada 0,001284
Salida 0,1021
Entrada 155,1
[ ]
Salida 81,2
Entrada 4898
[ ]
Salida 4235
182
FIGURA 7.11 TRAYECTORIA DE LA PARTCULA
FIGURA 7.12 TRAYECTORIA DE LA PARTCULA

(ACERCAMIENTO)
Ntese como algunos picos inferiores en la trayectoria de la
partcula atraviesan al eje Z generando un desplazamiento negativo.
Esto es un error que se introduce debido a que la discretizacin en

183
los clculos fue hecha en base al tiempo de cada, con pasos
temporales fijos, y no a la altura promedio de cada. Una manera de
corregir esto, si es que es meritorio, es con un paso ms pequeo.
FIGURA 7.13 PERFILES DE HUMEDAD (MADERA)
FIGURA 7.14 PERFILES DE TEMPERATURA (MADERA)

184
TABLA 4
RESUMEN DE RESULTADOS PARA ZINC
Regin B C D
( ) [ ] 40,56 39,91 39,50
[ ] 15,56 14,88 14,05
[ ] 14,88 14,05 12,53
[ ] 12,43 12,43 12,43
[ ] 9,36 12,41 12,43
12 17 31
[ ] 48,93 38,82 50,77
[ ] 2,877 2,501 2,945
[ ] 0,3142 0,3142 0,3142
[ ] 1,24 1,43 3,12
Entrada 58,90 59,03 59,28
Salida 59,03 59,28 61,27
Entrada 0,000312 0,0005522 0,00242
Salida 0,0005522 0,002420 0,05857
Entrada 423,9 360,9 274,9

[ ]
Salida 360,9 274,9 139,0
Entrada 55794 52729 47646

[ ]
Salida 52729 47646 38586
185
FIGURA 7.15 TRAYECTORIA DE LA PARTCULA (REGIN B)
FIGURA 7.16 TRAYECTORIA DE LA PARTCULA (REGIN C)

186
FIGURA 7.17 TRAYECTORIA DE LA PARTCULA (REGIN D)
B
C
D
FIGURA 7.18 PERFILES DE HUMEDAD (ZINC)

187
FIGURA 7.19 PERFILES DE TEMPERATURA (ZINC)
En las curvas de humedad de la partcula y temperatura del aire
(figuras 7.18 y 7.19) se pueden observar unos quiebres en las
transiciones de la regin B a la C y de la C a la D. El mismo
fenmeno se observa en la figura 7.7, para el caso del perfil de
humedad de los slidos obtenido a partir del modelo B sin factor de
correccin de rea.
La causa que provoca este comportamiento es que, en el
modelado, se consider un cambio brusco en dimetro de partcula
de regin a regin.
188
7.3 Comparaciones.
Una vez simulados en EES los procesos de secado de partculas de
madera y concentrados de zinc, se compararon los perfiles de
humedad y temperatura obtenidos mediante el modelo lagrangiano,
contra aquellos generados por los modelos A y B. En general, se
puede verificar que existe una buena respuesta del modelo
propuesto en este estudio.
De acuerdo a la figura 7.20, la curva de humedad en los slidos
obtenida del modelo lagrangiano presenta un comportamiento
similar a la del modelo A; pareciera que existe proporcionalidad
entre ellas. A la salida del secador, se tiene una diferencia relativa
de 6.36% entre ambas estimaciones. En este caso particular, el
modelo lagrangiano es ms conservador que el modelo A en cuanto
a la prediccin del contenido de humedad del producto en cualquier
seccin del secador.

189
FIGURA 7.20 HUMEDAD DE PARTCULA VS POSICIN AXIAL
FIGURA 7.21 TEMPERATURA DEL AIRE VS POSICIN AXIAL
Una situacin parecida se tiene en la figura 7.21, donde la
temperatura del aire estimada por el modelo lagrangiano es
ligeramente mayor que la que predice el modelo A para cualquier
valor de Z. Esto es congruente con las grficas de humedad

190
presentadas anteriormente; el secado es ms efectivo en el modelo
A y por ende el aire pierde ms energa. La diferencia relativa entre
ambos valores tericos de temperatura es de 9.94% a la salida del
secador.
FIGURA 7.22 TEMPERATURA DE PARTCULA VS POSICIN
Kamke valid su modelo mediante pruebas experimentales en un
secador con paletas internas. Luego simul el comportamiento del
mismo secador, operado bajo condiciones similares, utilizando
solamente paletas perifricas. Los resultados, presentados en la
figura 7.5, muestran que la temperatura de los slidos no es
afectada por la presencia de paletas internas. Considerando
aquello, y de acuerdo a lo mostrado en la figura 7.22, los datos

191
experimentales de temperatura convergen hacia el perfil dibujado
por el modelo lagrangiano en la direccin positiva de Z.
La figura 7.23 presenta los datos tericos y experimentales sobre el
contenido de humedad en los concentrados de zinc. A la salida del
secador, se tiene un error relativo de 4.62% para el modelo
lagrangiano y de 6.44% para el modelo B. Si se considera toda la
longitud del secador, se puede notar que el modelo lagrangiano se
apega ms hacia los datos experimentales. El error se hace ms
evidente en la ltima regin del secador.
FIGURA 7.23 HUMEDAD DE PARTCULA VS POSICIN AXIAL

192
Para ajustar su modelo de evaporacin, Ajayi utiliz distintos
factores de correccin de rea (FCA) para cada regin del secador.
Luego de obtener resultados satisfactorios, procedi a simular los
perfiles de temperatura del aire y del producto, sin corroborarlos
experimentalmente. Estos perfiles se presentan a continuacin,
junto a los que se obtuvieron del modelo lagrangiano.
FIGURA 7.24 TEMPERATURA DEL AIRE VS POSICIN

193
Suponiendo que la curva del modelo B sea correcta, el modelo
lagrangiano predice bien la temperatura del aire dentro de las dos
primeras regiones de secador analizadas. En la ltima etapa de
secado se genera una divergencia pronunciada debido al cambio
brusco en el tamao promedio de partcula.
Esta disminucin de dimetro no es producida por el encogimiento
de las partculas, sino ms bien por su rompimiento debido al
impacto; es como si la curva de distribucin de tamaos se sesgara
a la izquierda. Los concentrados de sulfuro de zinc son poco
resistentes al impacto.
194
Una situacin muy particular ocurre con el perfil de temperatura de
los slidos. La brecha existente entre el modelo B y el modelo
lagrangiano es relativamente grande; la diferencia es de unos 18C.
Lamentablemente, no existe informacin experimental que
dictamine cul de los dos modelos se apega ms a la situacin real.
Sin embargo, s se puede asegurar que el cambio en el tamao de
partcula tiene poca o ninguna influencia sobre la evolucin de
temperatura superficial en el producto. Para demostrar lo anterior,
se ejecut nuevamente el programa de secado considerando esta
vez un dimetro de partcula constante e igual a 15 mm; ste es el
dimetro promedio que corresponde a la regin B. Los resultados
obtenidos se presentan a continuacin.
FIGURA 7.26 HUMEDAD DE PARTCULA VS POSICIN

(DIMETRO DE PARTCULA CONSTANTE)
195
FIGURA 7.27 TEMPERATURA DEL AIRE VS POSICIN


196
Con la utilizacin de un dimetro constante de partcula, el perfil de
humedad de los slidos mejor notablemente en relacin a los
datos experimentales. Lo mismo se puede decir del perfil de
temperatura del aire; ya no existen quiebres en la curva y la
diferencia con respecto al modelo B es aproximadamente constante
a lo largo de todo el secador.
Por otro lado, la curva de temperatura superficial de los slidos se
mantuvo prcticamente intacta a pesar de la modificacin
efectuada.
7.4 Anlisis sobre el Modelo Lagrangiano.
Los tres modelos despliegan caractersticas diferentes. El modelo A
toma en cuenta la naturaleza del comportamiento de secado de la
madera e introduce una ecuacin que relaciona su contenido de
humedad en el equilibrio con la humedad relativa del aire. Tambin
considera la transferencia de calor en el interior de los slidos y la
prdida de calor a travs de la pared del cilindro.
A pesar de estas consideraciones, y de acuerdo a la informacin
experimental disponible, el modelo lagrangiano predice de mejor

197
manera la temperatura de los slidos en el secador. No obstante,
existe una ligera variacin en cuanto a los perfiles de humedad del
producto y temperatura del aire predichos por ambos modelos. A
continuacin se analizan los factores que podran provocar estas
diferencias.
Si se consideraran las prdidas de calor, por ejemplo, el perfil de
temperatura del aire que resulta del modelo lagrangiano se
acercara ms hacia la curva del modelo A; disminuira su potencial
de secado. Sin embargo, se tendra un efecto contrario al deseado
en la curva de humedad de los slidos. Aun cuando el modelo A s
considera las prdidas de calor, ste predice un secado ms
eficiente.
Tambin se descarta la posibilidad de que los mecanismos de
transferencia de calor y de masa en el interior de los slidos tengan
mayor influencia sobre el secado bajo las condiciones de operacin
dadas. Para este caso particular, con una humedad inicial en la
madera de 58.39% y una humedad final de 39.28%, y considerando
que la humedad relativa mxima del aire predicha con el modelo
lagrangiano es de solamente 10%, es muy poco probable que se
alcance el contenido de humedad en el equilibrio en la madera;
referirse a la figura 1.2. Es decir, hay suficiente humedad libre como

198
para que la evaporacin se de manera continua. Es difcil que la
difusin de humedad en el interior de los slidos pudiese llegar a
controlar la tasa de evaporacin en su superficie en este caso.
El transporte de las partculas y el efecto que pueda llegar a tener el
arrastre del aire sobre su movimiento, tiene gran influencia sobre la
longitud final del secador; no tanto en los parmetros de secado.
Por otro lado, el rea efectiva de transferencia de calor y de masa s
es un factor que podra alterar la forma de las curvas de humedad y
temperatura de manera apreciable. En este caso, el rea est dada
a partir del dimetro medio en la distribucin de tamaos,
asumiendo la forma de partcula esfrica. Esta suposicin fue muy
ligera considerando la forma irregular que poseen las astillas de
madera (vase la figura 7.1). Tal vez el trabajar con partculas de
forma cilndrica, o utilizando un dimetro de partcula de rea
equivalente, arroje un rea efectiva mayor, lo que aumentara la
tasa total de evaporacin. Es probable que este factor de correccin
de rea sea el que compense la desproporcionalidad leve que
existe entre los perfiles de humedad y temperatura en las figuras
7.20 y 7.21.
Otro factor que puede influenciar el rea efectiva de transferencia
de calor es el nmero real de partculas que se encuentran en

199
contacto ntimo con el aire de secado durante la produccin de la
cascada. Aqu se ha supuesto que todas las partculas estn
completamente envueltas por el aire y que, adems, no existe
interaccin entre ellas. En el caso del secador de madera, existe
una buena dispersin en la seccin transversal del secador debido
a que el ngulo dinmico de reposo es elevado.
El secador analizado con el modelo B es un caso especial porque
tanto la altura promedio de cada como el tamao promedio de
partcula son diferentes en las distintas regiones del secador.
La altura promedio de caa afecta directamente al tiempo de
secado.
Por otro lado, la principal consecuencia del cambio de dimetro de
partcula es el aumento del flujo promedio de partculas y la
disminucin del rea total de transferencia de calor y de masa. Si el
dimetro promedio disminuye en cierta proporcin , el flujo de
partculas aumenta en mientras que el rea total disminuye en
. Se cree que esto provoc la divergencia pronunciada en la
ltima regin de la curva de temperatura mostrada en la figura 7.24.
Se podra decir abiertamente que el efecto neto del cambio de

200
dimetro es igual a , pero no es tan simple porque existen otros
parmetros, afectados tambin por el cambio de dimetro, que
influyen sobre la tasa de evaporacin. Tal es el caso del coeficiente
promedio de transferencia de calor.
Al igual que en el secador de partculas de madera, aqu tambin se
descarta la posibilidad de que los mecanismos de transferencia de
calor y de masa en el interior de los slidos tengan mayor influencia
sobre el secado, ya que la curva de humedad en el equilibrio para el
xido de zinc es casi vertical (ver figura 1.2).
Cuando el modelo B no utiliza ningn factor de correccin de rea,
el modelo lagrangiano predice con un poco ms de precisin el
contenido de humedad de los slidos. El error es significativamente
menor cuando se considera el mismo dimetro promedio de
partcula en las tres regiones analizadas.
La temperatura del aire empleada para el secado de zinc fue mucho
mayor que la que se utiliz para el secado de madera. Ajayi report
temperaturas del casco sin aislamiento cercanas a los 300C
durante la prueba experimental (12). Entonces, aqu las prdidas de
calor hacia el entorno s pueden llegar a ser significativas en el
desempeo del secador. Puede ser que ste sea el motivo de la

201
diferencia entre las temperaturas de los slidos estimadas por el
modelo B y el lagrangiano.
El modelo B, aparte de incluir las prdidas de calor, considera que
la energa es transferida hacia los slidos por conveccin y
radiacin desde el aire de secado. Sin embargo, no existe mayor
diferencia entre las curvas de humedad estimadas por dicho modelo
y el modelo aqu propuesto, el cual no considera la transferencia de
calor por radiacin.
Aspectos relacionados con extender los resultados de una
partcula hacia todo el flujo de partculas
En este trabajo, para calcular la tasa total de evaporacin se
multiplic el agua evaporada por una partcula en un tiempo dado,
por el flujo total de partculas, y en base a este resultado se
corrigieron las propiedades del aire de secado. El flujo total de
partculas, por otro lado, se determin dividiendo el flujo msico de
producto hmedo para la masa de una partcula de tamao igual al
tamao promedio de la muestra, antes de su ingreso al secador.
Todo esto fundamentado en la suposicin de que una partcula es
representativa de todas las dems.
Es posible que estas simplificaciones, necesarias para el desarrollo
del modelo, traigan consecuencias sobre la precisin en las

202
predicciones de los perfiles de humedad y temperatura, tanto del
producto como del aire, a lo largo del secador.
La transferencia de masa y calor, as como el arrastre del aire, se
estimaron bajo la suposicin de que la partcula no es influenciada
por la presencia de las dems. La validez de esta aproximacin
depender de cun bien dispersos estn los slidos en la seccin
transversal del tambor; aqu interviene la calidad en el diseo de las
paletas elevadoras y la velocidad del aire. La dispersin de las
partculas no solamente influye sobre el coeficiente de transferencia
de calor y el coeficiente de arrastre, sino tambin en la cantidad real
de partculas desnudas en contacto con el aire.
Cuando un secador rotativo se encuentra sobrecargado (flujo
msico de producto mayor al flujo de diseo), es posible que cierta
cantidad de material se desplace a travs del mismo por efecto
kilning: el producto que no es recogido por la paleta se balancea en
la parte inferior del tambor y se desplaza hacia adelante debido a la
rotacin e inclinacin del mismo. Entonces, esto puede dar lugar a
que el flujo de partculas est un poco sobre estimado con respecto
a aquel que est siempre en contacto con el aire.
Los choques entre partculas, as como el arrastre del aire sobre la
cortina de slidos, afecta el desplazamiento de las mismas. Por lo

203
tanto, es posible que exista cierta diferencia entre el desplazamiento
estimado de la partcula en el modelo lagrangiano y aquel
desplazamiento medio de la cortina de slidos. Asimismo, el tiempo
de residencia estimado por el modelo lagrangiano tal vez difiera de
la media en la DTR del secador debido a los factores antes
mencionados.
Una modificacin a la propuesta de este estudio podra ser trabajar
con paquetes de partculas en lugar de hacerlo con una sola
partcula, esperando mejoras en la prediccin del comportamiento
de secado.
El enfoque lagrangiano, o modelado de trayectoria, proporciona una
descripcin directa del flujo de partculas al llevar un registro del
movimiento de partculas individuales, por lo que es un mtodo
conveniente para revelar la naturaleza discreta y transiente del
movimiento de las partculas. Sin embargo, cuando se trabaja con
paquetes de partculas, es necesario incluir tratamientos
estadsticos basados en tcnicas tales como el mtodo numrico de
Monte Carlo (24), que requieren grandes capacidades de cmputo
para efectuar promedios estadsticos y clculos de la interaccin
entre las fases slida y gaseosa. Esta es la principal limitacin de
este enfoque.
204
Anlisis paramtrico de humedad
Hasta el momento se ha considerado que los parmetros de
operacin que gobiernan el desempeo de los secadores se
mantienen constantes durante su funcionamiento. Sin embargo,
sera interesante observar cmo cambia la curva de humedad de
los slidos con la variacin de estos parmetros.
A continuacin se presentan unas grficas que estudian lo anterior,
tomando como ejemplo al secador de partculas de madera. Las
variables de control consideradas fueron dimetro promedio de
partcula, flujo de producto hmedo, flujo msico de aire,
temperatura de entrada del aire y contenido inicial de humedad en
el producto. Las condiciones de referencia son las mismas que
aparecen en la tabla 1. Se han tomado variaciones de
aproximadamente el 10% y 20% en cada parmetro, por encima y
por debajo de la referencia.

205
FIGURA 7.29 EFECTO DE LA VARIACIN DEL TAMAO DE

PARTCULA
FIGURA 7.30 EFECTO DE LA VARIACIN DEL FLUJO DE

PRODUCTO HMEDO
206
FIGURA 7.31 EFECTO DE LA VARIACIN DEL FLUJO MSICO

DE AIRE DE SECADO
FIGURA 7.32 EFECTO DE LA VARIACIN DE LA

TEMPERATURA INICIAL DEL AIRE
207
FIGURA 7.33 EFECTO DE LA VARIACIN DEL CONTENIDO DE

HUMEDAD INICIAL DEL PRODUCTO
208
CAPTULO 8
8. CONCLUSIONES Y RECOMENDACIONES
CONCLUSIONES
I. Se ha desarrollado un modelo de simulacin simple y verstil,
aplicable a secadores rotativos de diferentes dimensiones,
operando en flujo paralelo o en contraflujo. El modelo est basado
en principios bsicos de termodinmica, cintica del movimiento,
dinmica de fluidos y transferencia de calor y masa. No se han
utilizado relaciones empricas ligadas a estudios previos sobre
secadores rotativos.
II. En base a los resultados experimentales y tericos obtenidos por
otros investigadores, y a las comparaciones establecidas en este
trabajo, se puede decir que el modelo propuesto, para los dos casos
estudiados, es capaz de predecir los perfiles de humedad para los
slidos y temperatura del aire en el secador con un porcentaje bajo
de error.
209
III. Considerando la similitud en las curvas de humedad obtenidas
mediante los modelos A, B, y lagrangiano, el efecto de la difusin
interna de calor y masa en los slidos, considerada en los primeros
dos modelos, es de menor importancia al predecir el
comportamiento de secado bajo las condiciones examinadas.
IV. Bajo la suposicin de que los perfiles de humedad y temperatura
obtenidos por los modelos A y B son correctos, el haber ignorado
los efectos de radiacin y prdidas de calor en el secador no
introdujo mayores errores en la simulacin. Vale la pena indicar que
el modelo A s considera las prdidas de calor, mientras que el
modelo B considera ambas, las prdidas de calor y la transferencia
de calor por radiacin.
V. El ngulo dinmico de reposo, caracterstica inherente del producto,
tiene un efecto considerable sobre la altura promedio de cada, la
que a su vez repercute en el tiempo de secado y en el nmero de
ciclos necesarios para alcanzar determinada condicin de humedad
a la salida del secador.
VI. El perfil de humedad del producto es bastante sensible a los
cambios de flujo msico de producto, flujo de aire, y temperatura del
aire a la salida del calentador. Cuando el contenido de humedad
inicial de los slidos cambia, la curva de humedad versus posicin
axial en el tambor se desplaza verticalmente. Pequeas variaciones

210
en el tamao promedio de partcula no generan mayores
repercusiones en el perfil de humedad.
VII. En este trabajo se ha probado, en base a las comparaciones
establecidas contra los dos casos de estudio presentados, que es
viable analizar un secador rotativo desde una perspectiva
lagrangiana, rastreando a una partcula representativa a lo largo de
su paso por el secador.
En este caso, la representatividad de la partcula se midi en base
al tamao promedio de la muestra. A no ser que la distribucin de
tamaos cambie de manera significativa a lo largo del secador, ya
sea por desmoronamiento o por encogimiento del producto, la
precisin en la prediccin del perfil de humedad a lo largo del
secador depende ms del control en los parmetros de operacin
que de la representatividad que pueda tener o no la partcula. Las
figuras 7.29 a 7.33 demuestran lo anterior.

211
RECOMENDACIONES
a. Para ejecutar el programa de simulacin solamente se requiere
conocer la geometra de la seccin del tambor, y tres caractersticas
del producto (ngulo dinmico de reposo, densidad real, y dimetro
promedio de partcula), sin embargo, la integracin de los
resultados resulta un poco complicada para alguien que no est
familiarizado con la secuencia de anlisis o la estructura de
programacin, por lo que se recomienda mejorar estos aspectos
para que el software adquiera una mayor aplicabilidad.
b. Es importante analizar la naturaleza de los slidos antes de aplicar
el modelo propuesto. Por ejemplo, se debe conocer cmo vara el
ngulo dinmico de reposo con el contenido de humedad del
producto. Esto est relacionado directamente con la pegajosidad del
producto.
c. Se debe tener cuidado en seleccionar el dimetro de partcula ms
representativo de la muestra, que compute no solamente la media
en la distribucin de tamaos, sino tambin la forma de la partcula
o esfericidad. Lo anterior tiene un efecto directo sobre el flujo total
de partculas y sobre el rea efectiva de transferencia de calor y de
masa.
El modelo planteado puede no ser tan efectivo en el caso de
materiales que se rompan con el impacto o que se encojan durante

212
el secado. El software s permite trabajar con un dimetro variable
de partcula, pero hay que ser cauto al momento de cuantificar
dicha variacin.
d. El coeficiente de arrastre tiene un efecto considerable sobre el
desplazamiento de la partcula y, por ende, sobre la longitud final
del secador. Es importante efectuar pruebas experimentales de
arrastre para cada tipo de material, donde se tenga presente la
forma de las partculas y el efecto de la interaccin entre ellas.
e. Es necesario llevar a cabo pruebas experimentales que verifiquen,
de manera cuantitativa, hasta qu punto o para qu rango de
valores es aplicable el modelo y qu restricciones presenta. Se
debera estudiar la naturaleza de los materiales utilizados, y la
respuesta del modelo propuesto a cambios en los parmetros de
operacin.
213
APNDICES
APNDICE A
GRUPOS ADIMENSIONALES DE TRANSFERENCIA DE

CALOR Y MASA
Grupos adimensionales seleccionados de transferencia de calor y masa
Grupo Definicin Interpretacin
Nmero de Lewis Razn de las difusividades
( ) trmica y de masa
Nmero de Nusselt Gradiente de temperatura
( ) adimensional en la superficie
Nmero de Prandtl Razn de difusividades de
( ) momento y trmica
Nmero de Reynolds Razn de fuerzas de inercia y
( ) viscosas
Nmero de Schmidt Razn de difusividades de
( ) momento y de masa
Nmero de Sherwood Gradiente de concentracin
( ) adimensional en la superficie
Nmero de Stanton
Nmero de Nusselt modificado
( )
Nmero de Stanton
Nmero de Sherwood
para transferencia de
modificado
masa ( )
APNDICE B
RESULTADOS OBTENIDOS EN EES PARA EL SECADOR DE

PARTCULAS DE MADERA
APNDICE C
RESULTADOS OBTENIDOS EN EES PARA EL SECADOR DE

CONCENTRADOS DE ZINC
REFERENCIAS BIBLIOGRFICAS
1. Woodcock, C.R. Cohesion and adhesion. Bulk Solids Handling. Suffolk :

Chapman & Hall, 1987, pp. 30-41.
2. Shih Fan, Liang and Zhu, Chao. Size and Properties of Particles. [ed.] Arvind
Varma. Principles of Gas-Solid Flows. New York : Cambridge University Press,
1998, 1, pp. 3-10.
3. Treybal, Robert. Secado. Operaciones de transferencia de masa. Mxico :

McGraw-Hill, 1980, pp. 723-758.
4. Moyers, Charles. Psychrometry, Evaporative Cooling, and Solids Drying. [book

auth.] Robert Perry. [ed.] Robert H. Perry. Perry's Chemical Engineers' Handbook.
Seventh. New York : McGraw-Hill, 1997, 12, pp. 25-90.
5. Mujumdar, Arun. Rotary Drying. [book auth.] Taylor & Francis Group. Handbook
of Industrial Drying. Florida : CRC Press, 2006, pp. 151-171.
6. Studies in rotary drying. Friedman, S. J. and Marshall, W. R. 45, s.l. : American

Institute of Chemical Engineers, 1949, Chemical Engineering Progress, pp. 484,
573.
7. Behie, Leo. The Movement of Solids Through Flighted Rotating Drums. Part I:
Model Formulation. The Canadian Journal of Chemical Engineering. Junio 1993, pp.
337-346.
8. Barrozo, M.A.S. A study about particle motion in rotary dryers. Brazilian Journal
of Chemical Engineering. Febrero 2007, pp. 365-374.
9. Rotary Driers and Coolers for Granular Fertilizers. Schofield, F. R. and Glikin, P.
G. 40, 1962, Trans. Inst. Chem. Eng, pp. 183-190.
10. Modelado de Secaderos Rotatorios en Isocorriente. Castao, F., Rubio, F. R.

and Ortega, M. G. 2009, Revista Iberoamericana de Automtica e Informtica
Industrial, pp. 32-43.
11. Kamke, Frederick A. Engineering Analysis of a Rotary Dryer: Drying of Wood

Particles. 1983. pp. 43-44, 89-130, PhD Thesis.
12. Ajayi, Olukayode Oludamilola. Multiscale modelling of industrial flighted rotary

dryers. James Cook University. 2011. pp. 1-3,22-40,145-153, PhD Thesis.
13. Glikin, P. G. 56, 1978, Trans. Inst. Chem. Eng., p. 120.
14. Modeling and experimental study of hydrodynamic and drying characteristics of

an industrial rotary dryer. Barrozo, M. A. S. 2, Uberlandia : s.n., Septiembre 26,
2009, Brazilian Journal of Chemical Engineering, Vol. XXVI, pp. 331-341. ISSN
0104-6632.
15. Barrozo, M.A.S. and Arruda, E.B. Acerca de nosotros: IOPscience. sitio web de
IOP Science. [Online] 2008. [Cited: 11 2, 2013.] http://iopscience.iop.org/1742-
6596/135/1/012063.
16. Stroshine, Richard. Friction of solids and flow of granular solids. Physical
Properties of Agricultural Materials and Food Products. 2004, pp. 239-254.
17. The Design of Flights in Rotary Dryers. Revol, D., Briens, C. L. E. and
Chabagno, J. M. 2001, Powder Technology, Vol. 121, pp. 230-238.
18. The kinetic angle of repose of powders. Kelly, J. J. 21, Dublin : Society of
Costmetic Chemists of Great Britain, 1970, Journal of the society of cosmetic
chemists, pp. 37-51.
19. White, Frank M. Flujo alrededor de cuerpos. Mecnica de Fluidos. Quinta.

Mxico : McGraw Hill, 2008, pp. 449-501.
20. Incropera, Frank. Flujo externo. Fundamentos de Transferencia de Calor.

Mxico : Prentice Hall, 1999, pp. 346-394.
21. . Introduccin a la conveccin. Fundamentos de Transferencia de Calor.

Mxico : Prentice Hall, 1999, pp. 284-332.
22. Nellis, Gregory and Klein, Sanford. Mass Transfer. Heat Transfer. New York :
Cambridge University Press, 2008, pp. 974-978.
23. Leu, Shao Yuan. Robust enzymatic saccharification of a Douglas-fir forest

harvest residue by SPORL. Biomass and Energy. August 1, 2013, p. 5.
24. Shih Fan Llian, Zhu Chao. Basic Equations. [ed.] Varma A. Principles of Gas-
Solid Flows. New York : Cambridge University Press, 1998, 5, pp. 164-166.
25. Incropera, Frank. Transferencia de masa por difusin. Fundamentos de

Transferencia de Calor. Mxico : Prentice Hall, 1999, pp. 784-817.
26. Sper, Jos and Torres, Ivn. Proyecto de Graduacin. Escuela Superior
Politcnica del Litoral, Guayaquil : 2009.
27. Behie, Leo. The Movement of Solids Trough Flighted Rotating Drums. Part II:
Solids-Gas Interaction and Model Validation. The Canadian Journal of Chemical
Engineering. Abril 1994, pp. 240-248.
28. Wang, T. J. Granule transport and mean residence time in horizontal drum with
inclined flights. Powder Technology. Octubre 12, 2005, pp. 50-58.
29. Amanifard, N. and Haghi, A.K. On the drying process of a solid porous body:
some important criteria, trends and developments. [book auth.] G.E. Zaikov.
Progress in Chemical, and Biochemical Physics, Kinetics and Thermodynamics. New
York : Nova, 2008.
30. Malode, Parag. DEM analysis of the role of lifter design on performance on
rotary cooler. Jamshedpur. 2010, pp. 391-397.
31. Gas velocity effects on heat transfer in direct heat rotary dryers. McCormick, P.
Y. New York : American Institute of Chemical Engineers, 1962, pp. 57-61. ISBN.
32. Cengel, Yanus A. and J., Ghajar Afshin. Radiation Heat Transfer. Heat and
Mass Transfer: Fundamentals & Applications. Fourth. New York : McGraw Hill, 2011.
33. Price, James F. Woods Hole. [Online] June 7, 2007. [Cited: 10 10, 2013.]
http://www.whoi.edu/science/PO/people/jprice. 02543.
34. Barrozo, M.A.S. and Arruda, E.B. Experimental and Simulation Study of
Fertilizer Drying in a Rotary Dryer. Escuela de Ingeniera Qumica, Universidad
Federal de Uberlndia. Minas Gerair : s.n., 2011. Paper.
35. Modelling of Fertilizer Drying in a Rotary Dryer: Parametric Sensitivity Analysis.

Barrozo, M.A.S. and Arruda, E.B. 2, Minas Gerais : s.n., Septiembre 22, 2011,
Brazilian Journal of Chemical Engineering, Vol. XXIX, pp. 359-369. ISSN 0104-6632.
36. Kamke, F. A. and B., Wilson J. Computer Simulation of a Rotary Dryer. AlChE
Journal. February 1986, Vol. 32, 2, pp. 269-275.
37. Rastikian, K., Capart, R. and Benchimol, J. Modelling of sugar drying in a

countercurrent cascading rotary dryer from stationary profiles of temperature and
moisture. Journal of Food Engineering. 1999, 41, pp. 193-201.
38. Konidis, Jordan. Design of Direct Heated Rotary Dryers. Department of

Mechanical Engineering, Concordia University. Montreal : s.n., 1984. Technical
Report.
39. Drying and Heating Modelling of Granular Flow: Application to the Mix-Asphalt
Processes. LeGuen, L. 2, Bouguenais : s.n., 2011, Journal of Applied Fluid
Mechanics, Vol. 4, pp. 71-80. ISSN 1735-3572.
40. Price, James F. Lagrangian and Eulerian Representation of Fluid Flow:

Kinematics and the Equations of Motion. Woods Hole Oceanographic Institution.
Woods Hole : s.n., 2006. Essay.
41. Fakuori, Alireza. A new experimental method for flight design in rotary dryers.
Department of Mechanical Engineering, URMIA University of Technology.
42. Lee, Andrew. Modelling the Solids Transport Phenomena Within Flighted Rotary
Dryers. James Cook University. Queendland : s.n., 2008. PhD Thesis.
43. Mathematical Modeling and Computer Simulation of Continuous Rotary Dryer.

Hashemabadi, S. H. Narmak : s.n., 2004. 9th Iranian Chemical Engineering
Congress.
44. On modelling and control of a rotary sugar dryer. Pierce, Robert. San Diego :
Pergamon, 2001, Control Engineering Practice, pp. 249-266.
45. Cairo, Nadia, Colangelo, Gianpiero and Starace, Giuseppe. Performance
analysis of two industrial dryers (cross flow and rotary) for ligno-cellulosic biomass
desiccation. Dipartimento di Ingegneria dell'Innovazione, Universit del Salento.
Lecce : s.n., 2012.
46. Process analysis for an alfalfa rotary dryer using an improved dimensional
analysis method. Xianzhe, Zheng. 3, Harbin : s.n., 2009, Int J Agric & Biol Eng, Vol.
2, pp. 76-82.
47. Shih Fan, Liang and Zhu, Chao. Size and Properties of Particles. [ed.] Arvind
Varma. Principles of Gas-Solid Flows. New York : Cambridge University Press,
1998, 1, pp. 3-10.
48. Bostock, Pierce. Acerca de nosotros: Thermal Energy Systems. sitio web de
Thermal Enery Systems. [En lnea] 22 de Abril de 2009. [Citado el: 10 de Noviembre
de 2012.] http://www.thermalenergysystems.com/bsst/papers/2009agm3.pdf.

Tesis Carlos Durán

Cargado por

Copyright:

Formatos disponibles

Tesis Carlos Durán

Cargado por

Información del documento

Descripción original:

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Tesis Carlos Durán

Cargado por

Copyright:

Formatos disponibles

ESCUELA SUPERIOR POLITCNICA DEL LITORAL

Facultad de Ingeniera en Mecnica y Ciencias de la

Desarrollo de una Tcnica Lagrangiana para el Dimensionamiento

Previa a la obtencin del Ttulo de:

Carlos Geovanny Durn Oscuez

Al Ing. Vicente Adum Gilbert,

realizacin de este trabajo.

La responsabilidad del contenido de esta

Tesis de Grado, me corresponde

exclusivamente; y el patrimonio intelectual

de la misma a la ESCUELA SUPERIOR

(Reglamento de Graduacin de la ESPOL)

En este trabajo se pretende probar la hiptesis de que es factible analizar

secadores rotativos desde una perspectiva lagrangiana, la cual combine la

transferencia de calor y masa con el movimiento de una partcula individual a

travs del secador.

Un secador rotativo es un cilindro metlico giratorio dispuesto de forma

humedad del producto. En el interior del cilindro se encuentran montadas

unas paletas que recogen el material granular desde el fondo y lo descargan

hacia la corriente de aire para formar una cascada de partculas.

En la primera parte del estudio, se desarrolla un anlisis termodinmico del

para secar cierta cantidad de producto hmedo. Luego se establece un

ngulo de cada representativo de toda la cortina de slidos, con el cual se

estima la altura promedio de cada. A partir de all, el anlisis se concentra en

repetir el ciclo de levantamiento y cada a lo largo del secador.

Durante su recorrido, se determina la tasa de transferencia de calor y masa

en pasos de tiempo definidos, registrando las variaciones en las propiedades

de la partcula y del aire, y sumando los desplazamientos sucesivos de la

partcula hasta obtener la longitud de tambor correspondiente al nivel de

humedad deseado en el producto.

Para medir la validez de este planteamiento, se dise un software que

permiti aplicar la metodologa propuesta a dos secadores rotativos

previamente analizados con otros modelos existentes, uno para partculas de

madera y otro para concentrados de zinc. De este modo, se compararon los

datos tericos obtenidos con el modelo propuesto contra aquellos de los

modelos existentes, y contra la informacin experimental disponible.

del modelo propuesto y las del modelo existente para el contenido de

humedad final en la madera. Asimismo, se encontr un buen ajuste entre el

perfil de humedad terico generado por el modelo propuesto y los datos

experimentales de humedad para concentrados de zinc, en funcin de la

posicin axial en el secador, obtenindose un error relativo de 4.62% en el

contenido de humedad final del producto.

Aceleracin de la gravedad m/s2

RO Distancia desde el centro del tambor hasta m

Figura 5.2 Intercambio de calor latente y sensible en una interfaz gas-

Figura 7.26 Humedad de partcula vs Posicin (Dimetro de Partcula

Los secadores continuos de tipo rotativo son ampliamente utilizados para

secar productos granulados en la industria alimenticia, farmacutica, mineral

y de fertilizantes. Su desempeo est gobernado por tres factores:

movimiento, transferencia de masa y transferencia de calor. Estos factores

La literatura presenta modelos matemticos basados en la determinacin

experimental de transferencia de masa y coeficientes globales de

transferencia de calor para obtener los perfiles de humedad y temperatura

aplicables solamente a ciertos tipos de secadores, con configuraciones y

parmetros de operacin particulares.

Un modelo simple, rpidamente ajustable a diferentes condiciones de trabajo,

le permitira al diseador contar con una herramienta para analizar el

comportamiento de secado en funcin de los parmetros de operacin. Con

ello, eventualmente se podra optimizar el diseo para disminuir los costos