Ingeniería e Investigación

Universidad Nacional de Colombia

[email protected]
ISSN (Versión impresa): 0120-5609
COLOMBIA

2008
Everaldo Joaquín Montes Montes / Ramiro Torres Gallo / Ricardo David Andrade Pizarro /
Omar Andrés Pérez Sierra / José Luis Marimon Escobar / Isabel Inés Meza Herazo
MODELADO DE LA CINÉTICA DE SECADO DE ÑAME (DIOSCOREA ROTUNDATA) EN
CAPA DELGADA
Ingeniería e Investigación, agosto, año/vol. 28, número 002
Universidad Nacional de Colombia
Bogotá, Colombia
pp. 45-52

Red de Revistas Científicas de América Latina y el Caribe, España y Portugal

Universidad Autónoma del Estado de México

http://redalyc.uaemex.mx
REVISTA INGENIERÍA E INVESTIGACIÓN VOL. 28 No. 2, AGOSTO DE 2008 (45-52)

Modelado de la cinética de secado de ñame (dioscorea

rotundata) en capa delgada
Modelling the kinetics of thin-layer yam (dioscorea rotundata) drying
Everaldo Joaquín Montes Montes1, Ramiro Torres Gallo2, Ricardo David Andrade Pizarro3,
Omar Andrés Pérez Sierra4, José Luis Marimon Escobar5 e Isabel Inés Meza Herazo6

RESUMEN
El secado de ñame en capa delgada, de las variedades 9811-089 y 9811-091, fue evaluado en un secador de la-
boratorio tipo bandeja a temperaturas de 45, 55 y 70 °C y velocidad de aire promedio de 1 m/s, en el Laboratorio
de Ingeniería Aplicada de la Universidad de Córdoba; las muestras fueron adecuadas en geometría de rodajas de
radio 3,19 cm y espesor de 0,5 cm y láminas de 0,5x3x5 cm. Se evaluó el efecto de la temperatura, variedad y
geometría en el tiempo de secado, empleando un diseño completamente al azar en arreglo factorial con tres facto-
res: temperatura (45, 55 y 70 °C), geometría (rodajas y láminas) y variedad (9811-089 y 9811-091), realizando
tres repeticiones por tratamientoutilizando, dando como resultado una disminución en el tiempo de secado del
28,15% para la temperatura de 70 °C. Con los resultados obtenidos fueron construidas las curvas de secado, en
las cuales se observó que este proceso tuvo lugar en el periodo decreciente, lo que evidencia que la difusión es el
mecanismo que gobierna el secado de Dioscorea rotundata para las variedades estudiadas. Las curvas de secado
fueron ajustadas a los modelos matemáticos de Page, difusión, Thompsom, Newton, Page modificado, Henderson
y Pabis y el modelo logarítmico, determinando mediante el coeficiente de determinación (R2), error medio estimado
(SE) y desvío medio relativo (%P), que el modelo logarítmico es el que mejor describe el proceso de secado (R2 ≥
99.17 y SE ≤ 0.0299).
Palabras clave: coeficiente de difusión, Page, modelos matemáticos.

ABSTRACT
Thin-layer yam (varieties 9811-089 and 9811-091) drying was evaluated in a laboratory-type dish dryer at 45°C,
55°C and 70°C and 1 m/s average air speed in the Universidad de Cordoba’s Applied Engineering laboratory. The
samples were 3.19 cm long, 0.5 cm thick, in 0.5x3x5 cm slices. The effects of temperature, variety and geometry on
drying-time were evaluated using a completely random factorial adjustment design: temperature (45°C, 55°C and
70°C), geometry (slices and fillets) and variety (9811-089 and 9811-091). Three repetitions were made per treat-
ment, producing a 28.15% reduction in drying time at 70°C. Drying curves were constructed from the obtained re-
sults; it was observed that drying took place during the decreasing period, proving that diffusion was the mechanism
involved in Discorea rotundata drying for the studied varieties. Drying curves were adjusted to Page, diffusion,
Thompson, Newton, modified Page, Henderson and Pabis mathematical models. The logarithmic model (determi-
ned by determination coefficient (R2)) estimated mean error (SE) and relative mean deviation (%P), showing that the
logarithmic model better described the drying process (R2 ≥ 99.17 and SE ≤ 0.0299).
Keywords: diffusion coefficient, Page, mathematical model.

Recibido: febrero 11 de 2008

Aceptado: junio 27 de 2008

