Este trabalho tem por objetivo o desenvolvimento de procedimentos hierárquicos de cálculo para os parâmetros de avaliação de trocadores de calor tipo casco e tubos com precisão adequada e baixo custo computacional, facilitando a simulação e otimização de processos, implementados dentro do simulador EMSO. Com os recursos da linguagem de programação orientada a objetos disponível no simulador EMSO, foi possível estruturar a modelagem de equipamentos de troca térmica aproveitando ao máximo os conceitos de composição e herança e facilitando o desenvolvimento de novos modelos de trocadores de calor. Os modelos criados foram incorporados à biblioteca de modelos do simulador, atendendo os trocadores do tipo casco simples, casco duplo e trocadores casco e tubos multipasses conforme a designação das normas TEMA para esses equipamentos. O método de avaliação para o lado do casco é o método de Bell- Delaware amplamente descrito na literatura sobre o assunto. Para o lado dos tubos a transferência de calor e a perda de carga foram estimadas usando correlações existentes. Os procedimentos de cálculo também incluem a possibilidade de discretização do trocador para uma maior acurácia na avaliação dos coeficientes de troca térmica e na variação das propriedades físicas dos fluidos dentro do trocador. Para a validação dos modelos, várias simulações foram realizadas em uma bateria de trocadores de calor de uma unidade de destilação atmosférica de uma refinaria de petróleo. Os resultados encontrados foram satisfatórios podendo o usuário ter fácil acesso aos principais parâmetros que influenciam no desempenho do trocador de calor, ou incorporar novos procedimentos para o cálculo destes equipamentos. / The aim of this work is the development of hierarchical procedures for computing the performance parameters of shell and tube heat exchangers, with accuracy and low computational cost, facilitating the process simulation and optimization, implemented in the EMSO simulator. With the object oriented programming language resources available in the EMSO simulator, it was possible to structure the thermal exchange equipment modeling in a modular manner, making use of the composition and inheritance concepts. The models created have been incorporated into the library of simulator models, comprehending heat exchangers of type E Shell, F Shell, and multipass, in agreement with the TEMA standards for such equipment. The analysis method used for the shell side was the Bell-Delaware method, well described in the open literature. For the tube side, the heat transfer and pressure drop were estimated using the available correlations. The calculation procedures also include the discretization possibility of the exchanger for a better accuracy in the heat transfer coefficient evaluation for thermal exchange and the physical property variations of fluids inside the exchanger. For model validation, several simulations have been carried out in a series of heat exchangers of an atmospheric distillation unit of an oil refinery. The obtained results were satisfactory and the models also provide to the user an easy access to the main parameters of the heat exchanger and a flexible way to incorporate new procedures for the calculation of these pieces of equipment.

