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.

Modelos de atividade baseados em COSMO se apresentam como uma interessante alternativa para a predição do comportamento de substâncias em mistura. Estes modelos dependem apenas de informações das substâncias puras, como volume e área superficial e outra informação conhecida como perfil sigma, que podem ser determinados através de cálculos de química quântica computacional. O presente trabalho tem por objetivo criar e disponibilizar gratuitamente um banco de dados de perfis sigma para um vasto conjunto de moléculas utilizando o software GAMESS. Para a criação desse banco de dados com informações confiáveis de perfis sigma, foram testados diferentes níveis de teoria para o cálculo da estrutura eletrônica das moléculas, bem como conjuntos de bases distintos, considerando funções de polarização e efeitos difusos. Para realizar a validação dos perfis sigma obtidos, foram estimados valores de coeficientes de atividade em diluição infinita, obtidos através do modelo COSMO-SAC para um conjunto de 689 misturas, as quais não apresentam ligação de hidrogênio. A correlação com os valores de coeficiente de atividade em diluição infinita (IDAC) experimentais, foram construídos gráficos do logaritmo natural de IDAC experimentais e aqueles calculados pelo modelo. O tempo computacional total para obter os perfis sigma também foi avaliado com o intuito de escolher a melhor combinação de teoria e função de base. Dentre os métodos analisados, a teoria de Hartree-Fock (HF) com a função de base 6-311G(d,p) apresentou um bom custo benefício no que se refere ao tempo computacional e a precisão dos resultados frente a dados experimentais, apresentando um coeficiente de determinação superior a 0,94. A partir desses resultados, foi possível criar outros perfis sigma para um número maior de moléculas e disponibilizar essas informações em um banco de dados criado na plataforma do GitHub. / COSMO-based activity coefficient models are widely used to predict non-ideality in mixtures. They require only pure substance properties, such as the -profile, a representation of the screening charge density along the molecular surface. The present work aims to create and freely distribute a -profile database which can be easily extended by the community using the free GAMESS quantum chemistry (QC) package. Different combinations of QC methods and basis sets were tested in the calculation of the -profiles, aiming to define an efficient approach to create the database. The accuracy of COSMO-SAC predictions using these -profiles was evaluated by comparison with experimental infinite dilution activity coefficients and phase equilibrium diagrams. Among the studied alternatives, the Hartree-Fock method with the 6-311G(d,p) basis set led to results with a good balance between accuracy and computational cost. Our database was created using this approach, and it is already available to be used and extended by the community without requiring any commercial software.

O processo de destilação extrativa com N-formilmorfolina (NFM) pode ser utilizado para separação do benzeno de uma corrente contendo hidrocarbonetos não aromáticos com seis carbonos. Cerca de 20 componentes distintos constituem a corrente de interesse. Neste trabalho, primeiramente, buscou-se o conjunto de dados experimentais disponíveis de equilíbrio de fases e coeficiente de atividade em diluição infinita para estimar os parâmetros de interação binária do modelo NRTL capaz de representar o processo de separação de uma planta existente. Os parâmetros faltantes foram estimados a partir do modelo de predição baseado em contribuição de grupos, UNIFAC. O processo de destilação extrativa foi simulado utilizando-se um simulador comercial e uma boa correlação com os dados de planta foi encontrada. A influência na operação de diversos parâmetros, como relação solvente/carga e temperatura do solvente foram estudadas. Os valores ótimos teóricos foram obtidos e podem ser implementados para melhoria do processo. Identificou-se uma extrema sensibilidade do perfil da coluna com o calor fornecido pelo refervedor, indicando que pode ser uma fonte de instabilidades. Além disso, foi identificada baixa eficiência na seção de recheio, que também dificulta a especificação dos produtos de interesse. / Extractive distillation processes with N-formylmorpholine (NFM) are industrially used to separate benzene from non-aromatics hydrocarbons of six carbons. In the process studied in this work, the stream of interest consists of nearly 20 different hydrocarbons. A new set of NRTL parameters was correlated based on literature experimental data. Both vapor-liquid equilibrium data as well as infinite dilution activity coefficient were taken into account; missing parameters were estimated with the UNIFAC group contribution model. The extractive distillation process was simulated using a commercial simulator. Very good agreement with plant data was obtained. The influence of the main parameters in operation, solvent to feed ratio and solvent temperature were studied. The theoretical optimum values were obtained and can be implemented to improve the industrial process. Extreme static sensitivity with respect to reboiler heat was observed indicating that this can be the source of instabilities. In addition a low efficiency was identified at the packing section, contributing to the products specification problem as well.