1
Ingeniero químico, Universidad del Atlántico, Colombia. Maestría, en Ciencia y tecnología de alimentos, Universidad federal de viÇosa, Brasil.
Profesor titular, Universidad de Córdoba. [email protected], [email protected]
2
Ingeniero químico, Universidad del Atlántico, Colombia. Especialista, Ciencias Química, Universidad de Córdoba, Colombia. Profesor asistente, U-
niversidad de Córdoba, Colombia. [email protected]
3
Ingeniero químico, Universidad del Atlántico, Colombia. Especialista, Ciencia y tecnología de alimentos, Universidad del Magdalena, Colombia.
Especialista, Ingeniería de procesos de alimentos y biomateriales, UNAD, Colombia. Profesor asistente, Universidad de Córdoba. randrade@sinu.
unicordoba.edu.co
4
Ingeniero químico, Universidad del Atlántico, Colombia. Maestría, Ingeniería Química, Universidad Industrial del Santander, Colombia. Doctora-
do, Ciencia y tecnología de alimentos, Universidad federal de viÇosa, Brasil. Profesor asistente, Universidad de Córdoba, Colombia. operez@sinu.
unicordoba.edu.co
5
Ingeniero de alimentos, Universidad de Córdoba, Colombia. Ingeniero, La Soberana. [email protected]
6
Ingeniera de alimentos, Universidad de Córdoba, Colombia. Ingenieria, Superalmacenes Olímpica (SAO). [email protected]

45
MODELADO DE LA CINÉTICA DE SECADO DE ÑAME (DIOSCOREA ROTUNDATA) EN CAPA DELGADA

Introducción Desarrollo experimental

