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.

Este trabalho visa o estudo da viabilidade da redução do uso de água no processo de dessalgação em refinarias de petróleo. Em uma primeira fase, foi necessário o estudo teórico da separação das emulsões água/óleo. Em seguida, foi desenvolvido um modelo matemático baseado nas forças atuantes nas gotas de água, o que possibilitou a determinação do tempo entre as colisões de pares de gotas e o estabelecimento do critério para que ocorra o fenômeno de coalescência. Esse modelo foi empregado em um sistema desenvolvido com base em autômatos celulares, o qual possibilitou o acompanhamento do processo micro e macroscópico, através do cálculo para o conjunto das gotas, e o acompanhamento visual até a separação da fase contínua. Os experimentos de laboratório, para os quais foi usado equipamento ótico para a medição da intensidade de luz transmitida ou espalhada pelas gotas, possibilitaram avaliar a influência da qualidade da água de mistura no tempo de separação das emulsões. Na unidade industrial, foram realizados testes que permitiram analisar o desempenho das dessalgadoras em diferentes situações operacionais. Os resultados obtidos através dos experimentos de laboratório e da simulação usando o modelo matemático desenvolvido mostraram-se compatíveis com os dados obtidos nos testes na unidade industrial. O trabalho mostrou ser possível alterar os esquemas de usos de água nas dessalgadoras, aumentando a taxa de reciclagem e possibilitando a otimização do consumo de água fresca neste processo, o que resultaria em redução substancial no consumo geral de água na refinaria. / The aim of this work is the study of the reduction of water consumption in petroleum desalting processes. The study of the attraction forces acting on the droplets was necessary to know how the emulsion water/oil is separated. A mathematical model based upon these forces was built to calculate the time between each droplets collision and to establish criteria for their coalescence. This model was applied to a system developed based on cellular automata, which allows to follow the process micro and macroscopically. Computations were carried out to the ensemble of droplets and the visual progression, from the start of droplets separation of the continuous phase to the end of the process could be visualized. Laboratory experiments, in which optical equipment was used to measure the light intensity transmitted or scattered by the droplets, allowed to evaluate the influence of the type of mixing water in the separation time of the emulsions. Tests in the industrial unity allowed evaluating the performance of the desalting units at different operating conditions. Conclusions of the laboratory experiments and the results of the mathematical model were compared with results of the industrial tests, showing coherence between them. The work shows that it is possible to simulate the effect of the operating variables and to alter schemes of water use in desalting units, increasing the water recycling rate, allowing optimization of fresh water consumption in this process and reducing the total water consumption in the refinery.

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.

Dissertation submitted for obtain the degree of Doctor in Engineering Sciences, at the University of Portounder supervisory Professor Alírio E. Rodrigues

Dissertação apresentada para obtenção do grau de Doutor em Engenharia Química, na Faculdade de Engenharia da Universidade do Porto, sob a orientação do Professor Carlos Albino Veiga da Costa

