Orientadores: Antonio Claudio Correa, Maria Cristina de Castro Cunha / Dissertação (mestrado) - Universidade Estadual de Campinas, Instituto de Geociencias / Made available in DSpace on 2018-07-21T15:42:26Z (GMT). No. of bitstreams: 1 Almeida_SergioRibeirode_M.pdf: 1892513 bytes, checksum: 921d8f84806cd70ffa02c33ea865d2b2 (MD5) Previous issue date: 1995 / Resumo: Esta dissertação estuda a aplicação da teoria dos canais de fluxo associada a um modelo linear de dupa porosidade em reservatórios heterogêneos submetidos à injeção de um fluido incompressível molhante e conseqüente produção do fluido não-molhante. Durante todo o trabalho foi considerado que as linhas de fluxo utilizadas são constantes com o tempo, não havendo fluxo entre os canais. Desprezou-se as forças gravitacionais e o fenômeno da embebição foi modelado através de um modelo de tranferência de massa já utilizado no modelo linear. Soluções analfticas obtidas através da teoria do potencial complexo são apresentadas para os arranjos de poços "five-spot" e linha esconsa, bem como procedimentos para obtenção das linhas de fluxo. A generalização do modelo linear permitiu sua aplicação para o caso bidimensional. Um outro item abordado foi a obtenção do coeficiente de difusão em função da saturação a partir da curva de pressão capilar / Abstract: This work studies the application of strearntube theory joined with a two-porosity linear model in heterogeneous reservoirs under injection of a incompressible wetting fluid and the production of a non-wetting fluid. It was considered that the streamlines are fixed and there is no flow between each channel. Gravity was neglected and the imbibition phenomena was modeled by the use of a mass transfer model. Analytical solutions given by complex potential theory are shown for two well pattem, as well as procedures to obtain the streamlines. The generalized linear model allowed its use in the two-dimensional case. Another item studied was how to obtain a variable diffusion coefficient as a function of saturation from capillary pressure data / Mestrado / Mestre em Geoengenharia de Reservatorios

Orientador: Amauri Garcia / Dissertação (mestrado) - Universidade Estadual de Campinas, Faculdade de Engenharia Mecanica / Made available in DSpace on 2018-07-22T00:57:27Z (GMT). No. of bitstreams: 1 Santos_CarlosAlexandredos_M.pdf: 8099685 bytes, checksum: 057758332ab5622072aab31f5866c60a (MD5) Previous issue date: 1997 / Resumo: A estrutura, e consequentemente as propriedades de uma peça fundida, dependem das condições térmicas durante a solidificação, ou mais precisamente das taxas de resmamento, que por sua vez dependem das resistências térmicas metal/molde e molde/ambiente. O objetivo principal deste trabalho consiste na análise do processo de solidificação através de ferramentas matemáticas e, por meio de um confronto teórico/experimental de perfis de temperatura, caracterizar o fenômeno de transferência de calor nas interfaces metal/molde e molde/ambiente. Na análise experimental, utilizaram-se sistemas que permitiram a retirada de calor na forma unidimensional e bidimensional, utilizando-se como liga de trabalho Sn 10% Pb e moldes de aço e cobre com diferentes espessuras e acabamento de usinagem convencional em suas paredes internas. Os perfis de temperaturas gerados foram confrontados com os perfis simulados através de um modelo numérico de diferenças finitas, determinando assim a variação das condutâncias térmicas nas interfaces. Obtidos os valores dos coeficientes de transferência de calor das interfaces em questão, aplicou-se o modelo na simulação da solidificaçãode peças de geometria simples sob a condição simultânea de diferentes taxas de retirada de calor em direções distintas, confrontando os resultados simulados com os experimentais, demostrando, desta forma, a eficiência e maleabilidade do modelo utilizado / Abstract: The microstructure and properties of castings depend on cooling conditions during solidification,or more precisely, on cooling rates and consequently on the thermal resistance of metal/mould and mould/environment interfaces. The main objective of this work consists on analysing the solidification process by using mathematical tools, and on the characterisation of heat flow at these interfaces by fitting experimental and simulated temperature profiles. The experiments were carried out in unidirectional and bidimensional conditions, by using Sn 10% Pb as the standard material, and copper and steel moulds with variable thickness. The experimental temperature profiles were fitted to those simulated by a numerical model based on a Finite Difference Technique, and the variations of interface heat transfer coefficients were determined. By using these heat transfer coefficients simulations were performed by the numerical model in order to predict the progress of solidification under conditions of different directional cooling rates. The model has proved to be efficient and versatile / Mestrado / Materiais e Processos de Fabricação / Mestre em Engenharia Mecânica