El ñame pertenece a la familia Dioscoreacea, género Dio- Obtención y adecuación de la muestra
scorea; posee cerca de 600 especies, de las cuales D. Alata,
D. Cayenensis Lam y D. Rotundata tienen mayor importancia Los tubérculos de ñame de la variedad Dioscorea Rotundata
económica (Falde et al., 2006). Este tubérculo contiene fé- (9811-091; 9811- 089) fueron suministrados por el banco de
cula abundante (Rincón, 2000); y constituye un importante germoplasma del Programa de Biotecnología Agrícola (PBA)
ali-mento en las regiones tropicales (Guzmán y Buitrago, de la Facultad de Ciencias Agrícolas de la Universidad de
2000). Córdoba. Las muestras fueron sometidas a procesos de: se-
lección y clasificación (desechando los tubérculos en mal es-
La producción mundial de ñame está destinada casi en su tado y utilizando tubérculos de tamaños y formas similares),
totalidad a abastecer el mercado interno de los países pro- lavado superficial, pelado (retirando la cáscara mediante cu-
ductores y por lo general es comercializado sin ningún proce- chillos), corte (las muestras se cortaron en forma de láminas
samiento industrial (Floreze y Morini, 2000); este hecho, su- de 0,5x3x5 cm y rodajas de radio y espesor de 3,19 y 0,5
mado al alto contenido de humedad que posee el ñame cm. respectivamente) y luego se pesaron.
(entre 50 y 80%), provoca importantes pérdidas (20-30%)
durante la poscosecha (Guzmán y Buitrago, 2000; Sierra y Equipos y procedimiento
Munive, 2007). Se utilizó una cámara de secado experimental, provista de
El secado, definido como un proceso simultáneo de trans- un ventilador Blower de 1/4 HP para forzar el aire de seca-
ferencia de calor y masa entre el producto y el aire de se- do, un sistema de calefacción compuesto por dos resistencias
cado, consiste en la remoción de humedad excesiva conte- en alambre de 110V con recubrimiento en acrílico puestas
nida en el producto por medio de la evaporación, gene- en C, un sistema de control automático de temperatura
ralmente por convección forzada (Maldonado y Pacheco, (Makthermo Mc–17), un sistema de control de humedad
2003; Ertekin y Yaldiz, 2004); es un método clásico de pre- (Full Gauge Mt-530R) compuesto por una resistencia y una
servación de alimentos y se convierte en una excelente alter- unidad de refrigeración de 1/5 HP. Además el equipo cuenta
nativa para prolongar el periodo de vida útil del tubérculo, con dos bandejas de acero inoxidable de dimensiones 40x
así como para disminuir el peso para el transporte y reducir 30x2 cm.
el espacio requerido para su almacenamiento (Romero y Se utilizó también una balanza analítica (Mettler PS4000
Kieckush, 2003; Cardoso et al., 2004; Vega y Fito, 2005). con exactitud ± 0,01 g, un termohigrómetro (Durotherm con
Estudios en torno al secado de frutos y vegetales han deter- escala de -10 a 50 °C y exactitud ± 1°C; humedad relativa
minado que es un proceso difícil de describir, debido a la con exactitud ±1° %), una tajadora empleada para asegurar
complejidad de los fenómenos internos y externos que o- la homogeneidad del espesor de las muestras estudiadas y
curren durante él (Sandoval et al., 2006). No obstante, en la moldes de acero inoxidable circular con radio de 3,19 cm y
literatura se encuentran varios métodos propuestos para a- fiete de 5 cm de largo y semieje menor de 3 cm.
nalizar el secado en capa fina de productos alimenticios: mo- Determinación del contenido de humedad inicial
delos teóricos, semiteóricos y empíricos. Los modelos de se-
cado basados en la teoría de difusión de líquidos han sido Previo a los ensayos de secado, se determinó el peso seco de
preferidos por los investigadores en el área de secado de ali- la muestra mediante el método de estufa a 105 °C durante
mentos (Chaves et al., 2003). Entre estos se encuentran los 24 horas (AOAC, 1997), con el valor obtenido se calculó el
modelos de Newton, Page, Page modificado, Henderson y contenido de humedad inicial para ambas variedades.
Pabis modificado, Midilli et al. (Ertekin y Yaldiz, 2004; Vega y
Lemus, 2006). Estos modelos de secado permiten no sólo
Construcción de las curvas de secado
predecir el mejor proceso, sino que también ofrecen herra- Se evaluó el proceso de secado para Dioscorea rotundata en
mientas para predecir las condiciones de almacenamiento y las variedades 9811-089 y 9811-091, a temperatura de bul-
empaque; además, ayudan a establecer el contenido final de bo seco de 45, 55 y 70°C, velocidad de aire de secado de
humedad de los productos agrícolas y los requisitos del pro- 1m/s, y geometría de rodajas y láminas. Inicialmente se pro-
ceso de secado (Vega y Fito, 2005). cedió a tomar las dimensiones y peso de la bandeja, luego se
ajustaron las condiciones de trabajo del equipo (temperatura,
La aplicación del secado como una operación que garantice
humedad, velocidad de aire de secado) y se procedió a co-
la calidad del producto deshidratado y que además sea via-
locar la muestra previamente adecuada y pesada en la ban-
ble desde el punto de vista económico requiere de un es-
deja, verificando que esta abarcara por completo la base de
tudio previo que permita determinar las condiciones y pará-
la misma. Durante el proceso de secado las muestras fueron
metros a tener en cuenta para garantizarlo. En este trabajo se
pesadas a intervalos de 5 minutos durante la primera hora y
realizaron y modelaron las curvas de secado de Dioscorea ro-
luego, cada 15 minutos hasta obtener peso constante; este
tundata variedades 9811-089 y 9811-091.
procedimiento se realizó de igual forma en cada tratamiento.
Para la construcción de las curvas de secado se utilizó el mé-

46 REVISTA INGENIERÍA E INVESTIGACIÓN VOL. 28 No. 2, AGOSTO DE 2008 (45-52)

 MONTES, TORRES, ANDRADE, PÉREZ, MARIMON, MEZA

todo descrito por Geankoplis (1998). Con el fin de determi- el tiempo de secado y las constantes de los modelos evalua-
nar el modelo que mejor describe el comportamiento cinéti- dos. Todos los datos fueron procesados y evaluados a través
co de secado, se evaluó el ajuste de los datos experimentales del software Statistica 7.0 (StatSoft, Inc., 2004).
a 7 modelos matemáticos (Tabla 1), escogiendo el mejor mo-
delo que se ajusta a los datos experimentalmente. Resultados y discusión
Tabla 1. Modelos matemáticos para describir el comportamiento ciné- Curvas de rapidez de secado
tico de secado. Fuente: Ertekin y Yaldiz, 2004.
En las Figuras 1a, 1b, 2a y 2b se observa la in-
Nombre Ecuación fluencia de la temperatura del aire de secado so-
Newton MR = Exp (− kt ) bre el proceso de secado de Dioscorea rotundata
para las variedades 9811-089 y 9811-091 con
Page MR = Exp(− kt n ) geometrías de rodajas y láminas. Al aumentar la
Page Modificado [
MR = Exp (− kt )
n
] temperatura del aire de secado disminuye signifi-
cativamente el tiempo de secado, necesitando
Henderson y Pabis MR = a[Exp(− kt )] menos tiempo para finalizar el secado con la
temperatura de 70 °C (cerca de 28,15% menos)
Logarítmico MR = a[Exp (− kt )] + c que para las temperaturas de 55 y 45 °C; resul-
tados similares obtuvieron Park et al. (2002);
Thomson t = a ln (MR ) + b[ln (MR )]
2
MR = 1 + at + bt 2

Maldonado y Pacheco (2003); Vega y Fito

Difusión MR = a[Exp(− kt )] + (1 − a )[Exp(− kbt )] (2005); Simal et al., (2005); Doymaz (2005);
Mohapatra y Rao (2005); Vega y Lemus (2006);
Ocampo (2006) y Falade et al., (2006), trabajan-
Donde MR es la razón de humedad, t es el tiempo de secado do con otras especies, vegetales y cereales. También se ob-
-1
en horas, K, K0,K1 son las constantes de secado (h ), a,b,c,n serva que aproximadamente en las cuatro primeras horas de
los coeficientes de ajuste. secado (dependiendo de la geometría y variedad), que la hu-
Se utilizó el método de regresión no lineal para estimar las medad libre decrece linealmente (R2 > 0.90), esto es, dX/dt
constantes de los modelos. Los criterios usados para escoger es una constante, al igual que la velocidad de secado para
el modelo de mejor ajuste fueron el coeficiente de correla- este mismo periodo de tiempo. De igual forma, después de
ción, error medio estimado (Ecuación1) y desvío medio rela- la cuarta hora de secado la velocidad disminuye con mayor
tivo (Ecuación2). rapidez hasta llegar a su humedad de equilibrio. Esta etapa se
caracteriza porque en ella ocurre la mayor parte del secado,
⎛ ∑ (Vexp − V p )2 ⎞ pues la humedad se difunde con menor velocidad a través
SE = ⎜ ⎟ (1)
⎜ ⎟ del sólido. Por otra parte, se pudo observar que el tiempo
⎝ GRL ⎠ final de secado, a temperaturas de secado, fue aproximada-
mente igual para ambas variedades de Dioscorea rotundata,
100 n Vexp − V p (2) y mostró ser siempre menor para la geometría de láminas,
P= ∑
n i =1 Vexp resultado que concuerda con el obtenido por Fioreze y Mori-
ni (2000) en el secado de Dioscorea sp. con diferentes cortes
donde P es el desvío medio relativo (%), SE es el error medio y temperaturas.
estimado, Vexp es el valor obtenido experimentalmente, Vp es
el valor obtenido por el modelo, n es el número de datos El análisis de varianza, Anova (Tabla 2), evidencia una
experimentales. influencia altamente significativa de la temperatura sobre el
tiempo de secado de D. rotundata, y de igual manera, una
Para la construcción de las curvas de secado se empleó un influencia significativa de la geometría, en tanto que la
diseño completamente al azar en arreglo factorial con tres variedad no presentó influencia significativa.
factores; la temperatura en tres niveles (45, 55 y 70 °C), la
Tabla 2. Análisis de varianza (Anova) para el tiempo de secado (t) de
geometría en dos niveles (rodajas y láminas) y la variedad en
D. Rotundata
dos niveles (9811-089 y 9811-091), realizando tres repeticio- Fuente DF Tipo I SS Media F-Valor Pr>F
nes por tratamiento. Los resultados experimentales fueron so- Temperatura 2 324950.00 162475 302.28 <0.0001
metidos al test de Shapiro Wilk, mediante el cual se verificó Variedad 1 1406.25 1406.25 2.62 0.1188
la normalidad de los mismos. Se procedió a realizar un análi- Geometría 1 18906.25 18906.25 35.17 <0.0001
sis de varianza (Anova) al 95% de confiabilidad, usando un Temperatura*
2 1950.00 975.00 1.81 0.1847
modelo de tres vías de clasificación (variedad, temperatura y Variedad
Temperatura*
geometría) con interacciones. Para la comparación de me- 2 350.00 175.00 0.33 0.7252
geometría
dias se usó el test de Duncan con un nivel de significancia Variedad*
del 5%; de igual forma, el análisis de varianza permitió medir 1 306.25 306.25 0.57 0.4577
gesometría
la incidencia de la temperatura, geometría y variedad sobre Temp.*Var*Geom 2 200.00 100.00 0.19 0.8314

REVISTA INGENIERÍA E INVESTIGACIÓN VOL. 28 No. 2, AGOSTO DE 2008 (45-52) 47

MODELADO DE LA CINÉTICA DE SECADO DE ÑAME (DIOSCOREA ROTUNDATA) EN CAPA DELGADA

2,50
2,5
45ºC 45ºC
2,00 55ºC 55ºC
Humedad Libre (B.S)
2

Humedad Libre (b.s)

70ºC 70ºC
1,50
1,5

1,00
1

0,50
0,5

0,00
0
0 200 400 600 800
0 200 400 600 800
Tiempo de secado (min)
Tiempo de secado (min)
a. a.

2,50 2,50
45ºC
45ºC
55ºC
2,00 55ºC 2,00
70ºC
Humedad Libre (b.s)

Humedad Libre (b.s)

70ºC
1,50 1,50

1,00 1,00

0,50
0,50

0,00
0,00
0 200 400 600 800
0 200 400 600 800
Tiempo de Secado (min)
Tiempo de Secado (min)
b. b.

Figura 1. Curvas de rapidez de secado de D. rotundata 9811-089 a Figura 2. Curvas de rapidez de secado de D. rotundata 9811-091 a
diferentes temperaturas. a) Rodajas y b) Láminas diferentes temperaturas. a) Rodajas y b) Láminas.

