ARTÍCULOS

DOI: 10.22395/rium.v15n28a12

Evaluación del proceso de producción de detergente en polvo a partir de la simulación

Evaluation of the powder detergent production process from simulation

Junior Lorenzo Llanes^*; Lourdes Zumalacárregui de Cárdenas^**; Osney Pérez Ones^***

^* Ing., Profesor Asistente, Miembro del Grupo de Análisis de Procesos, Facultad de Ingeniería Química, Instituto Superior Politécnico José Antonio Echeverría, Cujae. Calle 114, No. 11901, e/ Ciclovía y Rotonda, Marianao, La Habana, Cuba, e-mail: junior@quimica.cujae.edu.cu

^** Ph. D., Profesora Titular, Consultante y de Mérito, Jefe del Grupo de Análisis de Procesos, Facultad de Ingeniería Química, Instituto Superior Politécnico José Antonio Echeverría, Cujae. Calle 114, No. 11901, e/ Ciclovía y Rotonda, Marianao, La Habana, Cuba, e-mail: lourdes@quimica.cujae.edu.cu

^*** Ph. D., Profesor Auxiliar, Miembro del Grupo de Análisis de Procesos, Facultad de Ingeniería Química, Instituto Superior Politécnico José Antonio Echeverría, Cujae. Calle 114, No. 11901, e/ Ciclovía y Rotonda, Marianao, La Habana, Cuba, e-mail: osney@quimica.cujae.edu.cu

Recibido: 13/01/2015

Aceptado: 11/12/2015

RESUMEN

En Cuba la producción de detergente en polvo se lleva a cabo por una sola fábrica. El presente trabajo tuvo como objetivo evaluar el proceso de producción de detergente en polvo a partir de la simulación. Para ello se empleó el simulador de procesos químicos HYSYS versión 3.2, siendo necesaria la determinación de propiedades de componentes hipotéticos que no forman parte de la base de datos de componentes puros. Los resultados del caso base mostraron ineficiencias en la operación del proceso, las más importantes, en la etapa de neutralización y en el calentamiento de los gases de secado. Los resultados del caso estudio 2 muestran que es posible reducir el pH de la pasta neutra en un 15,34 % respecto al caso base, para una producción de 848 kg_{pasta neutra}/h, sin afectar el porcentaje de ingrediente activo. Por otra parte, una reducción en el consumo de gas licuado de 1.500 L/d implica un aumento de las ganancias por valor de 227.301 $/año con respecto al caso base.

PALABRAS CLAVE

producción de detergente, simulación, HYSYS.

ABSTRACT

In Cuba, powder detergent production is carried out by a single factory. This work research was intended to evaluate the powder detergent production process from a simulation process. For this purpose, the HYSYS v. 3.2 chemical process simulator was employed and definition of properties of hypothetical components that are nor part of the pure components database was required. Results of the basic case showed inefficiencies of the process operation, the most important being during the neutralization stage and the drying gas heating stage. Results of Study Case 2 show that it is possible to reduce pH of the neutral paste to 15.34% when compared to the basic case, for a production of 848 kg neutral paste per hour, without affecting the percentage of the active ingredient. Additionally, a decreased consumption of liquefied gas of 1,500 L/d implies increased profit in the amount of $ 227,301 per year, compared to the basic case.

KEY WORDS

Detergent production; simulation; HYSYS.Key words: production of detergent, simulation, HYSYS

INTRODUCCIÓN

En Cuba, la producción de detergente en polvo se lleva a cabo por una sola fábrica. La misma se enfrenta a la necesidad de elevar su eficiencia y eficacia con el objetivo de aplicar los lineamientos de la política económica y social del Partido y la Revolución que le permitan insertarse exitosamente en el nuevo modelo económico cubano. Actualmente, la carencia de instrumentos de medición, así como de sistemas automatizados de control en dicha entidad impide una adecuada evaluación del mismo.

Lo anterior repercute de forma negativa en la toma de decisiones para incrementar los rendimientos. Por ello el presente trabajo tuvo como objetivo evaluar el proceso de producción de detergente en polvo a partir de la simulación. Según [1], con el propósito de entender el comportamiento de un sistema o de evaluar varias estrategias para la operación del mismo, se acude a la simulación.

