El Azúcar de Caña no centrifugada, (non-centrifugal cane sugar – NCS en inglés), también llamada “panela” (1), es obtenida de la extracción, purificación y concentración de los jugos de caña de azúcar (2). Colombia con el 14.9% de producción mundial, es el segundo productor, después de la India (3, 4), y presenta el mayor consumo per cápita en el mundo (24 kg/persona/año) (5).

En Colombia la elaboración de la panela se lleva a cabo en edificaciones denominadas “trapiches” (6), las cuales son plantas de procesamiento agroindustrial naturalmente ventiladas. Se estima que en Colombia hay cerca de 18.473 trapiches, su tamaño está sujeto a la capacidad de producción que incluyen unidades con rangos de capacidad pequeña (menor a 100 kg h^-1) , mediana (100 a 250 kg h^-1) y grande (mayor a 250 kg h^-1) (7-9).

En estas edificaciones se ubican los equipos necesarios para la transformación, especialmente, el horno, donde se lleva a cabo la evaporación de más del 90% del agua contenida en el jugo de caña (10). Lo anterior, hace que estas edificaciones emitan calor y vapor que pueden comprometer el confort y desempeño de los ocupantes en su interior.

Factores ambientales como la temperatura y humedad relativa contribuyen al confort, salud, rendimiento y productividad de los trabajadores (11). Si el proceso de diseño del edificio ignora los factores climáticos, se pueden generar afectaciones en las actividades de los ocupantes (12).

Una alternativa de diseño de edificaciones es la arquitectura bioclimática, la cual es una forma alternativa de diseño y construcción de edificaciones que involucra condiciones climáticas del entorno, el uso de recursos naturales y técnicas tradicionales y modernas de construcción (13), especialmente en el diseño de la envolvente del edificio (14), con el objetivo de proveer confort térmico a sus ocupantes (15), bajo parámetros de sostenibilidad y eficiencia energética.

La envolvente del edificio divide el ambiente interior de las condiciones climáticas externas (16), y está compuesta por cubiertas, paredes, ventanas, puertas y pisos (17). La envolvente interviene para regular el balance de masa (de aire) y energía entre la edificación y el exterior (18).

Entre las estrategias del diseño bioclimático pasivo se encuentra la ventilación natural, la cual tiene como objetivo, la sustitución del aire interior nocivo, por aire fresco, sin utilizar sistemas o equipos mecánicos (19, 20), esta puede ser térmica o mínima (higiénica), y con ella se busca reducir el consumo energético, incrementar el confort térmico y mejorar de la calidad del aire (21-23).

Colombia, es un país ubicado en la zona ecuatorial, donde predominan los climas cálidos y templados. Al interior de los trapiches se genera calor y vapor del proceso agroindustrial, por lo que estas edificaciones son naturalmente ventiladas. La ventilación natural en estas edificaciones es dinámica y térmica, producto la acción del viento, y el efecto chimenea que generan los gradientes de temperatura en su interior. Es de aclarar que la ventilación térmica en situaciones de calor, suplen las necesidades de la ventilación higiénica al tener tasas de renovación mayores a las requeridas para la calidad del aire por emisiones de gases.

Cuando la temperatura del ambiente interno es mayor que la temperatura exterior, el aire externo más frio y denso ingresa a la edificación a través de las aberturas bajas, y el aire interno más caliente y menos denso asciende a la parte superior de las cubiertas o sale por las aberturas elevadas; este fenómeno se conoce como efecto chimenea (24), el cual es fundamental en las edificaciones tipo trapiche, por la emisión de calor y vapor en su interior.

En la industria es necesario mantener confort térmico, ya que permite una mayor comodidad y productividad, al tiempo que reduce el riesgo de accidentes en el trabajo y fatiga laboral (25, 26). En ambientes cálidos la incomodidad puede ser causada por alta humedad en la piel generada por niveles altos de humedad del aire (27).

