Modelagem, simulação e otimização de uma unidade industrial de extração de óleos essenciais por arraste a vapor

Este trabalho aborda a modelagem matemática do processo de extração de óleos essenciais por arraste a vapor e a otimização deste processo em escala industrial. Utilizouse dados experimentais de óleos essenciais de alecrim e citronela obtidos em uma unidade industrial de extração por arraste a vapor. O modelo é constituído de um sistema de equações diferenciais parciais, que representam o balanço de massa unidimensional nas fases vapor e sólido. Para resolver o modelo desenvolvido, utilizouse o método das linhas com o uso da técnica de diferenças finitas para discretizar a coordenada espacial e suas derivadas, resultando em um sistema de equações diferenciais ordinárias implementado no simulador EMSO. Os parâmetros do modelo foram estimados pelo método dos poliedros flexíveis, utilizando dois conjuntos de dados experimentais e um terceiro para sua validação. O modelo ajustou-se bem aos dados experimentais industriais, resultando, para duas condições ambientais e geográficas das plantas aromáticas de alecrim, em coeficiente global de transferência de massa de 8,64 x 10 -4 s -1 e 9,10 x 10 -4 s -1, e a constante de equilíbrio de 2,72 x 10 -4 m3=kg e 1,07 x 10 -2m3/kg. Enquanto que para a citronela obteve-se coeficiente global de transferência de massa de 8,29 x 10 -4s -1 e 1,23 x 10 -3 s -1, e constante de equilíbrio de 1,25 x 10 -2 m3/kg e 1,52 x 10 -3 m3/kg. Foi proposto também um problema de otimização para determinar o tempo ótimo de extração e o número de bateladas por dia. Os resultados da otimização do processo indicam que se pode reduzir o tempo de extração, permitindo que outras extrações sejam realizadas por dia, maximizando o volume de óleo extraído ao fim do dia.This work has as objective the development of a mathematical model of essential oils extraction by steam distillation and the optimization of this process on an industrial scale. We used experimental data obtained from rosemary and citronella in an industrial steam distillation unit. The model consists a system of partial differential equations, which represent the one-dimensional mass balance in steam and solid phases. To solve the developed model, we used the method of lines with the finite differences technique to discretize the spatial coordinate and its derivatives, resulting in a system of ordinary differential equations implemented in the EMSO simulator. The model parameters were estimated by the method of flexible polyhedra, using two experimental runs and a third experimental run to test the model. The model fitted to the experimental data industry, resulting, for two environmental and geographical conditions of rosemary, in mass transfer coefficients of 8,64 x 10 -4 s -1 and 9,10 x 10 -4 s -1, and the equilibrium constant of 2,72 x 10 -4 m3/kg and 1,07 x 10 -2m3/kg. For citronella, the obtained mass transfer coefficient of 8,29 x 10 -4 s -1 and 1,23 x 10 -3 s -1 and equilibrium constant of 1,25 x 10 -2 m3/kg and 1,52 x 10 -3 m3/kg. It was also proposed an optimization problem to determine the optimum extraction time and number of batches per day. The optimization results indicate that the process can reduce the extraction time, allowing other extractions are performed each day, maximizing the essential oil volume extracted in the day