El simulador de procesos HYSYS versión 3.2 es uno de los simuladores más respetados en la industria. El mismo integra la simulación en estado estacionario y la simulación dinámica [2].

El proceso de producción de detergente en polvo consta de tres etapas fundamentales: neutralización, en la cual se neutraliza el ácido dodecil-benceno sulfónico lineal (ADBS) con hidróxido de sodio y se añade agua, obteniéndose una pasta neutra cuyo ingrediente activo es el dodecil-benceno sulfonato de sodio (LAS), que pasa a un tanque de maduración para la culminación de la reacción. Este tanque está conectado a una bomba de recirculación que mantiene la pasta en movimiento.

En este paso se añade agua, la cual debe ser controlada para no afectar la humedad de la pasta neutra [3].

Si la pasta se considera conforme, según las especificaciones del proceso, pasa a la siguiente etapa, la batición. Durante la batición se mezcla la pasta neutra con el resto de las materias primas que componen el detergente y que le dan las propiedades necesarias y mejoran su calidad. La batidora carga automáticamente los materiales, con excepción del aditivo sólido que lo añade manualmente el operador auxiliar que descarga las materias primas. Según las normas de buenas prácticas se deben tomar muestras de la pasta de batidora en el tanque de alimentación como mínimo dos veces en el turno de trabajo y enviarlas al laboratorio para determinar el porcentaje de humedad cuya especificación es 28-32 % [4].

La tercera etapa es el secado. En esta etapa la pasta húmeda, salida del proceso de batición, se atomiza empleando una bomba de alta presión. En estas condiciones se alimenta la torre de secado a contracorriente con aire caliente a 380°C aproximadamente. A la torre se alimenta un volumen de aire caliente suficiente para suministrar el calor necesario para completar la evaporación del agua. La transferencia de masa y energía se logra mediante el contacto directo del gas caliente con las gotas dispersas. A la salida de la torre, el gas que sale por la parte superior se enfría (para no dañar directamente la biodiversidad) y las partículas finas se separan del gas por medio de un filtro (para que las partículas finas no pasen a la atmósfera contaminándola). La torre de aspersión generalmente está construida con chapa de acero común; está provista de mirillas o puntos de observación e inyectores de pulverización. Una vez obtenido el polvo, se le atomiza perfume y se obtiene el polvo de detergente que se usa en los hogares, hospitales, industrias, etc.

1. MATERIALES Y MÉTODOS

Para llevar a cabo la evaluación del proceso de producción de detergente en polvo se empleó la simulación como herramienta. Para ello se utilizó el simulador de procesos HYSYS versión 3.2, dado que el mismo dispone de rigurosos modelos para la evaluación de propiedades físicas y termodinámicas y una base de datos de cerca de 1.500 componentes que permite el cálculo de propiedades de compuestos reales e hipotéticos; cuenta con un amplio conjunto de operaciones unitarias y otros utilitarios [2].

En la base de datos del simulador HYSYS versión 3.2 no aparecen varios de los compuestos que intervienen en el proceso; por ello fue necesario crearlos como componentes hipotéticos. Los componentes creados, así como las propiedades que se tuvieron en cuenta, se muestran en la tabla 1.

Para la determinación de las capacidades calóricas del ADBS y del LAS se utilizaron dos métodos de contribución de grupos dados por [8].

Método 1: Ley de Kopp's modificada por Werner en 1941 para sólidos. En este caso las capacidades calóricas del compuesto se obtienen a partir de la sumatoria del producto del número de átomos de cada elemento presente por un factor empírico para cada elemento, el cual depende del estado de agregación del compuesto en cuestión. En la tabla 2 se muestran los resultados de la aplicación del método.

Método 2: contribución de grupos de Chunch y Swanson (1973) para líquidos orgánicos. En este caso la capacidad calórica del compuesto se obtiene a partir de la sumatoria de la contribución de un factor reportado para cada grupo funcional (CH₃, CH₂, CH), por el número de veces que aparece dicho grupo en la estructura. Este factor depende del estado de agregación del compuesto en cuestión. En la tabla 3 se muestran los resultados de la aplicación del método.