Orientador: Rubens Maciel Filho / Dissertação (mestrado) - Universidade Estadual de Campinas, Faculdade de Engenharia Quimica / Made available in DSpace on 2018-07-22T05:33:25Z (GMT). No. of bitstreams: 1 Azeredo_MariaCristinaRibeirode_M.pdf: 5028037 bytes, checksum: 00e54e9e7dd3f6688242feda01cb51a8 (MD5) Previous issue date: 1996 / Resumo: A grande utilização dos reatores catalíticos de leito fixo na indústria química, associado com a necessidade de operação de alta performance requerendo o uso de técnicas de otimização e controle avançado tem motivado o desenvolvimento e o uso de modelos detalhados. As características heterogêneas, causadas devido a presença do sólido e ao fato do sistema ser de parâmetros distribuídos e as equações serem não-lineares, quando analisadas do ponto de vista dinâmico, podem levar a ocorrência da resposta inversa. Sem um modelo detalhado, este tipo de comportamento não pode ser predito de uma maneira segura, o que torna dificil a definição de políticas operacionais ótimas e de controle. A aplicação de um modelo heterogêneo poderia ser uma solução, mas o tempo computacional pode ser proibitivo para aplicação em tempo real. Um modelo dinâmico, bidimensional, pseudo-homogêneo modificado que procura levar em conta, de maneira efetiva, a presença do sólido, é apresentado. A reação de obtenção do acetaldeído através da oxidação do etanol sobre catalisador de Fe-Mo é utilizada como caso estudo, embora o modelo matemático determinístico para representar os fenômenos que ocorrem em um reator de leito fixo não dependa da equação da taxa. O fato da reação ser fortemente exotérmica exige a utilização de reatores resfriados por fluidos refrigerantes, e isto tem de ser considerado no desenvolvimento do modelo matemático do reator. Para a solução das equações do modelo dinâmico, propõe-se um método que consiste em aplicar o método da colocação ortogonal na direção radial e o método de Runge-Kutta de quarta ordem na direção axial e na variável tempo. Entretanto, o modelo dinâmico exige o conhecimento prévio do perfil das variáveis de estado no instante inicial. Deste modo, o modelo estacionário é apresentado como estudo para definir critérios para o modelo dinâmico. Uma extensiva análise sensitiva do reator é feita de modo a se definir critérios seguros de operação, e para isso é estudada a influência das variáveis de operação no comportamento dinâmico do reator / Abstract: The development and use of detailed models have been motivated by the great application of a fixed bed catalytic reactor in Chemical Industry and the need of high profile operations that require optimization and advanced controI. The heterogeneous features, caused by the presence of solid particles and due to the fact that this is a distributed parameter system and non-linearity of the equations, when analyzed from a dynamic point of view, may cause to inverse response. Without a detailed model, this behavior can not be predicted in an accurate fashion, therefore the optimum operational and control policies become difficult to be determined. The solution could be the application of a heterogeneous model, however, the running time might be prohibitive for a real time application. A dynamic, two-dimensional, modified pseudo-homogeneous model that considers, effectively, the presence of solids particles is presented. The reaction of obtention of acetaldehyde by the oxidation of ethanol over a Fe-Mo catalyst is used as a study case, although the deterministic mathematical mode! that represents the phenomena that occurs in a fixed bed catalytic reactor does not depend on rate equation. It is necessary to use cooling reactors because the reaction is highly exothermic. This fact was considered in the mathematical model of the reactor. The method proposed to solve the equation of the dynamic model consists on the application of orthogonal collocation for radial direction and fourth-order Runge-Kutta for the axial as well as for the time. However, the dynamic model requires the previous knowledge of state variables profile at the initial time. Therefore, the stationary model is presented as a study to define the criteria for the dynamic model. An extensive analysis of the reactor is performed in order to define accurate operational criteria, and therefore the influence of operational variables in the dynamic behavior in the reactor is studied. / Mestrado / Desenvolvimento de Processos Químicos / Mestre em Engenharia Química