Processos de polimerização de polietileno linear de baixa densidade em fase gás normalmente removem monômeros e outros hidrocarbonetos dissolvidos do polímero formado através de um sistema de stripping com nitrogênio ou vapor. A remoção é necessária para reduzir as perdas com monômero, assegurar baixos teores de componentes voláteis no polímero final exigidos pela legislação ambiental e reduzir riscos de explosão em silos e embalagens. Para construir uma metodologia de modelagem destes equipamentos especiais, os fenômenos de solubilidade e equilíbrio polímero-hidrocarboneto, bem como da difusividade multicomponente para polímeros são revistos e acoplados em um sistema de equações de transferência de massa e energia para um sistema de leito descendente com nitrogênio ou vapor em contracorrente. As equações diferenciais parciais para o modelo de leito descendente construído com as considerações acima são resolvidas através de diferenças finitas para as coordenadas axiais e por colocação ortogonal pelo método de elementos finitos para a coordenada radial utilizando o software gPROMS. Os resultados do modelamento são comparados com dados de equipamentos industriais utilizando nitrogênio e vapor. Para ajustar as condições industriais com o modelo proposto, o comprimento de difusão dentro da partícula teve que ser ajustado para considerar a influência da morfologia das partículas e da má distribuição de gás na transferência de massa. Logo, embora um entendimento do fenômeno de stripping de hidrocarbonetos pôde ser alcançado, uma melhor representação das ineficiências na transferência de massa é fundamental para melhorar a predição do modelo. Foi também observado que as condições de processo do reator de polimerização afetam a remoção de hidrocarbonetos, exigindo novas condições ótimas de stripping para cada produto. A metodologia formulada para equipamentos de stripping de polietileno pode ser utilizada para avaliação do impacto de novos polímeros no processo ou como linha base de dimensionamento de novos equipamentos já que há pouca literatura disponível neste campo. / Gas-phase linear low density polyethylene polymerization processes normally remove monomers and other dissolved hydrocarbons from the reacted polymer in steam or nitrogen stripping systems. The stripping is performed in order to decrease the monomer losses as well as to assure low volatile components in the final polymer to mitigate explosion risks and meet environmental standards. In order to increase the knowledge of these special equipments, the separation phenomena is reviewed and a mathematical model is developed for the counter-current steam or nitrogen stripping process to remove hydrocarbons from polyethylene under a moving-bed vessel. The modeling methodology takes into account the current polymer-solvent solubility and diffusion theories, the interphase equilibrium, the mass and heat transfer for a multicomponent diffusion-limited spherical polymer particle and the plug-flow and funnel-flow solid behavior. The partial differential equation (PDE) for the moving bed model built from the above considerations is solved by the method of lines using the finite differences method for the axial coordinates and orthogonal collocation in finite elements method for the radial coordinate using the gPROMS software. The modeling results are compared with industrial polyethylene stripping equipment data using nitrogen and steam. To fit the industrial conditions to the model, the diffusion length had to be adjusted to consider the particles morphology and gas distribution influence in the mass transfer. So, although a better stripping phenomena understanding could be reached, a correct representation of the mass transfer inefficiencies is crucial to increase the model prediction capability. It was also noted that the reactor process conditions affect the equipment performance asking new optimum stripping conditions for each reactor product. The formulated methodology for polyethylene stripping equipments can be used in the new process developments or as base line for new equipment sizing since a poor literature in this field is available.

A necessidade cada vez maior do uso de produtos poliméricos sintéticos, como para produção de embalagens, partes de eletrodomésticos, eletroeletrônicos, automóveis, etc., tem levado a indústria de polímeros a buscar cada vez mais a diminuição do desperdício e aumento da qualidade dos produtos. Para isso tem-se buscado entender melhor como as propriedades reológicas dos polímeros afetam seu processamento e a qualidade final dos produtos. Com o intuito de se obter resultados mais rápidos e com menor custo recorre-se cada vez mais a estudos de modelagem e simulação de processos de transformação de polímeros. Neste trabalho é apresentada uma nova ferramenta de CFD para simulações de escoamentos envolvendo fluidos viscoelásticos, o viscoelasticFluidFoam solver. A implementação do módulo foi feita no pacote de CFD OpenFOAM devido principalmente às vantagens oferecidas por esse software, como por exemplo, possibilidade de uso de geometrias complexas, malhas não-estruturadas, técnicas multigrid e paralelização do processamento de dados, além de ser um software gratuito e de código aberto. Foi feita a implementação do modelo de Maxwell, UCM, Oldroyd-B, Giesekus, FENE-P, FENE-CR, PTT na forma linear e exponencial, e DCPP, todos na forma multimodo. Dentre as várias metodologias disponíveis para resolver o problema da obtenção de soluções estáveis a altos valores deWeissenberg foi escolhida a DEVSS devido a sua estabilidade e aplicação a modelos complexos. Para se fazer a validação do solver desenvolvido foi feito a comparação com resultados numéricos e experimentais obtidos da literatura. É mostrada uma comparação entre vários modelos para obtenção da velocidade e diferença de tensões normais para um escoamento em uma contração plana abrupta 4:1. Os resultados obtidos foram satisfatórios sendo possível dar credibilidade ao solver implementado e garantir a disponibilidade de uma boa ferramenta para estudo de fluidos viscoelásticos para ser usada tanto no meio acadêmico como no setor industrial. / Synthetic polymer products are of great importance in several industrial sectors, such as for production of packaging, parts of appliances, electronics, and cars. Due to the increasing demand for this kind of material, reduction of waste and increase of quality has become a key issue in polymer industry. In this sense modeling and simulation of processing operations appears as a fundamental tool, leading to better understanding of how the rheological properties of polymers affect their processability and final product quality, and reducing time and costs related to the development of processes and products. This work presents a new Computational Fluid Dynamics (CFD) tool for the simulation of viscoelastic fluid flows, called viscoelasticFluidFoam solver, which consists of a viscoelastic fluid module to be used OpenFOAM CFD package. The advantages of using OpenFOAM as development platform include its characteristics with relation to flexibility to deal with complex geometries, unstructured and non orthogonal meshes, moving meshes, large variety of interpolation schemes and solvers for the linear discretized system, and the possibility of data processing parallelization. Linear Maxwell, Oldroyd-B, Giesekus, Phan-Thien-Tanner (PTT), the Finitely Extensible Nonlinear Elastic (FENE-P and FENE-CR), and DCPP constitutive equations have been implemented, in single and the multimode form. Among the various available methodologies to solve the problem of obtaining stable solutions to high Weissenberg values, the DEVSS was chosen due to its stability and application to complex models. The viscoelasticFluidFoam solver was tested by comparing its predictions with experimental and numerical data from literature for the analysis of a planar 4:1 contraction flow. These tests have shown the great potential of this solver for application both in academia and in industry.

