PRÁCTICA DE SECADO DE SÓLIDOS

MODELO DE CÁLCULO

1. Método experimental para el periodo de velocidad constante.

El secado se refiere a la eliminación parcial de agua (humedad) de un sólido por

medios difusivos y térmicos. Durante el proceso de reducción del contenido de
humedad en un sólido lo que se pretende es obtener el tiempo requerido para lograr el
grado de secado exigido, determinar el tamaño del equipo y las condiciones de
operación de humedad y temperatura del aire empleado. Por consiguiente, a partir de
los datos obtenidos en el laboratorio se recalculan los datos obtenidos por medio de la
ecuación (1) que define la humedad de una mezcla aire-vapor de agua.
W −W S
Xt= (1)
WS
Donde,
W : Sólido húmedo (kg sólido seco + kg de agua)
W S: Sólido seco (kg sólido seco)
X t : Humedad en el tiempo (kg H2O/kg sólido seco)
Una vez establecidas las condiciones de secado constante, se procede a calcular el
contenido de humedad en equilibrio, por medio del calculo de humedad libre para cada
valor obtenido anteriormente por medio de:
X =X t −X ¿ (2)
Donde,
X : Humedad libre (kg totales de agua/kg sólido seco)
X t : Humedad en el tiempo (kg H2O/kg sólido seco)
X ¿ : Humedad en equilibrio
Utilizando los datos de humedad en equilibrio es posible trazar la curva del contenido
de humedad en función del tiempo con el fin de medir las pendientes de las tangentes
a dicha curva y determinar así la velocidad de secado ( R ) y trazar la curva de
velocidad de secado, así:
−W S dX
R= (3)
A dt
Al tratarse de un periodo a velocidad constante,
dX
m= (4 )
dt
Donde,
A : Área superficial de la esfera (m 2)
2. Método experimental para el periodo de velocidad decreciente.
Para el periodo de velocidad decreciente, se obtienen una gráfica de humedad libre vs
tiempo, tal como en el periodo anterior y se realiza un ajuste por medio de regresión
lineal y así obtener la expresión de m derivando la ecuación de la curva.
Una vez encontrado el valor de m , reemplazar en la ecuación (3) para cada valor
obtenido y así analizar el comportamiento de la humedad en un periodo de velocidad
decreciente.
3. Coeficiente de transferencia de calor a partir de correlaciones empíricas.

Experimentalmente, este coeficiente se puede determinar como:

R c λ agua
h= (5)
T −T w
Donde,

Rc : Es la velocidad de secado para el periodo constante, en unidades de kg agua/s-

m2.
λ agua: Es la entalpia de vaporización del agua, en unidades de kJ/kg.
T : Es la temperatura de bulbo seco en K.
T w: Es la temperatura de bulbo húmedo en K.

Las temperaturas de bulbo seco y húmedo se obtienen experimentalmente, mientras

que el valor de Rc es el determinado por el promedio de la velocidad constante.
mientras que λ agua se determina de la siguiente forma:

λ agua=−2.4285T w +3166.2 (6)

La determinación teórica de este coeficiente se lleva a cabo por medio de la definición

de Nusselt:

Nu k
h= (7)
Lc
Donde,

Nu: relaciona la transferencia de calor por convección con la transferencia de calor por
conducción.
k : Es la conductividad termina del aire a T, en unidades de W/m-K.
Lc : Es la longitud característica de la partícula en m, en este caso será el diámetro de
la esfera.

Por lo dicho anteriormente, se define Reynolds para una superficie esférica como:

v ρD
ℜD = (8)
μ
Donde,

v: Representa la velocidad del aire en m/s.

ρ : Se refiera la densidad del aire evaluada a la temperatura de película en kg/m3.
D : Es el diámetro de la esfera en m.
μ: Es la viscosidad dinámica a la temperatura de película.

Para flujo externo en una superficie esférica existen numerosas correlaciones de

transferencia de calor propuestas, y Whitaker recomienda una expresión de la forma:

1 2 1 /4
´ D =2+ 0.4 ℜ D 2 + 0.06 ℜ D 3 Pr 0.4 μ
( )
Nu
μs ( ) (9)

Con las siguientes restricciones:

0.71< Pr <380

3.5< ℜD <7.6 x 104

1.0< ( μμ )<3.2
s

Todas las propiedades excepto μs se evalúan en T .͚

En esta correlación, se expresa el número adimensional de Prandtl, que es la relación

entre la velocidad de difusión de la cantidad de movimiento y la velocidad de difusión
de calor y se define como1:
v Cp μ
Pr= = (10)
α k

4. Coeficiente de transferencia de masa a partir de correlaciones empíricas