La búsqueda de buenos valores para propiedades físicas de los parámetros es la clave para una simulación exitosa. Lo anterior depende mucho de escoger el método correcto de estimación de propiedades. Los cuatro factores que se deben tener en cuenta al escoger los métodos de estimación de propiedades son: la naturaleza de las propiedades de interés; la composición de la mezcla; los intervalos de presión y temperaturas, y la disponibilidad de parámetros.

Dada la presencia de electrolitos en el proceso de producción de detergentes se hace necesaria la selección de un paquete de propiedades que se ajuste a esta situación. Los métodos más comunes para el trabajo con electrolitos son el Pitzer y el NRTL-modificado [9]. En el presente trabajo, teniendo en cuenta la presencia de compuestos orgánicos y agua, se empleó el paquete NRTL-modificado, ya que la ecuación de Pitzer no considera la interacción de las partículas con las moléculas de solvente [10].

En la tabla 4 se muestran los módulos empleados para simular las diferentes operaciones unitarias presentes en el proceso, así como los equipos asociados a estas.

A continuación se describen los parámetros que se fijaron para la simulación de cada una de las operaciones unitarias.

1.1.1 Neutralización

Para la simulación de esta operación se crearon tres corrientes de entrada con sus correspondientes variables de decisión, todas a 25°C y 101,3 kPa, que corresponden a los reactivos (ADBS y NaOH(ac)), así como una corriente de agua. Las especificaciones se muestran en la tabla 5. Por otra parte se creó una corriente de salida (LAS) con una temperatura de 55°C y presión de 101,3 kPa. También fue añadida una corriente de energía (Qperdido) de manera que el simulador determinase la energía cedida en el proceso. La corriente de salida del neutralizador constituye la entrada a un segundo reactor de conversión el cual es el equivalente al tanque intermedio descrito en el proceso. En este caso solo se fijó la presión a la corriente de salida.

Como uno de los parámetros a controlar en la etapa es el pH, fue necesario crear una hoja de cálculo (spreadsheet) para conocer su valor. El cálculo se realizó a partir de la diferencia entre la constante de disociación del agua, cuyo valor se tomó a la temperatura de 60°C que, según [11], tiene un valor de 9,6×10^-14 y el pOH de la pasta neutra, valor que se calcula a partir de la concentración de NaOH, dado que este es un electrolito fuerte y soluble. La concentración molar de NaOH se determina a partir del flujo molar de NaOH (n(NaOH)) y el flujo volumétrico (Vsalida) a la salida del reactor. Para ello se emplearon las ecuaciones desde la (1) hasta la (3). La ecuación 3 está afectada por el valor de la constante de disociación del agua a las condiciones de trabajo.

(1)

(2)

(3)

1.1.2 Batición

Para la simulación de esta operación se tienen dos corrientes de entrada. Una es la salida del tanque intermedio y otra se creó teniendo en cuenta los consumos de materias primas en la etapa (aditivos). La temperatura y presión de la corriente de Aditivos es 25°C y 101,3 kPa, respectivamente. En la tabla 6 se muestran las especificaciones de la corriente creada.

También se creó una corriente de salida (pasta batidora) a la cual se le fijó la temperatura de 45 °C valor que se obtiene de mediciones realizadas durante el proceso productivo; por otra parte, se conectó una corriente de energía (W) de manera que el simulador evaluase el trabajo entregado por el agitador para el mezclado.

1.1.3 Horno

Para la simulación de esta sección se crearon dos corrientes de entrada al reactor cuyas especificaciones se muestran en la tabla 7. Además de las corrientes de entrada se creó una de salida (gas secado) donde las únicas especificaciones hechas fueron las de temperatura y presión (400°C y 101,3 kPa). Al igual que en el resto de las operaciones descritas anteriormente se conectó al reactor una corriente de energía (Q) de manera que se calculara el calor cedido al medio.

1.1.4 Secado

En el proceso de secado están presentes las operaciones de transferencia de calor y masa. Para la simulación de la operación de intercambio de calor a partir del uso del módulo intercambiador de calor se crearon dos corrientes de materiales a la salida, pasta a secador y aire a secador, de las cuales solo a la corriente pasta a secador se le especificó la temperatura (90°C). Las corrientes de entrada fueron calculadas por el módulo que representa el horno y la bomba de alta.