Controladores PID são largamente utilizados no controle de processos industriais, tanto em sistemas monovariáveis como em sistemas multivariáveis. No entanto, muitos dos controladores encontrados na indústria são mal sintonizados. Um dos métodos mais simples de sintonia de controladores PID consiste em identificar algumas grandezas, as quais se relacionam com características do processo a controlar, e em seguida aplicar fórmulas para os parâmetros dos controladores baseadas nestas grandezas. Estas grandezas são o ganho e o período críticos do processo, os quais se relacionam diretamente com o limite de estabilidade do sistema. Uma característica bastante interessante deste método é que o mesmo pode facilmente ser implementado por um controle auto-ajustável. Por este fato, métodos de auto-ajuste destes controladores têm sido bastante utilizados em sistemas monovariáveis, através do uso do método do relé como procedimento de obtenção das grandezas críticas, necessárias à sintonia dos controladores. O método do relé consiste em aplicar um controle bang-bang em malha fechada com o processo do qual se quer identificar as grandezas críticas. Este procedimento, dentro de certas condições, provoca uma oscilação sustentada na saída do processo, da qual as grandezas críticas são obtidas. Tendo como objetivo o auto-ajuste de controladores PID em sistemas multivariáveis, o método do relé também pode ser utilizado para obtenção das grandezas críticas, devidamente adequado de maneira a obter as grandezas críticas multivariáveis do processo. Diferentes procedimentos de ensaios com relés podem ser aplicados em sistemas multivariáveis, porém somente um destes procedimentos é capaz de identificar as grandezas críticas multivariáveis do processo diretamente: o ensaio descentralizado com relés. Mesmo assim, a sintonia proposta na literatura é baseada nas fórmulas monovariáveis de Ziegler-Nichols e afins, mostrando-se, muitas vezes, inadequada. Desta forma, este trabalho tem o objetivo de apresentar um método de ajuste multivariável para controladores PID descentralizados, baseado nas grandezas críticas do processo. Este método estende o método do ponto crítico para sistemas monovariáveis aos sistemas multivariáveis, através da análise multivariável do problema. A análise do método do ponto crítico para sistemas monovariáveis mostra que um controlador PI ou PID sintonizado pelas fórmulas baseadas nas grandezas críticas do processo sempre irá deslocar o ponto crítico do processo para um outro ponto no plano complexo, determinado pelas fórmulas utilizadas. Da mesma forma, o método de sintonia proposto neste trabalho desloca o ponto crítico do processo para outro ponto no plano complexo, escolhido a priori, alterando a freqüência crítica do sistema. / PID controllers are widely used in process control, in singlevariable systems as well as in multivariable ones. Yet, many of the controllers found in industry are poorly tuned. One of the simplest tuning method of PID controllers consists in identifying some values which are related to the process characteristics, and simply apply some formulae based on these quantities to determine the parameters of the controllers. Theses quantities are the ultimate gain and the ultimate period of the process, which are directly related to the system stability limit. A very interesting characteristic of this method is that it is easily implemented by an auto-tuning control. Thus, auto-tuning methods of this kind of controllers have been largely used in singlevariable systems, using the relay feedback experiment in order to obtain the ultimate quantities, which are needed to tune the controllers. The relay feedback experiment consists in applying a bang-bang control to the process from which the ultimate quantities are to be identified. This procedure, under some conditions, provides a sustained oscillation in the process’ output, from which the ultimate quantities are obtained. Aiming at auto-tuning of PID controllers in multivariable systems, the relay feedback experiment can also be used in order to get the ultimate quantities. Different relay feedback procedures can be applied to multivariable processes, but only one of these can identify the real multivariable ultimate quantities, formally considering the multivariable nature of the process: the decentralized relay feedback (DRF). However, the tuning of the controllers proposed in the literature is based on Ziegler-Nichols like formulae, what seems to be, many times, inappropriate. This work presents a multivariable tuning method of decentralized PID controllers, based on the process’ ultimate quantities. This method extends the ultimate point method used in SISO systems to multivariable ones, through multivariable analysis of the problem. The analysis of the ultimate point method used in singlevariable systems shows that a PI or PID controller tuned through formulae based on the process’ ultimate quantities will always dislocate the ultimate point to another point in the complex plane, determined by the used formulae. The tuning method proposed in this work dislocates the process’ ultimate point of a multivariable process to another point in the complex plane, chosen a priori, modifying the system’s ultimate frequency.