Atualmente, pelo menos para fins de engenharia, os modelos preditivos de maior sucesso para coeficientes de atividade são os baseados em grupos funcionais, tais como UNIFAC e suas variantes. Enquanto esses modelos requerem grandes quantidades de dados experimentais, os baseados em COSMO-RS (COnductor-like Screening MOdel - for Real Solvents) requerem a calibração de um pequeno conjunto de parâmetros universais. No entanto, a precisão requerida por tarefas de engenharia, tais como a otimização de sistemas de separação, é maior do que a obtida por esta última categoria de modelos. Assim, um novo modelo é proposto neste trabalho, aqui chamado de F-SAC (Functional-Segment Activity Coefficient). Este novo modelo também é baseado no conceito de grupos funcionais, mas a energia de interação entre os grupos vem da teoria COSMO-RS. No presente trabalho, foram consideradas apenas misturas em que não há formação de ligação de hidrogênio ou quando esta pôde ser assumida negligenciável. Assim, foram necessários apenas três parâmetros para descrever cada grupo funcional. A princípio, uma vez ajustados os parâmetros de cada grupo, estes funcionariam para descrever a interação para qualquer par de grupos. Esta é a principal vantagem do modelo proposto. O número de parâmetros do modelo cresce proporcionalmente ao número de grupos funcionais, enquanto que no UNIFAC o número de parâmetros de interação cresce proporcionalmente ao quadrado do número de grupos. Para o banco de dados experimentais de coeficientes de atividade em diluição infinita considerado, a correlação do F-SAC apresentou um erro médio absoluto de 0,07 unidades de ln, enquanto que os modelos UNIFAC (Do) e COSMO-SAC apresentaram, respectivamente, erros de 0,12 e 0,21. O F-SAC foi também avaliado para mais de 1000 misturas binárias de um soluto dissolvido em líquido iônico com dados disponíveis na literatura. O modelo apresentou uma boa correlação aos dados experimentais, com erro médio absoluto de 0,17 unidades de ln, similar ao apresentado pelas misturas orgânicas. O poder de predição do novo modelo foi avaliado utilizando dados de equilíbrio líquido-vapor não considerados no procedimento de ajuste do modelo. Uma ótima concordância com os dados experimentais foi possível em toda a faixa de composição, bem como na predição de azeótropos. Esses resultados demonstram o potencial do modelo proposto. / At present, at least for engineering purposes, the most successful predictive models for activity coefficients are those based on functional groups, such as UNIFAC and its variants. While these models require large amounts of experimental data, the ones based on COSMO-RS require the calibration of a small set of universal parameters. However, the resolution required by engineering tasks, such as the optimization of separation systems, is higher than that obtained by COSMO-RS models. Thus, in this work a novel Functional-Segment Activity Coefficient (F-SAC) model is proposed. This new model is also based on the concept of functional groups, but the interaction energy between groups comes from the COSMO-RS theory. In this study, we considered only mixtures where there is no formation of hydrogen bonds or when they could be assumed negligible. Then, only three parameters were required to describe each functional group. In principle, once the parameters for each functional group are properly calibrated, they should work to describe the interactions with any other group. This is the main advantage of the proposed model, the number of model parameters grows linearly with the number of functional groups, whereas in UNIFAC the number of interaction parameters exhibit quadratic growth with respect to the number of groups. For the experimental database of infinite dilution activity coefficients considered, the correlation of F-SAC have shown a mean absolute error of 0.07 ln-unit. The UNIFAC (Do) and COSMO-SAC models presented errors of 0.12 and 0.21, respectively. F-SAC was also evaluated for more than 1000 binary mixtures of solute in ionic liquid with data available in the literature. Again, the model have shown good correlation to the experimental data, with mean absolute error of 0.17 ln-unit, similar to the performance with the organic mixtures. The predictive strength of the model was assessed by using vaporliquid equilibrium data not considered in the model fitting process. Very good agreement with experimental data was possible over the entire composition range, as well as in the prediction of azeotropes. These results demonstrate the potential of the proposed method.