Curvas de velocidad de secado   1,2

En términos generales, al aumentar la temperatura del aire 1,0 45ºC

Velocidad de secado

de secado aumenta la velocidad de este (Figuras 3a, 3b, 4a y 55ºC

0,8
4b). Por lo tanto, en función de los resultados obtenidos se 70ºC

deduce que la cinética de secado de D. rotundata es depen- 0,6

diente de la temperatura del aire de secado. Se observa que 0,4

en el secado de D. rotundata a la temperatura de 70 °C para
0,2
la variedad 9811-089, se genera una velocidad inicial pro-
medio de secado del orden de 0,8045 Kg de agua/h x m2 0,0
para la geometría de rodajas, que resulta significativamente 0,0 0,5 1,0 1,5 2,0 2,5

superior a la velocidad de secado alcanzada a una tempera- Hum edad Libre(B.S)

a.
tura 45 °C que es del orden de 0,3715 Kg de agua/h x m2 pa-
1,2
ra esta misma variedad y geometría. De igual forma, se logra
45ºC
Velocidad de secado

observar cómo los valores de velocidad de secado para las 1,0

55ºC
dos variedades de D. rotundata en estudio son similares para 0,8
70ºC
las mismas geometrías y temperaturas de aire de secado. Por 0,6
otra parte, se observan diferencias entre la velocidad de se- 0,4
cado de rodajas y láminas, lo cual se explica por que esta úl- 0,2
tima geometría posee mayor área de intercambio de energía 0,0
y masa con el aire de secado, resultado que coincide con 5 0,0 0,5 1,0 1,5 2,0 2,5
Fioreze y Morini (2000) en el secado de Dioscorea sp. con Hum edad Libre(B.S)
diferentes cortes y temperaturas. b.