La humedad relativa elevada reduce la capacidad del cuerpo para perder calor a través de la transpiración y la evaporación (28). Según (29), humedades relativas del 80 % al 95 %, pueden causar graves riesgos para la salud, además, puede causar daños en las superficies del edificio (30), por la condensación (31), situación que se presenta en las cubiertas de las instalaciones agroindustriales de fabricación de panela.

Para la clasificación de confort térmico, se han desarrollado diferentes índices, que agrupan en un parámetro el efecto conjunto de variables como: la temperatura de bulbo seco, humedad relativa, velocidad del aire y radiación térmica, y factores humanos como: el aislamiento de la ropa y el nivel de actividad metabólica (32).

Un índice utilizado para determinar el confort térmico en entornos de trabajo, es la temperatura efectiva ( T e ) (33-35), La cual está en función de la temperatura de bulbo seco, humedad relativa y velocidad del aire.

En este tipo de edificaciones, el proceso se lleva cabo en condiciones de alta humedad relativa que puede oscilar entre 80% y cerca de la saturación, y alta temperatura, que puede alcanzar entre 60 y 150°C (36) . Estas condiciones ambientales adversas, combinadas con la actividades físicas pesadas propias del proceso, pueden causar estrés térmico y problemas de salud, cambios en el comportamiento mental, disminución del rendimiento, aumento del riesgo de accidentes, y conducir a un deterioro progresivo del trabajador (37).

Varias investigaciones utilizando Dinámica de Fluidos Computacional (CFD), han estudiado el ambiente interno de edificaciones, incluyendo construcciones rurales, y han demostrado ser una herramienta valiosa y eficaz en este tipo de análisis (38), ayudando a reducir costos y tiempo en el diseño y evaluación bioclimática (39).

Relacionados al sector panelero se encuentran estudios utilizando CFD orientados al diseño del horno y las pailas o fondos que la integran (40-42), pero no se registran estudios en instalaciones de producción tradicionales tipo trapiche.

Este estudio tuvo como objetivo analizar el efecto de tratamientos con diferentes configuraciones de ventilación natural, combinadas con tres tipos de materiales de cubierta, sobre las condiciones de confort térmico de los ocupantes, dentro de una instalación para la producción de panela, utilizando CFD, con el fin de preservarlas y mejorarlas.

Este estudio se realizó en una edificación para la fabricación de panela, en el municipio de Santo Domingo - Antioquia - Colombia (N06°28′15″;W075°09′57″, 1440 msnm), el cual posee una temperatura media de 20°C, y humedad relativa de 80% (43). Su capacidad de producción es de 150 kg h^-1 (mediano), provista de hornos con cámara de combustión doble tipo Plana-CIMPA, que es un diseño mejorado de la cámara tradicional desarrollada por el Centro de Investigación para el Mejoramiento de la Agroindustria Panelera (CIMPA) (10). Las dimensiones de este trapiche son las aproximadas para una unidad con capacidad de producción de tamaño mediano y su distribución arquitectónica permite un flujo secuencial del proceso de producción.

En cada una de las tres áreas críticas de trabajo de la instalación (Hornilla, Apronte y Moldeo) se instaló un Datalogger de humedad y de temperatura marca EXTECH RHT10, rango de medida de -40ºC a 70ºC y de 0% a 100% de humedad relativa, precisión de 0.1ºC y 0.1% respectivamente. Se tomaron datos de 3 moliendas y 2 días de referencia para las 214 horas. En el área de bodega no se instaló dispositivo porque es utilizada como área de tránsito y almacenamiento.

Se simuló el ambiente interno del edificio para las operaciones que se llevan a cabo en el área del horno u hornilla, ya que, cuando esta está en funcionamiento se presenta la situación más crítica para las condiciones de operatividad.

La Figura 1 muestra la distribución arquitectónica de la instalación, la cual tiene un piso en forma escalonada con tres niveles, lo que facilita el transporte del jugo de caña por gravedad en la hornilla. En el primer nivel está el área de batido y moldeo, el área de apronte de caña, el área baja de la hornilla, en el segundo nivel está el área de bodega, el molino de caña, el prelimpiador, el área de apronte de caña, el área media de la hornilla, y en el nivel tres está la chimenea y el área alta de la hornilla.