Neste trabalho, o modelo F-SAC (Functional-Segment Activity Coefficient) foi revisado e estendido de modo a representar, simultaneamente, coeficiente de atividade à diluição infinita (IDAC), equilíbrio líquido-líquido (LLE) e equilíbrio líquido-vapor (VLE), em amplos intervalos de temperatura. Originalmente, o modelo F-SAC foi calibrado essencialmente com dados experimentais de IDAC para diversas misturas. Investigouse, portanto, a possibilidade de o modelo levar também em consideração sistemas que apresentem LLE. Para a presente investigação, dois conjuntos de misturas industrialmente interessantes foram selecionados: hidrocarbonetos (alcanos, cicloalcanos, alcenos, cicloalcenos e compostos aromáticos) + n-formilmorfolina (NFM) e hidrocarbonetos + água. Foram propostas duas modificações na equação que calcula a energia de interação entre os grupos: (i) inspirado pelo comportamento típico observado na expansão isobárica de líquidos, um novo parâmetro foi proposto de modo a diminuir o (pequeno) erro sistemático do modelo com a temperatura no termo sem associação; (ii) em misturas associativas, utilizou-se uma equação similar para representar a dependência da energia de ligação de hidrogênio entre dois segmentos com a temperatura. Na estimação de parâmetros do modelo foram utilizados 1130 pontos experimentais de LLE, para 105 sistemas binários e 16 ternários, juntamente com um 815 dados experimentais de IDAC, resultando numa nova tabela de parâmetros. Ao total, foram estimados 42 parâmetros dos grupos e subgrupos. Os resultados do F-SAC foram comparados, quando possível, com diferentes versões do UNIFAC encontradas na literatura. O modelo proposto foi capaz de correlacionar a maior parte dos dados de IDAC e LLE utilizados. Para os dados de IDAC, foram observados desvios médios absolutos de 0,139 e 0,608 unidades de ln para o F-SAC e o UNIFAC(PSRK), respectivamente. Para os dados de LLE, foram observados desvios semelhantes. Finalmente, o F-SAC foi utilizado para representar dados de VLE para os sistemas de compostos aromáticos + NFM. Observou-se desvios médios relativos de 0,038 e 0,082 no cálculo da pressão de bolha utilizando-se o F-SAC e o UNIFAC(PSRK), respectivamente. Estes resultados comprovam a boa capacidade preditiva do modelo, uma vez que nenhum dado de VLE foi utilizado nas estimações. Conclui-se, desta forma, que o modelo FSAC pode ser usado para representar, simultaneamente, IDAC, LLE e VLE de misturas com e sem ligação de hidrogênio, utilizando um único conjunto de parâmetros. / In this work, the F-SAC model (Functional -Segment Activity Coefficient) has been revised and extended in order to simultaneously represent activity coefficient at infinite dilution (IDAC), liquid-liquid equilibrium (LLE) and vapor-liquid equilibrium (VLE), in a wide range of temperature. Originally, the F- SAC model was calibrated essentially with IDAC experimental data for several mixtures. In the present work, the possibility of also taking into account LLE was investigated. For this study, two sets of industrially interesting mixtures were selected : hydrocarbons (alkanes, cycloalkanes, alkenes, cycloalkenes and aromatics) + n-formylmorpholine (NFM) and hydrocarbons + water. Modifications were suggested in the equation that computes the interaction energy for each contact between two segments: (i) inspired by the typical behavior observed in the isobaric expansion of liquids, a new parameter has been proposed in order to reduce the (small) systematic error with the temperature for non-associating systems; (ii) for associating mixtures, a similar equation was proposed to represent the temperature dependence of the hydrogen bond energy between two segments. In the parameterization of the model, 1130 LLE experimental points were used, for 105 binary and 16 ternary systems, along with 815 IDAC experimental points, resulting in a new parameter table. A total of 42 parameters for groups and subgroups were calibrated. The F-SAC results were compared, whenever possible, with different versions of the UNIFAC model, freely available in the literature. The proposed model was able to correlate most of the IDAC and LLE data. For IDAC data, absolute average deviations were 0.139 and 0.608 ln units using F-SAC and UNIFAC(PSRK), respectively. Similar deviations were observed for LLE data. Finally, F-SAC was used to represent VLE data for aromatics + NFM systems. The absolute relative deviations were 0.038 and 0.082 in the bubble pressure calculations for F-SAC and UNIFAC(PSRK), respectively. The VLE results indicate a strong predictive ability of the suggested model, since no VLE data was used in the parameter estimation procedure. The results also indicate that the F-SAC model can be used to represent, simultaneously, IDAC, LLE and VLE for associating and non-associating mixtures, with a single set of parameters.