Por otra parte, se consideró que no existía caída de presión ni por el lado de los tubos ni de la coraza. En el caso de la operación de transferencia de masa se empleó el módulo divisor de flujo, considerando la operación adiabática. Para ello se crearon dos corrientes de materiales (detergente y emisión) a la salida, así como una de energía (Q₁). Para determinar las composiciones de las corrientes de salida se tuvo en cuenta que la humedad del producto final debía ser de 2,5 % como máximo.

1.1.5 Bombeo

El sistema de bombeo está formado por dos bombas (alta y baja); para la simulación se tomó el módulo bomba. A cada uno de los módulos se le conectó una corriente de energía, potencia 1 y potencia 2, respectivamente, de manera que el simulador evaluase la potencia consumida. En la figura 1 se muestra el esquema de simulación del proceso de producción de detergente en polvo.

2. RESULTADOS DEL CASO BASE (CB)

2.1 Balance de materiales y energía en la sección de neutralización

En la tabla 8 se muestran los resultados del balance de materiales en el neutralizador y el cálculo del pH obtenido manualmente y a partir del simulador. Los cálculos manuales se realizaron con el objetivo de validar el modelo de simulación.

De este resultado se aprecia que el esquema de simulación reproduce con menos del 5 % de error relativo el balance de materiales en el neutralizador.

Para el cálculo de las pérdidas de energía en el neutralizador se estimó la energía cedida por convección y radiación (Qconv, Qrad). Estos valores se muestran en la tabla 9. No fue posible verificar el resultado de los balances de energía con el simulador, por no haber posibilidades de realizar mediciones experimentales que permitieran corroborar los estimados de propiedades energéticas de los componentes hipotéticos.

2.2 Balance de energía en la sección de batición

El resultado más importante de la etapa de batición es que la potencia requerida para realizar el mezclado de los componentes, en las condiciones de operación reportadas por la fábrica, a partir de HYSYS, es de 47,11 kW. Dado que la potencia máxima del motor del agitador es de 55 kW, operando al 86 % de eficiencia puede suministrar el valor requerido.

En la tabla 10 se presenta el resultado del balance de materiales en el horno.

De este resultado se aprecia que el esquema de simulación reproduce con menos del 5 % de error el balance de materiales en el horno.

Para el balance de energía, el cálculo del calor cedido al medio durante la combustión del gas en el horno es 4,46×10⁶ kJ/h, que representa que el 80 % está siendo cedido al medio lo que evidencia la baja eficiencia del sistema. Una manifestación de esta ineficiencia en el intercambio de calor en el horno es la temperatura de salida de los gases de combustión, muy por debajo de la media para una reacción de combustión que se desarrolle correctamente.

2.4 Balance de materiales y energía en la torre de secado

Para verificar la validez de representar la torre de secado por el módulo Divisor de flujo se calculó el agua evaporada en el proceso (973,73 kg/h) y se comparó con la reportada por el simulador (976,97 kg/h), de donde el error relativo en la estimación es menor que 5 %.

Adicionalmente se evaluó manualmente la energía requerida para el secado y el porcentaje de pérdidas en el secador, arrojando 0,77×10⁶ kJ/h y 23,32 % respectivamente.

En la tabla 11 se presenta en resultado de la evaluación del divisor de flujo según HYSYS.

Si se hubiera considerado un porcentaje de pérdidas de calor en la torre de secado igual al obtenido mediante el cálculo manual, la energía requerida para el secado hubiera sido de 2,77×10⁶ kJ/h, valor que difiere del calor absorbido hallado manualmente en un 9,3 %. Lo anterior indica que es posible emplear para la simulación este módulo de cálculo.

3. DISCUSIÓN DE CASOS DE ESTUDIO (CE)

A partir de los resultados del caso base se determinó realizar tres estudios de casos con el objetivo de valorar alternativas tecnológicas que conllevan al uso racional de materias primas y energía.

Los resultados del caso base muestran que el pH de la pasta neutra, como el porcentaje de componente activo (% LAS) se encuentran fuera de especificaciones; de ahí que se estudiaron dos variantes tecnológicas para mejorar estos resultados. La variante tecnológica analizada consistió en disminuir el flujo de la solución de hidróxido de sodio alimentada al neutralizador hasta 82 kg/h.