Orientador: Wagner Caroli do Amaral / Dissertação (mestrado) - Universidade Estadual de Campinas, Faculdade de Engenharia Eletrica e de Computação / Made available in DSpace on 2018-07-22T14:45:55Z (GMT). No. of bitstreams: 1 Campello_RicardoJoseGabrielliBarreto_M.pdf: 4071054 bytes, checksum: eb90b0948ad88acfde6bd0b7d4e10f07 (MD5) Previous issue date: 1997 / Resumo: Modelos nebulosos constituem uma importante ferramenta matemática na identificação de processos complexos onde técnicas convencionais normalmente não são eficientes. Uma classe particular é formada pelos modelos relacionais nebulosos, onde o mapeamento entre variáveis de entrada e saída é realizado através de equações relacionais nebulosas. Os principais procedimentos para modelagem de sistemas utilizando esses modelos são a determinação da estrutura, a estimação dos parâmetros da relação nebulosa e a construção das interfaces de entrada e saída quando são necessárias. Neste trabalho é proposto um algoritmo para estimação dos parâmetros da relação nebulosa em modelos relacionais de estrutura conhecida. A principal vantagem deste método, em comparação aos algoritmos existentes na literatura, é sua maior capacidade de adaptação. Essa característica evita a dependência de uma boa estimativa da relação nebulosa e torna o algoritmo adequado à modelagem de sistemas variantes no tempo. Outra característica é a utilização de um método de otimização simples que permite a identificação de sistemas nebulosos e não-nebulosos. Para construção das interfaces de entrada e saída utiliza-se um algoritmo baseado em otimização que reduz a influência das interfaces no modelo através da minimização dos erros nos mapeamentos entre grandezas numéricas e linguísticas, mantendo a integridade semântica dos conjuntos nebulosos. A otimização das interfaces de entrada e saída e a estimação dos parâmetros da relação nebulosa constituem a metodologia de identificação... Observação: O resumo, na íntegra, poderá ser visualizado no texto completo da tese digital / Abstract: Not informed. / Mestrado / Mestre em Engenharia Elétrica

Orientador: Rubens Maciel Filho / Tese (doutorado) - Universidade Estadual de Campinas, Faculdade de Engenharia Quimica / Made available in DSpace on 2018-07-23T09:09:51Z (GMT). No. of bitstreams: 1 Tomaz_Edson_D.pdf: 12844889 bytes, checksum: 3a2c4ca1dcf2aaa80c422e96c05065d6 (MD5) Previous issue date: 1998 / Resumo: Os fomos rotativos têm tido importantes aplicações em vários processos, especialmente para produção de cimento e incineração de resíduos sólidos perigosos. O Forno rotativo é o tipo de incinerador de resíduos perigosos mais empregado devido a sua capacidade de tratar materiais sólidos com uma ampla faixa de tamanho de partículas e proporcionar uma boa mistura do leito de sólido. Vários trabalhos sobre transferência de calor em fornos rotativos foram publicados com aplicações em engenharia metalúrgica e produção de cimento. Os trabalhos mais importantes foram desenvolvidos por Jenkin and Moles11, Gorog et all 13,14,15 e Silcox e Pershing6. Este último, desenvolveu o modelo mais realístico e constitue o ponto de partida para o presente trabalho. Neste trabalho, desenvolveu-se um modelo matemático para predizer os perfis axiais de temperatura e concentração dentro do forno rotativo e do pós-combustor do incinerador, em estado estacionário e não estacionário. Este modelo difere dos trabalhos anteriores principalmente porque considera a redução do volume do sólido devido à degradação térmica do resíduo, assim como, a combustão em fase gasosa dos compostos voláteis gerados. Além disso, a injeção de resíduo líquido no pós-combustor é permitida. O modelo é unidimensional e considera as trocas radiativas entre zonas imediatamente adjacentes. As correntes de fluxo de massa são divididas em jato ou chama, gás anular e sólido. O forno é dividido axialmente em zonas de iguais comprimentos, onde cada uma das correntes é considerada completamente misturada e caracterizada por uma única temperatura, única composição e única vazão mássica. A combustão dos voláteis oriundos do sólido em fase gasosa é considerada instantânea. A taxa de queima do resíduo é controlada pela transferência de massa dos voláteis a partir do sólido, expressa como taxa de perda de massa do sólido. Dentro de cada zona, o jato ou chama é representado como um cilindro de diâmetro constante e especificado para cada zona. O jato é dividido em três regiões: a região do estabelecimento do jato, a chama e a região pós-chama. Na região de estabelecimento do jato não ocorre combustão ou entrada de gás a partir da região anular. Na região pós-chama não há combustível auxiliar para combustão, mas pode haver a queima de voláteis oriundos do leito de sólido. Na chama, há a entrada de gás proveninente da região anular, cuja taxa é especificada para cada zona. A temperatura da última seção transversal do forno rotativo é calculada e não preestabelecida como nos trabalhos anteriores. As capacidades caloríficas, condutividades térmicas e viscosidades são calculadas levando-se em conta as dependência com a temperatura e a composição. Para resolver as equações de balanço de energia é adotado um procedimento matemático iterativo, e para a simulação numérica é desenvolvido um código FORTRAN. Uma análise de sensibilidade paramétrica é realizada para 21 parâmetros e os efeitos sobre os perfis de temperaturas são analisados, para variações negativas e positivas de 10%. As capacidades caloríficas do gás e do jato resultaram nos maiores efeitos. Devido à enorme variabilidade dos resíduos industriais, é analisada a influência de várias taxas de perda de massa sobre os perfis axiais de temperaturas. O trabalho também inclue a análise dos efeitos das taxa de entrada de gás na chama e da aproximação para a geometria da chama sobre o perfil de temperatura, assim como, um breve estudo sobre a redução do consumo de combustível por duas vias: pré-aquecimento do ar de combustão e enriquecimento com oxigênio / Abstract: Not informed / Doutorado / Desenvolvimento de Processos Químicos / Doutor em Engenharia Química