El coeficiente de transferencia de masa, se determina experimentalmente de la siguiente forma:

Rc
k y= (11)
y s− y
Donde,

y s : Representa la humedad de saturación.

y : Representa la humedad, en unidades de kg agua/kg aire.

Rc : Se determinado del promedio de velocidades constantes; mientras que y s e y se determinan
de la siguiente forma:

P sat
A
ys = sat (12)
P−P A
1
INCROPERA, Frank P. & DE WITT, David P. Fundamentos de transferencia de calor. Cuarta edición. Ciudad de México. Pearson,
1999 extraido : 3/10/2016
Donde,

Psat
A : Es la presión de vapor del agua pura a la temperatura establecida, en unidades de kPa.

P: Es la presión total en kPa, normalmente es la presión atmosférica del sistema.

Por otro lado, la definición de humedad solo depende de la presión parcial ( P A ) del vapor de
agua y de la presión total:
PA
y= (13)
P−P A
La determinación teórica del coeficiente se lleva a cabo teniendo en cuenta del número
adimensional análogo a Nusselt, denominado número de Sherwood:
Sh D AB
k y= (14)
Lc
Donde,

Sh: representa el cociente entre la transferencia de masa por convección y difusión.

D AB: Es la difusividad del agua en el aire, en unidades de m2/s.

Se define el número de Schmidt que es el cociente entre la difusividad de cantidad de

movimiento y la difusividad de masa como:

μ
Sc= (15)
ρ D AB
Como no hay convección formada, en la siguiente correlación el número de Sherwood es 2, por
lo tanto, la correlación se define de la siguiente forma:

Sh=2+C ℜm Sc 1/ 3 (16)

Donde C y m son constantes de la correlación.

ANÁLISIS DE RESULTADOS

Se establecieron condiciones de operación y especificaciones del sólido (Piedra

pómez), presentes en la tabla 1.

Tabla 1. Condiciones de operación

%H Entrada 0.109
T aire entrada (K) 324.25
desfera(m) 0.3243
Aesfera(m2) 0.0066
W sólido seco (Kg) 0.0194
Una vez establecidas las condiciones de operación se registraron lecturas de
variaciones de masa y porcentaje de humedad relativa a intervalos de 20 segundos
por 20 minutos (Tabla 2) y a intervalos de 4 minutos (Tabla 3) hasta no obtener
variaciones en el peso.

Tabla 2. Datos experimentales intervalos de 20s

Tiempo(h) W: Masa (kg) %HR a la salida T aire salida(K)

0 0.03613 0.114 325.15
0.00555556 0.0361 0.115 325.05
0.01111111 0.03606 0.116 324.95
0.01666667 0.03596 0.118 324.85
0.02222222 0.0359 0.119 324.75
0.02777778 0.03581 0.12 324.75
0.03333333 0.03577 0.121 324.65
0.03888889 0.03568 0.122 324.55
0.04444444 0.03567 0.123 324.45
0.05 0.0356 0.123 324.45
0.05555556 0.03558 0.123 324.35
0.06111111 0.03553 0.124 324.25
0.06666667 0.03544 0.124 324.15
0.07222222 0.03538 0.125 324.05
0.07777778 0.03533 0.125 323.95
0.08333333 0.03526 0.125 323.95
0.08888889 0.0352 0.126 323.85
0.09444444 0.03513 0.126 323.75
0.1 0.03509 0.127 323.65
0.10555556 0.03503 0.128 323.55
0.11111111 0.03498 0.128 323.55
0.11666667 0.03489 0.129 323.45
0.12222222 0.03481 0.129 323.45
0.12777778 0.03481 0.129 323.35
0.13333333 0.03477 0.13 323.35
0.13888889 0.03466 0.13 323.25
0.14444444 0.03461 0.13 323.25
0.15 0.03458 0.131 323.15
0.15555556 0.03451 0.132 323.05
0.16111111 0.03447 0.132 323.05
0.16666667 0.03437 0.132 323.05
0.17222222 0.03422 0.133 323.05
0.17777778 0.03418 0.133 322.95
0.18333333 0.03412 0.133 322.95
0.18888889 0.03406 0.133 322.95
0.19444444 0.034 0.133 322.95
 0.2 0.03392 0.133 322.95