A modelagem e simulação de reatores trifásicos representa um grande desafio para simuladores de processo. Muitos estudos acadêmicos têm sido desenvolvidos com o objetivo de descrever e prever o comportamento desses equipamentos, quer seja em seu estado estacionário, quer seja as respostas dinâmicas. Várias são as abordagens aplicadas para tentar melhor refletir os fenômenos envolvidos nos reatores trifásicos, ponderando sempre o custo computacional de cada abordagem. Neste trabalho, um reator de hidrogenação seletiva de nafta de pirólise foi simulado baseado em dados publicados na literatura. O reator é do tipo leito gotejante (trifásico), e as reações consideradas são de pseudo-primeira ordem. A técnica aplicada foi a de modelagem matemática por células, onde os leitos catalíticos foram subdivididos em reatores tipo CSTR dinâmicos associados em série. A cada célula, um cálculo de flash foi associado, aperfeiçoando os balanços de massa e energia comumente empregados em reatores de leito gotejante. A abordagem termodinâmica utilizada para prever o equilíbrio líquido-vapor foi a i - i, com a equação de estado SRK associada a parâmetros de interação binária específicos para a solubilidade de hidrogênio em gasolina de pirólise. Os parâmetros de interação binária são provenientes do pacote termodinâmico do simulador iiSE. O modelo do reator foi implementado no software EMSO (Environment for Modeling, Simulation and Optimization) consistindo em cerca de 9000 equações e variáveis. Os resultados obtidos com o modelo contruído foram similares aos reportados na literatura. A aplicação de modelagem por células mostrou-se não só aplicável mas também mais robusta do que as abordagens tradicionais que utilizam equações diferenciais ordinárias. A utilização da ferramenta EMSO para a modelagem por células mostrou-se ainda mais vantajosa ao permitir a avaliação do comportamento dinâmico do reator em algumas situações hipotéticas, mas que são bem comuns na indústria. / Modeling and simulation of three-phase reactors is a challenge for process simulators. Many academic studies have been developed in order to describe and predict the behavior of these equipments for both steady states or dynamic responses. There are several approaches applied trying to reflect better the phenomena involved in three-phase reactors, always considering the computational cost of each approach. In this work a selective hydrogenation reactor for pyrolysis naphtha was simulated based on published literature data.The reactor is trickle bed type and the reactions are considered pseudofirst order. The technique used is cell-network modeling, where the catalyst beds were subdivided into dynamics CSTR reactors linked in series. To each cell there is a calculation Flash associated, enhancing mass and energy balances commonly applied in TBR. The thermodynamic approach used to predict the vapor-liquid equilibrium was i- i with SRK equation of state, associated to specific binary interaction parameters for the solubility of hydrogen in pygas from the thermodynamic simulator package iiSE. The binary interaction parameters come from the thermodynamic package of IISE simulator. The reactor model was implemented in software EMSO (Environment for Modeling, Simulation and Optimization) and consists of 9,000 equations and variables approximately. The results obtained with the model constructed here were similar to those reported in the literature. The use of cell-network modeling proved to be not only applicable but also more robust than the traditional approaches that use ordinary differential equations. The use of EMSO tool for cell-network modeling proved to be even more advantageous because allows the evaluation of some dynamic behavior of the reactor for hypothetical situations but quite common in the industry.

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.