Las dimensiones de este edificio son, 16.00 m de ancho x 30.00 m de largo x 3.40 m de alto en promedio (Figura 2).

La edificación cuenta con una cubierta a dos aguas discontinua en dos niveles, y pendiente de 20%. Cada discontinuidad configura una ventana cenital, específicamente tiene una ventana cenital en la cumbrera de 27.70 m de largo × 0.65 m de ancho sentido este - oeste, una segunda y tercera ventana cenital de techo discontinuas por cada una de las dos aguas de la cubierta de 5.15 m de largo x 0.65 m de ancho sentido norte - sur costado norte, y una de 6.15 m de largo x 0.65 m de ancho sentido norte - sur y costado sur respectivamente, y una cuarta ventana cenital de 18.00 m de largo x 0.65 m de ancho en el costado sur, sentido este - oeste, (Figura 2).

El edificio tiene cerramiento en ladrillo perforado en sus caras laterales. La cubierta posee, tejas de fibrocemento, cerchas y perfiles metálicos de 0.08 x 0.04 m. El piso está compuesto por una losa de contrapiso en concreto reforzado.

En la Figura 3 se muestran las geometrías 3D de la instalación en formato CAD, la primera es una condición con cerramiento en ladrillo tipo perforado para evaluar tres distribuciones constructivas, y la segunda es una condición base sin cerramiento perimetral, utilizada para evaluar nueve distribuciones constructivas.

Se simularon diferentes opciones de aberturas de ventilación natural con el fin de determinar la influencia del área de ventilación, y la ubicación de las aberturas en la temperatura y humedad relativa dentro de la instalación. Los diferentes tratamientos se pueden ver en la Tabla 1 y Figura 4 , donde el tratamiento T1 es la geometría original (tratamiento testigo) que fue validada y simulada.

Posteriormente, se realizó la malla computacional tetraédrica con ayuda del software ANSYS-ICEM-CFD con 128.197 nodos, y 629.596 elementos para el modelo con cerramiento en ladrillo perforado (Figura 5a), y 154.197 nodos y 776.124 elementos para los modelos sin cerramiento perimetral y con cerramiento perimetral con ladrillo macizo (Figura 5b). Cuando las geometrías son complejas, o el rango de longitudes de escala es grande, como se presentó para los dos modelos propuestos, una malla tetraédrica es apropiada, porque requiere un número menor de elementos (44). Se obtuvo diferentes tamaños de elementos después de varias pruebas de malla donde el refinamiento mayor se encuentra en las aberturas, debido a la necesidad de mayor detalle de análisis fluidodinámico por ventilación natural.

Una vez se ha generado la malla computacional, se importó a la plataforma ANSYS-CFX-Pre, en este punto se insertaron las condiciones de frontera del modelo: propiedades de los materiales, condiciones del ambiente externo, componente energético de la hornilla, entre otros. Las propiedades térmicas de los materiales de construcción utilizados como condiciones de frontera se muestran en la Tabla 2.

Se realizó un modelo transitorio a partir de las condiciones de la actividad crítica durante el proceso de fabricación de panela: flujo másico de vapor producido en promedio por la hornilla es de 0.15 kg de agua s^-1, temperatura promedio de las pailas o fondos y jugos en la hornilla alcanzan un valor de 132°C, temperatura del suelo de 20°C, temperatura media del exterior de 20°C, humedad relativa media externa de 80%, conductividad térmica del aire es 0.02514 W m ^-1 K ^-1 y la difusividad térmica del aire es 0.02514 W m ^-1 K ^-1 y la difusividad térmica del aire de de 2.074 × 10^-5 m² s^-1 (31).

La humedad relativa se insertó como una variable adicional a partir de la ecuación [1] (46), en función de la presión de vapor del aire y la presión de vapor de saturación del mismo representada en la ecuación [2] de Magnus-Tetens (47).