Figura 3. Curvas de velocidad de secado de D. rotundata 9811-089,

para tres temperaturas. a) Rodajas y b) Filetes.

Además se distinguen tres etapas de secado; en la primera se

observa un corto periodo de inducción, en la que el produc-

48 REVISTA INGENIERÍA E INVESTIGACIÓN VOL. 28 No. 2, AGOSTO DE 2008 (45-52)

 MONTES, TORRES, ANDRADE, PÉREZ, MARIMON, MEZA

to se encuentra a una temperatura inferior a la que tendrá en donde MR es la razón de humedad, X es la humedad en base
el proceso, y la velocidad de secado va en aumento. Este seca (Kg agua/Kg s.s.) a cada tiempo, Xe es la humedad de e-
periodo es llamado también creciente de velocidad (Sando- quilibrio (Kg agua/Kg s.s.), Xo es la humedad inicial (Kg agua
val et al., 2006). En la segunda etapa se observa un periodo /Kg s.s.), t es el tiempo en hora, l es el semiespesor de la
de velocidad constante, en la que la temperatura del alimen- lámina (m) y Df es el coeficiente de difusión expresado en m2
to ha alcanzado su valor de equilibrio, la superficie del pro- s-1. El coeficiente difusional (Df) para cada temperatura se cal-
ducto está muy mojada al principio debido a que se encuen- culó a partir de la pendiente de la línea recta obtenida al gra-
tra cubierta de un mucílago característico de la especie r- ficar el Ln (MR) en función del tiempo de secado en horas
otundata; sobre este se forma una película continua de agua, (ver Figura 5).
esta capa es agua no combinada y el secado se da como si el
sólido no estuviera presente. La tercera etapa corresponde a Los valores de difusividad aumentaron con el incremento de
un periodo de velocidad decreciente en la que se lleva a la temperatura del aire de secado (Tabla 3), lo que muestra
una disminución de las resistencias internas del sólido con el
cabo la mayor parte del proceso, esto demuestra que la di-
acrecentamiento de la temperatura de secado. Una explica-
fusión es el mecanismo físico dominante que gobierna el mo-
ción de esa adición: puede ser debida al menor encogimien-
vimiento de la humedad en el secado de las muestras de
to observado en las temperaturas más altas; menor encogi-
ñame. Similares resultados fueron reportados para judías (Ro-
miento implica una mayor porosidad y consecuentemente u-
sello et al., 1997), okra (Doymaz, 2004; Gogus y Maskan,
na mayor área de transferencia de masa. Resultados similares
1999), berenjena (Ertekin y Yaldiz, 2004) y D. rotundata y D.
se obtuvieron en piña (Park et al., 2002); papaya chilena (Ve-
alata (Falde et al., 2006).
ga y Lemus, 2006); ñame (Falade et al., 2006) y en mango
(Ocampo, 2006). Por otra parte, el análisis de varianza (Tabla
1,2 4), indica que existe influencia significativa de la temperatura
45ºC
del aire de secado, y el análisis de comparación de medias
Velocidad de secado

1,0
55ºC de Duncan evidenció diferencia estadística significativa para
0,8
70ºC las difusividades calculadas para las tres temperaturas (45, 55
0,6
y 70 °C).
0,4

0,2 Tiempo de Secado (h)

0,0 0,00 5,00 10,00 15,00
0,0 0,5 1,0 1,5 2,0 2,5 0,5
Hum edad Libre(B.S)
(a )
0 45ºC
55ºC
Ln (X - X e)/(X o - X e)

1,2 -0,5
45º C 70ºC
Velocidad de secado

1,0 -1 -- calcu
55º C
0,8
70º C -1,5
0,6

0,4 -2

0,2 -2,5
0,0
-3
0,0 0,5 1,0 1,5 2,0 2,5
Hum edad Libre(B.S)
(b) -3,5

Figura 4. a) Curvas de velocidad de secado de D. rotundata 9811- Figura 5. Correlación logarítmica para las curvas de secado a 45, 55
091, para tres temperaturas. a) Rodajas y b) Filetes y 70°C de D. rotundata 9811-089. Geometría de rodajas