Este trabalho tem por objetivo o desenvolvimento de procedimentos hierárquicos de cálculo para os parâmetros de avaliação de trocadores de calor tipo casco e tubos com precisão adequada e baixo custo computacional, facilitando a simulação e otimização de processos, implementados dentro do simulador EMSO. Com os recursos da linguagem de programação orientada a objetos disponível no simulador EMSO, foi possível estruturar a modelagem de equipamentos de troca térmica aproveitando ao máximo os conceitos de composição e herança e facilitando o desenvolvimento de novos modelos de trocadores de calor. Os modelos criados foram incorporados à biblioteca de modelos do simulador, atendendo os trocadores do tipo casco simples, casco duplo e trocadores casco e tubos multipasses conforme a designação das normas TEMA para esses equipamentos. O método de avaliação para o lado do casco é o método de Bell- Delaware amplamente descrito na literatura sobre o assunto. Para o lado dos tubos a transferência de calor e a perda de carga foram estimadas usando correlações existentes. Os procedimentos de cálculo também incluem a possibilidade de discretização do trocador para uma maior acurácia na avaliação dos coeficientes de troca térmica e na variação das propriedades físicas dos fluidos dentro do trocador. Para a validação dos modelos, várias simulações foram realizadas em uma bateria de trocadores de calor de uma unidade de destilação atmosférica de uma refinaria de petróleo. Os resultados encontrados foram satisfatórios podendo o usuário ter fácil acesso aos principais parâmetros que influenciam no desempenho do trocador de calor, ou incorporar novos procedimentos para o cálculo destes equipamentos. / The aim of this work is the development of hierarchical procedures for computing the performance parameters of shell and tube heat exchangers, with accuracy and low computational cost, facilitating the process simulation and optimization, implemented in the EMSO simulator. With the object oriented programming language resources available in the EMSO simulator, it was possible to structure the thermal exchange equipment modeling in a modular manner, making use of the composition and inheritance concepts. The models created have been incorporated into the library of simulator models, comprehending heat exchangers of type E Shell, F Shell, and multipass, in agreement with the TEMA standards for such equipment. The analysis method used for the shell side was the Bell-Delaware method, well described in the open literature. For the tube side, the heat transfer and pressure drop were estimated using the available correlations. The calculation procedures also include the discretization possibility of the exchanger for a better accuracy in the heat transfer coefficient evaluation for thermal exchange and the physical property variations of fluids inside the exchanger. For model validation, several simulations have been carried out in a series of heat exchangers of an atmospheric distillation unit of an oil refinery. The obtained results were satisfactory and the models also provide to the user an easy access to the main parameters of the heat exchanger and a flexible way to incorporate new procedures for the calculation of these pieces of equipment.