Donde, HR es la humedad relativa (%), es la presión parcial de vapor (hPa), es la presión de vapor de saturación (hPa), es la temperatura de bulbo seco (°C).

Las ecuaciones de transporte son ecuaciones de derivadas parciales utilizadas para describir fenómenos de transporte de energía, masa, y cantidad de movimiento (48), y se expresan en términos generales en las ecuaciones [3], [4] y [5] (41), donde, ρ ¯ es la densidad media de Reynolds en kg m^-3; ũ_j, ũ_i, ũ_k son los vectores de velocidad promedio de Favre en m s^-1; xi, xi, xk son las posiciones del vector; u i " , u j " son los vectores de velocidad fluctuante de Favre en m s^-1; gi es el vector aceleración de la gravedad en m s^-2; μ es la viscosidad dinámica en Pa s; h ~ t Favre promediado de la entalpía total en kJ kg^-1; Pr Número de Prandt; S_r término promediado debido a la transferencia de calor por radiación en W m^-3; p es la presión estática media en Pa.

Posteriormente, en la plataforma SOLVER se simuló el modelo previamente realizado en ANSYS-CFX-Pre, para ello se utilizó como criterio de convergencia un error máximo de 10^-4 para todos los tratamientos.

Para el modelo de turbulencia se utilizó el método de turbulencia k-ε estándar, que evalúa la viscosidad (μτ) a partir de una relación entre la energía cinética turbulenta (k) y la disipación de energía cinética turbulenta (ε) (44). El modelo k-ε estándar se ha utilizado en estudios anteriores en aplicaciones industriales e instalaciones agrícolas (49, 50).

A partir de los resultados obtenidos en la modelación computacional, se utilizaron los valores obtenidos de temperatura de bulbo seco y humedad relativa dentro de la instalación, ubicando puntos en la edificación de las áreas de las actividades de molienda y procesamiento de jugos (clarificación, evaporación y concentración de jugos), distribuidos en cuatro áreas: hornilla, cuarto de batido moldeo, apronte de caña y bodega (Figura 2).

En la Figura 6, se indican las posiciones de los puntos ubicados en él edificio, en el plano horizontal los puntos están equidistantes a 3.00 m y en el plano vertical a tres alturas 0.80m, 1.50m y 2.20m sobre el nivel del suelo, para formar una malla de 180 puntos, identificando cada área con puntos de color. El objetivo de estas tres alturas fue determinar medidas de temperatura y humedad relativa a la altura de las mesas 0.80m, a la altura promedio de la cabeza del trabajador 1.50m y por encima de la cabeza del trabajador 2.20m.

En cada uno de los 180 puntos ubicados en el edificio se determinó la temperatura efectiva (T_e) (34), el cual está en función de la temperatura, la humedad relativa y la velocidad del aire, con el fin de evaluar las condiciones bioclimáticas de los trabajadores.

Por medio psicrométria y los datos de temperatura y humedad relativa, se determinaron los valores de la temperatura de bulbo húmedo. La ecuación [6], define la temperatura efectiva en un frente de trabajo (34):

Donde, T_e es la temperatura efectiva en °C, T_bs es la temperatura bulbo seco en °C, T_bh es la temperatura de bulbo húmedo en °C, V_a es la velocidad de la corriente del aire en m s^-1.

La temperatura efectiva considera las variables ambientales de temperatura y humedad relativa, con el fin de evaluar los resultados obtenidos en la simulación, se compararon los tiempos de permanecía en los frentes de trabajo establecidos en el Reglamento (34).

Los resultados del modelo CFD para el T1 testigo (geometría original), fueron validados por medio del error cuadrático medio normalizado (NMSE) recomendado por la ASTM (51) por medio de la ecuación [7].

Donde:

Ypi: Valor predicho

Ymi: Valor medido

n: número de datos