O modelo F–SAC (Functional–Segment Activity Coefficient), recentemente introduzido, combina a ideia de contribuição de grupos com a teoria de superfícies de contato COSMO–RS. Este modelo tem uma dependência reduzida dos parâmetros de interação binária quando comparado com as variantes do modelo UNIFAC e precisão melhorada quando comparada com modelos baseados em COSMO. No presente trabalho uma modificação na formulação do modelo F–SAC é proposta para a inclusão de interações dispersivas. Para testar esta modificação, foram considerados os dados experimentais de alcanos, ciclo–alcanos, alcenos, aromáticos e perfluorocarbonetos. O modelo proposto foi capaz de correlacionar entalpias de vaporização de substâncias puras, não consideradas em versões anteriores do modelo. Uma vez que a capacidade calorífica de líquidos está intimamente relacionada com a entalpia de vaporização, o modelo também pode prever a capacidade calorífica de substâncias puras. Em relação ao coeficiente de atividade em diluição infinita e dados de equilíbrio líquido–vapor, o modelo modificado manteve o bom desempenho do modelo original, também semelhante a outros modelos similares, como variantes do modelo UNIFAC. Além disso, o modelo modificado pode agora calcular valores consistentes para a entalpia e entropia de excesso para sistemas onde as interações são principalmente dispersivas. Para estes sistemas, a maioria dos modelos de coeficiente de atividade prevê entalpia de excesso zero, contrabalançando valores de energia de Gibbs de excesso confiáveis com entropias de excesso distorcidas. / The recently introduced F–SAC (Functional–Segment Activity Coefficient) model combines the group contribution idea with a COSMO–RS surface contacting theory. This model has a reduced dependency on binary interaction parameters when compared to classical UNIFAC type models and improved accuracy when compared with COSMO based models. In the present work a modified F–SAC formulation is proposed for including dispersive interactions. For testing the modification, experimental data of alkanes, cycloalkanes, alkenes aromatics and perfluorocarbons were considered. The proposed model was able to correlate pure compound enthalpies of vaporization, not considered in previous versions of the model. Since the heat capacity of liquids is closely related to the enthalpies of vaporization, the model also can predict pure compound heat capacity of liquids. Regarding mixture infinite dilution activity coefficient and vapor–liquid equilibrium data, the modified model maintained the good performance of the original model, also similar to other competing models such as UNIFAC variants. Additionally, the modified model now can compute consistent values for the excess enthalpy and entropy for systems where the interactions are mainly dispersive. For these systems, most activity coefficient models predict zero excess enthalpies, counterbalancing that with distorted excess entropies.

Neste trabalho apresenta-se a modelagem de chamas difusivas na forma de jato, para baixo número de Mach e elevado número de Damköhler. O modelo é baseado na solução das equações na forma flamelet para a parte química e na fração de mistura para o fluxo. Este modelo descreve bem o comportamento de chamas difusivas, exceto na sua extremidade (ponta), onde geralmente surgem instabilidades. Resultados numéricos são apresentados para uma cinética química de uma e multietapas, utilizando a técnica LES (Large-Eddy Simulation) com o modelo de Smagorinsky para a viscosidade turbulenta. A discretização das equações governantes é feita em diferenças finitas, com a aplicação da técnica TVD (Total Variation Diminishing). Além disso, apresentamse mecanismos reduzidos multietapas para o metano, o metanol e o etanol, visando obter resultados realistas. A modelagem de chamas de metanol e etanol diferencia-se da modelagem de chamas de metano por ocorrer uma mudança de fase antes da combustão. Modela-se o efeito global das gotas usando uma descrição Lagrangeana que é incorporada à descrição Euleriana do escoamento, via termos fonte. Testes numéricos foram realizados para chamas difusivas de metano, metanol e etanol, e os resultados estão em concordância com os dados encontrados na literatura. / This work presents a model for a jet diffusion flame, for low Mach and high Damköhler numbers. The model is based on the solution of the flamelet equations for the chemistry and on the mixture fraction for the flow. This model describes well the behavior of diffusion flames, except at the flame tip, where instabilities can often occur. Numerical results are presented for an one-step and multi-step chemical kinetic models, using the LES (Large-Eddy Simulation) technique with the Smagorinsky model for the turbulent viscosity. The discretization of the governing equations follows the finite difference method, with the application of the TVD (Total Variation Diminishing) technique. Besides, multi-step reduced mechanisms for the methane, the methanol and the ethanol are employed, obtaining realistic results. The flame modeling of methanol and ethanol differs from the modeling of methane flames because of a phase change occurs before the combustion. The droplets global effect is modeled based on a Lagrangian description, which is incorporated into the Eulerian description of the flow through source terms. Numerical tests were carried out for methane, methanol and ethanol diffusion flames, and the results compare well with data in the literature.