Processos de polimerização de polietileno linear de baixa densidade em fase gás normalmente removem monômeros e outros hidrocarbonetos dissolvidos do polímero formado através de um sistema de stripping com nitrogênio ou vapor. A remoção é necessária para reduzir as perdas com monômero, assegurar baixos teores de componentes voláteis no polímero final exigidos pela legislação ambiental e reduzir riscos de explosão em silos e embalagens. Para construir uma metodologia de modelagem destes equipamentos especiais, os fenômenos de solubilidade e equilíbrio polímero-hidrocarboneto, bem como da difusividade multicomponente para polímeros são revistos e acoplados em um sistema de equações de transferência de massa e energia para um sistema de leito descendente com nitrogênio ou vapor em contracorrente. As equações diferenciais parciais para o modelo de leito descendente construído com as considerações acima são resolvidas através de diferenças finitas para as coordenadas axiais e por colocação ortogonal pelo método de elementos finitos para a coordenada radial utilizando o software gPROMS. Os resultados do modelamento são comparados com dados de equipamentos industriais utilizando nitrogênio e vapor. Para ajustar as condições industriais com o modelo proposto, o comprimento de difusão dentro da partícula teve que ser ajustado para considerar a influência da morfologia das partículas e da má distribuição de gás na transferência de massa. Logo, embora um entendimento do fenômeno de stripping de hidrocarbonetos pôde ser alcançado, uma melhor representação das ineficiências na transferência de massa é fundamental para melhorar a predição do modelo. Foi também observado que as condições de processo do reator de polimerização afetam a remoção de hidrocarbonetos, exigindo novas condições ótimas de stripping para cada produto. A metodologia formulada para equipamentos de stripping de polietileno pode ser utilizada para avaliação do impacto de novos polímeros no processo ou como linha base de dimensionamento de novos equipamentos já que há pouca literatura disponível neste campo. / Gas-phase linear low density polyethylene polymerization processes normally remove monomers and other dissolved hydrocarbons from the reacted polymer in steam or nitrogen stripping systems. The stripping is performed in order to decrease the monomer losses as well as to assure low volatile components in the final polymer to mitigate explosion risks and meet environmental standards. In order to increase the knowledge of these special equipments, the separation phenomena is reviewed and a mathematical model is developed for the counter-current steam or nitrogen stripping process to remove hydrocarbons from polyethylene under a moving-bed vessel. The modeling methodology takes into account the current polymer-solvent solubility and diffusion theories, the interphase equilibrium, the mass and heat transfer for a multicomponent diffusion-limited spherical polymer particle and the plug-flow and funnel-flow solid behavior. The partial differential equation (PDE) for the moving bed model built from the above considerations is solved by the method of lines using the finite differences method for the axial coordinates and orthogonal collocation in finite elements method for the radial coordinate using the gPROMS software. The modeling results are compared with industrial polyethylene stripping equipment data using nitrogen and steam. To fit the industrial conditions to the model, the diffusion length had to be adjusted to consider the particles morphology and gas distribution influence in the mass transfer. So, although a better stripping phenomena understanding could be reached, a correct representation of the mass transfer inefficiencies is crucial to increase the model prediction capability. It was also noted that the reactor process conditions affect the equipment performance asking new optimum stripping conditions for each reactor product. The formulated methodology for polyethylene stripping equipments can be used in the new process developments or as base line for new equipment sizing since a poor literature in this field is available.