Orientadores: Sergio Persio Ravagnani, Lucia H. I. Mei / Dissertação (mestrado) - Universidade Estadual de Campinas, Faculdade de Engenharia Quimica / Made available in DSpace on 2018-07-23T10:54:20Z (GMT). No. of bitstreams: 1 Moncao_AnaPaulados_M.pdf: 2013063 bytes, checksum: c53abd549d973c08134f44a47800b23e (MD5) Previous issue date: 1997 / Resumo: Neste trabalho, hidrogéis de poli(2-hidroxietil metacrilato) ¿ PHEMA foram preparados por reticulação química, usando N,N¿ metileno ¿ bis ¿ acrilamida como agente reticulante. Hidrogéis de diferentes graus de reticulação foram obtidos variando-se a concentração do reticulante de 0,12% a 3% (p/p). Os filmes de PHEMA obtidos foram caracterizados através de análise de calorimetria diferencial de varredura ¿ DSC, microscopia eletrônica de varredura ¿ MEV e ensaio de tração. Dentre essas análises, o ensaio de tração foi o que mais mostrou as alterações das propriedades dos hidrogéis obtidos. O intumescimento foi estimado a partir da variação ocorrida no peso das amostras quando atingiram o equilíbrio, imersas em água destilada e solução tampão a pH=1, à temperatura ambiente. Observou-se que há influencia do meio externo no intumescimento do hidrogel. Os estudos da difusão foram realizados utilizando-se Rodamina B como espécie difusiva. As análises de difusão foram realizadas a ¿37 GRAUS¿C, em solução tampão de pH=1 e água destilada, pH~7. o coeficiente de difusão do soluto (D) foi calculado através da '2 POT. a¿ lei de Fick. Os dados experimentais obtidos e a elaboração de um modelo matemático permitiram a avaliação do coeficiente de difusão, utilizando-se métodos numéricos iterativos, os resultados obtidos mostraram que há reticulação física na formação do hidrogel de PHEMA, permitindo a sua utilização sem agente reticulante. Com a adição do agente reticulante, observou-se uma significativa alteração do coeficiente de difusão / Abstract: In this work, hydrogels of poly (2 ¿ hydroxyethyl methacrylate) ¿ PHEMA were prepared by chemical crosslinks using N, N¿methylene ¿ bis ¿ acrylamide as the crosslinking agent. Hydrogels with different crosslinking degrees were obtained by varying the crosslinker concentration from 0.12% up to 3.0% (w/w). the films of PHEMA obtained were characterized by differential scanning calorimetry ¿ DSC, scanning electron microscopy ¿ MEV and mechanical properties. The mechanical tensile test was the method of characterization that showed more properties variation of the hydrogels obtained. The swelling was estimated by the variation in the samples weight when they reached the equilibrium, after have been immersed into distilled water and buffer solution at room temperature. It was observed that there is influence of the pH in the hydrogel swelling. The diffusion studies were carried out using Rhodamine B as the diffusive model specimen. The diffusion analysis were conducted at ¿37 GRAUS¿C, in buffer solution at pH=1 and distilled water at pH~7. The Fick¿s Second Law has been used to calculate the solute diffusion coefficient (D). The experimental data obtained and the mathematical model elaborated permited the evaluation of diffusion coefficient using iterative numerical methods. The results obtained have shown that there is physical crosslinking in the PHEMA hydrogel formation, so its utilization without crosslinker was permited. A significative variation of the diffusion coefficient was observed when the crosslinking agent was added. / Mestrado / Ciencia e Tecnologia de Materiais / Mestre em Engenharia Química

