Spanish (pdf)
Article in xml format
Article references
Automatic translation
Send this article by e-mail
Cited by SciELO 
Cited by Google 
Similars in
SciELO 
Similars in Google 
Permalink
INTRODUCCIÓN
Uno de los problemas que presentan los caficultores para la comercialización del café pergamino seco es garantizar el contenido de humedad óptimo entre 10 y 12 % base húmeda, lo cual conlleva problemas de calidad y pérdidas al caficultor [1]. Las estructuras de secado solar permiten aumentar la temperatura y reducir la humedad relativa del aire con respecto a las condiciones ambientales, lo cual asegura la viabilidad del dispositivo para el secado de productos [2].
Varios estudios del comportamiento de la temperatura y velocidad del viento en sistemas de secado solar con aplicación en productos agroalimentarios han utilizado dinámica de fluidos computacional (CFD) con resultados numéricos comparables a los datos experimentales [3]. Algunas investigaciones realizaron la optimización de estructuras de secado con modificaciones geométricas [5-6]. Los modelos de elementos finitos se pueden usar para simular el secado de café pergamino y para un diseño óptimo de los secadores de café [7]. Por tanto, la simulación de secadores solares con CFD es útil para el diseño, en pos de garantizar la calidad en el secado y aprovechar al máximo la energía solar, dicha tecnología presenta un uso incipiente en el sector agroalimentario de Colombia.
Con el secador solar en estudio, se presentan tiempos de secado que se pueden prolongar hasta nueve días durante la época de lluvias, comportamiento presente en temporadas con mayores niveles de producción en la zona cafetera colombiana, o que vuelve ineficiente este tipo de sistemas de secado. Por tanto, el objetivo del presente estudio es simular térmicamente un secador solar de café con la herramienta CFD, mediante la modificación de la configuración geométrica en la apertura de ventanas.
1. MATERIALES Y MÉTODOS
1.1 Descripción del secador solar experimental tipo túnel
El secador solar consta de una estructura en madera, una cubierta plástica y ventanas frontales con cortinas plegables que permiten cambiar la apertura en la entrada y salida del aire. Se encuentra ubicado en el municipio de Fredonia, departamento de Antioquia, localizado en la coordenada 5°52'9,68" N; 75°40'6,22" O. Las condiciones climaticas medias del sitio corresponden a las presentes en las zonas cafeteras de Colombia con temperatura de 20,2 °C, humedad relativa de 74,0 %, y precipitación de 2.674 mm/ año [7]. El secador opera con ventilación natural, tiene un área de secado de 3,85 m x 9,00 m que permite una capacidad para 409 kg de café pergamino seco con capa de 3 cm y una cubierta plástica con forma curva tipo túnel, que permite la incidencia de la radiación solar durante las horas de luz del día, geometría típica en secadores utilizados por la caficultura colombiana (figura 1 y figura 2). El plástico utilizado fue polietileno PE calibre 6 con una emisividad de 0,1.
1.2 Adquisición de datos experimentales
La adquisición de datos se realizó con equipos marca Campbell Scientific, mediante la conexión de sensores de temperatura, humedad relativa, velocidad y dirección del viento al datalogger CR1000 [8]. Se instaló un anemómetro WindSonic4 en la ventana a la entrada del aire que permitió medir la velocidad con una resolución de 0,01 m/s en un rango de operación entre 0 y 60 m/s, la dirección del viento se midió con una resolución de 1°. Para la temperatura y humedad relativa del aire en la entrada y salida del aire se utilizó el sensor CS215, para medir temperaturas entre 5 y 40 °C con una precisión de ± 0,4 °C y humedades relativas entre 0 y 100 % con una precisión de ± 4 %. La temperatura del plástico y del café se midieron con el sensor de temperatura infrarroja IRR-P con precisión de ± 0,2 °C entre -10 y 65 °C. La tasa de agua evaporada desde la capa de café se determinó con el método Gravimet [10], mediante una balanza de precisión modelo PL150E con capacidad maxima para 1.520 g y precisión de ± 0,007 g.
Las condiciones climáticas externas al secador solar se midieron con una estación climática marca Decagon Device serie Em50G, ubicada a una distancia de 10 m y a una altura de 2 m, la medición de temperatura tiene una precisión de ± 0,5 °C entre 5 y 40 °C, para la humedad relativa la precisión es de ± 2 % entre 90 y 20 %.
1.3 Dominio computacional
El modelamiento se realizó con una geometría y malla tridimensional no estructurada, se refinaron nodos en los primeros 10 cm de las zonas adyacentes a la cubierta plástica y capa de café. Para garantizar calidad en los resultados se realizó una prueba de malla con el programa Ansys ICEM® para evaluar el comportamiento térmico del secador con 8.963, 19.328, 60.975 y 121.276 nodos. Los resultados del perfil de velocidad en el centro del secador obtenidos con 60.975 y 121.276 nodos presentaron una tendencia similar, mientras que con 8.963 y 19.328 nodos se presentó mayor dispersión de los resultados (figura 2). No obstante, la malla de 121.276 nodos incrementó considerablemente el tiempo de cómputo, por tanto, se seleccionó la malla generada con 60.975 nodos por presentar la precisión y el costo computacional esperado (figura 3). La computadora utilizada fue una Workstation con procesador Intel Xeon, memoria RAM de 16 GB, disco duro de 1 TB y tarjeta de video Nvidia.
El análisis se elaboró con las mediciones experimentales realizadas el 25 de abril del 2016, día en el cual predominaron condiciones climáticas desfavorables en el proceso de secado. En la tabla 1 se registran los valores de entrada constantes obtenidos de manera experimental y utilizados en las simulaciones, se operó el secador con la configuración geométrica propuesta T1 (figura 3), correspondiente a la operación habitual del secador. Los valores de temperatura, humedad relativa y velocidad fueron las propiedades del aire en la entrada del secador, mientras el flujo másico fue el agua evaporada desde la capa de café. Todas las propuestas fueron simuladas en estado estacionario a las 2:00 p.m., cuando se presentaron las máximas tasas evaporativas de agua, a las 8:00 a.m. y 6:00 p.m., se realizó la apertura y el cierre de ventanas en la operación habitual del secador, respectivamente.
1.4 Configuración del modelo computacional
La simulación se realizó con el paquete comercial Ansys CFX® versión 14 para resolver las ecuaciones de Navier-Stokes que discretizan el dominio del fluido con el método de volúmenes finitos [11]. Las ecuaciones de masa, momentum y energía también se conocen como las ecuaciones de Navier-Stokes y se expresan en términos generales de la forma en la que se presenta en las ecuaciones 1, 2 y 3.
Donde: Cp es el calor específico; [J kg-1 °C-1]. k es la conductividad térmica; [W m-1 °C-1]. p es la presión; [N m2]. U es el vector velocidad; [m s1]. p es la densidad; [kg m-3]. g es la gravedad; [m s 2]. µ T es la viscosidad dinámica; [kg m-1 s1]. Tes la temperatura; [°C]. S h incluye la fuente de calor debido a la incidencia de la radiación solar.
La simulación del secador solar se consideró con el modelo de turbulencia Reynolds Stress BSL por presentar mayor precisión respecto al modelo k - ε. El modelo BSL permite analizar la dinámica de fluidos alrededor de las paredes según el número de Reynolds del aire, al agregar una tensión adicional a la viscosidad dinámica µ T [11].
1.5 Condiciones de frontera
Las condiciones de frontera en el modelo CFD se definieron con los siguientes criterios:
La humedad relativa se definió como una variable adicional al utilizar las ecuaciones de Tetens 4 y 5 [12].
Donde: e v es la presión de vapor de agua en la atmósfera; [kPa]. e s es la presión de vapor de agua a saturación; [kPa]. ω es la humedad absoluta del aire; [g vapor de agua/kg aire seco]. T es la temperatura del aire, [°C].
2. La transferencia de calor por radiación se consideró al acoplar el método de Monte Carlo con el término fuente S h de la ecuación de energía de Navier-Stokes [14].
3. La tasa evaporativa de agua desde la capa de café se determinó con el método Gravimet, mediante la ecuación 6.
Donde: m ev , es la tasa evaporativa de agua en la capa de café [kg/s)]; ∆m es el diferencia de masa de café humedo [kg]; ∆t es el intervalo de tiempo en la medición de la masa de café humedo [s].
4. La entrada de aire a través de la ventana y la tasa evaporativa de agua en la capa de café se consideraron uniformes espacialmente.
1.6. Propuestas geométricas en el secador solar
El tiempo de secado depende del diseño geométrico del secador solar [14]. La optimización de un secador para frutas mediante modificaciones geométricas con CFD permitió aumentar la temperatura del aire y lograr una distribución uniforme del flujo [5].
La figura 4 muestra las seis propuestas con diferentes configuraciones geométricas en la apertura de ventanas con ventilación natural que se evaluaron en el presente estudio. La geometría T1 corresponde a la configuración original del secador con ventanas en la entrada y salida del aire de igual dimensión y apertura del 100 %, las geometrías T2, T3 y T4 solo tienen modificaciones en la ventana a la salida del aire con aperturas triangulares del 55 %, 40 % y 25 % con relación a la apertura original de la geometría T1, respectivamente. La geometría T5 tiene la salida del aire a traves de una sección rectangular en la parte inferior, y la geometría T6 corresponde a la apertura máxima de ambas ventanas, al suprimir el material plástico de los costados frontales.
1.7 Análisis estadístico
Para analizar estadísticamente las diferencias entre las configuraciones geométricas propuestas respecto a la temperatura y la humedad relativa en un plano a 10 cm sobre la capa de café, se realizó una prueba de comparación de medias con una prueba T-student con un nivel de significancia del 5 %.
2. RESULTADOS Y DISCUSIÓN
2.1 Comportamiento térmico en el interior y exterior del secador con operación original T1
La viabilidad del secador solar se determina mediante la diferencia de temperatura y de humedad relativa con respecto a las condiciones ambientales [2]. Los datos experimentales medidos en el interior y exterior del secador solar, con la configuración geométrica T1, evidencian que la temperatura del aire en el interior fue superior a la temperatura ambiental: la diferencia de temperatura fue 1,9 °C, 3,6 °C y 2,5 °C a las 8:00 a.m., 2:00 p.m., y 6:00 p.m., respectivamente, mientras que la humedad relativa disminuyó 0,5 %, 10,4 % y 6,8 % a las 8:00 a.m., 2:00 p.m., y 6:00 p.m. Janjai et al. [16], obtuvieron un incremento de temperatura de 8 °C aproximadamente, en un secador solar tipo túnel con ventilación forzada ubicado en un ambiente tropical.
2.2 Análisis de temperatura y humedad relativa entre configuraciones geométricas
El análisis de temperatura y humedad relativa en las seis configuraciones geométricas propuestas se realizaron con los resultados de las simulaciones en un plano a 10 cm sobre la capa de café Z= 10 cm, al considerar que la tasa de evaporación de agua, el tiempo de secado y la calidad en taza dependen de las condiciones térmicas del aire adyacente al grano, este hecho ha sido ampliamente estudiado en secado mecánico [16].
En la tabla 2 se indica que hay diferencias significativas (p<0,05) en la temperatura y humedad relativa de la geometría T6 durante las horas objeto del estudio, lo cual evidencia que la máxima apertura de ventanas no generó aumento de la temperatura, y una disminución de la humedad relativa útil para el proceso de secado. No obstante, se presenta una diferencia significativa en la temperatura de la geometría T4 a las 2:00 p.m., y 6:00 p.m., que puede aumentar entre 3,6 °C y 4,2 °C en el plano Z= 10 cm con respecto a la temperatura ambiental. El análisis anterior logra predecir mayores temperaturas cuando la ventana a la salida del aire tiene una apertura del 25 % (geometría T4), mientras la humedad relativa no presentó un aumento que afecte el proceso, lo cual permite aumentar la tasa de evaporación y disminuir el tiempo de secado del café [2].
Tabla 2 Valores medios e intervalos de confianza para la temperatura y humedad relativa en el plano Z= 10 cm sobre la capa de secado a las 8:00 a.m., 2:00 p.m., y 6:00 p.m.
Nota: valores con la misma letra significa que no existe diferencia significativa (p<0.05)
Fuente: elaboración propia
Un análisis estadístico previo de la pérdida de peso indicó un secado uniforme de la capa de café con la geometría T1, lo cual evidencia que la distribución de las condiciones térmicas al interior del secado solar no genera diferencias significativas en la pérdida de peso entre diferentes posiciones con la geometría T1. Por consiguiente, al comparar la desviación estándar de la geometría T1 con las demás geometrías registradas en la Tabla 2, se deduce que el comportamiento térmico de las geometrías propuestas en los tratamientos no afecta negativamente la uniformidad del secado, toda vez que las desviaciones estándar con la geometría T1 son similares a las desviaciones de las geometrías T2, T3, T4, T5 y T6, y fueron calculadas con datos de temperatura y humedad relativa obtenidos en el proceso de simulación en 27 posiciones uniformemente distribuidas sobre el plano Z=10 cm.
En las figuras 5, 6 y 7 se observa la distribución espacial de temperatura y humedad relativa en el plano Z=10 cm a las 8:00 a.m., 2:00 p.m., y 6:00 p.m., respectivamente. En todas las geometrías e instantes de tiempo se da un aumento de temperatura durante la circulación del aire desde la entrada hacia la salida, la humedad relativa presenta un comportamiento similar producto de la evaporación de agua desde la capa de café.
Fuente: elaboración propia
Figura 5 Comportamiento de la temperatura y humedad relativa en el plano Z=10 cm a las 8:00 a.m.
Fuente: elaboración propia
Figura 6 Comportamiento de la temperatura y humedad relativa en el plano Z=10 cm a las 2:00 p.m.
Fuente: elaboración propia
Figura 7 Comportamiento de la temperatura y humedad relativa en el plano Z=10 cm a las 6:00 p.m.
La temperatura de la capa de café y del plástico fueron superiores a la temperatura del aire que ingresa al secador solar, lo que generó los mayores diferenciales térmicos en las zonas adyacentes a dichos materiales, producto de la transferencia de calor por convección.
La geometría T6 presenta mayor uniformidad de la temperatura y humedad relativa por tener apertura en todo el ancho de la ventana a la salida del aire, no obstante, el incremento de temperatura es mínima con respecto a la temperatura atmosférica por la alta tasa de ventilación interna o renovación del aire generada con la apertura total.
En los instantes de tiempo simulados, las máximas temperaturas en el plano Z=10 cm se presentaron en la geometría T4; la apertura del 25 % de la ventana en la salida del aire restringe la transferencia de calor hacia el exterior del secador solar, además, la radiación en onda larga o infrarroja reflejada entre la capa de café y el plástico calienta el aire interno e intensifica el efecto invernadero [14].
3. CONCLUSIONES
La dinámica de fluidos computacional CFD fue una herramienta versátil y confiable para la comparación térmica en diferentes configuraciones geométricas para la apertura de ventanas del secador solar.
Los resultados de la simulación muestran que la geometría T4 con apertura del 25 % en la ventana a la salida del aire fue adecuada para aumentar la temperatura de manera uniforme sobre la capa de café; no se recomienda implementar la geometría T6 con la apertura máxima de ambas ventanas en el secado solar por el poco incremento de temperatura, hasta de 0,3 °C en épocas lluviosas.
La apertura máxima de las ventanas en la salida y entrada del aire correspondiente a la geometría T6 permitió disminuir hasta 0,7 % la humedad relativa del aire con respecto a las demás geometrías, producto de la alta tasa de renovación del aire.