Coeficiente difusivo (Df) y energía de activación (Ea) En general, la influencia de la temperatura de secado sobre la
difusividad del agua en D. rotundata sigue una tendencia ti-
Para modelar el periodo de velocidad decreciente en el seca-
po Arrhenius
do de D. rotundata se utilizó el modelo difusional, método
que ha sido utilizado en el estudio de la transferencia de ma- ⎛ E ⎞
D f = D0 .Exp⎜ − a ⎟ (4)
sa de productos agrícolas (Vega y Fito, 2005). Se utilizó la e-
cuación integrada de la segunda ley de Fick (Ecuación 3) pa- ⎝ RT ⎠
ra tiempos largos y geometría plana en una dimensión. donde Do es factor de Arrhenius, Ea es la energía de activa-
ción, R es la constante universal de los gases (8,314KJ mol-1) y
X − Xe 8 ⎛ − Df π t ⎞ 2
(3)
MR = = 2 exp⎜ ⎟ T es la temperatura absoluta.
X0 − Xe π ⎜ 4l 2 ⎟
⎝ ⎠

REVISTA INGENIERÍA E INVESTIGACIÓN VOL. 28 No. 2, AGOSTO DE 2008 (45-52) 49

MODELADO DE LA CINÉTICA DE SECADO DE ÑAME (DIOSCOREA ROTUNDATA) EN CAPA DELGADA

Tabla 3. Difusividad efectiva (Df), energía de activación (Ea) y factor Modelación matemática
de Arrhenius (D0), para las variedades de D. rotundata 9811-089 y
9811-091 Los modelos de Newton, Page, Page modificado, Hen-
derson y Pabis, logarítmico, difusional y Thomson des-
Difusividad efectiva (m2/s) criben satisfactoriamente las curvas de secado de Dios-
45ºC 55ºC 70ºC Ea(Kj/mol) D0(m2/s) corea rotundata a las temperaturas estudiadas, una vez
9811-089 Rodajas 2.67x10 -06
3.91x10 -06
4.91x10 -06
21.49 9.55x10-3 que sus valores de coeficiente de determinación (R2)
9811-089 Filetes 2.77x10 -06
3.99x10 -06
4.71x10 -06
18.62 3.45x10-3 fueron superiores a 90% y los errores medios relativos
9811-091 Rodajas 2.82x10-06 3.72x10-06 4.63x10-06 17.66 2.31x10-3 (P) inferiores al 10%. Sin embargo, el modelo que me-
9811-091 Filetes 2.98x10 -06
3.45x10 -06
4.81x10 -06
17.50 2.19x10-3 jor describe el comportamiento de los datos de secado
de Dioscorea rotundata variedades 9811-089 y 9811-
091 en geometrías de rodajas y láminas para las tem-
Tabla 4. Análisis de varianza para el coeficiente difusivo de la hume- peraturas de 45, 55 y 70 °C, es el modelo logarítmico
dad (Df)
(Tabla 5).

Fuente DF Tipo I SS Media F-Valor Pr>F En la Tabla 6 se puede observar como el parámetro
Temperatura 2 0.1784 0.0892 85.89 <0.0001 K del modelo logarítmico aumenta su valor con el
Variedad 1 0.0075 0.0075 7.25 0.0127 incremento de la temperatura, coincidiendo este
Geometría 1 0.01087 0.0108 10.47 0.0035 resultado con el obtenido por Ertekin y Yaldiz,
Temperatura*Variedad 2 0.0066 0.0033 3.20 0.0585 2004; Correa et al. 2004; Babilis y Belessiotis,
Temperatura*geometría 2 0.00775 0.00387 3.74 0.0387 2004; Vega y Fito, 2005; Vega y Lemus, 2006;
Variedad*geometría 1 0.000025 0.0000253 0.02 0.8772 Falade et al., 2006, para el secado de berenjena,
Temp.*Var*Geom 2 0.00134 0.000669 0.64 0.5339 cebolla, higo, pimiento rojo, papaya chilena y ñame
(variedades D. alata y D. rotundata), respectivamen-
te. De igual forma, los parámetros a y b varían para
Al representar Ln (Df) en función de 1/T se obtiene una línea las geometrías, siendo estos resultados confirmados con el a-
recta (Figura 6) de cuya pendiente se obtiene Ea, y de la or- nálisis de varianza (Anova) para estas constantes a un nivel de
denada en el origen, el factor de Arrhenius (D0). confiabilidad de 95%.

-12 Conclusiones
2,90E-03 3,00E-03 3,10E-03 3,20E-03
-12,1
Expc
La temperatura de aire de secado tuvo una gran influencia
-12,2 Lineal (Expc) en el tiempo de secado, evidenciándose menores tiempos
-12,3 con cada incremento de esta; de la misma manera, la geo-
-12,4
metría de láminas permitió obtener menor tiempo de secado
Ln Df

para ambas variedades en comparación a las rodajas, en tan-

-12,5
to las dos variedades de Discorea rotundata presentaron
-12,6 tiempos de secado estadísticamente iguales bajo un nivel de
-12,7 significancia del 5%.
-12,8 Las curvas de secado permitieron identificar claramente un
-12,9
LN(Df) = -2585,4/T - 4,6515
2
periodo de inducción, uno de velocidad de secado constante
R = 0,9442
y otro decreciente para cada una de las variedades, tempera-
-13
turas y geometrías estudiadas.
1/T