Para ello a la entrada del neutralizador se debe regular el flujo de solución de hidróxido de sodio en el cabezal de la bomba triple que corresponde a este reactivo. En la variante se analiza cómo la variable independiente (flujo másico de NaOH) influye sobre las variables dependientes (pH y fracción másica de componente activo). En este caso la variable independiente varía en un intervalo cuyo límite inferior es de 82 kg/h y el superior de 92 kg/h, tomándose un paso de 1,00 kg/h. Los resultados de este caso de estudio se muestran en la figura 2.

A partir del análisis de la figura 2 se puede observar que al disminuir el flujo de hidróxido de sodio desde 92 kg/h hasta 82 kg/h el pH de la pasta neutra disminuye desde 12,15, valor que se encuentra fuera de especificación según la norma de neutralización [4], hasta 10,60, lo que representa una disminución en el valor del pH de un 15,63 %. A partir de esta modificación el exceso de hidróxido de sodio disminuye desde 13,25 % hasta 0,3 %. Lo anterior se debe a que al disminuir la concentración de hidróxido de sodio a la salida el valor del pOH aumenta lo que provoca una disminución del pH de la pasta neutra.

Por otra parte, se observa un aumento lineal del porcentaje de componente activo (C₁₈H₂₉SO₃Na) desde un 39,40 % hasta un 39,85 %, resultado que se corresponde con el balance de materiales. No obstante ambos valores se hallan fuera de especificación que según [4] debe ser (38,2 % ± 1 %). Para un flujo de 85 kg/h de hidróxido de sodio se obtiene una mejor condición de operación, manteniendo constante el flujo de ácido al neutralizador.

Al igual que en el CE1 en la variante que se analiza a continuación se pretende mejorar el resultado obtenido para el pH y el porcentaje de componente activo. En este caso se disminuyeron los flujos de reactivos. La variante tecnológica analizada consistió en disminuir el flujo de la solución de hidróxido de sodio hasta 75 kg/h, así como el de ADBS hasta 305 kg/h.

Para ello, a la entrada del neutralizador se deben regular los flujos de solución de hidróxido de sodio y ADBS en los cabezales de la bomba triple que corresponden a estos reactivos. En la variante se analiza cómo las variables independientes (flujo másico de NaOH(ac) y de ADBS) influyen sobre las variables dependientes pH y fracción másica de componente activo. En este caso las variables independientes varían en los siguientes intervalos: flujo másico de NaOH de 70 kg/h-90 kg/h y flujo másico de ADBS entre 300 kg/h-330 kg/h; para ambos el paso es 5,00 kg/h.

En este caso se produce una reducción del pH de la pasta neutra desde 12,25 hasta 10,37, lo que representa una reducción de un 15,34 %. Por otra parte, el menor flujo de ácido e hidróxido provoca una disminución del porcentaje de componente activo en la pasta neutra desde un valor de 39,23 %, fuera de especificación, hasta 37,98 % valor que se encuentra dentro de la norma.

Al comparar estos resultados con los del CE1 se ve que el pH se reduce en ambos casos prácticamente en el mismo porcentaje, aunque en el CE2 adicionalmente se disminuye el consumo de ácido. La producción de pasta neutra para las mejores condiciones del CE2 (75 kg/h de hidróxido y 305 kg/h de ácido) es de 848 kg_{pasta neutra}/h, mientras que en las condiciones del CB es de 885 kg_{pasta neutra}/h, con una disminución de un 4,2 %. En la tabla 12 se muestra el comportamiento de las variables en el CE2.

3.3 Caso de estudio 3 (CE3)

El elevado valor de las pérdidas energéticas en el horno hace necesario el estudio de alguna alternativa que permita disminuir este valor.

Variante tecnológica analizada: disminuir el flujo de gas licuado para la combustión en el horno hasta 2.000 L/d.

En este caso se analiza la influencia de la variable independiente flujo de gas licuado sobre la variable dependiente calor perdido (Qp H). Al igual que en el caso base se mantiene constante la temperatura de salida de los gases de combustión. Para la simulación de este caso de estudio se varía el valor de la variable independiente desde los 0,095 kmol/h hasta los 0,22 kmol/h, tomándose un paso de 0,01kmol/h.

