INTRODUCCION
Colombia lidera actualmente la producción de ornamentales y flores de corte en América del Sur (Mena-Vásconez et al., 2016). La rosa (Rosa sp.) es una de las especies más cultivadas y comercializadas en el mercado internacional. Esta especie se encuentra establecida en invernaderos de plástico, que desde hace varias décadas se convirtieron en la opción tecnológica predominante para el cultivo protegido en climas templados y cálidos (McCartney et al., 2018). En Colombia existe una amplia gama de diseños y tamaños de invernaderos que han evolucionado a partir del invernadero tradicional colombiano, siendo uno de éstos el invernadero multitúnel tipo millenium, aunque se debe resaltar que mantiene la característica de estructuras de bajo costo y ventiladas pasivamente (Villagran et al., 2018).
El microclima generado en estos invernaderos es dependiente del clima exterior. Al no contar con sistema de control de clima activo, el control de las variables ambientales relevantes para la producción agrícola se limita al intercambio de aire entre el exterior y el interior del invernadero a través de las áreas de ventilación dispuestas en la estructura. Esto ha generado que en las últimas décadas existan numerosos trabajos dedicados a la optimización y evaluación del fenómeno de ventilación natural en invernaderos (Baeza et al., 2009; Villagran et al., 2012; He et al., 2015).
El estudio del comportamiento microclimático de un invernadero puede realizarse a través de modelos empíricos de balance de energía, experimentación en campo, termografía digital o mediante simulación numérica aplicando técnicas como la dinámica de fluidos computacional (CFD, por sus siglas en inglés) (Chen, 2009). Esta herramienta permite realizar el análisis cualitativo y cuantitativo del comportamiento de variables de interés, como temperatura, humedad y concentración de especies en estructuras agrícolas o pecuarias bajo diferentes escenarios de simulación (Bournet and Boulard, 2010). Un estudio mediante modelos CFD se compone de tres fases de desarrollo: (i) el pre-proceso, fase donde se genera la geometría y el mallado del invernadero; se establecen las condiciones de frontera y se seleccionan las principales propiedades físicas del fluido en estudio; (ii) la fase del proceso, en la cual se definen los modelos numéricos a utilizar y la forma de discretización de las ecuaciones parciales no lineales que se involucran en este tipo de estudios, y es la fase donde se busca la convergencia del problema alcanzando así una solución; y (iii) la fase de postproceso, donde los usuarios pueden observar los resultados numéricos y grá cos obtenidos en la geometría estudiada; esta fase suministra la información necesaria para realizar la calibración y validación del modelo CFD utilizado (Espinel-Montes et al., 2015).
El objetivo de este trabajo consistió en desarrollar y ajustar un modelo numérico CFD 3D, con el fin de estudiar el comportamiento de las variables temperatura y humedad y su relación con la ventilación natural en un invernadero multitúnel tipo millenium colombiano bajo condiciones climáticas de la sabana de Bogotá, usado para la producción de rosa.
MATERIALES Y MÉTODOS
Sitio y descripción del invernadero
El trabajo experimental se desarrolló en un invernadero comercial de plástico tipo multitúnel millenium de 5610 m2 de área cubierta, perteneciente a una finca productora de rosas (4e35'32,44" N; 74e12'35,14" W; 2510 msnm) ubicada en el municipio de El Rosal, Cundinamarca, Colombia. El invernadero evaluado estaba compuesto de 15 naves, cada una con una luz de 6,9 m (Fig.1A), las alturas mínimas y máximas bajo canal fueron de 3,0 y 8,2 m, respectivamente, la distancia longitudinal del invernadero era de 55 m y estaba orientado en sentido este-oeste (E-W). Cada nave disponía de un área de ventilación cenital fija de 0,53 m de apertura; el área total de ventilación natural se complementaba con aperturas de 2,0 m en los cuatro costados de la estructura, para un total de área de ventilación de 856 m2.
Se evaluaron las condiciones meteorológicas exteriores para las 24 horas del día, dentro del periodo comprendido entre las 00:00 horas del 05 de junio y las 23:00 horas de 06 de julio del año 2016. El registro de las variables se realizó por medio de una estación meteorológica (Vantage Pro2 Plus, Davis Instruments, Hayward, California, USA), que integraba sensores de radiación global, temperatura, humedad relativa, precipitación, velocidad y dirección del viento; la frecuencia de registro era cada 10 minutos.
Simulación numérica
Las ecuaciones gobernantes de un flujo pueden ser representadas como ecuaciones de convección-difusión de un fluido para tres leyes de conservación, que incluyen las ecuaciones de momento, energía y transporte de un fluido compresible y en un campo tridimensional (3D) en estado estacionario:
donde p es la densidad del fluido(kg nr3), V es
el operador nabla que es un operador vectorial que puede ser aplicado a un escalar o a un vector y representa su gradiente en el espacio, Ø representa la variable de concentración, ∆ es el vector de velocidad (m s-1), es el coeficiente de difusión (m2 s-1) y S representa el término fuente (Pisciaet al., 2012).
El comportamiento de la humedad y la conservación de masa fue modelado a partir de la siguiente oeprestón (Rojanoet al., 20d4):
donde Y es la fracción local de masa de cada especie a través de la solución de la ecuación convección-difusión, E es el termino fuente de humedad, J es el flujo de difusión de las especies, y Ri es la tasa de producción de especies por componente i. Debido a que las velocidades de viento serán bajas en algunas zonas en el interior del invernadero, los efectos de flotabilidad influenciados por el cambio de densidad del aire estarán presentes (Espinel et al., 2015) y pueden ser modelados a través del modelo de Boussinesq de acuerdo con el trabajo desarrollado por Kacira et al. (2004), que se describe mediante la siguiente expresión:
Así mismo todos los escenarios simulados consideraron la ecuación energía y el modelo de radiación seleccionado fue el de ordenadas discretas (DO) con discretización angular. Este modelo permite realizar el análisis de clima en condiciones de periodo nocturno, simulando y resolviendo el fenómeno de radiación desde el suelo del invernadero hacia el ambiente exterior. Para tal fin se consideró el cielo como un cuerpo negro con una temperatura equivalente (TC) de 10°C para un escenario predominante de noches húmedas y nubadas (Iglesias et al., 2009).
Estos conjuntos de ecuaciones fueron resueltos mediante la utilización del software comercial ANSYS Fluent (v. 17.0) (Ansys Corporation, 2016), usando una malla no estructurada y el método de volúmenes finitos. La naturaleza turbulenta del flujo de aire se simuló utilizando el modelo de turbulencia estándar k-ε, modelo ampliamente usado y validado en estudios enfocados a invernaderos demostrando un ajuste y precisión adecuada (Piscia et al., 2015).
Dominio computacional y generación de la malla
El software de pre-proceso ANSYS ICEM CFD (v. 17.0) se utilizó para generar un gran dominio computacional compuesto por el invernadero y sus alrededores, esto con el fin de garantizar la no afectación de la solución numérica del campo de flujo fuera del invernadero, y para permitir una definición apropiada de la capa límite atmosférica (Rico-Garcia, 2011). Las dimensiones del dominio computacional fueron de 333,5; 195 y 60 m para los ejes x, y, z, respectivamente (Fig. 1C). Este tamaño se determinó siguiendo las pautas para el cálculo de CFD del entorno eólico alrededor de los edificios (Tominaga et al., 2008).
El dominio computacional se dividió en 56,190,701 volúmenes discretizados en el espacio, número de elementos que se obtuvo luego de verificar la independencia de las soluciones numéricas del flujo de aire a tamaños de malla con un número superior e inferior de elementos de acuerdo al procedimiento reportado por He et al. (2017). Un criterio fundamental para establecer la precisión de las soluciones obtenidas mediante CFD consiste en evaluar la calidad de la malla, este procedimiento se realiza desde la interfaz del software de pre-proceso (Fig. 1B). Los parámetros de calidad evaluados fueron el tamaño de las celdas y la variación del tamaño de celda a celda encontrando que un 90,5% de las celdas de la malla estaban dentro del intervalo de calidad alta (0,95-1). También se evaluó el ortogonal quality, donde el valor mínimo obtenido fue de 0,88, resultados que se clasifican dentro del rango de alta calidad (Flores-Velázquez et al., 2015).
El software comercial de procesamiento CFD ANSYS FLUENT (v. 17.0) se utilizó para realizar las simulaciones. El método de solución semi implícito para la ecuación de presión (SIMPLE) y el algoritmo de resolución de partición de presión implícita con la división de operadores (PISO) fue aplicado para resolver el campo de flujo del fluido simulado. Los criterios de convergencia del modelo fueron establecidos en 10-8 para la ecuación de energía y en 10-6 para las ecuaciones de continuidad, momento y turbulencia (Baxevanou et al., 2017).
Condiciones de frontera
El límite superior del dominio y las superficies paralelas al flujo fueron fijados con condiciones de frontera de propiedades simétricas para no generar pérdidas de fricción del flujo de aire en contacto con estas superficies. Las simulaciones consideraron las características atmosféricas del sitio experimental, tales como presión atmosférica de 74993 Pa, viscosidad del aire igual a 1,7E-05 kg m-1 s-1 y gravedad de 9,81 m s-2. Al límite inferior y las paredes del invernadero se les fijó una condición de frontera de pared antideslizante. Otras propiedades del polietileno y del suelo agrícola como calor especifico (Cp), conductividad térmica (k) y densidad (p) fueron establecidas de acuerdo con Villagrán et al. (2012)
Los límites perimetrales del dominio computacional se establecieron como límites de entrada de aire o salida de presión según el caso evaluar. Se consideró un perfil uniforme de velocidad del viento evaluando velocidades con valores entre 0,32 y 1,22 m s-1 y valores de temperatura medios para cada hora del día. Dichos valores fueron establecidos a partir de la información climática recopilada en el periodo de medición experimental.
Otros parámetros de entrada que son necesarios para alimentar el modelo de radiación son las propiedades ópticas del material de cubierta establecidas para el fenómeno de radiación infrarroja de onda larga. Los valores para el polietileno transparente fueron coeficiente absorción (a) de 0,69, coeficiente de transmisión (7) de 0,19 y coeficiente de reflexión (p) de 0,11. El modelo no incluyó cultivo alguno y adicionalmente se asumió una hermeticidad máxima del invernadero. Estas simplificaciones son válidas puesto que son aplicadas a cada una de los casos simulados y los errores que puedan derivarse tendrán el mismo grado de magnitud para cada escenario.
Escenarios considerados
La evaluación del comportamiento micro climático se realizó estableciendo como condiciones de entrada al modelo CFD, el valor de los parámetros meteorológicos medios obtenidos para cada hora evaluada (Tabla 1) y la configuración de ventilación estándar del invernadero estudiado, para este caso, fue ventilación lateral combinada con ventilación de techo para horas diurnas y ventilación de techo para las horas nocturnas, para cada simulación se calculó el Índice de renovación (IR) de acuerdo con la siguiente expresión:
donde IR es el índice de renovación horario (Vol h-1), Q es el caudal de flujo de aire por unidad de tiempo (m3 s-1) en las áreas de ventilación, es el volumen del invernadero (m3).
Validación del modelo desarrollado
La temperatura de bulbo seco y de bulbo húmedo del aire del invernadero se registró por medio de treinta termopares tipo T (cobre-Constantán) conectados a un número igual de registradores de datos (Cox-Tracer Junior, Escort DLS, Edison, New Jersey, USA). Los conjuntos de termopares se desplegaron uniformemente a lo largo del eje transversal y longitudinal del invernadero a una altura de 1,5 m sobre el nivel del suelo y se programaron para registrar datos con una frecuencia de 10 minutos, datos que posteriormente fueron analizados y organizados en valores medios para cada hora, y mediante el uso de ecuaciones psicrométricas se determinó el valor de la humedad relativa. Los datos obtenidos en forma experimental serán comparados con los obtenidos mediante simulación CFD.
RESULTADOS Y DISCUSIÓN
Validación del modelo
La evolución temporal de la temperatura y la humedad relativa en el interior del invernadero se puede observar en la Fig. 2. Los órdenes en magnitud de los valores de temperatura simulados estuvieron de acuerdo con los medidos experimentalmente, se obtuvo un error cuadrático medio (MSE) de 0,17°C, un error medio absoluto (MAE) de 0,36°C y un error porcentual absoluto medio (MAPE) de 2,54%. Para el caso de la humedad relativa se encontró una tendencia similar entre los datos simulados y medidos, obteniendo un valor de MSE igual a 5,70%, el MAE fue de 2,16% y un MAPE de 2,67%. En general se puede concluir que el ajuste del modelo CFD es adecuado ya que los errores encontrados están dentro de los rangos aceptables para los modelos de clima bajo invernadero (Baptista, 2007). Por lo tanto, este modelo CFD permitirá analizar de forma fiable el comportamiento de la temperatura y la humedad relativa del invernadero estudiado.
Periodo diurno Campo de flujo
En la Fig. 3 se presentan los patrones de flujo de aire en el interior del invernadero para las 08:00, 12:00 y 16:00 horas. Estos flujos de aire son generados por el fenómeno de ventilación natural influenciados por la convección libre vía flotabilidad o efecto térmico y la convección forzada generada por el efecto eólico (Baeza et al., 2012). Para la hora 08:00 se observa un flujo que ingresa por la ventana frontal de barlovento con una velocidad media de 0,31 m s-1 (Fig. 3A); este flujo de aire fresco desplaza el aire cálido del interior del invernadero generando un movimiento convectivo que se dirige horizontalmente en el eje longitudinal del invernadero buscando la
salida por la ventana frontal de sotavento, y en el eje vertical en dirección de la aperturas fijas de techo; este efecto combinado genera un Índice de enovación (IR) de 7,3 vol h-1. Para las horas 12:00 y 16:00 se observa un comportamiento similar con elocidades de 1,09 y 0,38 m s-1, respectivamente, generando unos IR de 26,31 y 11,31 vol h-1 para cada una de estas horas (Fig. 3B, C).
En la Tabla 2 se presentan los valores de IR para el periodo comprendido entre las 06:00 y las 18:00 horas. Se debe mencionar que ara ninguna de estas horas evaluadas el IR cumple con el valor mínimo recomendado para invernaderos ventilados naturalmente, que debe ser igual a 40 vol h-1 (Montero et al., 2006). Los valores de IR moderadamente bajos pueden generar condiciones microclimáticas subóptimas y heterogéneas en el interior del invernadero. Adicionalmente se debe mencionar que los IR presentan una relación directa con los valores de velocidad de viento exterior y los diferenciales térmicos entre el interior y el exterior del invernadero (∆T= Tmedia interior-Tmedia exterior),, relación v media interior media exterior'' que fue previamente reportada por López et al. (2011).
Comportamiento térmico
La Fig. 4 muestra la distribución de la temperatura del aire en una vista isométrica a 1,5 m sobre el nivel del suelo para las 08:00, 12:00 y 16:00 horas. Las simulaciones permitieron observar que el valor de la temperatura para
la hora 08:00 presenta un comportamiento heterogéneo con dos zonas marcadas desde la zona central del invernadero y los costados de este. La temperatura alcanzó un valor medio de 16,2°C y la diferencia entre el valor de temperatura máximo y mínimo fue de aproximadamente 0,93°C (Fig. 4A). El comportamiento es similar para las 12:00 y 16:00 horas (Fig. 4B, C) donde la temperatura media alcanzó valores de 19,41 y 17,41°C, respectivamente, y las diferencias entre los valores máximos y mínimos de temperatura para cada caso son de 1,98 y 1,26°C.
En la Tabla 2 se presentan los valores de temperatura media interior y AT para las horas del periodo diurno; se debe destacar que estos valores son altamente dependientes de las condiciones presentadas en el ambiente exterior (Tabla 1) y principalmente del valor
de la radiación. En general se observa que el invernadero exhibe condiciones térmicas aceptables para la producción de rosas en el intervalo comprendido entre las 09:00 y 16:00 horas, donde la temperatura se mantiene dentro de los rangos óptimos recomendados para este cultivo, los cuales deben oscilar entre 17 y 25°C (Kinfe et al., 2014) y preferiblemente entre 21 y 24°C (Yong, 2004). Lo anterior se traduce que bajo las condiciones presentadas en el invernadero los ciclos vegetativos pueden alargarse debido a que los valores de temperatura no superan los 21°C. Esta situación puede ser compensada a través de un manejo adecuado de las áreas de ventilación, limitando el grado de apertura lo que permitiría aumentar la ganancia térmica del invernadero hasta generar las condiciones deseadas, aunque debe mencionarse que esta actividad debe ser realizada teniendo en cuenta los valores de humedad y CO2 presentes en el invernadero.
Comportamiento de la humedad relativa El comportamiento de la Humedad Relativa (HR) se puede observar en la Fig. 5. La distribución espacial de HR a 1,5 m sobre el nivel del suelo muestra un comportamiento heterogéneo para las 08:00, 12:00 y 16:00 horas. El valor medio de HR obtenido para la hora 8:00 fue 83,2%, bajo esta condición se encontraron zonas con valores máximos de HR de 87,3% y valores mínimos de 73,4% (Fig. 5A), al igual que en la variable temperatura este comportamiento no es recomendado para la producción agrícola, ya que se somete a las plantas a condiciones diferenciadas en un ambiente donde el manejo de las plantas se hace de forma homogénea sobre toda el área de cultivo establecida. Esto se traduce en una producción final heterogénea tanto en cantidad como en calidad. Por otro lado, el ataque de patógenos es más frecuente en las zonas donde las condiciones microclimáticas son más críticas. La HR muestra una relación inversa con la variable temperatura ya que las zonas de altas y bajas temperaturas concuerdan con las zonas de bajas y altas HR. Es importante mencionar que este comportamiento heterogéneo puede estar influenciado por los valores bajos de IR presentados para este invernadero. Para las horas 12:00 y 16:00 el comportamiento encontrado es similar al de la hora 08:00 y los valores medios de HR observados son de 70,5 y 69,1% respectivamente (Fig. 5B, C).
En la Tabla 2 se observan los valores de HR para las horas diurnas. Se puede mencionar que estos valores se encuentran en unos rangos óptimos para la producción de rosa en el periodo comprendido entre las 09:00 y 18:00 horas. Para las horas 06:00 y 07:00 la HR presenta valores superiores al 85%, estas mismas horas presentan temperaturas inferiores a 15°C; este comportamiento no es recomendado ya que puede inducir enfermedades en el follaje de la planta tales como mildiu velloso (Peronospora sparsa) y Botritys cinerea, enfermedades que en el contexto colombiano son frecuentes y generan una alta pérdida económica, debido al uso intensivo de fungicidas para el manejo y control de dichas enfermedades (Castillo et al., 2010; Restrepo, 2010).
Campo de flujo
La Fig. 6 muestra el comportamiento de los flujos de aire para las 20:00, 00:00 y 04:00 horas, bajo estas condiciones donde la velocidad exterior del viento es baja (Tabla 1) y las aperturas laterales de ventilación se encuentran cerradas. El movimiento del aire es inducido por el efecto térmico de la ventilación natural. Estos movimientos se caracterizan por ser de tipo convectivo rotatorio en dirección desde el suelo hacia la cubierta del invernadero. Este fenómeno ocurre debido a almacenamiento de calor por parte del suelo en las horas diurnas y su posterior liberación a la atmósfera en las horas nocturnas (Mesmoudi et al., 2017). Para la hora 20:00 pueden diferenciarse tres áreas de flujos convectivos donde la velocidad media del flujo de aire es de 0,28 m s-1 (Fig. 6A). Para las horas 00:00 y 04:00, las velocidades de los flujos de aire están por debajo de los 0,1 m s-1 y se logran identificar movimientos convectivos en cada una de las naves que componen el invernadero (Fig. 6B, C).
En la Tabla 3 se presentan los IR para las horas nocturnas evaluadas. En general estos IR son bajos, ya que el intercambio de aire se realiza a través de las aperturas de techo la cuales permanecen abiertas y permiten el intercambio de aire entre el interior y el exterior del invernadero. Este movimiento de aire está fuertemente influenciado por el fenómeno de flotabilidad.
Comportamiento térmico
La temperatura bajo los escenarios simulados de condición nocturna presentó un comportamiento espacial homogéneo (Fig. 7). El valor medio de la temperatura para las 20:00, 00:00 y 04:00 horas, fue de 13,15; 11,13; y 9,85°C, respectivamente. Esta tendencia decreciente del valor de la temperatura es la que ocurre típicamente en los invernaderos pasivos que carecen de sistemas de calefacción.
En la Tabla 3 se presentan los valores medios de temperatura en el interior del invernadero para el periodo comprendido entre las 18:00-23:00 y las 00:00-05:00 horas. En general se observa que estos valores están por debajo del valor mínimo recomendado para la producción de rosas, que es de 16°C; esta condición particular puede ocasionar desordenes fisiológicos en las plantas y flores de baja calidad comercial (Yong, 2004).
En la Tabla 3 se pueden encontrar los valores de AT. Los resultados reflejan que en el periodo comprendido entre las 23:00 y las 05:00 horas el valor de AT es negativo. Esto es un indicativo que el invernadero se encuentra bajo condiciones de
inversión térmica, fenómeno donde el ambiente en el interior del invernadero presenta un valor térmico inferior al del ambiente exterior. Esto se genera principalmente por una pérdida acelerada de calor a través de la cubierta y del flujo de calor extraído del interior del invernadero a través de las aperturas cenitales fijas (Montero et al., 2013). La presencia del fenómeno de inversión térmica no es recomendable en invernaderos usados para la producción agrícola, ya que las plantas estarán sometidas a condiciones térmicas subóptimas, que pueden limitar de forma parcial o total la producción de rosa. Por lo tanto, se sugiere la búsqueda de estrategias que permitan mejorar la hermeticidad del invernadero evaluado.
Comportamiento de la humedad relativa
La evolución espacial de la HR en el interior del invernadero para las horas 20:00, 00:00 y 04:00 se muestra en la Fig. 8. Se puede observar un comportamiento heterogéneo de esta variable climática; para la hora 08:00 el valor medio de la HR fue de 85,9%, se observan zonas con valores mínimos y máximos de HR de 81 y 89%,
respectivamente (Fig. 8A). El valor medio de la HR para la hora 00:00 es igual a 94,1%, con algunas regiones con valores mínimos y máximos de HR de 91 y 96% (Fig. 8B). Para la hora 04:00 se observa un comportamiento similar con un valor medio de HR de 98,2%.
Los valores para el comportamiento temporal de la HR se presentan en la Tabla 3. En general se puede observar que para el periodo evaluado el valor de HR es superior al 80% con una tendencia creciente a medida que trascurre la noche. Para el periodo comprendido entre las 1:00 y 5:00 horas, el valor de la HR presentó valores superiores al 94%, esto es indicativo de un aire altamente saturado lo cual es una característica común del ambiente exterior en la sabana de Bogotá. Esta condición no es adecuada para la producción de rosas, ya que los valores altos de HR en combinación con el fenómeno de inversión térmica y déficit de presión de vapor bajos, suelen generar condensación del vapor de agua en la cubierta del invernadero y en la superficie foliar de las plantas, creando las condiciones favorables para la aparición de enfermedades fúngicas (Baptista et al., 2012). Adicionalmente en los invernaderos pasivos cuando aparece el fenómeno de condensación, éste se extiende hasta las primeras horas del periodo diurno, lo cual tiende a reducir la tasa de transpiración de las plantas, generando desórdenes fisiológicos debido a la translocación de algunos nutrientes como el calcio (Wang et al., 2012).
CONCLUSIONES
El modelo numérico CFD-3D fue validado mediante la toma de datos experimentales, que permitieron realizar una comparación entre los datos simulados y medidos. Esta comparación mostró que el modelo utilizado tiene una alta capacidad para representar el comportamiento real de las variables estudiadas. El grado de ajuste es adecuado y los errores de simulación están por debajo del 10%.
El análisis desarrollado permitió establecer que el comportamiento de las variables temperatura y HR para el periodo diurno, presentó valores adecuados para la producción de rosas, caso contrario sucedió para el periodo nocturno donde el comportamiento de estas variables mostró valores subóptimos para su producción. Adicionalmente, se observó la presencia frecuente de inversión térmica, fenómeno que, complementado con la presencia de altas humedades relativas, favorece la condensación del vapor de agua, lo cual genera ambientes propicios para el ataque de enfermedades fúngicas y el desarrollo de desórdenes fisiológicos en las plantas. El enfoque de este trabajo puede ser utilizado en investigaciones futuras, aunque se recomienda incluir el cultivo y observar las variaciones que puedan ocasionar las plantas al microclima generado.