A modelagem termodinâmica do equilíbrio de fases e a determinação de propriedades termodinâmicas de sistemas em que as moléculas apresentam associação é um problema desafiador na indústria química. Embora o uso de equações de estado tenha sido durante muito tempo restrito aos sistemas de fluidos simples, há uma crescente demanda por modelos que também sejam adequados a compostos complexos e macromoleculares. O avanço teórico mais aparente para equações de estado com tais capacidades foi conseguido por aplicação dos princípios da mecânica estatística, o qual considera as forças de associação entre as moléculas em sua teoria molecular. Um exemplo desta categoria é a equação de estado PC-SAFT (Perturbed Chain-Statistical Associating Fluid Theory), cuja formulação requer apenas um pequeno conjunto de parâmetros para cada componente e estas interações são computadas através de regras de mistura. Desta forma, os objetivos deste trabalho foram a implementação do termo de associação desta equação e o desenvolvimento de ferramentas para o cálculo de propriedades termodinâmicas de sistemas contendo fluidos associados e não-associados, bem como para o cálculo do equilíbrio de fases. Todas as funções foram desenvolvidas no software MATLAB e a verificação da implementação foi feita através de comparação entre os resultados obtidos com estas ferramentas e dados experimentais e/ou outras equações de estado presentes na literatura. Através da comparação com dados de equilíbrio de fases disponíveis na literatura, foi analisado o desempenho da equação PC-SAFT com associação para diferentes sistemas multicomponentes e múltiplas fases envolvendo substâncias polares e apolares. Os resultados obtidos foram satisfatórios quando comparados a dados experimentais e aqueles obtidos através de outras equações de estado. / Modeling the phase equilibrium and thermodynamic properties of systems in which molecules exhibit associating interactions remains a challenging problem in chemical industry. Although the use of equations of state has for a long time been restricted to systems of simple fluids, there is an increasing demand for models that are also suitable for complex and macromolecular compounds. Progress toward this goal has been achieved by applying principles of statistical mechanics, which takes molecular theories into account. One example of such approach is the PC-SAFT (Perturbed Chain-Statistical Associating Fluid Theory) equation of state, which requires only a small set of parameters for each component and whose interactions are obtained by mixing rules. The goals of this work have been the implementation of the PC-SAFT association term and the development of tools for thermodynamic properties calculation of associating and non-associating systems, as well as for phase equilibria calculation. All the functions have been developed using MATLAB software and their implementation analysis has been made by comparing obtained results to experimental data and those obtained by other equations of state available in literature. Through these comparisons, PC-SAFT performance has been analyzed for different multicomponent and multiphase systems involving polar and non-polar substances. The obtained results have shown good correlation to both experimental data and those obtained by other equations of state.

A modelagem termodinâmica do equilíbrio de fases e a determinação de propriedades termodinâmicas de sistemas em que as moléculas apresentam associação é um problema desafiador na indústria química. Embora o uso de equações de estado tenha sido durante muito tempo restrito aos sistemas de fluidos simples, há uma crescente demanda por modelos que também sejam adequados a compostos complexos e macromoleculares. O avanço teórico mais aparente para equações de estado com tais capacidades foi conseguido por aplicação dos princípios da mecânica estatística, o qual considera as forças de associação entre as moléculas em sua teoria molecular. Um exemplo desta categoria é a equação de estado PC-SAFT (Perturbed Chain-Statistical Associating Fluid Theory), cuja formulação requer apenas um pequeno conjunto de parâmetros para cada componente e estas interações são computadas através de regras de mistura. Desta forma, os objetivos deste trabalho foram a implementação do termo de associação desta equação e o desenvolvimento de ferramentas para o cálculo de propriedades termodinâmicas de sistemas contendo fluidos associados e não-associados, bem como para o cálculo do equilíbrio de fases. Todas as funções foram desenvolvidas no software MATLAB e a verificação da implementação foi feita através de comparação entre os resultados obtidos com estas ferramentas e dados experimentais e/ou outras equações de estado presentes na literatura. Através da comparação com dados de equilíbrio de fases disponíveis na literatura, foi analisado o desempenho da equação PC-SAFT com associação para diferentes sistemas multicomponentes e múltiplas fases envolvendo substâncias polares e apolares. Os resultados obtidos foram satisfatórios quando comparados a dados experimentais e aqueles obtidos através de outras equações de estado. / Modeling the phase equilibrium and thermodynamic properties of systems in which molecules exhibit associating interactions remains a challenging problem in chemical industry. Although the use of equations of state has for a long time been restricted to systems of simple fluids, there is an increasing demand for models that are also suitable for complex and macromolecular compounds. Progress toward this goal has been achieved by applying principles of statistical mechanics, which takes molecular theories into account. One example of such approach is the PC-SAFT (Perturbed Chain-Statistical Associating Fluid Theory) equation of state, which requires only a small set of parameters for each component and whose interactions are obtained by mixing rules. The goals of this work have been the implementation of the PC-SAFT association term and the development of tools for thermodynamic properties calculation of associating and non-associating systems, as well as for phase equilibria calculation. All the functions have been developed using MATLAB software and their implementation analysis has been made by comparing obtained results to experimental data and those obtained by other equations of state available in literature. Through these comparisons, PC-SAFT performance has been analyzed for different multicomponent and multiphase systems involving polar and non-polar substances. The obtained results have shown good correlation to both experimental data and those obtained by other equations of state.