A necessidade cada vez maior do uso de produtos poliméricos sintéticos, como para produção de embalagens, partes de eletrodomésticos, eletroeletrônicos, automóveis, etc., tem levado a indústria de polímeros a buscar cada vez mais a diminuição do desperdício e aumento da qualidade dos produtos. Para isso tem-se buscado entender melhor como as propriedades reológicas dos polímeros afetam seu processamento e a qualidade final dos produtos. Com o intuito de se obter resultados mais rápidos e com menor custo recorre-se cada vez mais a estudos de modelagem e simulação de processos de transformação de polímeros. Neste trabalho é apresentada uma nova ferramenta de CFD para simulações de escoamentos envolvendo fluidos viscoelásticos, o viscoelasticFluidFoam solver. A implementação do módulo foi feita no pacote de CFD OpenFOAM devido principalmente às vantagens oferecidas por esse software, como por exemplo, possibilidade de uso de geometrias complexas, malhas não-estruturadas, técnicas multigrid e paralelização do processamento de dados, além de ser um software gratuito e de código aberto. Foi feita a implementação do modelo de Maxwell, UCM, Oldroyd-B, Giesekus, FENE-P, FENE-CR, PTT na forma linear e exponencial, e DCPP, todos na forma multimodo. Dentre as várias metodologias disponíveis para resolver o problema da obtenção de soluções estáveis a altos valores deWeissenberg foi escolhida a DEVSS devido a sua estabilidade e aplicação a modelos complexos. Para se fazer a validação do solver desenvolvido foi feito a comparação com resultados numéricos e experimentais obtidos da literatura. É mostrada uma comparação entre vários modelos para obtenção da velocidade e diferença de tensões normais para um escoamento em uma contração plana abrupta 4:1. Os resultados obtidos foram satisfatórios sendo possível dar credibilidade ao solver implementado e garantir a disponibilidade de uma boa ferramenta para estudo de fluidos viscoelásticos para ser usada tanto no meio acadêmico como no setor industrial. / Synthetic polymer products are of great importance in several industrial sectors, such as for production of packaging, parts of appliances, electronics, and cars. Due to the increasing demand for this kind of material, reduction of waste and increase of quality has become a key issue in polymer industry. In this sense modeling and simulation of processing operations appears as a fundamental tool, leading to better understanding of how the rheological properties of polymers affect their processability and final product quality, and reducing time and costs related to the development of processes and products. This work presents a new Computational Fluid Dynamics (CFD) tool for the simulation of viscoelastic fluid flows, called viscoelasticFluidFoam solver, which consists of a viscoelastic fluid module to be used OpenFOAM CFD package. The advantages of using OpenFOAM as development platform include its characteristics with relation to flexibility to deal with complex geometries, unstructured and non orthogonal meshes, moving meshes, large variety of interpolation schemes and solvers for the linear discretized system, and the possibility of data processing parallelization. Linear Maxwell, Oldroyd-B, Giesekus, Phan-Thien-Tanner (PTT), the Finitely Extensible Nonlinear Elastic (FENE-P and FENE-CR), and DCPP constitutive equations have been implemented, in single and the multimode form. Among the various available methodologies to solve the problem of obtaining stable solutions to high Weissenberg values, the DEVSS was chosen due to its stability and application to complex models. The viscoelasticFluidFoam solver was tested by comparing its predictions with experimental and numerical data from literature for the analysis of a planar 4:1 contraction flow. These tests have shown the great potential of this solver for application both in academia and in industry.

A modelagem e simulação de reatores trifásicos representa um grande desafio para simuladores de processo. Muitos estudos acadêmicos têm sido desenvolvidos com o objetivo de descrever e prever o comportamento desses equipamentos, quer seja em seu estado estacionário, quer seja as respostas dinâmicas. Várias são as abordagens aplicadas para tentar melhor refletir os fenômenos envolvidos nos reatores trifásicos, ponderando sempre o custo computacional de cada abordagem. Neste trabalho, um reator de hidrogenação seletiva de nafta de pirólise foi simulado baseado em dados publicados na literatura. O reator é do tipo leito gotejante (trifásico), e as reações consideradas são de pseudo-primeira ordem. A técnica aplicada foi a de modelagem matemática por células, onde os leitos catalíticos foram subdivididos em reatores tipo CSTR dinâmicos associados em série. A cada célula, um cálculo de flash foi associado, aperfeiçoando os balanços de massa e energia comumente empregados em reatores de leito gotejante. A abordagem termodinâmica utilizada para prever o equilíbrio líquido-vapor foi a i - i, com a equação de estado SRK associada a parâmetros de interação binária específicos para a solubilidade de hidrogênio em gasolina de pirólise. Os parâmetros de interação binária são provenientes do pacote termodinâmico do simulador iiSE. O modelo do reator foi implementado no software EMSO (Environment for Modeling, Simulation and Optimization) consistindo em cerca de 9000 equações e variáveis. Os resultados obtidos com o modelo contruído foram similares aos reportados na literatura. A aplicação de modelagem por células mostrou-se não só aplicável mas também mais robusta do que as abordagens tradicionais que utilizam equações diferenciais ordinárias. A utilização da ferramenta EMSO para a modelagem por células mostrou-se ainda mais vantajosa ao permitir a avaliação do comportamento dinâmico do reator em algumas situações hipotéticas, mas que são bem comuns na indústria. / Modeling and simulation of three-phase reactors is a challenge for process simulators. Many academic studies have been developed in order to describe and predict the behavior of these equipments for both steady states or dynamic responses. There are several approaches applied trying to reflect better the phenomena involved in three-phase reactors, always considering the computational cost of each approach. In this work a selective hydrogenation reactor for pyrolysis naphtha was simulated based on published literature data.The reactor is trickle bed type and the reactions are considered pseudofirst order. The technique used is cell-network modeling, where the catalyst beds were subdivided into dynamics CSTR reactors linked in series. To each cell there is a calculation Flash associated, enhancing mass and energy balances commonly applied in TBR. The thermodynamic approach used to predict the vapor-liquid equilibrium was i- i with SRK equation of state, associated to specific binary interaction parameters for the solubility of hydrogen in pygas from the thermodynamic simulator package iiSE. The binary interaction parameters come from the thermodynamic package of IISE simulator. The reactor model was implemented in software EMSO (Environment for Modeling, Simulation and Optimization) and consists of 9,000 equations and variables approximately. The results obtained with the model constructed here were similar to those reported in the literature. The use of cell-network modeling proved to be not only applicable but also more robust than the traditional approaches that use ordinary differential equations. The use of EMSO tool for cell-network modeling proved to be even more advantageous because allows the evaluation of some dynamic behavior of the reactor for hypothetical situations but quite common in the industry.