Orientador: João Frederico A. C. Meyer / Dissertação (mestrado) - Universidade Estadual de Campinas, Instituto de Matematica, Estatistica e Computação Cientifica / Made available in DSpace on 2018-07-23T10:36:46Z (GMT). No. of bitstreams: 1 Bernardes_Mateus_M.pdf: 2439026 bytes, checksum: 665f76c0d359c6a0bb45d5d3dfa1217d (MD5) Previous issue date: 1998 / Resumo: Não informado. / Abstract: Not informed. / Mestrado / Mestre em Matemática Aplicada

Orientador: Marcio Luiz de Souza-Santos / Dissertação (mestrado) - Universidade Estadual de Campinas, Faculdade de Engenharia Mecanica / Made available in DSpace on 2018-07-23T14:24:57Z (GMT). No. of bitstreams: 1 Costa_AlexandreMarconideSouzada_M.pdf: 847784 bytes, checksum: 7bc8fb092d72a04962a7f74ea4c89e9e (MD5) Previous issue date: 1998 / Mestrado

Orientador: Denis Jose Schiozer / Tese (doutorado) - Universidade Estadual de Campinas, Faculdade de Engenharia Mecanica / Made available in DSpace on 2018-07-23T15:11:14Z (GMT). No. of bitstreams: 1 Guedes_SolangedaSilva_D.pdf: 13102319 bytes, checksum: e8ce55e85a643438e655bad9ce693bd9 (MD5) Previous issue date: 1998 / Resumo: Os avanços atuais na área de caracterização de reservatórios são capazes de fornecer descrições de reservatórios que resultam em malhas de discretização com centenas de milhares a dezenas de milhões de blocos. Durante a simulação numérica de fluxo multifásico, estas informações podem ser consideradas através de técnicas de mudança de escala. O procedimento usual considera pseudo-funções multifásicas e o método proposto por Kyte e Berry (1975) é largamente utilizado. Esta tese apresenta um modelo computacional multiescala para fluxo multifásico que trata o processo de mudança escala de forma implícita e sem o uso das pseudo-funções. Na metodologia proposta, urna seqüência de problemas locais são resolvidos em subdomínios da malha fina, obtendo assim urna melhor aproximação na malha grossa, mas com um custo computacional menor do que a resposta da malha fina. A principal vantagem deste novo procedimento de mudança de escala é tomar tal tarefa um processo implícito durante a etapa de solução do problema, contornando alguns problemas de ordem prática durante aplicações das pseudo-funções tradicionais. Os resultados apresentados tratam problemas bidimensionais em meios homogêneos e heterogêneos e duas fases: óleo e água. Alguns exemplos são utilizados para comparar este procedimento e o programa PSEUDO, um módulo de mudança de escala baseado em pseudo-funções e disponível no simulador comercial ECLIPSE / Abstract: Advances in petroleum reservoir descriptions have provided an amount of data that can not be handled directly during numerical simulations. This detailed geological information must be incorporated into a coarser model during multiphase fluid flow simulations by means of some upscaling technique. The most used approach is the pseudo relative permeabilities and the more widely used is the Kyte and Berry method (1975). In this work, it is proposed a multi-scale computational model for multiphase tlow that implicitly treats the upscaling without using pseudo functions. By solving a sequence of local problems on subdomains of the refined scale it is possible to achieve results with a coarser grid without expensive computations of a fine grid model. The main advantage of this new procedure is to treat the upscaling step implicitly in the solution process, overcoming some practical difficulties related the use of traditional pseudo functions. Results of bidimensional two phase flow simulations considering homogeneous and heterogeneous porous media are presented. Some examples compare the results of this approach and the commercial upscaling program PSEUDO, a module of the reservoir simulation software ECLIPSE / Doutorado / Doutor em Engenharia de Petróleo