Una disminución del flujo de gas licuado desde 3 500 L/d (0,219 m³/h) hasta 2.000 L/d (0,125 m³/h) provoca una disminución en el calor cedido al medio por el horno, manteniendo constante el flujo de aire al horno, así como la temperatura de salida de los gases de combustión, desde 4,5×10⁶ kJ/h hasta 1,7×10⁶ kJ/h, lo que equivale a una reducción de un 59,57 %. No obstante, las pérdidas generadas continúan siendo grandes representando el 31,66 % del calor total que se evoluciona durante la combustión.

4. ANÁLISIS ECONÓMICO

Para llevar a cabo el análisis económico se determinó para cada alternativa la ganancia económica (GE), así como la relación beneficio costo (B/C) tomando como punto de partida los costos por concepto de materias primas y utilidades, así como los ingresos por concepto de detergente. Los resultados se muestran en la tabla 13.

Para profundizar en el análisis económico, en la alternativa de mayor ventaja, se realizó un análisis de sensibilidad de la ganancia ante la variación de los ingresos por venta de detergente y los egresos por concepto de ADBS e hidróxido de sodio. Para ello se definió como variable dependiente la ganancia y se consideraron cinco valores de cada variable independiente: el valor esperado (0 %), un valor optimista (+5 %), otro pesimista (-5 %) y valores extremos optimista (+10 %) y extremos pesimista (-10 %). Puede observarse en la figura 3, que la ganancia es sensible a todas las variables independientes analizadas, pero la mayor sensibilidad aparece ante la variable ingresos por venta de detergente que obtengan.

5. CONCLUSIONES

Se comprobó que el esquema de simulación desarrollado con el simulador de procesos químicos HYSYS arroja resultados que se corresponden con las condiciones reales de operación en el proceso de producción de detergentes en polvo, desde el punto de vista de los balances de materiales, lo que permite su uso para el análisis de alternativas de operación.

Por otro lado, los resultados obtenidos de los casos de estudio 1 y 2 muestran que es posible disminuir el pH de la pasta neutra en un 15,5 %, aproximadamente, aunque con la variante aplicada en el caso de estudio 2 se reduce la producción en un 4,2 %, mientras que se logra que el porcentaje de componente activo sea el exigido en la especificación. Lo anterior permite reajustar las condiciones de operación de modo que se eleve la productividad del proceso.

REFERENCIAS

1 González, C. E. ''Simulador de una columna de destilación binaria''. Tesis en opción del título de Ingeniero Industrial y en Sistemas Organizacionales, Escuela de Ingeniería, Universidad de La Salle. México DF. 2001.

2 Hyprotech. ''Hyprotech: simulation software for industry'' Disponible en internet: http://www.ucalgary.ca/community/research/hyprotech. [Citado el 23 de febrero de 2013].

3 Sardiñas, D. J. ''Desarrollo de una base de conocimientos para el entrenamiento de los operarios y tecnólogos de la planta de detergentes de la Empresa Suchel Debón y propuesta de un sistema de control para el proceso de neutralización''. [Tesis MSc], Instituto Superior Politécnico José Antonio Echeverría. 30 pp. 2012.

4 Reglamento técnico de la empresa, Pto: 8226-01, Detergentes en polvo multiusos. Neutralización y batición. Proceso tecnológico operativo. Empresa Suchel Debón, 2013.

5 Initial assessment report for 20th SIAM, UNEP Publications. Disponible en internet: jheinze@johnadams.com [Citado el 25 de febrero de 2013].

6 Adalid, G. N. ''Diseño de una planta para la fabricación de detergente'', Proyecto de fin de carrera, Universidad de Barcelona. 2001.

7 Perry, H. R; Green, W. D. ''Perry's Chemical Engineers' Handbook''. McGraw-Hill. New York. Section 2. 2008.

8 Sinnott, R. K. ''Chemical Engineering Design''. Vol 6, 4th Edition. New Jersey. Elsevier. 1031pp. 2005.

9 Carlson, E. C. ''Don't gamble with physical properties for simulations''. Chemical engineering progress, October pp35. 1996.

10 Aspentech. ''Aspen HYSYS Simulation Basis'' version 7. Aspen Technology, Inc. USA. 561 pp. 2008.

11 Petrucci, R. H; Harwood, W. S. ''General Chemistry: principles and modern applications''. Prentice-Hall. New Jersey. 587 pp. 1999.