0.20555556 0.03388 0.133 322.95
0.21111111 0.03385 0.133 322.95
0.21666667 0.03378 0.133 322.95
0.22222222 0.03372 0.133 322.95
0.22777778 0.03366 0.132 322.95
0.23333333 0.03364 0.132 322.95
0.23888889 0.03354 0.131 322.95
0.24444444 0.03347 0.131 322.95
0.25 0.03345 0.131 322.95
0.25555556 0.03338 0.131 322.95
0.26111111 0.0333 0.131 322.95
0.26666667 0.03325 0.131 322.95
0.27222222 0.03311 0.13 322.95
0.27777778 0.0331 0.13 322.95
0.28333333 0.03309 0.13 322.95
0.28888889 0.03302 0.13 322.95
0.29444444 0.03299 0.13 322.95
0.3 0.03289 0.13 323.05
0.30555556 0.03285 0.13 323.05
0.31111111 0.0328 0.129 323.05
0.31666667 0.03272 0.128 323.05
0.32222222 0.03268 0.128 323.05
0.32777778 0.03264 0.127 323.15
0.33333333 0.0326 0.127 323.15

Tabla 2. Datos experimentales intervalos de 4min

Tiempo(h) W: Masa (kg) %HR a la salida T aire salida(K)

0.4 0.02899 0.13 323.15
0.46666667 0.02872 0.129 323.05
0.53333333 0.02843 0.127 323.15
0.6 0.02819 0.119 323.65
0.66666667 0.02797 0.111 324.35
0.73333333 0.02775 0.101 325.35
0.8 0.02755 0.113 325.35
0.86666667 0.02739 0.122 324.65
0.93333333 0.02719 0.125 323.95
1 0.02705 0.126 323.85
1.06666667 0.02689 0.125 323.75
1.13333333 0.02673 0.117 324.05
A partir de los datos obtenidos se calcula el contenido de humedad ( X t ) y la humedad
libre ( X ) para obtener la gráfica de humedad libre vs tiempo.

En la gráfica 1 es posible interpretar el comportamiento de la humedad del sistema

respecto al tiempo, estos datos fueron tomados cada 20 segundos durante 20 minutos.
Con esta gráfica se puede observar que la humedad disminuye con respecto al
tiempo, esto se debe a que el túnel de viento retira la humedad y por ende a medida
que el tiempo transcurre el sólido esférico contiene cada vez menos agua. Se ajustan
los datos a un polinomio de orden 3.

Humedad vs tiempo (20 s)

0.5
f(x) = 1.23 x³ − 0.59 x² − 0.48 x + 0.48
0.45 R² = 1

0.4
Xlibre

0.35

0.3

0.25

0.2
0 0.05 0.1 0.15 0.2 0.25 0.3 0.35
t(s)

Gráfico 1. Humedad libre vs tiempo

En la gráfica 2 se puede interpretar el comportamiento de la humedad respecto al

tiempo para el segundo periodo en el cual se miden las variaciones de masa, y
humedad relativa cada 4 min hasta que no se obtienen variaciones importantes de
peso. También se ajustan los datos a un polinomio de orden 3.
 Humedad vs tiempo (4min)
0.14

0.12
f(x) = − 0.06 x³ + 0.22 x² − 0.38 x + 0.24
0.1 R² = 1

0.08
xllibre

0.06

0.04

0.02

0
0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 1.1 1.2
t(s)

Gráfico 2. Humedad libre vs tiempo

Teniendo en cuenta el cambio de humedad libre en el tiempo, se determinó la

velocidad de secado de la piedra pómez en función del tiempo. En la gráfica 3 se
puede observar el cambio de la velocidad de secado respecto a la humedad libre.

Velocidad de secado vs Xlibre

8.00E-05
7.00E-05
6.00E-05
R(Kg agua/m2S)

5.00E-05
4.00E-05
3.00E-05
2.00E-05
1.00E-05
0.00E+00
0 10 20 30 40 50 60 70
X libre(kg H2O/ kg ss)

Gráfica 3. Velocidad de secado vs Xmedia

El comportamiento de la velocidad tal como se observa en la gráfica 3 es inicialmente

creciente, luego llega a un periodo donde la velocidad es constante y luego empieza a
decrecer. Al principio la velocidad de secado incrementa debido a que tiene el mayor
porcentaje de humedad en la superficie del sólido, después de este periodo creciente
llega el periodo de velocidad constante tal como se esperaba del gráfico velocidad
constante vs Xlibre, se tiene un periodo de velocidad constante entre X = 0.31612-
0.4287 kgH2O/kgSS, , en este intervalo la superficie del solido se encuentra húmeda y el
agua superficial se evapora como si el sólido no hiciera parte del sistema, lo que
provoca que se presente un secado a velocidad constante. La humedad crítica se
presenta en X= 0.3757 kgH2O/kgSS, esto significa que a partir de ese momento la
superficie del solido se encuentra seca y el plano de evaporación se adentra en el
sólido lo que hace más difícil el proceso de secado y ocasiona que la velocidad de
secado disminuya.