Para esto, se realizó el promedio de los datos de temperatura y humedad relativa cada hora para un periodo de 24 horas simulados para las áreas de Hornilla, Apronte y Moldeo. Los valores con un NMSE inferior a 0.25 se consideran como buenos indicadores de concordancia; a medida que este valor se acerca a cero, la concordancia entre los valores medidos y predichos es mayor.

Posteriormente se realizó un análisis de varianza, acompañado de un análisis de medias mediante la prueba de Tukey, con el fin de observar si existe diferencia estadística entre los diferentes tratamientos o configuraciones de la envolvente, además, se realizaron diagramas de caja para las variables temperatura, humedad relativa y temperatura efectiva.

Con fines de validación del modelo, se realizó el cálculo del NMSE para el modelo CFD y los datos experimentales del tratamiento T1 o testigo (Tabla 3), para temperatura y humedad relativa, indicando buena concordancia debido a que todos los valores del NMSE son menores a 0.25, por lo que se concluye que puede ser utilizado para predecir el comportamiento de la temperatura y la humedad relativa al interior de la instalación.

La Tabla 4, muestra el test de Tukey, para análisis estadístico medias de temperatura, humedad relativa y temperatura efectiva para los doce tratamientos, en el cual para los grupos unidos con la misma letra no existe diferencia estadística significativa.

Para los grupos unidos con la misma letra no hay una diferencia estadística significativa por la prueba de Tukey (Temperatura P = <0,001 y F= 227.803; Humedad relativa P = <0.001 y F= 766.009; Temperatura efectiva P = <0,001 y F= 227.803).

Es de destacar que no hubo diferencia estadística de temperatura, humedad relativa y temperatura efectiva en los tres tipos de cubierta para cada tipo de cerramiento, lo que indica que en este tipo de edificaciones que emiten gran cantidad de calor y vapor, es más importante el cerramiento y la ventilación natural que el tipo de cubierta, para el intercambio de masa de aire caliente y húmedo con el exterior (52).

Los tratamientos T4, T5 y T6 (perímetro abierto y ventana cenital abierta, con cubierta de fibrocemento, zinc y barro respectivamente), con un valor alrededor de los 23°C, humedad relativa de 61.6%, y temperatura efectiva de 19.5°C, presentaron las mejores y únicas condiciones ambientales saludables para el trabajo, de acuerdo al (33), el cual sugiere rangos de temperatura en el ambiente laboral entre 14°C y 25°C y rangos de humedad relativa entre 30% y 70%.

Por otra parte, la situación más crítica se presentó en los tratamientos T10, T11 y T12 (cerramiento en muros de ladrillo macizo y ventana cenital abierta, con cubierta de fibrocemento, zinc y barro respectivamente), con valores de temperatura entorno de 86°C, humedad relativa cerca de la saturación, y valores de temperatura efectiva alrededor de 85°C.

Para los tratamientos T1, T2 y T3 (cerramiento de muros de ladrillo perforado y ventana cenital abierta con cubiertas de fibrocemento, zinc y barro respectivamente), no hubo diferencia estadística en las tres variables. La temperatura media se encuentra alrededor de 44°C, humedad relativa cercana a 73%, y temperatura efectiva alrededor de 41°C.

Aunque los tratamientos T1, T2 y T3 registran valores de temperatura efectiva aproximadamente de 4°C por encima de los tratamientos T7, T8 y T9, no se puede asegurar que exista diferencia estadística significativa entre estos seis tratamientos.

En el área circundante a la hornilla se presentan los mayores valores de temperatura y humedad relativa, ya que esta zona es la mayor fuente de calor y vapor,

Es importante considerar alternativas constructivas que faciliten la evacuación al ambiente externo del vapor generado por la hornilla, especialmente en países tropicales donde los niveles óptimos de humedad relativa para el confort se sitúan entre el 30% y el 70% (53).

La Tabla 5 muestra los tiempos de permanencia en frentes de trabajo. Para desarrollar labores sin límite de permanencia, los valores de una temperatura efectiva deben estar por debajo de 28°C, los tratamientos T4, T5 y T6, cumplen con esta condición y los requerimientos mínimos de confort térmico, los demás tratamientos no se recomiendan como ambientes saludables para los trabajadores.