Este trabalho constitui-se do estudo da cinética e termodinâmica da reação de cura de um sistema bi-epoxi ( Araldit F + anidrido ftálico + trietilamina utilizando a Análise Térmica. através da técnica do DSC. o sistema foi analisado segundo métodos não-isotérmicos com a finalidade de verificar os efeitos da concentração da trietilamina bem como o da velocidade de aquecimento. Os métodos utilizados foram o da regressão multi-linear31, de freeman-carrol de Kissinger 49,50 e de. Osawa51. Os resultados obtidos demonstraram Que a reação de polimerização-reticulação do sistema segue a equação de Arrhenius, indicando um processo aproximadamente de primeira ordem (1. 5 l em relação a amina, e que os valores da energia de ativação sofrem um aumento com a adição da amina ao meio reacional. Comparando-se o estudo!54 realizado anteriormente Araldit F + anidrido succínico + trietilaminal no Laboratório de Instrumentação e Dinâmica Molecular, deste Instituto, e o sistema em estudo ( Araldit F + anidrido ftálico + trietilamina ), verifica- se um comportamento cinético termodinâmico similar com diferenciação nos valores da energia de ativação e no fator pré-exponencial. sem diferenças nos valores da ordem de reação. A presença da amina causa um aumento considerável nos valores da energia de ativação nos dois sistemas. Este aumento da energia de ativação nos leva a considerar a amina como iniciador no processo de polimerização e não como um catalisador como comumente é considerada. / This work includes the Kinetic and Thermodynamical study by Thermal Analvsis, using DSC technics, of the cure reaction: Araldite F + Phtalic Anhydride + Triethylamine. The system was analysed by non-isothermal methods with the purpose to study the effects Df triethylamine concentration. as well the rate Df heating. The theoretical - mathematical treatment used here 31 was according to the multi linear regressl.on 43 Freeman-Carrol, kissinger 49,50 an d Osawa 51 techniques. The results demonstrated that the polymerization network-formation reaction Df the system studied followed the Arrhenius equation, indicating approximately a first order process (1.5) with refference to amine. The values Df the Activation Energy increased with addition of amine to reactional medium. The system studied in this work (Araldite F + Phtalic Anhydride + Triethylamine represents a Kinetic and Thermodynamical behaviour similar to the behaviour of the reaction system (Araldite F + Succinic Anhydride + Triethylamine) studied,4 before. We observed different values of Activation Energy as well as Pre-Exponencial Factor, but same reaction order for the two differents svstems. The amine addition leads to a remarkable increase Df the Activation Energv in the two systems. This increase in Activation Energy permited us to consider the amine as initiator Df the polymer-ization process, but not as catalvst as it is usuallv considerated.

Microalgas são micro-organismos fotossintetizantes que convertem CO2 majoritariamente em lipídios, proteínas e carboidratos. A produção de microalgas é indicada atualmente como alternativa para biomitigação de CO2 e para geração de biocombustíveis, a partir da conversão da biomassa através de processos químicos e biotecnológicos. As vantagens das microalgas estão na sua maior velocidade de crescimento em relação a vegetais oleaginosos superiores e em seu teor lipídico, que para algumas espécies pode chegar a valores maiores que 50%. Apesar do potencial envolvendo as microalgas, muitos desafios ainda precisam ser superados para tornar viável a aplicação desses micro-organismos para fins energéticos. São necessários desenvolvimentos nas mais diversas áreas, abordando os seguintes aspectos: seleção e melhoramento genético de espécies, alcance de maior eficiência fotossintética, desenvolvimento de sistemas de produção e de seu escalonamento, desenvolvimento de sistemas e processos de colheita, extração e processamento, aproveitamento de nutrientes provenientes de resíduos, otimização de condições operacionais, dentre outros. Este trabalho teve por objetivo determinar as melhores condições de crescimento da microalga marinha Nannochloropsis oculata em um fotobiorreator airlift. Os experimentos foram realizados seguindo um delineamento composto central rotacional com temperatura, (19 a 29 °C), concentração de nitrato no meio de cultivo (f/2) (25 a 125 mg.l-1) e intensidade luminosa (3636 a 10364 lux) como variáveis de estudo. Uma planta laboratorial com 12 fotobiorreatores foi construída, permitindo realizar os experimentos do referido planejamento em dois blocos. As variáveis de resposta estudadas e os respectivos melhores resultados experimentais obtidos foram de 482,7 mg.l-1 para a máxima concentração de biomassa, nas condições de 21 °C, 105 mg.l-1 de NaNO3 e 9000 lux, taxa instantânea de crescimento no período exponencial de 0,5624 d-1, nas condições de 24 °C, 75 mg.l-1 de NaNO3 e 7000 lux e percentual de lipídios em biomassa liofilizada de 30,36%, nas condições de 21 °C, 45 mg.l-1 de NaNO3 e 5000 lux. / Microalgae are photosynthetic microorganisms that convert CO2 mainly into lipids, proteins and carbohydrates. Microalgae production is currently pointed out as an alternative for CO2 biomitigation and generation of renewable biofuels, from biomass conversion by chemical and biotechnological processes. Microalgae advantages are fastest growth comparing with oil crops and their higher lipid content, which for some species can reach values above 50%. Despite the potential involving microalgae, many challenges remain to be overcome to make feasible the application of these microorganisms for energy purposes. We need to develop several areas, addressing the following issues: selection and genetic improvement of species, reaching higher photosynthetic efficiency, development of production systems and their scale up, development of systems and procedures for harvest, extraction and processing, use of nutrients from flue gases and wastewater, optimization of operating conditions, etc. This study aimed to determine the best conditions for growth of the marine microalgae Nannochloropsis oculata in airlift phootobioreactor. Experiments were carried out following a central composite design of the following variables: temperature (19 to 29 ° C), nitrate concentration in the culture medium (f/2) (25 to 125 mg.l-1) and irradiance (3636 to 10364 lux). A laboratory plant was built with 12 photobioreactors, allowing perform the central composite design in two blocks. The response variables studied and the best experimental results obtained were 482,7 mg.l-1 for maximum biomass concentration under conditions of 21 °C, NaNO3 concentration 105 mg.l-1 and 9000 lux, instantaneous rate of increase in exponential period equal to 0.5624 d-1 under the conditions of 24 °C, NaNO3 concentration 75 mg.l-1 and 7000 lux and percentual lipid content in lyophilized biomass of 30.36% under the conditions of 21 °C, NaNO3 concentration 45 mg.l-1 and 5000 lux.

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.

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.