Tabla 1. Error promedio en el cálculo de ambas constantes.

%Error Ky-Treybal 99.85%

%Error Ky-Froesling 106%
%Error h 13381.26%

El cálculo de los coeficientes de transferencia presenta un error relativo demasiado

alto, que se puede atribuir a la lectura experimental de los datos en el laboratorio y
posiblemente a que los modelos teóricos correlacionados no se ajusten
adecuadamente a las condiciones experimentales.
Respecto al coeficiente de transferencia de calor(h), es el resultado del intercambio de
calor entre el sólido esférico y el flujo de aire caliente que pasa dentro del túnel. El
comportamiento del coeficiente de transferencia de calor es similar al de la velocidad
de secado, por lo que experimentalmente es el modelo se ajusta, sin embargo, por el
modelo teórico se presenta una oscilación del valor de este. Esto corrobora que las
correlaciones no se ajustan a las condiciones experimentales. Sucede un fenómeno
similar con el coeficiente de transferencia de masa(Ky), sin embargo, la oscilación del
valor teórico es mucho menor lo que se reflejado directamente en el porcentaje de
error relativo dado que este es menor comparado con el error relativo del coeficiente
de transferencia de calor con un porcentaje de error de 99.85%.
Estos errores se deben a que para realizar los cálculos de los valores teóricos para
cada uno de los coeficientes se utilizan correlaciones que utilizan valores reportados
en la literatura para las diferentes propiedades de los fluidos y condiciones necesarias
propias de cada correlación que muchas veces no se cumple durante la práctica.

Otra posible causa de error es que durante el registro de datos experimentales, se

esperaría que al llegar al contenido de humedad de equilibrio, la masa del sólido y por
ende la humedad se volvieran constantes, ya que no es posible retirar más humedad
por medio de secado; sin embargo, se observa que esto no ocurre lo que indica que si
se registran más datos para el secado se podría llegar a un valor diferente para la
humedad de equilibrio y cambiar el valor de Rc, lo que modificaría los resultados
experimentales de los coeficientes de trasferencia de calor y masa.
CONCLUSIONES

 Se conocieron las operaciones unitarias del secado directo mediante el uso de

un secador discontinuo de túnel de viento en condiciones constantes de flujo
de aire y temperatura.

 Se construyeron las curvas de secado (humedad contra tiempo, velocidad de

secado contra humedad) para el caso de la piedra pómez.

 Se determinaron los coeficientes convectivos de transferencia de calor (hc) y

masa (Ky) para el período antecrítico y se compararon con los resultados
obtenidos por correlaciones empíricas.

 Se determinaron experimentalmente los coeficientes convectivos de

transferencia de calor (h) con un valor de 31.35 W/m2K y masa Ky (Kg As/m2
S)-Treybalde 0.413, Ky-Froesling 5.06 E-5 para el período antecrítico.
Obteniendo porcentajes de error de 13381.26% y 99.85% y 106%
respectivamente.

 Se determinó el mecanismo de secado en el período poscrítico.

 Es necesario que el sólido alcance la humedad de equilibrio para tener mayor

precisión en la determinación de cada uno de los resultados.

 Se logró modelar el secado de un sólido identificando todas las variables

relacionadas y características del proceso en función del tiempo.

CRITERIOS DE DESEMPEÑO

Diseño Comunicación Investigación Mejoramiento continuo

Se hizo uso Se realiza una Se procuró evitar
de Excel Identificando los previa investigación errores en la
como objetivos de la de los conceptos presentación de
herramienta práctica, se expresó básico de secado informes escritos que
computaciona de manera clara la (humedad, hayan sido cometidos
l para resolver resolución y la velocidad y tiempos con anterioridad y se
de forma adaptación de los críticos) de sólidos identificaron las
adecuada las modelos y se buscó variables del secado
ecuaciones experimentales y con implementar las que pueden ser
científicas base a los reportes de teorías. tomadas en
expresadas literatura se estimaron suposiciones y consideración para una
según los posibles errores que datos que podrían futura propuesta
principios se analizaron. relacionar errores ingenieril.
establecidos. en la práctica.

REFERENCIAS

- GEANKOPLIS, Christie J. Procesos de transporte y operaciones unitarias. Tercera

edición. Ciudad de México. Compañía Editorial Continental, S.A. 1998.
- Frank Incropera, David DeWitt. Fundamentos de Transferencia de Calor.
Prentice Hall. México 1999 (4° Ed). ISBN 0-471-30460-3. 888 pp.

- Welty, J. R., Wicks, C. E., Wilson, R. E., & Rorrer, G. L. (2007). Fundamentals
of momentum, heat, and mass transfer. Oregon: Wiley.