Como ya se ha mencionado, el área más crítica de la edificación es la circundante a la hornilla, los puntos tomados alrededor de esta se muestran en la Figura 6, identificados con puntos de color amarillo. Por la relevancia de esta área, se realizó un análisis de medias de la temperatura efectiva para ella, la cual se muestra en la Tabla 6.

Se puede observar que al igual que lo mencionado anteriormente, en los tres tipos de cubierta para cada tipo de cerramiento, pero si en el tipo de envolvente. También se puede observar que los tratamientos T4, T5 y T6, de acuerdo a los rangos presentados en la Tabla 6, para estos tratamientos los valores de la media de la temperatura efectiva para el área hornilla, se encuentran por debajo de 28°C, que corresponde a un tiempo de permanencia sin limitaciones, de lo que se concluye, las aberturas en los laterales y la abertura en la ventana cenital ofrece las condiciones higrotérmicas adecuadas al interior de la instalación para los trabajadores, independiente del material de la cubierta.

Los demás tratamientos superan el límite máximo de 32°C para permanecía según el índice de temperatura efectiva, lo que sugiere, que no es recomendable que los trabajadores laboren con las condiciones, y que se deben adoptar medidas para mejorar la extracción de calor y vapor del interior, sea por medios pasivos (modificación de la envolvente), o por medios activos como el uso de extractores.

Resultado de la modelación CFD, se obtuvieron los perfiles de temperatura y humedad relativa, para los tratamientos T1, T4, T7 y T10, siendo los tratamientos representativos de la modelación computacional, debido a que no hay diferencia estadística según los materiales de la cubierta.

En el tratamiento T1 (Figura 7) las temperaturas oscilan entre los valores de 60°C y 80°C (333.15K y 353.15K) principalmente en el área que colinda con la emisión de vapor generada en los fondos o pailas de la hornilla. En las áreas que no hay contacto directo con la hornilla se observan temperaturas entre 20°C y 40°C (293.15K y 313.15K) presentándose las temperaturas más altas en el área que corresponde al costado derecho de la hornilla, y las temperaturas más bajas al costado izquierdo. El ladrillo perforado disminuye el área de ventilación, pero permite el ingreso del aire lo que ayuda a evacuar el aire caliente.

Con respecto a la humedad relativa en el tratamiento T1, los valores oscilan entre el 80% y 90%, al reducirse el área de ventilación en las laterales, lo que en este caso no es suficiente para ofrecer valores de humedad relativa por debajo del 70% y ofrecer un ambiente no sofocante.

Se encontró que el tratamiento T4 (Figura 8), exceptuando el área que corresponde a los límites con la hornilla, la temperatura osciló entre 20°C y 25°C (293.15K y 298.15K) en casi toda la instalación. Se encontraron temperaturas de 60°C y 80°C (333.15K y 353.8K) en la superficie interna de la cubierta, esto por efecto del ascenso de vapor caliente que genera el proceso, y por la acumulación de calor en esta región por el efecto chimenea.

Por otro lado, se observa que la humedad relativa oscila entre el 60% y 65% en la mayoría del espacio de la instalación, valores aceptables para ofrecer un confort térmico a los trabajadores.

En el tratamiento T7 (Figura 9) al igual que en el tratamiento T1, la temperatura oscilan entre 60°C y 80°C(333.15K y 353.15K) en casi toda la instalación, esto se debe que al cerrarse o no poseer la ventana cenital, no hay salida del vapor por efecto chimenea. Se genera un volumen de aire caliente y húmedo en la parte superior de la edificación, con humedad relativa entre 70% y el 85%, y cerca de la ventana cenital de 90% a 100%, con peligros de condensación en la cubierta.