Figura 6. Correlación exponencial de la difusividad (Df) a 45, 55 y 70 La difusión es el mecanismo físico que gobierna la migración
°C de D. rotundata 9811-089, geometría de rodajas de humedad de la D. rotundata por cuanto los datos expe-
rimentales se ajustaron correctamente a la ecuación de Fick.
La energía de activación (Tabla 3) presenta una magnitud
que se encuentra dentro del rango de valores encontrados en Los modelos cinéticos de Page, Thompsom, Newton, Page
la literatura para materiales alimenticios (12,7 a 110 KJ/mol): modificado, Henderson y Pabis, difusional y logarítmico, se
28,39 KJ/mol para zanahoria (Doymaz, 2004) y 20KJ/mol pa- ajustan bien a los datos de secado bajo las condiciones es-
ra papa (Rosello et al.,1997). Estas magnitudes de energía de tudiadas; no obstante, es este último el que mejor describe el
activación representan la energía requerida para iniciar la proceso, una vez que presenta los valores más altos de coe-
difusión de la humedad del D. rotundata durante el secado y ficiente de determinación (R2) y los más bajos de error medio
dan un indicio de la influencia de la temperatura en el pro- estimado (SE) y desvío medio relativo (P%) para todas las
ceso (2001). temperaturas (45, 55 y 70 °C) y geometrías empleadas en es-
te estudio (rodajas y láminas).

50 REVISTA INGENIERÍA E INVESTIGACIÓN VOL. 28 No. 2, AGOSTO DE 2008 (45-52)

 MONTES, TORRES, ANDRADE, PÉREZ, MARIMON, MEZA

Tabla 5. Criterios estadísticos (R2 (%), P (%) y SE) de ajuste de modelos Doymaz, I., Pretreatment effect on sun drying of mulberry
para las curvas experimentales de secado de D. rotundata 9811-089 Fruits (Morusalba L)., Journal of food engineering, Vol.
y 9811-091 y dos geometrías a 45°C 65, No. 2, 2004, pp. 205-209.
Modelo Doymaz, I., Drying characteristics and
Geometría

kinetics of Okra., Journal of food engi-

Variedad

Estadístico
Page Henderson neering, Vol. 69, No. 3, 2005, pp.
Newton Page Logarítmico Difusional Thomson 275-279.
modificado y Pabis
2
R (%) 98.54 98.71 98.54 98.55 99.17 99.08 98.36
Ertekin, C., Yaldiz, O., Draying of e-
ggplant and selection of a suitable thin
Rodajas

SE 0.0381 0.0358 0.0381 0.0380 0.0289 0.0303 0.0404

layer drying model., Journal of food
9811-089

P (%) 32.608 30.711 32.463 32.735 23.556 24.943 31.644 engineering, Vol. 63, No. 3, 2004, pp.
2
R (%) 98.28 99.51 98.28 98.73 99.58 99.40 99.65 349-359.
Filetes

SE 0.0444 0.0236 0.0444 0.0383 0.0221 0.0264 0.0201 Falade, K., Olurin, T., Ike, E., Aworh,
P (%) 39.101 19.824 39.101 33.220 17.791 34.856 15.684 O., Effec of pre-treatment and tempe-
2
R (%) 98.33 99.04 98.33 98.65 99.19 99.15 98.86 rature on air – Drying of Discorea alata
Rodajas

SE 0.0428 0.0325 0.0429 0.0385 0.0299 0.0307 0.0354 and Discorea rotundata slices., Journal
of food engineering, Vol. 80, No. 4,
9811-091

P (%) 37.663 23.193 0.0377 33.671 23.529 23.865 27.239

2
2007, pp. 1002-1010.
R (%) 98.63 99.80 98.63 99.10 99.86 99.82 99.82
Fioreze, R., Morini, B. 2000. Yam
Filetes

SE 0.0395 0.0150 0.0395 0.0321 0.0125 0.0145 0.0142

(Dioscorea sp.) drying with different
P (%) 34.943 13.169 34.943 27.993 0.9723 11.313 11.347
cuts and temperatures: experimental
and simulated results. Ciênc. Tecnol.
Aliment. Vol. 20, No. 2, 2000, pp.
Tabla 6. Constantes del modelo logarítmico para las dos variedades,
262-266.
dos geometrías y tres temperaturas.
Geankoplis, C., Procesos de trans-
Variedad 9811-089 Variedad 9811-091 porte y operaciones unitarias., 3a.
Constantes

edición. México: Compañía Edi-