A pesquisa acerca da modelagem de colunas de destilação não é mais uma novidade nos dias de hoje. Os primeiros trabalhos, com metodologias propostas para a solução de sistemas de separação modelados prato a prato, surgiram na década de 30. Porém, só a partir da década de 50, com o advento do computador digital, foram realizados investimentos sólidos no desenvolvimento de novos algoritmos e simuladores. Apesar deste investimento, apenas modelos simplificados eram utilizados devido à baixa capacidade de processamento. A partir da década de 70, os primeiros simuladores comerciais começaram a ser introduzidos na indústria e o desenvolvimento de modelos rigorosos não parou mais. Dentro deste contexto, este trabalho tem como objetivo o desenvolvimento de um modelo rigoroso de coluna de destilação, genérico, dinâmico e de fácil aplicação nos mais variados tipos de estudos, desde simples simulações de operação até otimizações e previsões de comportamentos de parada e partida de plantas. Ou seja, um modelo projetado para situações onde uma representação fiel do comportamento dinâmico é necessária. Os modelos gerados neste estudo foram implementados na linguagem de modelagem do simulador EMSO e fazem parte da biblioteca EML (EMSO Model Library). Esta biblioteca é distribuída no conceito de software livre, disponibilizando todos os modelos via internet e sem custo. Para a validação dos modelos desenvolvidos foram utilizados dois exemplos: uma coluna de destilação reativa com dados provenientes da literatura e uma coluna deisobutanizadora industrial da PETROBRAS com 80 pratos que separa o isobutano de uma mistura de 13 componentes. Este último problema resulta em um sistema com mais de 6300 variáveis, provando que o simulador utilizado é muito eficiente para a resolução de problemas de grande escala. Os resultados da validação foram satisfatórios, assim como o desempenho do simulador EMSO quanto à eficiência computacional e à robustez na solução do problema. / Nowadays, modeling and simulation of distillation columns are not a new research subject anymore. It was around the 1930’s that the first methodologies were proposed to solve staged separation systems. However, it was in 1950’s, with the advent of digital computers that most important investments were made for the development of new algorithms and simulators. In spite of these investments, only very simple models were used due to the low speed of the computers available at that time. In the 1970’s, the industry started to use the first commercial simulators and the development of rigorous models for distillation columns has not stopped since then. Inside this context, the main goal of this work is the development of a rigorous, generic, dynamic and easy-to-use model of distillation column to be applied on diferent targets, from operation simulations to dynamic optimizations and the prediction of the start up and shut down plant behaviour. To do so, the model must represent the dynamic performance of the unit. The models generated were implemented in the EMSO simulator environment with its own modeling language and are part of the EMSO Model Library (EML), which is distributed under the Free Software concept. In order to analyze the model performance, two examples were chosen: a reactive distillation column from the literature and an 80 tray deisobutanizer column from PETROBRAS which separates isobutane of a 13 component mixture. The solution of the last problem generates a 6300 variable system, that is used to confirm the great efficiency of EMSO to solve large scale problems. The validation results have been as satisfactory as the model prediction and the EMSO simulator performance in computational efficiency and robustness.