O polietileno de baixa densidade (PEBD) é uma resina termoplástica fabricada através de processos a altas pressões, com reatores do tipo autoclave ou tubular. O cenário atual é de crescente demanda por PEBD produzido a custos competitivos, através de plantas de alto desempenho e capacidade, de modo que o desenvolvimento de processo e produto torna-se um fator chave para as empresas produtoras, sendo os simuladores de processo ferramentas de grande importância para esta finalidade. O objetivo desta dissertação é a implementação de um modelo matemático que represente reatores autoclave de produção de PEBD e a avaliação de sua capacidade preditiva através da comparação dos seus resultados numéricos com dados de literatura e dados de uma planta industrial. O modelo desenvolvido deverá ser capaz de descrever o comportamento do reator, através de seu perfil de temperatura, vazões de iniciador e conversão, e de algumas das principais propriedades do polímero produzido, através de suas massas molares médias. A modelagem dos reatores é composta de compartimentos genéricos de tanque agitado, implementada no simulador dinâmico de processos EMSO, utilizando cinética de literatura. Os balanços de massa e energia foram resolvidos simultaneamente para a obtenção de condições de operação factíveis, sendo necessário o uso de controladores para manter o ponto de operação no estado estacionário desejado. As propriedades são obtidas através da utilização do método dos momentos. É considerada no modelo a possibilidade de formação de duas fases. Parâmetros do modelo cinético e de mistura são ajustados para os dados de planta através de uma rotina de estimação, considerando duas resinas com características de processo e produto distintas. As predições do modelo para perfil de temperatura, vazões de iniciador, conversão e massas molares médias mostraram boa concordância com os dados de literatura e com os dados de planta industrial, para as duas resinas consideradas. / The low density polyethylene (LDPE) is a thermoplastic resin manufactured by high pressure processes, with autoclave or tubular type reactors. The present scenario is of growing demand for LDPE produced at competitive costs, through high performance and capacity plants, so that the process and product development becomes a key factor for the producers, and process simulators are high importance tools for this purpose. The objective of this dissertation is the implementation of a mathematical model which represents autoclave reactors for LDPE production and the evaluation of its predictive capacity through the comparison between their numerical results with literature data and with an industrial plant data. The model developed should be able to describe the reactor behavior, through its temperature profile, initiator flows and conversion, and the behavior of some of the main properties of the produced polymer, through their average molecular weights. The reactors modeling is composed of stirred tank generic compartments, implemented in the process dynamic simulator EMSO, using literature kinetics. The mass and energy balances were solved simultaneously in order to obtain feasible operating conditions, which requires the use of controllers to maintain the operating point at the desired steady state. The properties are obtained through the use of the method of moments. It is considered in the model the possibility of formation of two phases. Kinetics and mixture model parameters are fitted to the plant data through an estimation routine, considering two grades with distinct process and product characteristics. The model predictions for temperature profile, initiator flows, conversion and average molecular weights presented good agreement with the literature data and with the industrial plant data, for the two considered grades.