1. INTRODUCCIÓN
El sector textil, en conjunto con el sector de la confección, son responsables del 10% de los impactos ambientales del ciclo de vida europeo, alcanzando la cuarta posición en el ranking de categorías de productos que causan los mayores impactos ambientales [1]. Los impactos negativos asociados con este sector de forma general se pueden resumir en cinco principales categorías: (1) consumo de energía en la producción de materias primas (fibras sintéticas), fabricación de hilados, acabado de hilados y tejidos, lavado y secado de ropa, así como procesos de comercialización y venta; (2) consumo de agua asociado con el cultivo de materias primas, productos químicos de pretratamiento de la etapa de producción, tintes y acabados, así como lavandería; (3) uso de productos químicos (pretratamiento húmedo, teñido, acabado y lavado) y liberación en agua; (4) residuos sólidos que surgen de la fabricación, envasado y disposición final de hilados, tejidos y productos finales de productos al final de su vida útil; y (5) emisiones directas de CO2, particularmente las relacionadas con los procesos de transporte [2]. Debido a estos impactos negativos es importante destacar el control del consumo de energía asociado al funcionamiento de la planta valorando la eficiencia tanto de equipos principales, auxiliares, iluminación y climatización [3]. Para la climatización se debe estudiar el efecto de las características termo conductivas de los materiales de cerramiento en el consumo de energía siendo que esto aporta a los costos de operación de las áreas productivas del proyecto [4]. También, identificar costos asociados al consumo energético de los equipos seleccionaos, el costo asociado al mantenimiento, la demanda esperada y la vida útil de los componentes [5]. Además, de los cuidados de gestión energética existen estudios enfocados a la gestión de los residuos textiles postindustriales posteriores, los cuales ayudan a desarrollar materiales de aislamiento acústico utilizando desechos mixtos de algodón / poliéster y caucho natural [6].
Desde el punto de vista económico, algunos trabajos muestran cómo lograr una comprensión profunda de las demandas ecológicas de los clientes objetivos y una referencia para el diseño y producción de productos de moda sostenibles que ahorren energía. Esto se logra mediante un procedimiento sistemático que identifique los factores energéticos para proyectos sostenibles y el desarrollo de escenarios factibles y prácticos [7].
Por lo anterior expuesto reducir la huella medioambiental para la industria es necesaria, para ello la sostenibilidad y las plantas textiles deben formar parte de una economía circular, alargando la vida de los productos mediante la conservación de las fibras textiles dentro de un circuito cerrado. Esto se logra clasificando la ropa para reciclaje, donde usualmente se hace de manera manual, una propuesta de una máquina que analice los espectros infrarrojos de las muestras textiles y aplicando métodos estadísticos multivariados adecuados puede lograr la clasificación precisa permitiendo el reciclaje de los tejidos [8].
En industrias como la agrícola existen técnicas y sistemas concernientes a los sistemas de enfriamiento pasivo aplicados en la construcción ganadera, estos principios son de la misma forma aplicables a la industria textil [9]. Propuestas como la de un sistema de enfriamiento de tratamiento de aire híbrido, que incluye un sistema de tratamiento de aire “plug-and-play”, promueve la eficiencia energética mejorando la calidad del aire interior, mediante una unidad de tratamiento de oxidación basada en ozono y un dispositivo de depuración de aire. Los resultados de esta propuesta han indicado que un menor caudal de aire exterior conduce a una reducción significativa del consumo total de electricidad para refrigeración [10].
Trabajos en los que se estudia el consumo energético medio, específico y por equipos de la industria textil, revelan un potencial ahorro de energía calculando las emisiones de CO2 para la producción de una pieza de prenda, lo que demostró que la participación considerable en el consumo de energía la producción de vapor, los compresores y los equipos de iluminación es relevante [11]. No obstante, las aplicaciones de la industria textil no se limitan solo al confeccionadito de ropa, sino que también se usan para recubrir paredes en el campo de la ingeniería civil, por ejemplo, se han utilizado TRC (textiles de carbono, matriz cementosa) para reforzar o reforzar los elementos estructurales de antiguas obras de construcción donde se analizan el efecto de la temperatura elevada sobre el comportamiento del TRC y el rendimiento del TRC [12].
Existen zonas con condiciones climatológicas que causan problemas a los sistemas de refrigeración como puede ser la India en una época del año, esto causa una sobrecarga el serpentín de enfriamiento y una fuerte disminución en el rendimiento del condensador enfriado por aire en un clima tan extremo. Por lo tanto, se da lugar a un consumo excesivo de energía por parte del compresor del acondicionador de aire convencional. La combinando del efecto del serpentín de enfriamiento convencional junto con la carga de enfriamiento total de enfriamiento por evaporación sobre los serpentines se puede disminuir y simultáneamente mejorar la tasa de rechazo de calor en el condensador [13]. Es así que, un sistema de recuperación de calor residual presenta una gran oportunidad para recuperar energía de los sistemas de refrigeración in situ y convertirla en agua caliente útil, dependiendo de la aplicación se pueden obtener ahorros [14]. Estos sistemas de recuperación de calor además se pueden adecuar a distintos sistemas como el sistema con enfriamiento evaporativo de ciclo M [15].
Una inclusión interesante en la evaluación de la recuperación de la energía residual es abordada en este trabajo de investigación. Primeramente, describimos a las plantas textiles para realizar el estudio de caso que será abordado. Segundo, realizamos la metodología numérica para evaluar los estudios de casos, considerando en este apartado el dominio de análisis, las condiciones de frontera, el modelo de turbulencia, discretización del dominio y los parámetros de estudio. Este estudio es realizado en la ciudad de Lima, Perú, en las instalaciones de una planta textil industrial dedicada al proceso de producción de hilados de algodón, con una producción promedio de 12 toneladas por periodo de 24 horas de operación. Para evaluar el rendimiento del consumo energético en una planta textil no solo se toma en cuenta el consumo de equipos, sino que también costes de iluminación y climatización por lo que una modificación en este tipo de sistema no solo puede mejorar el consumo energético, sino que también mejorar el confort. Por último, verificamos los resultados numéricos para el caso de estudio mencionado mediante el análisis de un ambiente comparando al diseño de una planta textil que hace recircular el aire por el entretecho y controla la temperatura del ambiente de trabajo, aliviando de esta forma la carga térmica que se genera por los diferentes equipos de hilandería.
2. DESCRIPCIÓN DE LAS PLANTAS TEXTILES
2.1. Unidad de análisis
El estudio se realiza sobre una planta industrial textil de procesos de hilandería de algodón y mezclas, ubicada en Lima, Perú. El diseño considera la investigación documental y de campo para el caso materia del estudio.
El propósito del estudio es la posibilidad de implementar la propuesta para mejorar la viabilidad económica de una nave usada en la industria textil, esto dado que se necesita que el ambiente interior en todas las dependencias de la planta sea el adecuado, tanto para las personas como para las máquinas o elementos de la producción, temperatura, humedad relativa, ruido, composición del aire, etc.
En todos los casos, el ingreso de aire del entretecho a la sala principal es indeseable, por lo que es un factor bastante controlado en la planta. En condiciones de baja temperatura, la perdida de calor a través del cielo raso puede representar puntos o zonas frías que influyen de forma negativa en el proceso productivo; en condiciones de alta temperatura, el ingreso de calor del entretecho representa una sobrecarga para el sistema de climatización y su consecuente costo en energía. En ambos casos el ingreso de aire del entretecho puede representar además ingreso de polvo, suciedad, pelusa, fibras cortas, trazas de lana de vidrio del aislamiento, o cualquier otro elemento que se pudiera encontrar, por todo lo precedente, para el efecto de la simulación, resulta valido considerar que existe un sellado lo suficientemente confiable del cielo raso, de forma que no permite ni permitiría el paso de aire desde el entretecho a la sala, ni en la condición real ni en la simulada, siendo el cielo raso capaz de manejar el flujo y la presión adicional que se generará en el entretecho por el paso del aire residual
2.2. Consideración funcional
El estado natural del sub ambiente correspondiente al entretecho es de aire estanco, ya que no presenta sistemas de climatización ni de circulación de aire, además de encontrarse idealmente cerrada en todas sus superficies. En la práctica, por el efecto de la climatización de la sala principal que trabaja en presión menor a la atmosférica y por la reducida inversión en el sellado del cerramiento lateral del entretecho, el entretecho se convierte en un pulmón abastecedor de aire a la espera de la menor abertura en el cielo raso para ingresar a la sala principal. Ese efecto es factible de ser controlado asegurando un adecuado nivel de sellado del cielo raso, adicionalmente es usual la instalación de “extractores” eólicos que contribuyen en menor medida a reducir la presión del aire del entretecho sobre su intención de ingresar a la sala principal.
2.3. Propuesta
Analizar el sistema de ahorro de energía para plantas de procesos textiles aprovechando la energía contenida en el aire residual del sistema de climatización. Esto permitirá reducir la demanda de energía eléctrica necesaria para controlar las condiciones de temperatura del ambiente en las instalaciones materia del estudio, significando reducción de costos en la operación del sistema de climatización, sistemas auxiliares y consecuentemente en la operación de la planta, esto se muestra en la Fig. 1.
3. METODOLOGÍA NUMÉRICA PARA PLANTAS TEXTILES
A partir de la descripción problemática es identificado el caso de estudio de una planta textil con aprovechamiento de la energía residual y es elaborado la metodología numérica para una sección de la planta textil con la condición real y la propuesta. Además, se establece que la “nave industrial” o sección de la planta textil, consta de dos sub ambientes, el principal que es la nave propiamente, que inicia a nivel del suelo y en el que se encuentran el personal, las máquinas y se realiza el proceso productivo, y un segundo sub ambiente, que corresponde al que se encuentra entre el cielo raso o falso techo y el techo. A seguir describiremos con énfasis la propuesta numérica para la realización de la simulación de estos casos de estudio.
3.1 Primer caso de modelamiento
Este primer caso define una condición real o existente del sistema de climatización sin aprovechamiento de la energía residual. El modelamiento considera una sección de la planta de hilandería, la cual consta de los equipos de planta y los ductos de salida de aire y rejillas de ventilación, ver Fig. 2. Además, no es considerada toda la planta de hilandería sino este modelo observado en la Fig. 2 es catalogada como la sección representativa de la planta.
La Fig. 2 muestra la sección representativa de la geometría de la planta de hilandería, considerando ductos de salida de aire de forma externa con rejillas y equipos textiles igualmente espaciados.
3.2 Segundo caso de modelamiento
Corresponde a la condición propuesta simulada en el modelamiento, en la que, el aire residual del sistema de climatización, y que normalmente es expulsado al exterior, sea forzado a través del segundo sub ambiente que se encuentra en el área entre el cielo raso o falso techo y el techo. Es importante también considerar que, dentro del espacio entre los techos, se encuentran los conductos que trasladan el aire climatizado a inyectar en la sala, por lo que la superficie de los conductos (que son de plancha metálica de e<1mm revestida en lana de vidrio de e<25mm), representa una considerable área de transferencia de calor entre el aire que conduce y el circundante en el entretecho. Es importante considerar que, debido a la configuración del sistema de climatización, es usual que el aire expulsado del sistema, se dispone en la parte alta de la planta, sobre la línea de altura del cielo raso o, en el extremo, un poco por encima de la misma, lo que deja el aire residual y la energía residual térmica y cinética que éste contiene justo a disposición para ingresar o ser direccionada al entretecho. La Fig. 3 muestra una configuración similar a la planta convencional, pero agrega un ducto de recirculación del aire y dos salidas.
3.3 Condiciones de frontera
De acuerdo con las condiciones de estudio son establecidas algunas condiciones de contorno que ayudarán a determinar las condiciones reales a las cuales la planta de hilandería opera. La tabla I y tabla II muestran estás condiciones de contorno. Además, la Fig. 4 y Fig. 5 especifican de manera esquemática las condiciones de contorno de los casos de estudio propuestos.
Superficie seleccionada | Tipo de contorno | Etiqueta |
Entrada de aire fresco | Ingreso | INLET |
Salida de aire desde el sistema de | Salida | OUTLET |
Superficie intermedia del dominio (plano Y-Z) | Simetría | SYMMETRY1 |
Superficie intermedia del dominio (plano Y-Z) | Simetría | SYMMETRY2 |
Superficie seleccionada | Tipo de contorno | Etiqueta |
Entrada de aire fresco | Ingreso | INLET |
Salida de aire desde el sistema de | Salida | OUTLET1 |
Salida de aire desde el entretecho | Salida | OUTLET2 |
Superficie intermedia del dominio (plano Y-Z) | Simetría | SYMMETRY1 |
Superficie intermedia del dominio (plano Y-Z) | Simetría | SYMMETRY2 |
3.4 Modelo de turbulencia
Los modelos de turbulencia de las ecuaciones permiten la determinación tanto de una longitud turbulenta y una escala de tiempo resolviendo dos ecuaciones de transporte separadas. El modelo estándar en Ansys Fluent se encuentra dentro de esta clasificación y se ha convertido en uno de los más usados en el campo de la ingeniería, desde que fue propuesto por Launder y Spalding [16]. Particularmente, el modelo estándar está fundamentado en las ecuaciones de transporte, para la energía cinética turbulenta y su tasa de disipación [16].
La energía cinética turbulenta, y su tasa de disipación, se obtienen a partir de la ecuación (1) y la ecuación (2) de transporte.
Donde, es la generación de energía cinética turbulenta debido a los gradientes de la velocidad media, es la generación de energía cinética turbulenta debido a flotabilidad, , es la contribución a la disipación de la turbulencia debido a la fluctuación en la dilatación de turbulencia compresible, son constantes del modelo, son los números de Prandtl para k y . son los términos fuente de ambos parámetros.
La viscosidad turbulenta para este modelo es descrita por la combinación de k y como se muestra en la ecuación (3).
3.5 Discretización del dominio
Para calcular las propiedades de los campos involucrados en el fenómeno estudiado se discretiza el dominio en volúmenes de control, donde en cada una de las celdas generadas se resuelven las ecuaciones de conservación de la masa, cantidad de movimiento y energía.
La relación de aspecto en Fluent es una medida del estiramiento de una celda. Se calcula como la relación entre el valor máximo y el valor mínimo de cualquiera de las siguientes distancias: las distancias normales entre el centroide de la celda y el centroide de la cara, y las distancias entre el centroide de la celda y los nodos, esto se ve en la Fig. 6, [17].
Para los elementos del tipo poliedros la ortogonalidad se mide de la siguiente manera:
El producto escalar normalizado del vector de área de una cara y un vector desde el centroide de la celda al centroide de esa cara
El producto escalar normalizado del vector de área de una cara y un vector desde el centroide de la celda al centroide de la celda adyacente que comparte esa cara
El valor mínimo que resulta de calcular la ecuación (4) y ecuación (5) para todas las caras se define entonces como la ortogonalidad de la celda. Por lo tanto, las peores celdas tendrán una ortogonalidad más cercana al valor de cero y las mejores celdas tendrán una ortogonalidad más cercana al valor de uno. Los vectores usados para este cálculo se aprecian en la Fig. 7.
La discretización del dominio consta de un tipo de malla poliédrica ya que proporciona un mayor grado de precisión en la obtención de los campos de temperatura y velocidades del flujo circundante en todo el sistema de climatización. En los modelos propuestos se obtienen 832,498 elementos poliédricos para la malla correspondiente al módulo de análisis del sistema de climatización convencional, y 837,794 elementos poliédricos para la malla del módulo del sistema de climatización propuesto. La Fig. 8 y Fig. 9 muestran el mallado de los sistemas de climatización de las plantas considerando elementos poliédricos.
Se debe tener un equilibrio entre la cantidad de elementos utilizados y el tiempo computacional empleado, en el cual influye el tipo de elementos usados y el tamaño de los mismos en el enmallado. Sin embargo, existen métricas a revisar sobre esta discretización para obtener buenos resultados, estas son la ortogonalidad y la relación de aspecto. Las métricas mínimas de ortogonalidad se aprecian en la Fig. 10 y Fig. 11, mientras que la cantidad de elementos y datos de la relación de aspecto de los elementos se muestran en la Tabla III y Tabla IV.
Numero de celdas | Tipo de elementos | Ortogonalidad mínima | Ratio de aspecto máximo |
832498 | Poliedros | 0.13655291 | 48.654888 |
Numero de celdas | Tipo de elementos | Ortogonalidad mínima | Ratio de aspecto máximo |
837794 | Poliedros | 0.05002873 | 88.102013 |
3.6 Parámetros para estudio
Debido a la complejidad de los casos de estudio, fueron considerados los fenómenos de convección y radiación, las cuales son dadas para las condiciones de contorno mencionadas en la Tabla V. También, las propiedades térmicas como la lana de vidrio dependen de la humedad que esta posea, entre otras variables [18], para este caso de estudio se toma de referencia el valor de 0.03 W/m-K en las zonas que la contengan. La Tabla V, muestra los parámetros insertados en el modelo numérico para representar los fenómenos mencionados. La Tabla VI, muestra las velocidades y propiedades establecidas en la entrada y salida para ambos casos de estudios.
Frontera | Temperatura (°C) | Espesor (m) | Coeficiente de conducción (W/m-K) | Coeficiente de convección (W/m2-K) | Emisividad | Generación de calor (W/m2) |
Ductos | 24 | 0.013 | 0.03 | |||
Rejillas | 32 | 0.013 | 0.03 | |||
Paredes | 32 | 0.12 | 17 | |||
Techo | 40 | 0.001 | 202.4 | 20 | 0.7 | |
Máquinas | 2.8 |
4. RESULTADOS
A partir de la metodología numérica propuesta anteriormente se realizó x’mediante simulaciones numéricas el comportamiento del sistema de climatización del aire convencional y propuesto para una planta de hilandería, son obtenidos los campos de temperaturas, denominadas iso-temperaturas y líneas de corriente para representar las velocidades del flujo de aire.
La Fig. 12 muestra el cerramiento o cobertura del techo de la planta industrial con un valor de temperatura de 40 °C uniformemente representada en toda su superficie, que corresponde a una condición real en temporadas de verano, y que constituye una de las situaciones propuestas en el modelamiento numérico.
La Fig. 13 evidencia un cambio en las temperaturas de la cobertura del techo para el modelamiento propuesto, este cambio es ligero mostrando el cambio del campo de temperaturas en gran medida en la superficie del techo de la planta textil. Complementariamente, en la Fig. 14 y Fig. 15 se muestran las líneas de corriente con sus respectivas velocidades, donde se evidencia el aire que se hace recircular, tanto para el modelo real y propuesto, respectivamente.
CONCLUSIONES
Se puede aprovechar la energía residual cinética y térmica contenida en el aire residual como un factor de reducción de demanda de energía eléctrica del sistema de climatización de la planta materia del estudio.
El modelamiento permite, para el caso de estudio, identificar el comportamiento del aire climatizado en la planta industrial, que por las velocidades en la entrada muestra un comportamiento turbulento.
La implementación de la propuesta permitirá reducir la carga térmica y la magnitud del rango de temperatura de operación de los materiales que conforman el cerramiento, el aislamiento y otros activos que se pudieran encontrar en el entretecho de la nave industrial estudiada.
El uso de elementos poliédricos es factible en simulaciones de ambientes cerrados y geometrías parecidas el caso estudiado, dado que la calidad del mallado en general es muy buena y no incrementa la cantidad de elementos requeridos como si sería si se usasen tetraedros.
El ducto de recirculación de aire cumple su función y existe un aprovechamiento de la energía residual por el flujo que ingresa al entretecho.