Rodajas Laminas Rodajas Laminas torial Continen-tal S. A., 1998, pp.
45ºC 55ºC 70ºC 45ºC 55ºC 70ºC 45ºC 55ºC 70ºC 45ºC 55ºC 70ºC 1022.
K 0.1815 0.1883 0.272 0.185 0.269 0.306 0.177 0.245 0.310 0.183 0.214 0.329 Gögus, F., Maskan, M., Water
a 10.984 11.989 10.780 11.520 11.135 10.735 11.230 10.822 10.645 11.504 11.285 10.727 adsorption and drying characteris-
b -0.1205 -0.176 -0.094 -0.127 -0.078 -0.075 -0.102 -0.067 -0.055 -0.123 -0.107 -0.066 tics of Okra (Hibis cus esculentas
L)., Drying Tech-nology. Vol. 17,
No. 4, 1999, pp. 883-894.
Bibliografía Guzmán, M., Buitrago, G., ÑAME: Producción de Semilla
por Biotecnología., Santa fe de Bogotá: Unibiblos, 2000,
AOAC., Assocciation of Oficial Analytical Chemist. Oficial pp. 17-20.
methods of análisis., 16th edition. Washington, 1997, pp. Maldoado, R., Pacheco, T., Curvas de deshidratación del
3005. Brócoli (Brassica Oleracea L. var Italica Plenk) y Coliflor
Babilis, S., Belessiotis, V. Influence of the drying constants (Brassica oleraceae L. var Botritis)., Revista de la Facultad
and moisture diffusivity during the thin-layer drying of figs., de Agronomia, Universidad Central de Venezuela. Vol.
Journal of food Engineering, Vol. 65, No. 3, 2004, pp. 20, 2003, pp. 306-319.
449-458. Mohapatra, D., Rao, P., A thin layer drying model of
Cardoso, F., Palmeira, J., Rodrigues, M., Honorato, F., parboiled wheat., Journal of food engineering. Vol. 66,
ALVES, S. Comparación de modelos matemáticos de calor No. 4, 2005, pp. 513-518.
isotérmico de desorción en pulpa de Guayaba., Revista Nicoleti, F., Cinética de secado y encogimiento de rodajas
Brasileira de almacenamiento, Vol. 29, No. 1, 2004, pp. de piña, Universidad estatal paulista. Campus de San José
28-34. de Río prieto-Brasil. 2001
Chaves, M., Sgroppo, S., Avanza, J., Cinética de Secado de Ocampo, A., Modelo cinético del secado de la Pulpa de
Berenjena (Solanum melongena L.)., Universidad nacional mango., Revista Escuela de Ingeniería de Antioquia EIA,
del nordeste. Comunicaciones científicas y tecnológicas. Vol. 5, 2006, pp. 119-128.
Corrientes–Argentina. Facultad de ciencias exactas y Park, K., Vohnikova, Z., Reisbros, F., Evaluation of driying
agrimensura, 2003, [Citado Oct-06-2006], Resumen E- parameters and desorption isotherms of garden mint
60, Disponible en Internet: http:www.exa.unne.edu.ar. leaves (Mentha crispa L.)., Journal foo Engineering, Vol.
Correa, P., Da Silva, P., Almeida, L., Estudo das 51, 2002, pp. 193-199.
propriedades físicas e de transporte na secagem de CE- Rincón, A., Evaluación del posible uso tecnológico de al-
bola (Allium cepa L.) em camada delgada., Ciência y gunos tubérculos de las dioscoreas: Ñame congo (Diosco-
tecnología de los alimentos, Vol. 24, No. 3, 2004. rea bulbífera) y mapuey (Dioscorea trífida)., Archivos

REVISTA INGENIERÍA E INVESTIGACIÓN VOL. 28 No. 2, AGOSTO DE 2008 (45-52) 51

MODELADO DE LA CINÉTICA DE SECADO DE ÑAME (DIOSCOREA ROTUNDATA) EN CAPA DELGADA

Latinoamericanos de Nutrición. Vol. 50, No. 3, 2000, pp. ñame (D. alata 9506-021, 9506-027), para la obtención
286-290. del mayor contenido de almidón., Trabajo de grado para
Romero, L., Kieckbush, T., Influência de Condições de optar al titulo de Ingeniero de alimentos, Universidad de
Secagem na Qualidade de Fatias de Tomate., Brazilian Córdoba. Montería – Córdoba, 2007.
Journal of food Technology, Vol. 6, No. 1, 2003, pp. 69- Simal, S., Femenia, M., Rosello, C., Use of exponential,
76. Page an diffusional models to simulate the drying kinetics
Rosello, C., Simal, S., San Juan, N., Nonostropic mass of Kiwi fruit., Journal of food Engineering. Vol. 43, 2005,
Transfer model for green Bean Drying., Journal of pp. 109-114.
Agricultural and food Chemistry. Vol. 45, 1997, pp. 337- StatSoft, Inc., STATISTICA (data analysis software system).,
342. 7.0, 2004, www.statsoft.com.
Sandoval, S., Rodriguez, J., Méndez, L., Sánchez, J., Rapi- Vega, A., Fito, P., Modelado de la cinética de secado del
dez de secado reducida: Una aplicación al secado con- pimiento Rojo (capsicum annuum L. cv Lamuyo)., In-
vectivo de plátano Roatán., Revista Mexicana de formación tecnológica, Vol. 16, No. 6, 2005, pp. 3-11.
Ingeniería química, Vol. 5, No. 1, 2006, pp. 35-38. Vega, A., Lemus, R., Modelado de la cinética de secado de
Sierra, O., Munive, J., Determinación de las mejores con- papaya chilena (basconcellea pubescens)., Información
diciones del proceso de secado de dos variedades de tecnológica. Vol. 17, N° 3, 2006, pp. 23-31.

52 REVISTA INGENIERÍA E INVESTIGACIÓN VOL. 28 No. 2, AGOSTO DE 2008 (45-52)