Orientador: Maria Regina Wolf Maciel / Dissertação (mestrado) - Universidade Estadual de Campinas, Faculdade de Engenharia Quimica / Made available in DSpace on 2018-07-23T15:23:43Z (GMT). No. of bitstreams: 1 Perioto_FabianoRomero_M.pdf: 2598126 bytes, checksum: 3dc3c43bc391145773afca1ccf1f5e4d (MD5) Previous issue date: 1998 / Resumo: o presente trabalho representa urna contribuição à síntese do processo de destilação azeotrópica, o qual é utilizado para a separação de misturas líquidas com comportamento nãoideal, englobando a modelagem de coluna de destilação azeotrópica trifásica no estado estacionário. Os seguintes desenvolvimentos fazem parte desta tese: . Modelagem e Simulação de Separadores Simples Trifásicos : desenvolvimento da modelagem e simulação de separadores simples para sistemas multicomponentes, à temperatura e pressão fixas, englobando cálculos de equilíbrios líquido-líquido, líquidovapor e líquido-líquido-vapor, além de possibilitar também cálculos de pontos de bolha multifásicos. Esta metodologia é aplicável à caracterização terrnodinâmica de sistemas multicomponentes através das superficies de equilíbrio, sendo totalmente independente do fornecimento de estimativas iniciais; . Metodologia de Cálculo de Azeótropos Homogêneos e Heterogêneos para Sistemas Multicomponentes : desenvolvimento de nova metodologia de cálculo de pontos de azeotropia para sistemas multicomponentes englobando os casos homogêneos e heterogêneos. Esta metodologia é totalmente independente do fornecimento de.estimativas iniciais. Tal ferramenta é de suma importância na escolha de processos de separação adequados, particularmente orientando a escolha de solventes para a destilação azeotrópica; . Modelagem e Simulação de Colunas de Destilação Azeotrópica: desenvolvimento da modelagem e simulação de colunas de destilação azeotrópica aplicadas à separação de misturas multicomponentes, possibilitando o cálculo de equilíbrio líquido-líquido-vapor nos estágios internos da coluna. O programa apresenta urna metodologia de geração de estimativas dos perfis das variáveis de interesse e posterior identificação das regiões de equilíbrios líquido-líquido e líquido-líquido-vapor a partir do emprego de simuladores cOmerCiaIs / Abstract: This work represents a contribution to the synthesis of azeotropic distillation process, that is used to the separation of non-idealliquid mixtures, with the modelling of Steady-State Three-Phase Azeotropic Distillation Column. The following developments were carried out in this thesis: . Modelling and Simulation of Three-Phase Flash: development of the modelling and simulation of Three-Phase Flash to multicomponent systems, at fixed temperature and pressure, with liquid-liquid, liquid-vapor, liquid-liquid-vapor equilibrium calculations and multiphase bubbIe point calcuIations. This methodology is applicable in the thermodynamic characterization of the multicomponent systems through equilibrium surfaces, being totally independent of the initial estimates suppIy; . Methodology of Calculation of Heterogeneous and Homogeneous Azeotropes for MuIticomponent Systems: development of a new methodology of azeotropes calculations to multicomponents systems, with homogeneous and heterogeneous cases, independently of initial estimates supply. This tool is. very important for choosing adequate separation processes, especially to guide the choice of solvents to the azeotropic distillation application; . Modelling and Simulation of Azeotropic Distillation Columns: development of the modelling and simulation of the azeotropic distillation columns applicabIe to multicomponent systems separation, enabling the liquid-liquid-vapor equilibrium calculation in the internal stages of the column. The program presents a methodology to generation of estimates to the profiles of the interest variables and identification of the liquid-liquid and liquid-liquid-vapor equilibrium regions through commercial simulators / Mestrado / Mestre em Engenharia Química