El tratamiento T10 (Figura 10), presenta la condición más crítica, ya que al cerrarse las paredes laterales se limita la entrada de aire fresco (más frío y denso), por lo que se incrementa la temperatura y la humedad relativa al interior a valores que oscilan entre 80°C y 97°C (353.15K y 370.15K) y humedad de 90% a 100%, por la limitación de expulsión de los vapores calientes generados en el proceso. Esta condición de ambiente sofocante, puede afectar considerablemente el confort y la salud de los trabajadores (33), y generara condensaciones de agua en las superficies internas del edificio, con implicaciones negativas para los ocupantes, el edificio y el producto.

Por lo tanto, el tratamiento T4, presenta las condiciones más adecuadas de temperatura y humedad relativa para el confort higrotérmico de los trabajadores, cumpliendo con los rangos de temperatura y humedad relativa establecidos por (33).

La humedad del aire puede llegar a la saturación cuando se cierran las paredes del edificio, o se prescinde de la ventilación cenital, lo que también puede afectar el confort térmico humano. Con el fin de evitar el uso de equipos de climatización y mejorar las condiciones ambientales internas, se sugiere potenciar el efecto chimenea en este tipo de edificaciones.

Al favorecer el movimiento del aire dentro de la instalación mediante el efecto chimenea, se logra reducir la temperatura efectiva (T_e) mejorando la sensación térmica del trabajador, esto se debe a que la temperatura efectiva está directamente relacionada con la velocidad del aire.

La Figura 11 presenta un diagrama de cajas de la temperatura efectiva para los 12 tratamientos, donde se puede observar la mediana y la dispersión de datos de cada uno de los tratamientos estudiados.

Los tratamientos, T1, T2 y T3 con valores de mediana de 37.08°C, 37.03°C y 37.62°C, Mientras que T4, T5 y T6 presenta valores de mediana de 16.38°C, siendo los valores más bajos, cumpliendo con los requerimientos de confort; los datos atípicos (de mayor valor) se presentan justo encima de la hornilla donde se produce el vapor caliente.

Los tratamientos, T7, T8 y T9 con valores de mediana de 35°C, 35.47°C y 35.82°C respectivamente, presentan un comportamiento similar a T1, T2 y T3, pero con mayor dispersión de datos.

Los tratamientos T10, T11 y T12 con valores de mediana de 85.28°C, 85.62°C y 86.73°C respectivamente, representan la condición más crítica en la instalación para los trabajadores.

En Colombia es común que este tipo de instalaciones agroindustriales presentan envolventes con perímetro abierto, pero no con ventana cenital o chimenea en la cubierta, por lo anterior, se recomienda con fines de eficiencia energética y confort, adaptar estrategias pasivas de climatización como adoptar diseños con perímetro abierto, y ventanas cenitales o chimenea para favorecer el efecto chimenea, y con ello el intercambio de calor y humedad al exterior, favoreciendo las condiciones adecuadas para el ocupante, el edificio y el producto.

Con base en los resultados obtenidos, se concluyó que la simulación en CFD se puede utilizar en el diseño y evaluación para mejorar las condiciones de confort de edificaciones para la producción de panela tipo trapiche, teniendo en cuenta el diseño bioclimático.

De los doce tratamientos que se propusieron y se evaluaron, se observó que los mejores comportamientos bioclimáticos se presentan en los tratamientos con perímetro abierto y ventana cenital abierta (T4, T5 y T6), envolvente que favorece la ventilación natural, especialmente el efecto chimenea, lo que permite ofrecer condiciones higrotérmicas adecuadas para los trabajadores de esta agroindustria.

Se evidenció la importancia de promover y potenciar el efecto chimenea en la ventilación natural como estrategia pasiva de climatización, con el fin de evacuar vapores calientes de este proceso, y evitar la utilización de equipos de climatización, mayores gastos y reducir huella de carbono en la utilización de los mismos.

Los resultados del índice de temperatura efectiva para los tratamientos de perímetro abierto y ventana cenital abierta (T4, T5 y T6) en el área total de la edificación se encuentran dentro del rango adecuado con valores de 22.58, 22.51 y 22.56 para T4, T5 y T6 respectivamente, y el tiempo de permanencia no presenta límites.